The genomic and enzymatic basis for iridoid biosynthesis in cat thyme (Teucrium marum).

Iridoids are non-canonical monoterpenoids produced by both insects and plants. An example is the cat-attracting and insect-repelling volatile iridoid nepetalactone, produced by Nepeta sp. (catmint) and aphids. Recently, both nepetalactone biosynthetic pathways were elucidated, showing a remarkable convergent evolution. The iridoid, dolichodial, produced by Teucrium marum (cat thyme) and multiple insect species, has highly similar properties to nepetalactone but its biosynthetic origin remains unknown. We set out to determine the genomic, enzymatic, and evolutionary basis of iridoid biosynthesis in T. marum. First, we generated a de novo chromosome-scale genome assembly for T. marum using Oxford Nanopore Technologies long reads and proximity-by-ligation Hi-C reads. The 610.3 Mb assembly spans 15 pseudomolecules with a 32.9 Mb N50 scaffold size. This enabled identification of iridoid biosynthetic genes, whose roles were verified via activity assays. Phylogenomic analysis revealed that the evolutionary history of T. marum iridoid synthase, the iridoid scaffold-forming enzyme, is not orthologous to typical iridoid synthases but is derived from its conserved paralog. We discovered an enzymatic route from nepetalactol to diverse iridoids through the coupled activity of an iridoid oxidase cytochrome P450 and acetyltransferases, via an inferred acylated intermediate. This work provides a genomic resource for specialized metabolite research in mints and demonstration of the role of acetylation in T. marum iridoid diversity. This work will enable future biocatalytic or biosynthetic production of potent insect repellents, as well as comparative studies into iridoid biosynthesis in insects.

Single-cell multi-omics in the medicinal plant Catharanthus roseus.

Advances in omics technologies now permit the generation of highly contiguous genome assemblies, detection of transcripts and metabolites at the level of single cells and high-resolution determination of gene regulatory features. Here, using a complementary, multi-omics approach, we interrogated the monoterpene indole alkaloid (MIA) biosynthetic pathway in Catharanthus roseus, a source of leading anticancer drugs. We identified clusters of genes involved in MIA biosynthesis on the eight C. roseus chromosomes and extensive gene duplication of MIA pathway genes. Clustering was not limited to the linear genome, and through chromatin interaction data, MIA pathway genes were present within the same topologically associated domain, permitting the identification of a secologanin transporter. Single-cell RNA-sequencing revealed sequential cell-type-specific partitioning of the leaf MIA biosynthetic pathway that, when coupled with a single-cell metabolomics approach, permitted the identification of a reductase that yields the bis-indole alkaloid anhydrovinblastine. We also revealed cell-type-specific expression in the root MIA pathway.

Eleven biosynthetic genes explain the majority of natural variation in carotenoid levels in maize grain.

Vitamin A deficiency remains prevalent in parts of Asia, Latin America, and sub-Saharan Africa where maize (Zea mays) is a food staple. Extensive natural variation exists for carotenoids in maize grain. Here, to understand its genetic basis, we conducted a joint linkage and genome-wide association study of the US maize nested association mapping panel. Eleven of the 44 detected quantitative trait loci (QTL) were resolved to individual genes. Six of these were correlated expression and effect QTL (ceeQTL), showing strong correlations between RNA-seq expression abundances and QTL allelic effect estimates across six stages of grain development. These six ceeQTL also had the largest percentage of phenotypic variance explained, and in major part comprised the three to five loci capturing the bulk of genetic variation for each trait. Most of these ceeQTL had strongly correlated QTL allelic effect estimates across multiple traits. These findings provide an in-depth genome-level understanding of the genetic and molecular control of carotenoids in plants. In addition, these findings provide a roadmap to accelerate breeding for provitamin A and other priority carotenoid traits in maize grain that should be readily extendable to other cereals.

An updated gene atlas for maize reveals organ-specific and stress-induced genes.

Maize (Zea mays L.), a model species for genetic studies, is one of the two most important crop species worldwide. The genome sequence of the reference genotype, B73, representative of the stiff stalk heterotic group was recently updated (AGPv4) using long-read sequencing and optical mapping technology. To facilitate the use of AGPv4 and to enable functional genomic studies and association of genotype with phenotype, we determined expression abundances for replicated mRNA-sequencing datasets from 79 tissues and five abiotic/biotic stress treatments revealing 36 207 expressed genes. Characterization of the B73 transcriptome across six organs revealed 4154 organ-specific and 7704 differentially expressed (DE) genes following stress treatment. Gene co-expression network analyses revealed 12 modules associated with distinct biological processes containing 13 590 genes providing a resource for further association of gene function based on co-expression patterns. Presence-absence variants (PAVs) previously identified using whole genome resequencing data from 61 additional inbred lines were enriched in organ-specific and stress-induced DE genes suggesting that PAVs may function in phenological variation and adaptation to environment. Relative to core genes conserved across the 62 profiled inbreds, PAVs have lower expression abundances which are correlated with their frequency of dispersion across inbreds and on average have significantly fewer co-expression network connections suggesting that a subset of PAVs may be on an evolutionary path to pseudogenization. To facilitate use by the community, we developed the Maize Genomics Resource website (maize.plantbiology.msu.edu) for viewing and data-mining these resources and deployed two new views on the maize electronic Fluorescent Pictograph Browser (bar.utoronto.ca/efp_maize).

Quantitative trait loci and differential gene expression analyses reveal the genetic basis for negatively associated ß-carotene and starch content in hexaploid sweetpotato [Ipomoea batatas (L.) Lam.].

KEY MESSAGE: ß-Carotene content in sweetpotato is associated with the Orange and phytoene synthase genes; due to physical linkage of phytoene synthase with sucrose synthase, ß-carotene and starch content are negatively correlated. In populations depending on sweetpotato for food security, starch is an important source of calories, while ß-carotene is an important source of provitamin A. The negative association between the two traits contributes to the low nutritional quality of sweetpotato consumed, especially in sub-Saharan Africa. Using a biparental mapping population of 315 F1 progeny generated from a cross between an orange-fleshed and a non-orange-fleshed sweetpotato variety, we identified two major quantitative trait loci (QTL) on linkage group (LG) three (LG3) and twelve (LG12) affecting starch, ß-carotene, and their correlated traits, dry matter and flesh color. Analysis of parental haplotypes indicated that these two regions acted pleiotropically to reduce starch content and increase ß-carotene in genotypes carrying the orange-fleshed parental haplotype at the LG3 locus. Phytoene synthase and sucrose synthase, the rate-limiting and linked genes located within the QTL on LG3 involved in the carotenoid and starch biosynthesis, respectively, were differentially expressed in Beauregard versus Tanzania storage roots. The Orange gene, the molecular switch for chromoplast biogenesis, located within the QTL on LG12 while not differentially expressed was expressed in developing roots of the parental genotypes. We conclude that these two QTL regions act together in a cis and trans manner to inhibit starch biosynthesis in amyloplasts and enhance chromoplast biogenesis, carotenoid biosynthesis, and accumulation in orange-fleshed sweetpotato. Understanding the genetic basis of this negative association between starch and ß-carotene will inform future sweetpotato breeding strategies targeting sweetpotato for food and nutritional security.

The unique epigenetic features of Pack-MULEs and their impact on chromosomal base composition and expression spectrum.

Acquisition and rearrangement of host genes by transposable elements (TEs) is an important mechanism to increase gene diversity as exemplified by the ∼3000 Pack-Mutator-like TEs in the rice genome which have acquired gene sequences (Pack-MULEs), yet remain enigmatic. To identify signatures of functioning Pack-MULEs and Pack-MULE evolution, we generated transcriptome, translatome, and epigenome datasets and compared Pack-MULEs to genes and other TE families. Approximately 40% of Pack-MULEs were transcribed with 9% having translation evidence, clearly distinguishing them from other TEs. Pack-MULEs exhibited a unique expression profile associated with specificity in reproductive tissues that may be associated with seed traits. Expressed Pack-MULEs resemble regular protein-coding genes as exhibited by a low level of DNA methylation, association with active histone marks and DNase I hypersensitive sites, and an absence of repressive histone marks, suggesting that a substantial fraction of Pack-MULEs are potentially functional in vivo. Interestingly, the expression capacity of Pack-MULEs is independent of the local genomic environment, and the insertion and expression of Pack-MULEs may have altered the local chromosomal expression pattern as well as counteracted the impact of recombination on chromosomal base composition, which has profound consequences on the evolution of chromosome structure.

Genome diversity of tuber-bearing Solanum uncovers complex evolutionary history and targets of domestication in the cultivated potato.

Cultivated potatoes (Solanum tuberosum L.), domesticated from wild Solanum species native to the Andes of southern Peru, possess a diverse gene pool representing more than 100 tuber-bearing relatives (Solanum section Petota). A diversity panel of wild species, landraces, and cultivars was sequenced to assess genetic variation within tuber-bearing Solanum and the impact of domestication on genome diversity and identify key loci selected for cultivation in North and South America. Sequence diversity of diploid and tetraploid Stuberosum exceeded any crop resequencing study to date, in part due to expanded wild introgressions following polyploidy that captured alleles outside of their geographic origin. We identified 2,622 genes as under selection, with only 14-16% shared by North American and Andean cultivars, showing that a limited gene set drove early improvement of cultivated potato, while adaptation of upland (Stuberosum group Andigena) and lowland (S. tuberosum groups Chilotanum and Tuberosum) populations targeted distinct loci. Signatures of selection were uncovered in genes controlling carbohydrate metabolism, glycoalkaloid biosynthesis, the shikimate pathway, the cell cycle, and circadian rhythm. Reduced sexual fertility that accompanied the shift to asexual reproduction in cultivars was reflected by signatures of selection in genes regulating pollen development/gametogenesis. Exploration of haplotype diversity at potato's maturity locus (StCDF1) revealed introgression of truncated alleles from wild species, particularly Smicrodontum in long-day-adapted cultivars. This study uncovers a historic role of wild Solanum species in the diversification of long-day-adapted tetraploid potatoes, showing that extant natural populations represent an essential source of untapped adaptive potential.

Genome sequence of M6, a diploid inbred clone of the high-glycoalkaloid-producing tuber-bearing potato species Solanum chacoense, reveals residual heterozygosity.

Cultivated potato (Solanum tuberosum L.) is a highly heterozygous autotetraploid that presents challenges in genome analyses and breeding. Wild potato species serve as a resource for the introgression of important agronomic traits into cultivated potato. One key species is Solanum chacoense and the diploid, inbred clone M6, which is self-compatible and has desirable tuber market quality and disease resistance traits. Sequencing and assembly of the genome of the M6 clone of S. chacoense generated an assembly of 825 767 562 bp in 8260 scaffolds with an N50 scaffold size of 713 602 bp. Pseudomolecule construction anchored 508 Mb of the genome assembly into 12 chromosomes. Genome annotation yielded 49 124 high-confidence gene models representing 37 740 genes. Comparative analyses of the M6 genome with six other Solanaceae species revealed a core set of 158 367 Solanaceae genes and 1897 genes unique to three potato species. Analysis of single nucleotide polymorphisms across the M6 genome revealed enhanced residual heterozygosity on chromosomes 4, 8 and 9 relative to the other chromosomes. Access to the M6 genome provides a resource for identification of key genes for important agronomic traits and aids in genome-enabled development of inbred diploid potatoes with the potential to accelerate potato breeding.

Gene Discovery in Gelsemium Highlights Conserved Gene Clusters in Monoterpene Indole Alkaloid Biosynthesis.

Genome mining is a routine technique in microbes for discovering biosynthetic pathways. In plants, however, genomic information is not commonly used to identify novel biosynthesis genes. Here, we present the genome of the medicinal plant and oxindole monoterpene indole alkaloid (MIA) producer Gelsemium sempervirens (Gelsemiaceae). A gene cluster from Catharanthus roseus, which is utilized at least six enzymatic steps downstream from the last common intermediate shared between the two plant alkaloid types, is found in G. sempervirens, although the corresponding enzymes act on entirely different substrates. This study provides insights into the common genomic context of MIA pathways and is an important milestone in the further elucidation of the Gelsemium oxindole alkaloid pathway.

Genome Reduction Uncovers a Large Dispensable Genome and Adaptive Role for Copy Number Variation in Asexually Propagated Solanum tuberosum.

Clonally reproducing plants have the potential to bear a significantly greater mutational load than sexually reproducing species. To investigate this possibility, we examined the breadth of genome-wide structural variation in a panel of monoploid/doubled monoploid clones generated from native populations of diploid potato (Solanum tuberosum), a highly heterozygous asexually propagated plant. As rare instances of purely homozygous clones, they provided an ideal set for determining the degree of structural variation tolerated by this species and deriving its minimal gene complement. Extensive copy number variation (CNV) was uncovered, impacting 219.8 Mb (30.2%) of the potato genome with nearly 30% of genes subject to at least partial duplication or deletion, revealing the highly heterogeneous nature of the potato genome. Dispensable genes (>7000) were associated with limited transcription and/or a recent evolutionary history, with lower deletion frequency observed in genes conserved across angiosperms. Association of CNV with plant adaptation was highlighted by enrichment in gene clusters encoding functions for environmental stress response, with gene duplication playing a part in species-specific expansions of stress-related gene families. This study revealed unique impacts of CNV in a species with asexual reproductive habits and how CNV may drive adaption through evolution of key stress pathways.

A root-expressed L-phenylalanine:4-hydroxyphenylpyruvate aminotransferase is required for tropane alkaloid biosynthesis in Atropa belladonna.

The tropane alkaloids, hyoscyamine and scopolamine, are medicinal compounds that are the active components of several therapeutics. Hyoscyamine and scopolamine are synthesized in the roots of specific genera of the Solanaceae in a multistep pathway that is only partially elucidated. To facilitate greater understanding of tropane alkaloid biosynthesis, a de novo transcriptome assembly was developed for Deadly Nightshade (Atropa belladonna). Littorine is a key intermediate in hyoscyamine and scopolamine biosynthesis that is produced by the condensation of tropine and phenyllactic acid. Phenyllactic acid is derived from phenylalanine via its transamination to phenylpyruvate, and mining of the transcriptome identified a phylogenetically distinct aromatic amino acid aminotransferase (ArAT), designated Ab-ArAT4, that is coexpressed with known tropane alkaloid biosynthesis genes in the roots of A. belladonna. Silencing of Ab-ArAT4 disrupted synthesis of hyoscyamine and scopolamine through reduction of phenyllactic acid levels. Recombinant Ab-ArAT4 preferentially catalyzes the first step in phenyllactic acid synthesis, the transamination of phenylalanine to phenylpyruvate. However, rather than utilizing the typical keto-acid cosubstrates, 2-oxoglutarate, pyruvate, and oxaloacetate, Ab-ArAT4 possesses strong substrate preference and highest activity with the aromatic keto-acid, 4-hydroxyphenylpyruvate. Thus, Ab-ArAT4 operates at the interface between primary and specialized metabolism, contributing to both tropane alkaloid biosynthesis and the direct conversion of phenylalanine to tyrosine.

Characterization of the Xanthomonas translucens Complex Using Draft Genomes, Comparative Genomics, Phylogenetic Analysis, and Diagnostic LAMP Assays.

Prevalence of Xanthomonas translucens, which causes cereal leaf streak (CLS) in cereal crops and bacterial wilt in forage and turfgrass species, has increased in many regions in recent years. Because the pathogen is seedborne in economically important cereals, it is a concern for international and interstate germplasm exchange and, thus, reliable and robust protocols for its detection in seed are needed. However, historical confusion surrounding the taxonomy within the species has complicated the development of accurate and reliable diagnostic tools for X. translucens. Therefore, we sequenced genomes of 15 X. translucens strains representing six different pathovars and compared them with additional publicly available X. translucens genome sequences to obtain a genome-based phylogeny for robust classification of this species. Our results reveal three main clusters: one consisting of pv. cerealis, one consisting of pvs. undulosa and translucens, and a third consisting of pvs. arrhenatheri, graminis, phlei, and poae. Based on genomic differences, diagnostic loop-mediated isothermal amplification (LAMP) primers were developed that clearly distinguish strains that cause disease on cereals, such as pvs. undulosa, translucens, hordei, and secalis, from strains that cause disease on noncereal hosts, such as pvs. arrhenatheri, cerealis, graminis, phlei, and poae. Additional LAMP assays were developed that selectively amplify strains belonging to pvs. cerealis and poae, distinguishing them from other pathovars. These primers will be instrumental in diagnostics when implementing quarantine regulations to limit further geographic spread of X. translucens pathovars.

Genome-guided investigation of plant natural product biosynthesis.

The medicinal plant Madagascar periwinkle, Catharanthus roseus (L.) G. Don, produces hundreds of biologically active monoterpene-derived indole alkaloid (MIA) metabolites and is the sole source of the potent, expensive anti-cancer compounds vinblastine and vincristine. Access to a genome sequence would enable insights into the biochemistry, control, and evolution of genes responsible for MIA biosynthesis. However, generation of a near-complete, scaffolded genome is prohibitive to small research communities due to the expense, time, and expertise required. In this study, we generated a genome assembly for C. roseus that provides a near-comprehensive representation of the genic space that revealed the genomic context of key points within the MIA biosynthetic pathway including physically clustered genes, tandem gene duplication, expression sub-functionalization, and putative neo-functionalization. The genome sequence also facilitated high resolution co-expression analyses that revealed three distinct clusters of co-expression within the components of the MIA pathway. Coordinated biosynthesis of precursors and intermediates throughout the pathway appear to be a feature of vinblastine/vincristine biosynthesis. The C. roseus genome also revealed localization of enzyme-rich genic regions and transporters near known biosynthetic enzymes, highlighting how even a draft genome sequence can empower the study of high-value specialized metabolites.

Gene amplification of 5-enol-pyruvylshikimate-3-phosphate synthase in glyphosate-resistant Kochia scoparia.

MAIN CONCLUSION: Field-evolved resistance to the herbicide glyphosate is due to amplification of one of two EPSPS alleles, increasing transcription and protein with no splice variants or effects on other pathway genes. The widely used herbicide glyphosate inhibits the shikimate pathway enzyme 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS). Globally, the intensive use of glyphosate for weed control has selected for glyphosate resistance in 31 weed species. Populations of suspected glyphosate-resistant Kochia scoparia were collected from fields located in the US central Great Plains. Glyphosate dose response verified glyphosate resistance in nine populations. The mechanism of resistance to glyphosate was investigated using targeted sequencing, quantitative PCR, immunoblotting, and whole transcriptome de novo sequencing to characterize the sequence and expression of EPSPS. Sequence analysis showed no mutation of the EPSPS Pro106 codon in glyphosate-resistant K. scoparia, whereas EPSPS genomic copy number and transcript abundance were elevated three- to ten-fold in resistant individuals relative to susceptible individuals. Glyphosate-resistant individuals with increased relative EPSPS copy numbers had consistently lower shikimate accumulation in leaf disks treated with 100 µM glyphosate and EPSPS protein levels were higher in glyphosate-resistant individuals with increased gene copy number compared to glyphosate-susceptible individuals. RNA sequence analysis revealed seven nucleotide positions with two different expressed alleles in glyphosate-susceptible reads. However, one nucleotide at the seven positions was predominant in glyphosate-resistant sequences, suggesting that only one of two EPSPS alleles was amplified in glyphosate-resistant individuals. No alternatively spliced EPSPS transcripts were detected. Expression of five other genes in the chorismate pathway was unaffected in glyphosate-resistant individuals with increased EPSPS expression. These results indicate increased EPSPS expression is a mechanism for glyphosate resistance in these K. scoparia populations.

Targeted genotyping-by-sequencing of potato and data analysis with R/polyBreedR.

Mid-density targeted genotyping-by-sequencing (GBS) combines trait-specific markers with thousands of genomic markers at an attractive price for linkage mapping and genomic selection. A 2.5K targeted GBS assay for potato (Solanum tuberosum L.) was developed using the DArTag technology and later expanded to 4K targets. Genomic markers were selected from the potato Infinium single nucleotide polymorphism (SNP) array to maximize genome coverage and polymorphism rates. The DArTag and SNP array platforms produced equivalent dendrograms in a test set of 298 tetraploid samples, and 83% of the common markers showed good quantitative agreement, with RMSE (root mean squared error) <0.5. DArTag is suited for genomic selection candidates in the clonal evaluation trial, coupled with imputation to a higher density platform for the training population. Using the software polyBreedR, an R package for the manipulation and analysis of polyploid marker data, the RMSE for imputation by linkage analysis was 0.15 in a small half-diallel population (N = 85), which was significantly lower than the RMSE of 0.42 with the random forest method. Regarding high-value traits, the DArTag markers for resistance to potato virus Y, golden cyst nematode, and potato wart appeared to track their targets successfully, as did multi-allelic markers for maturity and tuber shape. In summary, the potato DArTag assay is a transformative and publicly available technology for potato breeding and genetics.

Advances in plant genome sequencing.

The study of plant biology in the 21st century is, and will continue to be, vastly different from that in the 20th century. One driver for this has been the use of genomics methods to reveal the genetic blueprints for not one but dozens of plant species, as well as resolving genome differences in thousands of individuals at the population level. Genomics technology has advanced substantially since publication of the first plant genome sequence, that of Arabidopsis thaliana, in 2000. Plant genomics researchers have readily embraced new algorithms, technologies and approaches to generate genome, transcriptome and epigenome datasets for model and crop species that have permitted deep inferences into plant biology. Challenges in sequencing any genome include ploidy, heterozygosity and paralogy, all which are amplified in plant genomes compared to animal genomes due to the large genome sizes, high repetitive sequence content, and rampant whole- or segmental genome duplication. The ability to generate de novo transcriptome assemblies provides an alternative approach to bypass these complex genomes and access the gene space of these recalcitrant species. The field of genomics is driven by technological improvements in sequencing platforms; however, software and algorithm development has lagged behind reductions in sequencing costs, improved throughput, and quality improvements. It is anticipated that sequencing platforms will continue to improve the length and quality of output, and that the complementary algorithms and bioinformatic software needed to handle large, repetitive genomes will improve. The future is bright for an exponential improvement in our understanding of plant biology.

Chromosome-scale assembly of the Verbenaceae species Queen's Wreath (Petrea volubilis L.).

OBJECTIVES: Petrea volubilis, a member of the Order Lamiales and the Verbenaceae family, is an important horticultural species that has been used in traditional folk medicine. To provide a genome sequence for comparative studies within the Order Lamiales that includes important families such as Lamiaceae (mints), we generated a long-read, chromosome-scale genome assembly of this species. DATA DESCRIPTION: Using a total of 45.5 Gb of Pacific Biosciences long read sequence, we generated a 480.2 Mb assembly of P. volubilis, of which, 93% is chromosome anchored. Representation of genic regions was robust with 96.6% of the Benchmarking of Universal Single Copy Orthologs present in the genome assembly. A total of 57.8% of the genome was annotated as a repetitive sequence. Using a gene annotation pipeline that included refinement of gene models using transcript evidence, 30,982 high confidence genes were annotated. Access to the P. volubilis genome will facilitate evolutionary studies in the Lamiales, a key order of Asterids that includes significant crop and medicinal plant species.

Genome Report: Genome sequence of 1S1, a transformable and highly regenerable diploid potato for use as a model for gene editing and genetic engineering.

Availability of readily transformable germplasm, as well as efficient pipelines for gene discovery are notable bottlenecks in the application of genome editing in potato. To study and introduce traits such as resistance against biotic and abiotic factors, tuber quality traits and self-fertility, model germplasm that is amenable to gene editing and regeneration is needed. Cultivated potato is a heterozygous autotetraploid and its genetic redundancy and complexity makes studying gene function challenging. Genome editing is simpler at the diploid level, with fewer allelic variants to consider. A readily transformable diploid potato would be further complemented by genomic resources that could aid in high throughput functional analysis. The heterozygous Solanum tuberosum Group Phureja clone 1S1 has a high regeneration rate, self-fertility, desirable tuber traits and is amenable to Agrobacterium-mediated transformation. We leveraged its amenability to Agrobacterium-mediated transformation to create a Cas9 constitutively expressing line for use in viral vector-based gene editing. To create a contiguous genome assembly, a homozygous doubled monoploid of 1S1 (DM1S1) was sequenced using 44 Gbp of long reads generated from Oxford Nanopore Technologies (ONT), yielding a 736 Mb assembly that encoded 31,145 protein-coding genes. The final assembly for DM1S1 represents a nearly complete genic space, shown by the presence of 99.6% of the genes in the Benchmarking Universal Single Copy Orthologs (BUSCO) set. Variant analysis with Illumina reads from 1S1 was used to deduce its alternate haplotype. These genetic and genomic resources provide a toolkit for applications of genome editing in both basic and applied research of potato.