Search NASA⌕ Search

SEARCH · Search NASA

Results for “pathway enrichment 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.

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

At least 19 records

Interpreting omics data with pathway enrichment analysis

Pathway enrichment analysis is indispensable for interpreting omics datasets and generating hypotheses. However, the foundations of enrichment analysis remain elusive to many biologists. Here, in this study, we discuss best practices in interpreting different types of omics data using pathway enrichment analysis and highlight the importance of considering intrinsic features of various types of omics data. We further explain major components that influence the outcomes of a pathway enrichment analysis, including defining background sets and choosing reference annotation databases. To improve reproducibility, we describe how to standardize reporting methodological details in publications. This article aims to serve as a primer for biologists to leverage the wealth of omics resources and motivate bioinformatics tool developers to enhance the power of pathway enrichment analysis.

60 APPLIED LIFE SCIENCES↗

RNA Splicing Events in Circulation Distinguish Individuals With and Without New-onset Type 1 Diabetes

Context: Alterations in RNA splicing may influence protein isoform diversity that contributes to or reflects the pathophysiology of certain diseases. Whereas specific RNA splicing events in pancreatic islets have been investigated in models of inflammation in vitro, how RNA splicing in the circulation correlates with or is reflective of type 1 diabetes (T1D) disease pathophysiology in humans remains unexplored. Objective: To use machine learning to investigate if alternative RNA splicing events differ between individuals with and without new-onset T1D and to determine if these splicing events provide insight into T1D pathophysiology. Methods: RNA deep sequencing was performed on whole blood samples from 2 independent cohorts: a training cohort consisting of 12 individuals with new-onset T1D and 12 age- and sex-matched nondiabetic controls and a validation cohort of the same size and demographics. Machine learning analysis was used to identify specific isoforms that could distinguish individuals with T1D from controls. Results: Distinct patterns of RNA splicing differentiated participants with T1D from unaffected controls. Notably, certain splicing events, particularly involving retained introns, showed significant association with T1D. Machine learning analysis using these splicing events as features from the training cohort demonstrated high accuracy in distinguishing between T1D subjects and controls in the validation cohort. Gene Ontology pathway enrichment analysis of the retained intron category showed evidence for a systemic viral response in T1D subjects. Conclusion: Alternative RNA splicing events in whole blood are significantly enriched in individuals with new-onset T1D and can effectively distinguish these individuals from unaffected controls. Further, our findings also suggest that RNA splicing profiles offer the potential to provide insights into disease pathogenesis.

60 APPLIED LIFE SCIENCES↗

3′ RNA-seq is superior to standard RNA-seq in cases of sparse data but inferior at identifying toxicity pathways in a model organism

The application of RNA-sequencing has led to numerous breakthroughs related to investigating gene expression levels in complex biological systems. Among these are knowledge of how organisms, such as the vertebrate model organism zebrafish (Danio rerio), respond to toxicant exposure. Recently, the development of 3' RNA-seq has allowed for the determination of gene expression levels with a fraction of the required reads compared to standard RNA-seq. While 3' RNA-seq has many advantages, a comparison to standard RNA-seq has not been performed in the context of whole organism toxicity and sparse data. Here, we examined samples from zebrafish exposed to perfluorobutane sulfonamide (FBSA) with either 3' or standard RNA-seq to determine the advantages of each with regards to the identification of functionally enriched pathways. We found that 3' and standard RNA-seq showed specific advantages when focusing on annotated or unannotated regions of the genome. We also found that standard RNA-seq identified more differentially expressed genes (DEGs), but that this advantage disappeared under conditions of sparse data. We also found that standard RNA-seq had a significant advantage in identifying functionally enriched pathways via analysis of DEG lists but that this advantage was minimal when identifying pathways via gene set enrichment analysis of all genes. These results show that each approach has experimental conditions where they may be advantageous. Our observations can help guide others in the choice of 3' RNA-seq vs standard RNA sequencing to query gene expression levels in a range of biological systems.

3’ RNA-seq↗

High School Citizen Scientists Use AI/ML to Predict Intra-Ocular Pressure From Gene Expression Data for Spaceflown Mice

Artificial Intelligence (AI) and Machine Learning (ML) have increasingly become pivotal in biological and biomedical research, largely due to the culture of open data sharing and its associated benefits. The methodologies inherent in AI/ML are particularly adept at identifying and forecasting biological phenotypes from the vast amounts of data generated by next-generation sequencing technologies. These techniques offer substantial promise for advancing research in space biosciences and for the development of automated systems for monitoring space health. Nevertheless, there are crucial aspects to consider when training, validating, and testing machine learning models in both biological research and clinical contexts. It is essential that Open Science principles, including data sharing and the availability of open-source code, are complemented by high-quality, publicly accessible training resources. These resources should focus on best practices and include modules based on real-world scientific cases and data to ensure that future AI/ML practitioners gain practical experience with genuine problems. Addressing this knowledge gap, we have designed, developed, and delivered both interactive and self-paced training programs for citizen scientists worldwide, enabling them to utilize AI/ML for space biology research. This initiative was made possible through generous funding from a Transformation to Open Science Training grant. The interactive training sessions, conducted this summer, utilized AI/ML techniques to analyze data from the Open Science Data Repository, specifically targeting the effects of spaceflight on ocular structure and function. The dataset OSD-583, from the Rodent Research 9 mission, provides experimental data detailing the ocular responses of mice subjected to a 35-day spaceflight, compared with ground control counterparts. Using OSD-583 as observational data, our summer training participants applied AI/ML methods to predict intraocular pressure from RNA-seq data and identify the genes most predictive of the observed responses. Further analysis through pathway enrichment and gene set enrichment revealed that these genes are involved in molecular and cellular processes contributing to retinal degeneration.

James Casaletto↗

Transcriptomic Analysis of Irradiated Mouse Retina Following Readaptation

Rodent models are used as analogs for studying the effects of spaceflight. NASA GeneLab provides access to relevant omics datasets generated from spaceflight and ground-based experiments allowing for additional retrospective analysis. In this study, we used GeneLab’s GLDS-203, a dataset generated by researchers at Loma Linda University to study the impact of prolonged unloading and/or low-dose radiation on mouse retina. We analyzed transcriptomics data from retina of mice irradiated with gamma-rays for 21 days followed by 7 days, 1 month, or 4 months of readaptation. We obtained raw gene counts from GeneLab and performed differential gene expression analysis after data normalization. For each of the three timepoints, we performed differential expression analysis to compare transcriptional profiles for retina from irradiated vs. non-irradiated (controls) mice, all exposed to gravity. We observed the highest number of differentially expressed genes at 7 days, followed by 1 month and 4 months. Enrichment analysis showed top pathways (adjusted p-value < 0.05) were related to transport along microtubule and photoreceptor cell development in the 7-day readaptation group. Fewer significantly enriched pathways were observed for the 1-month group and included mRNA metabolic processes and neuron differentiation. No significantly enriched pathways were found in the 4-month group. The Gene Ontology biological processes common between the 7 days and 1-month groups include visual perception, synapse organization, and perception of light stimulus. This analysis is part of a larger effort to characterize the molecular mechanisms involved in retinal readaptation following radiation exposure. Future analyses will include other related retina datasets in GeneLab repository to assess whether gene expression patterns are consistent across different study cohorts.

Prachi Kothiyal↗

Proteogenomic characterization of difficult-to-treat breast cancer with tumor cells enriched through laser microdissection

Abstract Background Breast cancer (BC) is the most commonly diagnosed cancer and the leading cause of cancer death among women globally. Despite advances, there is considerable variation in clinical outcomes for patients with non-luminal A tumors, classified as difficult-to-treat breast cancers (DTBC). This study aims to delineate the proteogenomic landscape of DTBC tumors compared to luminal A (LumA) tumors. Methods We retrospectively collected a total of 117 untreated primary breast tumor specimens, focusing on DTBC subtypes. Breast tumors were processed by laser microdissection (LMD) to enrich tumor cells. DNA, RNA, and protein were simultaneously extracted from each tumor preparation, followed by whole genome sequencing, paired-end RNA sequencing, global proteomics and phosphoproteomics. Differential feature analysis, pathway analysis and survival analysis were performed to better understand DTBC and investigate biomarkers. Results We observed distinct variations in gene mutations, structural variations, and chromosomal alterations between DTBC and LumA breast tumors. DTBC tumors predominantly had more mutations inTP53,PLXNB3, Zinc finger genes, and fewer mutations inSDC2,CDH1,PIK3CA,SVIL, andPTEN. Notably, Cytoband 1q21, which contains numerous cell proliferation-related genes, was significantly amplified in the DTBC tumors. LMD successfully minimized stromal components and increased RNA–protein concordance, as evidenced by stromal score comparisons and proteomic analysis. Distinct DTBC and LumA-enriched clusters were observed by proteomic and phosphoproteomic clustering analysis, some with survival differences. Phosphoproteomics identified two distinct phosphoproteomic profiles for high relapse-risk and low relapse-risk basal-like tumors, involving several genes known to be associated with breast cancer oncogenesis and progression, includingKIAA1522,DCK,FOXO3,MYO9B,ARID1A,EPRS,ZC3HAV1, andRBM14. Lastly, an integrated pathway analysis of multi-omics data highlighted a robust enrichment of proliferation pathways in DTBC tumors. Conclusions This study provides an integrated proteogenomic characterization of DTBC vs LumA with tumor cells enriched through laser microdissection. We identified many common features of DTBC tumors and the phosphopeptides that could serve as potential biomarkers for high/low relapse-risk basal-like BC and possibly guide treatment selections.

Oncology↗

Radiation Quality Effects on Transcriptome Profiles in 3-d Cultures After Particle Irradiation

In this work, we evaluate the differential effects of low- and high-LET radiation on 3-D organotypic cultures in order to investigate radiation quality impacts on gene expression and cellular responses. Reducing uncertainties in current risk models requires new knowledge on the fundamental differences in biological responses (the so-called radiation quality effects) triggered by heavy ion particle radiation versus low-LET radiation associated with Earth-based exposures. We are utilizing novel 3-D organotypic human tissue models that provide a format for study of human cells within a realistic tissue framework, thereby bridging the gap between 2-D monolayer culture and animal models for risk extrapolation to humans. To identify biological pathway signatures unique to heavy ion particle exposure, functional gene set enrichment analysis (GSEA) was used with whole transcriptome profiling. GSEA has been used extensively as a method to garner biological information in a variety of model systems but has not been commonly used to analyze radiation effects. It is a powerful approach for assessing the functional significance of radiation quality-dependent changes from datasets where the changes are subtle but broad, and where single gene based analysis using rankings of fold-change may not reveal important biological information. We identified 45 statistically significant gene sets at 0.05 q-value cutoff, including 14 gene sets common to gamma and titanium irradiation, 19 gene sets specific to gamma irradiation, and 12 titanium-specific gene sets. Common gene sets largely align with DNA damage, cell cycle, early immune response, and inflammatory cytokine pathway activation. The top gene set enriched for the gamma- and titanium-irradiated samples involved KRAS pathway activation and genes activated in TNF-treated cells, respectively. Another difference noted for the high-LET samples was an apparent enrichment in gene sets involved in cycle cycle/mitotic control. It is plausible that the enrichment in these particular pathways results from the complex DNA damage resulting from high-LET exposure where repair processes are not completed during the same time scale as the less complex damage resulting from low-LET radiation.

Patel, Z. S.↗

Multi-Omics Analysis of Mouse Retina Following Low Dose Radiation and/or Hindlimb Unloading

Rodent models have been used as analogs for studying the effects of spaceflight. NASA’s GeneLab provides access to omics datasets generated from spaceflight and ground-based experiments allowing for additional retrospective analysis. We used GeneLab’s GLDS-203, a dataset generated by researchers at Loma Linda University to study the impact of prolonged unloading and/or low-dose radiation on mouse retina. The purpose of this study was to understand the effect of gamma radiation and/or hindlimb unloading on mice retinas through a multi-omics analysis. In the experiment that generated the omics data, mice were irradiated with gamma-ray and/or subjected to hindlimb unloading for 21 days and multi-omics analysis was performed at 7 days, 1 month, or 4 months post exposure. In the current study, for each of the three timepoints, we compared epigenomic profiles for retinas from exposed mice against timepoint-matched controls. We identified a total of 5,271 differentially methylated loci (DML) and 321 differentially methylated regions (DMR; using a sliding window and step size of 500 bp) with methylation difference > 10% and q-value < 0.05 (sliding linear model corrected p-value) across the nine exposure groups. Highest correlation in methylation difference was seen for significant DMLs (q-value < 0.05) across different conditions at same post exposure timepoint (Figure 1).The location of DMLs and DMRs were characterized with respect to CpG islands and shores, putative promoters, gene body, and intergenic regions (Table 1). We analyzed RNA-seq counts and performed gene set enrichment analysis using differential expression results from comparing each exposure group to its timepoint-matched control group. Significant pathways (adjusted p-value <0.05) enriched in all three microgravity-only groups were related to morphogenesis of a branching epithelium, skeletal muscle cell differentiation, and response to fibroblast growth factor. Common processes across all timepoints in the radiation-only groups were retina homeostasis, synaptic vesicle exocytosis-endocytosis, and chemotaxis. In the combination groups, regulation of trans-synaptic signaling, and Rho protein signal transduction were enriched at all three timepoints. Processes related to purine nucleotide metabolism were enriched in all nine exposure groups, with activation at 1 month, and suppression at 7 days and 4 months. A total of 14 genes contained at least one DML and were differentially expressed at adjusted p-value < 0.05, including genes implicated in cataract development (Sipa1l3, Crybb3) and those involved in cytoskeletal organization (Plec, Flnb, Eef1a1). This analysis is part of a larger effort to understand the molecular mechanisms following spaceflight exposures that can help translate effects observed in animal models to human impacts.

Prachi Kothiyal↗

Radiation Quality Effects on Transcriptome Profiles in 3-D Cultures After Charged Particle Irradiation

In this work, we evaluated the differential effects of low- and high-LET radiation on 3-D organotypic cultures in order to investigate radiation quality impacts on gene expression and cellular responses. Current risk models for assessment of space radiation-induced cancer have large uncertainties because the models for adverse health effects following radiation exposure are founded on epidemiological analyses of human populations exposed to low-LET radiation. Reducing these uncertainties requires new knowledge on the fundamental differences in biological responses (the so-called radiation quality effects) triggered by heavy ion particle radiation versus low-LET radiation associated with Earth-based exposures. In order to better quantify these radiation quality effects in biological systems, we are utilizing novel 3-D organotypic human tissue models for space radiation research. These models hold promise for risk assessment as they provide a format for study of human cells within a realistic tissue framework, thereby bridging the gap between 2-D monolayer culture and animal models for risk extrapolation to humans. To identify biological pathway signatures unique to heavy ion particle exposure, functional gene set enrichment analysis (GSEA) was used with whole transcriptome profiling. GSEA has been used extensively as a method to garner biological information in a variety of model systems but has not been commonly used to analyze radiation effects. It is a powerful approach for assessing the functional significance of radiation quality-dependent changes from datasets where the changes are subtle but broad, and where single gene based analysis using rankings of fold-change may not reveal important biological information.

Patel, Zarana S.↗

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↗

Integrated transcriptomic and metabolomic investigation of the genes and metabolites involved in swine follicular cyst formation

Follicular cysts are a common reproductive disorder in mammals that is usually caused by stress. However, the pathogenesis of follicular cysts in sows remains unclear. To provide new insights into the mechanisms of follicular cyst formation in pigs, we conducted a combined transcriptomic and metabolomic analysis on theca interna and mural granulosa cells of follicular cysts and mature follicles. We identified 2,533 up-regulated and 1,355 down-regulated genes in follicular cysts, compared with mature follicles. These differentially expressed genes were mainly found in signaling pathways related to tumor formation and cortisol synthesis and secretion as shown by Ingenuity Pathway Analysis, which predicted 4,362 upstream regulatory factors. The combined gene expression and pathway analysis identified the following genes as potential biomarkers for porcine follicular cysts: cytochrome P450 family 2 subfamily C polypeptide 18, L-lactate dehydrogenase, carbamoyl-phosphate synthase, fibroblast growth factor 7, integrin binding sialoprotein, interleukin 23 receptor, prolactin receptor, epiregulin, interleukin 1 receptor type II, arginine vasopressin receptor 1A, fibroblast growth factor 10, claudin 7, G Protein Subunit Gamma 3, cholecystokinin B receptor and cytosolic phospholipase A2 . Metabolomics analysis found significant differences in 87 metabolites, which were enriched in unsaturated fatty acid biosynthesis, and sphingolipid signaling pathways. These results provide valuable information on the molecular mechanisms of follicular cyst formation, which may facilitate the development of new therapeutics to prevent and treat follicular cysts.

Dai, Jiage↗

Substrate Matters: Ionic Silver Alters Lettuce Growth, Nutrient Uptake, and Root Microbiome in a Hydroponics System

Ionic silver (Ag+) is being investigated as a residual biocide for use in NASA spacecraft potable water systems on future crewed missions. This water will be used to irrigate future spaceflight crop production systems. We have evaluated the impact of three concentrations (31 ppb, 125 ppb, and 500 ppb) of ionic silver biocide solutions on lettuce in an arcillite (calcinated clay particle substrate) and hydroponic (substrate-less) growth setup after 28 days. Lettuce plant growth was reduced in the hydroponic samples treated with 31 ppb silver and severely stunted for samples treated at 125 ppb and 500 ppb silver. No growth defects were observed in arcillite-grown lettuce. Silver was detectable in the hydroponic-grown lettuce leaves at each concentration but was not detected in the arcillite-grown lettuce leaves. Specifically, when 125 ppb silver water was applied to a hydroponics tray, Ag+ was detected at an average amount of 7 μg/g (dry weight) in lettuce leaves. The increase in Ag+ corresponded with a decrease in several essential elements in the lettuce tissue (Ca, K, P, S). In the arcillite growth setup, silver did not impact the plant root zone microbiome in terms of alpha diversity and relative abundance between treatments and control. However, with increasing silver concentration, the alpha diversity increased in lettuce root samples and in the water from the hydroponics tray samples. The genera in the hydroponic root and water samples were similar across the silver concentrations but displayed different relative abundances. This suggests that ionic silver was acting as a selective pressure for the microbes that colonize the hydroponic water. The surviving microbes likely utilized exudates from the stunted plant roots as a carbon source. Analysis of the root-associated microbiomes in response to silver showed enrichment of metagenomic pathways associated with alternate carbon source utilization, fatty-acid synthesis, and the ppGpp (guanosine 3′-diphosphate 5′-diphosphate) stringent response global regulatory system that operates under conditions of environmental stress. Nutrient solutions containing Ag+ in concentrations greater than 31 ppb in hydroponic systems lacking cation-exchange capacity can severely impact crop production due to stunting of plant growth.

lettuce↗

An emerging multi-omic understanding of the genetics of opioid addiction

Opioid misuse, addiction, and associated overdose deaths remain global public health crises. Despite the tremendous need for pharmacological treatments, current options are limited in number, use, and effectiveness. Fundamental leaps forward in our understanding of the biology driving opioid addiction are needed to guide development of more effective medication-assisted therapies. This Review focuses on the omics-identified biological features associated with opioid addiction. Recent GWAS have begun to identify robust genetic associations, including variants in OPRM1, FURIN, and the gene cluster SCAI/PPP6C/RABEPK. An increasing number of omics studies of postmortem human brain tissue examining biological features (e.g., histone modification and gene expression) across different brain regions have identified broad gene dysregulation associated with overdose death among opioid misusers. Drawn together by meta-analysis and multi-omic systems biology, and informed by model organism studies, key biological pathways enriched for opioid addiction–associated genes are emerging, which include specific receptors (e.g., GABAB receptors, GPCR, and Trk) linked to signaling pathways (e.g., Trk, ERK/MAPK, orexin) that are associated with synaptic plasticity and neuronal signaling. Studies leveraging the agnostic discovery power of omics and placing it within the context of functional neurobiology will propel us toward much-needed, field-changing breakthroughs, including identification of actionable targets for drug development to treat this devastating brain disease.

60 APPLIED LIFE SCIENCES↗

Seed coat transcriptomic profiling of 5-593, a genotype important for genetic studies of seed coat color and patterning in common bean ( Phaseolus vulgaris L.)

Common bean (Phaseolus vulgaris L.) market classes have distinct seed coat colors, which are directly related to the diverse flavonoids found in the mature seed coat. To understand and elucidate the molecular mechanisms underlying the regulation of seed coat color, RNA-Seq data was collected from the black bean 5-593 and used for a differential gene expression and enrichment analysis from four different seed coat color development stages. 5-593 carries dominant alleles for 10 of the 11 major genes that control seed coat color and expression and has historically been used to develop introgression lines used for seed coat genetic analysis. Pairwise comparison among the four stages identified 6,294 differentially expressed genes (DEGs) varying from 508 to 5,780 DEGs depending on the compared stages. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that phenylpropanoid biosynthesis, flavonoid biosynthesis, and plant hormone signal transduction comprised the principal pathways expressed during bean seed coat pigment development. Transcriptome analysis suggested that most structural genes for flavonoid biosynthesis and some potential regulatory genes were significantly differentially expressed. Further studies detected 29 DEGs as important candidate genes governing the key enzymatic flavonoid biosynthetic pathways for common bean seed coat color development. Additionally, four gene models, Pv5-593.02G016100, 593.02G078700, Pv5-593.02G090900, and Pv5-593.06G121300, encode MYB-like transcription factor family protein were identified as strong candidate regulatory genes in anthocyanin biosynthesis which could regulate the expression levels of some important structural genes in flavonoid biosynthesis pathway. These findings provide a framework to draw new insights into the molecular networks underlying common bean seed coat pigment development.

60 APPLIED LIFE SCIENCES↗

Fine scale sampling reveals early differentiation of rhizosphere microbiome from bulk soil in young Brachypodium plant roots

Abstract For a deeper and comprehensive understanding of the composition and function of rhizosphere microbiomes, we need to focus at the scale of individual roots in standardized growth containers. Root exudation patterns are known to vary along distinct parts of the root even in juvenile plants giving rise to spatially distinct microbial niches. To address this, we analyzed the microbial community from two spatially distinct zones of the developing primary root (tip and base) in young Brachypodium distachyon grown in natural soil using standardized fabricated ecosystems known as EcoFABs as well as in more conventional pot and tubes. 16S rRNA based community analysis showed a strong rhizosphere effect resulting in significant enrichment of several OTUs belonging to Actinobacteria, Bacteroidetes, Firmicutes and Proteobacteria. However, microbial community composition did not differ between root tips and root base or across different growth containers. Functional analysis of bulk metagenomics revealed significant differences between root tips and bulk soil. The genes associated with different metabolic pathways and root colonization were enriched in root tips. On the other hand, genes associated with nutrient-limitation and environmental stress were prominent in the bulk soil compared to root tips, implying the absence of easily available, labile carbon and nutrients in bulk soil relative to roots. Such insights into the relationships between developing root and microbial communities are critical for judicious understanding of plant-microbe interactions in early developmental stages of plants.

Acharya, Shwetha M.↗

Cross Kingdom Analysis of Data Within the GeneLab Repository Identifies a Potential Conserved Response of Life to the Stress Associated with Spaceflight

It is important to determine the health risks and potential survival for astronauts associated with long-term space missions. This entails not only understanding the impact the space environment will have on humans, but also how it will affect other organisms needed for humans to survive in space such as plants. In addition, it has been reported in the literature that hundreds of genes seem to be conserved and/or transferred between different organisms from bacteria, archaea, fungi, microorganisms, and plants to animals. Since space travel involves humans in a closed environment over a long period of time, we hypothesize that potential conserved biological factors will occur between the different organisms in that environment possibly due to transfer of genes. Determining the conserved factors that are commonly being regulated in space can shed insight into possible universal master regulators and also determine the symbiotic relationship between the organisms in space. Utilizing NASA's GeneLab Data Repository (a rapidly expanding, curated clustering of spaceflight-related ‘omics-level datasets for all organisms), we were able to uncover a novel pathway and factors that were commonly shared between humans, mice, plants, C. Elegans, and drosophilas. Through ChIP-Seq enrichment analysis techniques utilizing various GeneLab datasets from each species that were flown in space, we found the following factors to be conserved across all species: oxidative stress, DNA damage (through GABPA/NRFs and NFY), SIX5, GTF2B and glutamine synthetase. Such commonalities would likely reflect the effects of factors such as microgravity and the increased radiation exposure inherent in spaceflight on basic physical processes shared by all biological systems at the cellular level. Differences between organismal responses revealed by GeneLab's data should also help understand the unique reactions to life in space that arise from the very different lifestyles of microbes, animals and plants.

Barker, Richard↗

Enhanced Resistance Pines for Improved Renewable Biofuel and Chemical Production (Technical Report)

We completed phenotyping constitutive and inducible oleoresin flow across two seasons, constitutive resin canal number and density and wood terpene content in our ADEPT2 and CCLONES populations. We completed genetic association between 19 oleoresin phenotypes and a total of 523,192 SNP markers from ADEPT2 and 13,883 SNP markers in CCLONES using four mixed linear models. A total of 293 significant SNPs (FDR = 0.20) were identified. We used the MENTOR tool to mine mechanistic connections from a multiplex network constructed from poplar multi-omic data to construct a conceptual model for a subset of these significant SNPs. Our model contains 6 transcriptional regulators in addition to 3 monoterpene synthases. To generate more lines of evidence for these significant SNPs, we completed a time course RNAseq experiment after inducing vascular zone cells to differentiate into new resin canals with a methyl jasmonate treatment, a single nuclei RNAseq that identified differentiating resin canal epithelial cells and are completing analysis for a QTL study in a hybrid pine population. The time course identified 4634 significantly down and 1890 significantly up regulated transcripts after treatment with methyl jasmonate, an inducer of new resin canal formation in the vascular cambial meristem. To analyze this large set of differentially regulated genes, we created a predictive expression network and analyzed it with random walk restart using 6 seed genes coding for transcription factors regulating xylem differentiation in poplar. Of the top ranked 200 transcripts, 119 transcripts were significant differentially expressed supporting these transcripts as potential candidates regulating resin canal formation. Analysis of single nuclei sequencing of shoot tips that contain differentiating resin canals, identified 10 clusters. One cluster was highly enriched in transcripts coding for 9 of the enzymes in the MEP pathway 3 prenyl synthetases, and 3 monoterpene synthases strongly suggesting that this cluster represents resin canal epithelial cells. We are mining the additional transcripts to create a trajectory analysis. In summary, we have identified > 10 novel genes that are strongly supported candidates for further analysis in breeding lines and for genetic engineering over- and under- expressing lines to increase wood terpene content to improve resistance to insect and fungal pathogens while simultaneously increasing terpene supplies for renewable chemicals and biofuels.

59 BASIC BIOLOGICAL SCIENCES↗

Genetic variations and their interaction with thirdhand smoke exposure on anxiety and memory in Collaborative Cross mice

Thirdhand smoke (THS) is linked to adverse health effects, but the effect of genetic variations on behavioral outcomes is poorly understood. To investigate this, we assessed anxiety- and memory-related behaviors in 820 mice from 21 strains of the genetically diverse Collaborative Cross (CC) mouse that were exposed to THS from 4 through 10 weeks of age. Anxiety was evaluated with a light/dark box assay with a previously established risk score system. Females were generally more sensitive: THS reduced anxiety risk in strains CC013, CC019, and CC051, but increased risk in CC036 and CC061, while males showed no significant effects. Memory was tested using passive avoidance: impairments were observed in both sexes in CC016 and CC019, with sex-dependent effects in CC002 and CC051. A genome-wide association study identified 2,347 SNPs associated with anxiety and 1,568 SNPs with memory, with 32 and 85 SNPs, respectively, interacting with THS exposure. Enrichment analyses revealed distinct biological processes underlying susceptibility, including axonogenesis, synapse organization, cognition, and learning and memory. KEGG pathway analysis identified distinct genetic pathways, including GTPase binding and GTPase regulatory activity, that act as critical molecular switches in the brain that regulate synaptic plasticity, dendritic spine structure, and neuronal signaling, directly influencing anxiety-like behaviors and memory formation. These findings show that THS exposure affects neurobehavioral outcomes in a sex- and genotype-dependent manner, highlighting critical gene-environment interactions and providing a foundation for mechanistic insights into THS neurotoxicity

Anxiety↗