Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “RNA-seq”

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.

35 records · Page 2

Novel insights enabled by combining mouse muscle datasets from the Rodent Research-1 mission

Biological space experiments are often expensive and difficult to conduct. As such, it is critical to maximize the value of the data that is collected during these experiments. One way to do this is to combine multiple–previously separate–datasets. This can increase the number of replicates for the conditions of interest (and hence statistical power), allow new multi-factor questions to be asked, and potentially highlight new patterns that otherwise would not have been identified from single-dataset studies. However, the process of combining datasets introduces noise due to inherent technical variations between experiments. To better understand the insights that can be gained from multi-dataset analyses and the problems that may arise from joining multiple datasets, several mouse muscle RNA-Seq datasets from the Rodent Research-1 mission were first selected. Then, using the R package DESeq2, principal component analysis (PCA) plots and differentially expressed gene (DEG) lists between ground and flight muscle samples were generated for individual datasets and for different pairwise combinations of datasets. Several new DEGs were identified in the combined datasets, and patterns in the PCA plots were affected depending on which datasets were joined. Understanding the results of this work will be critical for future studies that seek to perform multi-dataset analyses.

spaceflight↗

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↗

Translational Line of Sight: a multi-omics longitudinal study of the murine retinal response to spaceflight hazard analogs

The space environment includes unique hazards like radiation and microgravity which adversely affect physiology and behavior of humans and rodent models. To better characterize the retinal response to spaceflight, we assessed a multi-omics NASA GeneLab dataset where 6-month-old female mice were gamma irradiated and/or hindlimb unloaded for 21 days followed by whole transcriptome shotgun sequencing (RNA-Seq) and reduced representation bisulfite sequencing (RRBS) of retina samples collected at 7 days, 1 month or 4 months post-exposure. We compared time-matched epigenomic and transcriptomic retinal profiles revealing a total of 4,178 differentially methylated loci or regions, and 457 differentially expressed genes. Highest correlation in methylation differences was seen across different conditions at the same time point (e.g., between radiation exposure and hindlimb unloaded at 7 days). Biological processes related to nucleotide metabolism were enriched in all groups with activation at 1 month and suppression at 7 days and 4 months. Genes and processes related to Notch and Wnt signaling showed alterations 4 months post-exposure. Interestingly, Notch3 and Lrg1 showed differential patterns in the NASA Twins Study in-flight samples and in response to stressors in the murine retina in the current study. A total of 23 genes were both differentially methylated and expressed, including genes involved in retinal disease or cataract development (Crybb3, Fgfr1, Pitpnm3, Sipa1l3, Sox9) and inflammatory response (B4galt6, Ppm1a, Sphk1). To our knowledge, the current multi-omics analysis is the first multi-omics study to interrogate the epigenomic and transcriptomic impacts of radiation and hindlimb unloading on the retina in isolation and in combination. The results provide an insight into the retinal response to individual spaceflight hazard analogs and their interplay at different post-exposure stages and contributes towards a mechanistic understanding of spaceflight-induced vision impairment using ground-based models.

Prachi Kothiyal↗

MULTI-OMICS ANALYSIS OF THE IMPACT OF CHRONIC LOW-DOSE RADIATION AND HINDLIMB SUSPENSION ON MURINE BRAIN AND RETINA

The space environment includes hazards like radiation and microgravity which can adversely affect biological systems. We assessed multi-omics multi-tissue NASA GeneLab datasets where 6-month-old female mice were gamma irradiated (IR) and/or hindlimb unloaded (HLU) for 21 days. Whole transcriptome shotgun sequencing (RNA-Seq) and reduced representation bisulfite sequencing (RRBS) of brain and retina samples collected at 4 months post-exposure was performed to better characterize the retinal and neurological responses to spaceflight. We compared epigenomic and transcriptomic profiles within each exposure group for both tissue types to identify correlation (Pearson’s correlation test; p-value < 0.05) between gene expression and DNA methylation levels that may be related to transcriptional regulation. We then obtained genes with methylation-expression correlation that also showed differences in mean expression or dispersion between exposed and control groups (adjusted p-value < 0.25; relaxed to denote ‘hypothesis’) in the brain (37 genes in HLU, 4 in IR, and 156 in HLU+IR) or retina (92 genes in HLU, 1 in IR, and 55 in HLU+IR). Enriched Gene Ontology (GO) terms for these genes are listed in Table 1 for HLU and HLU+IR for both tissue types. No enriched terms and only a few genes were detected with IR-only exposure in the brain (Chmp1a, Limd1, Rab40b, Ubc) and retina (retinoblastoma binding protein Rbbp7). Cellular components related to synapse were enriched in both tissue types. Previous analysis of differentially expressed genes in the retina after 1 month of HLU+IR showed enrichment in the somatodendritic compartment of the neuron, which was also observed in the brain 4 months post-exposure. Interestingly, genes related to ubiquitination showed correlation between methylation and expression and were differentially expressed or dispersed in different exposure groups indicating that this pathway may play an important role in multi-stressor response (Figure 1). The current multi-omics and multi-tissue analysis interrogates the epigenomic and transcriptomic impacts of radiation and hindlimb unloading, in isolation and in combination, on the retina and the brain. The results provide insight into the adaptive response to individual spaceflight hazard analogs and their interplay, as well as hypotheses to be further tested for understanding spaceflight-induced neurological and vision effects.

Prachi Kothiyal↗

Transcriptomics-based Machine Learning Analysis Predicts Space-Exposed Murine Livers

Limited sample sizes, high data dimensionality, and sensitivity to technical and biological variability of next generation sequencing (NGS), typically limits machine learning (ML) approaches in spaceflight studies that include radiation effects. However, pooling smaller studies while addressing intra- and inter-study variabilities allows for ML predictive modeling. Here, integration methods were applied to whole transcriptome shotgun sequencing (RNA-seq) data from six mouse liver GeneLab datasets (GLDS) (n ranging from 6 to 39 samples) from with a total of 81 spaceflight and ground-control samples to determine top features (i.e. genes) relevant to spaceflight including the effect of radiation exposure. RNASeq counts were normalized for each study, then merged and scaled across all datasets. Data dimensionality was reduced using a minimum redundancy maximum relevance (MRMR) methodology. Redundancy and relevance were computed using the Pearson correlation and F-statistic, respectively. The top 100 MRMR features were used to predict spaceflight vs. ground-control samples using Random Forest (RF), Support Vector Machine (SVM), and Linear Discriminant Analysis (LDA) classifiers with 5-fold cross validation (CV). Principal component analysis (PCA) on the complete feature set versus the MRMR features shows separation between spaceflight samples and ground controls (Figure 1A). The ML-based gene sets were compared against differential gene expression results obtained with DESeq2 from individual GLDS. Using all features or randomly sampled subsets at matching set sizes with MRMR, a maximum classifier accuracy of 69% was shown on the test set over 5 folds. For all classifiers, CV training using at least the top 30 MRMR genes show minimum 89% accuracy and 0.95 AUC value on the test set over 5 folds (Figure 1B). Baseline set analysis on differentially expressed genes (DEGs) identified using padj ≤ 0.05 show 295 DEGs that overlap at least two studies and 13 DEGs that overlap three studies (Figure 1C). Set analysis between the top 100 MRMR features and the DEGs showed 47 genes that overlap at least one study and 24 genes that overlap two studies. Over-representation analysis showed overlapping biological processes related to fatty acid and lipid metabolism which may indicate these processes in the response to spaceflight stressors. MRMR feature selection for the selected ML methods improve performance relative to a classifier built on all features or randomly sampled subsets. Permutation feature importance within the decorrelated MRMR features showed concordance in feature ranking between ML methods. A challenge of applying ML methods across heterogeneous NGS data is accounting for signal:noise. Here, signal validation across studies was shown by intersecting sets between top MRMR genes and DEGs from DESeq2 analysis. Non-intersecting sets introduce opportunity to explore genes relevant to differentiating space flight exposed groups and implementing ML methods across existing NGS datasets may overcome sample size limitations.

Machine Learning↗

A Multi-omics Longitudinal Study of the Murine Retinal Response to Chronic Low-dose Irradiation and/or Simulated Microgravity

The space environment includes unique hazards like radiation and microgravity which adversely affect physiology and behavior of humans and rodent models. To better characterize the retinal response to spaceflight, we assessed a multi-omics NASA GeneLab dataset where 6-month-old female mice were gamma irradiated and/or hindlimb unloaded for 21 days followed by whole transcriptome shotgun sequencing (RNA-Seq) and reduced representation bisulfite sequencing (RRBS) of retina samples collected at 7 days, 1 month or 4 months post-exposure. We compared time-matched epigenomic and transcriptomic retinal profiles revealing a total of 4,178 differentially methylated loci or regions, and 457 differentially expressed genes. Highest correlation in methylation differences was seen across different conditions at the same time point (e.g., between radiation exposure and hindlimb unloaded at 7 days). Biological processes related to nucleotide metabolism were enriched in all groups with activation at 1 month and suppression at 7 days and 4 months. Genes and processes related to Notch and Wnt signaling showed alterations 4 months post-exposure. Interestingly, Notch3 and Lrg1 showed differential patterns in the NASA Twins Study in-flight samples and in response to stressors in the murine retina in the current study. A total of 23 genes were both differentially methylated and expressed, including genes involved in retinal disease or cataract development (Crybb3, Fgfr1, Pitpnm3, Sipa1l3, Sox9) and inflammatory response (B4galt6, Ppm1a, Sphk1). To our knowledge, the current multi-omics analysis is the first multi-omics study to interrogate the epigenomic and transcriptomic impacts of radiation and hindlimb unloading on the retina in isolation and in combination. The results provide an insight into the retinal response to individual spaceflight hazard analogs and their interplay at different post-exposure stages and contributes towards a mechanistic understanding of spaceflight-induced vision impairment using ground-based models.

Prachi Kothiyal↗

Transcriptomics-based Machine Learning (ML) Analysis Predicts Space-Exposed Murine Livers

Limited sample sizes, high data dimensionality, and sensitivity to technical and biological variability of next generation sequencing (NGS), typically limits machine learning (ML) approaches in spaceflight studies that include radiation effects. However, pooling smaller studies while addressing intra- and inter-study variabilities allows for ML predictive modeling. Here, integration methods were applied to whole transcriptome shotgun sequencing (RNA-seq) data from six mouse liver GeneLab datasets (GLDS) (n ranging from 6 to 39 samples) from with a total of 81 spaceflight and ground-control samples to determine top features (i.e. genes) relevant to spaceflight including the effect of radiation exposure. RNASeq counts were normalized for each study, then merged and scaled across all datasets. Data dimensionality was reduced using a minimum redundancy maximum relevance (MRMR) methodology. Redundancy and relevance were computed using the Pearson correlation and F-statistic, respectively. The top 100 MRMR features were used to predict spaceflight vs. ground-control samples using Random Forest (RF), Support Vector Machine (SVM), and Linear Discriminant Analysis (LDA) classifiers with 5-fold cross validation (CV). Principal component analysis (PCA) on the complete feature set versus the MRMR features shows separation between spaceflight samples and ground controls (Figure 1A). The ML-based gene sets were compared against differential gene expression results obtained with DESeq2 from individual GLDS. Using all features or randomly sampled subsets at matching set sizes with MRMR, a maximum classifier accuracy of 69% on the test set over 5 folds. For all classifiers, CV training using at least the top 30 MRMR genes show minimum 89% accuracy and 0.95 AUC value on the test set over 5 folds (Figure 1B). Baseline set analysis on differentially expressed genes (DEGs) identified using padj ≤ 0.05 show 295 DEGs that overlap at least two studies and 13 DEGs that overlap three studies (Figure 1C). Set analysis between the top 100 MRMR features and the DEGs showed 47 genes that overlap at least one study and 24 genes that overlap two studies. Over-representation analysis showed overlapping biological processes related to fatty acid and lipid metabolism which may indicate these processes in the response to spaceflight stressors. MRMR feature selection for the selected ML methods improve performance relative to a classifier built on all features or randomly sampled subsets. Permutation feature importance within the decorrelated MRMR features showed concordance in feature ranking between ML methods. A challenge of applying ML methods across heterogeneous NGS data is accounting for signal:noise. Here, signal validation across studies was shown by intersecting sets between top MRMR genes and DEGs from DESeq2 analysis. Non-intersecting sets introduce opportunity to explore genes relevant to differentiating space flight exposed groups and implementing ML methods across existing NGS datasets may overcome sample size limitations.

Machine Learning↗

The NASA Twins Study: The Effect of One Year in Space on Long-Chain Fatty Acid Desaturases and Elongases

Background: To date, there is no clear understanding of the effect of long-duration spaceflight on the major enzymes that govern the metabolism of omega-6 and omega-3 fatty acids. To address this gap in knowledge, we used data from the NASA Twins Study, which includes a multi-scale omic investigation of the changes that occurred during a year-long (340 days) human spaceflight. Embedded within the NASA Twins data are specific analytes associated with fatty acid metabolism. Objectives: To examine the long-chain fatty acid desaturases and elongases in a single human during one year in space. Method: One male twin was on board the International Space Station (ISS) for one year, while his monozygotic twin served as a genetically matched ground control. Longitudinal assessments included the genome, epigenome, transcriptome, proteome, metabolome, microbiome, and immunome during the mission, as well as six months before and after. The gene-specific fatty acid desaturase and elongase transcriptome data (FADS1, FADS2, ELOVL2 and ELOVL5) were extracted from untargeted RNA-seq measurements derived from white blood cell fractions. Results: Most data from the elongases and desaturases exhibited relatively similar expression profiles (R2>0.6) over time for the CD8, CD19, and LD cell fractions, indicating overall conservation of function within and between the subjects. Both cell-type and temporal specificity was observed in some cases, and some differences were also apparent between the poly-adenylated fraction (polyA) of processed RNAs vs. the ribo-depleted (ribo-) fraction. The flight subject showed a stronger enrichment of the Fatty Acid Metabolic processes pathway across almost all cell types (columns, CD4, CD8, CPT, LD), most especially in the ribodepleted fraction of RNA, but also with the polyA+ fraction of RNA. GSEA enrichment measures across three related Fatty Acid Metabolism pathways showed a differential between the ground and flight subject. Conclusions: There appears to be no persistent alteration of desaturase and elongase gene expression associated with one year in space. However, these data provide evidence that cellular lipid metabolism can be responsive and dynamic to spaceflight, even though it appears cell-type- and context-specific, most notably in terms of the fraction of RNA measured and the collection protocols. These results also provide new evidence of mid-flight spikes in expression of selected genes, which may indicate transient responses to specific insults during spaceflight.

Elongase↗

Transcriptomics-based Machine Learning Analysis Predicts Space-Exposed Murine Livers

Limited sample sizes, high data dimensionality, and sensitivity to technical and biological variability of next generation sequencing (NGS), typically limits machine learning (ML) approaches in spaceflight studies that include radiation effects. However, pooling smaller studies while addressing intra- and inter-study variabilities allows for ML predictive modeling. Here, integration methods were applied to whole transcriptome shotgun sequencing (RNA-seq) data from six mouse liver GeneLab datasets (GLDS) (n ranging from 6 to 39 samples) from with a total of 81 spaceflight and ground-control samples to determine top features (i.e. genes) relevant to spaceflight including the effect of radiation exposure. RNASeq counts were normalized for each study, then merged and scaled across all datasets. Data dimensionality was reduced using a minimum redundancy maximum relevance (MRMR) methodology. Redundancy and relevance were computed using the Pearson correlation and F-statistic, respectively. The top 100 MRMR features were used to predict spaceflight vs. ground-control samples using Random Forest (RF), Support Vector Machine (SVM), and Linear Discriminant Analysis (LDA) classifiers with 5-fold cross validation (CV). Principal component analysis (PCA) on the complete feature set versus the MRMR features shows separation between spaceflight samples and ground controls (Figure 1A). The ML-based gene sets were compared against differential gene expression results obtained with DESeq2 from individual GLDS. Using all features or randomly sampled subsets at matching set sizes with MRMR, a maximum classifier accuracy of 69% was shown on the test set over 5 folds. For all classifiers, CV training using at least the top 30 MRMR genes show minimum 89% accuracy and 0.95 AUC value on the test set over 5 folds (Figure 1B). Baseline set analysis on differentially expressed genes (DEGs) identified using padj ≤ 0.05 show 295 DEGs that overlap at least two studies and 13 DEGs that overlap three studies (Figure 1C). Set analysis between the top 100 MRMR features and the DEGs showed 47 genes that overlap at least one study and 24 genes that overlap two studies. Over-representation analysis showed overlapping biological processes related to fatty acid and lipid metabolism which may indicate these processes in the response to spaceflight stressors. MRMR feature selection for the selected ML methods improve performance relative to a classifier built on all features or randomly sampled subsets. Permutation feature importance within the decorrelated MRMR features showed concordance in feature ranking between ML methods. A challenge of applying ML methods across heterogeneous NGS data is accounting for signal:noise. Here, signal validation across studies was shown by intersecting sets between top MRMR genes and DEGs from DESeq2 analysis. Non-intersecting sets introduce opportunity to explore genes relevant to differentiating space flight exposed groups and implementing ML methods across existing NGS datasets may overcome sample size limitations.

Machine Learning↗

MULTI-OMICS STUDY OF THE EFFECT OF REDOX-ACTIVE METALLOPORPHYRIN ON MURINE RETINA DURING SPACEFLIGHT

Astronauts returning from spaceflight have experienced eye problems, which may decrease retinal performance and lead to long-term effects on visual acuity. This study leverages the collected data from spaceflown murine retinas that were treated with redox-active metalloporphyrin (BuOE) to mitigate spaceflight-induced changes and respective ground controls. 10-week-old adult C57BL/6 male mice (n=5 in each of BuOE treated and saline control groups for spaceflown and ground control samples) were flown on Space-X 24 to the ISS national lab, kept in low earth orbit for 35 days and returned to Earth alive. Our multi-omics analysis of RNA-sequencing and reduced representation bisulfite sequencing (RRBS) data generated from subsequent murine retina tissues uncovered genes, pathways, and epigenetic modifications consistent with therapeutic potential of BuOE. From RNA-Seq analysis of spaceflown murine samples, the treatment group show differentially expressed genes relative to saline controls that reached significance (adjusted p-value < 0.05) and included genes Gpx3 and Crhbp, which are related to protection against cell oxidative damage and cellular response to organonitrogen compounds. Ranked fold-changes from the same contrast were used for gene set enrichment analysis, which showed biological processes reaching significance (adjusted p-value < 0.05) including glutathione metabolic processes and cellular response to xenobiotic stimulus. RRBS data of the spaceflown murine samples found 139 hyper or hypo differentially methylated sites spread across chromosomes 1-19 (20% promoters, 21% exons, 43% introns | 20 CpG islands, 7 CpG shores) with a 10% methylation difference (q-value < 0.05).The findings from this investigation have the potential to provide valuable insights into the molecular mechanisms underlying conditions like spaceflight associated neuro-ocular syndrome and assess the effectiveness of BuOE as a countermeasure for astronauts experiencing neuro-ophthalmic abnormalities, which can lead to long-term effects on visual acuity.

Biostatistics↗

Transcriptomics Processing Pipelines for Space Biology: An Open Source and Consensus-Driven Approach

Transcriptomics holds significant value in elucidating the relationship between gene expression, experimental factors, biological factors, and various types of omics data. Enhancing our understanding of these connections is paramount for foundational biology, which plays a pivotal role in devising solutions for challenges pertinent to both space travel and terrestrial life. The NASA GeneLab project, part of the Open Science Data Repository (OSDR.nasa.gov), seeks to accelerate space biology research through cataloging and democratizing ‘omics data, including transcriptomics. Since raw omics data are largely inaccessible to non-bioinformaticians, GeneLab works with the scientific community via the Open Science Analysis Working Groups (AWGs) to develop standard processing pipelines to generate and publish processed data. Unlike raw data, processed data have greater immediate value to diverse users with varying technical backgrounds and computational capabilities. Standardizing processing workflows is essential to match the pace of raw data generation, ensure reproducibility, and enable standardized processed data for comparison across datasets. As of June 2023, transcriptomics studies comprise over half of GeneLab datasets hosted on the OSDR, including data from bulk RNA-seq and Affymetrix or Agilent 1-Channel DNA microarray assays. In collaboration with the AWGs, GeneLab developed consensus processing pipelines for these transcriptomics data types that includes quality control, background correction (microarray only), data normalization and quantification, culminating in the detection and annotation of differentially expressed genes. The work presented here describes Nextflow implementations of GeneLab’s consensus transcriptomics pipelines that automates and accelerates processing of these datasets. In addition to the core data processing, these workflows also include raw data staging and a robust verification and validation program to identify errors in real-time, stop additional downstream computation, and preserve computational resources. These workflows are used to generate GeneLab processed data hosted on the OSDR, and are publicly available as open source software for others to use at: https://github.com/nasa/GeneLab_Data_Processing.

Jonathan Oribello↗

Multi-Omics Study of the Effect of Redox-Active Metalloporphyrin on Murine Retina During Spaceflight

Astronauts returning from spaceflight have experienced eye problems, which may decrease retinal performance and lead to long-term effects on visual acuity. This study leverages the collected data from spaceflown murine retinas that were treated with redox-active metalloporphyrin (BuOE) to mitigate spaceflight-induced changes and respective ground controls. 10-week-old adult C57BL/6 male mice (n=5 in each of BuOE treated and saline control groups for spaceflown and ground control samples) were flown on Space-X 24 to the ISS national lab, kept in low earth orbit for 35 days and returned to Earth alive. Our multi-omics analysis of RNA-sequencing and reduced representation bisulfite sequencing (RRBS) data generated from subsequent murine retina tissues uncovered genes, pathways, and epigenetic modifications consistent with therapeutic potential of BuOE. From RNA-Seq analysis of spaceflown murine samples, the treatment group show differentially expressed genes relative to saline controls that reached significance (adjusted p-value < 0.05) and included genes Gpx3 and Crhbp, which are related to protection against cell oxidative damage and cellular response to organonitrogen compounds. Ranked fold-changes from the same contrast were used for gene set enrichment analysis, which showed biological processes reaching significance (adjusted p-value < 0.05) including glutathione metabolic processes and cellular response to xenobiotic stimulus. RRBS data of the spaceflown murine samples found 139 hyper or hypo differentially methylated sites spread across chromosomes 1-19 (20% promoters, 21% exons, 43% introns | 20 CpG islands, 7 CpG shores) with a 10% methylation difference (q-value < 0.05).The findings from this investigation have the potential to provide valuable insights into the molecular mechanisms underlying conditions like spaceflight associated neuro-ocular syndrome and assess the effectiveness of BuOE as a countermeasure for astronauts experiencing neuro-ophthalmic abnormalities, which can lead to long-term effects on visual acuity.

Biostatistics↗

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↗

Metabolic Stress in Space: ROS-Induced Mutations in Mice Hint at A New Path to Cancer

Long-duration spaceflight beyond Earth's magnetosphere poses serious health risks, including muscle atrophy, bone loss, liver and kidney damage, and the Spaceflight-Associated Neuro-ocular Syndrome (SANS). RNA-seq of mice aboard the International Space Station (ISS) for 37 days revealed extraordinary hypermutation in tissue-specific genes, with guanine base conversion predominating, potentially contributing to spaceflight-associated health risks. Our results suggest that the genome-wide accelerated mutation that we measured, seemingly independent of radiation dose, was induced by oxidative damage from higher atmospheric carbon dioxide (CO 2 ) levels and increased reactive oxygen species (ROS) on the ISS. This accelerated mutation, faster via RNA transcription than replication and more numerous than by radiation alone, unveils novel hotspots in the mammalian proteome. Notably, these hotspots correlate with commonly mutated genes across various human cancers, highlighting the ISS as a crucial platform for studying accelerated mutation, genome instability, and the induction of disease-causing mutations in model organisms. Our results suggest that metabolic processes can contribute to somatic mutation, and thus may play a role in the development of cancer. A metabolic link to genetic instability potentially has far-reaching implications for various diseases, with implications for human health on Earth and in space.

mice↗

Evaluation of Correction Methods for NASA GeneLab Transcriptomic Datasets

Conducting space biology experiments aboard the International Space Station, particularly those utilizing complex model organisms like mice, is expensive and difficult due to limited crew availability, hardware, and space. As a result, sample numbers from these studies are low, reducing the statistical power of any one experiment. Aggregating spaceflight datasets serves as a method to increase sample numbers, allowing for novel insights through bioinformatic analysis of ‘omics data from merged datasets. However, aggregating datasets can introduce unwanted variation including 1) differences in sample handling, processing, and sequencing platforms between datasets (technical variation) as well as 2) differences in experimental design between datasets such as sex or age of the model organism used. In the present study, NASA GeneLab-hosted RNAseq datasets from rodent liver tissues were used to evaluate several statistical methods to correct for this unwanted variation through two approaches, reference-based and standard. The following correction algorithms were applied with (reference-based) and/or without (standard) considering Universal Mouse RNA Reference samples: ComBat and ComBat_seq from the SVA package, median polish, empirical Bayes, and ANOVA-based algorithms from the MBatch package, and negative binomial regression normalization in the DESeq2 package. For each approach, after the correction algorithm was applied, differential gene expression (DGE) analysis of flight and ground control samples was performed with the combined data. The robustness of each tool was evaluated using BatchQC, to determine statistical differences between datasets before and after correction, Principal Component Analysis, to evaluate global gene expression in samples before and after correction, and by comparing DGE analysis of individual datasets and combined datasets before and after correction. The results showed that the reference-based approach introduced several additional (and likely artificial) DEGs when compared with the standard approach. Thus, the most robust standard correction will be implemented in the GeneLab Visualization 2.0 platform when datasets are combined.

GeneLab, RNA-seq, Batch Correction↗

Evaluation of Correction Methods for NASA GeneLab Transcriptomic Datasets

Conducting space biology experiments aboard the International Space Station, particularly those utilizing complex model organisms like mice, is expensive and difficult due to limited crew availability, hardware, and space. As a result, sample numbers from these studies are low, reducing the statistical power of any one experiment. Aggregating spaceflight datasets serves as a method to increase sample numbers, allowing for novel insights through bioinformatic analysis of ‘omics data from merged datasets. However, aggregating datasets can introduce unwanted variation including 1) differences in sample handling, processing, and sequencing platforms between datasets (technical variation) as well as 2) differences in experimental design between datasets. In the present study, NASA GeneLab-hosted RNAseq datasets from mouse liver tissues were used to evaluate several statistical methods to correct for this unwanted variation through two approaches, reference-based and standard. The following correction algorithms were applied with (reference-based) and/or without (standard) considering Universal Mouse RNA Reference samples: ComBat and ComBat_seq from the SVA package, median polish, empirical Bayes, and ANOVA-based algorithms from the MBatch package, and negative binomial regression normalization in the DESeq2 package. For each approach, after the correction algorithm was applied, differential gene expression (DGE) analysis of flight and ground control samples was performed with the combined data. The robustness of each tool was evaluated using BatchQC to determine statistical differences between datasets before and after correction, Principal Component Analysis to evaluate global gene expression in samples before and after correction, and by comparing DGE analysis of individual datasets and combined datasets before and after correction. The results showed that the reference-based approach introduced several additional (and likely artificial) DEGs when compared with the respective standard approach. Of the methods tested, standard ComBat and DESeq2 were identified as the most robust correction methods for combining spaceflight mouse liver RNAseq datasets hosted on GeneLab.

GeneLab↗

Regulatory response to a hybrid ancestral nitrogenase in Azotobacter vinelandii

Biological nitrogen fixation, the microbial reduction of atmospheric nitrogen to bioavailable ammonia, represents both a major limitation on biological productivity and a highly desirable engineering target for synthetic biology. However, the engineering of nitrogen fixation requires an integrated understanding of how the gene regulatory dynamics of host diazotrophs respond across sequence-function space of its central catalytic metalloenzyme, nitrogenase. Here, we interrogate this relationship by analyzing the transcriptome of Azotobacter vinelandii engineered with a phylogenetically inferred ancestral nitrogenase protein variant. The engineered strain exhibits reduced cellular nitrogenase activity but recovers wild-type growth rates following an extended lag period. We find that expression of genes within the immediate nitrogen fixation network is resilient to the introduced nitrogenase sequence-level perturbations. Rather the sustained physiological compatibility with the ancestral nitrogenase variant is accompanied by reduced expression of genes that support trace metal and electron resource allocation to nitrogenase. Our results spotlight gene expression changes in cellular processes adjacent to nitrogen fixation as productive engineering considerations to improve compatibility between remodeled nitrogenase proteins and engineered host diazotrophs.

nitrogen fixation↗