Search NASA⌕ Search

SEARCH · Search NASA

Results for “sequence divergence”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Ancestral archaea expanded the genetic code with pyrrolysine

The pyrrolysyl-tRNA synthetase (PylRS) facilitates the cotranslational installation of the 22nd amino acid pyrrolysine. Owing to its tolerance for diverse amino acid substrates, and its orthogonality in multiple organisms, PylRS has emerged as a major route to install noncanonical amino acids into proteins in living cells. Recently, a novel class of PylRS enzymes was identified in a subset of methanogenic archaea. Enzymes within this class (ΔPylSn) lack the N-terminal tRNA-binding domain that is widely conserved amongst PylRS enzymes, yet remain active and orthogonal in bacteria and eukaryotes. In this study, we use biochemical and in vivo UAG-readthrough assays to characterize the aminoacylation efficiency and substrate spectrum of a ΔPylSn class PylRS from the archaeon Candidatus Methanomethylophilus alvus. We show that, compared with the full-length enzyme from Methanosarcina mazei, the Ca. M. alvus PylRS displays reduced aminoacylation efficiency but an expanded amino acid substrate spectrum. To gain insight into the evolution of ΔPylSn enzymes, we performed molecular phylogeny using 156 PylRS and 105 pyrrolysine tRNA (tRNA Pyl ) sequences from diverse archaea and bacteria. This analysis suggests that the PylRS•tRNA Pyl pair diverged before the evolution of the three domains of life, placing an early limit on the evolution of the Pyl-decoding trait. Furthermore, our results document the coevolutionary history of PylRS and tRNA Pyl and reveal the emergence of tRNA Pyl sequences with unique A73 and U73 discriminator bases. The orthogonality of these tRNA Pyl species with the more common G73-containing tRNA Pyl will enable future efforts to engineer PylRS systems for further genetic code expansion.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Atypical Divergence of SARS-CoV-2 Orf8 from Orf7a within the Coronavirus Lineage Suggests Potential Stealthy Viral Strategies in Immune Evasion

Orf8, one of the most puzzling genes in the SARS lineage of coronaviruses, marks a unique and striking difference in genome organization between SARS-CoV-2 and SARS-CoV-1. Here, using sequence comparisons, we unequivocally reveal the distant sequence similarities between SARS-CoV-2 Orf8 with its SARS-CoV-1 counterparts and the X4-like genes of coronaviruses, including its highly divergent “paralog” gene Orf7a, whose product is a potential immune antagonist of known structure. Supervised sequence space walks unravel identity levels that drop below 10% and yet exhibit subtle conservation patterns in this novel superfamily, characterized by an immunoglobulin-like beta sandwich topology. We document the high accuracy of the sequence space walk process in detail and characterize the subgroups of the superfamily in sequence space by systematic annotation of gene and taxon groups. While SARS-CoV-1 Orf7a and Orf8 genes are most similar to bat virus sequences, their SARS-CoV-2 counterparts are closer to pangolin virus homologs, reflecting the fine structure of conservation patterns within the SARS-CoV-2 genomes. The divergence between Orf7a and Orf8 is exceptionally idiosyncratic, since Orf7a is more constrained, whereas Orf8 is subject to rampant change, a peculiar feature that may be related to hitherto-unknown viral infection strategies. Despite their common origin, the Orf7a and Orf8 protein families exhibit different modes of evolutionary trajectories within the coronavirus lineage, which might be partly attributable to their complex interactions with the mammalian host cell, reflected by a multitude of functional associations of Orf8 in SARS-CoV-2 compared to a very small number of interactions discovered for Orf7a.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1994

Transcription factor IID in the Archaea: sequences in the Thermococcus celer genome would encode a product closely related to the TATA-binding protein of eukaryotes

The first step in transcription initiation in eukaryotes is mediated by the TATA-binding protein, a subunit of the transcription factor IID complex. We have cloned and sequenced the gene for a presumptive homolog of this eukaryotic protein from Thermococcus celer, a member of the Archaea (formerly archaebacteria). The protein encoded by the archaeal gene is a tandem repeat of a conserved domain, corresponding to the repeated domain in its eukaryotic counterparts. Molecular phylogenetic analyses of the two halves of the repeat are consistent with the duplication occurring before the divergence of the archael and eukaryotic domains. In conjunction with previous observations of similarity in RNA polymerase subunit composition and sequences and the finding of a transcription factor IIB-like sequence in Pyrococcus woesei (a relative of T. celer) it appears that major features of the eukaryotic transcription apparatus were well-established before the origin of eukaryotic cellular organization. The divergence between the two halves of the archael protein is less than that between the halves of the individual eukaryotic sequences, indicating that the average rate of sequence change in the archael protein has been less than in its eukaryotic counterparts. To the extent that this lower rate applies to the genome as a whole, a clearer picture of the early genes (and gene families) that gave rise to present-day genomes is more apt to emerge from the study of sequences from the Archaea than from the corresponding sequences from eukaryotes.

NASA Discipline Exobiology↗
DOE OSTI2023

Ecological divergence of syntopic marine bacterial species is shaped by gene content and expression

Identifying mechanisms by which bacterial species evolve and maintain genomic diversity is particularly challenging for the uncultured lineages that dominate the surface ocean. A longitudinal analysis of bacterial genes, genomes, and transcripts during a coastal phytoplankton bloom revealed two co-occurring, highly related Rhodobacteraceae species from the deeply branching and uncultured NAC11-7 lineage. These have identical 16S rRNA gene amplicon sequences, yet their genome contents assembled from metagenomes and single cells indicate species-level divergence. Moreover, shifts in relative dominance of the species during dynamic bloom conditions over 7 weeks confirmed the syntopic species’ divergent responses to the same microenvironment at the same time. Genes unique to each species and genes shared but divergent in per-cell inventories of mRNAs accounted for 5% of the species’ pangenome content. These analyses uncover physiological and ecological features that differentiate the species, including capacities for organic carbon utilization, attributes of the cell surface, metal requirements, and vitamin biosynthesis. Such insights into the coexistence of highly related and ecologically similar bacterial species in their shared natural habitat are rare.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2001

Evolution of hematopoiesis: Three members of the PU.1 transcription factor family in a cartilaginous fish, Raja eglanteria

T lymphocytes and B lymphocytes are present in jawed vertebrates, including cartilaginous fishes, but not in jawless vertebrates or invertebrates. The origins of these lineages may be understood in terms of evolutionary changes in the structure and regulation of transcription factors that control lymphocyte development, such as PU.1. The identification and characterization of three members of the PU.1 family of transcription factors in a cartilaginous fish, Raja eglanteria, are described here. Two of these genes are orthologs of mammalian PU.1 and Spi-C, respectively, whereas the third gene, Spi-D, is a different family member. In addition, a PU.1-like gene has been identified in a jawless vertebrate, Petromyzon marinus (sea lamprey). Both DNA-binding and transactivation domains are highly conserved between mammalian and skate PU.1, in marked contrast to lamprey Spi, in which similarity is evident only in the DNA-binding domain. Phylogenetic analysis of sequence data suggests that the appearance of Spi-C may predate the divergence of the jawed and jawless vertebrates and that Spi-D arose before the divergence of the cartilaginous fish from the lineage leading to the mammals. The tissue-specific expression patterns of skate PU.1 and Spi-C suggest that these genes share regulatory as well as structural properties with their mammalian orthologs.

NASA Discipline Evolutionary Biology↗
DOE OSTI2018

Crystal structure of the RxLR effector PcRxLR12 from Phytophthora capsici

RxLR genes are a prominent class of effectors in oomycetes, and almost half of these proteins contain a conserved sequence motif termed the WY domain, that may exist singly, or as divergent tandem repeats in different effectors. Here we describe the crystal structure of PcRxLR12 (63-488) from Phytophthora capsici at 3.0 Å resolution. The structure consists of five tandemly arrayed WY-domains linked to each other by short connecting helices. Superposition of the WY-2 domain on the other four domains of PcRxLR12, show that the first α-helix termed the K motif, and Loop 3 which connects α{sub 3} and α{sub 4} are the key regions of structural divergence between the WY domains. A similar pattern was observed when WY-2 was superposed on the 11 WY domains from other oomycete effectors. We also note that an added connecting helix between WY domains in some RXLR effectors, ensures that the WY domains are oriented in the same direction.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Deciphering the evolution of flavin-dependent monooxygenase stereoselectivity using ancestral sequence reconstruction

Controlling the selectivity of a reaction is critical for target-oriented synthesis. Accessing complementary selectivity profiles enables divergent synthetic strategies, but is challenging to achieve in biocatalytic reactions given enzymes’ innate preferences of a single selectivity. Thus, it is critical to understand the structural features that control selectivity in biocatalytic reactions to achieve tunable selectivity. Here, we investigate the structural features that control the stereoselectivity in an oxidative dearomatization reaction that is key to making azaphilone natural products. Crystal structures of enantiocomplementary biocatalysts guided the development of multiple hypotheses centered on the structural features that control the stereochemical outcome of the reaction; however, in many cases, direct substitutions of active site residues in natural proteins led to inactive enzymes. Ancestral sequence reconstruction (ASR) and resurrection were employed as an alternative strategy to probe the impact of each residue on the stereochemical outcome of the dearomatization reaction. These studies suggest that two mechanisms are active in controlling the stereochemical outcome of the oxidative dearomatization reaction: one involving multiple active site residues in AzaH and the other dominated by a single Phe to Tyr switch in TropB and AfoD. Moreover, this study suggests that the flavin-dependent monooxygenases (FDMOs) adopt simple and flexible strategies to control stereoselectivity, which has led to stereocomplementary azaphilone natural products produced by fungi. This paradigm of combining ASR and resurrection with mutational and computational studies showcases sets of tools for understanding enzyme mechanisms and provides a solid foundation for future protein engineering efforts.

Science & Technology - Other Topics↗
DOE OSTI2025

Leaky ribosomal scanning enables tunable translation of bicistronic ORFs in green algae

Advances in sequencing technology have unveiled examples of nucleus-encoded polycistrons, once considered rare. Exclusively polycistronic transcripts are prevalent in green algae, although the mechanism by which multiple polypeptides are translated from a single transcript is unknown. Here, we used bioinformatic and in vivo mutational analyses to evaluate competing mechanistic models for translation of bicistronic mRNAs in green algae. High-confidence manually curated datasets of bicistronic loci from two divergent green algae, Chlamydomonas reinhardtii and Auxenochlorella protothecoides, revealed a preference for weak Kozak-like sequences for ORF 1 and an underrepresentation of potential initiation codons before the ORF 2 start codon, which are suitable conditions for leaky ribosome scanning to allow ORF 2 translation. We used mutational analysis in A. protothecoides to test the mechanism. In vivo manipulation of the ORF 1 Kozak-like sequence and start codon altered reporter expression at ORF 2, with a weaker Kozak-like sequence enhancing expression and a stronger one diminishing it. A synthetic bicistronic dual reporter demonstrated inversely adjustable activity of green fluorescent protein expressed from ORF 1 and luciferase from ORF 2, depending on the strength of the ORF 1 Kozak-like sequence. Our findings demonstrate that translation of multiple ORFs in green algal bicistronic transcripts is consistent with episodic leaky scanning of ORF 1 to allow translation at ORF 2. This work has implications for the potential functionality of upstream open reading frames (uORFs) found across eukaryotic genomes and for transgene expression in synthetic biology applications.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1992

How close is close: 16S rRNA sequence identity may not be sufficient to guarantee species identity

16S rRNA (genes coding for rRNA) sequence comparisons were conducted with the following three psychrophilic strains: Bacillus globisporus W25T (T = type strain) and Bacillus psychrophilus W16AT, and W5. These strains exhibited more than 99.5% sequence identity and within experimental uncertainty could be regarded as identical. Their close taxonomic relationship was further documented by phenotypic similarities. In contrast, previously published DNA-DNA hybridization results have convincingly established that these strains do not belong to the same species if current standards are used. These results emphasize the important point that effective identity of 16S rRNA sequences is not necessarily a sufficient criterion to guarantee species identity. Thus, although 16S rRNA sequences can be used routinely to distinguish and establish relationships between genera and well-resolved species, very recently diverged species may not be recognizable.

Non-NASA Center↗
DOE OSTI2019

Rhizosphere microbiomes diverge among Populus trichocarpa plant-host genotypes and chemotypes, but it depends on soil origin

Abstract Background Plants have developed defense strategies for phytopathogen and herbivore protection via coordinated metabolic mechanisms. Low-molecular weight metabolites produced within plant tissues, such as salicylic acid, represent one such mechanism which likely mediates plant – microbe interactions above and below ground. Salicylic acid is a ubiquitous phytohormone at low levels in most plants, yet are concentrated defense compounds in Populus , likely acting as a selective filter for rhizosphere microbiomes. We propagated twelve Populus trichocarpa genotypes which varied an order of magnitude in salicylic acid (SA)-related secondary metabolites, in contrasting soils from two different origins. After four months of growth, plant properties (leaf growth, chlorophyll content, and net photosynthetic rate) and plant root metabolomics specifically targeting SA metabolites were measured via GC-MS. In addition, rhizosphere microbiome composition was measured via Illumina MiSeq sequencing of 16S and ITS2 rRNA-genes. Results Soil origin was the primary filter causing divergence in bacterial/archaeal and fungal communities with plant genotype secondarily influential. Both bacterial/archaeal and fungal evenness varied between soil origins and bacterial/archaeal diversity and evenness correlated with at least one SA metabolite (diversity: populin; evenness: total phenolics). The production of individual salicylic acid derivatives that varied by host genotype resulted in compositional differences for bacteria /archaea (tremuloidin) and fungi (salicylic acid) within one soil origin (Clatskanie) whereas soils from Corvallis did not illicit microbial compositional changes due to salicylic acid derivatives. Several dominant bacterial (e.g., Betaproteobacteria, Acidobacteria, Verrucomicrobia, Chloroflexi, Gemmatimonadete, Firmicutes ) and one fungal phyla ( Mortierellomycota ) also correlated with specific SA secondary metabolites; bacterial phyla exhibited more negative interactions (declining abundance with increasing metabolite concentration) than positive interactions. Conclusions These results indicate microbial communities diverge most among soil origin. However, within a soil origin, bacterial/archaeal communities are responsive to plant SA production within greenhouse-based rhizosphere microbiomes. Fungal microbiomes are impacted by root SA-metabolites, but overall to a lesser degree within this experimental context. These results suggest plant defense strategies, such as SA and its secondary metabolites, may partially drive patterns of both bacterial/archaeal and fungal taxa-specific colonization and assembly.

16S rRNA↗
DOE OSTI2018

Structure and specificity of several triclocarban-binding single domain camelid antibody fragments

The variable VHH domains of camelid single chain antibodies have been useful in numerous biotechnology applications due to their simplicity, biophysical properties, and abilities to bind to their cognate antigens with high affinities and specificity. Their interactions with proteins have been well-studied, but considerably less work has been done to characterize their ability to bind haptens. Here, a high-resolution structural study of three nanobodies (T4, T9, and T10) which have been shown to bind triclocarban (TCC, 3-(4-chlorophenyl)-1-(3,4-dichlorophenyl)urea) with near-nanomolar affinity shows that binding occurs in a tunnel largely formed by CDR1 rather than a surface or lateral binding mode seen in other nanobody-hapten interactions. Additional significant interactions are formed with a non-hypervariable loop, sometimes dubbed “CDR4”. A comparison of apo and holo forms of T9 and T10 shows that the binding site undergoes little conformational change upon binding of TCC. Structures of three nanobody-TCC complexes demonstrated there was not a standard binding mode. T4 and T9 have a high degree of sequence identity and bind the hapten in a nearly identical manner, while the more divergent T10 binds TCC in a slightly displaced orientation with the urea moiety rotated approximately 180° along the long axis of the molecule. In addition to methotrexate, this is the second report of haptens binding in a tunnel formed by CDR1, suggesting that compounds with similar hydrophobicity and shape could be recognized by nanobodies in analogous fashion. Structure-guided mutations failed to improve binding affinity for T4 and T9 underscoring the high degree of natural optimization.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Benefits and Limits of Phasing Alleles for Network Inference of Allopolyploid Complexes

Abstract Accurately reconstructing the reticulate histories of polyploids remains a central challenge for understanding plant evolution. Although phylogenetic networks can provide insights into relationships among polyploid lineages, inferring networks may be hindered by the complexities of homology determination in polyploid taxa. We use simulations to show that phasing alleles from allopolyploid individuals can improve phylogenetic network inference under the multispecies coalescent by obtaining the true network with fewer loci compared with haplotype consensus sequences or sequences with heterozygous bases represented as ambiguity codes. Phased allelic data can also improve divergence time estimates for networks, which is helpful for evaluating allopolyploid speciation hypotheses and proposing mechanisms of speciation. To achieve these outcomes in empirical data, we present a novel pipeline that leverages a recently developed phasing algorithm to reliably phase alleles from polyploids. This pipeline is especially appropriate for target enrichment data, where the depth of coverage is typically high enough to phase entire loci. We provide an empirical example in the North American Dryopteris fern complex that demonstrates insights from phased data as well as the challenges of network inference. We establish that our pipeline (PATÉ: Phased Alleles from Target Enrichment data) is capable of recovering a high proportion of phased loci from both diploids and polyploids. These data may improve network estimates compared with using haplotype consensus assemblies by accurately inferring the direction of gene flow, but statistical nonidentifiability of phylogenetic networks poses a barrier to inferring the evolutionary history of reticulate complexes.

Evolutionary Biology↗
DOE OSTI2017

Comparative genomics of Mortierella elongata and its bacterial endosymbiont Mycoavidus cysteinexigens

Endosymbiosis of bacteria by eukaryotes is a defining feature of cellular evolution. In addition to well-known bacterial origins for mitochondria and chloroplasts, multiple origins of bacterial endosymbiosis are known within the cells of diverse animals, plants and fungi. Early-diverging lineages of terrestrial fungi harbor endosymbiotic bacteria belonging to the Burkholderiaceae. Furthermore, we sequenced the metagenome of the soil-inhabiting fungus Mortierella elongata and assembled the complete circular chromosome of its endosymbiont, Mycoavidus cysteinexigens, which we place within a lineage of endofungal symbionts that are sister clade to Burkholderia. The genome of M. elongata strain AG77 features a core set of primary metabolic pathways for degradation of simple carbohydrates and lipid biosynthesis, while the M. cysteinexigens (AG77) genome is reduced in size and function. Experiments using antibiotics to cure the endobacterium from the host demonstrate that the fungal host metabolism is highly modulated by presence/ absence of M. cysteinexigens. In independent comparative phylogenomic analyses of fungal and bacterial genomes we find that they are consistent with an ancient origin for M. elongata M. cysteinexigens symbiosis, most likely over 350 million years ago and concomitant with the terrestrialization of Earth and diversification of land fungi and plants.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

cis -Regulatory Elements in Plant Development, Adaptation, and Evolution

cis-Regulatory elements encode the genomic blueprints that ensure the proper spatiotemporal patterning of gene expression necessary for appropriate development and responses to the environment. Accumulating evidence implicates changes to gene expression as a major source of phenotypic novelty in eukaryotes, including acute phenotypes such as disease and cancer in mammals. Moreover, genetic and epigenetic variation affecting cis-regulatory sequences over longer evolutionary timescales has become a recurring theme in studies of morphological divergence and local adaptation. Here, we discuss the functions of and methods used to identify various classes of cis-regulatory elements, as well as their role in plant development and response to the environment. We highlight opportunities to exploit cis-regulatory variants underlying plant development and environmental responses for crop improvement efforts. Although a comprehensive understanding of cis-regulatory mechanisms in plants has lagged behind that in animals, we showcase several breakthrough findings that have profoundly influenced plant biology and shaped the overall understanding of transcriptional regulation in eukaryotes.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1976

Solution to a gene divergence problem under arbitrary stable nucleotide transition probabilities

A nucleic acid chain, L nucleotides in length, with the specific base sequence B(1)B(2) ... B(L) is defined by the L-dimensional vector B = (B(1), B(2), ..., B(L)). For twelve given constant non-negative transition probabilities that, in a specified position, the base B is replaced by the base B' in a single step, an exact analytical expression is derived for the probability that the position goes from base B to B' in X steps. Assuming that each base mutates independently of the others, an exact expression is derived for the probability that the initial gene sequence B goes to a sequence B' = (B'(1), B'(2), ..., B'(L)) after X = (X(1), X(2), ..., X(L)) base replacements. The resulting equations allow a more precise accounting for the effects of Darwinian natural selection in molecular evolution than does the idealized (biologically less accurate) assumption that each of the four nucleotides is equally likely to mutate to and be fixed as one of the other three. Illustrative applications of the theory to some problems of biological evolution are given.

Holmquist, R.↗
DOE OSTI2023

The molecular basis of the neutralization breadth of the RBD-specific antibody CoV11

SARS-CoV-2, the virus behind the COVID-19 pandemic, has changed over time to the extent that the current virus is substantially different from what originally led to the pandemic in 2019–2020. Viral variants have modified the severity and transmissibility of the disease and continue do so. How much of this change is due to viral fitness versus a response to immune pressure is hard to define. One class of antibodies that continues to afford some level of protection from emerging variants are those that closely overlap the binding site for angiotensin-converting enzyme 2 (ACE2) on the receptor binding domain (RBD). Some members of this class that were identified early in the course of the pandemic arose from the VH 3-53 germline gene (IGHV3-53*01) and had short heavy chain complementarity-determining region 3s (CDR H3s). Here, we describe the molecular basis of the SARS-CoV-2 RBD recognition by the anti-RBD monoclonal antibody CoV11 isolated early in the COVID-19 pandemic and show how its unique mode of binding the RBD determines its neutralization breadth. CoV11 utilizes a heavy chain VH 3-53 and a light chain VK 3-20 germline sequence to bind to the RBD. Two of CoV11’s four heavy chain changes from the VH 3-53 germline sequence, $Thr^{28}_{FWRH1}$ to Ile and $Ser^{31}_{CDRH1}$ to Arg, and some unique features in its CDR H3 increase its affinity to the RBD, while the four light chain changes from the VK 3-20 germline sequence sit outside of the RBD binding site. Antibodies of this type can retain significant affinity and neutralization potency against variants of concern (VOCs) that have diverged significantly from original virus lineage such as the prevalent omicron variant. We also discuss the mechanism by which VH 3-53 encoded antibodies recognize spike antigen and show how minimal changes to their sequence, their choice of light chain, and their mode of binding influence their affinity and impact their neutralization breadth.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Novel Toxin Biosynthetic Gene Cluster in Harmful Algal Bloom-Causing Heteroscytonema crispum : Insights into the Origins of Paralytic Shellfish Toxins

Caused by both eukaryotic dinoflagellates and prokaryotic cyanobacteria, harmful algal blooms are events of severe ecological, economic, and public health consequence, and their incidence has become more common of late. Despite coordinated research efforts to identify and characterize the genomes of harmful algal bloom-causing organisms, the genomic basis and evolutionary origins of paralytic shellfish toxins produced by harmful algal blooms remain at best incomplete. The paralytic shellfish toxin saxitoxin has an especially complex genomic architecture and enigmatic phylogenetic distribution, spanning dinoflagellates and multiple cyanobacterial genera. Using filtration and extraction techniques to target the desired cyanobacteria from nonaxenic culture, coupled with a combination of short- and long-read sequencing, we generated a reference-quality hybrid genome assembly for Heteroscytonema crispum UTEX LB 1556, a freshwater, paralytic shellfish toxin-producing cyanobacterium thought to have the largest known genome in its phylum. We report a complete, novel biosynthetic gene cluster for the paralytic shellfish toxin saxitoxin. Leveraging this biosynthetic gene cluster, we find support for the hypothesis that paralytic shellfish toxin production has appeared in divergent Cyanobacteria lineages through widespread and repeated horizontal gene transfer. This work demonstrates the utility of long-read sequencing and metagenomic assembly toward advancing our understanding of paralytic shellfish toxin biosynthetic gene cluster diversity and suggests a mechanism for the origin of paralytic shellfish toxin biosynthetic genes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Geographic patterns of genomic diversity and structure in the C4 grass Panicum hallii across its natural distribution

Abstract Geographic patterns of within-species genomic diversity are shaped by evolutionary processes, life history and historical and contemporary factors. New genomic approaches can be used to infer the influence of such factors on the current distribution of infraspecific lineages. In this study, we evaluated the genomic and morphological diversity as well as the genetic structure of the C4 grass Panicum hallii across its complex natural distribution in North America. We sampled extensively across the natural range of P. hallii in Mexico and the USA to generate double-digestion restriction-associated DNA (ddRAD) sequence data for 423 individuals from 118 localities. We used these individuals to study the divergence between the two varieties of P. hallii, P. hallii var. filipes and P. hallii var. hallii as well as the genetic diversity and structure within these groups. We also examined the possibility of admixture in the geographically sympatric zone shared by both varieties, and assessed distribution shifts related with past climatic fluctuations. There is strong genetic and morphological divergence between the varieties and consistent genetic structure defining seven genetic clusters that follow major ecoregions across the range. South Texas constitutes a hotspot of genetic diversity with the co-occurrence of all genetic clusters and admixture between the two varieties. It is likely a recolonization and convergence point of populations that previously diverged in isolation during fragmentation events following glaciation periods.

54 ENVIRONMENTAL SCIENCES↗