Search NASA⌕ Search

SEARCH · Search NASA

Results for “Sequence Analysis, RNA”

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 91 records · Page 5

Transcriptomic and functional analyses uncover a conserved effector driving genotype-dependent virulence in the Sphaerulina musiva-Populus trichocarpa interaction

The introduction of invasive microbes compromises the structure, biodiversity, and function of naïve ecosystems. Sphaerulina musiva, a hemibiotrophic pathogen that causes leaf spot and stem cankers in Populus species, exemplifies an invasive fungal pathogen spread by human activities. However, the genetic mechanisms of pathogenicity and virulence are poorly understood, impeding mitigation strategies. We utilized RNA sequencing to identify fungal effectors linked to stem canker formation, informing the development of future strategies for effective disease management. Our analysis revealed 70 genes differentially expressed at 2 weeks and 110 genes at 3 weeks between inoculated trees and controls. Notably, the gene with the highest expression at 2 weeks and the second highest at 3 weeks was homologous to Extracellular protein 2 (Ecp2). Complementary genome-wide association studies linked sequence polymorphisms in this locus to phenotypic variation in disease severity. Infiltration of S. musiva Ecp2 into Populus trichocarpa leaves induced necrosis in susceptible genotypes. Gene disruption using a CRISPR-Cas9 RNP system resulted in a genotype-dependent reduction of stem canker and disease severity. Tracing the evolutionary history of this effector across the fungal kingdom, we uncovered clade-specific gene-family expansions and orthologs in new species. These findings raise questions about the function and adaptive significance of these gene families in fungal lifestyles. Our study provides the first tractable target for breeding resistant poplar genotypes, addressing the challenges of managing S. musiva and uncovering mechanisms that drive its virulence, and provides deeper insights into the evolutionary dynamics of a conserved small-secreted protein with a diversity of functions.

Sondreli, Kelsey L [Oregon State University]↗

Combined analysis of microbial metagenomic and metatranscriptomic sequencing data to assess in situ physiological conditions in the premature infant gut

Microbes alter their transcriptomic profiles in response to the environment. The physiological conditions experienced by a microbial community can thus be inferred using meta-transcriptomic sequencing by comparing transcription levels of specifically chosen genes. However, this analysis requires accurate reference genomes to identify the specific genes from which RNA reads originate. In addition, such an analysis should avoid biases in transcript counts related to differences in organism abundance. In this study we describe an approach to address these difficulties. Sample-specific meta-genomic assembled genomes (MAGs) were used as reference genomes to accurately identify the origin of RNA reads, and transcript ratios of genes with opposite transcription responses were compared to eliminate biases related to differences in organismal abundance, an approach hereafter named the “diametric ratio” method. We used this approach to probe the environmental conditions experienced by Escherichia spp. in the gut of 4 premature infants, 2 of whom developed necrotizing enterocolitis (NEC), a severe inflammatory intestinal disease. We analyzed twenty fecal samples taken from four premature infants (4–6 time points from each infant), and found significantly higher diametric ratios of genes associated with low oxygen levels in samples of infants later diagnosed with NEC than in samples without NEC. We also show this method can be used for examining other physiological conditions, such as exposure to nitric oxide and osmotic pressure. These study results should be treated with caution, due to the presence of confounding factors that might also distinguish between NEC and control infants. Nevertheless, together with benchmarking analyses, we show here that the diametric ratio approach can be applied for evaluating the physiological conditions experienced by microbes in situ. Results from similar studies can be further applied for designing diagnostic methods to detect NEC in its early developmental stages.

59 BASIC BIOLOGICAL SCIENCES↗

Genomics Characterization of an Engineered Corynebacterium glutamicum in Bioreactor Cultivation Under Ionic Liquid Stress

Corynebacterium glutamicum is an ideal microbial chassis for production of valuable bioproducts including amino acids and next generation biofuels. Here we resequence engineered isopentenol (IP) producing C. glutamicum BRC-JBEI 1.1.2 strain and assess differential transcriptional profiles using RNA sequencing under industrially relevant conditions including scale transition and compare the presence vs absence of an ionic liquid, cholinium lysinate ([Ch][Lys]). Analysis of the scale transition from shake flask to bioreactor with transcriptomics identified a distinct pattern of metabolic and regulatory responses needed for growth in this industrial format. These differential changes in gene expression corroborate altered accumulation of organic acids and bioproducts, including succinate, acetate, and acetoin that occur when cells are grown in the presence of 50 mM [Ch][Lys] in the stirred-tank reactor. This new genome assembly and differential expression analysis of cells grown in a stirred tank bioreactor clarify the cell response of an C. glutamicum strain engineered to produce IP.

59 BASIC BIOLOGICAL SCIENCES↗

Metastable multimeric G-quadruplex 2′FY-RNA aptamers that selectively bind pyoverdines

Two 2′FY-RNA aptamers with distinct sequences were selected for specific binding to pyoverdine-Pf5 (PVD-Pf5), increasing chromophore fluorescence upon binding. They also recognized the peptide portion of pyoverdines, as shown by their differential specificity for related variants. Computational analysis and experimental data (NMM binding, CD spectra) identified G-quadruplex structures that were thermally metastable but reformed in the presence of PVD-Pf5. Further structural studies mainly with one aptamer revealed imino proton peaks in 1D H-NMR and pressure stability up to 2 kbar. Electrophoretic evidence identified dimeric G-quadruplexes formed by the 2′FY-RNA aptamers and their RNA equivalents. While cations were necessary for PVD-Pf5 binding, they were not required for G-quadruplex formation. Given the established role of G-quadruplexes as protein interaction sites, multimeric G-quadruplexes offer a potential framework for structure-based regulatory mechanisms in cellular RNAs. In addition to previously characterized multimeric G-quadruplexes, these aptamers contribute novel sequences that expand the repertoire of known multimeric G-quadruplexes.

2′FY-RNA↗

Learning perturbation-inducible cell states from observability analysis of transcriptome dynamics

Abstract A major challenge in biotechnology and biomanufacturing is the identification of a set of biomarkers for perturbations and metabolites of interest. Here, we develop a data-driven, transcriptome-wide approach to rank perturbation-inducible genes from time-series RNA sequencing data for the discovery of analyte-responsive promoters. This provides a set of biomarkers that act as a proxy for the transcriptional state referred to as cell state. We construct low-dimensional models of gene expression dynamics and rank genes by their ability to capture the perturbation-specific cell state using a novel observability analysis. Using this ranking, we extract 15 analyte-responsive promoters for the organophosphate malathion in the underutilized host organism Pseudomonas fluorescens SBW25. We develop synthetic genetic reporters from each analyte-responsive promoter and characterize their response to malathion. Furthermore, we enhance malathion reporting through the aggregation of the response of individual reporters with a synthetic consortium approach, and we exemplify the library’s ability to be useful outside the lab by detecting malathion in the environment. The engineered host cell, a living malathion sensor, can be optimized for use in environmental diagnostics while the developed machine learning tool can be applied to discover perturbation-inducible gene expression systems in the compendium of host organisms.

59 BASIC BIOLOGICAL SCIENCES↗

Deep structural clustering for single-cell RNA-seq data jointly through autoencoder and graph neural network

Abstract Single-cell RNA sequencing (scRNA-seq) permits researchers to study the complex mechanisms of cell heterogeneity and diversity. Unsupervised clustering is of central importance for the analysis of the scRNA-seq data, as it can be used to identify putative cell types. However, due to noise impacts, high dimensionality and pervasive dropout events, clustering analysis of scRNA-seq data remains a computational challenge. Here, we propose a new deep structural clustering method for scRNA-seq data, named scDSC, which integrate the structural information into deep clustering of single cells. The proposed scDSC consists of a Zero-Inflated Negative Binomial (ZINB) model-based autoencoder, a graph neural network (GNN) module and a mutual-supervised module. To learn the data representation from the sparse and zero-inflated scRNA-seq data, we add a ZINB model to the basic autoencoder. The GNN module is introduced to capture the structural information among cells. By joining the ZINB-based autoencoder with the GNN module, the model transfers the data representation learned by autoencoder to the corresponding GNN layer. Furthermore, we adopt a mutual supervised strategy to unify these two different deep neural architectures and to guide the clustering task. Extensive experimental results on six real scRNA-seq datasets demonstrate that scDSC outperforms state-of-the-art methods in terms of clustering accuracy and scalability. Our method scDSC is implemented in Python using the Pytorch machine-learning library, and it is freely available at https://github.com/DHUDBlab/scDSC.

Gan, Yanglan↗

Machine Learning of All Mycobacterium tuberculosis H37Rv RNA-seq Data Reveals a Structured Interplay between Metabolism, Stress Response, and Infection

Mycobacterium tuberculosis is one of the most consequential human bacterial pathogens, posing a serious challenge to 21st century medicine. A key feature of its pathogenicity is its ability to adapt its transcriptional response to environmental stresses through its transcriptional regulatory network (TRN). While many studies have sought to characterize specific portions of the M. tuberculosis TRN, and some studies have performed system-level analysis, few have been able to provide a network-based model of the TRN that also provides the relative shifts in transcriptional regulator activity triggered by changing environments. Here, we compiled a compendium of nearly 650 publicly available, high quality M. tuberculosis RNA-sequencing data sets and applied an unsupervised machine learning method to obtain a quantitative, top-down TRN. It consists of 80 independently modulated gene sets known as “iModulons,” 41 of which correspond to known regulons. These iModulons explain 61% of the variance in the organism’s transcriptional response. We show that iModulons (i) reveal the function of poorly characterized regulons, (ii) describe the transcriptional shifts that occur during environmental changes such as shifting carbon sources, oxidative stress, and infection events, and (iii) identify intrinsic clusters of regulons that link several important metabolic systems, including lipid, cholesterol, and sulfur metabolism. This transcriptome-wide analysis of the M. tuberculosis TRN informs future research on effective ways to study and manipulate its transcriptional regulation and presents a knowledge-enhanced database of all published high-quality RNA-seq data for this organism to date.

59 BASIC BIOLOGICAL SCIENCES↗

A genetic toolbox for metabolic engineering of Issatchenkia orientalis

The nonconventional yeast Issatchenkia orientalis can grow under highly acidic conditions and has been explored for production of various organic acids. Yet, its broader application is hampered by the lack of efficient genetic tools to enable sophisticated metabolic manipulations. We recently constructed an episomal plasmid based on the autonomously replicating sequence (ARS) from Saccharomyces cerevisiae (ScARS) in I. orientalis and developed a CRISPR/Cas9 system for multiplexed gene deletions. Here we report three additional genetic tools including: (1) identification of a 0.8 kb centromere-like (CEN-L) sequence from the I. orientalis genome by using bioinformatics and functional screening; (2) discovery and characterization of a set of constitutive promoters and terminators under different culture conditions by using RNA-Seq analysis and a fluorescent reporter; and (3) development of a rapid and efficient in vivo DNA assembly method in I. orientalis, which exhibited ~100% fidelity when assembling a 7 kb-plasmid from seven DNA fragments ranging from 0.7 kb to 1.7 kb. As evidence of concept, we used these genetic tools to rapidly construct a functional xylose utilization pathway in I. orientalis.

59 BASIC BIOLOGICAL SCIENCES↗

EvoNet: A phylogenomic and systems biology approach to identify genes underlying plant survival in marginal, low‐N soils

The DOE‐BER “EvoNet” project investigates the genetic and molecular basis of plant resilience in extreme environments. We do this by identifying key genes that enable “extreme survivor” species to thrive in the nitrogen-poor soils of Chile’s hyper-arid Atacama Desert. Our collections focus on 32 Atacama extremophile species, including seven grass species with potential biofuel applications. To identify genes-of-importance to survival we compared genomic and transcriptomic profiles of extremophile species that thrive in the Atacama to those of closely related “sister” species from nitrogen-rich arid and mesic regions of California. Deep RNA sequencing and de novo transcriptome assembly across these triplet species sets supported a phylogenomic framework for identifying positively selected genes associated with adaptive divergence. Our integrative analysis combined ecological and environmental data, metagenomics, evolutionary and systems biology, and metabolomics. This enabled us to create an unprecedented framework for systematically understanding how non-model plants have adapted to survive in extreme conditions. Our resulting database of positively selected ortholog groups in the extremophile plants offers promising targets for engineering crop and biofuel species with enhanced resilience to drought and extreme weather. Additionally, our newest dataset explores and exploits a complementary metabolomic approach. This new aspect provides innovative strategies to manipulate plant cell metabolism, further supporting efforts to improve agricultural productivity in the face of extreme climates. Importantly, our combined evolutionary- and metabolomic-based strategies focused on convergent patterns of adaptation, providing a genetic and metabolomic toolkit for improving crop and biofuel resilience across diverse plant species. Finally, our novel exploration of ecological and evolutionary dynamics delivered to the community a phylogenomic computational pipeline called “PhyloGeneious.” Our continued adaptations of this pipeline are publicly available to expedite evolutionary genomic research for future scientific discoveries. In total, our DOE-BER has provided genomic, metabolomic, and computational strategies to understand how extremophile plants provide evolutionary and physiological targets for improving agricultural and biofuel production.

59 BASIC BIOLOGICAL SCIENCES↗

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

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

60 APPLIED LIFE SCIENCES↗

Space flight associated changes in astronauts’ plasma-derived small extracellular vesicle microRNA: Biomarker identification

This pilot study suggests relatively short (median 12 days long) low-Earth orbit (LEO) spaceflight induces changes in circulating plasma small extracellular vesicle (sEV) microRNA expression. Normalization of small RNA sequencing (sRNAseq) data and quantitative polymerase chain reaction (qPCR) validation confirmed miR-4732-3p is significantly upregulated up to 3 days post-landing, and enrichment analysis suggests this miRNA is expressed in various central nervous system tissues and hematopoietic cells and may be linked to different organ disorders.

60 APPLIED LIFE SCIENCES↗

Identifying strengths and weaknesses of methods for computational network inference from single-cell RNA-seq data

Single-cell RNA-sequencing (scRNA-seq) offers unparalleled insight into the transcriptional programs of different cellular states by measuring the transcriptome of thousands of individual cells. An emerging problem in the analysis of scRNA-seq is the inference of transcriptional gene regulatory networks and a number of methods with different learning frameworks have been developed to address this problem. Here, we present an expanded benchmarking study of eleven recent network inference methods on seven published scRNA-seq datasets in human, mouse, and yeast considering different types of gold standard networks and evaluation metrics. We evaluate methods based on their computing requirements as well as on their ability to recover the network structure. We find that, while most methods have a modest recovery of experimentally derived interactions based on global metrics such as Area Under the Precision Recall curve, methods are able to capture targets of regulators that are relevant to the system under study. Among the top performing methods that use only expression were SCENIC, PIDC, MERLIN or Correlation. Addition of prior biological knowledge and the estimation of transcription factor activities resulted in the best overall performance with the Inferelator and MERLIN methods that use prior knowledge outperforming methods that use expression alone. We found that imputation for network inference did not improve network inference accuracy and could be detrimental. Comparisons of inferred networks for comparable bulk conditions showed that the networks inferred from scRNA-seq datasets are often better or at par with the networks inferred from bulk datasets. Our analysis should be beneficial in selecting methods for network inference. At the same time, this highlights the need for improved methods and better gold standards for regulatory network inference from scRNAseq datasets.

59 BASIC BIOLOGICAL SCIENCES↗

CABO-16S—a Combined Archaea, Bacteria, Organelle 16S rRNA database framework for amplicon analysis of prokaryotes and eukaryotes in environmental samples

Abstract Identification of both prokaryotic and eukaryotic microorganisms in environmental samples is currently challenged by the need for additional sequencing to obtain separate 16S and 18S ribosomal RNA (rRNA) amplicons or the constraints imposed by “universal” primers. Organellar 16S rRNA sequences are amplified and sequenced along with prokaryote 16S rRNA and provide an alternative method to identify eukaryotic microorganisms. CABO-16S combines bacterial and archaeal sequences from the SILVA database with 16S rRNA sequences of plastids and other organelles from the PR2 database to enable identification of all 16S rRNA sequences. Comparison of CABO-16S with SILVA 138.2 results in equivalent taxonomic classification of mock communities and increased classification of diverse environmental samples. In particular, identification of phototrophic eukaryotes in shallow seagrass environments, marine waters, and lake waters was increased. The CABO-16S framework allows users to add custom sequences for further classification of underrepresented clades and can be easily updated with future releases of reference databases. Addition of sequences obtained from Sanger sequencing of methane seep sediments and curated sequences of the polyphyletic SEEP-SRB1 clade resulted in differentiation of syntrophic and non-syntrophic SEEP-SRB1 in hydrothermal vent sediments. CABO-16S highlights the benefit of combining and amending existing training sets when studying microorganisms in diverse environments.

Eitel, Eryn M. (ORCID:0009000723919297)↗

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

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

3’ RNA-seq↗

Armillaria root rot fungi host single-stranded RNA viruses

Species of Armillaria are distributed globally and include some of the most important pathogens of forest and ornamental trees. Some of them form large long-living clones that are considered as one of the largest organisms on earth and are capable of long-range spore-mediated transfer as well as vegetative spread by drought-resistant hyphal cords called rhizomorphs. However, the virus community infecting these species has remained unknown. In this study we used dsRNA screening and high-throughput sequencing to search for possible virus infections in a collection of Armillaria isolates representing three different species: Armillaria mellea from South Africa, A. borealis from Finland and Russia (Siberia) and A. cepistipes from Finland. Our analysis revealed the presence of both negative-sense RNA viruses and positive-sense RNA viruses, while no dsRNA viruses were detected. The viruses included putative new members of virus families Mymonaviridae, Botourmiaviridae and Virgaviridae and members of a recently discovered virus group tentatively named “ambiviruses” with ambisense bicistronic genomic organization. We demonstrated that Armillaria isolates can be cured of viruses by thermal treatment, which enables the examination of virus effects on host growth and phenotype using isogenic virus-infected and virus-free strains.

59 BASIC BIOLOGICAL SCIENCES↗

Quantitative trait locus (QTL) mapping and transcriptome profiling identify QTLs and candidate genes associated with heat stress response during reproductive development in Camelina sativa

Camelina sativa (L.) Crantz is a low-input oilseed crop that has great potential in providing sustainable feedstock for biofuels and bioproducts. Climate change is threatening production of camelina with rising global temperatures. Elucidating the genetic response to high temperatures is essential for successful breeding of heat-tolerant camelina varieties. Here, we report a combinatorial approach to identifying candidate genes associated with heat stress by quantitative trait locus (QTL) mapping and comparative transcriptome profiling. A population of recombinant inbred lines (RILs) was grown in a controlled growth chamber under the high-temperature regimes for 14 days beginning at the onset of the reproductive stage. Several traits related to seed production were evaluated at maturity. The QTL analysis identified several regions with co-located traits on chromosomes 8, 10, and 12. Two RILs with contrasting phenotypic responses to heat stress were chosen for gene expression profiling via RNA sequencing. Multiple pathways and genes were found to be strongly affected by heat stress, and many genes expressed differently between the two RILs. Several genes identified within the QTL regions were considered strong candidates that may control heat tolerance during reproduction in camelina. These studies provide resources for future studies that may assist in improving the heat tolerance of camelina.

60 APPLIED LIFE SCIENCES↗

Metatranscriptomics analysis reveals a novel transcriptional and translational landscape during Middle East respiratory syndrome coronavirus infection

Among all RNA viruses, coronavirus RNA transcription is the most complex and involves a process termed “discontinuous transcription” that results in the production of a set of 3'-nested, co-terminal genomic and subgenomic RNAs during infection. While the expression of the classic canonical set of subgenomic RNAs depends on the recognition of a 6- to 7-nt transcription regulatory core sequence (TRS), here, we use deep sequence and metagenomics analysis strategies and show that the coronavirus transcriptome is even more vast and more complex than previously appreciated and involves the production of leader-containing transcripts that have canonical and noncanonical leader-body junctions. Moreover, by ribosome protection and proteomics analyses, we show that both positive- and negative-sense transcripts are translationally active. The data support the hypothesis that the coronavirus proteome is much vaster than previously noted in the literature.

59 BASIC BIOLOGICAL SCIENCES↗

Data from A Genetic Toolbox for Metabolic Engineering of Issatchenkia orientalis

The nonconventional yeast Issatchenkia orientalis can grow under highly acidic conditions and has been explored for production of various organic acids. However, its broader application is hampered by the lack of efficient genetic tools to enable sophisticated metabolic manipulations. We recently constructed an episomal plasmid based on the autonomously replicating sequence (ARS) from Saccharomyces cerevisiae (ScARS) in I. orientalis and developed a CRISPR/Cas9 system for multiplexed gene deletions. Here we report three additional genetic tools including: (1) identification of a 0.8 kb centromere-like (CEN-L) sequence from the I. orientalis genome by using bioinformatics and functional screening; (2) discovery and characterization of a set of constitutive promoters and terminators under different culture conditions by using RNA-Seq analysis and a fluorescent reporter; and (3) development of a rapid and efficient in vivo DNA assembly method in I. orientalis , which exhibited ~100% fidelity when assembling a 7 kb-plasmid from seven DNA fragments ranging from 0.7 kb to 1.7 kb. As proof of concept, we used these genetic tools to rapidly construct a functional xylose utilization pathway in I. orientalis .

Conversion↗