Search NASA⌕ Search

SEARCH · Search NASA

Results for “co-expression networks”

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.

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

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Mapping of flumioxazin tolerance in a snap bean diversity panel leads to the discovery of a master genomic region controlling multiple stress resistance genes

Effective weed management tools are crucial for maintaining the profitable production of snap bean (Phaseolus vulgaris L.). Preemergence herbicides help the crop to gain a size advantage over the weeds, but the few preemergence herbicides registered in snap bean have poor waterhemp (Amaranthus tuberculatus) control, a major pest in snap bean production. Waterhemp and other difficult-to-control weeds can be managed by flumioxazin, an herbicide that inhibits protoporphyrinogen oxidase (PPO). However, there is limited knowledge about crop tolerance to this herbicide. We aimed to quantify the degree of snap bean tolerance to flumioxazin and explore the underlying mechanisms. We investigated the genetic basis of herbicide tolerance using genome-wide association mapping approach utilizing field-collected data from a snap bean diversity panel, combined with gene expression data of cultivars with contrasting response. The response to a preemergence application of flumioxazin was measured by assessing plant population density and shoot biomass variables. Snap bean tolerance to flumioxazin is associated with a single genomic location in chromosome 02. Tolerance is influenced by several factors, including those that are indirectly affected by seed size/weight and those that directly impact the herbicide's metabolism and protect the cell from reactive oxygen species-induced damage. Transcriptional profiling and co-expression network analysis identified biological pathways likely involved in flumioxazin tolerance, including oxidoreductase processes and programmed cell death. Transcriptional regulation of genes involved in those processes is possibly orchestrated by a transcription factor located in the region identified in the GWAS analysis. Several entries belonging to the Romano class, including Bush Romano 350, Roma II, and Romano Purpiat presented high levels of tolerance in this study. The alleles identified in the diversity panel that condition snap bean tolerance to flumioxazin shed light on a novel mechanism of herbicide tolerance and can be used in crop improvement.

60 APPLIED LIFE 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 OSTI2023

Altered Expression of Two Small Secreted Proteins (ssp4 and ssp6) Affects the Degradation of a Natural Lignocellulosic Substrate by Pleurotus ostreatus

Pleurotus ostreatus is a white-rot fungus that can degrade lignin in a preferential manner using a variety of extracellular enzymes, including manganese and versatile peroxidases (encoded by the vp1-3 and mnp1-6 genes, respectively). This fungus also secretes a family of structurally related small secreted proteins (SSPs) encoded by the ssp1-6 genes. Using RNA sequencing (RNA-seq), we determined that ssp4 and ssp6 are the predominant members of this gene family that were expressed by P. ostreatus during the first three weeks of growth on wheat straw. Downregulation of ssp4 in a strain harboring an ssp RNAi construct (KDssp1) was then confirmed, which, along with an increase in ssp6 transcript levels, coincided with reduced lignin degradation and the downregulation of vp2 and mnp1. In contrast, we observed an increase in the expression of genes related to pectin and side-chain hemicellulose degradation, which was accompanied by an increase in extracellular pectin-degrading capacity. Genome-wide comparisons between the KDssp1 and the wild-type strains demonstrated that ssp silencing conferred accumulated changes in gene expression at the advanced cultivation stages in an adaptive rather than an inductive mode of transcriptional response. Based on co-expression networking, crucial gene modules were identified and linked to the ssp knockdown genotype at different cultivation times. Based on these data, as well as previous studies, we propose that P. ostreatus SSPs have potential roles in modulating the lignocellulolytic and pectinolytic systems, as well as a variety of fundamental biological processes related to fungal growth and development.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Rhythmic Mechanisms Governing CAM Photosynthesis in Kalanchoe fedtschenkoi : High-Resolution Temporal Transcriptomics

Crassulacean acid metabolism (CAM) is a specialized photosynthetic pathway that enhances water-use efficiency by temporally separating nocturnal CO 2 uptake from daytime decarboxylation and carbon fixation. To uncover the regulatory mechanisms coordinating these temporal dynamics, we generated high-resolution, 48 h time-course transcriptomes for the CAM model Kalanchoe fedtschenkoi under both 12 h/12 h light/dark (LD) cycles and continuous light (LL). A rhythmicity analysis revealed that diel light cues are the dominant driver of transcript oscillations: 16,810 genes (54.3% of annotated genes) exhibited rhythmic expression only under LD, whereas just 399 genes (1.3%) remained rhythmic under LL. A smaller set of 3009 genes (9.7%) oscillated in both conditions, indicating that the intrinsic circadian clock sustains rhythmicity for a limited subset of the transcriptome. A gene co-expression network analysis revealed extensive integration between circadian clock components, core CAM pathway enzymes, and stomatal regulators, defining regulatory modules that coordinate metabolic and physiological timing. Notably, key hub genes associated with post-translational and post-transcriptional regulation, including the E3 ubiquitin ligase HUB2 and several pentatricopeptide repeat (PPR) proteins, act as central nodes in CAM-associated networks. This discovery implicates epigenetic and organellar regulation as previously unrecognized critical tiers of control in CAM. Together, our results support a regulatory model in which CAM rhythmicity is governed by both external light/dark cues and the endogenous circadian clock through multi-level control spanning transcriptional and protein-level regulation. To support community exploration, we also provide an interactive eFP (electronic Fluorescent Pictograph) browser for visualizing time-resolved gene expression profiles.

09 BIOMASS FUELS↗
DOE OSTI2024

Transcriptomic and Proteomic Insights into Host Immune Responses in Pediatric Severe Malarial Anemia: Dysregulation in HSP60-70-TLR2/4 Signaling and Altered Glutamine Metabolism

Severe malarial anemia (SMA, Hb < 6.0 g/dL) is a leading cause of childhood morbidity and mortality in holoendemic Plasmodium falciparum transmission zones. This study explored the entire expressed human transcriptome in whole blood from 66 Kenyan children with non-SMA (Hb ≥ 6.0 g/dL, n = 41) and SMA (n = 25), focusing on host immune response networks. RNA-seq analysis revealed 6862 differentially expressed genes, with equally distributed up-and down-regulated genes, indicating a complex host immune response. Deconvolution analyses uncovered leukocytic immune profiles indicative of a diminished antigenic response, reduced immune priming, and polarization toward cellular repair in SMA. Weighted gene co-expression network analysis revealed that immune-regulated processes are central molecular distinctions between non-SMA and SMA. A top dysregulated immune response signaling network in SMA was the HSP60-HSP70-TLR2/4 signaling pathway, indicating altered pathogen recognition, innate immune activation, stress responses, and antigen recognition. Validation with high-throughput gene expression from a separate cohort of Kenyan children (n = 50) with varying severities of malarial anemia (n = 38 non-SMA and n = 12 SMA) confirmed the RNA-seq findings. Proteomic analyses in 35 children with matched transcript and protein abundance (n = 19 non-SMA and n = 16 SMA) confirmed dysregulation in the HSP60-HSP70-TLR2/4 signaling pathway. Additionally, glutamine transporter and glutamine synthetase genes were differentially expressed, indicating altered glutamine metabolism in SMA. This comprehensive analysis underscores complex immune dysregulation and novel pathogenic features in SMA.

Microbiology↗
NASA NTRS2000

Genetic Network Inference: From Co-Expression Clustering to Reverse Engineering

Advances in molecular biological, analytical, and computational technologies are enabling us to systematically investigate the complex molecular processes underlying biological systems. In particular, using high-throughput gene expression assays, we are able to measure the output of the gene regulatory network. We aim here to review datamining and modeling approaches for conceptualizing and unraveling the functional relationships implicit in these datasets. Clustering of co-expression profiles allows us to infer shared regulatory inputs and functional pathways. We discuss various aspects of clustering, ranging from distance measures to clustering algorithms and multiple-duster memberships. More advanced analysis aims to infer causal connections between genes directly, i.e., who is regulating whom and how. We discuss several approaches to the problem of reverse engineering of genetic networks, from discrete Boolean networks, to continuous linear and non-linear models. We conclude that the combination of predictive modeling with systematic experimental verification will be required to gain a deeper insight into living organisms, therapeutic targeting, and bioengineering.

Dhaeseleer, Patrik↗
DOE OSTI2026

A tissue‐resolved, network‐based transcriptomic framework for abiotic stress responses in sorghum

Developing climate‐resilient crops requires a detailed understanding of stress‐induced gene expression dynamics, as maladaptive responses can compromise their productivity and survival. Sorghum, a globally important cereal with exceptional tolerance to multiple abiotic stresses, provides a powerful system for investigating these dynamics. However, how stress type, tissue specificity, and temporal progression jointly shape transcriptomic responses in crops remains poorly understood. Here, we present a comparative, time‐resolved transcriptomic atlas of sorghum responses to drought, heat, and salinity stress across shoot and root tissues. Integrative analyses revealed that tissue specificity is the dominant determinant of abiotic stress‐induced gene reprogramming across all three stresses. Building on these global comparisons, we focused on heat stress, as it elicited the most coherent and pronounced transcriptional and regulatory responses, enabling deeper network‐level interrogation. Co‐expression network analysis identified tissue‐specific modules enriched for phytohormone‐responsive genes, while gene regulatory network (GRN) mapping and cistrome analyses uncovered transcription factors (TFs) controlling key hub genes within these modules. Together, this study provides a foundational transcriptomic and network‐based resource for dissecting the regulatory architecture of abiotic stress responses in sorghum and offers prioritized candidates for future functional validation and engineering of climate‐resilient crops.

abiotic stress↗
DOE OSTI2024

Microbial inoculum effects on the rumen epithelial transcriptome and rumen epimural metatranscriptome in calves

Manipulation of the rumen microbial ecosystem in early life may affect ruminal fermentation and enhance the productive performance of dairy cows. The objective of this experiment was to evaluate the effects of dosing three different types of microbial inoculum on the rumen epithelium tissue (RE) transcriptome and the rumen epimural metatranscriptome (REM) in dairy calves. For this objective, 15 Holstein bull calves were enrolled in the study at birth and assigned to three different intraruminal inoculum treatments dosed orally once weekly from three to six weeks of age. The inoculum treatments were prepared from rumen contents collected from rumen fistulated lactating cows and were either autoclaved (control; ARF), processed by differential centrifugation to create the bacterial-enriched inoculum (BE), or through gravimetric separation to create the protozoal-enriched inoculum (PE). Calves were fed 2.5 L/d pasteurized waste milk 3x/d from 0 to 7 weeks of age and texturized starter until euthanasia at 9 weeks of age, when the RE tissues were collected for transcriptome and microbial metatranscriptome analyses, from four randomly selected calves from each treatment. The different types of inoculum altered the RE transcriptome and REM. Compared to ARF, 9 genes were upregulated in the RE of BE and 92 in PE, whereas between BE and PE there were 13 genes upregulated in BE and 114 in PE. Gene ontology analysis identified enriched GO terms in biological process category between PE and ARF, with no enrichment between BE and ARF. The RE functional signature showed different KEGG pathways related to BE and ARF, and no specific KEGG pathway for PE. We observed a lower alpha diversity index for RE microbiome in ARF (observed genera and Chao1 (p < 0.05)). Five microbial genera showed a significant correlation with the changes in host gene expression: Roseburia (25 genes), Entamoeba (two genes); Anaerosinus, Lachnospira, and Succiniclasticum were each related to one gene. sPLS-DA analysis showed that RE microbial communities differ among the treatments, although the taxonomic and functional microbial profiles show different distributions. Co-expression Differential Network Analysis indicated that both BE and PE had an impact on the abundance of KEGG modules related to acyl-CoA synthesis, type VI secretion, and methanogenesis, while PE had a significant impact on KEGGs related to ectoine biosynthesis and D-xylose transport. Our study indicated that artificial dosing with different microbial inocula in early life alters not only the RE transcriptome, but also affects the REM and its functions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Tripogon loliiformis tolerates rapid desiccation after metabolic and transcriptional priming during initial drying

Abstract Crop plants and undomesticated resilient species employ different strategies to regulate their energy resources and growth. Most crop species are sensitive to stress and prioritise rapid growth to maximise yield or biomass production. In contrast, resilient plants grow slowly, are small, and allocate their resources for survival in challenging environments. One small group of plants, termed resurrection plants, survive desiccation of their vegetative tissue and regain full metabolic activity upon watering. However, the precise molecular mechanisms underlying this extreme tolerance remain unknown. In this study, we employed a transcriptomics and metabolomics approach, to investigate the mechanisms of desiccation tolerance in Tripogon loliiformis , a modified desiccation-tolerant plant, that survives gradual but not rapid drying. We show that T. loliiformis can survive rapid desiccation if it is gradually dried to 60% relative water content (RWC). Furthermore, the gene expression data showed that T. loliiformis is genetically predisposed for desiccation in the hydrated state, as evidenced by the accumulation of MYB, NAC, bZIP, WRKY transcription factors along with the phytohormones, abscisic acid, salicylic acid, amino acids (e.g., proline) and TCA cycle sugars during initial drying. Through network analysis of co-expressed genes, we observed differential responses to desiccation between T. loliiformis shoots and roots. Dehydrating shoots displayed global transcriptional changes across broad functional categories, although no enrichment was observed during drying. In contrast, dehydrating roots showed distinct network changes with the most significant differences occurring at 40% RWC. The cumulative effects of the early stress responses may indicate the minimum requirements of desiccation tolerance and enable T. loliiformis to survive rapid drying. These findings potentially hold promise for identifying biotechnological solutions aimed at developing drought-tolerant crops without growth and yield penalties.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Protocol for applying a network-enabled gene discovery pipeline to non-model plant species

Identifying upstream regulators of key genes is essential for understanding gene regulatory mechanisms and translating these insights into functional targets. Here, we present a protocol for applying the network-enabled gene discovery pipeline (NEEDLE) to non-model plant species. We describe steps for environment setup, data preparation, computational analysis, expected outputs, and parameter considerations. NEEDLE integrates RNA sequencing (RNA-seq) processing, weighted gene co-expression analysis (WGCNA), Gene Network Inference with Ensemble of trees (GENIE3), and promoter conservation analysis to prioritize candidate transcriptional regulators.

Plant Sciences↗
DOE OSTI2023

Improved Biofuel Production through Discovery and Engineering of Terpene Metabolism in Switchgrass

Project Objectives - Of the myriad specialized metabolites that plants deploy to adapt to environmental challenges, terpenes form the largest group. In many major crops, unique terpene blends serve as key stress defenses that directly impact plant fitness and yield. In addition, terpenes, such as bisabolene and pinene, are used for producing renewable biofuels. Essential to advancing a broader use of terpenes for biofuel feedstock engineering is a system-wide knowledge of the diverse biosynthetic machinery and defensive potential of often species-specific terpene blends. The proposed project would merge genome-wide enzyme discovery with comparative –omics, protein structural and plant microbiome studies to define the biosynthesis and stress-defensive functions of the switchgrass (Panicum virgatum) terpene network. These insights would be combined with developing and applying non-transgenic genome editing tools to design plants with desirable terpene blends for higher productivity and biofuel production on marginal lands. As a dedicated lignocellulosic feedstock for U.S. biofuel production with high net energy yield, stress tolerance, and available genome resources, switchgrass is well-suited for devising new avenues for biofuel production. Project Description – The diversity of plant terpene defenses is governed by species-specific families of terpene synthase (TPS) and cytochrome P450 monooxygenase (P450) enzymes. Mining of the switchgrass genome (genotype Alamo) identified ~100 TPS and P450 candidate genes, and combinatorial biochemical analysis of synthesized TPSs and P450s revealed more than a dozen enzymes with common and novel activities. In addition, several identified terpene metabolites and the corresponding transcripts were up-regulated in response to abiotic stressors. These findings demonstrate a unique switchgrass terpene network with probable importance to abiotic stress tolerance, thus providing a large chemical portfolio for optimizing crop resistance, yield, and biofuel composition. Leveraging these preliminary data, we propose to generate a genome-wide map of the switchgrass terpene metabolic network through multi-gene co-expression analyses that allow the efficient cross-validation of TPS and P450 functions. Key enzymes would further be applied to structure-function studies via X-ray protein crystallography, homology modeling and site-directed mutagenesis to gain mechanistic insight into the catalytic specificity of switchgrass terpene metabolism and provide gene and amino acid targets for genome editing. In tandem with terpene pathway discovery, system-wide metabolomics, transcriptomics and proteomics studies in switchgrass accessions of contrasting drought tolerance would define the role of switchgrass terpene metabolism in conferring abiotic stress resilience. Metabolic changes would be assessed in a combined approach of targeted (terpenes) and untargeted metabolite profiling using a high-resolution LC-MS/MS approach, differential gene expression analyses through multiplexed Illumina RNA sequencing, and quantitative analysis of high-priority pathway enzymes using multiple reaction monitoring (MRM). Drawing on these insights, knock-down/out mutants of stress-associated pathway nodes would be generated by optimizing transient virus-induced gene silencing (VIGS) and CRISPR/Cas9 systems under control of the Tobacco Rattle Virus (TRV). The resulting mutant lines would then be analyzed for stress susceptibility and the impact on the root microbiome to define gene functions in planta. Knowledge of terpene pathways, enzyme mechanisms and bioactivities would be applied to enhance switchgrass stress resilience and to tailor-make terpene blends for biofuel production. Here, TRV-enabled CRISPR/Cas9 genome editing, including allele-specific knock-out of redundant genes, engineering of enzyme specificity via structure-guided point mutations, and overexpression of terpene genes relevant to stress-protection or biofuel production, would be used to increase metabolic flux toward desired pathways. Broader Impacts - Integrating the system-wide discovery, mechanistic analysis and non-transgenic genome engineering of the switchgrass terpene network aligns the required steps to unlock the chemical potential of this important metabolite class to generate crops that are more resistant to stress and provide advanced biofuel production in light of rising climate pressures as foreseeable challenges for bioenergy crop cultivation. The proposed project would further offer interdisciplinary student training through active involvement in the project and integration of research concepts and outcomes into newly-developed graduate and undergraduate courses on Plant Biotechnology.

09 BIOMASS FUELS↗
DOE OSTI2024

Evolution of the regulatory subunits for the heteromeric acetyl-CoA carboxylase

The committed step for de novo fatty acid (FA) synthesis is the ATP-dependent carboxylation of acetyl-coenzyme A catalysed by acetyl-CoA carboxylase (ACCase). In most plants, ACCase is a multi-subunit complex orthologous to prokaryotes. However, unlike prokaryotes, the plant and algal orthologues are comprised both catalytic and additional dedicated regulatory subunits. Novel regulatory subunits, biotin lipoyl attachment domain-containing proteins (BADC) and carboxyltransferase interactors (CTI) (both three-gene families inArabidopsis) represent new effectors specific to plants and certain algal species. The evolutionary history of these genes in autotrophic eukaryotes remains elusive, making it an ongoing area of research. Analyses of potential protein–protein and co-occurrence interactions, informed by gene network patterns using the STRING database, inArabidopsis thalianaandChlamydomonas reinhardtiiunveil intricate gene associations with ACCase, suggesting a complex interplay between FA synthesis and other cellular processes. Among both species, a higher number of co-expressed genes was identified inArabidopsis, indicating a wider potential regulatory network of ACCase in plants. This review investigates the extent to which these genes arose in autotrophic eukaryotes and provides insights into their evolutionary trajectory. This article is part of the theme issue ‘The evolution of plant metabolism’.

Life Sciences & Biomedicine - Other Topics↗
DOE OSTI2022

Exploring Camelina sativa lipid metabolism regulation by combining gene co-expression and DNA affinity purification analyses

Camelina (Camelina sativa) is an annual oilseed plant that is gaining momentum as a biofuel cover crop. Understanding gene regulatory networks is essential to deciphering plant metabolic pathways, including lipid metabolism. Furthermore, we take advantage of a growing collection of gene expression datasets to predict transcription factors (TFs) associated with the control of Camelina lipid metabolism. We identified approximately 350 TFs highly co-expressed with lipid-related genes (LRGs). These TFs are highly represented in the MYB, AP2/ERF, bZIP, and bHLH families, including a significant number of homologs of well-known Arabidopsis lipid and seed developmental regulators. After prioritizing the top 22 TFs for further validation, we identified DNA-binding sites and predicted target genes for 16 out of the 22 TFs tested using DNA affinity purification followed by sequencing (DAP-seq). Enrichment analyses of targets supported the co-expression prediction for most TF candidates, and the comparison to Arabidopsis revealed some common themes, but also aspects unique to Camelina. Within the top potential lipid regulators, we identified CsaMYB1, CsaABI3AVP1-2, CsaHB1, CsaNAC2, CsaMYB3, and CsaNAC1 as likely involved in the control of seed fatty acid elongation and CsaABI3AVP1-2 and CsabZIP1 as potential regulators of the synthesis and degradation of triacylglycerols (TAGs), respectively. Altogether, the integration of co-expression data and DNA-binding assays permitted us to generate a high-confidence and short list of Camelina TFs involved in the control of lipid metabolism during seed development.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Metabolome and RNA-seq Analysis of Responses to Nitrogen Deprivation and Resupply in Tea Plant (Camellia sinensis) Roots

Nitrogen (N) is an important contributor in regulating plant growth and development as well as secondary metabolites synthesis, so as to promote the formation of tea quality and flavor. Theanine, polyphenols, and caffeine are important secondary metabolites in tea plant. In this study, the responses of Camellia sinensis roots to N deprivation and resupply were investigated by metabolome and RNA-seq analysis. N deficiency induced content increase for most amino acids (AAs) and reduction for the remaining AAs, polyphenols, and caffeine. After N recovery, the decreased AAs and polyphenols showed a varying degree of recovery in content, but caffeine did not. Meanwhile, theanine increased in content, but its related synthetic genes were down-regulated, probably due to coordination of the whole N starvation regulatory network. Flavonoids-related pathways were relatively active following N stress according to KEGG enrichment analysis. Gene co-expression analysis revealed TCS2 , AMT1;1 , TAT2 , TS , and GOGAT as key genes, and TFs like MYB, bHLH, and NAC were also actively involved in N stress responses in C. sinensis roots. These findings facilitate the understanding of the molecular mechanism of N regulation in tea roots and provide genetic reference for improving N use efficiency in tea plant.

Xu, Wenluan↗
NASA NTRS2000

From Coexpression to Coregulation: An Approach to Inferring Transcriptional Regulation Among Gene Classes from Large-Scale Expression Data

We provide preliminary evidence that existing algorithms for inferring small-scale gene regulation networks from gene expression data can be adapted to large-scale gene expression data coming from hybridization microarrays. The essential steps are (I) clustering many genes by their expression time-course data into a minimal set of clusters of co-expressed genes, (2) theoretically modeling the various conditions under which the time-courses are measured using a continuous-time analog recurrent neural network for the cluster mean time-courses, (3) fitting such a regulatory model to the cluster mean time courses by simulated annealing with weight decay, and (4) analysing several such fits for commonalities in the circuit parameter sets including the connection matrices. This procedure can be used to assess the adequacy of existing and future gene expression time-course data sets for determining transcriptional regulatory relationships such as coregulation.

Mjolsness, Eric↗
DOE OSTI2023

Ectomycorrhizal symbiosis prepares its host locally and systemically for abiotic cue signaling

Tree growth and survival are dependent on their ability to perceive signals, integrate them, and trigger timely and fitted molecular and growth responses. While ectomycorrhizal symbiosis is a predominant tree-microbe interaction in forest ecosystems, little is known about how and to what extent it helps trees cope with environmental changes. We hypothesized that the presence of Laccaria bicolor influences abiotic cue perception by Populus trichocarpa and the ensuing signaling cascade. We submitted ectomycorrhizal or non-ectomycorrhizal P. trichocarpa cuttings to short-term cessation of watering or ozone fumigation to focus on signaling networks before the onset of any physiological damage. Poplar gene expression, metabolite levels, and hormone levels were measured in several organs (roots, leaves, mycorrhizas) and integrated into networks. We discriminated the signal responses modified or maintained by ectomycorrhization. Ectomycorrhizas buffered hormonal changes in response to short-term environmental variations systemically prepared the root system for further fungal colonization and alleviated part of the root abscisic acid (ABA) signaling. Furthermore, the presence of ectomycorrhizas in the roots also modified the leaf multi-omics landscape and ozone responses, most likely through rewiring of the molecular drivers of photosynthesis and the calcium signaling pathway. In conclusion, P. trichocarpa-L. bicolor symbiosis results in a systemic remodeling of the host's signaling networks in response to abiotic changes. In addition, ectomycorrhizal, hormonal, metabolic, and transcriptomic blueprints are maintained in response to abiotic cues, suggesting that ectomycorrhizas are less responsive than non-mycorrhizal roots to abiotic challenges.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Transcriptional networks underlying a primary ovarian insufficiency disorder in alligators naturally exposed to EDCs

Interactions between the endocrine system and environmental contaminants are responsible for impairing reproductive development and function. Despite the taxonomic diversity of affected species and attendant complexity inherent to natural systems, the underlying signaling pathways and cellular consequences are mostly studied in lab models. To resolve the genetic and endocrine pathways that mediate affected ovarian function in organisms exposed to endocrine disrupting contaminants in their natural environments, we assessed broad-scale transcriptional and steroidogenic responses to exogenous gonadotropin stimulation in juvenile alligators (Alligator missippiensis) originating from a lake with well-documented pollution (Lake Apopka, FL) and a nearby reference site (Lake Woodruff, FL). Here, we found that individuals from Lake Apopka are characterized by hyperandrogenism and display hyper-sensitive transcriptional responses to gonadotropin stimulation when compared to individuals from Lake Woodruff. Site-specific transcriptomic divergence appears to be driven by wholly distinct subsets of transcriptional regulators, indicating alterations to fundamental genetic pathways governing ovarian function. Consistent with broad-scale transcriptional differences, ovaries of Lake Apopka alligators displayed impediments to folliculogenesis, with larger germinal beds and decreased numbers of late-stage follicles. After resolving the ovarian transcriptome into clusters of co-expressed genes, most site-associated modules were correlated to ovarian follicule phenotypes across individuals. However, expression of two site-specific clusters were independent of ovarian cellular architecture and are hypothesized to represent alterations to cell-autonomous transcriptional programs. Collectively, our findings provide high resolution mapping of transcriptional patterns to specific reproductive function and advance our mechanistic understanding regarding impaired reproductive health in an established model of environmental endocrine disruption.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Chromosome-level genome assembly of Quercus variabilis provides insights into the molecular mechanism of cork thickness

Quercus variabilis is a deciduous woody species with high ecological and economic value and is a major source of cork in East Asia. Cork from thick softwood sheets have higher commercial value than those from thin sheets. It is extremely difficult to genetically improve Q. variabilis to produce high quality softwood due to the lack of genomic information. Here, we present a high-quality chromosomal genome assembly for Q. variabilis with length of 791,89 Mb and 54,606 predicted genes. Comparative analysis of protein sequences of Q. variabilis with 11 other species revealed that specific and expanded gene families were significantly enriched in the "fatty acid biosynthesis" pathway in Q. variabilis, which may contribute to the formation of its unique cork. Additionally, based on weighted correlation network analysis of time-course (i.e., five important developmental ages) gene expression data in thick-cork versus thin-cork genotypes of Q. variabilis, we identified one co-expression gene module associated with the thick-cork trait. Within this co-expression gene module, 10 hub genes were associated with suberin biosynthesis. Furthermore, we identified a total of 198 suberin biosynthesis-related new candidate genes that were up-regulated in trees with a thick cork layer relative to those with a thin cork layer. Also, we found that some genes related to cell expansion and cell division were highly expressed in trees with a thick cork layer. Collectively, our results revealed that two metabolic pathways (i.e., suberin biosynthesis, fatty acid biosynthesis), along with other genes involved in cell expansion, cell division, and transcriptional regulation, were associated with the thick-cork trait in Q. variabilis, providing insights into the molecular basis of cork development and knowledge for informing genetic improvement of cork thickness in Q. variabilis and closely related species.

59 BASIC BIOLOGICAL SCIENCES↗