Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “genomic mutations”

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 199 records · Page 11

AlphaBeta: computational inference of epimutation rates and spectra from high-throughput DNA methylation data in plants

Stochastic changes in DNA methylation (i.e., spontaneous epimutations) contribute to methylome diversity in plants. Here, we describe AlphaBeta, a computational method for estimating the precise rate of such stochastic events using pedigree-based DNA methylation data as input. We demonstrate how AlphaBeta can be employed to study transgenerationally heritable epimutations in clonal or sexually derived mutation accumulation lines, as well as somatic epimutations in long-lived perennials. Application of our method to published and new data reveals that spontaneous epimutations accumulate neutrally at the genome-wide scale, originate mainly during somatic development and that they can be used as a molecular clock for age-dating trees.

59 BASIC BIOLOGICAL SCIENCES↗

Comparative genomics of the Liberibacter genus reveals widespread diversity in genomic content and positive selection history

‘Candidatus Liberibacter’ is a group of bacterial species that are obligate intracellular plant pathogens and cause Huanglongbing disease of citrus trees and Zebra Chip in potatoes. Here, we examined the extent of intra- and interspecific genetic diversity across the genus using comparative genomics. Our approach examined a wide set of Liberibacter genome sequences including five pathogenic species and one species not known to cause disease. By performing comparative genomics analyses, we sought to understand the evolutionary history of this genus and to identify genes or genome regions that may affect pathogenicity. With a set of 52 genomes, we performed comparative genomics, measured genome rearrangement, and completed statistical tests of positive selection. We explored markers of genetic diversity across the genus, such as average nucleotide identity across the whole genome. These analyses revealed the highest intraspecific diversity amongst the ‘Ca. Liberibacter solanacearum’ species, which also has the largest plant host range. We identified sets of core and accessory genes across the genus and within each species and measured the ratio of nonsynonymous to synonymous mutations (dN/dS) across genes. We identified ten genes with evidence of a history of positive selection in the Liberibacter genus, including genes in the Tad complex, which have been previously implicated as being highly divergent in the ‘Ca. L. capsica’ species based on high values of dN.

59 BASIC BIOLOGICAL SCIENCES↗

Haploidy and aneuploidy in switchgrass mediated by misexpression of CENH3

Cross bred species such as switchgrass may benefit from advantageous breeding strategies requiring inbred lines. Doubled haploid production methods offer several ways that these lines can be produced that often involve uniparental genome elimination as the rate limiting step. We have used a centromere-mediated genome elimination strategy in which modified CENH3 is expressed to induce the process. Transgenic tetraploid switchgrass lines coexpressed Cas9, a poly-cistronic tRNA-gRNA tandem array containing eight guide RNAs that target two CENH3 genes, and different chimeric versions of CENH3 with alterations to the N-terminal tail region. Genotyping of CENH3 genes in transgenics identified edits including frameshift mutations and deletions in one or both copies of the two CENH3 genes. Flow cytometry of T 1 seedlings identified two T 0 lines that produced five haploid individuals representing an induction rate of 0.5% and 1.4%. Eight different T0 lines produced aneuploids at rates ranging from 2.1 to 14.6%. A sample of aneuploid lines were sequenced at low coverage and aligned to the reference genome, revealing missing chromosomes and chromosome arms.

59 BASIC BIOLOGICAL SCIENCES↗

Pressurized DNA state inside herpes capsids—A novel antiviral target

Drug resistance in viruses represents one of the major challenges of healthcare. As part of an effort to provide a treatment that avoids the possibility of drug resistance, we discovered a novel mechanism of action (MOA) and specific compounds to treat all nine human herpesviruses and animal herpesviruses. The novel MOA targets the pressurized genome state in a viral capsid, “turns off” capsid pressure, and blocks viral genome ejection into a cell nucleus, preventing viral replication. This work serves as a proof-of-concept to demonstrate the feasibility of a new antiviral target—suppressing pressure-driven viral genome ejection —that is likely impervious to developing drug resistance. This pivotal finding presents a platform for discovery of a new class of broad-spectrum treatments for herpesviruses and other viral infections with genome-pressure-dependent replication. A biophysical approach to antiviral treatment such as this is also a vital strategy to prevent the spread of emerging viruses where vaccine development is challenged by high mutation rates or other evasion mechanisms.

59 BASIC BIOLOGICAL SCIENCES↗

Identification of a broadly conserved family of enzymes that hydrolyze (p)ppApp

Bacteria produce a variety of nucleotide second messengers to adapt to their surroundings. Although chemically similar, the nucleotides guanosine penta- and tetraphosphate [(p)ppGpp] and adenosine penta- and tetraphosphate [(p)ppApp] have distinct functions in bacteria. (p)ppGpp mediates survival under nutrient-limiting conditions and its intracellular levels are regulated by synthetases and hydrolases belonging to the RelA-SpoT homolog (RSH) family of enzymes. By contrast, (p)ppApp is not known to be involved in nutrient stress responses and is synthesized by RSH-resembling toxins that inhibit the growth of bacterial cells. However, it remains unclear whether there exists a family of hydrolases that specifically act on (p)ppApp to reverse its toxic effects. Here, we present the structure and biochemical characterization of adenosine 3'-pyrophosphohydrolase 1 (Aph1), the founding member of a monofunctional (p)ppApp hydrolase family of enzymes. Our work reveals that Aph1 adopts a histidine-aspartate (HD)-domain fold characteristic of phosphohydrolase metalloenzymes and its activity mitigates the growth inhibitory effects of (p)ppApp-synthesizing toxins. Using an informatic approach, we identify over 2,000 putative (p)ppApp hydrolases that are widely distributed across bacterial phyla and found in diverse genomic contexts, and we demonstrate that 12 representative members hydrolyze ppApp. In addition, our in silico analyses reveal a unique molecular signature that is specific to (p)ppApp hydrolases, and we show that mutation of two residues within this signature broadens the specificity of Aph1 to promiscuously hydrolyze (p)ppGpp in vitro. Overall, our findings indicate that like (p)ppGpp hydrolases, (p)ppApp hydrolases are widespread in bacteria and may play important and underappreciated role(s) in bacterial physiology.

(p)ppApp↗

MAL33 drives natural variation in maltose metabolism in Saccharomyces eubayanus

Maltose is one of the most abundant sugars in brewer’s wort, and its efficient utilization is critical for successful fermentation. However, maltose consumption varies naturally among Saccharomyces eubayanus strains isolated from different host trees, such as Quercus and Nothofagus. To identify the genetic determinants underlying these phenotypic differences, we performed bulk segregant analysis (BSA) and quantitative trait loci (QTL) mapping using an F 2 offspring derived from QC18 (Quercus-associated) and CL467.1 (Nothofagus-associated) strains. QTL mapping identified two significant genomic regions on subtelomeric loci of chromosomes V-R and XVI-L, each containing complete MAL loci composed of MAL32 (encoding maltase), MAL31 (transporter), and MAL33 (transcriptional activator) genes. Comparative polymorphism analyses identified mutations in MAL32 and MAL33 of QC18, including frameshift mutations resulting in premature stop codons. Functional validation demonstrated that the heterologous expression of MAL33 ChrV from CL467.1 fully restored maltose utilization in QC18, indicating the functional presence of MAL33 cis-regulatory sequences and MAL32 and MAL31 genes in QC18. While structural protein predictions identified truncation and impaired functionality in the maltose-responsive activation domain of Mal33p from QC18, overexpression of QC18’s own MAL33 ChrV allele also improved maltose metabolism, suggesting dosage-dependent transcriptional limitations rather than complete functional loss. These results indicate that allelic variations in the maltose-responsive activation domain of Mal33p result in differences in maltose consumption between strains. Here, we hypothesized that reduced maltose metabolism in QC18 is an adaptive response to the distinct sugar composition in Quercus robur bark, contrasting with the starch-rich environment of Nothofagus pumilio. These findings highlight subtelomeric MAL gene diversity as a reservoir of genetic variation, representing a key evolutionary mechanism that influences maltose adaptation among natural Saccharomyces isolates.

evolutionary plasticity↗

Proteogenomic characterization of difficult-to-treat breast cancer with tumor cells enriched through laser microdissection

Abstract Background Breast cancer (BC) is the most commonly diagnosed cancer and the leading cause of cancer death among women globally. Despite advances, there is considerable variation in clinical outcomes for patients with non-luminal A tumors, classified as difficult-to-treat breast cancers (DTBC). This study aims to delineate the proteogenomic landscape of DTBC tumors compared to luminal A (LumA) tumors. Methods We retrospectively collected a total of 117 untreated primary breast tumor specimens, focusing on DTBC subtypes. Breast tumors were processed by laser microdissection (LMD) to enrich tumor cells. DNA, RNA, and protein were simultaneously extracted from each tumor preparation, followed by whole genome sequencing, paired-end RNA sequencing, global proteomics and phosphoproteomics. Differential feature analysis, pathway analysis and survival analysis were performed to better understand DTBC and investigate biomarkers. Results We observed distinct variations in gene mutations, structural variations, and chromosomal alterations between DTBC and LumA breast tumors. DTBC tumors predominantly had more mutations inTP53,PLXNB3, Zinc finger genes, and fewer mutations inSDC2,CDH1,PIK3CA,SVIL, andPTEN. Notably, Cytoband 1q21, which contains numerous cell proliferation-related genes, was significantly amplified in the DTBC tumors. LMD successfully minimized stromal components and increased RNA–protein concordance, as evidenced by stromal score comparisons and proteomic analysis. Distinct DTBC and LumA-enriched clusters were observed by proteomic and phosphoproteomic clustering analysis, some with survival differences. Phosphoproteomics identified two distinct phosphoproteomic profiles for high relapse-risk and low relapse-risk basal-like tumors, involving several genes known to be associated with breast cancer oncogenesis and progression, includingKIAA1522,DCK,FOXO3,MYO9B,ARID1A,EPRS,ZC3HAV1, andRBM14. Lastly, an integrated pathway analysis of multi-omics data highlighted a robust enrichment of proliferation pathways in DTBC tumors. Conclusions This study provides an integrated proteogenomic characterization of DTBC vs LumA with tumor cells enriched through laser microdissection. We identified many common features of DTBC tumors and the phosphopeptides that could serve as potential biomarkers for high/low relapse-risk basal-like BC and possibly guide treatment selections.

Oncology↗

Deep Sequencing of MHC-Adapted Viral Lines Reveals Complex Recombinational Exchanges With Endogenous Retroviruses Leading to High-Frequency Variants

Experimental evolution (serial passage) of Friend virus complex (FVC) in mice demonstrates phenotypic adaptation to specific host major histocompatibility complex (MHC) genotypes. These evolved viral lines show increased fitness and virulence in their host-genotype-of-passage, but display fitness and virulence tradeoffs when infecting unfamiliar host MHC genotypes. Here, we deep sequence these viral lines in an attempt to discover the genetic basis of FVC adaptation. The principal prediction for genotype-specific adaptation is that unique mutations would rise to high frequency in viral lines adapted to each host MHC genotype. This prediction was not supported by our sequencing data as most observed high-frequency variants were present in each of our independently evolved viral lines. However, using a multi-variate approach to measure divergence between viral populations, we show that populations of replicate evolved viral lines from the same MHC congenic mouse strain were more similar to one another than to lines derived from different MHC congenic mouse strains, suggesting that MHC genotype does predictably act on viral evolution in our model. Sequence analysis also revealed rampant recombination with endogenous murine leukemia virus sequences (EnMuLVs) that are encoded within the BALB/c mouse genome. The highest frequency variants in all six lines contained a 12 bp insertion from a recombinant EnMuLV source, suggesting such recombinants were either being favored by selection or were contained in a recombinational hotspot. Interestingly, they did not reach fixation, as if they are low fitness. The amount of background mutations linked to FVC/EnMuLV variable sites indicated that FVC/EnMuLV recombinants had not reached mutation selection equilibrium and thus, that EnMuLV sequences are likely continuously introgressing into the replicating viral population. These discoveries raise the question: is the expression of EnMuLV sequences in mouse splenocytes that permit recombination with exogenous FVC a pathogen or host adaptation?

59 BASIC BIOLOGICAL SCIENCES↗

A view of the pan‐genome of domesticated Cowpea ( Vigna unguiculata [L.] Walp.)

Abstract Cowpea, Vigna unguiculata L . Walp., is a diploid warm‐season legume of critical importance as both food and fodder in sub‐Saharan Africa. This species is also grown in Northern Africa, Europe, Latin America, North America, and East to Southeast Asia. To capture the genomic diversity of domesticates of this important legume, de novo genome assemblies were produced for representatives of six subpopulations of cultivated cowpea identified previously from genotyping of several hundred diverse accessions. In the most complete assembly (IT97K‐499‐35), 26,026 core and 4963 noncore genes were identified, with 35,436 pan genes when considering all seven accessions. GO terms associated with response to stress and defense response were highly enriched among the noncore genes, while core genes were enriched in terms related to transcription factor activity, and transport and metabolic processes. Over 5 million single nucleotide polymorphisms (SNPs) relative to each assembly and over 40 structural variants >1 Mb in size were identified by comparing genomes. Vu10 was the chromosome with the highest frequency of SNPs, and Vu04 had the most structural variants. Noncore genes harbor a larger proportion of potentially disruptive variants than core genes, including missense, stop gain, and frameshift mutations; this suggests that noncore genes substantially contribute to diversity within domesticated cowpea.

59 BASIC BIOLOGICAL SCIENCES↗

A cell type-aware framework for nominating non-coding variants in Mendelian regulatory disorders

Abstract Unsolved Mendelian cases often lack obvious pathogenic coding variants, suggesting potential non-coding etiologies. Here, we present a single cell multi-omic framework integrating embryonic mouse chromatin accessibility, histone modification, and gene expression assays to discover cranial motor neuron (cMN)cis-regulatory elements and subsequently nominate candidate non-coding variants in the congenital cranial dysinnervation disorders (CCDDs), a set of Mendelian disorders altering cMN development. We generate single cell epigenomic profiles for ~86,000 cMNs and related cell types, identifying ~250,000 accessible regulatory elements with cognate gene predictions for ~145,000 putative enhancers. We evaluate enhancer activity for 59 elements using an in vivo transgenic assay and validate 44 (75%), demonstrating that single cell accessibility can be a strong predictor of enhancer activity. Applying our cMN atlas to 899 whole genome sequences from 270 genetically unsolved CCDD pedigrees, we achieve significant reduction in our variant search space and nominate candidate variants predicted to regulate known CCDD disease genesMAFB, PHOX2A, CHN1, andEBF3– as well as candidates in recurrently mutated enhancers through peak- and gene-centric allelic aggregation. This work delivers non-coding variant discoveries of relevance to CCDDs and a generalizable framework for nominating non-coding variants of potentially high functional impact in other Mendelian disorders.

Science & Technology - Other Topics↗

High allelic diversity in Arabidopsis NLRs is associated with distinct genomic features

Plants rely on Nucleotide-binding, Leucine-rich repeat Receptors (NLRs) for pathogen recognition. Highly variable NLRs (hvNLRs) show remarkable intraspecies diversity, while their low-variability paralogs (non-hvNLRs) are conserved between ecotypes. At a population level, hvNLRs provide new pathogen-recognition specificities, but the association between allelic diversity and genomic and epigenomic features has not been established. Our investigation of NLRs in Arabidopsis Col-0 has revealed that hvNLRs show higher expression, less gene body cytosine methylation, and closer proximity to transposable elements than non-hvNLRs. hvNLRs show elevated synonymous and nonsynonymous nucleotide diversity and are in chromatin states associated with an increased probability of mutation. Diversifying selection maintains variability at a subset of codons of hvNLRs, while purifying selection maintains conservation at non-hvNLRs. How these features are established and maintained, and whether they contribute to the observed diversity of hvNLRs is key to understanding the evolution of plant innate immune receptors.

59 BASIC BIOLOGICAL SCIENCES↗

Expanding the genetic toolkit: adenine and cytosine base editors for gene disruption in Aspergillus niger

Despite revolutionizing fungal genetic engineering, conventional CRISPR/Cas9-mediated knockouts rely on DNA double-strand breaks (DSBs), which can cause unwanted insertions and deletions, chromosomal abnormalities, and cytotoxicity. Base editors such as adenine base editors (ABEs), which convert A‧T to G‧C, and cytosine base editors (CBEs), which convert C‧G to T‧A, offer a safer alternative by enabling predictable, target-specific single-nucleotide changes without introducing DSBs. To overcome the limitations of traditional genome editing in filamentous fungi, we developed efficient base-editing systems in Aspergillus niger . For the first time, we constructed an ABE in A. niger , achieving up to 80% editing efficiency and inducing predictable A-to-G mutations at the intended intron sites, disrupting gene function through mRNA mis-splicing. We also developed a highly efficient CBE system, capable of introducing premature stop codons with 50–100% efficiency. To broaden the editing scope, we implemented a Cas9-NG variant recognizing a relaxed PAM sequence requiring only a single guanine (G), enabling editing at start codons and splice sites. Leveraging this expanded scope, we established gene disruption approaches by targeting start codons via ABE-mediated A-to-G conversions (ATG-to-GTG and ATG-to-ACG) and CBE-mediated C-to-T conversion (ATG-to-ATA). Additionally, our base-editing systems enable multiplex gRNA delivery and marker-free editing of multiple genes. Collectively, the scope-expanding strategies increase the number of genes targetable for disruption by base-editing in A. niger by 26.3% and enable near-complete coverage of 96% of the coding genes. Overall, this work demonstrates the potential of ABE and CBE systems as versatile, efficient, and safer alternatives to DSBs-based gene disruption in filamentous fungi.

Aspergillus↗

Proteomic and phosphoproteomic measurements enhance ability to predict ex vivo drug response in AML

Acute Myeloid Leukemia (AML) affects 20,000 patients in the US annually with a five-year survival rate of approximately 25%. One reason for the low survival rate is the high prevalence of clonal evolution that gives rise to heterogeneous sub-populations of leukemic cells with diverse mutation spectra, which eventually leads to disease relapse. This genetic heterogeneity drives the activation of complex signaling pathways that is reflected at the protein level. This diversity makes it difficult to treat AML with targeted therapy, requiring custom patient treatment protocols tailored to each individual’s leukemia. Toward this end, the Beat AML research program prospectively collected genomic and transcriptomic data from over 1000 AML patients and carried out ex vivo drug sensitivity assays to identify genomic signatures that could predict patient-specific drug responses. However, there are inherent weaknesses in using only genetic and transcriptomic measurements as surrogates of drug response, particularly the absence of direct information about phosphorylation-mediated signal transduction. As a member of the Clinical Proteomic Tumor Analysis Consortium, we have extended the molecular characterization of this cohort by collecting proteomic and phosphoproteomic measurements from a subset of these patient samples (38 in total) to evaluate the hypothesis that proteomic signatures can improve the ability to predict response to 26 drugs in AML ex vivo samples. In this work we describe our systematic, multi-omic approach to evaluate proteomic signatures of drug response and compare protein levels to other markers of drug response such as mutational patterns. We explore the nuances of this approach using two drugs that target key pathways activated in AML: quizartinib (FLT3) and trametinib (Ras/MEK), and show how patient-derived signatures can be interpreted biologically and validated in cell lines. In conclusion, this pilot study demonstrates strong promise for proteomics-based patient stratification to assess drug sensitivity in AML.

60 APPLIED LIFE SCIENCES↗

Comparative functional genomics identifies an iron-limited bottleneck in a Saccharomyces cerevisiae strain with a cytosolic-localized isobutanol pathway

Metabolic engineering strategies have been successfully implemented to improve the production of isobutanol, a next-generation biofuel, in Saccharomyces cerevisiae. Here, we explore how two of these strategies, pathway relocalization and redox cofactor-balancing, affect the performance and physiology of isobutanol producing strains. We equipped yeast with isobutanol cassettes which had either a mitochondrial or cytosolic localized isobutanol pathway and used either a redox-imbalanced (NADPH-dependent) or redox-balanced (NADHdependent) ketol-acid reductoisomerase enzyme. We then conducted transcriptomic, proteomic and metabolomic analyses to elucidate molecular differences between the engineered strains. Pathway localization had a large effect on isobutanol production with the strain expressing the mitochondrial-localized enzymes producing 3.8-fold more isobutanol than strains expressing the cytosolic enzymes. Cofactor-balancing did not improve isobutanol titers and instead the strain with the redox-imbalanced pathway produced 1.5-fold more isobutanol than the balanced version, albeit at low overall pathway flux. Functional genomic analyses suggested that the poor performances of the cytosolic pathway strains were in part due to a shortage in cytosolic Fe–S clusters, which are required cofactors for the dihydroxyacid dehydratase enzyme. We then demonstrated that this cofactor limitation may be partially recovered by disrupting iron homeostasis with a fra2 mutation, thereby increasing cellular iron levels. The resulting isobutanol titer of the fra2 null strain harboring a cytosolic-localized isobutanol pathway outperformed the strain with the mitochondrial-localized pathway by 1.3-fold, demonstrating that both localizations can support flux to isobutanol.

59 BASIC BIOLOGICAL SCIENCES↗

Molecular and Epidemiological Investigation of Fluconazole-resistant Candida parapsilosis —Georgia, United States, 2021

Abstract Background Reports of fluconazole-resistant Candida parapsilosis bloodstream infections are increasing. We describe a cluster of fluconazole-resistant C parapsilosis bloodstream infections identified in 2021 on routine surveillance by the Georgia Emerging Infections Program in conjunction with the Centers for Disease Control and Prevention. Methods Whole-genome sequencing was used to analyze C parapsilosis bloodstream infections isolates. Epidemiological data were obtained from medical records. A social network analysis was conducted using Georgia Hospital Discharge Data. Results Twenty fluconazole-resistant isolates were identified in 2021, representing the largest proportion (34%) of fluconazole-resistant C parapsilosis bloodstream infections identified in Georgia since surveillance began in 2008. All resistant isolates were closely genetically related and contained the Y132F mutation in the ERG11 gene. Patients with fluconazole-resistant isolates were more likely to have resided at long-term acute care hospitals compared with patients with susceptible isolates (P = .01). There was a trend toward increased mechanical ventilation and prior azole use in patients with fluconazole-resistant isolates. Social network analysis revealed that patients with fluconazole-resistant isolates interfaced with a distinct set of healthcare facilities centered around 2 long-term acute care hospitals compared with patients with susceptible isolates. Conclusions Whole-genome sequencing results showing that fluconazole-resistant C parapsilosis isolates from Georgia surveillance demonstrated low genetic diversity compared with susceptible isolates and their association with a facility network centered around 2 long-term acute care hospitals suggests clonal spread of fluconazole-resistant C parapsilosis. Further studies are needed to better understand the sudden emergence and transmission of fluconazole-resistant C parapsilosis.

Misas, Elizabeth (ORCID:0000000162437716)↗

Novel mutation leading to splice donor loss in a conserved site of DMD gene causes Duchenne muscular dystrophy with cryptorchidism

Background As one of the most common congenital abnormalities in male births, cryptorchidism has been found to have a polygenic aetiology according to previous studies of common variants. However, little is known about genetic predisposition of rare variants for cryptorchidism, since rare variants have larger effective size on diseases than common variants. Methods In this study, a cohort of 115 Chinese probands with cryptorchidism was analysed using whole-genome sequencing, alongside 19 parental controls and 2136 unaffected men. Additionally, CRISPR-Cas9 editing of a conserved variant was performed in a mouse model, with MRI screening used to observe the phenotype. Results In 30 of 115 patients (26.1%), we identified four novel genes ( ARSH , DMD , MAGEA4 and SHROOM2 ) affecting at least five unrelated patients and four known genes ( USP9Y , UBA1 , BCORL1 and KDM6A ) with the candidate rare pathogenic variants affecting at least two cases. Burden tests of rare variants revealed the genome-wide significances for newly identified genes (p<2.5×10 −6 ) under the Bonferroni correction. Surprisingly, novel and known genes were mainly found on X chromosome (seven on X and one on Y) and all rare X-chromosomal segregating variants exhibited a maternal inheritance rather than de novo origin. CRISPR-Cas9 mouse modelling of a splice donor loss variant in DMD (NC_000023.11:g.32454661C>G), which resides in a conserved site across vertebrates, replicated bilateral cryptorchidism phenotypes, confirmed by MRI at 4 and 10 weeks. The movement tests further revealed symptoms of Duchenne muscular dystrophy (DMD) in transgenic mice. Conclusion Our results revealed the role of the DMD gene mutation in causing cryptorchidism. The results also suggest that maternal-X inheritance of pathogenic defects could have a predominant role in the development of cryptorchidism.

Chen, Jianhai (ORCID:0000000300932003)↗

CO2 response screen in grass Brachypodium reveals the key role of a MAP kinase in CO2-triggered stomatal closure

Abstract Plants respond to increased CO2 concentrations through stomatal closure, which can contribute to increased water use efficiency. Grasses display faster stomatal responses than eudicots due to dumbbell-shaped guard cells flanked by subsidiary cells working in opposition. However, forward genetic screening for stomatal CO2 signal transduction mutants in grasses has yet to be reported. The grass model Brachypodium distachyon is closely related to agronomically important cereal crops, sharing largely collinear genomes. To gain insights into CO2 control mechanisms of stomatal movements in grasses, we developed an unbiased forward genetic screen with an EMS-mutagenized B. distachyon M5 generation population using infrared imaging to identify plants with altered leaf temperatures at elevated CO2. Among isolated mutants, a “chill1” mutant exhibited cooler leaf temperatures than wild-type Bd21-3 parent control plants after exposure to increased CO2. chill1 plants showed strongly impaired high CO2-induced stomatal closure despite retaining a robust abscisic acid-induced stomatal closing response. Through bulked segregant whole-genome sequencing analyses followed by analyses of further backcrossed F4 generation plants and generation and characterization of sodium azide and CRISPR-cas9 mutants, chill1 was mapped to a protein kinase, Mitogen-Activated Protein Kinase 5 (BdMPK5). The chill1 mutation impaired BdMPK5 protein-mediated CO2/HCO3− sensing together with the High Temperature 1 (HT1) Raf-like kinase in vitro. Furthermore, AlphaFold2-directed structural modeling predicted that the identified BdMPK5-D90N chill1 mutant residue is located at the interface of BdMPK5 with the BdHT1 Raf-like kinase. BdMPK5 is a key signaling component that mediates CO2-induced stomatal movements and is proposed to function as a component of the primary CO2 sensor in grasses.

Lopez, Bryn N. K. (ORCID:0009000937288216)↗