Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “gene 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 19 records

Evaluating the performance of random forest and iterative random forest based methods when applied to gene expression data

Gene-to-gene networks, such as Gene Regulatory Networks (GRN) and Predictive Expression Networks (PEN) capture relationships between genes and are beneficial for use in downstream biological analyses. There exists multiple network inference tools to produce these gene-to-gene networks from matrices of gene expression data. Random Forest-Leave One Out Prediction (RF-LOOP) is a method that has been shown to be efficient at producing these gene-to-gene networks, frequently known as GEne Network Inference with Ensemble of trees (GENIE3). Random Forest can be replaced in this process by iterative Random Forest (iRF), which performs variable selection and boosting. Here we validate that iterative Random Forest-Leave One Out Prediction (iRF-LOOP) produces higher quality networks than GENIE3 (RF-LOOP). We use both synthetic and empirical networks from the Dialogue for Reverse Engineering Assessment and Methods (DREAM) Challenges by Sage Bionetworks, as well as two additional empirical networks created from Arabidopsis thaliana and Populus trichocarpa expression data.

59 BASIC BIOLOGICAL SCIENCES↗

Disease resistance in coral is mediated by distinct adaptive and plastic gene expression profiles

Infectious diseases are an increasing threat to coral reefs, resulting in altered community structure and hindering the functional contributions of disease-susceptible species. We exposed seven reef-building coral species from the Caribbean to white plague disease and determined processes involved in (i) lesion progression, (ii) within-species gene expression plasticity, and (iii) expression-level adaptation among species that lead to differences in disease risk. Gene expression networks enriched in immune genes and cytoskeletal arrangement processes were correlated to lesion progression rates. Whether or not a coral developed a lesion was mediated by plasticity in genes involved in extracellular matrix maintenance, autophagy, and apoptosis, while resistant coral species had constitutively higher expression of intracellular protein trafficking. This study offers insight into the process involved in lesion progression and within- and between-species dynamics that lead to differences in disease risk that is evident on current Caribbean reefs.

59 BASIC BIOLOGICAL SCIENCES↗

Gene co-expression network analysis in zebrafish reveals chemical class specific modules

Zebrafish is a popular animal model used for high-throughput screening of chemical hazards, however, investigations of transcriptomic mechanisms of toxicity are still needed. Here, our goal was to identify genes and biological pathways that Aryl Hydrocarbon Receptor 2 (AHR2) Activators and flame retardant chemicals (FRCs) alter in developing zebrafish. Taking advantage of a compendium of phenotypically-anchored RNA sequencing data collected from 48-h post fertilization (hpf) zebrafish, we inferred a co-expression network that grouped genes based on their transcriptional response. Genes responding to the FRCs and AHR2 Activators localized to distinct regions of the network, with FRCs inducing a broader response related to neurobehavior. AHR2 Activators centered in one region related to chemical stress responses. We also discovered several highly co-expressed genes in this module, including cyp1a, and we subsequently show that these genes are definitively within the AHR2 signaling pathway. Systematic removal of the two chemical types from the data, and analysis of network changes identified neurogenesis associated with FRCs, and regulation of vascular development associated with both chemical classes. We also identified highly connected genes responding specifically to each class that are potential biomarkers of exposure. Overall, we created the first zebrafish chemical-specific gene co-expression network illuminating how chemicals alter the transcriptome relative to each other. In addition to our conclusions regarding FRCs and AHR2 Activators, our network can be leveraged by other studies investigating chemical mechanisms of toxicity.

59 BASIC BIOLOGICAL SCIENCES↗

MicroRNA Expression Profiles in Cultured Human Fibroblasts in Space

Microgravity, or an altered gravity environment from the static 1g, has been shown to influence global gene expression patterns and protein levels in living organisms. However, it is unclear how these changes in gene and protein expressions are related to each other or are related to other factors regulating such changes. A different class of RNA, the small non-coding microRNA (miRNA), can have a broad effect on gene expression networks by mainly inhibiting the translation process. Previously, we investigated changes in the expression of miRNA and related genes under simulated microgravity conditions on the ground using the NASA invented bioreactor. In comparison to static 1 g, simulated microgravity altered a number of miRNAs in human lymphoblastoid cells. Pathway analysis with the altered miRNAs and RNA expressions revealed differential involvement of cell communication and catalytic activity, as well as immune response signaling and NGF activation of NF-kB pathways under simulated microgravity condition. The network analysis also identified several projected networks with c- Rel, ETS1 and Ubiquitin C as key factors. In a flight experiment on the International Space Station (ISS), we will investigate the effects of actual spaceflight on miRNA expressions in nondividing human fibroblast cells in mostly G1 phase of the cell cycle. A fibroblast is a type of cell that synthesizes the extracellular matrix and collagen, the structural framework for tissues, and plays a critical role in wound healing and other functions. In addition to miRNA expressions, we will investigate the effects of spaceflight on the cellular response to DNA damages from bleomycin treatment.

Wu, Honglu↗

Evaluating the Performance of Random Forest and Iterative Random Forest Based Methods when Applied to Gene Expression Data

Gene-to-gene networks, such as Gene Regulatory Networks (GRN) and Predictive Expression Networks (PEN) capture relationships between genes and are beneficial for use in downstream biological analyses. There exists multiple network inference tools to produce these gene-to-gene networks from matrices of gene expression data. Random Forest-Leave One Out Prediction (RF-LOOP) is a method that has been shown to be efficient at producing these gene-to-gene networks, frequently known as GEne Network Inference with Ensemble of trees (GENIE3). Here we validate that iterative Random Forest-Leave One Out Prediction (iRF-LOOP) produces higher quality networks than GENIE3. We use both synthetic and empirical networks from the Dialogue for Reverse Engineering Assessment and Methods (DREAM) Challenges by Sage Bionetworks, as well as two additional empirical networks created from Arabidopsis thaliana and Populus trichocarpa expression data.

iRF-Loop, expression network, Populus Trichocarpa↗

Constitutive and inducible oleoresin defenses share genetic architectures and mechanisms in Pinus taeda

The oleoresin defense system of loblolly pine (Pinus taeda) protects trees from insects and pathogens and is an important source of renewable biofuels and chemicals, but the genetic basis of oleoresin production is poorly understood. We characterized the genetic architecture of oleoresin flow, resin canal number, stem wood terpene content, and monoterpene composition in two clonal populations of P. taeda. We used quantitative genetic analyses, genome-wide association studies (GWASs), multiplex network learning, and gene expression profiling to elucidate shared gene networks underlying defense traits and to identify high-quality candidates for breeding and engineering loblolly pine. Genetic analyses revealed polygenic inheritance and trait-to-trait correlations provide strong evidence for shared genes regulating constitutive and induced oleoresin flow. We identified 236 single nucleotide polymorphisms associated with oleoresin flow, resin canal number, and terpene composition and highlight candidate genes likely involved in terpene biosynthesis, cambial meristem reprogramming, and pathogen perception and immune signaling. Fourteen GWAS candidates were methyl jasmonate-responsive in tissues where resin canals initiate and terpene production occurs. Integrating quantitative genetics, GWAS, gene expression, and multiplex network analyses enabled the prioritization of high-quality candidate genes. This work advances the development of more resilient loblolly pine optimized for ecological performance, renewable chemical, and biofuel production.

genome-wide association study↗

A Sorghum F-Box Protein Induces an Oxidative Burst in the Defense Against Colletotrichum sublineola

The hemibiotrophic fungal pathogen Colletotrichum sublineola is the causal agent of anthracnose in sorghum ( Sorghum bicolor), resulting in leaf blight, stalk rot, and head blight in susceptible genotypes, with yield losses of up to 50%. The development of anthracnose-resistant cultivars can reduce reliance on fungicides and provide a more sustainable and economical means for disease management. A previous genome-wide association study of the sorghum association panel identified the candidate resistance gene Sobic.005G172300 encoding an F-box protein. To better understand the role of this gene in the defense against C. sublineola, gene expression following infection with C. sublineola was monitored by RNA sequencing in seedlings of sorghum accession SC110, which harbored the resistance allele, and three accessions that harbored a susceptible allele. Only in SC110 did the expression of Sobic.005G172300 increase during the biotrophic phase of infection. Subsequent transcriptome analysis, gene co-expression networks, and gene regulatory networks of inoculated and mock-inoculated seedlings of resistant and susceptible accessions suggest that the increase in expression of Sobic.005G172300 induces an oxidative burst by lowering the concentration of ascorbic acid during the biotrophic phase of infection. Based on gene regulatory network analysis, the protein encoded by Sobic.005G172300 is proposed to target proteins involved in the biosynthesis of ascorbic acid for polyubiquitination through the SCF E3 ubiquitin ligase, causing their degradation via the proteasome.

Plant Sciences↗

Concentration-response gene expression analysis in zebrafish reveals phenotypically-anchored transcriptional responses to retene

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous environmental contaminants and are associated with human disease. Canonically, many PAHs induce toxicity via activation of the aryl hydrocarbon receptor (AHR) pathway. While the interaction between PAHs and the AHR is well-established, understanding which AHR-regulated transcriptional effects directly result in observable phenotypes and which are adaptive or benign is important to better understand PAH toxicity. Retene is a frequently detected PAH in environmental sampling and has been associated with AHR2-dependent developmental toxicity in zebrafish, though its mechanism of toxicity has not been fully elucidated. To interrogate transcriptional changes causally associated with retene toxicity, we conducted whole-animal RNA sequencing at 48 hours post-fertilization after exposure to eight retene concentrations. The concentrations were selected to produce effects ranging from no phenotype to mortality and malformations in 100% of animals at 5 days post-fertilization. We identified a concentration-response relationship between retene teratogenicity and differential gene expression in both number of DEGs and magnitude of expression change. Elevated expression of cyp1a at retene concentrations below the threshold for teratogenicity suggested that while cyp1a expression is a sensitive biomarker of AHR activation, it may be too sensitive to serve as a biomarker of AHR-dependent teratogenicity. Genes differentially expressed at only non-teratogenic concentrations were enriched for transforming growth factor-ß (TGF-ß) signaling pathway disruption while DEGs identified at only teratogenic concentrations were significantly enriched for response to xenobiotic stimulus and reduction-oxidation reaction activity. DEGs which spanned both non-teratogenic and teratogenic concentrations showed similar disrupted biological processes to those unique to teratogenic concentrations, indicating these processes were disrupted at low exposure concentrations. Gene co-expression network analysis identified several gene modules, including those associated with PAHs and AHR2 activation. One, Module 7, was strongly enriched for AHR2-associated genes and contained the strongest responses to retene. Benchmark concentration (BMC) of Module 7 genes identified a median BMC of 7.5 µM, nearly the highest retene concentration with no associated teratogenicity, supporting the hypothesis that Module 7 genes are largely responsible for retene toxicity.

Toxin, Zebrafish, Retene, Transcriptomics, network↗

Machine-learning from Pseudomonas putida KT2440 transcriptomes reveals its transcriptional regulatory network

Bacterial gene expression is orchestrated by numerous transcription factors (TFs). Elucidating how gene expression is regulated is fundamental to understanding bacterial physiology and engineering it for practical use. In this study, a machine-learning approach was applied to uncover the genome-scale transcriptional regulatory network (TRN) in Pseudomonas putida KT2440, an important organism for bioproduction. We performed independent component analysis of a compendium of 321 high-quality gene expression profiles, which were previously published or newly generated in this study. We identified 84 groups of independently modulated genes (iModulons) that explain 75.7% of the total variance in the compendium. With these iModulons, we (i) expand our understanding of the regulatory functions of 39 iModulon associated TFs (e.g., HexR, Zur) by systematic comparison with 1993 previously reported TF-gene interactions; (ii) outline transcriptional changes after the transition from the exponential growth to stationary phases; (iii) capture group of genes required for utilizing diverse carbon sources and increased stationary response with slower growth rates; (iv) unveil multiple evolutionary strategies of transcriptome reallocation to achieve fast growth rates; and (v) define an osmotic stimulon, which includes the Type VI secretion system, as coordination of multiple iModulon activity changes. Taken together, this study provides the first quantitative genome-scale TRN for P. putida KT2440 and a basis for a comprehensive understanding of its complex transcriptome changes in a variety of physiological states.

09 BIOMASS FUELS↗

Time‐series multi‐omics analysis of micronutrient stress in Sorghum bicolor reveals iron and zinc crosstalk and regulatory network conservation

Micronutrient stress impacts growth, biomass production, and grain yield in crops. Multi-omics studies are valuable resources in identifying genes for functional studies and trait improvement, such as accumulation of Fe or Zn under deficient or excess conditions for bioenergy or grain agriculture. We conducted transcriptomics and ionomics analyses on Sorghum bicolor BTx623, grown under Fe and Zn limited and excess conditions over a 21-day period. To identify early and late transcriptional response in roots and leaves, 180 RNAseq libraries were sequenced for differential expression and co-expression network analyses. Fe and Zn accumulation was measured using ICP-MS at each time point, and a fluorometer was used to estimate chlorophyll content in leaves. Among the four treatments, Fe limitation and Zn excess resulted in the largest phenotypic effects and transcriptional response in roots and leaves. Several of the reduction (Strategy I) and chelation (Strategy II) strategy genes that improve bioavailability of Fe and Zn in plant roots often used by non-grass and grass species, respectively, were differentially expressed. Gene regulatory network (GRN) analysis of roots revealed enrichment of genes from Fe limiting and Zn excess which strongly connect to homologues of SbFIT, SbPYE, and SbBTS as hub genes. The GRN for leaf responses showed homologues of SbPYE and SbBTS as hubs connecting genes for chloroplast biosynthesis, Fe-S cluster assembly, photosynthesis, and ROS scavenging. Expression analyses suggest sorghum uses Strategy II genes for Fe and Zn uptake, as expected, but can also utilize Strategy I genes, which may be advantageous in variable moisture environments. We found strong overlap between Fe and Zn responsive GRNs, indicative of micronutrient crosstalk. We also found conservation of root and leaf GRNs, and known homologous genes suggest strong constraints on homeostasis networks in plants. These data will provide a resource for functional genetics to enhance micronutrient transport in sorghum, and opportunities to conduct further comparative GRN analysis across diverse crops species.

59 BASIC BIOLOGICAL SCIENCES↗

Argonaute and Dicer are essential for communication between Trichoderma atroviride and fungal hosts during mycoparasitism

ABSTRACT Trichoderma species are known for their mycoparasitic activity against phytopathogenic fungi that cause significant economic losses in agriculture. During mycoparasitism, Trichoderma spp. recognize molecules produced by the host fungus and release secondary metabolites and hydrolytic enzymes to kill and degrade the host’s cell wall. Here, we explored the participation of the Trichoderma atroviride RNAi machinery in the interaction with six phytopathogenic fungi of economic importance. We determined that both Argonaute-3 and Dicer-2 play an essential role during mycoparasitism. Using an RNA-Seq approach, we identified that perception, detox, and cell wall degradation depend on the T. atroviride -RNAi when interacting with Alternaria alternata, Rhizoctonia solani AG2 , and R. solani AG5. Furthermore, we constructed a gene co-expression network that provides evidence of two gene modules regulated by RNAi, which play crucial roles in essential processes during mycoparasitism. In addition, based on small RNA-seq, we conclude that siRNAs regulate amino acid and carbon metabolism and communication during the Trichoderma -host interaction. Interestingly, our data suggest that siRNAs might regulate allorecognition ( het ) and transport genes in a cross-species manner. Thus, these results reveal a fine-tuned regulation in T. atroviride dependent on siRNAs that is essential during the biocontrol of phytopathogenic fungi, showing a greater complexity of this process than previously established. IMPORTANCE There is an increasing need for plant disease control without chemical pesticides to avoid environmental pollution and resistance, and the health risks associated with the application of pesticides are increasing. Employing Trichoderma species in agriculture to control fungal diseases is an alternative plant protection strategy that overcomes these issues without utilizing chemical fungicides. Therefore, understanding the biocontrol mechanisms used by Trichoderma species to antagonize other fungi is critical. Although there has been extensive research about the mechanisms involved in the mycoparasitic capability of Trichoderma species, there are still unsolved questions related to how Trichoderma regulates recognition, attack, and defense mechanisms during interaction with a fungal host. In this work, we report that the Argonaute and Dicer components of the RNAi machinery and the small RNAs they process are essential for gene regulation during mycoparasitism by Trichoderma atroviride .

59 BASIC BIOLOGICAL SCIENCES↗

Transcriptomic Network Analysis of Cyanobacterial-Methylotroph Interactions in Coculture and Axenic Conditions

A previous study demonstrated the potential for Cyanobacterial-Methylotroph cocultures to facilitate biogas processing as well as to be used in other biotechnological applications. To advance this technology, we investigated potential interactions between Cyanobacterium stanieri HL-69 (HL69) and Methylotuvimicrobium alkaliphilum 20Z (20Z) by inferring and analyzing gene co-expression networks under co-culture and axenic conditions. Five different co-expression networks were examined. These networks were inferred using gene expression profiles for 20Z axenic condition, HL-69 axenic, HL-69, 20Z coculture, HL-69 coculture, and cross-species HL-69-20Z coculture. Through the analysis of node (gene) betweenness and node normalized degree values in all five network cases, we compared adjustments in gene expression between growth conditions (axenic vs co-culture) as well as identify biological functions relevant to interspecies interactions. This analysis was done to distinguish between gene interactions within an organism and gene interactions between two organisms. Moreover, for all five cases we investigated two different network cutoff levels of 3,000 and 10,000. By shedding light on inter- and intra- species interactions, we hope to gain a better understanding of how these two organisms interact. This research will allow the investigation of further biotechnological applications of coculture systems and optimization of such applications for biotechnological purposes.

59 BASIC BIOLOGICAL SCIENCES↗

Metabolic and transcriptomic study of pennycress natural variation identifies targets for oil improvement

Pennycress (Thlaspi arvense L.), a member of the Brassicaceae family, produces seed oil high in erucic acid, suitable for biodiesel and aviation fuel. Although pennycress, a winter annual, could be grown as a dedicated bioenergy crop, an increase in its seed oil content is required to improve its economic competitiveness. The success of crop improvement relies upon finding the right combination of biomarkers and targets, and the best genetic engineering and/or breeding strategies. In this work, we combined biomass composition with metabolomic and transcriptomic studies of developing embryos from 22 pennycress natural variants to identify targets for oil improvement. The selected accession collection presented diverse levels of fatty acids at maturity ranging from 29% to 41%. Pearson correlation analyses, weighted gene co‐expression network analysis and biomarker identifications were used as complementary approaches to detect associations between metabolite level or gene expression and oil content at maturity. The results indicated that improving seed oil content can lead to a concomitant increase in the proportion of erucic acid without affecting the weight of embryos. Processes, such as carbon partitioning towards the chloroplast, lipid metabolism, photosynthesis, and a tight control of nitrogen availability, were found to be key for oil improvement in pennycress. Besides identifying specific targets, our results also provide guidance regarding the best timing for their modification, early or middle maturation. Thus, this work lays out promising strategies, specific for pennycress, to accelerate the successful development of lines with increased seed oil content for biofuel applications.

59 BASIC BIOLOGICAL SCIENCES↗

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↗

A Path-Based Analysis of Infected Cell Line and COVID-19 Patient Transcriptome Reveals Novel Potential Targets and Drugs Against SARS-CoV-2

Most transcriptomic studies of SARS-CoV-2 infection have focused on differentially expressed genes, which do not necessarily reveal the genes mediating the transcriptomic changes. In contrast, exploiting curated biological network, our PathExt tool identifies central genes from the differentially active paths mediating global transcriptomic response. Here we apply PathExt to multiple cell line infection models of SARS-CoV-2 and other viruses, as well as to COVID-19 patient-derived PBMCs. The central genes mediating SARS-CoV-2 response in cell lines were uniquely enriched for ATP metabolic process, G1/S transition, leukocyte activation and migration. In contrast, PBMC response reveals dysregulated cell-cycle processes. In PBMC, the most frequently central genes are associated with COVID-19 severity. Importantly, relative to differential genes, PathExt-identified genes show greater concordance with several benchmark anti-COVID-19 target gene sets. We propose six novel anti-SARS-CoV-2 targets ADCY2, ADSL, OCRL, TIAM1, PBK, and BUB1, and potential drugs targeting these genes, such as Bemcentinib, Phthalocyanine, and Conivaptan.

59 BASIC BIOLOGICAL SCIENCES↗

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↗

Generating Co-expression Networks for Three Cyanobacteria: Synechococcus sp. PCC 7942, Synechococcus sp. PCC 7002, Synechocystis sp. PCC 6803

Cyanobacteria are photosynthetic organisms capable of high growth rate and represent a promising bioplatform for harnessing the sun’s energy to make biofuel. Additionally, the process of photosynthesis absorbs CO2 from the environment. Understanding the metabolic processes involved in photosynthesis could lead to solutions to the recent rise of CO2 concentration in Earth’s atmosphere and the associated climate change. More research on the transcriptional regulation of these cells is needed to learn how to harness the untapped potential of cyanobacteria for these applications. Transcriptional analysis via RNA-seq provides an understanding of how gene expression changes at the mRNA level under diverse growing conditions. I systematically collected and analyzed RNA-Seq data obtained under a variety of conditions and available on the NCBI database for three cyanobacteria model organisms: Synechococcus elongatus sp. PCC 7942, Synechococcus sp. PCC 7002, and Synechocystis sp. PCC 6803. For each organism, the data was mapped to a reference genome to characterize the RNA expression profile. Samples were checked for quality based on the number of reads and the correlation of the expression profile between labeled replicates. All samples were transformed into transcripts per million reads, followed by a log transformation to account for the wide range of sample sizes. Gene co-expression networks were generated and analyzed for each species using Cytoscape. These networks provide a base level of gene expression for each species. The network topology and high-betweeness nodes of these networks need to be analyzed further to provide insight on potential ways to harness cyanobacteria genetics. Additionally, these datasets can be used together to form a core genome network analysis- one that includes only the genes that are homologous between the three species. This project has prepared the way for a more in-depth study on photosynthetic microbes on a genetic level.

59 BASIC BIOLOGICAL SCIENCES↗