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

EvoNet: A phylogenomic and systems biology approach to identify genes underlying plant survival in marginal, low‐N soils

The DOE‐BER “EvoNet” project investigates the genetic and molecular basis of plant resilience in extreme environments. We do this by identifying key genes that enable “extreme survivor” species to thrive in the nitrogen-poor soils of Chile’s hyper-arid Atacama Desert. Our collections focus on 32 Atacama extremophile species, including seven grass species with potential biofuel applications. To identify genes-of-importance to survival we compared genomic and transcriptomic profiles of extremophile species that thrive in the Atacama to those of closely related “sister” species from nitrogen-rich arid and mesic regions of California. Deep RNA sequencing and de novo transcriptome assembly across these triplet species sets supported a phylogenomic framework for identifying positively selected genes associated with adaptive divergence. Our integrative analysis combined ecological and environmental data, metagenomics, evolutionary and systems biology, and metabolomics. This enabled us to create an unprecedented framework for systematically understanding how non-model plants have adapted to survive in extreme conditions. Our resulting database of positively selected ortholog groups in the extremophile plants offers promising targets for engineering crop and biofuel species with enhanced resilience to drought and extreme weather. Additionally, our newest dataset explores and exploits a complementary metabolomic approach. This new aspect provides innovative strategies to manipulate plant cell metabolism, further supporting efforts to improve agricultural productivity in the face of extreme climates. Importantly, our combined evolutionary- and metabolomic-based strategies focused on convergent patterns of adaptation, providing a genetic and metabolomic toolkit for improving crop and biofuel resilience across diverse plant species. Finally, our novel exploration of ecological and evolutionary dynamics delivered to the community a phylogenomic computational pipeline called “PhyloGeneious.” Our continued adaptations of this pipeline are publicly available to expedite evolutionary genomic research for future scientific discoveries. In total, our DOE-BER has provided genomic, metabolomic, and computational strategies to understand how extremophile plants provide evolutionary and physiological targets for improving agricultural and biofuel production.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Genetic stability of foot-and-mouth disease virus during long-term infections in natural hosts

Foot-and-mouth disease (FMD) is a severe infection caused by a picornavirus that affects livestock and wildlife. Persistence in ruminants is a well-documented feature of Foot-andmouth disease virus (FMDV) pathogenesis and a major concern for disease control. Persistently infected animals harbor virus for extended periods, providing a unique opportunity to study within-host virus evolution. This study investigated the genetic dynamics of FMDV during persistent infections of naturally infected Asian buffalo. Using next-generation sequencing (NGS) we obtained 21 near complete FMDV genome sequences from 12 sub-clinically infected buffalo over a period of one year. Four animals yielded only one virus isolate and one yielded two isolates of different serotype suggesting a serial infection. Seven persistently infected animals yielded more than one virus of the same serotype showing a longterm intra-host viral genetic divergence at the consensus level of less than 2.5%. Quasi-species analysis showed few nucleotide variants and non-synonymous substitutions of progeny virus despite intra-host persistence of up to 152 days. Phylogenetic analyses of serotype Asia-1 VP1 sequences clustered all viruses from persistent animals with Group VII viruses circulating in Pakistan in 2011, but distinct from those circulating on 2008–2009. Furthermore, signature amino acid (aa) substitutions were found in the antigenically relevant VP1 of persistent viruses compared with viruses from 2008–2009. Intra-host purifying selective pressure was observed, with few codons in structural proteins undergoing positive selection. However, FMD persistent viruses did not show a clear pattern of antigenic selection. Our findings provide insight into the evolutionary dynamics of FMDV populations within naturally occurring subclinical and persistent infections that may have implications to vaccination strategies in the region.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Genome evolution in the allotetraploid frog Xenopus laevis

To explore the origins and consequences of tetraploidy in the African clawed frog, we sequenced the Xenopus laevis genome and compared it to the related diploid X. tropicalis genome. We characterize the allotetraploid origin of X. laevis by partitioning its genome into two homoeologous subgenomes, marked by distinct families of 'fossil' transposable elements. On the basis of the activity of these elements and the age of hundreds of unitary pseudogenes, we estimate that the two diploid progenitor species diverged around 34 million years ago (Ma) and combined to form an allotetraploid around 17-18 Ma. More than 56% of all genes were retained in two homoeologous copies. Protein function, gene expression, and the amount of conserved flanking sequence all correlate with retention rates. The subgenomes have evolved asymmetrically, with one chromosome set more often preserving the ancestral state and the other experiencing more gene loss, deletion, rearrangement, and reduced gene expression.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Structural basis of divergent substrate recognition and inhibition of human neurolysin

A zinc metallopeptidase neurolysin (Nln) processes diverse bioactive peptides to regulate signaling in the mammalian nervous system. To understand how Nln interacts with various peptides with dissimilar sequences, we determined crystal structures of Nln in complex with diverse peptides including dynorphins, angiotensin, neurotensin, and bradykinin. The structures show that Nln binds these peptides in a large dumbbell-shaped interior cavity constricted at the active site, making minimal structural changes to accommodate different peptide sequences. The structures also show that Nln readily binds similar peptides with distinct registers, which can determine whether the peptide serves as a substrate or a competitive inhibitor. We analyzed the activities and binding of Nln toward various forms of dynorphin A peptides, which highlights the promiscuous nature of peptide binding and shows how dynorphin A (1–13) potently inhibits the Nln activity while dynorphin A (1–8) is efficiently cleaved. Our work provides insights into the broad substrate specificity of Nln and may aid in the future design of small molecule modulators for Nln.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Expanding the substrate scope of pyrrolysyl-transfer RNA synthetase enzymes to include non-$α$-amino acids in vitro and in vivo

The absence of orthogonal aminoacyl-transfer RNA (tRNA) synthetases that accept non-L-α-amino acids is a primary bottleneck hindering the in vivo translation of sequence-defined hetero-oligomers and biomaterials. Here we report that pyrrolysyl-tRNA synthetase (PylRS) and certain PylRS variants accept α-hydroxy, α-thio and N-formyl-L-α-amino acids, as well as α-carboxy acid monomers that are precursors to polyketide natural products. These monomers are accommodated and accepted by the translation apparatus in vitro; those with reactive nucleophiles are incorporated into proteins in vivo. High-resolution structural analysis of the complex formed between one PylRS enzyme and a m-substituted 2-benzylmalonic acid derivative revealed an active site that discriminates prochiral carboxylates and accommodates the large size and distinct electrostatics of an α-carboxy substituent. This work emphasizes the potential of PylRS-derived enzymes for acylating tRNA with monomers whose α-substituent diverges substantially from the α-amine of proteinogenic amino acids. These enzymes or derivatives thereof could synergize with natural or evolved ribosomes and/or translation factors to generate diverse sequence-defined non-protein heteropolymers.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Targeted rRNA depletion enables efficient mRNA sequencing in diverse bacterial species and complex co-cultures

Bacterial mRNA sequencing is inefficient due to the abundance of ribosomal RNA that is challenging to deplete. While commercial kits target rRNA from common bacterial species, they are frequently inefficient when applied to divergent species, including those from environmental isolates. Similarly, other methods typically employ large probe sets that tile the entire length of rRNAs; however, such approaches are infeasible when applied to many species. Therefore, we present EMBR-seq+, which requires fewer than 10 oligonucleotides per rRNA by combining rRNA blocking primers with RNase H-mediated depletion to achieve rRNA removal efficiencies of up to 99% in diverse bacterial species. Furthermore, in more complex microbial co-cultures between Fibrobacter succinogenes strain UWB7 and anaerobic fungi, EMBR-seq+ depleted both bacterial and fungal rRNA, with a fourfold improvement in bacterial rRNA depletion compared with a commercial kit, thereby demonstrating that the method can be applied to non-model microbial mixtures. Notably, for microbes with unknown rRNA sequences, EMBR-seq+ enables rapid iterations in probe design without requiring to start experiments from total RNA. Finally, efficient depletion of rRNA enabled systematic quantification of the reprogramming of the bacterial transcriptome when cultured in the presence of the anaerobic fungi Anaeromyces robustus or Caecomyces churrovis. We observed that F. succinogenes strain UWB7 downregulated several lignocellulose-degrading carbohydrate-active enzymes in the presence of anaerobic gut fungi, suggesting close interactions between two cellulolytic species that specialize in different aspects of biomass breakdown. Thus, EMBR-seq+ enables efficient, cost-effective, and rapid quantification of the transcriptome to gain insights into non-model microbial systems.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Pests, diseases, and aridity have shaped the genome of Corymbia citriodora

Corymbia citriodora is a member of the predominantly Southern Hemisphere Myrtaceae family, which includes the eucalypts (Eucalyptus, Corymbia and Angophora; ~800 species). Corymbia is grown for timber, pulp and paper, and essential oils in Australia, South Africa, Asia, and Brazil, maintaining a high-growth rate under marginal conditions due to drought, poor-quality soil, and biotic stresses. To dissect the genetic basis of these desirable traits, we sequenced and assembled the 408 Mb genome of Corymbia citriodora, anchored into eleven chromosomes. Comparative analysis with Eucalyptus grandis reveals high synteny, although the two diverged approximately 60 million years ago and have different genome sizes (408 vs 641 Mb), with few large intra-chromosomal rearrangements. C. citriodora shares an ancient whole-genome duplication event with E. grandis but has undergone tandem gene family expansions related to terpene biosynthesis, innate pathogen resistance, and leaf wax formation, enabling their successful adaptation to biotic/abiotic stresses and arid conditions of the Australian continent.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

Identification and Characterization of the Major Porin of Desulfovibrio vulgaris Hildenborough

ABSTRACT Due in large part to their ability to facilitate the diffusion of a diverse range of solutes across the outer membrane (OM) of Gram-negative bacteria, the porins represent one of the most prominent and important bacterial membrane protein superfamilies. Notably, for the Gram-negative bacterium Desulfovibrio vulgaris Hildenborough, a model organism for studies of sulfate-reducing bacteria, no genes for porins have been identified or proposed in its annotated genome. Results from initial biochemical studies suggested that the product of the DVU0799 gene, which is one of the most abundant proteins of the D. vulgaris Hildenborough OM and purified as a homotrimeric complex, was a strong porin candidate. To investigate this possibility, this protein was further characterized biochemically and biophysically. Structural analyses via electron microscopy of negatively stained protein identified trimeric particles with stain-filled depressions and structural modeling suggested a β-barrel structure for the monomer, motifs common among the known porins. Functional studies were performed in which crude OM preparations or purified DVU0799 was reconstituted into proteoliposomes and the proteoliposomes were examined for permeability against a series of test solutes. The results obtained establish DVU0799 to be a pore-forming protein with permeability properties similar to those observed for classical bacterial porins, such as those of Escherichia coli . Taken together, these findings identify this highly abundant OM protein to be the major porin of D. vulgaris Hildenborough. Classification of DVU0799 in this model organism expands the database of functionally characterized porins and may also extend the range over which sequence analysis strategies can be used to identify porins in other bacterial genomes. IMPORTANCE Porins are membrane proteins that form transmembrane pores for the passive transport of small molecules across the outer membranes of Gram-negative bacteria. The present study identified and characterized the major porin of the model sulfate-reducing bacterium Desulfovibrio vulgaris Hildenborough, observing its preference for anionic sugars over neutral ones. Its predicted architecture appears to be novel for a classical porin, as its core β-barrel structure is of a type typically found in solute-specific channels. Broader use of the methods employed here, such as assays for channel permeability and electron microscopy of purified samples, is expected to help expand the database of confirmed porin sequences and improve the range over which sequence analysis-based strategies can be used to identify porins in other Gram-negative bacteria. Functional characterization of these critical gatekeeping proteins from divergent Desulfovibrio species should offer an improved understanding of the physiological features that determine their habitat range and supporting activities.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Anatomy of a Summertime Convective Event over the Arabian Region

This study investigates the structure and evolution of a summertime convective event that occurred on 14 July 2015 over the Arabian region. We use the WRF Model with 1-km horizontal grid spacing and test three PBL parameterizations: the Mellor–Yamada–Nakanishi–Niino (MYNN) scheme; the Asymmetrical Convective Model, version 2, (ACM2) scheme; and the quasi-normal scale-elimination (QNSE) scheme. Convection initiates near the Al Hajar Mountains of northern Oman at around 1100 local time (LT; 0700 UTC) and propagates northwestward. A nonorographic convective band along the west coast of the United Arab Emirates (UAE) develops after 1500 LT as a result of the convergence of cold pools with the sea breeze from the Arabian Gulf. The model simulation employing the QNSE scheme simulates the convection initiation and propagation well. Although the MYNN and ACM2 simulations show convective initiation near the Al Hajar Mountains, they fail to simulate the development of the convective band along the UAE west coast. The MYNN run simulates colder near-surface temperatures and a weaker sea breeze, whereas the ACM2 run simulates a stronger sea breeze but a drier lower troposphere. Sensitivity simulations using horizontal grid spacings of 9 and 3 km show that lower-resolution runs develop broader convective structures and weaker cold pools and horizontal wind divergence, affecting the development of convection along the west coast of the UAE. The 1-km run using the QNSE PBL scheme realistically captures the sequence of events that leads to the moist convection over the UAE and adjacent mountains.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Diversity, taxonomy, and evolution of archaeal viruses of the class Caudoviricetes

The archaeal tailed viruses (arTV), evolutionarily related to tailed double-stranded DNA (dsDNA) bacteriophages of the class Caudoviricetes , represent the most common isolates infecting halophilic archaea. Only a handful of these viruses have been genomically characterized, limiting our appreciation of their ecological impacts and evolution. Here, we present 37 new genomes of haloarchaeal tailed virus isolates, more than doubling the current number of sequenced arTVs. Analysis of all 63 available complete genomes of arTVs, which we propose to classify into 14 new families and 3 orders, suggests ancient divergence of archaeal and bacterial tailed viruses and points to an extensive sharing of genes involved in DNA metabolism and counterdefense mechanisms, illuminating common strategies of virus–host interactions with tailed bacteriophages. Coupling of the comparative genomics with the host range analysis on a broad panel of haloarchaeal species uncovered 4 distinct groups of viral tail fiber adhesins controlling the host range expansion. The survey of metagenomes using viral hallmark genes suggests that the global architecture of the arTV community is shaped through recurrent transfers between different biomes, including hypersaline, marine, and anoxic environments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Turicibacter bilis sp. nov., a novel bacterium isolated from the chicken eggshell and swine ileum

Three novel, anaerobic, Gram-positive bacteria were isolated from the eggshell of two separate white leghorn chicken flocks and the ileum of a healthy pig, and designated MMM721 T , ISU324 and PIG517 respectively. Cells were pleomorphic and capable of forming long chains of rods or coccoid clusters. Phylogenetic analysis of the 16S rRNA gene sequences identified these strains to be within the genus Turicibacter, of which only one species, Turicibacter sanguinis, has been formally described. However, whole genome sequencing of novel isolates returned a digital DNA–DNA hybridization value of 22.5% and average nucleotide identity (ANI) values of 76.4% (ANIb) and 86.0% (ANIm), indicating divergence between the type strain MMM721 T and T. sanguinis, suggesting the strains represented a novel species. The major fatty acid methyl esters of strain MMM721 T were C 16:0 , C 18:1 ω7c and C 18:0 . The strains mainly produced the volatile fatty acid lactate, along with smaller amounts of acetate and butyrate. Together, these data indicate that MMM721 T , along with ISU324 and PIG517, represent a novel species within the genus Turicibacter. We propose the name Turicibacter bilis sp. nov. for the species. The type strain is MMM721 T (=ATCC TSD-238 T =CCUG 74757 T ).

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

Protein abundances can distinguish between naturally-occurring and laboratory strains of Yersinia pestis , the causative agent of plague

Adaptive processes in bacterial species can occur rapidly in laboratory culture, leading to genetic divergence between naturally occurring and laboratory-adapted strains. Differentiating wild and closely-related laboratory strains is clearly important for biodefense and bioforensics; however, DNA sequence data alone has thus far not provided a clear signature, perhaps due to lack of understanding of how diverse genome changes lead to adapted phenotypes. Protein abundance profiles from mass spectrometry-based proteomics analyses are a molecular measure of phenotype. Proteomics data contains sufficient information that powerful statistical methods can uncover signatures that distinguish wild strains of Yersinia pestis from laboratory-adapted strains.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Genome Sequencing of Cladobotryum protrusum Provides Insights into the Evolution and Pathogenic Mechanisms of the Cobweb Disease Pathogen on Cultivated Mushroom

Cladobotryum protrusum is one of the mycoparasites that cause cobweb disease on cultivated edible mushrooms. However, the molecular mechanisms of evolution and pathogenesis of C. protrusum on mushrooms are largely unknown. Here, we report a high-quality genome sequence of C. protrusum using the single-molecule, real-time sequencing platform of PacBio and perform a comparative analysis with closely related fungi in the family Hypocreaceae. The C. protrusum genome, the first complete genome to be sequenced in the genus Cladobotryum, is 39.09 Mb long, with an N50 of 4.97 Mb, encoding 11,003 proteins. The phylogenomic analysis confirmed its inclusion in Hypocreaceae, with its evolutionary divergence time estimated to be ~170.1 million years ago. The genome encodes a large and diverse set of genes involved in secreted peptidases, carbohydrate-active enzymes, cytochrome P450 enzymes, pathogen–host interactions, mycotoxins, and pigments. Moreover, C. protrusum harbors arrays of genes with the potential to produce bioactive secondary metabolites and stress response-related proteins that are significant for adaptation to hostile environments. Knowledge of the genome will foster a better understanding of the biology of C. protrusum and mycoparasitism in general, as well as help with the development of effective disease control strategies to minimize economic losses from cobweb disease in cultivated edible mushrooms.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Evolutionary characteristics of intergenic transcribed regions indicate rare novel genes and widespread noisy transcription in the Poaceae

Abstract Extensive transcriptional activity occurring in intergenic regions of genomes has raised the question whether intergenic transcription represents the activity of novel genes or noisy expression. To address this, we evaluated cross-species and post-duplication sequence and expression conservation of intergenic transcribed regions (ITRs) in four Poaceae species. Among 43,301 ITRs across the four species, 34,460 (80%) are species-specific. ITRs found across species tend to be more divergent in expression and have more recent duplicates compared to annotated genes. To assess if ITRs are functional (under selection), machine learning models were established in Oryza sativa (rice) that could accurately distinguish between phenotype genes and pseudogenes (area under curve-receiver operating characteristic = 0.94). Based on the models, 584 (8%) and 4391 (61%) rice ITRs are classified as likely functional and nonfunctional with high confidence, respectively. ITRs with conserved expression and ancient retained duplicates, features that were not part of the model, are frequently classified as likely-functional, suggesting these characteristics could serve as pragmatic rules of thumb for identifying candidate sequences likely to be under selection. This study also provides a framework to identify novel genes using comparative transcriptomic data to improve genome annotation that is fundamental for connecting genotype to phenotype in crop and model systems.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Genome-scale model development and genomic sequencing of the oleaginous clade Lipomyces

The Lipomyces clade contains oleaginous yeast species with advantageous metabolic features for biochemical and biofuel production. Limited knowledge about the metabolic networks of the species and limited tools for genetic engineering have led to a relatively small amount of research on the microbes. Here, a genome-scale metabolic model (GSM) of Lipomyces starkeyi NRRL Y-11557 was built using orthologous protein mappings to model yeast species. Phenotypic growth assays were used to validate the GSM (66% accuracy) and indicated that NRRL Y-11557 utilized diverse carbohydrates but had more limited catabolism of organic acids. The final GSM contained 2,193 reactions, 1,909 metabolites, and 996 genes and was thus named iLst996. The model contained 96 of the annotated carbohydrate-active enzymes. iLst996 predicted a flux distribution in line with oleaginous yeast measurements and was utilized to predict theoretical lipid yields. Twenty-five other yeasts in the Lipomyces clade were then genome sequenced and annotated. Sixteen of the Lipomyces species had orthologs for more than 97% of the iLst996 genes, demonstrating the usefulness of iLst996 as a broad GSM for Lipomyces metabolism. Pathways that diverged from iLst996 mainly revolved around alternate carbon metabolism, with ortholog groups excluding NRRL Y-11557 annotated to be involved in transport, glycerolipid, and starch metabolism, among others. Overall, this study provides a useful modeling tool and data for analyzing and understanding Lipomyces species metabolism and will assist further engineering efforts in Lipomyces .

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

The Kalanchoe genome provides insights into convergent evolution and building blocks of crassulacean acid metabolism

Crassulacean acid metabolism (CAM) is a specialized photosynthetic adaptation to arid environments, found predominantly in diverse eudicotyledonous (eudicot) and monocotyledonous (monocot) lineages that diverged approximately 135 million years ago. To test whether convergent evolution underpins the independent emergences of CAM, we present de novo genome assembly and gene expression data for Kalanchoë fedtschenkoi, an obligate CAM species that was shown by multigene phylogenetic analysis to represent one of the earliest-diverging lineages of core eudicots. Our combined analysis of K. fedtschenkoi and two monocot CAM species (Ananas comosus and Phalaenopsis equestris) identified signatures of convergence in protein sequence and in the diel re-scheduling of genes involved in metabolism and signaling. Lastly, our results provide significant insight into CAM evolution, facilitating CAM-into-C 3 engineering for enhancing drought tolerance in crops.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Genome resources for three modern cotton lines guide future breeding efforts

Cotton ( Gossypium hirsutum L.) is the key renewable fibre crop worldwide, yet its yield and fibre quality show high variability due to genotype-specific traits and complex interactions among cultivars, management practices and environmental factors. Modern breeding practices may limit future yield gains due to a narrow founding gene pool. Precision breeding and biotechnological approaches offer potential solutions, contingent on accurate cultivar-specific data. Here we address this need by generating high-quality reference genomes for three modern cotton cultivars (‘UGA230’, ‘UA48’ and ‘CSX8308’) and updating the ‘TM-1’ cotton genetic standard reference. Despite hypothesized genetic uniformity, considerable sequence and structural variation was observed among the four genomes, which overlap with ancient and ongoing genomic introgressions from ‘Pima’ cotton, gene regulatory mechanisms and phenotypic trait divergence. Differentially expressed genes across fibre development correlate with fibre production, potentially contributing to the distinctive fibre quality traits observed in modern cotton cultivars. These genomes and comparative analyses provide a valuable foundation for future genetic endeavours to enhance global cotton yield and sustainability.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Generation of a chromosome-scale genome assembly of the insect-repellent terpenoid-producing Lamiaceae species, Callicarpa americana

Plants exhibit wide chemical diversity due to the production of specialized metabolites that function as pollinator attractants, defensive compounds, and signaling molecules. Lamiaceae (mints) are known for their chemodiversity and have been cultivated for use as culinary herbs, as well as sources of insect repellents, health-promoting compounds, and fragrance. We report the chromosome-scale genome assembly of Callicarpa americana L. (American beautyberry), a species within the early-diverging Callicarpoideae clade of Lamiaceae, known for its metallic purple fruits and use as an insect repellent due to its production of terpenoids. Using long-read sequencing and Hi-C scaffolding, we generated a 506.1-Mb assembly spanning 17 pseudomolecules with N50 contig and N50 scaffold sizes of 7.5 and 29.0 Mb, respectively. In all, 32,164 genes were annotated, including 53 candidate terpene synthases and 47 putative clusters of specialized metabolite biosynthetic pathways. Our analyses revealed 3 putative whole-genome duplication events, which, together with local tandem duplications, contributed to gene family expansion of terpene synthases. Kolavenyl diphosphate is a gateway to many of the bioactive terpenoids in C. americana; experimental validation confirmed that CamTPS2 encodes kolavenyl diphosphate synthase. Syntenic analyses with Tectona grandis L. f. (teak), a member of the Tectonoideae clade of Lamiaceae known for exceptionally strong wood resistant to insects, revealed 963 collinear blocks and 21,297 C. americana syntelogs. Access to the C. americana genome provides a road map for rapid discovery of genes encoding plant-derived agrichemicals and a key resource for understanding the evolution of chemical diversity in Lamiaceae.

59 BASIC BIOLOGICAL SCIENCES↗