Search NASA⌕ Search

SEARCH · Search NASA

Results for “RNA sequencing”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

ACTINN: automated identification of cell types in single cell RNA sequencing

Abstract Motivation Cell type identification is one of the major goals in single cell RNA sequencing (scRNA-seq). Current methods for assigning cell types typically involve the use of unsupervised clustering, the identification of signature genes in each cluster, followed by a manual lookup of these genes in the literature and databases to assign cell types. However, there are several limitations associated with these approaches, such as unwanted sources of variation that influence clustering and a lack of canonical markers for certain cell types. Here, we present ACTINN (Automated Cell Type Identification using Neural Networks), which employs a neural network with three hidden layers, trains on datasets with predefined cell types and predicts cell types for other datasets based on the trained parameters. Results We trained the neural network on a mouse cell type atlas (Tabula Muris Atlas) and a human immune cell dataset, and used it to predict cell types for mouse leukocytes, human PBMCs and human T cell sub types. The results showed that our neural network is fast and accurate, and should therefore be a useful tool to complement existing scRNA-seq pipelines. Availability and implementation The codes and datasets are available at https://figshare.com/articles/ACTINN/8967116. Tutorial is available at https://github.com/mafeiyang/ACTINN. All codes are implemented in python. Supplementary information Supplementary data are available at Bioinformatics online.

Ma, Feiyang (ORCID:0000000262600787)↗
NASA NTRS1995

Structurally complex and highly active RNA ligases derived from random RNA sequences

Seven families of RNA ligases, previously isolated from random RNA sequences, fall into three classes on the basis of secondary structure and regiospecificity of ligation. Two of the three classes of ribozymes have been engineered to act as true enzymes, catalyzing the multiple-turnover transformation of substrates into products. The most complex of these ribozymes has a minimal catalytic domain of 93 nucleotides. An optimized version of this ribozyme has a kcat exceeding one per second, a value far greater than that of most natural RNA catalysts and approaching that of comparable protein enzymes. The fact that such a large and complex ligase emerged from a very limited sampling of sequence space implies the existence of a large number of distinct RNA structures of equivalent complexity and activity.

Non-NASA Center↗
NASA NTRS1997

Equally parsimonious pathways through an RNA sequence space are not equally likely

An experimental system for determining the potential ability of sequences resembling 5S ribosomal RNA (rRNA) to perform as functional 5S rRNAs in vivo in the Escherichia coli cellular environment was devised previously. Presumably, the only 5S rRNA sequences that would have been fixed by ancestral populations are ones that were functionally valid, and hence the actual historical paths taken through RNA sequence space during 5S rRNA evolution would have most likely utilized valid sequences. Herein, we examine the potential validity of all sequence intermediates along alternative equally parsimonious trajectories through RNA sequence space which connect two pairs of sequences that had previously been shown to behave as valid 5S rRNAs in E. coli. The first trajectory requires a total of four changes. The 14 sequence intermediates provide 24 apparently equally parsimonious paths by which the transition could occur. The second trajectory involves three changes, six intermediate sequences, and six potentially equally parsimonious paths. In total, only eight of the 20 sequence intermediates were found to be clearly invalid. As a consequence of the position of these invalid intermediates in the sequence space, seven of the 30 possible paths consisted of exclusively valid sequences. In several cases, the apparent validity/invalidity of the intermediate sequences could not be anticipated on the basis of current knowledge of the 5S rRNA structure. This suggests that the interdependencies in RNA sequence space may be more complex than currently appreciated. If ancestral sequences predicted by parsimony are to be regarded as actual historical sequences, then the present results would suggest that they should also satisfy a validity requirement and that, in at least limited cases, this conjecture can be tested experimentally.

NASA Discipline Exobiology↗
DOE OSTI2025

Long-read RNA sequencing atlas of human microglia isoforms elucidates disease-associated genetic regulation of splicing

Microglia, the innate immune cells of the central nervous system, have been genetically implicated in multiple neurodegenerative diseases. Mapping the genetics of gene expression in human microglia has identified several loci associated with disease-associated genetic variants in microglia-specific regulatory elements. However, identifying genetic effects on splicing is challenging because of the use of short sequencing reads. Here, we present the isoform-centric microglia genomic atlas (isoMiGA), which leverages long-read RNA sequencing to identify 35,879 novel microglia isoforms. We show that these isoforms are involved in stimulation response and brain region specificity. We then quantified the expression of both known and novel isoforms in a multi-ancestry meta-analysis of 555 human microglia short-read RNA sequencing samples from 391 donors, and found associations with genetic risk loci in Alzheimer’s and Parkinson’s disease. We nominate several loci that may act through complex changes in isoform and splice-site usage.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Analysis of Transcriptome and Epitranscriptome in Plants Using PacBio Iso-Seq and Nanopore-Based Direct RNA Sequencing

Nanopore sequencing from Oxford Nanopore Technologies (ONT) and Pacific BioSciences (PacBio) single-molecule real-time (SMRT) long-read isoform sequencing (Iso-Seq) are revolutionizing the way transcriptomes are analyzed. These methods offer many advantages over most widely used high-throughput short-read RNA sequencing (RNA-Seq) approaches and allow a comprehensive analysis of transcriptomes in identifying full-length splice isoforms and several other post-transcriptional events. In addition, direct RNA-Seq provides valuable information about RNA modifications, which are lost during the PCR amplification step in other methods. Here, we present a comprehensive summary of important applications of these technologies in plants, including identification of complex alternative splicing (AS), full-length splice variants, fusion transcripts, and alternative polyadenylation (APA) events. Furthermore, we discuss the impact of the newly developed nanopore direct RNA-Seq in advancing epitranscriptome research in plants. Additionally, we summarize computational tools for identifying and quantifying full-length isoforms and other co/post-transcriptional events and discussed some of the limitations with these methods. Sequencing of transcriptomes using these new single-molecule long-read methods will unravel many aspects of transcriptome complexity in unprecedented ways as compared to previous short-read sequencing approaches. Analysis of plant transcriptomes with these new powerful methods that require minimum sample processing is likely to become the norm and is expected to uncover novel co/post-transcriptional gene regulatory mechanisms that control biological outcomes during plant development and in response to various stresses.

alternative polyadenylation↗
DOE OSTI2019

A comparison between whole transcript and 3’ RNA sequencing methods using Kapa and Lexogen library preparation methods

Background: 3’ RNA sequencing provides an alternative to whole transcript analysis. However, we do not know a priori the relative advantage of each method. Thus, a comprehensive comparison between the whole transcript and the 3' method is needed to determine their relative merits. To this end, we used two commercially available library preparation kits, the KAPA Stranded mRNA-Seq kit (traditional method) and the Lexogen QuantSeq 3’ mRNASeq kit (3' method), to prepare libraries from mouse liver RNA. We then sequenced and analyzed the libraries to determine the advantages and disadvantages of these two approaches. Results: We found that the traditional whole transcript method and the 3’ RNA-Seq method had similar levels of reproducibility. As expected, the whole transcript method assigned more reads to longer transcripts, while the 3' method assigned roughly equal numbers of reads to transcripts regardless of their lengths. We found that the 3’ RNA-Seq method detected more short transcripts than the whole transcript method. With regard to differential expression analysis, we found that the whole transcript method detected more differentially expressed genes, regardless of the level of sequencing depth. Conclusions: The 3’ RNA-Seq method was better able to detect short transcripts, while the whole transcript RNA-Seq was able to detect more differentially expressed.

3’ RNA-Seq↗
DOE OSTI2018

Global regulatory function of the low oxygen-induced transcriptional regulator LoiA in Salmonella Typhimurium revealed by RNA sequencing

Salmonella enterica serovar Typhimurium (S. Typhimurium) is a major intestinal pathogen that can infect both humans and a variety of animals. LoiA, a novel virulence-regulating protein encoded in Salmonella pathogenicity island (SPI)-14, has been shown to be induced under low oxygen conditions and contribute to S. Typhimurium invasion into intestinal epithetical cells by activating the SPI-1 invasion genes. However, the global regulatory network of LoiA remains unknown. Here, we used high-throughput RNA sequencing (RNA-seq) technology to investigate the regulatory function of LoiA in S. Typhimurium under low oxygen conditions. A total of 1250 genes were differentially expressed between the loiA mutant and the wild-type strain; 413 genes were up-regulated and 837 were down-regulated. SPI-1 gene expression was down-regulated in the loiA mutant, consistent with previous results. SPI-2 gene expression was not affected by deletion of loiA; the expression of most genes involved in flagellar basal body and hook biosynthesis was up-regulated in the loiA mutant, while the expression of genes associated with flagellin, motility, and chemotaxis was down-regulated; the expression of lon, encoding an ATP-dependent protease, was up-regulated in the mutant. This study indicates that LoiA regulates a variety of virulence-associated genes in S. Typhimurium. The negative regulation of Lon protease by LoiA indicates that LoiA can regulates several virulence-associated genes in S. Typhimurium via the Lon protease.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Optimal Bayesian supervised domain adaptation for RNA sequencing data

Abstract Motivation When learning to subtype complex disease based on next-generation sequencing data, the amount of available data is often limited. Recent works have tried to leverage data from other domains to design better predictors in the target domain of interest with varying degrees of success. But they are either limited to the cases requiring the outcome label correspondence across domains or cannot leverage the label information at all. Moreover, the existing methods cannot usually benefit from other information available a priori such as gene interaction networks. Results In this article, we develop a generative optimal Bayesian supervised domain adaptation (OBSDA) model that can integrate RNA sequencing (RNA-Seq) data from different domains along with their labels for improving prediction accuracy in the target domain. Our model can be applied in cases where different domains share the same labels or have different ones. OBSDA is based on a hierarchical Bayesian negative binomial model with parameter factorization, for which the optimal predictor can be derived by marginalization of likelihood over the posterior of the parameters. We first provide an efficient Gibbs sampler for parameter inference in OBSDA. Then, we leverage the gene-gene network prior information and construct an informed and flexible variational family to infer the posterior distributions of model parameters. Comprehensive experiments on real-world RNA-Seq data demonstrate the superior performance of OBSDA, in terms of accuracy in identifying cancer subtypes by utilizing data from different domains. Moreover, we show that by taking advantage of the prior network information we can further improve the performance. Availability and implementation The source code for implementations of OBSDA and SI-OBSDA are available at the following link. https://github.com/SHBLK/BSDA. Supplementary information Supplementary data are available at Bioinformatics online.

Biochemistry & Molecular Biology↗
NASA NTRS1985

Phylogenetic origins of the plant mitochondrion based on a comparative analysis of 5S ribosomal RNA sequences

The complete nucleotide sequences of 5S ribosomal RNAs from Rhodocyclus gelatinosa, Rhodobacter sphaeroides, and Pseudomonas cepacia were determined. Comparisons of these 5S RNA sequences show that rather than being phylogenetically related to one another, the two photosynthetic bacterial 5S RNAs share more sequence and signature homology with the RNAs of two nonphotosynthetic strains. Rhodobacter sphaeroides is specifically related to Paracoccus denitrificans and Rc. gelatinosa is related to Ps. cepacia. These results support earlier 16S ribosomal RNA studies and add two important groups to the 5S RNA data base. Unique 5S RNA structural features previously found in P. denitrificans are present also in the 5S RNA of Rb. sphaeroides; these provide the basis for subdivisional signatures. The immediate consequence of obtaining these new sequences is that it is possible to clarify the phylogenetic origins of the plant mitochondrion. In particular, a close phylogenetic relationship is found between the plant mitochondria and members of the alpha subdivision of the purple photosynthetic bacteria, namely, Rb. sphaeroides, P. denitrificans, and Rhodospirillum rubrum.

Villanueva, E.↗
DOE OSTI2018

RNAcentral: a hub of information for non-coding RNA sequences

RNAcentral is a comprehensive database of non-coding RNA (ncRNA) sequences, collating information on ncRNA sequences of all types from a broad range of organisms. We have recently added a new genome mapping pipeline that identifies genomic locations for ncRNA sequences in 296 species. We have also added several new types of functional annotations, such as tRNA secondary structures, Gene Ontology annotations, and miRNA-target interactions. A new quality control mechanism based on Rfam family assignments identifies potential contamination, incomplete sequences, and more. The RNAcentral database has become a vital component of many workflows in the RNA community, serving as both the primary source of sequence data for academic and commercial groups, as well as a source of stable accessions for the annotation of genomic and functional features. These examples are facilitated by an improved RNAcentral web interface, which features an updated genome browser, a new sequence feature viewer, and improved text search functionality.

Sweeney, Blake A.↗
NASA NTRS2002

Identification of characteristic oligonucleotides in the bacterial 16S ribosomal RNA sequence dataset

MOTIVATION: The phylogenetic structure of the bacterial world has been intensively studied by comparing sequences of 16S ribosomal RNA (16S rRNA). This database of sequences is now widely used to design probes for the detection of specific bacteria or groups of bacteria one at a time. The success of such methods reflects the fact that there are local sequence segments that are highly characteristic of particular organisms or groups of organisms. It is not clear, however, the extent to which such signature sequences exist in the 16S rRNA dataset. A better understanding of the numbers and distribution of highly informative oligonucleotide sequences may facilitate the design of hybridization arrays that can characterize the phylogenetic position of an unknown organism or serve as the basis for the development of novel approaches for use in bacterial identification. RESULTS: A computer-based algorithm that characterizes the extent to which any individual oligonucleotide sequence in 16S rRNA is characteristic of any particular bacterial grouping was developed. A measure of signature quality, Q(s), was formulated and subsequently calculated for every individual oligonucleotide sequence in the size range of 5-11 nucleotides and for 15mers with reference to each cluster and subcluster in a 929 organism representative phylogenetic tree. Subsequently, the perfect signature sequences were compared to the full set of 7322 sequences to see how common false positives were. The work completed here establishes beyond any doubt that highly characteristic oligonucleotides exist in the bacterial 16S rRNA sequence dataset in large numbers. Over 16,000 15mers were identified that might be useful as signatures. Signature oligonucleotides are available for over 80% of the nodes in the representative tree.

NASA Discipline Life Sciences Technologies↗
DOE OSTI2021

Single-cell RNA sequencing of batch Chlamydomonas cultures reveals heterogeneity in their diurnal cycle phase

The photosynthetic unicellular alga Chlamydomonas (Chlamydomonas reinhardtii) is a versatile reference for algal biology because of its ease of culture in the laboratory. Genomic and systems biology approaches have previously described transcriptome responses to environmental changes using bulk data, thus representing the average behavior from pools of cells. Here, we apply single-cell RNA sequencing (scRNA-seq) to probe the heterogeneity of Chlamydomonas cell populations under three environments and in two genotypes differing by the presence of a cell wall. First, we determined that RNA can be extracted from single algal cells with or without a cell wall, offering the possibility to sample natural algal communities. Second, scRNA-seq successfully separated single cells into nonoverlapping cell clusters according to their growth conditions. Cells exposed to iron or nitrogen deficiency were easily distinguished despite a shared tendency to arrest photosynthesis and cell division to economize resources. Notably, these groups of cells not only recapitulated known patterns observed with bulk RNA-seq but also revealed their inherent heterogeneity. A substantial source of variation between cells originated from their endogenous diurnal phase, although cultures were grown in constant light. We exploited this result to show that circadian iron responses may be conserved from algae to land plants. We document experimentally that bulk RNA-seq data represent an average of typically hidden heterogeneity in the population.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

A robust method of nuclei isolation for single-cell RNA sequencing of solid tissues from the plant genus Populus

Single-cell transcriptome analysis has been extensively applied in humans and animal models to uncover gene expression heterogeneity between the different cell types of a tissue or an organ. It demonstrated its capability to discover key regulatory elements that determine cell fate during developmental programs. Single-cell analysis requires the isolation and labeling of the messenger RNA (mRNA) derived from each cell. These challenges were primarily addressed in mammals by developing microfluidic-based approaches. For plant species whose cells contain cell walls, these approaches have generally required the generation of isolated protoplasts. Many plant tissues’ secondary cell wall hinders enzymatic digestion required for individual protoplast isolation, resulting in an unequal representation of cell types in a protoplast population. This limitation is especially critical for cell types located in the inner layers of a tissue or the inner tissues of an organ. Consequently, single-cell RNA sequencing (scRNA-seq) studies using microfluidic approaches in plants have mainly been restricted to Arabidopsis roots, for which well-established procedures of protoplast isolation are available. Here we present a simple alternative approach to generating high-quality protoplasts from plant tissue by characterizing the mRNA extracted from individual nuclei instead of whole cells. We developed the protocol using two different plant materials with varying cellular complexity levels and cell wall structure, Populus shoot apices, and more lignified stems. Using the 10× Genomics Chromium technology, we show that this procedure results in intact mRNA isolation and limited leakage, with a broad representation of individual cell transcriptomes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Single cell RNA sequencing reveals shifts in cell maturity and function of endogenous and infiltrating cell types in response to acute intervertebral disc injury

Intervertebral disc (IVD) degeneration contributes to disabling back pain. Degeneration can be initiated by injury and progressively leads to an irreversible loss of cells and function. IVD function restoration through cell replacement therapies have had limited success due to knowledge gaps in the critical cell populations important for repair. Here, in this study, we used single cell RNA sequencing to identify the transcriptional changes of IVD resident and infiltrating cell populations from Control and Injured coccygeal IVDs extracted from 12-week-old female C57BL/6J mice 7 days post injury. Clustering, gene ontology, and pseudotime trajectory analyses determined transcriptomic divergences with injury, flow cytometry identified they types of infiltrating immune cells, and immunofluorescence was utilized to define mesenchymal stem cell (MSC) localization. We identified 11 distinct clusters that included IVD, immune, vascular cells, and MSCs. Differential gene expression analysis determined that Outer Annulus Fibrosus, Neutrophils, Saa2-High MSCs, Macrophages, and Krt18 + Nucleus Pulposus (NP) cells were the major drivers of transcriptomic differences between Control and Injured cells. Gene ontology revealed that the most upregulated biological pathways were angiogenesis and T cell-related while wound healing and ECM regulation were downregulated. Pseudotime trajectory analyses revealed that IVD injury directed cells towards increased differentiation in all clusters, except for Krt18 + NP cells which remained in a less mature cell state. Saa2-High and Grem1-High MSCs populations shifted towards more differentiated IVD cells profiles with injury and localized distinctly within the IVD. This study revealed novel MSC populations with the potential to be leveraged for future IVD repair studies.

Cartilage↗
NASA NTRS2019

A Standard RNA Sequencing Assay for Space Biology

Given the limited opportunities for biological experimentation in space, it is often desirable to compare results across experiments to gain additional insights into the effects of spaceflight on biological systems. However, this approach is made difficult by a multitude of confounding factors including differences in strain, hardware configuration, and sample processing. To help harmonize datasets, the NASA GeneLab Project has developed consistent sample and data processing protocols for the generation of raw and processed RNA-sequencing data from various mouse tissues. We will present these and discuss how they can be used to make novel discoveries from these precious samples.

Galazka, Jonathan↗
DOE OSTI2024

Dual-RNA-sequencing to elucidate the interactions between sorghum and Colletotrichum sublineola

In warm and humid regions, the productivity of sorghum is significantly limited by the fungal hemibiotrophic pathogen Colletotrichum sublineola , the causal agent of anthracnose, a problematic disease of sorghum ( Sorghum bicolor (L.) Moench) that can result in grain and biomass yield losses of up to 50%. Despite available genomic resources of both the host and fungal pathogen, the molecular basis of sorghum− C. sublineola interactions are poorly understood. By employing a dual-RNA sequencing approach, the molecular crosstalk between sorghum and C. sublineola can be elucidated. In this study, we examined the transcriptomes of four resistant sorghum accessions from the sorghum association panel (SAP) at varying time points post-infection with C. sublineola . Approximately 0.3% and 93% of the reads mapped to the genomes of C. sublineola and Sorghum bicolor , respectively. Expression profiling of in vitro versus in planta C. sublineola at 1-, 3-, and 5-days post-infection (dpi) indicated that genes encoding secreted candidate effectors, carbohydrate-active enzymes (CAZymes), and membrane transporters increased in expression during the transition from the biotrophic to the necrotrophic phase (3 dpi). The hallmark of the pathogen-associated molecular pattern (PAMP)-triggered immunity in sorghum includes the production of reactive oxygen species (ROS) and phytoalexins. The majority of effector candidates secreted by C. sublineola were predicted to be localized in the host apoplast, where they could interfere with the PAMP-triggered immunity response, specifically in the host ROS signaling pathway. The genes encoding critical molecular factors influencing pathogenicity identified in this study are a useful resource for subsequent genetic experiments aimed at validating their contributions to pathogen virulence. This comprehensive study not only provides a better understanding of the biology of C. sublineola but also supports the long-term goal of developing resistant sorghum cultivars.

Vela, Saddie↗
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 NTRS1987

Ribosomal RNA sequence suggest microsporidia are extremely ancient eukaryotes

A comparative sequence analysis of the 18S small subunit ribosomal RNA (rRNA) of the microsporidium Vairimorpha necatrix is presented. The results show that this rRNA sequence is more unlike those of other eukaryotes than any known eukaryote rRNA sequence. It is concluded that the lineage leading to microsporidia branched very early from that leading to other eukaryotes.

Vossbrinck, C. R.↗