The genome and population genomics of allopolyploid Coffea arabica reveal the diversification history of modern coffee cultivars.

Coffea arabica, an allotetraploid hybrid of Coffea eugenioides and Coffea canephora, is the source of approximately 60% of coffee products worldwide, and its cultivated accessions have undergone several population bottlenecks. We present chromosome-level assemblies of a di-haploid C. arabica accession and modern representatives of its diploid progenitors, C. eugenioides and C. canephora. The three species exhibit largely conserved genome structures between diploid parents and descendant subgenomes, with no obvious global subgenome dominance. We find evidence for a founding polyploidy event 350,000-610,000 years ago, followed by several pre-domestication bottlenecks, resulting in narrow genetic variation. A split between wild accessions and cultivar progenitors occurred ~30.5 thousand years ago, followed by a period of migration between the two populations. Analysis of modern varieties, including lines historically introgressed with C. canephora, highlights their breeding histories and loci that may contribute to pathogen resistance, laying the groundwork for future genomics-based breeding of C. arabica.

A cytochrome P450 CYP87A4 imparts sterol side-chain cleavage in digoxin biosynthesis.

Digoxin extracted from the foxglove plant is a widely prescribed natural product for treating heart failure. It is listed as an essential medicine by the World Health Organization. However, how the foxglove plant synthesizes digoxin is mostly unknown, especially the cytochrome P450 sterol side chain cleaving enzyme (P450scc), which catalyzes the first and rate-limiting step. Here we identify the long-speculated foxglove P450scc through differential transcriptomic analysis. This enzyme converts cholesterol and campesterol to pregnenolone, suggesting that digoxin biosynthesis starts from both sterols, unlike previously reported. Phylogenetic analysis indicates that this enzyme arises from a duplicated cytochrome P450 CYP87A gene and is distinct from the well-characterized mammalian P450scc. Protein structural analysis reveals two amino acids in the active site critical for the foxglove P450scc's sterol cleavage ability. Identifying the foxglove P450scc is a crucial step toward completely elucidating digoxin biosynthesis and expanding the therapeutic applications of digoxin analogs in future work.

A well-characterized polycistronic-like gene expression system in yeast.

Efficient expression of multiple genes is critical to yeast metabolic engineering for the bioproduction of bulk and fine chemicals. A yeast polycistronic expression system is of particular interest because one promoter can drive the expression of multiple genes. 2A viral peptides enable the cotranslation of multiple proteins from a single mRNA by ribosomal skipping. However, the wide adaptation of 2A viral peptides for polycistronic-like gene expression in yeast awaits in-depth characterizations. Additionally, a one-step assembly of such a polycistronic-like system is highly desirable. To this end, we have developed a modular cloning (MoClo) compatible 2A peptide-based polycistronic-like system capable of expressing multiple genes from a single promoter in yeast. Characterizing the bi-, tri-, and quad-cistronic expression of fluorescent proteins showed high cleavage efficiencies of three 2A peptides: E2A from equine rhinitis B virus, P2A from porcine teschovirus-1, and O2A from Operophtera brumata cypovirus-18. Applying the polycistronic-like system to produce geraniol, a valuable industrial compound, resulted in comparable or higher titers than using conventional monocistronic constructs. In summary, this highly-characterized polycistronic-like gene expression system is another tool to facilitate multigene expression for metabolic engineering in yeast.

Machine-learning guided elucidation of contribution of individual steps in the mevalonate pathway and construction of a yeast platform strain for terpenoid production.

The production of terpenoids from engineered microbes contributes markedly to the bioeconomy by providing essential medicines, sustainable materials, and renewable fuels. The mevalonate pathway leading to the synthesis of terpenoid precursors has been extensively targeted for engineering. Nevertheless, the importance of individual pathway enzymes to the overall pathway flux and final terpenoid yield is less known, especially enzymes that are thought to be non-rate-limiting. To investigate the individual contribution of the five non-rate-limiting enzymes in the mevalonate pathway, we created a combinatorial library of 243 Saccharomyces cerevisiae strains, each having an extra copy of the mevalonate pathway integrated into the genome and expressing the non-rate-limiting enzymes from a unique combination of promoters. High-throughput screening combined with machine learning algorithms revealed that the mevalonate kinase, Erg12p, stands out as the critical enzyme that influences product titer. ERG12 is ideally expressed from a medium-strength promoter which is the 'sweet spot' resulting in high product yield. Additionally, a platform strain was created by targeting the mevalonate pathway to both the cytosol and peroxisomes. The dual localization synergistically increased terpenoid production and implied that some mevalonate pathway intermediates, such as mevalonate, isopentyl pyrophosphate (IPP), and dimethylallyl pyrophosphate (DMAPP), are diffusible across peroxisome membranes. The platform strain resulted in 94-fold, 60-fold, and 35-fold improved titer of monoterpene geraniol, sesquiterpene α-humulene, and triterpene squalene, respectively. The terpenoid platform strain will serve as a chassis for producing any terpenoids and terpene derivatives.

Rapid Assembly of Multi-Gene Constructs using Modular Golden Gate Cloning.

The Golden Gate cloning method enables the rapid assembly of multiple genes in any user-defined arrangement. It utilizes type IIS restriction enzymes that cut outside of their recognition sites and create a short overhang. This modular cloning (MoClo) system uses a hierarchical workflow in which different DNA parts, such as promoters, coding sequences (CDS), and terminators, are first cloned into an entry vector. Multiple entry vectors then assemble into transcription units. Several transcription units then connect into a multi-gene plasmid. The Golden Gate cloning strategy is of tremendous advantage because it allows scar-less, directional, and modular assembly in a one-pot reaction. The hierarchical workflow typically enables the facile cloning of a large variety of multi-gene constructs with no need for sequencing beyond entry vectors. The use of fluorescent protein dropouts enables easy visual screening. This work provides a detailed, step-by-step protocol for assembling multi-gene plasmids using the yeast modular cloning (MoClo) kit. We show optimal and suboptimal results of multi-gene plasmid assembly and provide a guide for screening for colonies. This cloning strategy is highly applicable for yeast metabolic engineering and other situations in which multi-gene plasmid cloning is required.

Genome-wide SNP identification and characterization in two soybean cultivars with contrasting Mungbean Yellow Mosaic India Virus disease resistance traits.

Mungbean yellow mosaic India virus (MYMIV) is a bipartite Geminivirus, which causes severe yield loss in soybean (Glycine max). Considering this, the present study was conducted to develop large-scale genome-wide single nucleotide polymorphism (SNP) markers and identify potential markers linked with known disease resistance loci for their effective use in genomics-assisted breeding to impart durable MYMIV tolerance. The whole-genome re-sequencing of MYMIV resistant cultivar 'UPSM-534' and susceptible Indian cultivar 'JS-335' was performed to identify high-quality SNPs and InDels (insertion and deletions). Approximately 234 and 255 million of 100-bp paired-end reads were generated from UPSM-534 and JS-335, respectively, which provided ~98% coverage of reference soybean genome. A total of 3083987 SNPs (1559556 in UPSM-534 and 1524431 in JS-335) and 562858 InDels (281958 in UPSM-534 and 280900 in JS-335) were identified. Of these, 1514 SNPs were found to be present in 564 candidate disease resistance genes. Among these, 829 non-synonymous and 671 synonymous SNPs were detected in 266 and 286 defence-related genes, respectively. Noteworthy, a non-synonymous SNP (in chromosome 18, named 18-1861613) at the 149th base-pair of LEUCINE-RICH REPEAT RECEPTOR-LIKE PROTEIN KINASE gene responsible for a G/C transversion [proline (CCC) to alanine(GCC)] was identified and validated in a set of 12 soybean cultivars. Taken together, the present study generated a large-scale genomic resource such as, SNPs and InDels at a genome-wide scale that will facilitate the dissection of various complex traits through construction of high-density linkage maps and fine mapping. In the present scenario, these markers can be effectively used to design high-density SNP arrays for their large-scale validation and high-throughput genotyping in diverse natural and mapping populations, which could accelerate genomics-assisted MYMIV disease resistance breeding in soybean.