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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Ex vivo to in vivo model of malignant peripheral nerve sheath tumors for precision oncology

Malignant peripheral nerve sheath tumors (MPNST) are aggressive soft tissue sarcomas that often develop in patients with neurofibromatosis type 1 (NF1). To address the critical need for novel therapeutics in MPNST, we aimed to establish an ex vivo 3D platform that accurately captured the genomic diversity of MPNST and could be utilized in a medium-throughput manner for drug screening studies to be validated in vivo using patient-derived xenografts (PDX). Genomic analysis was performed on all PDX-tumor pairs. Selected PDX were harvested for assembly into 3D microtissues. Based on prior work in our labs, we evaluated drugs (trabectedin, olaparib, and mirdametinib) ex vivo and in vivo. For 3D microtissue studies, cell viability was the endpoint as assessed by Zeiss Axio Observer. For PDX drug studies, tumor volume was measured twice weekly. Bulk RNA sequencing was performed to identify pathways enriched in cells. We developed 13 NF1-associated MPNST-PDX and identified mutations or structural abnormalities in NF1 (100%), SUZ12 (85%), EED (15%), TP53 (15%), CDKN2A (85%), and chromosome 8 gain (77%). We successfully assembled PDX into 3D microtissues, categorized as robust (>90% viability at 48 h), good (>50%), or unusable (<50%). We evaluated drug response to “robust” or “good” microtissues, namely MN-2, JH-2-002, JH-2-079-c, and WU-225. Drug response ex vivo predicted drug response in vivo, and enhanced drug effects were observed in select models. These data support the successful establishment of a novel 3D platform for drug discovery and MPNST biology exploration in a system representative of the human condition.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Phaseolus vulgaris SUT1.1 is a high affinity sucrose–proton co–transporter

Plant sucrose transporters are required for phloem loading, and therefore are essential for plant growth and development. In common beans (Phaseolus vulgaris) there are only two sucrose transporters functionally characterized. Through a previous RNA-seq study, we identified a putative sucrose transporter in common bean, which we hypothesize to function in import of sucrose into plant cells. In silico analysis revealed that PvSUT1.1 is a putative sucrose-proton co-transporter distinct from other characterized sucrose transporters in common bean indicating that this is a previously undescribed transporter protein in beans. Further analysis revealed that PvSUT1.1 shares high protein sequence homology to the phloem loader Arabidopsis SUC2; both have 12 transmembrane domains, a typical characteristic of plant sucrose transporters. Heterologous expression in yeast further showed PvSUT1.1 to be functional and it imported sucrose into yeast cells with a K m of 0.7 mM sucrose. Import of sucrose through PvSUT1.1 is also pH-dependent with highest uptake at pH 4.0, and activity is lost in the presence of the uncoupler carbonyl cyanide 3-chlorophenylhydrazone. Consistent with identification of PvSUT1.1 as a Type I transporter, PvSUT1.1 also transports esculin. Finally, PvSUT1.1 showed expression in multiple tissues and the protein was localized to the plasma membrane. The results show that PvSUT1.1 is a sucrose transporter that is probably involved in the uptake of sucrose into source and sink cells. The potential role of PvSUT1.1 in leaf phloem loading of sucrose in common beans and its importance in heat tolerance of reproductive tissues are further discussed.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Insights of Molecular Mechanism of Xylem Development in Five Black Poplar Cultivars

Black poplar (Populus deltoides, P. nigra, and their hybrids) is the main poplar cultivars in China. It offers interesting options of large-scale biomass production for bioenergy due to its rapid growth and high yield. Poplar wood properties were associated with chemical components and physical structures during wood formation. In this study, five poplar cultivars, P. euramericana ‘Zhonglin46’ (Pe1), P. euramericana ‘Guariento’ (Pe2), P. nigra ‘N179’ (Pn1), P. deltoides ‘Danhong’ (Pd1), and P. deltoides ‘Nanyang’ (Pd2), were used to explore the molecular mechanism of xylem development. We analyzed the structural differences of developing xylem in the five cultivars and profiled the transcriptome-wide gene expression patterns through RNA sequencing. The cross sections of the developing xylem showed that the cell wall thickness of developed fiber in Pd1 was thickest and the number of xylem vessels of Pn1 was the least. A total of 10,331 differentially expressed genes were identified among 10 pairwise comparisons of the five cultivars, most of them were related to programmed cell death and secondary cell wall thickening. K-means cluster analysis and Gene Ontology enrichment analysis showed that the genes highly expressed in Pd1 were related to nucleotide decomposition, metabolic process, transferase, and microtubule cytoskeleton; whereas the genes highly expressed in Pn1 were involved in cell wall macromolecule decomposition and polysaccharide binding processes. Based on a weighted gene co-expression network analysis, a large number of candidate regulators for xylem development were identified. And their potential regulatory roles to cell wall biosynthesis genes were validated by a transient overexpression system. This study provides a set of promising candidate regulators for genetic engineering to improve feedstock and enhance biofuel conversion in the bioenergy crop Populus.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Regional epithelial cell diversity in the small intestine of pigs

Abstract Understanding regional distribution and specialization of small intestinal epithelial cells is crucial for developing methods to control appetite, stress, and nutrient uptake in swine. To establish a better understanding of specific epithelial cells found across different regions of the small intestine in pigs, we utilized single-cell RNA sequencing (scRNA-seq) to recover and analyze epithelial cells from duodenum, jejunum, and ileum. Cells identified included crypt cells, enterocytes, BEST4 enterocytes, goblet cells, and enteroendocrine (EE) cells. EE cells were divided into two subsets based on the level of expression of the EE lineage commitment gene, NEUROD1. NEUROD1hi EE cells had minimal expression of hormone-encoding genes and were dissimilar to EE cells in humans and mice, indicating a subset of EE cells unique to pigs. Recently discovered BEST4 enterocytes were detected in both crypts and villi throughout the small intestine via in situ staining, unlike in humans, where BEST4 enterocytes are found only in small intestinal villi. Proximal-to-distal gradients of expression were noted for hormone-encoding genes in EE cells and nutrient transport genes in enterocytes via scRNA-seq, demonstrating regional specialization. Regional gene expression in EE cells and enterocytes was validated via quantitative PCR (qPCR) analysis of RNA isolated from epithelial cells of different small intestinal locations. Though many genes had similar patterns of regional expression when assessed by qPCR of total epithelial cells, some regional expression was only detected via scRNA-seq, highlighting advantages of scRNA-seq to deconvolute cell type-specific regional gene expression when compared to analysis of bulk samples. Overall, results provide new information on regional localization and transcriptional profiles of epithelial cells in the pig small intestine.

Agriculture↗
DOE OSTI2025

Comparative transcriptomic and phenotypic analysis of monoclonal and polyclonal Populus deltoides genotypes

Populus species are highly valued for bioenergy and bioproducts due to their rapid growth and productivity. Polyclonal plantings, or mixtures of Populus clones, have shown the potential to enhance resource utilization and productivity, likely due to phenotypic differences arising from niche differentiation. In this study, we investigated gene expression and productivity in monoclonal and polyclonal stands of P. deltoides . Phenotypic results showed that polyclonal plots exhibited higher leaf area index (LAI; p < 0.01, 2.96 ± 0.057 m 2 ) and total biomass (p < 0.01, 2.74 ± 0.06) compared to monoclonal plots, indicating superior productivity. RNA sequencing revealed upregulation of key genes such as exocyst subunit exo70 family protein H7 (EXO70H7), NDH-dependent cyclic electron flow 5 (NDF5) , and expansin-like A3 (EXLA3) . We also observed enrichment in phenylalanine metabolism and other secondary metabolic pathways in clone S7C8. Phenotypic results, upregulated genes and enriched biological pathways identified in this study may explain the enhanced productivity, increased nitrate content, and expanded canopy in polyclonal plantings. Overall, this study provides a foundation for future research to enhance forest productivity by linking molecular mechanisms to practical applications in field plantings.

differential gene expression↗
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↗
DOE OSTI2024

Discovering Innovations in Stress Tolerance through Comparative Gene Regulatory Network Analysis and Cell-Type Specific Expression Maps (Final Technical Report with Cover Page)

Through this grant, we developed a comparative framework to elucidate the mechanisms behind variations in environmental stress responses among a diverse group of species within the Brassicaceae family. Our focus was on the differences in physiological and transcriptomic responses to abscisic acid (ABA), a hormone associated with water stress. We examined the differential growth responses of four Brassicaceae species, finding that most exhibited reduced root growth correlated with smaller meristem size. In contrast, Schrenkiella parvula showed accelerated growth due to increased root cell elongation. We employed RNA sequencing to analyze the transcriptional responses to ABA across these species, and innovative bioinformatics techniques were used to pinpoint biological pathways with significant divergence. Additionally, we utilized DAP-seq to map the gene regulatory networks associated with ABAresponsive transcription factors, revealing that variations in the regulation of growth hormone biosynthesis play a critical role in the distinct ABA effects on root growth among the species. This research sets a new standard for comparative physiology by integrating comparative genomics and transcriptomics to uncover pathway divergences.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Single-molecule epitranscriptomic analysis of full-length HIV-1 RNAs reveals functional roles of site-specific m6As

Abstract Although the significance of chemical modifications on RNA is acknowledged, the evolutionary benefits and specific roles in human immunodeficiency virus (HIV-1) replication remain elusive. Most studies have provided only population-averaged values of modifications for fragmented RNAs at low resolution and have relied on indirect analyses of phenotypic effects by perturbing host effectors. Here we analysed chemical modifications on HIV-1 RNAs at the full-length, single RNA level and nucleotide resolution using direct RNA sequencing methods. Our data reveal an unexpectedly simple HIV-1 modification landscape, highlighting three predominant N 6 -methyladenosine (m 6 A) modifications near the 3′ end. More densely installed in spliced viral messenger RNAs than in genomic RNAs, these m 6 As play a crucial role in maintaining normal levels of HIV-1 RNA splicing and translation. HIV-1 generates diverse RNA subspecies with distinct m 6 A ensembles, and maintaining multiple of these m 6 As on its RNAs provides additional stability and resilience to HIV-1 replication, suggesting an unexplored viral RNA-level evolutionary strategy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Differentiation strategies of soil rare and abundant microbial taxa in response to changing climatic regimes

Despite the important roles of soil microbes, especially the most diverse rare taxa in maintaining community diversity and multifunctionality, how different climate regimes alter the stability and functions of the rare microbial biosphere remains unknown. We reciprocally transplanted field soils across a latitudinal gradient to simulate climate change and sampled the soils annually after harvesting the maize over the following 6 years (from 2005 to 2011). By sequencing microbial 16S ribosomal RNA gene amplicons, we found that changing climate regimes significantly altered the composition and dynamics of soil microbial communities. A continuous succession of the rare and abundant communities was observed. Rare microbial communities were more stable under changing climatic regimes, with lower variations in temporal dynamics, and higher stability and constancy of diversity. More nitrogen cycling genes were detected in the rare members than in the abundant members, including amoA, napA, nifH, nirK, nirS, norB and nrfA. Random forest analysis and receiver operating characteristics analysis showed that rare taxa may act as potential contributors to maize yield under changing climatics. The study indicates that the taxonomically and functionally diverse rare biosphere has the potential to increase functional redundancy and enhance the ability of soil communities to counteract environmental disturbances. With ongoing global climate change, exploring the succession process and functional changes of rare taxa may be important in elucidating the ecosystem stability and multifunctionality that are mediated by microbial communities.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Evaluation of co-circulating pathogens and microbiome from COVID-19 infections

Co-infections or secondary infections with SARS-CoV-2 have the potential to affect disease severity and morbidity. Additionally, the potential influence of the nasal microbiome on COVID-19 illness is not well understood. In this study, we analyzed 203 residual samples, originally submitted for SARS-CoV-2 testing, for the presence of viral, bacterial, and fungal pathogens and non-pathogens using a comprehensive microarray technology, the Lawrence Livermore Microbial Detection Array (LLMDA). Eighty-seven percent of the samples were nasopharyngeal samples, and 23% of the samples were oral, nasal and oral pharyngeal swabs. We conducted bioinformatics analyses to examine differences in microbial populations of these samples, as a proxy for the nasal and oral microbiome, from SARS-CoV-2 positive and negative specimens. We found 91% concordance with the LLMDA relative to a diagnostic RT-qPCR assay for detection of SARS-CoV-2. Sixteen percent of all the samples (32/203) revealed the presence of an opportunistic bacterial or frank viral pathogen with the potential to cause co-infections. The two most detected bacteria, Streptococcus pyogenes and Streptococcus pneumoniae , were present in both SARS-CoV-2 positive and negative samples. Human metapneumovirus was the most prevalent viral pathogen in the SARS-CoV-2 negative samples. Sequence analysis of 16S rRNA was also conducted to evaluate bacterial diversity and confirm LLMDA results.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Spatiotemporal analysis of lung immune dynamics in lethal Coccidioides posadasii infection

Coccidioidomycosis, or Valley fever, is a lung disease caused by inhalation of Coccidioides fungi, prevalent in the Southwestern United States, Mexico, and parts of Central and South America. Annually, the United States reports 10,000–20,000 cases, although those numbers are expected to increase as climate change expands the fungal geographic range. While 60% of infections are asymptomatic, 40% symptomatic infections are often misdiagnosed due to similarities with bronchitis or pneumonia. A small subset of infection progress to severe illness, necessitating a better understanding of immune responses during lethal infection. Using single-cell RNA sequencing and spatial transcriptomics, we characterized lung responses during Coccidioides infection. We identified monocyte-derived Spp1-expressing macrophages as potential mediators of tissue remodeling and fibrosis, marked by high expression of profibrotic and proinflammatory transcripts. These macrophages showed elevated TGF-β and IL-6 signaling, pathways involved in fibrosis pathogenesis. Additionally, we observed significant neutrophil infiltration and defective lymphocyte responses, indicating severe adaptive immunity dysregulation in lethal, acute infection. These findings enhance our understanding of Coccidioides infection and suggest new therapeutic targets.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Depth-resolved sagebrush root metabolomics, rhizosphere microbial communities, and geochemistry at the East River Watershed

This data set consists of results from soil nutrient profile, untargeted metabolomics, mass spec imaging, and amplicon sequencing. Data for soil nutrient profile includes common cations (Ca, Mg, Na, and K etc.) extracted from 3 digesting steps – ammonia acetate (for exchangeable cations), nitric acid (for acid dissolved fraction), and hydrofluoric acid/perchloric acid (HF/HClO4) for whole soil digestion. It also includes concentration of organic carbon, inorganic nitrogen (ammonia and nitrate) and phosphorus (Bray-1 P and nitric acid extract), and total nitrogen and phosphorus. Data for untargeted metabolomics includes metabolomic profile for root exudate/tissues and soil extracts from depths at surface soil to saprolite, that were measured using gas chromatography – mass spectrometry (GC-MS), and liquid chromatography – tandem mass spectrometry (LC-MS/MS). Data for mass spec imaging includes spatial distribution of metabolites that were detected and annotated with Fourier transformation ion cyclotron resonance mass spectrometer (FTICR-MS). Data for amplicon sequencing includes the base paired 16S and ITS ribosomal RNA sequences from Miseq Illumina sequencing. All samples were collected from 2 sampling campaign October 2022 and June 2023. Collectively, these datasets enable a mechanistic evaluation of how nutrient acquisition, especially nitrogen and phosphorus, differs between shallow roots operating in soil and deep roots functioning within the fractured bedrock zone. All files are provided as comma-separated values (CSV) fies (.csv) and (GZIP) file (.gz). The compressed .gz FASTQ files can be read directly in R using the dada2 package as part of the amplicon sequence analysis workflow. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231. This research was performed on a project award 60563 (https://dx.doi.org/10.46936/expl.proj.2022.60563/60008727) from the Environmental Molecular Sciences Laboratory, a DOE Office of Science User Facility sponsored by the Biological and Environmental Research program under Contract No. DE-AC05-76RL01830.

EARTH SCIENCE > AGRICULTURE > SOILS > CARBON↗
DOE OSTI2020

Single-Cell Transcriptomic Analysis of Tumor-Derived Fibroblasts and Normal Tissue-Resident Fibroblasts Reveals Fibroblast Heterogeneity in Breast Cancer

Cancer-associated fibroblasts (CAFs) are a prominent stromal cell type in solid tumors and molecules secreted by CAFs play an important role in tumor progression and metastasis. CAFs coexist as heterogeneous populations with potentially different biological functions. Although CAFs are a major component of the breast cancer stroma, molecular and phenotypic heterogeneity of CAFs in breast cancer is poorly understood. In this study, we investigated CAF heterogeneity in triple-negative breast cancer (TNBC) using a syngeneic mouse model, BALB/c-derived 4T1 mammary tumors. Using single-cell RNA sequencing (scRNA-seq), we identified six CAF subpopulations in 4T1 tumors including: 1) myofibroblastic CAFs, enriched for α-smooth muscle actin and several other contractile proteins; 2) ‘inflammatory’ CAFs with elevated expression of inflammatory cytokines; and 3) a CAF subpopulation expressing major histocompatibility complex (MHC) class II proteins that are generally expressed in antigen-presenting cells. Comparison of 4T1-derived CAFs to CAFs from pancreatic cancer revealed that these three CAF subpopulations exist in both tumor types. Interestingly, cells with inflammatory and MHC class II-expressing CAF profiles were also detected in normal breast/pancreas tissue, suggesting that these phenotypes are not tumor microenvironment-induced. This work enhances our understanding of CAF heterogeneity, and specifically targeting these CAF subpopulations could be an effective therapeutic approach for treating highly aggressive TNBCs.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Saccharomycopsis sp. strain UFMG-CM-Y6632 internal transcribed spacer 1, partial sequence; 5.8S ribosomal RNA gene and internal transcribed spacer 2, complete sequence; and large subunit ribosomal RNA gene, partial sequence

Three yeast isolates were obtained from soil and rotting wood samples collected in an Amazonian rainforest biome in Brazil. Comparison of the intergenic spacer 5.8S region and the D1/D2 domains of the large subunit rRNA gene showed that the isolates represent a novel species of the genus Saccharomycopsis. A tree inferred from the D1/D2 sequences placed the novel species near a subclade containing Saccharomycopsis lassenensis, Saccharomycopsis fermentans, Saccharomycopsis javanensis, Saccharomycopsis babjevae, Saccharomycopsis schoenii and Saccharomycopsis oosterbeekiorum, but with low bootstrap support. In terms of sequence divergence, the novel species had the highest identity in the D1/D2 domains with Saccharomycopsis capsularis, from which it differed by 36 substitutions. In contrast, a phylogenomic analysis based on 1061 single-copy orthologs for a smaller set of Saccharomycopsis species whose whole genome sequences are available indicated that the novel species represented by strain UFMG-CM-Y6991 is phylogenetically closer to Saccharomycopsis fodiens and Saccharomycopsis sp. TF2021a (=Saccharomycopsis phalluae). The novel yeast is homothallic and produces asci with one spheroidal ascospore with an equatorial or subequatorial ledge. The name Saccharomycopsis praedatoria sp. nov. is proposed to accommodate the novel species. The holotype of Saccharomycopsis praedatoria is CBS 16589(T). The MycoBank number is MB849369. S. praedatoria was able to kill cells of Saccharomyces cerevisiae by means of penetration with infection pegs, a trait common to most species of Saccharomycopsis.

Amazonian Forest↗
DOE OSTI2022

Saccharomycopsis sp. isolate UFMG-CM-Y6632 large subunit ribosomal RNA gene, partial sequence

Three yeast isolates were obtained from soil and rotting wood samples collected in an Amazonian rainforest biome in Brazil. Comparison of the intergenic spacer 5.8S region and the D1/D2 domains of the large subunit rRNA gene showed that the isolates represent a novel species of the genus Saccharomycopsis. A tree inferred from the D1/D2 sequences placed the novel species near a subclade containing Saccharomycopsis lassenensis, Saccharomycopsis fermentans, Saccharomycopsis javanensis, Saccharomycopsis babjevae, Saccharomycopsis schoenii and Saccharomycopsis oosterbeekiorum, but with low bootstrap support. In terms of sequence divergence, the novel species had the highest identity in the D1/D2 domains with Saccharomycopsis capsularis, from which it differed by 36 substitutions. In contrast, a phylogenomic analysis based on 1061 single-copy orthologs for a smaller set of Saccharomycopsis species whose whole genome sequences are available indicated that the novel species represented by strain UFMG-CM-Y6991 is phylogenetically closer to Saccharomycopsis fodiens and Saccharomycopsis sp. TF2021a (=Saccharomycopsis phalluae). The novel yeast is homothallic and produces asci with one spheroidal ascospore with an equatorial or subequatorial ledge. The name Saccharomycopsis praedatoria sp. nov. is proposed to accommodate the novel species. The holotype of Saccharomycopsis praedatoria is CBS 16589(T). The MycoBank number is MB849369. S. praedatoria was able to kill cells of Saccharomyces cerevisiae by means of penetration with infection pegs, a trait common to most species of Saccharomycopsis.

Amazonian Forest↗
DOE OSTI2024

Analysis of RNA viral function from metatranscriptomic reads in drought simulated soils

Soil viruses have received increased interest in recent years because of their abundance and unelucidated potential as drivers for organic carbon derived from fungal and bacterial necromass. With decreasing levels of rainfall and climbing instances of unseasonal conditions, it becomes crucial to our goals for negative carbon emissions to understand how soil carbon flux is affected. Here we leveraged metatranscriptomic data from a section of research plot of Hopland Research and Extension in northern California with a suite of viral detection tools: VirSorter, DRAMv, and VContact2. We aimed to cross analyze moisture regiments with viral populations. We find that viral diversity is driven by site, like bacteria. We found that viral diversity at face value showed some differentiation but did not significantly differ by legacy moisture content. Additionally, there was detection of auxiliary metabolic genes in a viral sequence which may further support virus’s direct input into the soil carbon pool.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

The Transcriptional Response of Soil Bacteria to Long-Term Warming and Short-Term Seasonal Fluctuations in a Terrestrial Forest

Terrestrial ecosystems are an important carbon store, and this carbon is vulnerable to microbial degradation with climate warming. After 30 years of experimental warming, carbon stocks in a temperate mixed deciduous forest were observed to be reduced by 30% in the heated plots relative to the controls. In addition, soil respiration was seasonal, as was the warming treatment effect. We therefore hypothesized that long-term warming will have higher expressions of genes related to carbohydrate and lipid metabolism due to increased utilization of recalcitrant carbon pools compared to controls. Because of the seasonal effect of soil respiration and the warming treatment, we further hypothesized that these patterns will be seasonal. We used RNA sequencing to show how the microbial community responds to long-term warming (~30 years) in Harvard Forest, MA. Total RNA was extracted from mineral and organic soil types from two treatment plots (+5°C heated and ambient control), at two time points (June and October) and sequenced using Illumina NextSeq technology. Treatment had a larger effect size on KEGG annotated transcripts than on CAZymes, while soil types more strongly affected CAZymes than KEGG annotated transcripts, though effect sizes overall were small. Although, warming showed a small effect on overall CAZymes expression, several carbohydrate-associated enzymes showed increased expression in heated soils (~68% of all differentially expressed transcripts). Further, exploratory analysis using an unconstrained method showed increased abundances of enzymes related to polysaccharide and lipid metabolism and decomposition in heated soils. Compared to long-term warming, we detected a relatively small effect of seasonal variation on community gene expression. Together, these results indicate that the higher carbohydrate degrading potential of bacteria in heated plots can possibly accelerate a self-reinforcing carbon cycle-temperature feedback in a warming climate.

54 ENVIRONMENTAL SCIENCES↗