Results for “RNA-seq”

NASA NTRSDate not supplied

A Pipeline for Assessing the Quality of Rna-Seq Datasets in GeneLab

Transcriptome profiling by RNA sequencing (RNA-seq) is a powerful approach to identify gene expression changes in organisms exposed to unique environments such as spaceflight. One of the challenges of evaluating RNA-seq data both within and across different space-relevant studies is the ability to control for technical differences, including the use of different library preparation kits, sequencing platforms, RNA yield, and person-to-person variation. To help address this issue, the National Institute of Standards and Technology (NIST, nist.gov) initiated a consortium, at the request of industry and academia, to develop a set of controls for gene expression measurements. The result was a set of 92 unlabeled, polyadenylated transcripts that range from 250 – 2,000 nucleotides in length to mimic natural eukaryotic mRNAs. These External RNA Controls Consortium (ERCC) genes can be used in any RNA-seq experiment, by adding known concentrations of the ERCC genes to samples after RNA extraction, to offer a standard measurement for data comparison. At NASA GeneLab, we employ these controls as part of our standard operating procedures for every in-house RNA-seq study to assess the limit of detection, dynamic range, and power of differential expression analysis both within and across experiments. Here we will discuss the use, benefits, and limitations of ERCC genes and other types of controls, such as universal RNA references, to generate quality control information for RNA-seq studies conducted at GeneLab.

GeneLab↗

NASA NTRS2019

The Omics of Stem Cell Mediated Regeneration: A Pilot Single Cell RNA-Seq Study of Mechanotransduction

Mechanical forces are potent modulators of stem cell based tissue regenerative mechanisms, inducing cell fate decisions and tissue specific commitment. A unique platform for investigating mechanotransduction is spaceflight, where microgravity and altered fluid mechanics provide a loading-null experimental condition. Seminal investigations of regenerative capacity in a wholly regenerative species, the newt model, and in a variety of totipotent and adult stem cell populations have demonstrated the detrimental effects of unloading on maintenance of stem cell based regeneration. Of particular interest is the observation that unloading interferes with the transition of stem cell pools from proliferative state to differentiation commitment. In this work we sought to test the hypothesis that gravity mechanotransduction regulates stem cell tissue regenerative processes by modulating stem cell proliferation and differentiation fates at specific cell cycle stages. To do this, clonally-derived ESCs were plated on a collagen matrix and expanded for 36 hours before re-plating on a non-adherent culture dish in the absence of leukemia inhibitory factor (LIF) to form spheroid aggregate EBs. After formation, the EBs were transferred to a collagen matrix coated culture dishes and given 4 days to allow implantation and outgrowth. In parallel, totipotent ESCs were plated 24 hours before mechanical stimulation on collagen matrix culture dishes in the presence of LIF to maintain totipotency and serve as un-differentiation committed controls. The EBs and ESCs were then subjected to either a 60 minute pulse of gravity (static loading) or 60 minutes of cyclic stretch (dynamic loading) mechanotransduction. Six hours post-stimulation, we used a 10X Genomics Single Cell controller to generate bar-coded single cell Illumina libraries and sequenced expressomes for 5,000 static loaded cells, representative of a change in gravity mechanotransduction, 5,000 dynamic loaded cells, representative of tissue loading associate with physiologic function, and 5,000 unstimulated 1g control cells. The comparison of these 3 libraries by cluster assignment based on like gene expression patterns show substantial alteration in cluster geometry due to mechanical loading. Specifically the mechanically loaded EB outgrowth cells to retain potency markers (PAX6, SOX2, CD34) and suppress early commitment markers (Dhh, VCAN, Igf1). Whereas the EBs cultured under the non-stimulated conditions display clear departure from the ESC expressome with lineage commitment markers upregulated and several tissue specific markers being expressed (BMP "early musculoskeletal development, Mesp1" early cardiovascular cell lineage). These markers are not seen in the mechano-stimulated cultures or the totipotent ESC cultures. Comparison of like clusters between our experimental conditions revealed an array of regenerative and stem cell genes are significantly mechano-regulated. Of particular importance CDKN1a/p21, a gene shown by previous investigation of our research team to be significantly upregulated in unloading, was suppressed in the static and dynamic loaded EBS. In addition to CDKN1a/p21 many genes related to cell cycle and transitory differentiation markers had elevated expression in the mechano-stimulated EBs, but surprisingly these trends were not observed in the ESC cultures. This study is the first of its kind investigating for mechano-signaling and mechano-regulated pathways, and has alre

Single cell RNA-sequecing↗

NASA NTRSDate not supplied

Combining RNA-SEQ Datasets from NASA GENELAB: An Evaluation of Correction Methods

Background: 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. Methods: 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, the 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. Results: The results showed that the reference-based approach introduced several additional (and likely artificial) differentially expressed genes when compared with the respective standard approach. Conclusions: Of the methods tested, standard ComBat_seq and DESeq2 were identified as the most robust correction methods for combining spaceflight mouse liver RNAseq datasets hosted on GeneLab.

Finsam Samson↗

NASA NTRSDate not supplied

Differential Responses to Mechanostimulation in Embryonic Stem Cells Versus the Embryoid Body Model of Development Assessed at Single Cell RNA-Seq Resolution

Mechanicalforces generated by gravity have shaped life on Earth and impact gene expression and morphogenesis during early development. In contrast disuse canreduce normal mechanical loading, resulting in altered cell and tissue function. Although loading in adult mammals is known to promote increased cell proliferation and differentiation, little is known about how cells respondto this stimulusduring early development. In this study we sought to understand, with single cell RNA-sequencing resolution, how a 60-minute pulse of 50xg hypergravity-generated 5kPa hydrostatic pressure, influences transcriptomic regulation of developmental processes in the Embryoid Body (EB) model. Our study included both day-9 EBs and progenitor mouse embryonic stem cells (ESCs) with or without the hydrostatic pressurepulse. Single cell tSNE mapping shows limited transcriptome shifts in response to thispulse in either ESCs or EBs; this pulse,however, induces greater positional shifts in EB mapping compared to ESCs, indicating the influence of mechanotransduction is more pronounced in later states of cell commitment within the developmental program.We assessed ESCs and EBs for differentially expressed (DE) genes with hydrostatic pressurepulse and found approximately 1/3 DE genes were shared. However, gene ontology (GO) pathway analysis show that EBs have choreographed responses associated with upregulation ofpathways formulticellular development, mechanical signal transduction, and DNA damage repair. Cluster transcriptome analysis of the EBs showsmechanostimulationpromotes maintenance of transitory cell phenotypes in early development,including EB cluster co-expression of markers for progenitor, post-implant epiblast and primitive endoderm phenotypes versus expression exclusivity in the non-pulsed clusters. Pseudotime analysisidentified three branching cell types susceptible tohydrostatic pressureinduction of cell fate decisions. In summary, this study provides novel evidence that ESC maintenance and EB development can be regulated by mechanostimulation,and that stem cells committed to a differentiation program are more sensitive to force-induced changes to their transcriptome.

Cassandra Juran↗

NASA NTRSDate not supplied

Space Flown Rodent Liver RNA Sequencing Data for Machine Learning in Space Biology Research

High-throughput nucleic acid sequencing (DNA-seq, RNA-seq) has become widespread in biomedical research due to the growing availability and affordability of these assays. Data analysis has been accelerated in recent years by the adoption of artificial intelligence (AI) and machine learning (ML) techniques by biomedical researchers. In space biology research, RNAseq datasets from space-flown experimental samples are critical for characterizing the gene expression aberrations associated with exposure to spaceflight stressors. However, space biological experiments tend to be very low sample size, so identifying proper AI/ML algorithms for sequencing data analysis is an ongoing challenge since these algorithms typically require large sample size. The NASA Science Mission Directorate (SMD) has started the “Benchmark Initiative for AI/ML”, focused on creating datasets meant for three main applications: 1) scientific benchmarking, which finds the best algorithm for a specific problem; 2) application benchmarking, which measures algorithm performance against a set of parameters; and 3) system benchmarking, which evaluates performance of hardware and software architecture. These scientific benchmarks consist of an AI-ready dataset and a reference implementation on a specific scientific question. In this work, we focused on generating standardized datasets to allow the scientific community to benchmark AI/ML algorithms in the domain of space biology. We present here a standardized, AI-ready, publicly available benchmark dataset for space biology RNA-seq data as a collaboration between the NASA AI4LS (Artificial Intelligence for Life Sciences) working group. and NASA’s SMD. This dataset consists of space-flown and ground control mouse liver found in the NASA GeneLab omics database. However, to amplify the small sample number (n=112 samples) for ML purposes, we employ Gaussian noise and a generative adversarial network to extend this dataset to 6,000 synthetic samples, matching the original gene expression characteristics.

James Casaletto↗

NASA NTRS2016

Transcriptomic Response of Drosophila Melanogaster Pupae Developed in Hypergravity

The metamorphosis of Drosophila is evolutionarily adapted to Earth's gravity, and is a tightly regulated process. Deviation from 1g to microgravity or hypergravity can influence metamorphosis, and alter associated gene expression. Understanding the relationship between an altered gravity environment and developmental processes is important for NASA's space travel goals. In the present study, 20 female and 20 male synchronized (Canton S, 2 to 3day old) flies were allowed to lay eggs while being maintained in a hypergravity environment (3g). Centrifugation was briefly stopped to discard the parent flies after 24hrs of egg laying, and then immediately continued until the eggs developed into P6-staged pupae (25 - 43 hours after pupation initiation). Post hypergravity exposure, P6-staged pupae were collected, total RNA was extracted using Qiagen RNeasy mini kits. We used RNA-Seq and qRT-PCR techniques to profile global transcriptomic changes in early pupae exposed to chronic hypergravity. During the pupal stage, Drosophila relies upon gravitational cues for proper development. Assessing gene expression changes in the pupa under altered gravity conditions helps highlight gravity dependent genetic pathways. A robust transcriptional response was observed in hypergravity-exposed pupae compared to controls, with 1,513 genes showing a significant (q < 0.05) difference in gene expression. Five major biological processes were affected: ion transport, redox homeostasis, immune response, proteolysis, and cuticle development. This outlines the underlying molecular changes occurring in Drosophila pupae in response to hypergravity.

RNASeq↗

NASA NTRS2018

GeneLab: Omics Database for Spaceflight Experiments

Motivation - To curate and organize expensive spaceflight experiments conducted aboard space stations and maximize the scientific return of investment, while democratizing access to vast amounts of spaceflight related omics data generated from several model organisms. Results - The GeneLab Data System (GLDS) is an open access database containing fully coordinated and curated "omics" (genomics, transcriptomics, proteomics, metabolomics) data, detailed metadata and radiation dosimetry for a variety of model organisms. GLDS is supported by an integrated data system allowing federated search across several public bioinformatics repositories. Archived datasets can be queried using full-text search (e.g., keywords, Boolean and wildcards) and results can be sorted in multifactorial manner using assistive filters. GLDS also provides a collaborative platform built on GenomeSpace for sharing files and analyses with collaborators. It currently houses 172 datasets and supports standard guidelines for submission of datasets, MIAME (for microarray), ENCODE Consortium Guidelines (for RNA-seq) and MIAPE Guidelines (for proteomics).

omics↗

NASA NTRS2019

Deciphering the Effects of Microgravity Cell by Cell

Forces generated by gravity have a profound impact on the behavior of cells in tissues affecting the course of the cell cycle and differentiation fate of progenitors in mammalian tissues. These cells are contributing to normal tissue regenerative health and defence against disease. In Human space exploration context, it is extremely important to determine spaceflight provoked changes in tissue's regenerational capabilities. Microgravity experienced during spaceflight causes unloading and mechanical disuse on all orthostatic support tissues, therefore impacting stem cell fate and lineage commitment decisions. Investigating how ESCs respond to mechanical stimulation is a platform for fundamental developmental and regeneration research applicable for spaceflight. However, the gene expression programs associatiated with early committment stem cell pathways in response to physical stimulation are not readily known. Single-cell RNA-seq technologies have recently revolutionized the world of molecular biology by providing the capability to assess gene expression pattern within a single cell. Our method isolates and separately barcodes mRNAs from thousands of single cells and sequences their expressomes. Understanding regenerative processes on a molecular level would not only help reduce long-term spaceflight impact on health, but also may enable the development of novel tissue regenerative approaches to tissue degeneration on Earth.

Single cell sequencing↗

NASA NTRS2019

GeneLab Analysis Working Group Pipelines

GeneLab must establish data processing pipelines for common data types including microarray, RNA-sequencing, and metagenomic profiling. Here we give an overview of current microarray and RNA-seq pipelines and discuss future pipelines including metagenomic profiling pipelines

Galazka, Jonathan M.↗

NASA NTRSDate not supplied

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↗

NASA NTRSDate not supplied

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↗

NASA NTRSDate not supplied

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↗

NASA NTRSDate not supplied

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↗

NASA NTRSDate not supplied

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↗

NASA NTRSDate not supplied

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↗

NASA NTRSDate not supplied

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↗

NASA NTRS2020

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↗

NASA NTRSDate not supplied

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↗