Non-Aqueous Isolation and Enrichment of Glandular Capitate Stalked and Sessile Trichomes from Cannabis sativa.

This paper presents a protocol for the convenient and high-throughput isolation and enrichment of glandular capitate stalked and sessile trichomes from Cannabis sativa. The biosynthetic pathways for cannabinoid and volatile terpene metabolism are localized primarily in the Cannabis trichomes, and isolated trichomes are beneficial for transcriptome analysis. The existing protocols for isolating glandular trichomes for transcriptomic characterization are inconvenient and deliver compromised trichome heads and a relatively low amount of isolated trichomes. Furthermore, they rely on expensive apparatus and isolation media containing protein inhibitors to avoid RNA degradation. The present protocol suggests combining three individual modifications to obtain a large amount of isolated glandular capitate stalked and sessile trichomes from C. sativa mature female inflorescences and fan leaves, respectively. The first modification involves substituting liquid nitrogen for the conventional isolation medium to facilitate the passage of trichomes through the micro-sieves. The second modification involves using dry ice to detach the trichomes from the plant source. The third modification involves passing the plant material consecutively through five micro-sieves of diminishing pore sizes. Microscopic imaging demonstrated the effectiveness of the isolation technique for both trichome types. In addition, the quality of RNA extracted from the isolated trichomes was appropriate for downstream transcriptomic analysis.

Pan-genome and multi-parental framework for high-resolution trait dissection in melon (Cucumis melo).

Linking genotype with phenotype is a fundamental goal in biology and requires robust data for both. Recent advances in plant-genome sequencing have expedited comparisons among multiple-related individuals. The abundance of structural genomic within-species variation that has been discovered indicates that a single reference genome cannot represent the complete sequence diversity of a species, leading to the expansion of the pan-genome concept. For high-resolution forward genetics, this unprecedented access to genomic variation should be paralleled and integrated with phenotypic characterization of genetic diversity. We developed a multi-parental framework for trait dissection in melon (Cucumis melo), leveraging a novel pan-genome constructed for this highly variable cucurbit crop. A core subset of 25 diverse founders (MelonCore25), consisting of 24 accessions from the two widely cultivated subspecies of C. melo, encompassing 12 horticultural groups, and 1 feral accession was sequenced using a combination of short- and long-read technologies, and their genomes were assembled de novo. The construction of this melon pan-genome exposed substantial variation in genome size and structure, including detection of ~300 000 structural variants and ~9 million SNPs. A half-diallel derived set of 300 F2 populations, representing all possible MelonCore25 parental combinations, was constructed as a framework for trait dissection through integration with the pan-genome. We demonstrate the potential of this unified framework for genetic analysis of various melon traits, including rind color intensity and pattern, fruit sugar content, and resistance to fungal diseases. We anticipate that utilization of this integrated resource will enhance genetic dissection of important traits and accelerate melon breeding.

2-oxoglutarate-dependent dioxygenases drive expansion of steroidal alkaloid structural diversity in the genus Solanum.

Solanum steroidal glycoalkaloids (SGAs) are renowned defence metabolites exhibiting spectacular structural diversity. Genes and enzymes generating the SGA precursor pathway, SGA scaffold and glycosylated forms have been largely identified. Yet, the majority of downstream metabolic steps creating the vast repertoire of SGAs remain untapped. Here, we discovered that members of the 2-OXOGLUTARATE-DEPENDENT DIOXYGENASE (2-ODD) family play a prominent role in SGA metabolism, carrying out three distinct backbone-modifying oxidative steps in addition to the three formerly reported pathway reactions. The GLYCOALKALOID METABOLISM34 (GAME34) enzyme catalyses the conversion of core SGAs to habrochaitosides in wild tomato S. habrochaites. Cultivated tomato plants overexpressing GAME34 ectopically accumulate habrochaitosides. These habrochaitoside enriched plants extracts potently inhibit Puccinia spp. spore germination, a significant Solanaceae crops fungal pathogen. Another 2-ODD enzyme, GAME33, acts as a desaturase (via hydroxylation and E/F ring rearrangement) forming unique, yet unreported SGAs. Conversion of bitter α-tomatine to ripe fruit, nonbitter SGAs (e.g. esculeoside A) requires two hydroxylations; while the known GAME31 2-ODD enzyme catalyses hydroxytomatine formation, we find that GAME40 catalyses the penultimate step in the pathway and generates acetoxy-hydroxytomatine towards esculeosides accumulation. Our results highlight the significant contribution of 2-ODD enzymes to the remarkable structural diversity found in plant steroidal specialized metabolism.

Extreme sugar accumulation in late fig ripening is accompanied by global changes in sugar metabolism and transporter gene expression.

Female fig (Ficus carica L.) fruit are characterized by a major increase in volume and sugar content during the final week of development. A detailed developmental analysis of water and dry matter accumulation during these final days indicated a temporal separation between the increase in volume due to increasing water content and a subsequent sharp increase in sugar content during a few days. The results present fig as an extreme example of sugar import and accumulation, with calculated import rates that are one order of magnitude higher than those of other sugar-accumulating sweet fruit species. To shed light on the metabolic changes occurring during this period, we followed the expression pattern of 80 genes encoding sugar metabolism enzymes and sugar transporter proteins identified in fig fruit. A parallel comparison with male fig fruits, which do not accumulate sugar during ripening, highlighted the genes specifically related to sugar accumulation. Tissue-specific analysis indicated that the expression of genes involved in sugar metabolism and transport undergoes a global transition.

QTL mapping and genomic analyses of earliness and fruit ripening traits in a melon Recombinant Inbred Lines population supported by de novo assembly of their parental genomes.

Earliness and ripening behavior are important attributes of fruits on and off the vine, and affect quality and preference of both growers and consumers. Fruit ripening is a complex physiological process that involves metabolic shifts affecting fruit color, firmness, and aroma production. Melon is a promising model crop for the study of fruit ripening, as the full spectrum of climacteric behavior is represented across the natural variation. Using Recombinant Inbred Lines (RILs) population derived from the parental lines "Dulce" (reticulatus, climacteric) and "Tam Dew" (inodorus, non-climacteric) that vary in earliness and ripening traits, we mapped QTLs for ethylene emission, fruit firmness and days to flowering and maturity. To further annotate the main QTL intervals and identify candidate genes, we used Oxford Nanopore long-read sequencing in combination with Illumina short-read resequencing, to assemble the parental genomes de-novo. In addition to 2.5 million genome-wide SNPs and short InDels detected between the parents, we also highlight here the structural variation between these lines and the reference melon genome. Through systematic multi-layered prioritization process, we identified 18 potential polymorphisms in candidate genes within multi-trait QTLs. The associations of selected SNPs with earliness and ripening traits were further validated across a panel of 177 diverse melon accessions and across a diallel population of 190 F1 hybrids derived from a core subset of 20 diverse parents. The combination of advanced genomic tools with diverse germplasm and targeted mapping populations is demonstrated as a way to leverage forward genetics strategies to dissect complex horticulturally important traits.

Underground heterosis for yield improvement in melon.

Heterosis, the superiority of hybrids over their parents, is a major genetic force associated with plant fitness and crop yield enhancement. We investigated root-mediated yield heterosis in melons (Cucumis melo) by characterizing a common variety grafted onto 190 hybrid rootstocks, resulting from crossing 20 diverse inbreds in a diallel-mating scheme. Hybrid rootstocks improved yield by more than 40% compared with their parents, and the best hybrid yield outperformed the reference commercial variety by 65% under both optimal and minimal irrigation treatments. To characterize the genetics of underground heterosis we conducted whole genome re-sequencing of the 20 founder lines, and showed that parental genetic distance was no predictor for the level of heterosis. Through inference of the 190 hybrid genotypes from their parental genomes, followed by genome-wide association analysis, we mapped multiple quantitative trait loci for root-mediated yield. Yield enhancement of the four best-performing hybrid rootstocks was validated in multiple experiments with four different scion varieties. Our grafting approach is complementary to the common roots genetic approach that focuses mainly on variation in root system architecture, and is a step towards discovery of candidate genes involved in root function and yield enhancement.

Prenyltransferases catalyzing geranyldiphosphate formation in tomato fruit.

Monoterpenes contribute either favorably or adversely to the flavor of tomato, yet modern tomato varieties generally lack monoterpenes in their fruit. The main immediate biosynthetic precursor of monoterpenes is geranyldiphosphate (GPP), produced by the action of GPP synthases (GPPSs). Plant GPPSs are often heteromeric enzymes consisting of a non-catalytic small subunit (GPPS.SSU) and a large subunit (GPPS.LSU), the latter similar to geranylgeranyldiphosphate synthases (GGPPSs) which generate longer prenylphosphate chains. We show here that LeGGPPS2, an enzyme previously reported to support carotenoid biosynthesis, can synthesize farnesyldiphosphate (FPP) and GPP in vitro, in addition to geranylgeranyldiphosphate, depending on the assay conditions. Moreover, GPP formation is favored in vitro by the interaction of LeGGPPS2 with GPPS.SSU from either Anthirrhinum majus (AmGPPS.SSU) or from a newly discovered GPPS.SSU ortholog present in the genome of M82 tomato. SlGPPS.SSU is not expressed in M82 tomato fruit but its orthologs are expressed in fruit of wild tomato relatives, such as Solanum pimpinelifollium and S. cheesmaniae that accumulate monoterpenes.

Sugar-regulated susceptibility of tomato fruit to Colletotrichum and Penicillium requires differential mechanisms of pathogenicity and fruit responses.

Colletotrichum gloeosporioides and Penicillium expansum cause postharvest diseases in tropical and deciduous fruit. During colonization, C. gloeosporioides and P. expansum secrete ammonia in hosts with low sugar content (LowSC) and gluconic acid in hosts with high sugar content (HighSC), respectively, as a mechanism to modulate enhanced pathogenicity. We studied the pathogens interactions with tomato lines of similar genetic background but differing in their sugar content. Colletotrichum gloeosporioides showed enhanced colonization of the LowSC line with differential expression response of 15% of its genes including enhanced relative expression of glycosyl hydrolases, glucanase and MFS-transporter genes. Enhanced colonization of P. expansum occurred in the HighSC line, accompanied by an increase in carbohydrate metabolic processes mainly phosphoenolpyruvate carboxykinase, and only 4% of differentially expressed genes. Gene response of the two host lines strongly differed depending on the sugar level. Limited colonization of HighSC line by C. gloeosporioides was accompanied by a marked alteration of gene expression compared the LowSC response to the same pathogen; while colonization by P. expansum resulted in a similar response of the two different hosts. We suggest that this differential pattern of fungal/host responses may be the basis for the differential of host range of both pathogens in nature.

Expression of Hexokinase in Stomata of Citrus Fruit Reduces Fruit Transpiration and Affects Seed Development.

The temporal formation and spatial distribution of stomata on the surface of citrus floral organs and, specifically, on the ovule from which the fruit develops, were analyzed using citrus plants that express green fluorescent protein (GFP) under the guard cell-specific KST1 promoter. Stomata are found on the style, sepal, and anther of the closed flower and on ovules from the stage of anthesis. It has previously been shown that hexokinase (HXK) mediates sugar-sensing in leaf guard cells and stimulates stomatal closure. The activity and response of citrus fruit stomata to sugar-sensing by HXK was examined using plants that express HXK under the KST1 promoter. Those plants are referred to as GCHXK plants. The transpiration of young green GCHXK citrus fruits was significantly reduced, indicating that their stomata respond to sugar similar to leaf stomata. Toward fruit maturation, fruit stomata are plugged and stop functioning, which explains why WT and GCHXK mature yellow fruits exhibited similar water loss. Seeds of the GCHXK plants were smaller and germinated more slowly than the WT seeds. We suggest that the stomata of young green citrus fruits, but not mature yellow fruits, respond to sugar levels via HXK and that fruit stomata are important for proper seed development.

Comparative Metabolomics and Molecular Phylogenetics of Melon (Cucumis melo, Cucurbitaceae) Biodiversity.

The broad variability of Cucumis melo (melon, Cucurbitaceae) presents a challenge to conventional classification and organization within the species. To shed further light on the infraspecific relationships within C. melo, we compared genotypic and metabolomic similarities among 44 accessions representative of most of the cultivar-groups. Genotyping-by-sequencing (GBS) provided over 20,000 single-nucleotide polymorphisms (SNPs). Metabolomics data of the mature fruit flesh and rind provided over 80,000 metabolomic and elemental features via an orchestra of six complementary metabolomic platforms. These technologies probed polar, semi-polar, and non-polar metabolite fractions as well as a set of mineral elements and included both flavor- and taste-relevant volatile and non-volatile metabolites. Together these results enabled an estimate of "metabolomic/elemental distance" and its correlation with the genetic GBS distance of melon accessions. This study indicates that extensive and non-targeted metabolomics/elemental characterization produced classifications that strongly, but not completely, reflect the current and extensive genetic classification. Certain melon Groups, such as Inodorous, clustered in parallel with the genetic classifications while other genome to metabolome/element associations proved less clear. We suggest that the combined genomic, metabolic, and element data reflect the extensive sexual compatibility among melon accessions and the breeding history that has, for example, targeted metabolic quality traits, such as taste and flavor.

High-density NGS-based map construction and genetic dissection of fruit shape and rind netting in Cucumis melo.

Melon is an important crop that exhibits broad variation for fruit morphology traits that are the substrate for genetic mapping efforts. In the post-genomic era, the link between genetic maps and physical genome assemblies is key for leveraging QTL mapping results for gene cloning and breeding purposes. Here, using a population of 164 melon recombinant inbred lines (RILs) that were subjected to genotyping-by-sequencing, we constructed and compared high-density sequence- and linkage-based recombination maps that were aligned to the reference melon genome. These analyses reveal the genome-wide variation in recombination frequency and highlight regions of disrupted collinearity between our population and the reference genome. The population was phenotyped over 3 years for fruit size and shape as well as rind netting. Four QTLs were detected for fruit size, and they act in an additive manner, while significant epistatic interaction was found between two neutral loci for this trait. Fruit shape displayed transgressive segregation that was explained by the action of four QTLs, contributed by alleles from both parents. The complexity of rind netting was demonstrated on a collection of 177 diverse accessions. Further dissection of netting in our RILs population, which is derived from a cross of smooth and densely netted parents, confirmed the intricacy of this trait and the involvement of major locus and several other interacting QTLs. A major netting QTL on chromosome 2 co-localized with results from two additional populations, paving the way for future study toward identification of a causative gene for this trait.

Low temperature upregulates cwp expression and modifies alternative splicing patterns, increasing the severity of cwp-induced tomato fruit cuticular microfissures.

The cwp (cuticular water permeability) gene controls the development of cuticular microfissuring and subsequent fruit dehydration in tomato. The gene underwent silencing in the evolution of the fleshy cultivated tomato but is expressed in the primitive wild tomato relatives. The introgression of the expressed allele from the wild S. habrochaites (cwp h ) into the cultivated tomato (Solanum lycopersicum) leads to the phenotype of fruit water loss during and following ripening. In this report, we show that low temperature impacts on the severity of the cuticular microfissure phenotype via a combination of effects on both expression and alternative splicing of cwp h . The cwp gene, comprising four exons and three introns, undergoes post-transcriptional alternative splicing processes, leading to seven alternative transcripts that differ in reading-frame lengths. Transgenic plants expressing each of the alternative transcripts identified the longest reading frame (VAR1) as the functional splice variant. Low temperature led to a strong upregulation of cwp h expression, compounded by an increase in the relative proportion of the functional VAR1 transcript, leading to increased severity of microfissuring of the cuticle. In summary, we demonstrate the molecular mechanism behind the horticultural phenomenon of the low-temperature effect on cuticular microfissures in the dehydrating tomato.

The multi-allelic APRR2 gene is associated with fruit pigment accumulation in melon and watermelon.

Color and pigment contents are important aspects of fruit quality and consumer acceptance of cucurbit crops. Here, we describe the independent mapping and cloning of a common causative APRR2 gene regulating pigment accumulation in melon and watermelon. We initially show that the APRR2 transcription factor is causative for the qualitative difference between dark and light green rind in both crops. Further analyses establish the link between sequence or expression level variations in the CmAPRR2 gene and pigment content in the rind and flesh of mature melon fruits. A genome-wide association study (GWAS) of young fruit rind color in a panel composed of 177 diverse melon accessions did not result in any significant association, leading to an earlier assumption that multiple genes are involved in shaping the overall phenotypic variation in this trait. Through resequencing of 25 representative accessions and allelism tests between light rind accessions, we show that multiple independent single nucleotide polymorphisms in the CmAPRR2 gene are causative of the light rind phenotype. The multi-haplotypic nature of this gene explains the lack of detection power obtained through genotyping by sequencing-based GWAS and confirms the pivotal role of this gene in shaping fruit color variation in melon. This study demonstrates the power of combining bi- and multi-allelic designs with deep sequencing, to resolve lack of power due to high haplotypic diversity and low allele frequencies. Due to its central role and broad effect on pigment accumulation in fruits, the APRR2 gene is an attractive target for carotenoid bio-fortification of cucurbit crops.

Cucurbit Genomics Database (CuGenDB): a central portal for comparative and functional genomics of cucurbit crops.

The Cucurbitaceae family (cucurbit) includes several economically important crops, such as melon, cucumber, watermelon, pumpkin, squash and gourds. During the past several years, genomic and genetic data have been rapidly accumulated for cucurbits. To store, mine, analyze, integrate and disseminate these large-scale datasets and to provide a central portal for the cucurbit research and breeding community, we have developed the Cucurbit Genomics Database (CuGenDB; http://cucurbitgenomics.org) using the Tripal toolkit. The database currently contains all available genome and expressed sequence tag (EST) sequences, genetic maps, and transcriptome profiles for cucurbit species, as well as sequence annotations, biochemical pathways and comparative genomic analysis results such as synteny blocks and homologous gene pairs between different cucurbit species. A set of analysis and visualization tools and user-friendly query interfaces have been implemented in the database to facilitate the usage of these large-scale data by the community. In particular, two new tools have been developed in the database, a 'SyntenyViewer' to view genome synteny between different cucurbit species and an 'RNA-Seq' module to analyze and visualize gene expression profiles. Both tools have been packed as Tripal extension modules that can be adopted in other genomics databases developed using the Tripal system.

Natural genetic variation for expression of a SWEET transporter among wild species of Solanum lycopersicum (tomato) determines the hexose composition of ripening tomato fruit.

The sugar content of Solanum lycopersicum (tomato) fruit is a primary determinant of taste and quality. Cultivated tomato fruit are characterized by near-equimolar levels of the hexoses glucose and fructose, derived from the hydrolysis of translocated sucrose. As fructose is perceived as approximately twice as sweet as glucose, increasing its concentration at the expense of glucose can improve tomato fruit taste. Introgressions of the FgrH allele from the wild species Solanum habrochaites (LA1777) into cultivated tomato increased the fructose-to-glucose ratio of the ripe fruit by reducing glucose levels and concomitantly increasing fructose levels. In order to identify the function of the Fgr gene, we combined a fine-mapping strategy with RNAseq differential expression analysis of near-isogenic tomato lines. The results indicated that a SWEET protein was strongly upregulated in the lines with a high fructose-to-glucose ratio. Overexpressing the SWEET protein in transgenic tomato plants dramatically reduced the glucose levels and increased the fructose : glucose ratio in the developing fruit, thereby proving the function of the protein. The SWEET protein was localized to the plasma membrane and expression of the SlFgr gene in a yeast line lacking native hexose transporters complemented growth with glucose, but not with fructose. These results indicate that the SlFgr gene encodes a plasma membrane-localized glucose efflux transporter of the SWEET family, the overexpression of which reduces glucose levels and may allow for increased fructose levels. This article identifies the function of the tomato Fgr gene as a SWEET transporter, the upregulation of which leads to a modified sugar accumulation pattern in the fleshy fruit. The results point to the potential of the inedible wild species to improve fruit sugar accumulation via sugar transport mechanisms.

Differential metabolism of L-phenylalanine in the formation of aromatic volatiles in melon (Cucumis melo L.) fruit.

Studies on the active pathways and the genes involved in the biosynthesis of L-phenylalanine-derived volatiles in fleshy fruits are sparse. Melon fruit rinds converted stable-isotope labeled L-phe into more than 20 volatiles. Phenylpropanes, phenylpropenes and benzenoids are apparently produced via the well-known phenylpropanoid pathway involving phenylalanine ammonia lyase (PAL) and being (E)-cinnamic acid a key intermediate. Phenethyl derivatives seemed to be derived from L-phe via a separate biosynthetic route not involving (E)-cinnamic acid and PAL. To explore for a biosynthetic route to (E)-cinnamaldehyde in melon rinds, soluble protein cell-free extracts were assayed with (E)-cinnamic acid, CoA, ATP, NADPH and MgSO4, producing (E)-cinnamaldehyde in vitro. In this context, we characterized CmCNL, a gene encoding for (E)-cinnamic acid:coenzyme A ligase, inferred to be involved in the biosynthesis of (E)-cinnamaldehyde. Additionally we describe CmBAMT, a SABATH gene family member encoding a benzoic acid:S-adenosyl-L-methionine carboxyl methyltransferase having a role in the accumulation of methyl benzoate. Our approach leads to a more comprehensive understanding of L-phe metabolism into aromatic volatiles in melon fruit.

Deciphering genetic factors that determine melon fruit-quality traits using RNA-Seq-based high-resolution QTL and eQTL mapping.

Combined quantitative trait loci (QTL) and expression-QTL (eQTL) mapping analysis was performed to identify genetic factors affecting melon (Cucumis melo) fruit quality, by linking genotypic, metabolic and transcriptomic data from a melon recombinant inbred line (RIL) population. RNA sequencing (RNA-Seq) of fruit from 96 RILs yielded a highly saturated collection of > 58 000 single-nucleotide polymorphisms, identifying 6636 recombination events that separated the genome into 3663 genomic bins. Bin-based QTL analysis of 79 RILs and 129 fruit-quality traits affecting taste, aroma and color resulted in the mapping of 241 QTL. Thiol acyltransferase (CmThAT1) gene was identified within the QTL interval of its product, S-methyl-thioacetate, a key component of melon fruit aroma. Metabolic activity of CmThAT1-encoded protein was validated in bacteria and in vitro. QTL analysis of flesh color intensity identified a candidate white-flesh gene (CmPPR1), one of two major loci determining fruit flesh color in melon. CmPPR1 encodes a member of the pentatricopeptide protein family, involved in processing of RNA in plastids, where carotenoid and chlorophyll pigments accumulate. Network analysis of > 12 000 eQTL mapped for > 8000 differentially expressed fruit genes supported the role of CmPPR1 in determining the expression level of plastid targeted genes. We highlight the potential of RNA-Seq-based QTL analysis of small to moderate size, advanced RIL populations for precise marker-assisted breeding and gene discovery. We provide the following resources: a RIL population genotyped with a unique set of SNP markers, confined genomic segments that harbor QTL governing 129 traits and a saturated set of melon eQTLs.

Genome-Wide Linkage-Disequilibrium Mapping to the Candidate Gene Level in Melon (Cucumis melo).

Cucumis melo is highly diverse for fruit traits providing wide breeding and genetic research opportunities, including genome-wide association (GWA) analysis. We used a collection of 177 accessions representing the two C. melo subspecies and 11 horticultural groups for detailed characterization of fruit traits variation and evaluation of the potential of GWA for trait mapping in melon. Through genotyping-by-sequencing, 23,931 informative SNPs were selected for genome-wide analyses. We found that linkage-disequilibrium decays at ~100 Kb in this collection and that population structure effect on association results varies between traits. We mapped several monogenic traits to narrow intervals overlapping with known causative genes, demonstrating the potential of diverse collections and GWA for mapping Mendelian traits to a candidate-gene level in melon. We further report on mapping of fruit shape quantitative trait loci (QTLs) and comparison with multiple previous QTL studies. Expansion of sample size and a more balanced representation of taxonomic groups might improve efficiency for simple traits dissection. But, as in other plant species, integrated linkage-association multi-allelic approaches are likely to produce better combination of statistical power, diversity capture and mapping resolution in melon. Our data can be utilized for selection of the most appropriate accessions for such approaches.