Search NASA⌕ Search

SEARCH · Search NASA

Results for “Genome-wide association studies”

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

Genome-wide profiling of histone (H3) lysine 4 (K4) tri-methylation (me3) under drought, heat, and combined stresses in switchgrass

Background: Switchgrass (Panicum virgatum L.) is a warm-season perennial (C4) grass identified as an important biofuel crop in the United States. It is well adapted to the marginal environment where heat and moisture stresses predominantly affect crop growth. However, the underlying molecular mechanisms associated with heat and drought stress tolerance still need to be fully understood in switchgrass. The methylation of H3K4 is often associated with transcriptional activation of genes, including stress-responsive. Therefore, this study aimed to analyze genome-wide histone H3K4-tri-methylation in switchgrass under heat, drought, and combined stress. Results: In total, ~ 1.3 million H3K4me3 peaks were identified in this study using SICER. Among them, 7,342; 6,510; and 8,536 peaks responded under drought (DT), drought and heat (DTHT), and heat (HT) stresses, respectively. Most DT and DTHT peaks spanned 0 to + 2000 bases from the transcription start site [TSS]. By comparing differentially marked peaks with RNA-Seq data, we identified peaks associated with genes: 155 DT-responsive peaks with 118 DT-responsive genes, 121 DTHT-responsive peaks with 110 DTHT-responsive genes, and 175 HT-responsive peaks with 136 HT-responsive genes. We have identified various transcription factors involved in DT, DTHT, and HT stresses. Gene Ontology analysis using the AgriGO revealed that most genes belonged to biological processes. Most annotated peaks belonged to metabolite interconversion, RNA metabolism, transporter, protein modifying, defense/immunity, membrane traffic protein, transmembrane signal receptor, and transcriptional regulator protein families. Further, we identified significant peaks associated with TFs, hormones, signaling, fatty acid and carbohydrate metabolism, and secondary metabolites. qRT-PCR analysis revealed the relative expressions of six abiotic stress-responsive genes (transketolase, chromatin remodeling factor-CDH3, fatty-acid desaturase A, transmembrane protein 14C, beta-amylase 1, and integrase-type DNA binding protein genes) that were significantly (P < 0.05) marked during drought, heat, and combined stresses by comparing stress-induced against un-stressed and input controls. Conclusion: Our study provides a comprehensive and reproducible epigenomic analysis of drought, heat, and combined stress responses in switchgrass. Significant enrichment of H3K4me3 peaks downstream of the TSS of protein-coding genes was observed. In addition, the cost-effective experimental design, modified ChIP-Seq approach, and analyses presented here can serve as a prototype for other non-model plant species for conducting stress studies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Diversity and scale: Genetic architecture of 2068 traits in the VA Million Veteran Program

One of the justifiable criticisms of human genetic studies is the underrepresentation of participants from diverse populations. Lack of inclusion must be addressed at-scale to identify causal disease factors and understand the genetic causes of health disparities. We present genome-wide associations for 2068 traits from 635,969 participants in the Department of Veterans Affairs Million Veteran Program, a longitudinal study of diverse United States Veterans. Systematic analysis revealed 13,672 genomic risk loci; 1608 were only significant after including non-European populations. Fine-mapping identified causal variants at 6318 signals across 613 traits. One-third (n = 2069) were identified in participants from non-European populations. This reveals a broadly similar genetic architecture across populations, highlights genetic insights gained from underrepresented groups, and presents an extensive atlas of genetic associations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

HighDimMixedModels.jl: Robust high-dimensional mixed-effects models across omics data

High-dimensional mixed-effects models are an increasingly important form of regression in which the number of covariates rivals or exceeds the number of samples, which are collected in groups or clusters. The penalized likelihood approach to fitting these models relies on a coordinate descent algorithm that lacks guarantees of convergence to a global optimum. Here, we empirically study the behavior of this algorithm on simulated and real examples of three types of data that are common in modern biology: transcriptome, genome-wide association, and microbiome data. Our simulations provide new insights into the algorithm’s behavior in these settings, and, comparing the performance of two popular penalties, we demonstrate that the smoothly clipped absolute deviation (SCAD) penalty consistently outperforms the least absolute shrinkage and selection operator (LASSO) penalty in terms of both variable selection and estimation accuracy across omics data. To empower researchers in biology and other fields to fit models with the SCAD penalty, we implement the algorithm in a Julia package, HighDimMixedModels.jl .

Gorstein, Evan↗
NASA NTRSDate not supplied

Somatic Mutation Analysis in Spaceflight: NASA Twins Genome Study

The NASA Twins Genome Study investigates the effects of spaceflight on somatic mutation accumulation by comparing genome-wide sequence data from a spaceflight astronaut and his Earth-bound twin. Utilizing advanced computational software on high performance computers, this study identifies and maps somatic mutations, with implications for understanding spaceflight-associated health risks, including cancer, neurodegeneration, and cardiovascular disease. The findings aim to bridge rodent and human space research, offering insights into tissue-specific pathophysiology, risk models, and potential therapeutic interventions.

somatic mutation↗
DOE OSTI2024

Genetic underpinnings of arthropod community distributions in Populus trichocarpa

Community genetics seeks to understand the mechanisms by which natural genetic variation in heritable host phenotypes can encompass assemblages of organisms such as bacteria, fungi, and many animals including arthropods. Prior studies that focused on plant genotypes have been unable to identify genes controlling community composition, a necessary step to predict ecosystem structure and function as underlying genes shift within plant populations. Here, we surveyed arthropods within an association population of Populus trichocarpa in three common gardens to discover plant genes that contributed to arthropod community composition. We analyzed our surveys with traditional single-trait genome-wide association analysis (GWAS), multitrait GWAS, and functional networks built from a diverse set of plant phenotypes. Plant genotype was influential in structuring arthropod community composition among several garden sites. Candidate genes important for higher level organization of arthropod communities had broadly applicable functions, such as terpenoid biosynthesis and production of dsRNA binding proteins and protein kinases, which may be capable of targeting multiple arthropod species. We have demonstrated the ability to detect, in an uncontrolled environment, individual genes that are associated with the community assemblage of arthropods on a host plant, further enhancing our understanding of genetic mechanisms that impact ecosystem structure.

Community genetics↗
DOE OSTI2026

Assessing the Application of a Genomic Network Analysis in Population Ecology: Inferring Patterns of Dispersal and Geographic Structure in the Emerging Pathogen, Coccidioides

A challenge in population ecology studies is identifying how to best group individuals into populations, especially when individual origin is unknown. Machine learning has improved upon traditional methods of identifying population structure and is more efficient at handling large, complex datasets. We demonstrate the applicability of a machine learning method to identify hierarchical population structure in an emerging pathogen, Coccidioides spp., the causative agent of Valley fever. We compared the network clusters to structure identified by traditional tools as a validation of the network performance. We used publicly available whole-genome data for 48 C. immitis and 102 C. posadasii, resulting in 168,211 genome-wide SNPs among the two species. The network analysis grouped samples into populations comparable to the literature for these species but also identified fine-scale geographic structure and travel-associated cases not reported thus far. Exploring different resolutions in the network made it easy to identify unique genotypes specific to California and possibly Nevada, as well as Phoenix- and Tucson-acquired infections in non-endemic areas, regardless of reported travel history. The present study provides a promising example of how a ML-based network analysis can improve our ability to understand pathogen ecology, group cases into populations and infer travel-associated infections.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2024

Metabolic Stress in Space: ROS-Induced Mutations in Mice Hint at A New Path to Cancer

Long-duration spaceflight beyond Earth's magnetosphere poses serious health risks, including muscle atrophy, bone loss, liver and kidney damage, and the Spaceflight-Associated Neuro-ocular Syndrome (SANS). RNA-seq of mice aboard the International Space Station (ISS) for 37 days revealed extraordinary hypermutation in tissue-specific genes, with guanine base conversion predominating, potentially contributing to spaceflight-associated health risks. Our results suggest that the genome-wide accelerated mutation that we measured, seemingly independent of radiation dose, was induced by oxidative damage from higher atmospheric carbon dioxide (CO 2 ) levels and increased reactive oxygen species (ROS) on the ISS. This accelerated mutation, faster via RNA transcription than replication and more numerous than by radiation alone, unveils novel hotspots in the mammalian proteome. Notably, these hotspots correlate with commonly mutated genes across various human cancers, highlighting the ISS as a crucial platform for studying accelerated mutation, genome instability, and the induction of disease-causing mutations in model organisms. Our results suggest that metabolic processes can contribute to somatic mutation, and thus may play a role in the development of cancer. A metabolic link to genetic instability potentially has far-reaching implications for various diseases, with implications for human health on Earth and in space.

mice↗
DOE OSTI2025

Editorial: Transcriptional and epigenetic landscapes of abiotic stress response in plants

In nature, plants constantly face various biotic and abiotic stresses that impact their growth, development, and productivity. Among these, abiotic stresses often have a more severe impact than biotic stresses. For instance, drought has been reported to cause greater yield losses than the combined impact of all plant pathogens (Gupta et al., 2020). Abiotic stresses are the immediate outcome of climate change, and the magnitude of these stresses has gradually increased every year with the rise in global temperatures. Thus, it has become imperative to study the impact of these stresses on plants and how plants respond to them at different levels to show resilient traits. This includes analysing the plants at morpho-physiological, biochemical, and molecular levels. Researchers often compare stressed plants to control (non-stressed) plants or evaluate contrasting genotypes, such as tolerant and sensitive lines, to elucidate the mechanisms underlying stress responses. While these studies have provided some insights, a comprehensive understanding of the intricate mechanisms governing plant responses to abiotic stress remains largely unknown. Recent advances in next-generation tools and technologies have enabled researchers to dissect the molecular basis of plant stress responses at genomic, transcriptomic, proteomic, metabolomic, epigenetic and epigenomic levels. Among these, knowledge of the transcriptional/epigenomic landscape of the trait-associated variations is limited. Given the importance of transcriptional changes and histone modifications in abiotic stress responses, this Research Topic was edited to collage the knowledge available on transcriptional and epigenetic landscapes of abiotic stress response in plants. The Research Topic features eight original research articles and one review, covering various aspects of transcriptome and epigenetic reprogramming in plants during abiotic stresses. Four of the research articles employ transcriptomics integrated with other omics approaches to explore transcriptome reprogramming, candidate gene identification, and the role of long non-coding RNA during different stresses. Two articles focus on the functional characterization of specific candidate genes involved in stress response, while another provides a genome-wide analysis of a stress-responsive gene family. Additionally, one study investigates genome-wide histone modifications, specifically H3K4me3 and H3K27me3, in response to abiotic stresses.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2019

Prenatal 2g Exposure Alters Placental Expression of Stress-Related Genes

Extra-terrestrial colonization is of growing interest to space agencies and private entities, emphasizing the importance of research on reproduction and development in the absence of Earth's 1G. Maternal stressors can modify offspring development, exerting significant lifespan and crossgenerational changes through prenatal programming. The space environment is stressful, therefore exposure to altered gravity during pregnancy may impact later life outcomes in offspring. In ground-based studies, we exposed pregnant rats to continuous +G (above Earth gravity), and observed overweight and elevated anxiety in adult male (but not female) offspring, common phenotypes associated with prenatal maternal stress. Here we hypothesize that exposure to increased gravity during pregnancy elicits changes in the expression of stress-related genes in placenta that may mediate emergence of later life outcomes. While the placenta transports maternal factors to the fetus and produces endogenous fetal hormones, stress-induced changes at the placental-uterine interface may also alter communication between mother and fetus, facilitating prenatal transmission of unfavorable later life outcomes and cross-generational epigenetic alterations. Maternal stress elevates maternal glucocorticoids however placental 11b-hydroxysteroid dehydrogenase type 2 (HSD11B2) buffers fetal exposure by converting cortisol/corticosterone into inactive metabolites. Maternal stress during pregnancy down-regulates this enzyme and can induce epigenetic changes in placental and fetal tissues accounting for heightened adult HPA reactivity. Past studies have shown a placenta-specific increase in DNA methyltransferase (DNMT3a) mRNA in stressed mothers, an effect with implications for genome-wide epigenetic changes that may account for diverse phenotypic outcomes following maternal stress. Here we exposed groups of pregnant rats to one of five gravity loads (1, 1.5, 1.75 and 2G) and analyzed placental samples during late gestation. We predicted a systematic dose-response relationship between gravity load and the expression of the HSD11B2 and DNMT3 genes, thereby linking maternal exposure to altered gravity during pregancy with maternal stress.

Kalotia, Simranjit K.↗
DOE OSTI2024

Variants in the MS4A cluster interact with soluble TREM2 expression on biomarkers of neuropathology

Recent evidence suggests that Alzheimer’s disease (AD) genetic risk variants (rs1582763 and rs6591561) of the MS4A locus are genome-wide significant regulators of soluble TREM2 levels such that the minor allele of the protective variant (rs1582763) is associated with higher sTREM2 and lower AD risk while the minor allele of (rs6591561) relates to lower sTREM2 and higher AD risk. Our group previously found that higher sTREM2 relates to higher Aβ 40 , worse blood–brain barrier (BBB) integrity (measured with the CSF/plasma albumin ratio), and higher CSF tau, suggesting strong associations with amyloid abundance and both BBB and neurodegeneration complicate interpretation. We expand on this work by leveraging these common variants as genetic tools to tune the interpretation of high CSF sTREM2, and by exploring the potential modifying role of these variants on the well-established associations between CSF sTREM2 as well as TREM2 transcript levels in the brain with AD neuropathology. Biomarker analyses leveraged data from the Vanderbilt Memory & Aging Project (n = 127, age = 72 ± 6.43) and were replicated in the Alzheimer’s Disease Neuroimaging Initiative (n = 399, age = 73 ± 7.39). Autopsy analyses were performed leveraging data from the Religious Orders Study and Rush Memory and Aging Project (n= 577, age = 89 ± 6.46). We found that the protective variant rs1582763 attenuated the association between CSF sTREM2 and Aβ 40 (β= -0.44, p-value= 0.017) and replicated this interaction in ADNI (β = -0.27, p = 0.017). We did not observe this same interaction effect between TREM2 mRNA levels and Aβ peptides in brain (Aβ total β = -0.14, p = 0.629; Aβ 1-38 , β = 0.11, p = 0.200). In contrast to the effects on Aβ, the minor allele of this same variant seemed to enhance the association with blood–brain barrier dysfunction (β = 7.0e-4, p = 0.009), suggesting that elevated sTREM2 may carry a much different interpretation in carriers vs. non-carriers of this allele. When evaluating the risk variant (rs6591561) across datasets, we did not observe a statistically significant interaction against any outcome in VMAP and observed opposing directions of associations in ADNI and ROS/MAP on Aβ levels. Together, our results suggest that the protective effect of rs1582763 may act by decoupling the associations between sTREM2 and amyloid abundance, providing important mechanistic insight into sTREM2 changes and highlighting the need to incorporate genetic context into the analysis of sTREM2 levels, particularly if leveraged as a clinical biomarker of disease in the future.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Impact of genotype-calling methodologies on genome-wide association and genomic prediction in polyploids

This dataset contains all data used in the paper "Impact of genotype-calling methodologies on genome-wide association and genomic prediction in polyploids". The dataset includes genotypes and phenotypic data from two autotetraploid species Miscanthus sacchariflorus and Vaccinium corymbosum that was used used for genome wide association studies and genomic prediction and the scripts used in the analysis. In this V2, 2 files have the raw data are added: "Miscanthus_sacchariflorus_RADSeq.vcf" is the VCF file with the raw SNP calls of the Miscanthus sacchariflorus data used for genotype calling using the 6 genotype calling methods. "Blueberry_data_read_depths.RData" is the a RData file with the read depth data that was used for genotype calling in the Blueberry dataset.

allelic dosage↗
DOE OSTI2023

Genome-wide fitness profiling reveals molecular mechanisms that bacteria use to interact with Trichoderma atroviride exometabolites

Trichoderma spp. are ubiquitous rhizosphere fungi capable of producing several classes of secondary metabolites that can modify the dynamics of the plant-associated microbiome. However, the bacterial-fungal mechanisms that mediate these interactions have not been fully characterized. Here, a random barcode transposon-site sequencing (RB-TnSeq) approach was employed to identify bacterial genes important for fitness in the presence of Trichoderma atroviride exudates. We selected three rhizosphere bacteria with RB-TnSeq mutant libraries that can promote plant growth: the nitrogen fixers Klebsiella michiganensis M5aI and Herbaspirillum seropedicae SmR1, and Pseudomonas simiae WCS417. As a non-rhizosphere species, Pseudomonas putida KT2440 was also included. From the RB-TnSeq data, nitrogen-fixing bacteria competed mainly for iron and required the siderophore transport system TonB/ExbB for optimal fitness in the presence of T . atroviride exudates. In contrast, P . simiae and P . putida were highly dependent on mechanisms associated with membrane lipid modification that are required for resistance to cationic antimicrobial peptides (CAMPs). A mutant in the Hog1-MAP kinase (Δ tmk 3) gene of T . atroviride showed altered expression patterns of many nonribosomal peptide synthetase (NRPS) biosynthetic gene clusters with potential antibiotic activity. In contrast to exudates from wild-type T . atroviride , bacterial mutants containing lesions in genes associated with resistance to antibiotics did not show fitness defects when RB-TnSeq libraries were exposed to exudates from the Δ tmk3 mutant. Unexpectedly, exudates from wild-type T . atroviride and the Δ tmk 3 mutant rescued purine auxotrophic mutants of H . seropedicae , K . michiganensis and P . simiae . Metabolomic analysis on exudates from wild-type T . atroviride and the Δ tmk 3 mutant showed that both strains excrete purines and complex metabolites; functional Tmk3 is required to produce some of these metabolites. This study highlights the complex interplay between Trichoderma -metabolites and soil bacteria, revealing both beneficial and antagonistic effects, and underscoring the intricate and multifaceted nature of this relationship.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Genome‐wide association identifies a BAHD acyltransferase activity that assembles an ester of glucuronosylglycerol and phenylacetic acid

SUMMARY Genome‐wide association studies (GWAS) are an effective approach to identify new specialized metabolites and the genes involved in their biosynthesis and regulation. In this study, GWAS of Arabidopsis thaliana soluble leaf and stem metabolites identified alleles of an uncharacterized BAHD‐family acyltransferase (AT5G57840) associated with natural variation in three structurally related metabolites. These metabolites were esters of glucuronosylglycerol, with one metabolite containing phenylacetic acid as the acyl component of the ester. Knockout and overexpression of AT5G57840 in Arabidopsis and heterologous overexpression in Nicotiana benthamiana and Escherichia coli demonstrated that it is capable of utilizing phenylacetyl‐CoA as an acyl donor and glucuronosylglycerol as an acyl acceptor. We, thus, named the protein Glucuronosylglycerol Ester Synthase (GGES). Additionally, phenylacetyl glucuronosylglycerol increased in Arabidopsis CYP79A2 mutants that overproduce phenylacetic acid and was lost in knockout mutants of UDP‐sulfoquinovosyl: diacylglycerol sulfoquinovosyl transferase, an enzyme required for glucuronosylglycerol biosynthesis and associated with glycerolipid metabolism under phosphate‐starvation stress. GGES is a member of a well‐supported clade of BAHD family acyltransferases that arose by duplication and neofunctionalized during the evolution of the Brassicales within a larger clade that includes HCT as well as enzymes that synthesize other plant‐specialized metabolites. Together, this work extends our understanding of the catalytic diversity of BAHD acyltransferases and uncovers a pathway that involves contributions from both phenylalanine and lipid metabolism.

09 BIOMASS FUELS↗
DOE OSTI2023

Genome‐wide association and genomic prediction for yield and component traits of Miscanthus sacchariflorus

Abstract Accelerating biomass improvement is a major goal of Miscanthus breeding. The development and implementation of genomic‐enabled breeding tools, like marker‐assisted selection (MAS) and genomic selection, has the potential to improve the efficiency of Miscanthus breeding. The present study conducted genome‐wide association (GWA) and genomic prediction of biomass yield and 14 yield‐components traits in Miscanthus sacchariflorus . We evaluated a diversity panel with 590 accessions of M. sacchariflorus grown across 4 years in one subtropical and three temperate locations and genotyped with 268,109 single‐nucleotide polymorphisms (SNPs). The GWA study identified a total of 835 significant SNPs and 674 candidate genes across all traits and locations. Of the significant SNPs identified, 280 were localized in mapped quantitative trait loci intervals and proximal to SNPs identified for similar traits in previously reported Miscanthus studies, providing additional support for the importance of these genomic regions for biomass yield. Our study gave insights into the genetic basis for yield‐component traits in M. sacchariflorus that may facilitate marker‐assisted breeding for biomass yield. Genomic prediction accuracy for the yield‐related traits ranged from 0.15 to 0.52 across all locations and genetic groups. Prediction accuracies within the six genetic groupings of M. sacchariflorus were limited due to low sample sizes. Nevertheless, the Korea/NE China/Russia ( N = 237) genetic group had the highest prediction accuracy of all genetic groups (ranging 0.26–0.71), suggesting that with adequate sample sizes, there is strong potential for genomic selection within the genetic groupings of M. sacchariflorus . This study indicated that MAS and genomic prediction will likely be beneficial for conducting population‐improvement of M. sacchariflorus .

09 BIOMASS FUELS↗
DOE OSTI2026

Identification of proteins influencing CRISPR-associated transposases for enhanced genome editing

CRISPR-associated transposases (CASTs) hold tremendous potential for microbial genome editing because of their ability to integrate large DNA cargos in a programmable, site-specific manner. However, their widespread application has been hindered by poorly understood host factor requirements for transposition. To address this gap, we conducted the first genome-wide screen for host factors affecting Vibrio cholerae CAST (VchCAST) activity using an Escherichia coli RB-TnSeq library and identified 15 genes affecting VchCAST transposition. Of these, seven factors were validated to improve VchCAST activity, and two were inhibitory. Guided by the identification of homologous recombination effectors, RecD and RecA, we tested the λ-Red recombineering system in our VchCAST editing vectors and increased editing efficiency by 55.2-fold in E. coli, 5.6-fold in Pseudomonas putida, and 10.8-fold in Klebsiella michiganensis while maintaining high target specificity and similar insertion arrangements. This study improves the understanding of factors affecting VchCAST activity and enhances its efficiency as a bacterial genome editor.

Song, Leo C T↗
DOE OSTI2024

Learning epistatic polygenic phenotypes with Boolean interactions

Detecting epistatic drivers of human phenotypes is a considerable challenge. Traditional approaches use regression to sequentially test multiplicative interaction terms involving pairs of genetic variants. For higher-order interactions and genome-wide large-scale data, this strategy is computationally intractable. Moreover, multiplicative terms used in regression modeling may not capture the form of biological interactions. Building on the Predictability, Computability, Stability (PCS) framework, we introduce the epiTree pipeline to extract higher-order interactions from genomic data using tree-based models. The epiTree pipeline first selects a set of variants derived from tissue-specific estimates of gene expression. Next, it uses iterative random forests (iRF) to search training data for candidate Boolean interactions (pairwise and higher-order). We derive significance tests for interactions, based on a stabilized likelihood ratio test, by simulating Boolean tree-structured null (no epistasis) and alternative (epistasis) distributions on hold-out test data. Finally, our pipeline computes PCS epistasis p-values that probabilisticly quantify improvement in prediction accuracy via bootstrap sampling on the test set. We validate the epiTree pipeline in two case studies using data from the UK Biobank: predicting red hair and multiple sclerosis (MS). In the case of predicting red hair, epiTree recovers known epistatic interactions surrounding MC1R and novel interactions, representing non-linearities not captured by logistic regression models. In the case of predicting MS, a more complex phenotype than red hair, epiTree rankings prioritize novel interactions surrounding HLA-DRB1 , a variant previously associated with MS in several populations. Taken together, these results highlight the potential for epiTree rankings to help reduce the design space for follow up experiments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Improved Biofuel Production through Discovery and Engineering of Terpene Metabolism in Switchgrass

Project Objectives - Of the myriad specialized metabolites that plants deploy to adapt to environmental challenges, terpenes form the largest group. In many major crops, unique terpene blends serve as key stress defenses that directly impact plant fitness and yield. In addition, terpenes, such as bisabolene and pinene, are used for producing renewable biofuels. Essential to advancing a broader use of terpenes for biofuel feedstock engineering is a system-wide knowledge of the diverse biosynthetic machinery and defensive potential of often species-specific terpene blends. The proposed project would merge genome-wide enzyme discovery with comparative –omics, protein structural and plant microbiome studies to define the biosynthesis and stress-defensive functions of the switchgrass (Panicum virgatum) terpene network. These insights would be combined with developing and applying non-transgenic genome editing tools to design plants with desirable terpene blends for higher productivity and biofuel production on marginal lands. As a dedicated lignocellulosic feedstock for U.S. biofuel production with high net energy yield, stress tolerance, and available genome resources, switchgrass is well-suited for devising new avenues for biofuel production. Project Description – The diversity of plant terpene defenses is governed by species-specific families of terpene synthase (TPS) and cytochrome P450 monooxygenase (P450) enzymes. Mining of the switchgrass genome (genotype Alamo) identified ~100 TPS and P450 candidate genes, and combinatorial biochemical analysis of synthesized TPSs and P450s revealed more than a dozen enzymes with common and novel activities. In addition, several identified terpene metabolites and the corresponding transcripts were up-regulated in response to abiotic stressors. These findings demonstrate a unique switchgrass terpene network with probable importance to abiotic stress tolerance, thus providing a large chemical portfolio for optimizing crop resistance, yield, and biofuel composition. Leveraging these preliminary data, we propose to generate a genome-wide map of the switchgrass terpene metabolic network through multi-gene co-expression analyses that allow the efficient cross-validation of TPS and P450 functions. Key enzymes would further be applied to structure-function studies via X-ray protein crystallography, homology modeling and site-directed mutagenesis to gain mechanistic insight into the catalytic specificity of switchgrass terpene metabolism and provide gene and amino acid targets for genome editing. In tandem with terpene pathway discovery, system-wide metabolomics, transcriptomics and proteomics studies in switchgrass accessions of contrasting drought tolerance would define the role of switchgrass terpene metabolism in conferring abiotic stress resilience. Metabolic changes would be assessed in a combined approach of targeted (terpenes) and untargeted metabolite profiling using a high-resolution LC-MS/MS approach, differential gene expression analyses through multiplexed Illumina RNA sequencing, and quantitative analysis of high-priority pathway enzymes using multiple reaction monitoring (MRM). Drawing on these insights, knock-down/out mutants of stress-associated pathway nodes would be generated by optimizing transient virus-induced gene silencing (VIGS) and CRISPR/Cas9 systems under control of the Tobacco Rattle Virus (TRV). The resulting mutant lines would then be analyzed for stress susceptibility and the impact on the root microbiome to define gene functions in planta. Knowledge of terpene pathways, enzyme mechanisms and bioactivities would be applied to enhance switchgrass stress resilience and to tailor-make terpene blends for biofuel production. Here, TRV-enabled CRISPR/Cas9 genome editing, including allele-specific knock-out of redundant genes, engineering of enzyme specificity via structure-guided point mutations, and overexpression of terpene genes relevant to stress-protection or biofuel production, would be used to increase metabolic flux toward desired pathways. Broader Impacts - Integrating the system-wide discovery, mechanistic analysis and non-transgenic genome engineering of the switchgrass terpene network aligns the required steps to unlock the chemical potential of this important metabolite class to generate crops that are more resistant to stress and provide advanced biofuel production in light of rising climate pressures as foreseeable challenges for bioenergy crop cultivation. The proposed project would further offer interdisciplinary student training through active involvement in the project and integration of research concepts and outcomes into newly-developed graduate and undergraduate courses on Plant Biotechnology.

09 BIOMASS FUELS↗
DOE OSTI2024

Climate adaptation and sustainability in switchgrass: exploring plant-microbe-soil interactions across continental scale environmental gradients

Less carbon-intensive energy sources are needed to reduce greenhouse gas emissions and their predicted role in climate change. There is growing interest in the potential of biofuels for meeting this need. A critical question is whether large-scale biofuel production can be sustainable over the time scales needed to mitigate our carbon debt from fossil fuel consumption. The carbon balance and ultimately the sustainability of biofuel feedstock production is the result of complex climate-coupled interactions between carbon fixation, sequestration, and release through combustion. Similarly, the long-term productivity of biofuels depends on the environmental factors limiting plant growth. These factors are often related to soil resources which involve complex interactions at the plant-microbe-soil interface impacting their availability and cycling. Our collaborative project addressed sustainable switchgrass (Panicum virgatum) production by exploring Plant Systems, Plant-Microbiome Interactions, and Ecosystem Processes through the integrating lens of Multi-Scale Modeling. Our research was based on detailed characterization of genetically diverse switchgrass genotypes planted in common gardens across a continental latitudinal gradient. The underlying theme of our Plant Systems research was the use of locally adapted plant material to explore plant function, to understand the mechanistic basis of environmental interactions, and to discover the plant genes important for adaptation and sustainability in the face of climate change. Our Plant-Microbiome Interaction project characterized the microbial communities associated with switchgrass using genomic tools. Our Ecosystem Processes research focused on carbon cycle responses at the ecosystem level using stand level plantings. Finally, our Multi-Scale Modeling helped to define conditions of a sustainable biofuel system and identify key tradeoffs between genetic diversity, productivity, and ecosystem services. Genome-wide association analyses were used to identify alleles that contribute to successful establishment and biomass production across North America. Together, our work provided a baseline analyses of the potential of switchgrass as a biofuel feedstock. Our project resulted in a number of successful outcomes. First, we were successful in collecting switchgrass germplasm across the species range, propagating the material, and establishing common garden experiments across the species range. In collaboration with DOE JGI, we successfully assembled the first tetraploid switchgrass genome and published this resource with an analyses of the genetic basis local adaptation from our gardens (Lowry et al. 2019, Lovell et al. 2021). The gardens were used to characterize the genetic architecture for a number of important plant phenotypes. Our project also conducted extensive sampling and sequencing to characterize the bacterial and fungal associates of switchgrass roots and leaves. We showed that host genotype, location, and harvesting practices can play a role in microbiome assembly (Singer et al. 2019 & 2022, Van Wallendael et al. 2020 & 2022, Edwards et al. 2023). Our ecosystem processes work created baseline dataset of carbon and nutrient cycling in realistic stand plantings of switchgrass. Data from this experiment provided new insight into the role of plant traits, phenology, and local environments in ecosystem processes like soil respiration, net-ecosystem exchange, and dynamics of soil and plant nutrients (Ricketts et al. 2023). Finally, our crop modelling experiments help to characterize the sensitivity of common modeling frameworks to parameters, identify key limiters of productivity across large geographic scales, and leverage patterns of local adaptation in prediction. Ultimately, these studies help to identify critical plant-microbe-soil traits that may be manipulated, through breeding or agronomic management, to improve the sustainability of biofuel feedstocks.

09 BIOMASS FUELS↗