Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “evolutionary biology”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Population Genomic Analyses Reveal Connectivity via Human-Mediated Transport across Populus Plantations in North America and an Undescribed Subpopulation of Sphaerulina musiva

Domestication of plant species has affected the evolutionary dynamics of plant pathogens in agriculture and forestry. A model system for studying the consequences of plant domestication on the evolution of an emergent plant disease is the fungal pathogen Sphaerulina musiva. This ascomycete causes leaf spot and stem canker disease of Populus spp. and their hybrids. A population genomics approach was used to determine the degree of population structure and evidence for selection on the North American population of S. musiva. In total, 122 samples of the fungus were genotyped identifying 120,016 single-nucleotide polymorphisms after quality filtering. In North America, S. musiva has low to moderate degrees of differentiation among locations. Three main genetic clusters were detected: southeastern United States, midwestern United States and Canada, and a new British Columbia cluster (BC2). Population genomics suggest that BC2 is a novel genetic cluster from central British Columbia, clearly differentiated from previously reported S. musiva from coastal British Columbia, and the product of a single migration event. Phenotypic measurements from greenhouse experiments indicate lower aggressiveness of BC2 on Populus trichocarpa. In summary, S. musiva has geographic structure across broad regions indicative of gene flow among clusters. The interconnectedness of the North American S. musiva populations across large geographic distances further supports the hypothesis of anthropogenic-facilitated transport of the pathogen.

Biochemistry & Molecular Biology↗
DOE OSTI2020

Mycena genomes resolve the evolution of fungal bioluminescence

Mushroom-forming fungi in the order Agaricales represent an independent origin of bioluminescence in the tree of life; yet the diversity, evolutionary history, and timing of the origin of fungal luciferases remain elusive. We sequenced the genomes and transcriptomes of five bonnet mushroom species ( Mycena spp.), a diverse lineage comprising the majority of bioluminescent fungi. Two species with haploid genome assemblies ~150 Mb are among the largest in Agaricales, and we found that a variety of repeats between Mycena species were differentially mediated by DNA methylation. We show that bioluminescence evolved in the last common ancestor of mycenoid and the marasmioid clade of Agaricales and was maintained through at least 160 million years of evolution. Analyses of synteny across genomes of bioluminescent species resolved how the luciferase cluster was derived by duplication and translocation, frequently rearranged and lost in most Mycena species, but conserved in the Armillaria lineage. Luciferase cluster members were coexpressed across developmental stages, with the highest expression in fruiting body caps and stipes, suggesting fruiting-related adaptive functions. Our results contribute to understanding a de novo origin of bioluminescence and the corresponding gene cluster in a diverse group of enigmatic fungal species.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Gene family expansions and transcriptome signatures uncover fungal adaptations to wood decay

Because they comprise some of the most efficient wood-decayers, Polyporales fungi impact carbon cycling in forest environment. Despite continuous discoveries on the enzymatic machinery involved in wood decomposition, the vision on their evolutionary adaptation to wood decay and genome diversity remains incomplete. We combined the genome sequence information from 50 Polyporales species, including 26 newly sequenced genomes and sought for genomic and functional adaptations to wood decay through the analysis of genome composition and transcriptome responses to different carbon sources. The genomes of Polyporales from different phylogenetic clades showed poor conservation in macrosynteny, indicative of genome rearrangements. We observed different gene family expansion/contraction histories for plant cell wall degrading enzymes in core polyporoids and phlebioids and captured expansions for genes involved in signalling and regulation in the lineages of white rotters. Furthermore, we identified conserved cupredoxins, thaumatin-like proteins and lytic polysaccharide monooxygenases with a yet uncharacterized appended module as new candidate players in wood decomposition. Given the current need for enzymatic toolkits dedicated to the transformation of renewable carbon sources, the observed genomic diversity among Polyporales strengthens the relevance of mining Polyporales biodiversity to understand the molecular mechanisms of wood decay.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Functional Divergence in a Multi-gene Family Is a Key Evolutionary Innovation for Anaerobic Growth in Saccharomyces cerevisiae

The amplification and diversification of genes into large multi-gene families often mark key evolutionary innovations, but this process often creates genetic redundancy that hinders functional investigations. When the model budding yeast Saccharomyces cerevisiae transitions to anaerobic growth conditions, the cell massively induces the expression of seven serine/threonine-rich anaerobically-induced cell wall mannoproteins (anCWMPs): TIP1, TIR1, TIR2, TIR3, TIR4, DAN1, and DAN4. Here, we show that these genes likely derive evolutionarily from a single ancestral anCWMP locus, which was duplicated and translocated to new genomic contexts several times both prior to and following the budding yeast whole genome duplication (WGD) event. Based on synteny and their phylogeny, we separate the anCWMPs into four gene subfamilies. To resolve prior inconclusive genetic investigations of these genes, we constructed a set of combinatorial deletion mutants to determine their contributions toward anaerobic growth in S. cerevisiae. We found that two genes, TIR1 and TIR3, were together necessary and sufficient for the anCWMP contribution to anaerobic growth. Overexpressing either gene alone was insufficient for anaerobic growth, implying that they encode non-overlapping functional roles in the cell during anaerobic growth. We infer from the phylogeny of the anCWMP genes that these two important genes derive from an ancient duplication that predates the WGD event, whereas the TIR1 subfamily experienced gene family amplification after the WGD event. Taken together, the genetic and molecular evidence suggests that one key anCWMP gene duplication event, several auxiliary gene duplication events, and functional divergence underpin the evolution of anaerobic growth in budding yeasts.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Major facilitator family transporters specifically enhance caffeyl alcohol uptake during C‐lignin biosynthesis

The mode of transport of lignin monomers to the sites of polymerization in the apoplast remains controversial. C-Lignin is a recently discovered form of lignin found in some seed coats that is composed exclusively of units derived from caffeyl alcohol. RNA-seq and proteome analyses identified a number of transporters co-expressed with C-lignin deposition in the seed coat of Cleome hassleriana. Cloning and influx/efflux analysis assays in yeast identified two low-affinity transporters, ChPLT3 and ChSUC1, that were active with caffeyl alcohol but not with the classical monolignols p-coumaryl, coniferyl, and sinapyl alcohols, consistent with molecular modeling and docking studies. Expression of ChPLT3 in Arabidopsis seedlings enhanced root growth in the presence of caffeyl alcohol, and expression of ChPLT3 and ChSUC1 correlated with lignin C-unit content in hairy roots of Medicago truncatula. We present a model, consistent with phylogenetic and evolutionary considerations, whereby passive caffeyl alcohol transport may be supplemented by hitchhiking on secondary active transporters to ensure the synthesis of C-lignin, and inhibition of synthesis of G-lignin, in the apoplast.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Functional and Evolutionary Integration of a Fungal Gene With a Bacterial Operon

Abstract Siderophores are crucial for iron-scavenging in microorganisms. While many yeasts can uptake siderophores produced by other organisms, they are typically unable to synthesize siderophores themselves. In contrast, Wickerhamiella/Starmerella (W/S) clade yeasts gained the capacity to make the siderophore enterobactin following the remarkable horizontal acquisition of a bacterial operon enabling enterobactin synthesis. Yet, how these yeasts absorb the iron bound by enterobactin remains unresolved. Here, we demonstrate that Enb1 is the key enterobactin importer in the W/S-clade species Starmerella bombicola. Through phylogenomic analyses, we show that ENB1 is present in all W/S clade yeast species that retained the enterobactin biosynthetic genes. Conversely, it is absent in species that lost the ent genes, except for Starmerella stellata, making this species the only cheater in the W/S clade that can utilize enterobactin without producing it. Through phylogenetic analyses, we infer that ENB1 is a fungal gene that likely existed in the W/S clade prior to the acquisition of the ent genes and subsequently experienced multiple gene losses and duplications. Through phylogenetic topology tests, we show that ENB1 likely underwent horizontal gene transfer from an ancient W/S clade yeast to the order Saccharomycetales, which includes the model yeast Saccharomyces cerevisiae, followed by extensive secondary losses. Taken together, these results suggest that the fungal ENB1 and bacterial ent genes were cooperatively integrated into a functional unit within the W/S clade that enabled adaptation to iron-limited environments. This integrated fungal-bacterial circuit and its dynamic evolution determine the extant distribution of yeast enterobactin producers and cheaters.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Targeting tRNA-synthetase interactions towards novel therapeutic discovery against eukaryotic pathogens

The development of chemotherapies against eukaryotic pathogens is especially challenging because of both the evolutionary conservation of drug targets between host and parasite, and the evolution of strain-dependent drug resistance. There is a strong need for new nontoxic drugs with broad-spectrum activity against trypanosome parasites such as Leishmania and Trypanosoma. A relatively untested approach is to target macromolecular interactions in parasites rather than small molecular interactions, under the hypothesis that the features specifying macromolecular interactions diverge more rapidly through coevolution. We computed tRNA Class-Informative Features in humans and independently in eight distinct clades of trypanosomes, identifying parasite-specific informative features, including base pairs and base mis-pairs, that are broadly conserved over approximately 250 million years of trypanosome evolution. Validating these observations, we demonstrated biochemically that tRNA:aminoacyl-tRNA synthetase (aaRS) interactions are a promising target for anti-trypanosomal drug discovery. From a marine natural products extract library, we identified several fractions with inhibitory activity toward Leishmania major alanyl-tRNA synthetase (AlaRS) but no activity against the human homolog. These marine natural products extracts showed cross-reactivity towards Trypanosoma cruzi AlaRS indicating the broad-spectrum potential of our network predictions. We also identified Leishmania major threonyl-tRNA synthetase (ThrRS) inhibitors from the same library. We discuss why chemotherapies targeting multiple aaRSs should be less prone to the evolution of resistance than monotherapeutic or synergistic combination chemotherapies targeting only one aaRS.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Intrinsic and environmental drivers of pairwise cohesion in wild Canis social groups

Animals within social groups respond to costs and benefits of sociality by adjusting the proportion of time they spend in close proximity to other individuals in the group (cohesion). Variation in cohesion between individuals, in turn, shapes important group-level processes such as subgroup formation and fission–fusion dynamics. Although critical to animal sociality, a comprehensive understanding of the factors influencing cohesion remains a gap in our knowledge of cooperative behavior in animals. We tracked 574 individuals from six species within the genus Canis in 15 countries on four continents with GPS telemetry to estimate the time that pairs of individuals within social groups spent in close proximity and test hypotheses regarding drivers of cohesion. Pairs of social canids (Canis spp.) varied widely in the proportion of time they spent together (5%–100%) during seasonal monitoring periods relative to both intrinsic characteristics and environmental conditions. The majority of our data came from three species of wolves (gray wolves, eastern wolves, and red wolves) and coyotes. For these species, cohesion within social groups was greatest between breeding pairs and varied seasonally as the nature of cooperative activities changed relative to annual life history patterns. Across species, wolves were more cohesive than coyotes. For wolves, pairs were less cohesive in larger groups, and when suitable, small prey was present reflecting the constraints of food resources and intragroup competition on social associations. Pair cohesion in wolves declined with increased anthropogenic modification of the landscape and greater climatic variability, underscoring challenges for conserving social top predators in a changing world. We show that pairwise cohesion in social groups varies strongly both within and across Canis species, as individuals respond to changing ecological context defined by resources, competition, and anthropogenic disturbance. Our work highlights that cohesion is a highly plastic component of animal sociality that holds significant promise for elucidating ecological and evolutionary mechanisms underlying cooperative behavior.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Host population diversity as a driver of viral infection cycle in wild populations of green sulfur bacteria with long standing virus-host interactions

Abstract Temperate phages are viruses of bacteria that can establish two types of infection: a lysogenic infection in which the virus replicates with the host cell without producing virions, and a lytic infection where the host cell is eventually destroyed, and new virions are released. While both lytic and lysogenic infections are routinely observed in the environment, the ecological and evolutionary processes regulating these viral dynamics are still not well understood, especially for uncultivated virus-host pairs. Here, we characterized the long-term dynamics of uncultivated viruses infecting green sulfur bacteria (GSB) in a model freshwater lake (Trout Bog Lake, TBL). As no GSB virus has been formally described yet, we first used two complementary approaches to identify new GSB viruses from TBL; one in vitro based on flow cytometry cell sorting, the other in silico based on CRISPR spacer sequences. We then took advantage of existing TBL metagenomes covering the 2005–2018 period to examine the interactions between GSB and their viruses across years and seasons. From our data, GSB populations in TBL were constantly associated with at least 2-8 viruses each, including both lytic and temperate phages. The dominant GSB population in particular was consistently associated with two prophages with a nearly 100% infection rate for >10 years. We illustrate with a theoretical model that such an interaction can be stable given a low, but persistent, level of prophage induction in low-diversity host populations. Overall, our data suggest that lytic and lysogenic viruses can readily co-infect the same host population, and that host strain-level diversity might be an important factor controlling virus-host dynamics including lytic/lysogeny switch.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

End-to-end learning of multiple sequence alignments with differentiable Smith–Waterman

Abstract Motivation Multiple sequence alignments (MSAs) of homologous sequences contain information on structural and functional constraints and their evolutionary histories. Despite their importance for many downstream tasks, such as structure prediction, MSA generation is often treated as a separate pre-processing step, without any guidance from the application it will be used for. Results Here, we implement a smooth and differentiable version of the Smith–Waterman pairwise alignment algorithm that enables jointly learning an MSA and a downstream machine learning system in an end-to-end fashion. To demonstrate its utility, we introduce SMURF (Smooth Markov Unaligned Random Field), a new method that jointly learns an alignment and the parameters of a Markov Random Field for unsupervised contact prediction. We find that SMURF learns MSAs that mildly improve contact prediction on a diverse set of protein and RNA families. As a proof of concept, we demonstrate that by connecting our differentiable alignment module to AlphaFold2 and maximizing predicted confidence, we can learn MSAs that improve structure predictions over the initial MSAs. Interestingly, the alignments that improve AlphaFold predictions are self-inconsistent and can be viewed as adversarial. This work highlights the potential of differentiable dynamic programming to improve neural network pipelines that rely on an alignment and the potential dangers of optimizing predictions of protein sequences with methods that are not fully understood. Availability and implementation Our code and examples are available at: https://github.com/spetti/SMURF. Supplementary information Supplementary data are available at Bioinformatics online.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Snekmer: a scalable pipeline for protein sequence fingerprinting based on amino acid recoding

Abstract Motivation The vast expansion of sequence data generated from single organisms and microbiomes has precipitated the need for faster and more sensitive methods to assess evolutionary and functional relationships between proteins. Representing proteins as sets of short peptide sequences (kmers) has been used for rapid, accurate classification of proteins into functional categories; however, this approach employs an exact-match methodology and thus may be limited in terms of sensitivity and coverage. We have previously used similarity groupings, based on the chemical properties of amino acids, to form reduced character sets and recode proteins. This amino acid recoding (AAR) approach simplifies the construction of protein representations in the form of kmer vectors, which can link sequences with distant sequence similarity and provide accurate classification of problematic protein families. Results Here, we describe Snekmer, a software tool for recoding proteins into AAR kmer vectors and performing either (i) construction of supervised classification models trained on input protein families or (ii) clustering for de novo determination of protein families. We provide examples of the operation of the tool against a set of nitrogen cycling families originally collected using both standard hidden Markov models and a larger set of proteins from Uniprot and demonstrate that our method accurately differentiates these sequences in both operation modes. Availability and implementation Snekmer is written in Python using Snakemake. Code and data used in this article, along with tutorial notebooks, are available at http://github.com/PNNL-CompBio/Snekmer under an open-source BSD-3 license. Supplementary information Supplementary data are available at Bioinformatics Advances online.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

The contributions from the progenitor genomes of the mesopolyploid Brassiceae are evolutionarily distinct but functionally compatible

The members of the tribe Brassiceae share a whole-genome triplication (WGT), and one proposed model for its formation is a two-step pair of hybridizations producing hexaploid descendants. However, evidence for this model is incomplete, and the evolutionary and functional constraints that drove evolution after the hexaploidy are even less understood. Here, we report a new genome sequence of Crambe hispanica , a species sister to most sequenced Brassiceae. Using this new genome and three others that share the hexaploidy, we traced the history of gene loss after the WGT using the Polyploidy Orthology Inference Tool (POInT). We confirm the two-step formation model and infer that there was a significant temporal gap between those two allopolyploidizations, with about a third of the gene losses from the first two subgenomes occurring before the arrival of the third. We also, for the 90,000 individual genes in our study, make parental subgenome assignments, inferring, with measured uncertainty, from which of the progenitor genomes of the allohexaploidy each gene derives. We further show that each subgenome has a statistically distinguishable rate of homoeolog losses. There is little indication of functional distinction between the three subgenomes: the individual subgenomes show no patterns of functional enrichment, no excess of shared protein-protein or metabolic interactions between their members, and no biases in their likelihood of having experienced a recent selective sweep. We propose a "mix and match" model of allopolyploidy, in which subgenome origin drives homoeolog loss propensities but where genes from different subgenomes function together without difficulty.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Quantifying microbial control of soil organic matter dynamics at macrosystem scales

Soil organic matter (SOM) stocks, decomposition and persistence are largely the product of controls that act locally. Yet the controls are shaped and interact at multiple spatiotemporal scales, from which macrosystem patterns in SOM emerge. Theory on SOM turnover recognizes the resulting spatial and temporal conditionality in the effect sizes of controls that play out across macrosystems, and couples them through evolutionary and community assembly processes. For example, climate history shapes plant functional traits, which in turn interact with contemporary climate to influence SOM dynamics. Selection and assembly also shape the functional traits of soil decomposer communities, but it is less clear how in turn these traits influence temporal macrosystem patterns in SOM turnover. Here, we review evidence that establishes the expectation that selection and assembly should generate decomposer communities across macrosystems that have distinct functional effects on SOM dynamics. Representation of this knowledge in soil biogeochemical models affects the magnitude and direction of projected SOM responses under global change. Yet there is high uncertainty and low confidence in these projections. To address these issues, we make the case that a coordinated set of empirical practices are required which necessitate (1) greater use of statistical approaches in biogeochemistry that are suited to causative inference; (2) long-term, macrosystem-scale, observational and experimental networks to reveal conditionality in effect sizes, and embedded correlation, in controls on SOM turnover; and (3) use of multiple measurement grains to capture local- and macroscale variation in controls and outcomes, to avoid obscuring causative understanding through data aggregation. Here, when employed together, along with process-based models to synthesize knowledge and guide further empirical work, we believe these practices will rapidly advance understanding of microbial controls on SOM and improve carbon cycle projections that guide policies on climate adaptation and mitigation.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1991

A model of onshore-offshore change in faunal diversity

Onshore-offshore patterns of faunal change occurred at many taxonomic scales during the Paleozoic Era, ranging from replacement of the Cambrian evolutionary fauna by the Paleozoic fauna to the environmental expansion of many orders and classes. A simple mathematical model is constructed to investigate such change. The environmental gradient across the marine shelf-slope is treated as a linear array of discrete habitats, each of which holds a set number of species, as observed in the fossil record. During any interval of time, some portion of the species in each habitat becomes extinct by background processes, with rates of extinction varying among both clades and habitats, as also observed in the record. After extinction, species are replaced from within the habitat and from immediately adjacent habitats, with proportions dependent on surviving species. This model leads to the prediction that extinction-resistant clades will always diversify at the expense of extinction-prone clades. But if extinction intensity is highest in nearshore habitats, extinction-resistant clades will expand preferentially in the onshore direction, build up diversity there, and then diversify outward toward the offshore. Thus, onshore-offshore patterns of diversification may be the expectation for faunal change quite independently of whether or not clades originate onshore. When the model is parameterized for Paleozoic trilobites and brachiopods, numerical solutions exhibit both a pattern of faunal change and a time span for diversification similar to that seen in the fossil record. They also generate structure similar to that seen in global diversification, including logistic patterns of growth, declining origination but constant extinction within clades through time, and declining overall extinction across clades through time.

NASA Discipline Number 52-40↗
DOE OSTI2026

Pyrodictium abyssi AbpX reveals a calcium-responsive family of microbial biomatrix proteins that form thermostable hydrogels

Evolutionary pressure on microbial communities propagating under extreme environmental conditions often results in unique structural adaptations to promote cell survival. In this work, we report an investigation of AbpX, a biomatrix protein identified in cultures of the hyperthermophilic archaeon Pyrodictium abyssi. Under ex vivo and in vitro conditions, AbpX assembles into a paracrystalline lattice composed of semiflexible fibrils. CryoEM analysis of recombinant AbpX fibrils reveals that the precursor protein polymerizes through donor strand complementation (DSC), a process previously reported for chaperone-usher fimbriae in Gram-negative bacteria. Unlike the latter DSC protein polymers, AbpX undergoes chaperone-free polymerization in the presence of calcium ions, which are sequestered at the donor strand-acceptor groove interface between protomers in the fibril. Using a combination of cryoEM and crystallographic information, a structural model is proposed for the AbpX lattice that provides insight into its potential role in biofilm formation. These findings suggest that calcium ion coordination may contribute to fibril assembly and preorganize fibrils for incorporation into the protein lattice. Bioinformatic analysis indicates that AbpX exemplifies a distinct and broadly distributed clade of calcium ion responsive biomatrix proteins within the TasA superfamily that can be fabricated into hydrogel biomaterials in vitro under environmentally benign conditions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Two major chromosome evolution events with unrivaled conserved gene content in pomegranate

Pomegranate has a unique evolutionary history given that different cultivars have eight or nine bivalent chromosomes with possible crossability between the two classes. Therefore, it is important to study chromosome evolution in pomegranate to understand the dynamics of its population. Here, we de novo assembled the Azerbaijani cultivar “Azerbaijan guloyshasi” (AG2017; 2n = 16) and re-sequenced six cultivars to track the evolution of pomegranate and to compare it with previously published de novo assembled and re-sequenced cultivars. High synteny was observed between AG2017, Bhagawa (2n = 16), Tunisia (2n = 16), and Dabenzi (2n = 18), but these four cultivars diverged from the cultivar Taishanhong (2n = 18) with several rearrangements indicating the presence of two major chromosome evolution events. Major presence/absence variations were not observed as >99% of the five genomes aligned across the cultivars, while >99% of the pan-genic content was represented by Tunisia and Taishanhong only. We also revisited the divergence between soft- and hard-seeded cultivars with less structured population genomic data, compared to previous studies, to refine the selected genomic regions and detect global migration routes for pomegranate. We reported a unique admixture between soft- and hard-seeded cultivars that can be exploited to improve the diversity, quality, and adaptability of local pomegranate varieties around the world. Our study adds body knowledge to understanding the evolution of the pomegranate genome and its implications for the population structure of global pomegranate diversity, as well as planning breeding programs aiming to develop improved cultivars.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Ancient Rapid Radiation Explains Most Conflicts Among Gene Trees and Well-Supported Phylogenomic Trees of Nostocalean Cyanobacteria

Prokaryotic genomes are often considered to be mosaics of genes that do not necessarily share the same evolutionary history due to widespread horizontal gene transfers (HGTs). Consequently, representing evolutionary relationships of prokaryotes as bifurcating trees has long been controversial. However, studies reporting conflicts among gene trees derived from phylogenomic data sets have shown that these conflicts can be the result of artifacts or evolutionary processes other than HGT, such as incomplete lineage sorting, low phylogenetic signal, and systematic errors due to substitution model misspecification. Here, we present the results of an extensive exploration of phylogenetic conflicts in the cyanobacterial order Nostocales, for which previous studies have inferred strongly supported conflicting relationships when using different concatenated phylogenomic data sets. We found that most of these conflicts are concentrated in deep clusters of short internodes of the Nostocales phylogeny, where the great majority of individual genes have low resolving power. We then inferred phylogenetic networks to detect HGT events while also accounting for incomplete lineage sorting. Our results indicate that most conflicts among gene trees are likely due to incomplete lineage sorting linked to an ancient rapid radiation, rather than to HGTs. Moreover, the short internodes of this radiation fit the expectations of the anomaly zone, i.e., a region of the tree parameter space where a species tree is discordant with its most likely gene tree. In this work, we demonstrated that concatenation of different sets of loci can recover up to 17 distinct and well-supported relationships within the putative anomaly zone of Nostocales, corresponding to the observed conflicts among well-supported trees based on concatenated data sets from previous studies. Our findings highlight the important role of rapid radiations as a potential cause of strongly conflicting phylogenetic relationships when using phylogenomic data sets of bacteria. We propose that polytomies may be the most appropriate phylogenetic representation of these rapid radiations that are part of anomaly zones, especially when all possible genomic markers have been considered to infer these phylogenies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Genotype–environment associations across spatial scales reveal the importance of putative adaptive genetic variation in divergence

Abstract Identifying areas of high evolutionary potential is a judicious strategy for developing conservation priorities in the face of environmental change. For wide‐ranging species occupying heterogeneous environments, the evolutionary forces that shape distinct populations can vary spatially. Here, we investigate patterns of genomic variation and genotype–environment associations in the hermit thrush ( Catharus guttatus ), a North American songbird, at broad (across the breeding range) and narrow spatial scales (at a hybrid zone). We begin by building a genoscape or map of genetic variation across the breeding range and find five distinct genetic clusters within the species, with the greatest variation occurring in the western portion of the range. Genotype–environment association analyses indicate higher allelic turnover in the west than in the east, with measures of temperature surfacing as key predictors of putative adaptive genomic variation rangewide. Since broad patterns detected across a species' range represent the aggregate of many locally adapted populations, we investigate whether our broadscale analysis is consistent with a finer scale analysis. We find that top rangewide temperature‐associated loci vary in their clinal patterns (e.g., steep clines vs. fixed allele frequencies) across a hybrid zone in British Columbia, suggesting that the environmental predictors and the associated candidate loci identified in the rangewide analysis are of variable importance in this particular region. However, two candidate loci exhibit strong concordance with the temperature gradient in British Columbia, suggesting a potential role for temperature‐related barriers to gene flow and/or temperature‐driven ecological selection in maintaining putative local adaptation. This study demonstrates how patterns identified at the broad (macrogeographic) scale can be validated by investigating genotype–environment correlations at the local (microgeographic) scale. Furthermore, our results highlight the importance of considering the spatial distribution of putative adaptive variation when assessing population‐level sensitivity to climate change and other stressors.

59 BASIC BIOLOGICAL SCIENCES↗