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.

Genome-wide association identifies a missing hydrolase for tocopherol synthesis in plants.

The tocopherol biosynthetic pathway, encoded by VTE genes 1 through 6, is highly conserved in plants but most large effect quantitative trait loci for seed total tocopherols (totalT) lack VTE genes, indicating other activities are involved. A genome-wide association study of Arabidopsis seed tocopherols showed five of seven significant intervals lacked VTE genes, including the most significant, which mapped to an uncharacterized, seed-specific, envelope-localized, alpha/beta hydrolase with esterase activity, designated AtVTE7. Atvte7 null mutants decreased seed totalT 55% while a leaky allele of the maize ortholog, ZmVTE7, decreased kernel and leaf totalT 38% and 49%, respectively. Overexpressing AtVTE7 or ZmVTE7 partially or fully complemented the Atvte7 seed phenotype and increased leaf totalT by 3.6- and 6.9-fold, respectively. VTE7 has the characteristics of an esterase postulated to provide phytol from chlorophyll degradation for tocopherol synthesis, but bulk chlorophyll levels were unaffected in vte7 mutants and overexpressing lines. Instead, levels of specific chlorophyll biosynthetic intermediates containing partially reduced side chains were impacted and strongly correlated with totalT. These intermediates are generated by a membrane-associated biosynthetic complex containing protochlorophyllide reductase, chlorophyll synthase, geranylgeranyl reductase (GGR) and light harvesting-like 3 protein, all of which are required for both chlorophyll and tocopherol biosynthesis. We propose a model where VTE7 releases prenyl alcohols from chlorophyll biosynthetic intermediates, which are then converted to the corresponding diphosphates for tocopherol biosynthesis.

The biosynthesis of thymol, carvacrol, and thymohydroquinone in Lamiaceae proceeds via cytochrome P450s and a short-chain dehydrogenase.

Thymol and carvacrol are phenolic monoterpenes found in thyme, oregano, and several other species of the Lamiaceae. Long valued for their smell and taste, these substances also have antibacterial and anti-spasmolytic properties. They are also suggested to be precursors of thymohydroquinone and thymoquinone, monoterpenes with anti-inflammatory, antioxidant, and antitumor activities. Thymol and carvacrol biosynthesis has been proposed to proceed by the cyclization of geranyl diphosphate to γ-terpinene, followed by a series of oxidations via p-cymene. Here, we show that γ-terpinene is oxidized by cytochrome P450 monooxygenases (P450s) of the CYP71D subfamily to produce unstable cyclohexadienol intermediates, which are then dehydrogenated by a short-chain dehydrogenase/reductase (SDR) to the corresponding ketones. The subsequent formation of the aromatic compounds occurs via keto-enol tautomerisms. Combining these enzymes with γ-terpinene in in vitro assays or in vivo in Nicotiana benthamiana yielded thymol and carvacrol as products. In the absence of the SDRs, only p-cymene was formed by rearrangement of the cyclohexadienol intermediates. The nature of these unstable intermediates was inferred from reactions with the γ-terpinene isomer limonene and by analogy to reactions catalyzed by related enzymes. We also identified and characterized two P450s of the CYP76S and CYP736A subfamilies that catalyze the hydroxylation of thymol and carvacrol to thymohydroquinone when heterologously expressed in yeast and N. benthamiana Our findings alter previous views of thymol and carvacrol formation, identify the enzymes involved in the biosynthesis of these phenolic monoterpenes and thymohydroquinone in the Lamiaceae, and provide targets for metabolic engineering of high-value terpenes in plants.

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.

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.

Impact of choice of future climate change projection on growth chamber experimental outcomes: a preliminary study in potato.

Understanding the impacts of climate change on agriculture is essential to ensure adequate future food production. Controlled growth experiments provide an effective tool for assessing the complex effects of climate change. However, a review of the use of climate projections in 57 previously published controlled growth studies found that none considered within-season variations in projected future temperature change, and few considered regional differences in future warming. A fixed, often arbitrary, temperature perturbation typically was applied for the entire growing season. This study investigates the utility of employing more complex climate change scenarios in growth chamber experiments. A case study in potato was performed using three dynamically downscaled climate change projections for the mid-twenty-first century that differ in terms of the timing during the growing season of the largest projected temperature changes. The climate projections were used in growth chamber experiments for four elite potato cultivars commonly planted in Michigan's major potato growing region. The choice of climate projection had a significant influence on the sign and magnitude of the projected changes in aboveground biomass and total tuber count, whereas all projections suggested an increase in total tuber weight and a decrease in specific gravity, a key market quality trait for potato, by mid-century. These results demonstrate that the use of more complex climate projections that extend beyond a simple incremental change can provide additional insights into the future impacts of climate change on crop production and the accompanying uncertainty.

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.

Chromosome-scale genome assembly of the 'Munstead' cultivar of Lavandula angustifolia.

OBJECTIVES: Lavandula angustifolia (English lavender) is commercially important not only as an ornamental species but also as a major source of fragrances. To better understand the genomic basis of chemical diversity in lavender, we sequenced, assembled, and annotated the 'Munstead' cultivar of L. angustifolia. DATA DESCRIPTION: A total of 80 Gb of Oxford Nanopore Technologies reads was used to assemble the 'Munstead' genome using the Canu genome assembler software. Following multiple rounds of error correction and scaffolding using Hi-C data, the final chromosome-scale assembly represents 795,075,733 bp across 25 chromosomes with an N50 scaffold length of 31,371,815 bp. Benchmarking Universal Single Copy Orthologs analysis revealed 98.0% complete orthologs, indicative of a high-quality assembly representative of genic space. Annotation of protein-coding sequences revealed 58,702 high-confidence genes encoding 88,528 gene models. Access to the 'Munstead' genome will permit comparative analyses within and among lavender accessions and provides a pivotal species for comparative analyses within Lamiaceae.

Identification of genes associated with abiotic stress tolerance in sweetpotato using weighted gene co-expression network analysis.

Sweetpotato, Ipomoea batatas (L.), a key food security crop, is negatively impacted by heat, drought, and salinity stress. The orange-fleshed sweetpotato cultivar "Beauregard" was exposed to heat, salt, and drought treatments for 24 and 48 h to identify genes responding to each stress condition in leaves. Analysis revealed both common (35 up regulated, 259 down regulated genes in the three stress conditions) and unique sets of up regulated (1337 genes by drought, 516 genes by heat, and 97 genes by salt stress) and down regulated (2445 genes by drought, 678 genes by heat, and 204 genes by salt stress) differentially expressed genes (DEGs) suggesting common, yet stress-specific transcriptional responses to these three abiotic stressors. Gene Ontology analysis of down regulated DEGs common to both heat and salt stress revealed enrichment of terms associated with "cell population proliferation" suggestive of an impact on the cell cycle by the two stress conditions. To identify shared and unique gene co-expression networks under multiple abiotic stress conditions, weighted gene co-expression network analysis was performed using gene expression profiles from heat, salt, and drought stress treated 'Beauregard' leaves yielding 18 co-expression modules. One module was enriched for "response to water deprivation," "response to abscisic acid," and "nitrate transport" indicating synergetic crosstalk between nitrogen, water, and phytohormones with genes encoding osmotin, cell expansion, and cell wall modification proteins present as key hub genes in this drought-associated module. This research lays the groundwork for exploring to a further degree, mechanisms for abiotic stress tolerance in sweetpotato.

Cold stress induces differential gene expression of retained homeologs in Camelina sativa cv Suneson.

Camelina sativa (L.) Crantz, a member of the Brassicaceae, has potential as a biofuel feedstock which is attributable to the production of fatty acids in its seeds, its fast growth cycle, and low input requirements. While a genome assembly is available for camelina, it was generated from short sequence reads and is thus highly fragmented in nature. Using long read sequences, we generated a chromosome-scale, highly contiguous genome assembly (644,491,969 bp) for the spring biotype cultivar 'Suneson' with an N50 contig length of 12,031,512 bp and a scaffold N50 length of 32,184,682 bp. Annotation of protein-coding genes revealed 91,877 genes that encode 133,355 gene models. We identified a total of 4,467 genes that were significantly up-regulated under cold stress which were enriched in gene ontology terms associated with "response to cold" and "response to abiotic stress". Coexpression analyses revealed multiple coexpression modules that were enriched in genes differentially expressed following cold stress that had putative functions involved in stress adaptation, specifically within the plastid. With access to a highly contiguous genome assembly, comparative analyses with Arabidopsis thaliana revealed 23,625 A. thaliana genes syntenic with 45,453 Suneson genes. Of these, 24,960 Suneson genes were syntenic to 8,320 A. thaliana genes reflecting a 3 camelina homeolog to 1 Arabidopsis gene relationship and retention of all three homeologs. Some of the retained triplicated homeologs showed conserved gene expression patterns under control and cold-stressed conditions whereas other triplicated homeologs displayed diverged expression patterns revealing sub- and neo-functionalization of the homeologs at the transcription level. Access to the chromosome-scale assembly of Suneson will enable both basic and applied research efforts in the improvement of camelina as a sustainable biofuel feedstock.

Transcriptome-wide association and prediction for carotenoids and tocochromanols in fresh sweet corn kernels.

Sweet corn (Zea mays L.) is consistently one of the most highly consumed vegetables in the United States, providing a valuable opportunity to increase nutrient intake through biofortification. Significant variation for carotenoid (provitamin A, lutein, zeaxanthin) and tocochromanol (vitamin E, antioxidants) levels is present in temperate sweet corn germplasm, yet previous genome-wide association studies (GWAS) of these traits have been limited by low statistical power and mapping resolution. Here, we employed a high-quality transcriptomic dataset collected from fresh sweet corn kernels to conduct transcriptome-wide association studies (TWAS) and transcriptome prediction studies for 39 carotenoid and tocochromanol traits. In agreement with previous GWAS findings, TWAS detected significant associations for four causal genes, ß-carotene hydroxylase (crtRB1), lycopene epsilon cyclase (lcyE), γ-tocopherol methyltransferase (vte4), and homogentisate geranylgeranyltransferase (hggt1) on a transcriptome-wide level. Pathway-level analysis revealed additional associations for deoxy-xylulose synthase2 (dxs2), diphosphocytidyl methyl erythritol synthase2 (dmes2), cytidine methyl kinase1 (cmk1), and geranylgeranyl hydrogenase1 (ggh1), of which, dmes2, cmk1, and ggh1 have not previously been identified through maize association studies. Evaluation of prediction models incorporating genome-wide markers and transcriptome-wide abundances revealed a trait-dependent benefit to the inclusion of both genomic and transcriptomic data over solely genomic data, but both transcriptome- and genome-wide datasets outperformed a priori candidate gene-targeted prediction models for most traits. Altogether, this study represents an important step toward understanding the role of regulatory variation in the accumulation of vitamins in fresh sweet corn kernels.

Biocatalytic routes to stereo-divergent iridoids.

Thousands of natural products are derived from the fused cyclopentane-pyran molecular scaffold nepetalactol. These natural products are used in an enormous range of applications that span the agricultural and medical industries. For example, nepetalactone, the oxidized derivative of nepetalactol, is known for its cat attractant properties as well as potential as an insect repellent. Most of these naturally occurring nepetalactol-derived compounds arise from only two out of the eight possible stereoisomers, 7S-cis-trans and 7R-cis-cis nepetalactols. Here we use a combination of naturally occurring and engineered enzymes to produce seven of the eight possible nepetalactol or nepetalactone stereoisomers. These enzymes open the possibilities for biocatalytic production of a broader range of iridoids, providing a versatile system for the diversification of this important natural product scaffold.

Combining GWAS and TWAS to identify candidate causal genes for tocochromanol levels in maize grain.

Tocochromanols (tocopherols and tocotrienols, collectively vitamin E) are lipid-soluble antioxidants important for both plant fitness and human health. The main dietary sources of vitamin E are seed oils that often accumulate high levels of tocopherol isoforms with lower vitamin E activity. The tocochromanol biosynthetic pathway is conserved across plant species but an integrated view of the genes and mechanisms underlying natural variation of tocochromanol levels in seed of most cereal crops remains limited. To address this issue, we utilized the high mapping resolution of the maize Ames panel of ∼1,500 inbred lines scored with 12.2 million single-nucleotide polymorphisms to generate metabolomic (mature grain tocochromanols) and transcriptomic (developing grain) data sets for genetic mapping. By combining results from genome- and transcriptome-wide association studies, we identified a total of 13 candidate causal gene loci, including 5 that had not been previously associated with maize grain tocochromanols: 4 biosynthetic genes (arodeH2 paralog, dxs1, vte5, and vte7) and a plastid S-adenosyl methionine transporter (samt1). Expression quantitative trait locus (eQTL) mapping of these 13 gene loci revealed that they are predominantly regulated by cis-eQTL. Through a joint statistical analysis, we implicated cis-acting variants as responsible for colocalized eQTL and GWAS association signals. Our multiomics approach provided increased statistical power and mapping resolution to enable a detailed characterization of the genetic and regulatory architecture underlying tocochromanol accumulation in maize grain and provided insights for ongoing biofortification efforts to breed and/or engineer vitamin E and antioxidant levels in maize and other cereals.

Leveraging prior biological knowledge improves prediction of tocochromanols in maize grain.

With an essential role in human health, tocochromanols are mostly obtained by consuming seed oils; however, the vitamin E content of the most abundant tocochromanols in maize (Zea mays L.) grain is low. Several large-effect genes with cis-acting variants affecting messenger RNA (mRNA) expression are mostly responsible for tocochromanol variation in maize grain, with other relevant associated quantitative trait loci (QTL) yet to be fully resolved. Leveraging existing genomic and transcriptomic information for maize inbreds could improve prediction when selecting for higher vitamin E content. Here, we first evaluated a multikernel genomic best linear unbiased prediction (MK-GBLUP) approach for modeling known QTL in the prediction of nine tocochromanol grain phenotypes (12-21 QTL per trait) within and between two panels of 1,462 and 242 maize inbred lines. On average, MK-GBLUP models improved predictive abilities by 7.0-13.6% when compared with GBLUP. In a second approach with a subset of 545 lines from the larger panel, the highest average improvement in predictive ability relative to GBLUP was achieved with a multi-trait GBLUP model (15.4%) that had a tocochromanol phenotype and transcript abundances in developing grain for a few large-effect candidate causal genes (1-3 genes per trait) as multiple response variables. Taken together, our study illustrates the enhancement of prediction models when informed by existing biological knowledge pertaining to QTL and candidate causal genes.

Phased, chromosome-scale genome assemblies of tetraploid potato reveal a complex genome, transcriptome, and predicted proteome landscape underpinning genetic diversity.

Cultivated potato is a clonally propagated autotetraploid species with a highly heterogeneous genome. Phased assemblies of six cultivars including two chromosome-scale phased genome assemblies revealed extensive allelic diversity, including altered coding and transcript sequences, preferential allele expression, and structural variation that collectively result in a highly complex transcriptome and predicted proteome, which are distributed across the homologous chromosomes. Wild species contribute to the extensive allelic diversity in tetraploid cultivars, demonstrating ancestral introgressions predating modern breeding efforts. As a clonally propagated autotetraploid that undergoes limited meiosis, dysfunctional and deleterious alleles are not purged in tetraploid potato. Nearly a quarter of the loci bore mutations are predicted to have a high negative impact on protein function, complicating breeder's efforts to reduce genetic load. The StCDF1 locus controls maturity, and analysis of six tetraploid genomes revealed that 12 allelic variants of StCDF1 are correlated with maturity in a dosage-dependent manner. Knowledge of the complexity of the tetraploid potato genome with its rampant structural variation and embedded deleterious and dysfunctional alleles will be key not only to implementing precision breeding of tetraploid cultivars but also to the construction of homozygous, diploid potato germplasm containing favorable alleles to capitalize on heterosis in F1 hybrids.

The tepary bean genome provides insight into evolution and domestication under heat stress.

Tepary bean (Phaseolus acutifolis A. Gray), native to the Sonoran Desert, is highly adapted to heat and drought. It is a sister species of common bean (Phaseolus vulgaris L.), the most important legume protein source for direct human consumption, and whose production is threatened by climate change. Here, we report on the tepary genome including exploration of possible mechanisms for resilience to moderate heat stress and a reduced disease resistance gene repertoire, consistent with adaptation to arid and hot environments. Extensive collinearity and shared gene content among these Phaseolus species will facilitate engineering climate adaptation in common bean, a key food security crop, and accelerate tepary bean improvement.

Genome sequencing of four culinary herbs reveals terpenoid genes underlying chemodiversity in the Nepetoideae.

Species within the mint family, Lamiaceae, are widely used for their culinary, cultural, and medicinal properties due to production of a wide variety of specialized metabolites, especially terpenoids. To further our understanding of genome diversity in the Lamiaceae and to provide a resource for mining biochemical pathways, we generated high-quality genome assemblies of four economically important culinary herbs, namely, sweet basil (Ocimum basilicum L.), sweet marjoram (Origanum majorana L.), oregano (Origanum vulgare L.), and rosemary (Rosmarinus officinalis L.), and characterized their terpenoid diversity through metabolite profiling and genomic analyses. A total 25 monoterpenes and 11 sesquiterpenes were identified in leaf tissue from the 4 species. Genes encoding enzymes responsible for the biosynthesis of precursors for mono- and sesqui-terpene synthases were identified in all four species. Across all 4 species, a total of 235 terpene synthases were identified, ranging from 27 in O. majorana to 137 in the tetraploid O. basilicum. This study provides valuable resources for further investigation of the genetic basis of chemodiversity in these important culinary herbs.

Construction of a chromosome-scale long-read reference genome assembly for potato.

BACKGROUND: Worldwide, the cultivated potato, Solanum tuberosum L., is the No. 1 vegetable crop and a critical food security crop. The genome sequence of DM1-3 516 R44, a doubled monoploid clone of S. tuberosum Group Phureja, was published in 2011 using a whole-genome shotgun sequencing approach with short-read sequence data. Current advanced sequencing technologies now permit generation of near-complete, high-quality chromosome-scale genome assemblies at minimal cost. FINDINGS: Here, we present an updated version of the DM1-3 516 R44 genome sequence (v6.1) using Oxford Nanopore Technologies long reads coupled with proximity-by-ligation scaffolding (Hi-C), yielding a chromosome-scale assembly. The new (v6.1) assembly represents 741.6 Mb of sequence (87.8%) of the estimated 844 Mb genome, of which 741.5 Mb is non-gapped with 731.2 Mb anchored to the 12 chromosomes. Use of Oxford Nanopore Technologies full-length complementary DNA sequencing enabled annotation of 32,917 high-confidence protein-coding genes encoding 44,851 gene models that had a significantly improved representation of conserved orthologs compared with the previous annotation. The new assembly has improved contiguity with a 595-fold increase in N50 contig size, 99% reduction in the number of contigs, a 44-fold increase in N50 scaffold size, and an LTR Assembly Index score of 13.56, placing it in the category of reference genome quality. The improved assembly also permitted annotation of the centromeres via alignment to sequencing reads derived from CENH3 nucleosomes. CONCLUSIONS: Access to advanced sequencing technologies and improved software permitted generation of a high-quality, long-read, chromosome-scale assembly and improved annotation dataset for the reference genotype of potato that will facilitate research aimed at improving agronomic traits and understanding genome evolution.