Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SnP”

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 OSTI2020

Genome-wide association analysis of the strength of the MAMP-elicited defense response and resistance to target leaf spot in sorghum

Abstract Plants have the capacity to respond to conserved molecular features known as microbe-associated molecular patterns (MAMPs). The goal of this work was to assess variation in the MAMP response in sorghum, to map loci associated with this variation, and to investigate possible connections with variation in quantitative disease resistance. Using an assay that measures the production of reactive oxygen species, we assessed variation in the MAMP response in a sorghum association mapping population known as the sorghum conversion population (SCP). We identified consistent variation for the response to chitin and flg22—an epitope of flagellin. We identified two SNP loci associated with variation in the flg22 response and one with the chitin response. We also assessed resistance to Target Leaf Spot (TLS) disease caused by the necrotrophic fungus Bipolaris cookei in the SCP. We identified one strong association on chromosome 5 near a previously characterized disease resistance gene. A moderately significant correlation was observed between stronger flg22 response and lower TLS resistance. Possible reasons for this are discussed.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Visualizing novel connections and genetic similarities across diseases using a network-medicine based approach

Understanding the genetic relationships between human disorders could lead to better treatment and prevention strategies, especially for individuals with multiple comorbidities. A common resource for studying genetic-disease relationships is the GWAS Catalog, a large and well curated repository of SNP-trait associations from various studies and populations. Some of these populations are contained within mega-biobanks such as the Million Veteran Program (MVP), which has enabled the genetic classification of several diseases in a large well-characterized and heterogeneous population. Here we aim to provide a network of the genetic relationships among diseases and to demonstrate the utility of quantifying the extent to which a given resource such as MVP has contributed to the discovery of such relations. We use a network-based approach to evaluate shared variants among thousands of traits in the GWAS Catalog repository. Our results indicate many more novel disease relationships that did not exist in early studies and demonstrate that the network can reveal clusters of diseases mechanistically related. Finally, we show novel disease connections that emerge when MVP data is included, highlighting methodology that can be used to indicate the contributions of a given biobank.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Pressurized formic acid dehydrogenation: an entropic spring replaces hydrogen compression cost

Formic acid is unique among liquid organic hydrogen carriers (LOHCs), because its dehydrogenation is highly entropically driven. This enables the evolution of high-pressure hydrogen at mild temperatures that is difficult to achieve with other LOHCs, conceptually by releasing the “spring” of energy stored entropically in the liquid carrier. Applications calling for hydrogen-on-demand, such as vehicle filling, require pressurized H 2 . Hydrogen compression dominates the cost for such applications, yet there are very few reports of selective, catalytic dehydrogenation of formic acid at elevated pressure. Herein, we show that homogenous catalysts with various ligand frameworks, including Noyori-type tridentate (PNP, SNS, SNP, SNPO), bidentate chelates (pyridyl)NHC, (pyridyl)phosphine, (pyridyl)sulfonamide, and their metallic precursors, are suitable catalysts for the dehydrogenation of neat formic acid under self-pressurizing conditions. Quite surprisingly, we discovered that their structural differences can be related to performance differences in their respective structural families, with some tolerant or intolerant of pressure and others that are significantly advantaged by pressurized conditions. Here, we further find important roles for H 2 and CO in catalyst activation and speciation. In fact, for certain systems, CO behaves as a healing reagent when trapped in a pressurizing reactor system, enabling extended life from systems that would be otherwise deactivated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Deep potential generation scheme and simulation protocol for the Li 10 GeP 2 S 12 -type superionic conductors

We report solid-state electrolyte materials with superior lithium ionic conductivities are vital to the next-generation Li-ion batteries. Molecular dynamics could provide atomic scale information to understand the diffusion process of Li-ion in these superionic conductor materials. Here, we implement the deep potential generator to set up an efficient protocol to automatically generate interatomic potentials for Li 10 GeP 2 S 12 -type solid-state electrolyte materials (Li 10 GeP 2 S 12 , Li 10 SiP 2 S 12 , and Li 10 SnP 2 S 12 ). The reliability and accuracy of the fast interatomic potentials are validated. With the potentials, we extend the simulation of the diffusion process to a wide temperature range (300 K–1000 K) and systems with large size (~1000 atoms). Important technical aspects such as the statistical error and size effect are carefully investigated, and benchmark tests including the effect of density functional, thermal expansion, and configurational disorder are performed. The computed data that consider these factors agree well with the experimental results, and we find that the three structures show different behaviors with respect to configurational disorder. Our work paves the way for further research on computation screening of solid-state electrolyte materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

NLR diversity and candidate fusiform rust resistance genes in loblolly pine

Resistance to fusiform rust disease in loblolly pine (Pinus taeda) is a classic gene-for-gene system. Early resistance gene mapping in the P. taeda family 10-5 identified RAPD markers for a major fusiform rust resistance gene, Fr1. More recently, single nucleotide polymorphism (SNP) markers associated with resistance were mapped to a full-length gene model in the loblolly pine genome encoding for a nucleotide-binding site leucine-rich repeat (NLR) protein. NLR genes are one of the most abundant gene families in plant genomes and are involved in effector-triggered immunity. Inter- and intraspecies studies of NLR gene diversity and expression have resulted in improved disease resistance. To characterize NLR gene diversity and discover potential resistance genes, we assembled de novo transcriptomes from 92 loblolly genotypes from across the natural range of the species. In these transcriptomes, we identified novel NLR transcripts that are not present in the loblolly pine reference genome and found significant geographic diversity of NLR genes providing evidence of gene family evolution. We designed capture probes for these NLRs to identify and map SNPs that stably cosegregate with resistance to the SC20-21 isolate of Cronartium quercuum f.sp. fusiforme (Cqf) in half-sib progeny of the 10-5 family. We identified 10 SNPs and 2 quantitative trait loci associated with resistance to SC20-21 Cqf. The geographic diversity of NLR genes provides evidence of NLR gene family evolution in loblolly pine. The SNPs associated with rust resistance provide a resource to enhance breeding and deployment of resistant pine seedlings.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

MINE: maximally informative next experiment—toward a new GWAS experimental design and methodology

Abstract The computational methodology of Genome Wide Association Studies (GWAS) currently has several limitations: (i) the number of observations (rows) on a quantitative trait tends to be smaller than the number of single nucleotide polymorphisms (SNPs) (columns) in the design matrix; (ii) each SNP is usually modeled separately, failing to acknowledge interaction between each other (ie epistasis); (iii) there is implicit linkage disequilibrium (LD) between neighboring SNPs due to their linkage. To overcome these issues, we developed a tool that uses ensemble methods to fit mixed linear models to GWAS data, and these ensemble methods include the development of a new experimental design approach in GWAS, which uses the resultant models and data to select the next informative experiment over time. This new adaptive and staged approach for GWAS experimental design was developed and tested in a 3 yr adaptive model-guided discovery experiment against a fixed classical design. In Sorghum bicolor a total of 79, 86, and 78 accessions were tested in years 1, 2, and 3, respectively out of 343 accessions available in the Bioenergy Association Panel (BAP) each identified for 232,303 SNPs, 1 every 2–3 kb in the genomes. We demonstrated the feasibility of MINE enacted with 8 people in the field per year over 3 yr vs in 1 large classical design enacted with 20 people in 1 yr. The MINE results for chromosomal regions identified controlling dry weight were confirmed against results from previous sorghum GWAS experiments and 1 large classical design for the BAP panel.

Genetics & Heredity↗
DOE OSTI2020

The success of particular Acinetobacter baumannii clones: accumulating resistance and virulence inside a sugary shield

Abstract Background In Portugal, carbapenem-resistant Acinetobacter baumannii (CRAB) has been associated with ST98, ST103 and ST208 (Oxford Scheme, Oxf) and a clone has usually been associated with a particular period of time. These clonal shifts were primarily explained by an increased antimicrobial resistance profile. Here we explore genomic and biochemical differences among these and more recent clones, which could further explain the diversity and evolution of this species. Methods A total of 116 CRAB isolates (2010–15), together with representatives of a previously described CRAB collection (4 isolates, 2001–06) were characterized by attenuated total reflection Fourier transform infrared spectroscopy (FTIR-ATR) and MLST. Representatives of different FTIR-ATR/MLST clusters were selected for WGS (n = 13), which allowed the in silico extraction of resistance and virulence genes, capsule locus and SNP analysis. Results A. baumannii clonal shifts of OXA-58-producing ST103Oxf (2001–04), OXA-40-producing ST98Oxf (2002–06), OXA-23-producing ST208Oxf (2006–10) and OXA-23-producing ST218Oxf (2010–15) were accompanied by an increase in AMR genes and virulence factors. FTIR-ATR clustering was congruent with sugar composition predicted from the capsular locus: a fucosamine cluster comprising ST98Oxf, ST103Oxf and a single ST218Oxf isolate; a pseudaminic acid cluster of ST208Oxf and ST1557Oxf isolates; and legionaminic acid, resembling the sialic acid from mammalian cells, in a cluster comprising ST218Oxf isolates. The whole-genome phylogenetic tree was congruent with MLST, with isolates presenting 5–28 938 SNPs. ST208Oxf and ST218Oxf presented ∼1900 SNPs while ST103Oxf and ST1557Oxf showed a greater number of SNPs (∼28 000). Conclusions Clonal shifts of CRAB were promoted, in our country, by consecutive virulence and AMR gene pool enlargement, together with features increasing pathogen–host adaptation. Worldwide dominance of ST218Oxf is supported by the combination of high AMR and virulence levels.

Silva, Liliana↗
DOE OSTI2021

Fitness cost of vancomycin-resistant Enterococcus faecium plasmids associated with hospital infection outbreaks

Abstract Background Vancomycin resistance is mostly associated with Enterococcus faecium due to Tn1546-vanA located on narrow- and broad-host plasmids of various families. This study’s aim was to analyse the effects of acquiring Tn1546-carrying plasmids with proven epidemicity in different bacterial host backgrounds. Methods Widespread Tn1546-carrying plasmids of different families RepA_N (n = 5), Inc18 (n = 4) and/or pHTβ (n = 1), and prototype plasmids RepA_N (pRUM) and Inc18 (pRE25, pIP501) were analysed. Plasmid transferability and fitness cost were assessed using E. faecium (GE1, 64/3) and Enterococcus faecalis (JH2-2/FA202/UV202) recipient strains. Growth curves (Bioscreen C) and Relative Growth Rates were obtained in the presence/absence of vancomycin. Plasmid stability was analysed (300 generations). WGS (Illumina-MiSeq) of non-evolved and evolved strains (GE1/64/3 transconjugants, n = 49) was performed. SNP calling (Breseq software) of non-evolved strains was used for comparison. Results All plasmids were successfully transferred to different E. faecium clonal backgrounds. Most Tn1546-carrying plasmids and Inc18 and RepA_N prototypes reduced host fitness (–2% to 18%) while the cost of Tn1546 expression varied according to the Tn1546-variant and the recipient strain (9%–49%). Stability of Tn1546-carrying plasmids was documented in all cases, often with loss of phenotypic resistance and/or partial plasmid deletions. SNPs and/or indels associated with essential bacterial functions were observed on the chromosome of evolved strains, some of them linked to increased fitness. Conclusions The stability of E. faecium Tn1546-carrying plasmids in the absence of selective pressure and the high intra-species conjugation rates might explain the persistence of vancomycin resistance in E. faecium populations despite the significant burden they might impose on bacterial host strains.

Tedim, Ana P.↗
DOE OSTI2023

Genomic description of human clinical Aspergillus fumigatus isolates, California, 2020

Abstract Aspergillus fumigatus, an environmental mold, causes life-threatening infections. Studies on the phylogenetic structure of human clinical A. fumigatus isolates are limited. Here, we performed whole genome sequencing of 24 A. fumigatus isolates collected from 18 patients in U.S. healthcare facilities in California. Single-nucleotide polymorphism (SNP) differences between patient isolates ranged from 187 to 70 829 SNPs. For five patients with multiple isolates, we calculated the within-host diversities. Three patients had a within-host diversity that ranged from 4 to 10 SNPs and two patients ranged from 2 to 16 977 SNPs. Findings revealed highly diverse A. fumigatus strains among patients and two patterns of diversity for isolates that come from the same patient, low and extremely high diversity.

Misas, Elizabeth (ORCID:0000000162437716)↗
DOE OSTI2025

Signatures of Selection for Resistance/Tolerance to Perkinsus olseni in Grooved Carpet Shell Clam ( Ruditapes decussatus ) Using a Population Genomics Approach

ABSTRACT The grooved carpet shell clam ( Ruditapes decussatus ) is a bivalve of high commercial value distributed throughout the European coast. Its production has suffered a decline caused by different factors, especially by the parasite Perkinsus olsenii . Improving production of R . decussatus requires genomic resources to ascertain the genetic factors underlying resistance/tolerance to P. olseni i . In this study, the first reference genome of R . decussatus was assembled through long‐ and short‐read sequencing (1677 contigs; 1.386 Mb) and further scaffolded at chromosome level with Hi‐C (19 superscaffolds; 95.4% of assembly). Repetitive elements were identified (32%) and masked for annotation of 38,276 coding‐ and 13,056 non‐coding genes. This genome was used as a reference to develop a 2bRAD‐Seq 13,438 SNP panel for a genomic screening on six shellfish beds distributed across the Atlantic Ocean and Mediterranean Sea. Beds were selected by perkinsosis prevalence and the infection level was individually evaluated in all the samples. Genetic diversity was significantly higher in the Mediterranean than in the Atlantic region. The main genetic breakage was detected between those regions (F ST = 0.224), being the Mediterranean more heterogeneous than the Atlantic. Several loci under divergent selection (394 outliers; 261 genomic windows) were detected across shellfish beds. Samples were also inspected to detect signals of selection for resistance/tolerance to P. olseni i by using infection‐level and population‐genomics approaches, and 90 common divergent outliers for resistance/tolerance to perkinsosis were identified and used for gene mining. Candidate genes and markers identified provide invaluable information for controlling perkinsosis and for improving production of the grooved carpet shell clam.

Sambade, Inés M. [Department of Zoology, Genetics ↗
DOE OSTI2025

Genetic Identification of Lamprey Genera and Anadromous Ecotypes in Watersheds of the Northeastern Pacific Ocean

ABSTRACT Nonparasitic, nonmigratory Western Brook Lamprey (WBL; Lampetra ayresii ), and parasitic, anadromous Western River Lamprey (WRL; L. ayresii ) are sympatric lampreys that likely represent different life history variations of a single species. Novel genetic tools are critical for differentiating WBL and WRL, whose larvae preclude morphological identification (ID) and will enable comprehensive assessment of imperiled native lampreys of the Northeastern Pacific (including WBL, WRL, and Pacific Lamprey, Entosphenus tridentatus ). We developed 47 candidate single nucleotide polymorphism (SNP) markers using whole genome resequencing of WBL ( N = 24) and WRL ( N = 15) from Ksi Ts'oohl Ts'ap Creek (Nass River, British Columbia, Canada) which are likely ecotypes distinguished by few divergent SNPs across multiple chromosomes. We used five novel candidate SNPs to perform genetic ID of WBL and WRL ecotypes in collections of mixed native lampreys from lower Columbia River tributaries ( N = 1474), Ksi Ts'oohl Ts'ap Creek ( N = 352), and ocean phase WRL from the Georgia Basin (Salish Sea, British Columbia, Canada; N = 91). Two previously published SNPs were used to ID genera, Entosphenus versus Lampetra . Morphological ID utilized photographs collected from a subset of genotyped lampreys, and high concordance was demonstrated between ID methods for genera (99%) and Lampetra ecotypes (> 98%). We characterized spatial and temporal composition of lamprey genera and ecotypes surveyed across NE Pacific tributaries under the expectation these compositions would be similar across nearby sites and across years at the same site. Proportions of lamprey genera were highly variable within regions and across years; however, Lampetra ecotypic proportions were spatially and temporally stable. WRL were rare in lower Columbia tributaries (~1% average rate among Lampetra ) and common further north (> 40% of Lampetra ). Genetic ID methods are powerful monitoring tools that create the novel ability to ascertain genera and ecotypes regardless of life stage, while increasing the efficiency of surveys by eliminating time‐intensive morphological data collection.

Silver, G. S. [Columbia River Inter‐Tribal Fish Co↗
DOE OSTI2023

Genomic variation across Chinook salmon populations reveals effects of a duplication on migration alleles and supports fine scale structure

Abstract The distribution of ecotypic variation in natural populations is influenced by neutral and adaptive evolutionary forces that are challenging to disentangle. This study provides a high‐resolution portrait of genomic variation in Chinook salmon ( Oncorhynchus tshawytscha ) with emphasis on a region of major effect for ecotypic variation in migration timing. With a filtered data set of ~13 million single nucleotide polymorphisms (SNPs) from low‐coverage whole genome resequencing of 53 populations (3566 barcoded individuals), we contrasted patterns of genomic structure within and among major lineages and examined the extent of a selective sweep at a major effect region underlying migration timing (GREB1L/ROCK1). Neutral variation provided support for fine‐scale structure of populations, while allele frequency variation in GREB1L/ROCK1 was highly correlated with mean return timing for early and late migrating populations within each of the lineages ( r 2 = .58–.95; p < .001). However, the extent of selection within the genomic region controlling migration timing was much narrower in one lineage (interior stream‐type) compared to the other two major lineages, which corresponded to the breadth of phenotypic variation in migration timing observed among lineages. Evidence of a duplicated block within GREB1L/ROCK1 may be responsible for reduced recombination in this portion of the genome and contributes to phenotypic variation within and across lineages. Lastly, SNP positions across GREB1L/ROCK1 were assessed for their utility in discriminating migration timing among lineages, and we recommend multiple markers nearest the duplication to provide highest accuracy in conservation applications such as those that aim to protect early migrating Chinook salmon. These results highlight the need to investigate variation throughout the genome and the effects of structural variants on ecologically relevant phenotypic variation in natural species.

Horn, Rebekah L.↗
DOE OSTI2024

The auxin efflux carrier PIN1a regulates vascular patterning in cereal roots

Barley (Hordeum vulgare) is an important global cereal crop and a model in genetic studies. Despite advances in characterising barley genomic resources, few mutant studies have identified genes controlling root architecture and anatomy, which plays a critical role in capturing soil resources. Our phenotypic screening of a TILLING mutant collection identified line TM5992 exhibiting a short-root phenotype compared with wild-type (WT) Morex background. Outcrossing TM5992 with barley variety Proctor and subsequent SNP array-based bulk segregant analysis, fine mapped the mutation to a cM scale. Exome sequencing pinpointed a mutation in the candidate gene HvPIN1a, further confirming this by analysing independent mutant alleles. Detailed analysis of root growth and anatomy in Hvpin1a mutant alleles exhibited a slower growth rate, shorter apical meristem and striking vascular patterning defects compared to WT. Expression and mutant analyses of PIN1 members in the closely related cereal brachypodium (Brachypodium distachyon) revealed that BdPIN1a and BdPIN1b were redundantly expressed in root vascular tissues but only Bdpin1a mutant allele displayed root vascular defects similar to Hvpin1a. We conclude that barley PIN1 genes have sub-functionalised in cereals, compared to Arabidopsis (Arabidopsis thaliana), where PIN1a sequences control root vascular patterning.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Soybean ( Glycine max ) Haplotype Map (GmHapMap): a universal resource for soybean translational and functional genomics

Here, we describe a worldwide haplotype map for soybean (GmHapMap) constructed using whole-genome sequence data for 1007 Glycine max accessions and yielding 14.9 million variants as well as 4.3 M tag single-nucleotide polymorphisms (SNPs). When sampling random subsets of these accessions, the number of variants and tag SNPs plateaued beyond approximately 800 and 600 accessions, respectively. This suggests extensive coverage of diversity within the cultivated soybean. GmHapMap variants were imputed onto 21 618 previously genotyped accessions with up to 96% success for common alleles. A local association analysis was performed with the imputed data using markers located in a 1-Mb region known to contribute to seed oil content and enabled us to identify a candidate causal SNP residing in the NPC1 gene. We determined gene-centric haplotypes (407 867 GCHs) for the 55 589 genes and showed that such haplotypes can help to identify alleles that differ in the resulting phenotype. Finally, we predicted 18 031 putative loss-of-function (LOF) mutations in 10 662 genes and illustrated how such a resource can be used to explore gene function. The GmHapMap provides a unique worldwide resource for applied soybean genomics and breeding.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

The first two chromosome‐scale genome assemblies of American hazelnut enable comparative genomic analysis of the genus Corylus

Summary The native, perennial shrub American hazelnut ( Corylus americana ) is cultivated in the Midwestern United States for its significant ecological benefits, as well as its high‐value nut crop. Implementation of modern breeding methods and quantitative genetic analyses of C. americana requires high‐quality reference genomes, a resource that is currently lacking. We therefore developed the first chromosome‐scale assemblies for this species using the accessions ‘Rush’ and ‘Winkler’. Genomes were assembled using HiFi PacBio reads and Arima Hi‐C data, and Oxford Nanopore reads and a high‐density genetic map were used to perform error correction. N50 scores are 31.9 Mb and 35.3 Mb, with 90.2% and 97.1% of the total genome assembled into the 11 pseudomolecules, for ‘Rush’ and ‘Winkler’, respectively. Gene prediction was performed using custom RNAseq libraries and protein homology data. ‘Rush’ has a BUSCO score of 99.0 for its assembly and 99.0 for its annotation, while ‘Winkler’ had corresponding scores of 96.9 and 96.5, indicating high‐quality assemblies. These two independent assemblies enable unbiased assessment of structural variation within C. americana , as well as patterns of syntenic relationships across the Corylus genus. Furthermore, we identified high‐density SNP marker sets from genotyping‐by‐sequencing data using 1343 C. americana , C. avellana and C. americana × C. avellana hybrids, in order to assess population structure in natural and breeding populations. Finally, the transcriptomes of these assemblies, as well as several other recently published Corylus genomes, were utilized to perform phylogenetic analysis of sporophytic self‐incompatibility (SSI) in hazelnut, providing evidence of unique molecular pathways governing self‐incompatibility in Corylus .

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Analysis of the ARTIC Version 3 and Version 4 SARS-CoV-2 Primers and Their Impact on the Detection of the G142D Amino Acid Substitution in the Spike Protein

The ARTIC Network provides a common resource of PCR primer sequences and recommendations for amplifying SARS-CoV-2 genomes. The initial tiling strategy was developed with the reference genome Wuhan-01, and subsequent iterations have addressed areas of low amplification and sequence drop out. Recently, a new version (V4) was released, based on new variant genome sequences, in response to the realization that some V3 primers were located in regions with key mutations. Herein, we compare the performance of the ARTIC V3 and V4 primer sets with a matched set of 663 SARS-CoV-2 clinical samples sequenced with an Illumina NovaSeq 6000 instrument. We observe general improvements in sequencing depth and quality, and improved resolution of the SNP causing the D950N variation in the spike protein. Importantly, we also find nearly universal presence of spike protein substitution G142D in Delta-lineage samples. Due to the prior release and widespread use of the ARTIC V3 primers during the initial surge of the Delta variant, it is likely that the G142D amino acid substitution is substantially underrepresented among early Delta variant genomes deposited in public repositories. In addition to the improved performance of the ARTIC V4 primer set, this study also illustrates the importance of the primer scheme in downstream analyses.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

The landscape of regulatory element evolution in a C4 perennial grass

Gene regulatory evolution is a well-known source of phenotypic diversity and adaptive evolution. Although cis-regulatory elements (CREs) play a vital role in gene expression evolution, the molecular evolution of CREs remains mostly unknown due to the difficulty in identifying and characterizing these functional elements. Comparative genomic analyses of noncoding DNA can be leveraged to identify conserved noncoding sequences (CNS), many of which may harbor functional CREs conserved by purifying selection. However, purely computational inference of CREs from putative CNS can be erroneous due to the complex genomic architecture in plants. One promising experimental approach to identify CREs is by profiling accessible chromatin regions (ACRs) that are often associated with the location of CREs. In this study, we use comparative genomics along with the profiling of ACRs to study the molecular evolution of putative functional noncoding regulatory regions in Panicoid grasses. We identified sets of CNS that varied in relationship to the degree of evolutionary divergence among the studied taxa, including identifying core-Panicoid-CNS. We augmented this analysis by profiling ACRs in Panicum hallii ecotypes using ATAC-seq. ACRs had low SNP density at the summit, harbored a high frequency of core-Panicoid-CNS, and were enriched with expression QTL. These data help to annotate the P. hallii genome for putative functional elements and suggest that a large proportion of these ACRs are evolving under purifying selection. Turnover in CNS and ACR between ecotypes of P. hallii identifies a small set of putatively divergent CREs that may underlie differences in gene regulation between genotypes from inland and coastal habitats. In summary, we profiled ACRs in Panicoid grasses and integrated this data with our putative CNS prediction framework, which provides unique insight into patterns of polymorphism and divergence in CREs in C4 perennial grasses.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Plastome evolution in annual Brachypodium species reveals widespread heteroplasmy and chloroplast capture, lineage-specific codon usage bias, and low positive selection

Comparative genomics and plastome phylogenomics have advanced significantly in recent years, highlighting the diversity, possible admixture, and non-neutral evolution of the predominantly considered non-recombinant chloroplast genomes in angiosperms. The grass genus Brachypodium serves as a powerful model for studying evolutionary processes in monocots. We analyzed 287 plastomes across the native circum-Mediterranean range of the three annual Brachypodium species ( B. distachyon, B. stacei, B. hybridum ), focusing on their structural variation, selection patterns and phylogenomic relationships. Our analyses confirmed the differentiation of the S and D plastomes, inherited respectively from the diploid progenitor species B. stacei and B. distachyon . We identified novel structural rearrangements and indels, and unique repeat motifs, along with widespread heteroplasmy, particularly in ancestral B. hybridum -D plastotypes. SNP diversity varied among plastotypes, reflecting population dynamics and evolutionary histories, with B. hybridum -D plastotypes showing the highest normalized diversity and B. hybridum -S the lowest. Positive selection was detected in 29 plastid genes by Tajima’s neutrality test, and in nine genes by site and branch-site evolutionary models, including matK, ndhF, rbcL, and rpoC2. Phylogenomic analyses revealed well-supported clades corresponding to the S and D plastome lineages, with frequent chloroplast capture events and long-distance dispersals shaping their evolutionary trajectories.

allopolyploidy↗