SEARCH · Search NASA

Results for “gene expression analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Stage-resolved gene regulatory network analysis reveals developmental reprogramming and genes with robust stem-preferred expression in sorghum

Sorghum bicolor is a deep-rooted, heat- and drought-tolerant crop that thrives on marginal lands and is increasingly valued for its applications in biofuel, bioenergy, and biopolymer production. The sorghum stem, which can reach 4–5 m in length, serves as the primary reservoir of both lignocellulosic biomass and soluble sugars, making it a promising bioenergy feedstock. Although recent advances in genetic, genomic, and transcriptomic resources have improved our understanding of sorghum biology, comprehensive genome-wide analyses of functional dynamics across diverse organ types and developmental stages remain limited. In particular, candidate genes with stem preferred expression pattern or their associated cis-regulatory elements, which may program key stem-related functions and enable organ- or tissue-specific engineering, have not yet been identified.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Genomic Analysis of Aspergillus Section Terrei Reveals a High Potential in Secondary Metabolite Production and Plant Biomass Degradation

Aspergillus terreus has attracted interest due to its application in industrial biotechnology, particularly for the production of itaconic acid and bioactive secondary metabolites. As related species also seem to possess a prosperous secondary metabolism, they are of high interest for genome mining and exploitation. Here, we present draft genome sequences for six species from Aspergillus section Terrei and one species from Aspergillus section Nidulantes. Whole-genome phylogeny confirmed that section Terrei is monophyletic. Genome analyses identified between 70 and 108 key secondary metabolism genes in each of the genomes of section Terrei, the highest rate found in the genus Aspergillus so far. The respective enzymes fall into 167 distinct families with most of them corresponding to potentially unique compounds or compound families. Moreover, 53% of the families were only found in a single species, which supports the suitability of species from section Terrei for further genome mining. Intriguingly, this analysis, combined with heterologous gene expression and metabolite identification, suggested that species from section Terrei use a strategy for UV protection different to other species from the genus Aspergillus. Section Terrei contains a complete plant polysaccharide degrading potential and an even higher cellulolytic potential than other Aspergilli, possibly facilitating additional applications for these species in biotechnology.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2023

Galactomannan utilization by Cellvibrio japonicus relies on a single essential α‐galactosidase encoded by the aga27A gene

Plant mannans are a component of lignocellulose that can have diverse compositions in terms of its backbone and side-chain substitutions. Consequently, the degradation of mannan substrates requires a cadre of enzymes for complete reduction to substituent monosaccharides that can include mannose, galactose, and/or glucose. One bacterium that possesses this suite of enzymes is the Gram-negative saprophyte Cellvibrio japonicus, which has 10 predicted mannanases from the Glycoside Hydrolase (GH) families 5, 26, and 27. Here we describe a systems biology approach to identify and characterize the essential mannan-degrading components in this bacterium. The transcriptomic analysis uncovered significant changes in gene expression for most mannanases, as well as many genes that encode carbohydrate active enzymes (CAZymes) when mannan was actively being degraded. A comprehensive mutational analysis characterized 54 CAZyme-encoding genes in the context of mannan utilization. Growth analysis of the mutant strains found that the man26C, aga27A, and man5D genes, which encode a mannobiohydrolase, α-galactosidase, and mannosidase, respectively, were important for the deconstruction of galactomannan, with Aga27A being essential. Our updated model of mannan degradation in C. japonicus proposes that the removal of galactose sidechains from substituted mannans constitutes a crucial step for the complete degradation of this hemicellulose.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Regional epithelial cell diversity in the small intestine of pigs

Abstract Understanding regional distribution and specialization of small intestinal epithelial cells is crucial for developing methods to control appetite, stress, and nutrient uptake in swine. To establish a better understanding of specific epithelial cells found across different regions of the small intestine in pigs, we utilized single-cell RNA sequencing (scRNA-seq) to recover and analyze epithelial cells from duodenum, jejunum, and ileum. Cells identified included crypt cells, enterocytes, BEST4 enterocytes, goblet cells, and enteroendocrine (EE) cells. EE cells were divided into two subsets based on the level of expression of the EE lineage commitment gene, NEUROD1. NEUROD1hi EE cells had minimal expression of hormone-encoding genes and were dissimilar to EE cells in humans and mice, indicating a subset of EE cells unique to pigs. Recently discovered BEST4 enterocytes were detected in both crypts and villi throughout the small intestine via in situ staining, unlike in humans, where BEST4 enterocytes are found only in small intestinal villi. Proximal-to-distal gradients of expression were noted for hormone-encoding genes in EE cells and nutrient transport genes in enterocytes via scRNA-seq, demonstrating regional specialization. Regional gene expression in EE cells and enterocytes was validated via quantitative PCR (qPCR) analysis of RNA isolated from epithelial cells of different small intestinal locations. Though many genes had similar patterns of regional expression when assessed by qPCR of total epithelial cells, some regional expression was only detected via scRNA-seq, highlighting advantages of scRNA-seq to deconvolute cell type-specific regional gene expression when compared to analysis of bulk samples. Overall, results provide new information on regional localization and transcriptional profiles of epithelial cells in the pig small intestine.

Agriculture↗

DOE OSTI2026

Compilation and utilization of a sorghum transcriptome compendium for gene regulatory network analysis and crop trait engineering

Sorghum bicolor (Sorghum) is a drought and heat tolerant C4 grass crop used to produce grain, forage, biofuels, and other bioproducts. Genetic improvement of sorghum hybrid crops is aided by a large and diverse germplasm, sorghum's diploid inbreeding genetics, and a relatively small genome that has facilitated genomic research. Over the past 20 years, the sorghum research community characterized the cytogenetic and recombinant landscapes of sorghum's 10 chromosomes, sequenced and annotated the sorghum genome, and used that information to identify genes/alleles that modulate flowering time, plant height, seed shattering, and other important traits. More recently, >1000 RNA-seq transcriptome profiles were collected from 15 sorghum genotypes to help understand the genetic basis of variation in growth and development of sorghum stems, tillers, roots, and leaves, and the regulation of biosynthetic pathways that produce epicuticular wax, dhurrin, and RFOs, compounds that contribute to sorghum's resilience. Transcriptome studies were designed to identify differentially expressed genes that are co-expressed during development or in response to a treatment to enable construction of gene regulatory networks. Co-expression and network analysis identified transcription factors and their cognate binding sites in target gene promoters and signaling pathways that modulate gene regulatory networks providing gene editing targets for further trait optimization. RNA-seq data from >20 experiments targeting sorghum organs, tissues, cell types, developmental stages, and responses to environmental conditions (i.e., diel, day-length, shading, water-deficit, temperature) has been compiled in a sorghum transcriptome compendium. The goal of this resource paper is to describe compendium content, accessibility, and a compendium data analysis pipeline and to illustrate the types of information that can be derived from the compendium with a focus on the elucidation of gene regulatory networks useful for guiding the improvement of sorghum traits through gene editing.

RNA-seq↗

DOE OSTI2022

Astronauts Plasma-Derived Exosomes Induced Aberrant EZH2-Mediated H3K27me3 Epigenetic Regulation of the Vitamin D Receptor

There are unique stressors in the spaceflight environment. Exposure to such stressors may be associated with adverse effects on astronauts' health, including increased cancer and cardiovascular disease risks. Small extracellular vesicles (sEVs, i.e., exosomes) play a vital role in intercellular communication and regulate various biological processes contributing to their role in disease pathogenesis. To assess whether spaceflight alters sEVs transcriptome profile, sEVs were isolated from the blood plasma of 3 astronauts at two different time points: 10 days before launch (L-10) and 3 days after return (R+3) from the Shuttle mission. AC16 cells (human cardiomyocyte cell line) were treated with L-10 and R+3 astronauts-derived exosomes for 24 h. Total RNA was isolated and analyzed for gene expression profiling using Affymetrix microarrays. Enrichment analysis was performed using Enrichr. Furthermore, transcription factor (TF) enrichment analysis using the ENCODE/ChEA Consensus TF database identified gene sets related to the polycomb repressive complex 2 (PRC2) and Vitamin D receptor (VDR) in AC16 cells treated with R+3 compared to cells treated with L-10 astronauts-derived exosomes. Further analysis of the histone modifications using datasets from the Roadmap Epigenomics Project confirmed enrichment in gene sets related to the H3K27me3 repressive mark. Interestingly, analysis of previously published H3K27me3–chromatin immunoprecipitation sequencing (ChIP-Seq) ENCODE datasets showed enrichment of H3K27me3 in the VDR promoter. Collectively, our results suggest that astronaut-derived sEVs may epigenetically repress the expression of the VDR in human adult cardiomyocytes by promoting the activation of the PRC2 complex and H3K27me3 levels.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

The JA‐responsive MYC2‐ BADH ‐ like transcriptional regulatory module in Poncirus trifoliata contributes to cold tolerance by modulation of glycine betaine biosynthesis

Summary Glycine betaine (GB) is known to accumulate in plants exposed to cold, but the underlying molecular mechanisms and associated regulatory network remain unclear. Here, we demonstrated that PtrMYC2 of Poncirus trifoliata integrates the jasmonic acid (JA) signal to modulate cold‐induced GB accumulation by directly regulating PtrBADH‐l , a betaine aldehyde dehydrogenase (BADH)‐like gene. PtrBADH‐l was identified based on transcriptome and expression analysis in P. trifoliata . Overexpression and VIGS (virus‐induced gene silencing)‐mediated knockdown showed that PtrBADH‐l plays a positive role in cold tolerance and GB synthesis. Yeast one‐hybrid library screening using PtrBADH‐l promoter as baits unraveled PtrMYC2 as an interacting candidate. PtrMYC2 was confirmed to directly bind to two G‐box cis ‐acting elements within PtrBADH‐l promoter and acts as a transcriptional activator. In addition, PtrMYC2 functions positively in cold tolerance through modulation of GB synthesis by regulating PtrBADH‐l expression. Interestingly, we found that GB accumulation under cold stress was JA‐dependent and that PtrMYC2 orchestrates JA‐mediated PtrBADH‐l upregulation and GB accumulation. This study sheds new light on the roles of MYC2 homolog in modulating GB synthesis. In particular, we propose a transcriptional regulatory module PtrMYC2‐PtrBADH‐l to advance the understanding of molecular mechanisms underlying the GB accumulation under cold stress.

Ming, Ruhong↗

DOE OSTI2022

Transcriptome Sequencing of Agave angustifolia Reveals Conservation and Diversification in the Expression of Cinnamyl Alcohol Dehydrogenase Genes in Agave Species

Agave angustifolia is an important crassulacean acid metabolism plant species, with wide applications in beverage and sisal fiber production. In this study, we carried out a transcriptome analysis of A. angustifolia leaves, generating a total of 58,482,436 clean reads through Illumina paired-end sequencing. De novo transcriptome assembly generated 67,314 unigenes, with about half of them having homologs in four public databases. In the Nr database, Asparagus officinalis was shown to be most closely related to agave, and the early angiosperm Amborella trichopoda was selected as an outgroup species. We further identified five, two, seven, seven, seven, six and six CAD genes in asparagus, amborella, A. deserti, A. tequilana, A. americana, A. hybrid H11648 and A. angustifolia, respectively. The maximum likelihood phylogenetic tree revealed the species-specific expansion of CAD genes in arabidopsis, rice and agave. The expression analysis indicated the conservatively expressed CAD1/2/4/6, providing candidate targets for manipulation to improve lignin traits. The species-specific expression of CAD3/5/7 indicates the existence of different regulatory mechanisms controlling the expression of these genes in agave species. This study presents the first transcriptome dataset of A. angustifolia, facilitating future studies on lignin biosynthesis in agave.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Omics-Lethal Human Viruses, Interferon-Stimulated Response Experiment IFNaHUH001

The purpose of this experiment was to evaluate the human host cellular response to treatment with and without interferon alpha/beta (IFNα/β) treatment. Sample time course data was obtained from human hepatoma carcinoma cells (HUH-7) expressing transcriptional activation protein VP30, stably transfected with the Ebola VP30 gene for enabling replication of VP30-deficient Ebola virus for transcriptome expression analysis. Secondary host-associated viral dataset downloads contain one or more statistically processed (normalization data transformation) quantitative dataset collections resulting in qualitative expression analyses of primary host-pathogen experimental study designs. Transcriptomics dataset download have a direct relationship to a primary sample submission corresponding to a specific interferon response to viral infection.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Phenotypically anchored transcriptomics across diverse agrichemicals reveals conserved pathways and unique gene expression signatures in zebrafish

Agrichemicals such as herbicides, fungicides, insecticides, and biocides are widely used in agriculture, yet some are associated with adverse effects in humans and the environment. While many of these chemicals have been extensively studied in vitro and are included in the EPA’s ToxCast program, comprehensive in vivo comparisons using RNA sequencing across structurally diverse agrichemicals, in a single screening platform, are lacking. In this study, we examined structurally diverse agrichemicals found in the U.S. Environmental Protection Agency’s (EPA) Toxcast Phase I and II library by statically exposing early life stage zebrafish at 6 h post fertilization (hpf) until 120 hpf at concentrations ranging from 0.25 to 100 µM. Morphological outcomes were assessed at 120 hpf across 10 endpoints, including yolk sac edema, craniofacial malformations, and axis abnormalities. Chemicals that produced robust concentration-response relationships were selected for transcriptomic profiling. For transcriptomic analysis, zebrafish were statically exposed to each chemical and sampled at 48 hpf, prior to the onset of morphological effects observed at 120 hpf. Differential expression analysis identified between 0 and 4,538 differentially expressed genes (DEGs) per chemical, with no clear correlation to morphological severity. Both DEG and co-expression network analyses revealed chemical-specific expression patterns that converged on shared biological pathways, including neurodevelopment and cytoskeletal organization. Key regulatory genes such as mylpfa and krt4 were identified within co-expression modules, suggesting their potential role in conserved toxicity mechanisms. Semantic similarity analysis of enriched gene ontology (GO) terms, when compared to existing datasets, highlighted gaps in the annotation of neurodevelopmental processes, indicating that some in vivo effects may not be fully captured by current curated resources. The results provide new insights into the modes of action of diverse agrichemicals and establish a framework for understanding how agrichemical structure relates to biological function in a vertebrate model.

agrichemical↗

DOE OSTI2022

Metatranscriptomics reveals a shift in microbial community composition and function during summer months in a coastal marine environment

Abstract Temperate coastal marine waters are often thermally stratified from spring through fall but can be dynamic and disrupted by tidal currents and wind‐driven upwelling. These mixing events introduce deeper, cooler water with a higher partial pressure of CO 2 (pCO2) and its associated microbial communities to the surface. Anecdotally, these events impact shellfish hatcheries and farms, warranting improved understanding of changes in composition and activity of marine microbial communities in relation to environmental processes. To characterize both compositional and functional changes associated with abiotic factors, here, we generate a reference metatranscriptome from the Strait of Georgia over representative seasons and analyze metatranscriptomic profiles of the microorganisms present within intake water containing different pCO 2 levels at a shellfish hatchery in British Columbia from June through October. Abiotic factors studied include pH, temperature, alkalinity, aragonite, calcite, and pCO 2 . Community composition changes were observed to occur at broad taxonomic levels and most notably to vary with temperature and pCO 2 . Functional gene expression profiles indicated a strong difference between early (June–July) and late summer (August–October) associated with viral activity. The taxonomic data suggest this could be due to the termination of cyanobacteria and phytoplankton blooms by viral lysis in the late season. Functional analysis indicated fewer differentially expressed transcripts associated with abiotic variables (e.g., pCO 2 ) than with the temporal effect. Microbial composition and activity in these waters vary with both short‐term effects observed alongside abiotic variation and long‐term effects observed across seasons. The analysis of both taxonomy and functional gene expression simultaneously in the same samples by environmental RNA (eRNA metatranscriptomics) provided a more comprehensive view for monitoring water bodies than either would in isolation.

Sutherland, Ben J. G.↗

DOE OSTI2024

Transcriptomic Analysis of Arachidonic Acid Pathway Genes Provides Mechanistic Insight into Multi-Organ Inflammatory and Vascular Diseases

Arachidonic acid (AA) metabolites have been associated with several diseases across various organ systems, including the cardiovascular, pulmonary, and renal systems. Lipid mediators generated from AA oxidation have been studied to control macrophages, T-cells, cytokines, and fibroblasts, and regulate inflammatory mediators that induce vascular remodeling and dysfunction. AA is metabolized by cyclooxygenase (COX), lipoxygenase (LOX), and cytochrome P450 (CYP) to generate anti-inflammatory, pro-inflammatory, and pro-resolutory oxidized lipids. As comorbid states such as diabetes, hypertension, and obesity become more prevalent in cardiovascular disease, studying the expression of AA pathway genes and their association with these diseases can provide unique pathophysiological insights. In addition, the AA pathway of oxidized lipids exhibits diverse functions across different organ systems, where a lipid can be both anti-inflammatory and pro-inflammatory depending on the location of metabolic activity. Therefore, we aimed to characterize the gene expression of these lipid enzymes and receptors throughout multi-organ diseases via a transcriptomic meta-analysis using the Gene Expression Omnibus (GEO) Database. In our study, we found that distinct AA pathways were expressed in various comorbid conditions, especially those with prominent inflammatory risk factors. Comorbidities, such as hypertension, diabetes, and obesity appeared to contribute to elevated expression of pro-inflammatory lipid mediator genes. Our results demonstrate that expression of inflammatory AA pathway genes may potentiate and attenuate disease; therefore, we suggest further exploration of these pathways as therapeutic targets to improve outcomes.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

A Novel Trichothecene Toxin Phenotype Associated with Horizontal Gene Transfer and a Change in Gene Function in Fusarium

Fusarium trichothecenes are among the mycotoxins of most concern to food and feed safety. Production of these mycotoxins and presence of the trichothecene biosynthetic gene (TRI) cluster have been confirmed in only two multispecies lineages of Fusarium: the Fusarium incarnatum-equiseti (Incarnatum) and F. sambucinum (Sambucinum) species complexes. Here, we identified and characterized a TRI cluster in a species that has not been formally described and is represented by Fusarium sp. NRRL 66739. This fungus is reported to be a member of a third Fusarium lineage: the F. buharicum species complex. Cultures of NRRL 66739 accumulated only two trichothecenes, 7-hydroxyisotrichodermin and 7-hydroxyisotrichodermol. Although these are not novel trichothecenes, the production profile of NRRL 66739 is novel, because in previous reports 7-hydroxyisotrichodermin and 7-hydroxyisotrichodermol were components of mixtures of 6–8 trichothecenes produced by several Fusarium species in Sambucinum. Heterologous expression analysis indicated that the TRI13 gene in NRRL 66739 confers trichothecene 7-hydroxylation. This contrasts the trichothecene 4-hydroxylation function of TRI13 in other Fusarium species. Phylogenetic analyses suggest that NRRL 66739 acquired the TRI cluster via horizontal gene transfer from a close relative of Incarnatum and Sambucinum. These findings provide insights into evolutionary processes that have shaped the distribution of trichothecene production among Fusarium species and the structural diversity of the toxins.

Food Science & Technology↗

DOE OSTI2020

Transcriptome Analysis of Drought-Resistant and Drought-Sensitive Sorghum ( Sorghum bicolor ) Genotypes in Response to PEG-Induced Drought Stress

Drought is a major limiting factor of crop yields. In response to drought, plants reprogram their gene expression, which ultimately regulates a multitude of biochemical and physiological processes. The timing of this reprogramming and the nature of the drought-regulated genes in different genotypes are thought to confer differential tolerance to drought stress. Sorghum is a highly drought-tolerant crop and has been increasingly used as a model cereal to identify genes that confer tolerance. Also, there is considerable natural variation in resistance to drought in different sorghum genotypes. Here, we evaluated drought resistance in four genotypes to polyethylene glycol (PEG)-induced drought stress at the seedling stage and performed transcriptome analysis in seedlings of sorghum genotypes that are either drought-resistant or drought-sensitive to identify drought-regulated changes in gene expression that are unique to drought-resistant genotypes of sorghum. Our analysis revealed that about 180 genes are differentially regulated in response to drought stress only in drought-resistant genotypes and most of these (over 70%) are up-regulated in response to drought. Among these, about 70 genes are novel with no known function and the remaining are transcription factors, signaling and stress-related proteins implicated in drought tolerance in other crops. This study revealed a set of drought-regulated genes, including many genes encoding uncharacterized proteins that are associated with drought tolerance at the seedling stage.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

A Genome-Scale Atlas Reveals Complex Interplay of Transcription and Translation in an Archaeon

The scale of post-transcriptional regulation and the implications of its interplay with other forms of regulation in environmental acclimation are underexplored for organisms of the domain Archaea. Here, we have investigated the scale of post-transcriptional regulation in the extremely halophilic archaeon Halobacterium salinarum NRC-1 by integrating the transcriptome-wide locations of transcript processing sites (TPSs) and SmAP1 binding, the genome-wide locations of antisense RNAs (asRNAs), and the consequences of RNase_2099C knockout on the differential expression of all genes. This integrated analysis has discovered that 54% of all protein-coding genes in the genome of this haloarchaeon are likely targeted by multiple mechanisms for putative post-transcriptional processing and regulation, with about 20% of genes likely being regulated by combinatorial schemes involving SmAP1, asRNAs, and RNase_2099C. Comparative analysis of mRNA levels (transcriptome sequencing [RNA-Seq]) and protein levels (sequential window acquisition of all theoretical fragment ion spectra mass spectrometry [SWATH-MS]) for 2,579 genes over four phases of batch culture growth in complex medium generated additional evidence for the conditional post-transcriptional regulation of 7% of all protein-coding genes. We demonstrate that post-transcriptional regulation may act to fine-tune specialized and rapid acclimation to stressful environments, e.g., as a switch to turn on gas vesicle biogenesis to promote vertical relocation under anoxic conditions and modulate the frequency of transposition by insertion sequence (IS) elements of the IS200/IS605, IS4, and ISH3 families. Findings from this study are provided as an atlas in a public Web resource (https://halodata.systemsbiology.net).

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2026

Data for The Stem Cell-Type Transcriptome of Bioenergy Sorghum Reveals the Spatial Regulation of Secondary Cell Wall Networks

Bioenergy sorghum is a low-input, drought-resilient, deep-rooting annual crop that has high biomass yield potential enabling the sustainable production of biofuels, biopower, and bioproducts. Bioenergy sorghum’s 4-5 m stems account for ~80% of the harvested biomass. Stems accumulate high levels of sucrose that could be used to synthesize bioethanol and useful biopolymers if information about stem cell-type gene expression and regulation was available to enable engineering. To obtain this information, Laser Capture Microdissection (LCM) was used to isolate and collect transcriptome profiles from five major cell types that are present in stems of the sweet sorghum Wray. Transcriptome analysis identified genes with cell-type specific and cell-preferred expression patterns that reflect the distinct metabolic, transport, and regulatory functions of each cell type. Analysis of cell-type specific gene regulatory networks (GRNs) revealed that unique TF families contribute to distinct regulatory landscapes, where regulation is organized through various modes and identifiable network motifs. Cell-specific transcriptome data was combined with a stem developmental transcriptome dataset to identify the GRN that differentially activates the secondary cell wall (SCW) formation in stem xylem sclerenchyma and epidermal cells. The cell-type transcriptomic dataset provides a valuable source of information about the function of sorghum stem cell types and GRNs that will enable the engineering of bioenergy sorghum stems.

Software↗

DOE OSTI2021

Integrative genomics reveals paths to sex dimorphism in Salix purpurea L

Sex dimorphism and gene expression were studied in developing catkins in 159 F 2 individuals from the bioenergy crop Salix purpurea , and potential mechanisms and pathways for regulating sex development were explored. Differential expression, eQTL, bisulfite sequencing, and network analysis were used to characterize sex dimorphism, detect candidate master regulator genes, and identify pathways through which the sex determination region (SDR) may mediate sex dimorphism. Eleven genes are presented as candidates for master regulators of sex, supported by gene expression and network analyses. These include genes putatively involved in hormone signaling, epigenetic modification, and regulation of transcription. eQTL analysis revealed a suite of transcription factors and genes involved in secondary metabolism and floral development that were predicted to be under direct control of the sex determination region. Furthermore, data from bisulfite sequencing and small RNA sequencing revealed strong differences in expression between males and females that would implicate both of these processes in sex dimorphism pathways. These data indicate that the mechanism of sex determination in Salix purpurea is likely different from that observed in the related genus Populus . This further demonstrates the dynamic nature of SDRs in plants, which involves a multitude of mechanisms of sex determination and a high rate of turnover.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Gene expression of functionally-related genes coevolves across fungal species: detecting coevolution of gene expression using phylogenetic comparative methods

Researchers often measure changes in gene expression across conditions to better understand the shared functional roles and regulatory mechanisms of different genes. Analogous to this is comparing gene expression across species, which can improve our understanding of the evolutionary processes shaping the evolution of both individual genes and functional pathways. One area of interest is determining genes showing signals of coevolution, which can also indicate potential functional similarity, analogous to co-expression analysis often performed across conditions for a single species. However, as with any trait, comparing gene expression across species can be confounded by the non-independence of species due to shared ancestry, making standard hypothesis testing inappropriate. We compared RNA-Seq data across 18 fungal species using a multivariate Brownian Motion phylogenetic comparative method (PCM), which allowed us to quantify coevolution between protein pairs while directly accounting for the shared ancestry of the species. Our work indicates proteins which physically-interact show stronger signals of coevolution than randomly-generated pairs. Interactions with stronger empirical and computational evidence also showing stronger signals of coevolution. We examined the effects of number of protein interactions and gene expression levels on coevolution, finding both factors are overall poor predictors of the strength of coevolution between a protein pair. Simulations further demonstrate the potential issues of analyzing gene expression coevolution without accounting for shared ancestry in a standard hypothesis testing framework. Furthermore, our simulations indicate the use of a randomly-generated null distribution as a means of determining statistical significance for detecting coevolving genes with phylogenetically-uncorrected correlations, as has previously been done, is less accurate than PCMs, although is a significant improvement over standard hypothesis testing. These methods are further improved by using a phylogenetically-corrected correlation metric. Our work highlights potential benefits of using PCMs to detect gene expression coevolution from high-throughput omics scale data. This framework can be built upon to investigate other evolutionary hypotheses, such as changes in transcription regulatory mechanisms across species.

59 BASIC BIOLOGICAL SCIENCES↗