Search NASA⌕ Search

SEARCH · Search NASA

Results for “Genetic prediction”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

The unique architecture of umbrella toxins permits a two-tiered molecular bet-hedging strategy for interbacterial antagonism

Bacteria exist in competitive and rapidly changing environments in which the nature of future threats cannot be easily predicted. Streptomyces coelicolor produces three antibacterial umbrella particles that harbor distinct polymorphic toxin domains and an overlapping set of six diversified lectins. Here, we show that the exquisite specificity of umbrella particles derives from lectin-mediated species-specific binding to previously undescribed hypervariable surface glycoconjugates. A cryo-electron microscopy (cryo-EM) structure of one such lectin in complex with its oligosaccharide substrate defines the molecular basis for targeting through the coordinated recognition of multiple glycan features. Biochemical and genetic studies of several target species, in conjunction with lectin-swapping experiments, support a model whereby S. coelicolor umbrella toxin diversification at the levels of lectin composition and toxin polymorphism represents a unique, two-tiered bet-hedging strategy. Bioinformatic analyses support this as a means by which the unusual architecture of umbrella toxins offers Streptomyces a generalizable strategy to antagonize an unpredictable array of competitors.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

High–throughput measurement of plant fitness traits with an object detection method using Faster R–CNN

Revealing the contributions of genes to plant phenotype is frequently challenging because loss-of-function effects may be subtle or masked by varying degrees of genetic redundancy. Such effects can potentially be detected by measuring plant fitness, which reflects the cumulative effects of genetic changes over the lifetime of a plant. However, fitness is challenging to measure accurately, particularly in species with high fecundity and relatively small propagule sizes such as Arabidopsis thaliana. An image segmentation-based method using the software ImageJ and an object detection-based method using the Faster Region-based Convolutional Neural Network (R-CNN) algorithm were used for measuring two Arabidopsis fitness traits: seed and fruit counts. The segmentation-based method was error-prone (correlation between true and predicted seed counts, r 2 = 0.849) because seeds touching each other were undercounted. By contrast, the object detection-based algorithm yielded near perfect seed counts (r 2 = 0.9996) and highly accurate fruit counts (r 2 = 0.980). Comparing seed counts for wild-type and 12 mutant lines revealed fitness effects for three genes; fruit counts revealed the same effects for two genes. Our study provides analysis pipelines and models to facilitate the investigation of Arabidopsis fitness traits and demonstrates the importance of examining fitness traits when studying gene functions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Hubbard-corrected oxide formation enthalpies without adjustable parameters

A density functional theory (DFT) approach to computing transition metal oxide heat of formation without adjustable parameters is presented. Different degrees of d-electron localization in oxides are treated within the DFT+U approach with site-dependent, first-principles Hubbard U-parameters obtained from linear response theory, and delocalized states in the metallic phases are treated without Hubbard corrections. Comparison of relative stabilities of these differently treated phases is enabled by a local d-electron density matrix-dependent model, which was found by genetic programming against experimental reference formation enthalpies. This mathematically simple model does not explicitly depend on the Hubbard-corrected ionic species and is shown to reproduce the heats of formation of the Mott insulators Ca 2 RuO 4 and Y 2 Ru 2 O 7 within ~3% of experimental results, where the experimental training data did not contain Ru oxides. This newly developed method thus absolves from the need for element-specific corrections fitted to experiments in existing Hubbard-corrected approaches to the prediction of reaction energies of transition metal oxides and metals. The absence of fitting parameters opens up here the possibility to predict relative thermodynamic stabilities and reaction energies involving d-states of varying degree of localization at transition metal oxide interfaces and defects, where site-dependent U-parameters will be particularly important and devising a fitting scheme against experimental data with predictive power would be exceedingly difficult.

transition metal oxides↗
DOE OSTI2023

From Optimization to Sampling Through Gradient Flows

Optimization and sampling algorithms play a central role in science and engineering as they enable finding optimal predictions, policies, and recommendations, as well as expected and equilibrium states of complex systems. The notion of “optimality” is formalized by the choice of an objective function, while the notion of an “expected” state is specified by a probabilistic model for the distribution of states. Optimizing rugged objective functions and sampling multimodal distributions is computationally challenging, especially in high-dimensional problems. Here, for this reason, many optimization and sampling methods have been developed by researchers working in disparate fields such as Bayesian statistics, molecular dynamics, genetics, quantum chemistry, machine learning, weather forecasting, econometrics, and medical imaging.

Trillos, N. García↗
DOE OSTI2023

RNAseq analysis of Cellvibrio japonicus during starch utilization differentiates between genes encoding carbohydrate active enzymes controlled by substrate detection or growth rate

ABSTRACT Bacterial utilization of starch is increasingly of interest as the importance and contributions of animal gut microbiomes become more defined. Consequently, identifying and characterizing the bacterial enzymes responsible for the degradation, transport, and metabolism of starch will enable developments in pharmaceutical, biotechnological, and culinary industries searching for novel prebiotics, carrier molecules, and low glycemic index sweeteners. The current challenge is that bacteria proficient at starch utilization often have hundreds of carbohydrate active enzymes, and it is unclear which are essential for starch utilization using only homology-based bioinformatics or computational methods. Complementary experimental data are also needed, especially to understand the regulation of bacterial starch utilization. We have completed an RNAseq analysis of the Gram-negative bacterium Cellvibrio japonicus and found that it has sophisticated regulation that includes substrate sensing and growth rate components for genes that encode starch-degrading enzymes. Among the 22 genes predicted to encode starch-active enzymes, C. japonicus has 10 alpha-amylases, 4 alpha-glucosidases, 2 pullulnases, and 2 cyclomaltodextrin glucanotransferases, 15 of which were up-regulated during exponential growth on starch and 8 up-regulated in stationary phase. Growth analyses with an enzyme secretion deficient mutant of C. japonicus suggested that secreted amylases are essential for this bacterium to degrade starch. Our approach of coupling a physiological growth assay with transcriptomic data provides a platform to identify targets for further genetic or biochemical analysis that can be broadly applied to other starch-utilizing bacteria. IMPORTANCE Understanding the bacterial metabolism of starch is important as this polysaccharide is a ubiquitous ingredient in foods, supplements, and medicines, all of which influence gut microbiome composition and health. Our RNAseq and growth data set provides a valuable resource to those who want to better understand the regulation of starch utilization in Gram-negative bacteria. These data are also useful as they provide an example of how to approach studying a starch-utilizing bacterium that has many putative amylases by coupling transcriptomic data with growth assays to overcome the potential challenges of functional redundancy. The RNAseq data can also be used as a part of larger meta-analyses to compare how C. japonicus regulates carbohydrate active enzymes, or how this bacterium compares to gut microbiome constituents in terms of starch utilization potential.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Enhanced resonant ultrasound spectroscopy for measurement of the elastic properties of multi-material systems

Understanding the elastic properties of materials is critical for their safe incorporation and predictable performance. Current methods of bulk elastic characterization often have notable limitations for in situ structural applications, with usage restricted to simple geometries and material distributions. To address these existing issues, this study sought to expand the capabilities of resonant ultrasound spectroscopy (RUS), an established nondestructive evaluation method, to include the characterization of isotropic multi-material samples. In this work, finite-element-based RUS analysis consisted of numerical simulations and experimental testing of composite samples comprised of material pairs with varying elasticity and density contrasts. Utilizing genetic algorithm inversion and mode matching, our results demonstrate that elastic properties of multi-material samples can be reliably identified within several percent of known or nominal values using a minimum number of identified resonance modes, given sample mass is held consistent. The accurate recovery of material properties for composite samples of varying material similarity and geometry expands the pool of viable samples for RUS and advances the method towards in situ inspection and evaluation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Li + Diffusion in Amorphous and Crystalline Al 2 O 3 for Battery Electrode Coatings

Al 2 O 3 is often applied protectively to lithium-ion battery anode and cathode materials to inhibit surface degradation, suppress dendrite formation, and relieve mechanical stresses. Given the very high intrinsic band gap and diffusion barrier of the material, the mechanism that allows Li diffusion through these coatings is not well understood, and widely varying laboratory results indicate that there may be dependencies on morphology and stoichiometry. Using nudged elastic band calculations and ab initio molecular dynamics, we perform a systematic investigation across Al 2 O 3 structures, both crystalline and amorphous, and at various concentrations of Li + to uncover the optimal parameters for maximally diffusive coatings. We find a correlation between the low proximity of Li + to Al 3+ and the low Li + migration barrier. Although barriers are the lowest in the highly diffusive one-dimensional channels of crystalline θ-Al 2 O 3 , the system is structurally delicate and subject to detrimental distortion as the Li + content is increased. The α-Al 2 O 3 lattice is, conversely, highly stable against distortion at all Li + concentrations but disadvantageous for Li + migration. In amorphous systems, unscreened Li + –Li + Coulomb repulsion and pre-emptive occupation of “trapping sites” combine to lower the energy barriers as a function of increasing concentration. One of our most important findings is that Al-deficient materials can sharply increase Li + movement, and we predict that an amorphous material with a combination of high Li + concentration and Al deficiency would enable highly Li + -conductive protective coatings for electrodes.

25 ENERGY STORAGE↗
DOE OSTI2022

One hundred years into the study of ecotypes, new advances are being made through large-scale field experiments in perennial plant systems

A hundred years after Turesson first clearly described how locally adaptive variation is distributed within species, plant biologists are making major breakthroughs in our understanding of mechanisms underlying adaptation from local populations to the scale of continents. Although the genetics of local adaptation has typically been studied in smaller reciprocal transplant experiments, here it is being evaluated with whole genomes in large-scale networks of common garden experiments with perennial switchgrass and poplar trees. These studies support the hypothesis that a complex combination of loci, both with and without adaptive trade-offs, underlies local adaptation and that hybridization and adaptive introgression play a key role in the evolution of these species. Future studies incorporating high-throughput phenotyping, gene expression, and modeling will be used to predict responses of these species to climate change.

Trade-off↗
DOE OSTI2021

Antimicrobial resistance and genomic characterization of Salmonella Dublin isolates in cattle from the United States

Salmonella enterica subspecies enterica serotype Dublin is a host-adapted serotype in cattle, associated with enteritis and systemic disease. The primary clinical manifestation of Salmonella Dublin infection in cattle, especially calves, is respiratory disease. While rare in humans, it can cause severe illness, including bacteremia, with hospitalization and death. In the United States, S . Dublin has become one of the most multidrug-resistant serotypes. The objective of this study was to characterize S . Dublin isolates from sick cattle by analyzing phenotypic and genotypic antimicrobial resistance (AMR) profiles, the presence of plasmids, and phylogenetic relationships. S . Dublin isolates (n = 140) were selected from submissions to the NVSL for Salmonella serotyping (2014–2017) from 21 states. Isolates were tested for susceptibility against 14 class-representative antimicrobial drugs. Resistance profiles were determined using the ABRicate with Resfinder and NCBI databases, AMRFinder and PointFinder. Plasmids were detected using ABRicate with PlasmidFinder. Phylogeny was determined using vSNP. We found 98% of the isolates were resistant to more than 4 antimicrobials. Only 1 isolate was pan-susceptible and had no predicted AMR genes. All S . Dublin isolates were susceptible to azithromycin and meropenem. They showed 96% resistance to sulfonamides, 97% to tetracyclines, 95% to aminoglycosides and 85% to beta-lactams. The most common AMR genes were: sulf2 and tetA (98.6%), aph(6)-Id (97.9%), aph(3’’)-Ib, (97.1%), floR (94.3%), and blaCMY-2 (85.7%). All quinolone resistant isolates presented mutations in gyr A. Ten plasmid types were identified among all isolates with IncA/C2, IncX1, and IncFII(S) being the most frequent. The S . Dublin isolates show low genomic genetic diversity. This study provided antimicrobial susceptibility and genomic insight into S . Dublin clinical isolates from cattle in the U.S. Further sequence analysis integrating food and human origin S . Dublin isolates may provide valuable insight on increased virulence observed in humans.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Bullying among children with heart conditions, National Survey of Children’s Health, 2018–2020

Abstract Children with chronic illnesses report being bullied by peers, yet little is known about bullying among children with heart conditions. Using 2018–2020 National Survey of Children’s Health data, the prevalence and frequency of being bullied in the past year (never; annually or monthly; weekly or daily) were compared between children aged 6–17 years with and without heart conditions. Among children with heart conditions, associations between demographic and health characteristics and being bullied, and prevalence of diagnosed anxiety or depression by bullying status were examined. Differences were assessed with chi-square tests and multivariable logistic regression using predicted marginals to produce adjusted prevalence ratios and 95% confidence intervals. Weights yielded national estimates. Of 69,428 children, 2.2% had heart conditions. Children with heart conditions, compared to those without, were more likely to be bullied (56.3% and 43.3% respectively; adjusted prevalence ratio [95% confidence interval] = 1.3 [1.2, 1.4]) and bullied more frequently (weekly or daily = 11.2% and 5.3%; p < 0.001). Among children with heart conditions, characteristics associated with greater odds of weekly or daily bullying included ages 9–11 years compared to 15–17 years (3.4 [2.0, 5.7]), other genetic or inherited condition (1.7 [1.0, 3.0]), ever overweight (1.7 [1.0, 2.8]), and a functional limitation (4.8 [2.7, 8.5]). Children with heart conditions who were bullied, compared to never, more commonly had anxiety (40.1%, 25.9%, and 12.8%, respectively) and depression (18.0%, 9.3%, and 4.7%; p < 0.01 for both). Findings highlight the social and psychological needs of children with heart conditions.

Cardiovascular System & Cardiology↗
DOE OSTI2023

Study of Tritium Diffusivity in Pure and Sn-Defective Zr: A First-Principles Density Functional Theory Approach

Zirconium alloys (e.g., zircaloy-4) are used as tritium ( 3 H) getter materials in tritium-producing burnable absorber rods (TPBARs) owing to their ability to capture 3 H and chemically convert them into metal hydrides. Understanding of 3 H diffusion mechanisms in zircaloy is crucial for the optimal design of material performance in nuclear technology. Here, for this work, we perform first-principles density functional theory calculations to study the 3 H diffusion mechanism in pure and impure Zr with a low concentration of tin (Sn) atoms to determine the impact of the presence of Sn on the movement of 3 H atoms through the material. First, we calculated the diffusion barriers for 3 H in pure Zr by taking different migration pathways. We then introduced a low concentration of Sn impurity and systematically explored the impurity effect on the diffusion barriers for 3 H. Using our calculated diffusion energy barriers, we further obtained the diffusion coefficients and analyzed the results by comparing them with the experimental and previously calculated values. A diffusion coefficient of the order of 10 –8 m 2 /s is predicted. We also found that the presence of a Sn impurity could reduce the diffusivity up to 4 orders of magnitude. In conclusion, our results could serve as guidelines for further experimental investigations.

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

Systems-Level Modeling for CRISPR-Based Metabolic Engineering

The CRISPR-Cas system has enabled the development of sophisticated, multigene metabolic engineering programs through the use of guide RNA-directed activation or repression of target genes. To optimize biosynthetic pathways in microbial systems, we need improved models to inform design and implementation of transcriptional programs. Recent progress has resulted in new modeling approaches for identifying gene targets and predicting the efficacy of guide RNA targeting. Genome-scale and flux balance models have successfully been applied to identify targets for improving biosynthetic production yields using combinatorial CRISPR-interference (CRISPRi) programs. Here, the advent of new approaches for tunable and dynamic CRISPR activation (CRISPRa) promises to further advance these engineering capabilities. Once appropriate targets are identified, guide RNA prediction models can lead to increased efficacy in gene targeting. Developing improved models and incorporating approaches from machine learning may be able to overcome current limitations and greatly expand the capabilities of CRISPR-Cas9 tools for metabolic engineering.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Genome-Wide Analysis of RNA Decay in the Cyanobacterium Synechococcus sp. Strain PCC 7002

ABSTRACT RNA degradation is an important process that influences the ultimate concentration of individual proteins inside cells. While the main enzymes that facilitate this process have been identified, global maps of RNA turnover are available for only a few species. Even in these cases, there are few sequence elements that are known to enhance or destabilize a native transcript; even fewer confer the same effect when added to a heterologous transcript. To address this knowledge gap, we assayed genome-wide RNA degradation in the cyanobacterium Synechococcus sp. strain PCC 7002 by collecting total RNA samples after stopping nascent transcription with rifampin. We quantified the abundance of each position in the transcriptome as a function of time using RNA-sequencing data and later analyzed the global mRNA decay map using machine learning principles. Half-lives, calculated on a per-ORF (open reading frame) basis, were extremely short, with a median half-life of only 0.97 min. Despite extremely rapid turnover of most mRNA, transcripts encoding proteins involved in photosynthesis were both highly expressed and highly stable. Upon inspection of these stable transcripts, we identified an enriched motif in the 3′ untranslated region (UTR) that had similarity to Rho-independent terminators. We built statistical models for half-life prediction and used them to systematically identify sequence motifs in both 5′ and 3′ UTRs that correlate with stabilized transcripts. We found that transcripts linked to a terminator containing a poly(U) tract had a longer half-life than both those without a poly(U) tract and those without a terminator. IMPORTANCE RNA degradation is an important process that affects the final concentration of individual mRNAs, affecting protein expression and cellular physiology. Studies of how RNA is degraded increase our knowledge of this fundamental process as well as enable the creation of genetic tools to manipulate RNA stability. By studying global transcript turnover, we searched for sequence elements that correlated with transcript (in)stability and used these sequences to guide tool design. This study probes global RNA turnover in a cyanobacterium, Synechococcus sp. strain PCC 7002, that both has a unique array of RNases that facilitate RNA degradation and is an industrially relevant strain that could be used to convert CO 2 and sunlight into useful products.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Predicting Field Efficiency of Round-Baling Operations in High-Yielding Biomass Crops

Model simulations for bioenergy harvest planning need to utilize equipment-capacity relationships for equipment operating under the high-yield conditions typical of a biomass crop. These performance assumptions have a direct bearing on the estimates of machine capacity, the number of machines required, and, therefore, the cost to fulfill the biorefinery plant demands for a given harvest window. Typically, two major issues in these models have been poorly understood: the available time required to complete the harvest operation (often called probability of workdays) and the capacity of the harvest equipment as impacted by yield. Simulations use annual yield estimates, which incorporate weather events, to demonstrate year-to-year effects. Some simulations also incorporate potential yield increases from genetically modified energy crops. There are limited field performance data for most current forage equipment used for harvesting high-yield biomass crops. Analysis shows that the impact of wrap/eject time for round balers resulted in a 50% reduction in achieved throughput capacity (Mg/h). After the maximum throughput is reached, the cost of the round bale operation (3.23 USD/Mg) is double that of the large-square baler (1.63 USD/Mg). The round baler achieved throughput capacity is 50% less (32.7 Mg/h compared to 71.0 Mg/h) than the large-square baler.

09 BIOMASS FUELS↗
DOE OSTI2022

The histone code of the fungal genus Aspergillus uncovered by evolutionary and proteomic analyses

Chemical modifications of DNA and histone proteins impact the organization of chromatin within the nucleus. Changes in these modifications, catalysed by different chromatin-modifying enzymes, influence chromatin organization, which in turn is thought to impact the spatial and temporal regulation of gene expression. While combinations of different histone modifications, the histone code, have been studied in several model species, we know very little about histone modifications in the fungal genus Aspergillus, whose members are generally well studied due to their importance as models in cell and molecular biology as well as their medical and biotechnological relevance. Here, we used phylogenetic analyses in 94 Aspergilli as well as other fungi to uncover the occurrence and evolutionary trajectories of enzymes and protein complexes with roles in chromatin modifications or regulation. We found that these enzymes and complexes are highly conserved in Aspergilli, pointing towards a complex repertoire of chromatin modifications. Nevertheless, we also observed few recent gene duplications or losses, highlighting Aspergillus species to further study the roles of specific chromatin modifications. SET7 (KMT6) and other components of PRC2 (Polycomb Repressive Complex 2), which is responsible for methylation on histone H3 at lysine 27 in many eukaryotes including fungi, are absent in Aspergilli as well as in closely related Penicillium species, suggesting that these lost the capacity for this histone modification. We corroborated our computational predictions by performing untargeted MS analysis of histone post-translational modifications in Aspergillus nidulans. This systematic analysis will pave the way for future research into the complexity of the histone code and its functional implications on genome architecture and gene regulation in fungi.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Modeling the Effects of Trimethylsilanol on Syngas Combustion Kinetics

A perturbation analysis modeling approach based on a genetic algorithm was used to identify possible reaction pathways to explain previous experimental observations of the strong acceleration of syngas auto-ignition by trimethylsilanol (TMSO). Organosilicon reactions were taken from an existing chemical kinetic mechanism for tetramethylsilane pyrolysis which also contained TMSO reactions. The experimental targets for optimization were ignition delay times and peak OH radical concentrations in the range 1010–1070 K, 5 atm with dilute (4%) ϕ = 0.1 syngas mixtures in N 2 /Ar containing 0, 200, and 1000 ppm TMSO. Hundreds of thousands of trials resulted in three models which accurately predicted the experimental auto-ignition times and peak OH concentrations. Rate of production analyses indicated that TMSO forms dimethylsilanediol, which then reacts creating two “catalytic loops.” The net effect of the loops is to accelerate the reaction H 2 + O = OH + H at a faster effective rate than H 2 /O 2 kinetics alone, thus enhancing the reactivity of syngas. The formation of the loops may be attributable to the high Si–C bond energies in siloxanols, reducing rates of thermal decomposition and allowing dimethylsilanediol to react with H 2 for a longer duration compared with its alkane analog. Here, the results of the current work provide valuable direction for fundamental studies of these new hypothesized reaction pathways.

Mansfield, Andrew B. [Eastern Michigan University,↗
DOE OSTI2022

Genetic and Structural Analysis of SARS-CoV-2 Spike Protein for Universal Epitope Selection

Evaluation of immunogenic epitopes for universal vaccine development in the face of ongoing SARS-CoV-2 evolution remains a challenge. Herein, we investigate the genetic and structural conservation of an immunogenically relevant epitope (C662–C671) of spike (S) protein across SARS-CoV-2 variants to determine its potential utility as a broad-spectrum vaccine candidate against coronavirus diseases. Comparative sequence analysis, structural assessment, and molecular dynamics simulations of C662–C671 epitope were performed. Mathematical tools were employed to determine its mutational cost. We found that the amino acid sequence of C662–C671 epitope is entirely conserved across the observed major variants of SARS-CoV-2 in addition to SARS-CoV. Its conformation and accessibility are predicted to be conserved, even in the highly mutated Omicron variant. Costly mutational rate in the context of energy expenditure in genome replication and translation can explain this strict conservation. These observations may herald an approach to developing vaccine candidates for universal protection against emergent variants of coronavirus.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Divergent selection and climate adaptation fuel genomic differentiation between sister species of Sphagnum (peat moss)

Abstract Background and Aims New plant species can evolve through the reinforcement of reproductive isolation via local adaptation along habitat gradients. Peat mosses (Sphagnaceae) are an emerging model system for the study of evolutionary genomics and have well-documented niche differentiation among species. Recent molecular studies have demonstrated that the globally distributed species Sphagnum magellanicum is a complex of morphologically cryptic lineages that are phylogenetically and ecologically distinct. Here, we describe the architecture of genomic differentiation between two sister species in this complex known from eastern North America: the northern S. diabolicum and the largely southern S. magniae. Methods We sampled plant populations from across a latitudinal gradient in eastern North America and performed whole genome and restriction-site associated DNA sequencing. These sequencing data were then analyzed computationally. Key Results Using sliding-window population genetic analyses we find that differentiation is concentrated within ‘islands’ of the genome spanning up to 400 kb that are characterized by elevated genetic divergence, suppressed recombination, reduced nucleotide diversity and increased rates of non-synonymous substitution. Sequence variants that are significantly associated with genetic structure and bioclimatic variables occur within genes that have functional enrichment for biological processes including abiotic stress response, photoperiodism and hormone-mediated signalling. Demographic modelling demonstrates that these two species diverged no more than 225 000 generations ago with secondary contact occurring where their ranges overlap. Conclusions We suggest that this heterogeneity of genomic differentiation is a result of linked selection and reflects the role of local adaptation to contrasting climatic zones in driving speciation. This research provides insight into the process of speciation in a group of ecologically important plants and strengthens our predictive understanding of how plant populations will respond as Earth’s climate rapidly changes.

58 GEOSCIENCES↗