Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “transcriptome”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Identification and functional analysis of strigolactone pathway genes regulating tillering traits in sugarcane

Abstract Saccharum officinarum and Saccharum spontaneum are two fundamental species of modern sugarcane cultivars, exhibiting divergent tillering patterns crucial for sugarcane architecture and yield. Strigolactones (SLs), a class of plant hormones, are considered to play a central role in shaping plant form and regulating tillering. Our study highlights the distinct tillering patterns observed between S. officinarum and S. spontaneum and implicates significant differences in SL levels in root exudates between the two species. Treatment with rac-GR24 (an artificial SL analog) suppressed tillering in S. spontaneum. Based on transcriptome analysis, we focused on two genes, TRANSCRIPTION ELONGATION FACTOR 1 (TEF1) and CIRCADIAN CLOCK ASSOCIATED1 (CCA1), which show higher expression in S. spontaneum or S. officinarum, respectively. While the overexpression of SoCCA1 did not lead to significant phenotypic differences, overexpression of SsTEF1 in rice stimulated tillering and inhibited plant height, demonstrating its role in tillering regulation. However, the overexpression of suggests that SoCCA1 may not be the key regulator of sugarcane tillering. Yeast one-hybrid assays identified four transcription factors (TFs) regulating SsTEF1 and four and five TFs regulating SsCCA1 and SoCCA1. This study provides a theoretical foundation for deciphering the molecular mechanisms underlying the different tillering behaviors between S. officinarum and S. spontaneum, providing valuable insights for the molecular-based design of sugarcane breeding strategies.

Qi, Yiying↗
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 OSTI2024

The N-terminal domains of NLR immune receptors exhibit structural and functional similarities across divergent plant lineages

Nucleotide-binding domain and leucine-rich repeat (NLR) proteins are a prominent class of intracellular immune receptors in plants. However, our understanding of plant NLR structure and function is limited to the evolutionarily young flowering plant clade. Here, we describe an extended spectrum of NLR diversity across divergent plant lineages and demonstrate the structural and functional similarities of N-terminal domains that trigger immune responses. We show that the broadly distributed coiled-coil (CC) and toll/interleukin-1 receptor (TIR) domain families of nonflowering plants retain immune-related functions through translineage activation of cell death in the angiosperm Nicotiana benthamiana. We further examined a CC subfamily specific to nonflowering lineages and uncovered an essential N-terminal MAEPL motif that is functionally comparable with motifs in resistosome-forming CC-NLRs. Consistent with a conserved role in immunity, the ectopic activation of CC MAEPL in the nonflowering liverwort Marchantia polymorpha led to profound growth inhibition, defense gene activation, and signatures of cell death. Moreover, comparative transcriptomic analyses of CC MAEPL activity delineated a common CC-mediated immune program shared across evolutionarily divergent nonflowering and flowering plants. Collectively, our findings highlight the ancestral nature of NLR-mediated immunity during plant evolution that dates its origin to at least ~500 million years ago.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

The DYRKP1 kinase regulates cell wall degradation in Chlamydomonas by inducing matrix metalloproteinase expression

Abstract The cell wall of plants and algae is an important cell structure that protects cells from changes in the external physical and chemical environment. This extracellular matrix, composed of polysaccharides and glycoproteins, must be constantly remodeled throughout the life cycle. However, compared to matrix polysaccharides, little is known about the mechanisms regulating the formation and degradation of matrix glycoproteins. We report here that a plant kinase belonging to the dual-specificity tyrosine phosphorylation-regulated kinase (DYRKP1) family present in all eukaryotes regulates cell wall degradation after mitosis of Chlamydomonas reinhardtii by inducing the expression of matrix metalloproteinases. Without DYRKP1, daughter cells cannot disassemble parental cell walls and remain trapped inside for more than 10 days. On the other hand, the dual-specificity tyrosine phosphorylation-regulated kinase complementation lines show normal degradation of the parental cell wall. Transcriptomic and proteomic analyses indicate a marked downregulation of MMP gene expression and accumulation, respectively, in the dyrkp1 mutants. The mutants deficient in matrix metalloproteinases retain palmelloid structures for a longer time than the background strain, like dyrkp1 mutants. Our findings show that dual-specificity tyrosine phosphorylation-regulated kinase, by ensuring timely MMP expression, enables the successful execution of the cell cycle. Altogether, this study provides insight into the life cycle regulation in plants and algae.

Kim, Minjae (ORCID:0000000223561295)↗
DOE OSTI2025

Too dim, too bright, and just right: Systems analysis of the Chlamydomonas diurnal program under limiting and excess light

Photosynthetic organisms coordinate their metabolism and growth with diurnal light, which can range in intensity from limiting to excessive. Little is known about how light intensity impacts the diurnal program in Chlamydomonas reinhardtii, or how diurnal rhythms in gene expression and metabolism shape photoprotective responses at different times of day. To address these questions, we performed a systems analysis of synchronized Chlamydomonas populations acclimated to low, moderate, and high diurnal light. Transcriptomic and proteomic data revealed that the Chlamydomonas rhythmic gene expression program is resilient to limiting and excess light: genome-wide, waves of transcripts, and proteins peak at the same times in populations acclimated to stressful light intensities as in populations acclimated to moderate light. Yet, diurnal photoacclimation gives rise to hundreds of gene expression changes, even at night. Time course measurements of photosynthetic efficiency and pigments responsive to excess light showed that high light-acclimated cells partially overcome photodamage in the latter half of the day prior to cell division. Although gene expression and photodamage are dynamic over the diurnal cycle, Chlamydomonas populations acclimated to low and high diurnal light maintain altered photosystem abundance, thylakoid architecture, and non-photochemical quenching capacity through the night phase. This suggests that cells remember or anticipate the light intensities that they have typically encountered during the day. The integrated data constitute an excellent resource for understanding photoacclimation in eukaryotes under environmentally relevant conditions.

Dupuis, Sunnyjoy↗
DOE OSTI2026

Gene and genome duplications have contrasting impacts on biosynthetic and flower developmental pathways in California poppy

Benzylisoquinoline alkaloids (BIAs) represent a vast group of specialized plant metabolites with diverse pharmaceutical applications, synthesized by a variety of gene families. Among the multiple plant lineages that produce BIAs, the most notable is the poppy family (Papaveraceae), with California poppy (Eschscholzia californica) emerging as a model organism. Here, we report a haplotype-resolved genome assembly, in combination with a high-density expression atlas, for California poppy. Genome analyses reveal recent diversification of BIA biosynthesis genes in poppy through localized duplications. Furthermore, we demonstrate that the degree of phylogenetic relatedness among paralogs within BIA biosynthesis-associated gene families correlates with similarities in gene expression. In contrast, gene families involved in carotenoid biosynthesis, which contributes to the intense orange petal pigmentation, are not phylogenetically clustered, and floral developmental regulators exhibit a high degree of retention of gene duplicates associated with ancient polyploidy events. These findings illustrate alternative roles for gene and genome duplications as drivers of trait evolution. Given the position of California poppy in the angiosperm phylogeny, the high-quality genomic resources generated for this work constitute a valuable resource for comparative genomic and transcriptomic analyses for poppies and flowering plants more generally.

Rössner, Le-Han [Justus-Liebig University, Giessen↗
DOE OSTI2021

Machine learning-enabled phenotyping for GWAS and TWAS of WUE traits in 869 field-grown sorghum accessions

Abstract Sorghum (Sorghum bicolor) is a model C4 crop made experimentally tractable by extensive genomic and genetic resources. Biomass sorghum is studied as a feedstock for biofuel and forage. Mechanistic modeling suggests that reducing stomatal conductance (gs) could improve sorghum intrinsic water use efficiency (iWUE) and biomass production. Phenotyping to discover genotype-to-phenotype associations remains a bottleneck in understanding the mechanistic basis for natural variation in gs and iWUE. This study addressed multiple methodological limitations. Optical tomography and a machine learning tool were combined to measure stomatal density (SD). This was combined with rapid measurements of leaf photosynthetic gas exchange and specific leaf area (SLA). These traits were the subject of genome-wide association study and transcriptome-wide association study across 869 field-grown biomass sorghum accessions. The ratio of intracellular to ambient CO2 was genetically correlated with SD, SLA, gs, and biomass production. Plasticity in SD and SLA was interrelated with each other and with productivity across wet and dry growing seasons. Moderate-to-high heritability of traits studied across the large mapping population validated associations between DNA sequence variation or RNA transcript abundance and trait variation. A total of 394 unique genes underpinning variation in WUE-related traits are described with higher confidence because they were identified in multiple independent tests. This list was enriched in genes whose Arabidopsis (Arabidopsis thaliana) putative orthologs have functions related to stomatal or leaf development and leaf gas exchange, as well as genes with nonsynonymous/missense variants. These advances in methodology and knowledge will facilitate improving C4 crop WUE.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Phenotyping stomatal closure by thermal imaging for GWAS and TWAS of water use efficiency-related genes

Abstract Stomata allow CO2 uptake by leaves for photosynthetic assimilation at the cost of water vapor loss to the atmosphere. The opening and closing of stomata in response to fluctuations in light intensity regulate CO2 and water fluxes and are essential for maintaining water-use efficiency (WUE). However, a little is known about the genetic basis for natural variation in stomatal movement, especially in C4 crops. This is partly because the stomatal response to a change in light intensity is difficult to measure at the scale required for association studies. Here, we used high-throughput thermal imaging to bypass the phenotyping bottleneck and assess 10 traits describing stomatal conductance (gs) before, during and after a stepwise decrease in light intensity for a diversity panel of 659 sorghum (Sorghum bicolor) accessions. Results from thermal imaging significantly correlated with photosynthetic gas exchange measurements. gs traits varied substantially across the population and were moderately heritable (h2 up to 0.72). An integrated genome-wide and transcriptome-wide association study identified candidate genes putatively driving variation in stomatal conductance traits. Of the 239 unique candidate genes identified with the greatest confidence, 77 were putative orthologs of Arabidopsis (Arabidopsis thaliana) genes related to functions implicated in WUE, including stomatal opening/closing (24 genes), stomatal/epidermal cell development (35 genes), leaf/vasculature development (12 genes), or chlorophyll metabolism/photosynthesis (8 genes). These findings demonstrate an approach to finding genotype-to-phenotype relationships for a challenging trait as well as candidate genes for further investigation of the genetic basis of WUE in a model C4 grass for bioenergy, food, and forage production.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Transcript and metabolite network perturbations in lignin biosynthetic mutants of Arabidopsis

Abstract Lignin, one of the most abundant polymers in plants, is derived from the phenylpropanoid pathway, which also gives rise to an array of metabolites that are essential for plant fitness. Genetic engineering of lignification can cause drastic changes in transcription and metabolite accumulation with or without an accompanying development phenotype. To understand the impact of lignin perturbation, we analyzed transcriptome and metabolite data from the rapidly lignifying stem tissue in 13 selected phenylpropanoid mutants and wild-type Arabidopsis (Arabidopsis thaliana). Our dataset contains 20,974 expressed genes, of which over 26% had altered transcript levels in at least one mutant, and 18 targeted metabolites, all of which displayed altered accumulation in at least one mutant. We found that lignin biosynthesis and phenylalanine supply via the shikimate pathway are tightly co-regulated at the transcriptional level. The hierarchical clustering analysis of differentially expressed genes (DEGs) grouped the 13 mutants into 5 subgroups with similar profiles of mis-regulated genes. Functional analysis of the DEGs in these mutants and correlation between gene expression and metabolite accumulation revealed system-wide effects on transcripts involved in multiple biological processes.

Plant Sciences↗
DOE OSTI2024

A maize semi-dwarf mutant reveals a GRAS transcription factor involved in brassinosteroid signaling

Brassinosteroids (BR) and gibberellins (GA) regulate plant height and leaf angle in maize (Zea mays). Mutants with defects in BR or GA biosynthesis or signaling identify components of these pathways and enhance our knowledge about plant growth and development. In this study, we characterized three recessive mutant alleles of GRAS transcription factor 42 (gras42) in maize, a GRAS transcription factor gene orthologous to the DWARF AND LOW TILLERING (DLT) gene of rice (Oryza sativa). These maize mutants exhibited semi-dwarf stature, shorter and wider leaves, and more upright leaf angle. Transcriptome analysis revealed a role for GRAS42 as a determinant of BR signaling. Analysis of the expression consequences from loss of GRAS42 in the gras42-mu1021149 mutant indicated a weak loss of BR signaling in the mutant, consistent with its previously demonstrated role in BR signaling in rice. Loss of BR signaling was also evident by the enhancement of weak BR biosynthetic mutant alleles in double mutants of nana plant1-1 and gras42-mu1021149. The gras42-mu1021149 mutant had little effect on GA-regulated gene expression, suggesting that GRAS42 is not a regulator of core GA signaling genes in maize. Single-cell expression data identified gras42 expressed among cells in the G2/M phase of the cell cycle consistent with its previously demonstrated role in cell cycle gene expression in Arabidopsis (Arabidopsis thaliana). Cis-acting natural variation controlling GRAS42 transcript accumulation was identified by expression genome-wide association study (eGWAS) in maize. Our results demonstrate a conserved role for GRAS42/SCARECROW-LIKE 28 (SCL28)/DLT in BR signaling, clarify the role of this gene in GA signaling, and suggest mechanisms of tillering and leaf angle control by BR.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Packaging “vegetable oils”: Insights into plant lipid droplet proteins

Abstract Plant neutral lipids, also known as “vegetable oils”, are synthesized within the endoplasmic reticulum (ER) membrane and packaged into subcellular compartments called lipid droplets (LDs) for stable storage in the cytoplasm. The biogenesis, modulation, and degradation of cytoplasmic LDs in plant cells are orchestrated by a variety of proteins localized to the ER, LDs, and peroxisomes. Recent studies of these LD-related proteins have greatly advanced our understanding of LDs not only as steady oil depots in seeds but also as dynamic cell organelles involved in numerous physiological processes in different tissues and developmental stages of plants. In the past 2 decades, technology advances in proteomics, transcriptomics, genome sequencing, cellular imaging and protein structural modeling have markedly expanded the inventory of LD-related proteins, provided unprecedented structural and functional insights into the protein machinery modulating LDs in plant cells, and shed new light on the functions of LDs in nonseed plant tissues as well as in unicellular algae. Here, we review critical advances in revealing new LD proteins in various plant tissues, point out structural and mechanistic insights into key proteins in LD biogenesis and dynamic modulation, and discuss future perspectives on bridging our knowledge gaps in plant LD biology.

Cai, Yingqi (ORCID:0000000203575809)↗
DOE OSTI2026

Metabolic flux, metabolite, and transcript analysis uncover reprogramming of metabolism toward higher seed oil

Overexpression of WRINKLED1 (WRI1), a master regulator of glycolysis and fatty acid biosynthesis, together with DIACYLGLYCEROL ACYLTRANSFERASE1 (DGAT1), which catalyzes the final step of triacylglycerol assembly, is a promising strategy for enhancing seed oil content. However, how these regulators coordinate system-wide metabolic reprogramming at the levels of gene expression, metabolite pools, and fluxes remains poorly understood. To address this, we performed 13 C-metabolic flux analysis, metabolomics, and transcriptomics on in vitro cultured pennycress (Thlaspi arvense L.) embryos overexpressing the native WRI1 and DGAT1 homologs. Here, in cultured embryos, WRI1/DGAT1 overexpression increased triacylglycerol accumulation by 28% while reducing protein content by 34%, relative to the wild type. Embryos showed ∼20-fold and 50-fold upregulation of WRI1 and DGAT1 along with induction of WRI1 target genes in glycolysis and fatty acid biosynthesis. Genes associated with photosynthesis and Calvin cycle functions were also upregulated, whereas genes encoding ribosomal proteins and seed storage proteins were strongly repressed, consistent with the observed lipid–protein tradeoff. Flux analysis revealed that enhanced triacylglycerol biosynthesis is supported by increased flux through the Rubisco shunt and cytosolic pyruvate kinase, while the oxidative pentose phosphate pathway and malic enzyme contributed little to NADPH or pyruvate supply. Metabolomic profiling revealed extensive perturbations in glycolytic intermediates, tricarboxylic acid cycle metabolites, and amino acids. In plant grown seeds, WRI1/DGAT1 lines also showed a modest but significant increase in total lipid content. Collectively, these findings reveal how WRI1 and DGAT1 reprogram central metabolism to enhance oil accumulation, with relevance to mature seeds.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Predicting oncology drug-induced cardiotoxicity with donor-specific iPSC-CMs—a proof-of-concept study with doxorubicin

Abstract Many oncology drugs have been found to induce cardiotoxicity in a subset of patients, which significantly limits their clinical use and impedes the benefit of lifesaving anticancer treatments. Human induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) carry donor-specific genetic information and have been proposed for exploring the interindividual difference in oncology drug-induced cardiotoxicity. Herein, we evaluated the inter- and intraindividual variability of iPSC-CM-related assays and presented a proof of concept to prospectively predict doxorubicin (DOX)-induced cardiotoxicity (DIC) using donor-specific iPSC-CMs. Our findings demonstrated that donor-specific iPSC-CMs exhibited greater line-to-line variability than the intraindividual variability in impedance cytotoxicity and transcriptome assays. The variable and dose-dependent cytotoxic responses of iPSC-CMs resembled those observed in clinical practice and largely replicated the reported mechanisms. By categorizing iPSC-CMs into resistant and sensitive cell lines based on their time- and concentration-related phenotypic responses to DOX, we found that the sensitivity of donor-specific iPSC-CMs to DOX may predict in vivo DIC risk. Furthermore, we identified a differentially expressed gene, DND microRNA-mediated repression inhibitor 1 (DND1), between the DOX-resistant and DOX-sensitive iPSC-CMs. Our results support the utilization of donor-specific iPSC-CMs in assessing interindividual differences in DIC. Further studies will encompass a large panel of donor-specific iPSC-CMs to identify potential novel molecular and genetic biomarkers for predicting DOX and other oncology drug-induced cardiotoxicity.

Toxicology↗
DOE OSTI2022

Quantitative Resistance Loci to Southern Rust Mapped in a Temperate Maize Diversity Panel

Southern rust is a severe foliar disease of maize (Zea mays) resulting from infection with the obligate biotrophic fungus Puccinia polysora. This disease reduces photosynthetic productivity, which in turn reduces yields, with the greatest yield losses (up to 50%) associated with earlier onset infections. P. polysora urediniospores overwinter only in tropical and subtropical regions but cause outbreaks when environmental conditions favor initial infection. Increased temperatures and humidity during the growing season combined with an increased frequency of moderate winters are likely to increase the frequency of severe southern rust outbreaks in the U.S. Corn Belt. In summer 2020, a severe outbreak of southern rust was observed in eastern Nebraska, United States. We scored a replicated maize association panel planted in Lincoln, NE for disease severity and found that disease incidence and severity showed significant variation among maize genotypes. Genome-wide association studies identified four loci associated with significant quantitative variation in disease severity. These loci were associated with candidate genes with plausible links to quantitative disease resistance. A transcriptome-wide association study identified additional genes associated with disease severity. Together, these results indicate that substantial diversity in resistance to southern rust exists among current temperate-adapted maize germplasm, including several candidate loci that may explain the observed variation in resistance to southern rust. [Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY 4.0 International license .

Plant Sciences↗
DOE OSTI2024

A Sorghum F-Box Protein Induces an Oxidative Burst in the Defense Against Colletotrichum sublineola

The hemibiotrophic fungal pathogen Colletotrichum sublineola is the causal agent of anthracnose in sorghum ( Sorghum bicolor), resulting in leaf blight, stalk rot, and head blight in susceptible genotypes, with yield losses of up to 50%. The development of anthracnose-resistant cultivars can reduce reliance on fungicides and provide a more sustainable and economical means for disease management. A previous genome-wide association study of the sorghum association panel identified the candidate resistance gene Sobic.005G172300 encoding an F-box protein. To better understand the role of this gene in the defense against C. sublineola, gene expression following infection with C. sublineola was monitored by RNA sequencing in seedlings of sorghum accession SC110, which harbored the resistance allele, and three accessions that harbored a susceptible allele. Only in SC110 did the expression of Sobic.005G172300 increase during the biotrophic phase of infection. Subsequent transcriptome analysis, gene co-expression networks, and gene regulatory networks of inoculated and mock-inoculated seedlings of resistant and susceptible accessions suggest that the increase in expression of Sobic.005G172300 induces an oxidative burst by lowering the concentration of ascorbic acid during the biotrophic phase of infection. Based on gene regulatory network analysis, the protein encoded by Sobic.005G172300 is proposed to target proteins involved in the biosynthesis of ascorbic acid for polyubiquitination through the SCF E3 ubiquitin ligase, causing their degradation via the proteasome.

Plant Sciences↗
DOE OSTI2024

BRAKER3: Fully automated genome annotation using RNA-seq and protein evidence with GeneMark-ETP, AUGUSTUS, and TSEBRA

Gene prediction has remained an active area of bioinformatics research for a long time. Still, gene prediction in large eukaryotic genomes presents a challenge that must be addressed by new algorithms. The amount and significance of the evidence available from transcriptomes and proteomes vary across genomes, between genes, and even along a single gene. User-friendly and accurate annotation pipelines that can cope with such data heterogeneity are needed. The previously developed annotation pipelines BRAKER1 and BRAKER2 use RNA-seq or protein data, respectively, but not both. A further significant performance improvement integrating all three data types was made by the recently released GeneMark-ETP. We here present the BRAKER3 pipeline that builds on GeneMark-ETP and AUGUSTUS, and further improves accuracy using the TSEBRA combiner. BRAKER3 annotates protein-coding genes in eukaryotic genomes using both short-read RNA-seq and a large protein database, along with statistical models learned iteratively and specifically for the target genome. We benchmarked the new pipeline on genomes of 11 species under an assumed level of relatedness of the target species proteome to available proteomes. BRAKER3 outperforms BRAKER1 and BRAKER2. The average transcript-level F1-score is increased by about 20 percentage points on average, whereas the difference is most pronounced for species with large and complex genomes. BRAKER3 also outperforms other existing tools, MAKER2, Funannotate, and FINDER. The code of BRAKER3 is available on GitHub and as a ready-to-run Docker container for execution with Docker or Singularity. Overall, BRAKER3 is an accurate, easy-to-use tool for eukaryotic genome annotation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

A Graphical Model for Fusing Diverse Microbiome Data

This paper develops a Bayesian graphical model for fusing disparate types of count data. The motivating application is the study of bacterial communities from diverse high-dimensional features, in this case, transcripts, collected from different treatments. In such datasets, there are no explicit correspondences between the communities and each corresponds to different factors, making data fusion challenging. We introduce a flexible multinomial-Gaussian generative model for jointly modeling such count data. This latent variable model jointly characterizes the observed data through a common multivariate Gaussian latent space that parameterizes the set of multinomial probabilities of the transcriptome counts. The covariance matrix of the latent variables induces a covariance matrix of co-dependencies between all the transcripts, effectively fusing multiple data sources. We present a computationally scalable variational Expectation-Maximization (EM) algorithm for inferring the latent variables and the parameters of the model. Here, the inferred latent variables provide a common dimensionality reduction for visualizing the data and the inferred parameters provide a predictive posterior distribution. In addition to simulation studies that demonstrate the variational EM procedure, we apply our model to a bacterial microbiome dataset.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Desiccation induces varied responses within a soil bacterial genus

Abstract Desiccation impacts a suite of physiological processes in microbes by elevating levels of damaging reactive oxygen species and inducing DNA strand breaks. In response to desiccation‐induced stress, microbes have evolved specialized mechanisms to help them survive. Here, we performed a 128‐day lab desiccation experiment on nine strains from three clades of an abundant soil bacterium, Curtobacterium . We sequenced RNA from each strain at three time points to investigate their response. Curtobacterium was highly resistant to desiccation, outlasting both Escherichia coli and a famously DNA damage‐resistant bacterium, Deinococcus radiodurans . However, within the genus, there were also 10‐fold differences in survival rates among strains. Transcriptomic profiling revealed responses shared within the genus including up‐regulation of genes involved in DNA damage repair, osmolyte production, and efflux pumps, but also up‐regulation of pathways and genes unique to the three clades. For example, trehalose synthesis gene otsB , the chaperone groEL , and the oxygen scavenger katA were all found in either one or two clades but not the third. Here, we provide evidence of considerable variation in closely related strains, and further elucidation of the phylogenetic conservation of desiccation tolerance remains an important goal for microbial ecologists.

59 BASIC BIOLOGICAL SCIENCES↗