A Genome-Wide Identification and Comparative Analysis of the Heavy-Metal-Associated Gene Family in Cucurbitaceae Species and Their Role in Cucurbita pepo under Arsenic Stress.

The heavy-metal-associated (HMA) proteins are a class of PB1-type ATPases related to the intracellular transport and detoxification of metals. However, due to a lack of information regarding the HMA gene family in the Cucurbitaceae family, a comprehensive genome-wide analysis of the HMA family was performed in ten Cucurbitaceae species: Citrullus amarus, Citrullus colocynthis, Citrullus lanatus, Citrullus mucosospermus, Cucumis melo, Cucumis sativus, Cucurbita maxima, Cucurbita moschata, Cucurbita pepo, and Legenaria siceraria. We identified 103 Cucurbit HMA proteins with various members, ranging from 8 (Legenaria siceraria) to 14 (Cucurbita pepo) across species. The phylogenetic and structural analysis confirmed that the Cucurbitaceae HMA protein family could be further classified into two major clades: Zn/Co/Cd/Pb and Cu/Ag. The GO-annotation-based subcellular localization analysis predicted that all HMA gene family members were localized on membranes. Moreover, the analysis of conserved motifs and gene structure (intron/exon) revealed the functional divergence between clades. The interspecies microsynteny analysis demonstrated that maximum orthologous genes were found between species of the Citrullus genera. Finally, nine candidate HMA genes were selected, and their expression analysis was carried out via qRT-PCR in root, leaf, flower, and fruit tissues of C. pepo under arsenic stress. The expression pattern of the CpeHMA genes showed a distinct pattern of expression in root and shoot tissues, with a remarkable expression of CpeHMA6 and CpeHMA3 genes from the Cu/Ag clade. Overall, this study provides insights into the functional analysis of the HMA gene family in Cucurbitaceae species and lays down the basic knowledge to explore the role and mechanism of the HMA gene family to cope with arsenic stress conditions.

Lessons learned from metabolic engineering in hairy roots: Transcriptome and metabolic profile changes caused by Rhizobium-mediated plant transformation in Cucurbitaceae species.

Cucurbitaceae species are used in traditional medicine around the world. Cucurbitacins are highly oxygenated triterpenoids found in Cucurbitaceae species and exhibit potent anticancer activity alone and in combination with other existing chemotherapeutic drugs. Therefore, increasing production of these specialized metabolites is of great relevance. We recently showed that hairy roots of Cucurbita pepo can be used as a platform for metabolic engineering of cucurbitacins to modify their structure and increase their production. To study the changes in cucurbitacin accumulation upon formation of hairy roots, an empty vector (EV) control and Cucurbitacin inducing bHLH transcription factor 1 (CpCUCbH1)-overexpressing hairy roots of C. pepo were compared to untransformed (WT) roots. Whilst CpCUCbH1-overexpression increased production of cucurbitacins I and B by 5-fold, and cucurbitacin E by 3-fold when compared to EV lines, this increase was not significantly different when compared to WT roots. This indicated that Rhizobium rhizogenes transformation lowered the cucurbitacins levels in hairy roots, but that increasing expression of cucurbitacin biosynthetic genes by CpCUCbH1-overexpression restored cucurbitacin production to WT levels. Subsequent metabolomic and RNA-seq analysis indicated that the metabolic profile and transcriptome of hairy roots was significantly changed when compared to WT roots. Interestingly, it was observed that 11% of the differentially expressed genes were transcription factors. It was noteworthy that the majority of transcripts showing highest Pearson correlation coefficients to the Rhizobium rhizogenes genes rolB, rolC and ORF13a, were predicted to be transcription factors. In summary, hairy roots are an excellent platform for metabolic engineering of plant specialized metabolites, but these extensive transcriptome and metabolic profile changes should be considered in subsequent studies.

Metabolic engineering of cucurbitacins in Cucurbita pepo hairy roots.

In this paper we show that metabolic engineering in Cucurbita pepo hairy roots can be used to both effectively increase and modify cucurbitacins. Cucurbitacins are highly-oxygenated triterpenoids originally described in the Cucurbitaceae family, but have since been found in 15 taxonomically distant plant families. Cucurbitacin B, D, E and I are the most widespread amongst the Cucurbitaceae and they have both important biological and pharmacological activities. In this study C. pepo hairy roots were used as a platform to boost production and alter the structures of the afore mentioned cucurbitacins by metabolic engineering to potentially provide new or more desirable bioactivities. We report that the ability to induce cucurbitacin biosynthesis by basic Helix-Loop-Helix transcription factors is partially conserved within the Cucurbitaceae and therefore can potentially be used as a biotechnological tool to increase cucurbitacins in several genera of this family. Additionally, overexpression of a novel acyltransferase from cucurbitacin producing Iberis amara generates a hitherto undescribed acetylation at the C3-hydroxyl group of the cucurbitadienol backbone. While overexpression of the cytochromes P450 CsCYP88L2 and McCYP88L7 from Cucumis sativus and Momordica charantia (respectively), results in accumulation of new spectral feature as revealed by High resolution liquid chromatography mass spectroscopy analysis; the m/z of the new peak supports it might be a cucurbitacin hydroxylated at the C19 position in C. pepo hairy roots. Finally, this paper is a case study of how hairy roots can be used to metabolically engineer and introduce novel modifications in metabolic pathways that have not been fully elucidated.

A GPCR-based yeast biosensor for biomedical, biotechnological, and point-of-use cannabinoid determination.

Eukaryotic cells use G-protein coupled receptors to sense diverse signals, ranging from chemical compounds to light. Here, we exploit the remarkable sensing capacity of G-protein coupled receptors to construct yeast-based biosensors for real-life applications. To establish proof-of-concept, we focus on cannabinoids because of their neuromodulatory and immunomodulatory activities. We construct a CB2 receptor-based biosensor, optimize it to achieve high sensitivity and dynamic range, and prove its effectiveness in three applications of increasing difficulty. First, we screen a compound library to discover agonists and antagonists. Second, we analyze 54 plants to discover a new phytocannabinoid, dugesialactone. Finally, we develop a robust portable device, analyze body-fluid samples, and confidently detect designer drugs like JWH-018. These examples demonstrate the potential of yeast-based biosensors to enable diverse applications that can be implemented by non-specialists. Taking advantage of the extensive sensing repertoire of G-protein coupled receptors, this technology can be extended to detect numerous compounds.

An Independent Evolutionary Origin for Insect Deterrent Cucurbitacins in Iberis amara.

Pieris rapae and Phyllotreta nemorum are Brassicaceae specialists, but do not feed on Iberis amara spp. that contain cucurbitacins. The cucurbitacins are highly oxygenated triterpenoid, occurring widespread in cucurbitaceous species and in a few other plant families. Using de novo assembled transcriptomics from I. amara, gene co-expression analysis and comparative genomics, we unraveled the evolutionary origin of the insect deterrent cucurbitacins in I. amara. Phylogenetic analysis of five oxidosqualene cyclases and heterologous expression allowed us to identify the first committed enzyme in cucurbitacin biosynthesis in I. amara, cucurbitadienol synthase (IaCPQ). In addition, two species-specific cytochrome P450s (CYP708A16 and CYP708A15) were identified that catalyze the unique C16 and C22 hydroxylation of the cucurbitadienol backbone, enzymatic steps that have not been reported before. Furthermore, the draft genome assembly of I. amara showed that the IaCPQ was localized to the same scaffold together with CYP708A15 but spanning over 100 kb, this contrasts with the highly organized cucurbitacin gene cluster in the cucurbits. These results reveal that cucurbitacin biosynthesis has evolved convergently via different biosynthetic routes in different families rather than through divergence from an ancestral pathway. This study thus provides new insight into the mechanism of recurrent evolution and diversification of a plant defensive chemical.

Transgenic Kalanchoë blossfeldiana, Containing Individual rol Genes and Open Reading Frames Under 35S Promoter, Exhibit Compact Habit, Reduced Plant Growth, and Altered Ethylene Tolerance in Flowers.

Reduced growth habit is a desirable trait for ornamental potted plants and can successfully be obtained through Rhizobium rhizogenes transformation in a stable and heritable manner. Additionally, it can also be obtained by transformation with Agrobacterium tumefaciens harboring specific genes from R. rhizogenes. The bacterial T-DNA harbors four root oncogenic loci (rol) genes and 14 less known open reading frames (ORFs). The four rol genes, i.e., rolA, rolB, rolC, and rolD, are conceived as the common denominator for the compact phenotype and the other less characterized ORFs seem auxiliary but present a potential breeding target for less aberrant and/or more tailored phenotypes. In this study, Kalanchoë blossfeldiana 'Molly' was transformed with individual rol genes and selected ORFs in 35S overexpressing cassettes to comprehensively characterize growth traits, gene copy and expression, and ethylene tolerance of the flowers. An association of reduced growth habit, e.g. height and diameter, was observed for rolB2 and ORF14-2 when a transgene single copy and high gene expression were detected. Chlorophyll content was reduced in overexpressing lines compared to wild type (WT), except for one ΔORF13a (a truncated ORF13a, where SPXX DNA-binding motif is absent). The flower number severely decreased in the overexpressing lines compared to WT. The anthesis timing showed that WT opened the first flower at 68.9 ± 0.9 days and the overexpressing lines showed similar or up to 24 days delay in flowering. In general, a single or low relative gene copy insertion was correlated to higher gene expression, ca. 3 to 5-fold, in rolB and ΔORF13a lines, while in ORF14 such relation was not directly linked. The increased gene expression observed in rolB2 and ΔORF13a-2 contributed to reducing plant growth and a more compact habit. Tolerance of detached flowers to 0.5 µl L-1 ethylene was markedly higher for ORF14 with 66% less flower closure at day 3 compared to WT. The subcellular localization of rolC and ΔORF13a was investigated by transient expression in Nicotiana benthamiana and confocal images showed that rolC and ΔORF13a are soluble and localize in the cytoplasm being able to enter the nucleus.

Plant triterpenoids with bond-missing skeletons: biogenesis, distribution and bioactivity.

Covering: up to December 2018 The polycyclic ABCD(E) framework of triterpenoids can miss a single endocyclic C-C bond as a result of a modification of the cyclization cascade that triggers their formation (interrupted- or diverted cascades), or can be the result of post-cyclization ring cleavage by late-stage oxidative modifications (seco-triterpenoids). Because of mechanistic and biogenetic differences, ring opening associated with loss of a skeletal fragment, as in nor-seco-triterpenoids (limonoids, quassinoids), will not be covered, nor will compounds where ring opening is part of a fragmentation cascade or of a multiple diversion from it. Even with these limitations, 342 bond-missing triterpenoids could be retrieved from the literature, with transversal distribution in the plant kingdom. Their structural diversity translates into a variety of biological targets, with dominance of potential applications in the realm of cancer, neuroprotection, and anti-infective therapy. In addition to the bioactivity and chemotaxonomic relevance of bond-missing triterpenoids, current knowledge on the genetic basis of interrupted- and diverted oxidosqualene cyclases will be summarized. This untapped source of enzymes could be useful to selectively modify triterpenoids by metabolic engineering, circumventing the bottlenecks of their isolation (poor yield or inadequate supply chain) to explore new areas of their chemical space.

Evolution of Structural Diversity of Triterpenoids.

Plants have evolved to produce a blend of specialized metabolites that serve functional roles in plant adaptation. Among them, triterpenoids are one of the largest subclasses of such specialized metabolites, with more than 14,000 known structures. They play a role in plant defense and development and have potential applications within food and pharma. Triterpenoids are cyclized from oxidized squalene precursors by oxidosqualene cyclases, creating more than 100 different cyclical triterpene scaffolds. This limited number of scaffolds is the first step towards creating the vast structural diversity of triterpenoids followed by extensive diversification, in particular, by oxygenation and glycosylation. Gene duplication, divergence, and selection are major forces that drive triterpenoid structural diversification. The triterpenoid biosynthetic genes can be organized in non-homologous gene clusters, such as in Avena spp., Cucurbitaceae and Solanum spp., or scattered along plant chromosomes as in Barbarea vulgaris. Paralogous genes organized as tandem repeats reflect the extended gene duplication activities in the evolutionary history of the triterpenoid saponin pathways, as seen in B. vulgaris. We review and discuss examples of convergent and divergent evolution in triterpenoid biosynthesis, and the apparent mechanisms occurring in plants that drive their increasing structural diversity within and across species. Using B. vulgaris' saponins as examples, we discuss the impact a single structural modification can have on the structure of a triterpenoid and how this affect its biological properties. These examples provide insight into how plants continuously evolve their specialized metabolome, opening the way to study uncharacterized triterpenoid biosynthetic pathways.

Co-expression of squalene epoxidases with triterpene cyclases boosts production of triterpenoids in plants and yeast.

Triterpene cyclases catalyze the first committed step in triterpene biosynthesis, by forming mono- to pentacyclic backbone structures from oxygenated C30 isoprenoid precursors. Squalene epoxidase precedes this cyclization by providing the oxygenated and activated substrate for triterpene biosynthesis. Three squalene epoxidases from Cucurbita pepo (CpSEs) were isolated and shown to have evolved under purifying selection with signs of sites under positive selection in their N- and C-termini. They all localize to the Endoplasmic Reticulum (ER) and produce 2,3-oxidosqualene and 2,3:22,23-dioxidosqualene when expressed in a yeast erg1 (squalene epoxidase) erg7 (lanosterol synthase) double mutant. Co-expression of the CpSEs with four different triterpene cyclases, either transiently in Nicotiana benthamiana or constitutively in yeast, showed that CpSEs boost triterpene production. CpSE2 was the best performing in this regard, which could reflect either increased substrate production or superior channeling of the substrate to the triterpene cyclases. Fluorescence Lifetime Imaging Microscopy (FLIM) analysis with C. pepo cucurbitadienol synthase (CpCPQ) revealed a specific interaction with CpSE2 but not with the other CpSEs. When CpSE2 was transformed into C. pepo hairy root lines, cucurbitacin E production was increased two folds compared to empty vector control lines. This study provides new insight into the importance of SEs in triterpene biosynthesis, suggesting that they may facilitate substrate channeling, and demonstrates that SE overexpression is a new tool for increasing triterpene production in plants and yeast.

A Single Oxidosqualene Cyclase Produces the Seco-Triterpenoid α-Onocerin.

8,14-seco-Triterpenoids are characterized by their unusual open C-ring. Their distribution in nature is rare and scattered in taxonomically unrelated plants. The 8,14-seco-triterpenoid α-onocerin is only known from the evolutionarily distant clubmoss genus Lycopodium and the leguminous genus Ononis, which makes the biosynthesis of this seco-triterpenoid intriguing from an evolutionary standpoint. In our experiments with Ononis spinosa, α-onocerin was detected only in the roots. Through transcriptome analysis of the roots, an oxidosqualene cyclase, OsONS1, was identified that produces α-onocerin from squalene-2,3;22,23-dioxide when transiently expressed in Nicotiana bethamiana In contrast, in Lycopodium clavatum, two sequential cyclases, LcLCC and LcLCD, are required to produce α-onocerin in the N. benthamiana transient expression system. Expression of OsONS1 in the lanosterol synthase knockout yeast strain GIL77, which accumulates squalene-2,3;22,23-dioxide, verified the α-onocerin production. A phylogenetic analysis predicts that OsONS1 branches off from specific lupeol synthases and does not group with the known L. clavatum α-onocerin cyclases. Both the biochemical and phylogenetic analyses of OsONS1 suggest convergent evolution of the α-onocerin pathways. When OsONS1 was coexpressed in N. benthamiana leaves with either of the two O. spinosa squalene epoxidases, OsSQE1 or OsSQE2, α-onocerin production was boosted, most likely because the epoxidases produce higher amounts of squalene-2,3;22,23-dioxide. Fluorescence lifetime imaging microscopy analysis demonstrated specific protein-protein interactions between OsONS1 and both O. spinosa squalene epoxidases. Coexpression of OsONS1 with the two OsSQEs suggests that OsSQE2 is the preferred partner of OsONS1 in planta. Our results provide an example of the convergent evolution of plant specialized metabolism.

Genome-wide Diversity and Association Mapping for Capsaicinoids and Fruit Weight in Capsicum annuum L.

Accumulated capsaicinoid content and increased fruit size are traits resulting from Capsicum annuum domestication. In this study, we used a diverse collection of C. annuum to generate 66,960 SNPs using genotyping by sequencing. The study identified 1189 haplotypes containing 3413 SNPs. Length of individual linkage disequilibrium (LD) blocks varied along chromosomes, with regions of high and low LD interspersed with an average LD of 139 kb. Principal component analysis (PCA), Bayesian model based population structure analysis and an Euclidean tree built based on identity by state (IBS) indices revealed that the clustering pattern of diverse accessions are in agreement with capsaicin content (CA) and fruit weight (FW) classifications indicating the importance of these traits in shaping modern pepper genome. PCA and IBS were used in a mixed linear model of capsaicin and dihydrocapsaicin content and fruit weight to reduce spurious associations because of confounding effects of subpopulations in genome-wide association study (GWAS). Our GWAS results showed SNPs in Ankyrin-like protein, IKI3 family protein, ABC transporter G family and pentatricopeptide repeat protein are the major markers for capsaicinoids and of 16 SNPs strongly associated with FW in both years of the study, 7 are located in known fruit weight controlling genes.

Genome-Wide Divergence and Linkage Disequilibrium Analyses for Capsicum baccatum Revealed by Genome-Anchored Single Nucleotide Polymorphisms.

Principal component analysis (PCA) with 36,621 polymorphic genome-anchored single nucleotide polymorphisms (SNPs) identified collectively for Capsicum annuum and Capsicum baccatum was used to characterize population structure and species domestication of these two important incompatible cultivated pepper species. Estimated mean nucleotide diversity (π) and Tajima's D across various chromosomes revealed biased distribution toward negative values on all chromosomes (except for chromosome 4) in cultivated C. baccatum, indicating a population bottleneck during domestication of C. baccatum. In contrast, C. annuum chromosomes showed positive π and Tajima's D on all chromosomes except chromosome 8, which may be because of domestication at multiple sites contributing to wider genetic diversity. For C. baccatum, 13,129 SNPs were available, with minor allele frequency (MAF) ≥0.05; PCA of the SNPs revealed 283 C. baccatum accessions grouped into 3 distinct clusters, for strong population structure. The fixation index (FST ) between domesticated C. annuum and C. baccatum was 0.78, which indicates genome-wide divergence. We conducted extensive linkage disequilibrium (LD) analysis of C. baccatum var. pendulum cultivars on all adjacent SNP pairs within a chromosome to identify regions of high and low LD interspersed with a genome-wide average LD block size of 99.1 kb. We characterized 1742 haplotypes containing 4420 SNPs (range 9-2 SNPs per haplotype). Genome-wide association study (GWAS) of peduncle length, a trait that differentiates wild and domesticated C. baccatum types, revealed 36 significantly associated genome-wide SNPs. Population structure, identity by state (IBS) and LD patterns across the genome will be of potential use for future GWAS of economically important traits in C. baccatum peppers.

Differential gene expression and alternative splicing between diploid and tetraploid watermelon.

The exploitation of synthetic polyploids for producing seedless fruits is well known in watermelon. Tetraploid progenitors of triploid watermelon plants, compared with their diploid counterparts, exhibit wide phenotypic differences. Although many factors modulate alternative splicing (AS) in plants, the effects of autopolyploidization on AS are still unknown. In this study, we used tissues of leaf, stem, and fruit of diploid and tetraploid sweet watermelon to understand changes in gene expression and the occurrence of AS. RNA-sequencing analysis was performed along with reverse transcription quantitative PCR and rapid amplification of cDNA ends (RACE)-PCR to demonstrate changes in expression and splicing. All vegetative tissues except fruit showed an increased level of AS in the tetraploid watermelon throughout the growth period. The ploidy levels of diploids and the tetraploid were confirmed using a ploidy analyser. We identified 5362 and 1288 genes that were up- and downregulated, respectively, in tetraploid as compared with diploid plants. We further confirmed that 22 genes underwent AS events across tissues, indicating possibilities of generating different protein isoforms with altered functions of important transcription factors and transporters. Arginine biosynthesis, chlorophyllide synthesis, GDP mannose biosynthesis, trehalose biosynthesis, and starch and sucrose degradation pathways were upregulated in autotetraploids. Phloem protein 2, chloroplastic PGR5-like protein, zinc-finger protein, fructokinase-like 2, MYB transcription factor, and nodulin MtN21 showed AS in fruit tissues. These results should help in developing high-quality seedless watermelon and provide additional transcriptomic information related to other cucurbits.

Identification of gene-specific polymorphisms and association with capsaicin pathway metabolites in Capsicum annuum L. collections.

Pepper (Capsicum annuum L.) is an economically important crop with added nutritional value. Production of capsaicin is an important quantitative trait with high environmental variance, so the development of markers regulating capsaicinoid accumulation is important for pepper breeding programs. In this study, we performed association mapping at the gene level to identify single nucleotide polymorphisms (SNPs) associated with capsaicin pathway metabolites in a diverse Capsicum annuum collection during two seasons. The genes Pun1, CCR, KAS and HCT were sequenced and matched with the whole-genome sequence draft of pepper to identify SNP locations and for further characterization. The identified SNPs for each gene underwent candidate gene association mapping. Association mapping results revealed Pun1 as a key regulator of major metabolites in the capsaicin pathway mainly affecting capsaicinoids and precursors for acyl moieties of capsaicinoids. Six different SNPs in the promoter sequence of Pun1 were found associated with capsaicin in plants from both seasons. Our results support that CCR is an important control point for the flux of p-coumaric acid to specific biosynthesis pathways. KAS was found to regulate the major precursors for acyl moieties of capsaicinoids and may play a key role in capsaicinoid production. Candidate gene association mapping of Pun1 suggested that the accumulation of capsaicinoids depends on the expression of Pun1, as revealed by the most important associated SNPs found in the promoter region of Pun1.