A physical and genetic map of Cannabis sativa identifies extensive rearrangements at the THC/CBD acid synthase loci.

Cannabis sativa is widely cultivated for medicinal, food, industrial, and recreational use, but much remains unknown regarding its genetics, including the molecular determinants of cannabinoid content. Here, we describe a combined physical and genetic map derived from a cross between the drug-type strain Purple Kush and the hemp variety "Finola." The map reveals that cannabinoid biosynthesis genes are generally unlinked but that aromatic prenyltransferase (AP), which produces the substrate for THCA and CBDA synthases (THCAS and CBDAS), is tightly linked to a known marker for total cannabinoid content. We further identify the gene encoding CBCA synthase (CBCAS) and characterize its catalytic activity, providing insight into how cannabinoid diversity arises in cannabis. THCAS and CBDAS (which determine the drug vs. hemp chemotype) are contained within large (>250 kb) retrotransposon-rich regions that are highly nonhomologous between drug- and hemp-type alleles and are furthermore embedded within ∼40 Mb of minimally recombining repetitive DNA. The chromosome structures are similar to those in grains such as wheat, with recombination focused in gene-rich, repeat-depleted regions near chromosome ends. The physical and genetic map should facilitate further dissection of genetic and molecular mechanisms in this commercially and medically important plant.

Mediator Complex Subunits MED2, MED5, MED16, and MED23 Genetically Interact in the Regulation of Phenylpropanoid Biosynthesis.

The phenylpropanoid pathway is a major global carbon sink and is important for plant fitness and the engineering of bioenergy feedstocks. In Arabidopsis thaliana, disruption of two subunits of the transcriptional regulatory Mediator complex, MED5a and MED5b, results in an increase in phenylpropanoid accumulation. By contrast, the semidominant MED5b mutation reduced epidermal fluorescence4-3 (ref4-3) results in dwarfism and constitutively repressed phenylpropanoid accumulation. Here, we report the results of a forward genetic screen for suppressors of ref4-3. We identified 13 independent lines that restore growth and/or phenylpropanoid accumulation in the ref4-3 background. Two of the suppressors restore growth without restoring soluble phenylpropanoid accumulation, indicating that the growth and metabolic phenotypes of the ref4-3 mutant can be genetically disentangled. Whole-genome sequencing revealed that all but one of the suppressors carry mutations in MED5b or other Mediator subunits. RNA-seq analysis showed that the ref4-3 mutation causes widespread changes in gene expression, including the upregulation of negative regulators of the phenylpropanoid pathway, and that the suppressors reverse many of these changes. Together, our data highlight the interdependence of individual Mediator subunits and provide greater insight into the transcriptional regulation of phenylpropanoid biosynthesis by the Mediator complex.

Identification of olivetolic acid cyclase from Cannabis sativa reveals a unique catalytic route to plant polyketides.

Δ(9)-Tetrahydrocannabinol (THC) and other cannabinoids are responsible for the psychoactive and medicinal properties of Cannabis sativa L. (marijuana). The first intermediate in the cannabinoid biosynthetic pathway is proposed to be olivetolic acid (OA), an alkylresorcinolic acid that forms the polyketide nucleus of the cannabinoids. OA has been postulated to be synthesized by a type III polyketide synthase (PKS) enzyme, but so far type III PKSs from cannabis have been shown to produce catalytic byproducts instead of OA. We analyzed the transcriptome of glandular trichomes from female cannabis flowers, which are the primary site of cannabinoid biosynthesis, and searched for polyketide cyclase-like enzymes that could assist in OA cyclization. Here, we show that a type III PKS (tetraketide synthase) from cannabis trichomes requires the presence of a polyketide cyclase enzyme, olivetolic acid cyclase (OAC), which catalyzes a C2-C7 intramolecular aldol condensation with carboxylate retention to form OA. OAC is a dimeric α+ß barrel (DABB) protein that is structurally similar to polyketide cyclases from Streptomyces species. OAC transcript is present at high levels in glandular trichomes, an expression profile that parallels other cannabinoid pathway enzymes. Our identification of OAC both clarifies the cannabinoid pathway and demonstrates unexpected evolutionary parallels between polyketide biosynthesis in plants and bacteria. In addition, the widespread occurrence of DABB proteins in plants suggests that polyketide cyclases may play an overlooked role in generating plant chemical diversity.

The hexanoyl-CoA precursor for cannabinoid biosynthesis is formed by an acyl-activating enzyme in Cannabis sativa trichomes.

The psychoactive and analgesic cannabinoids (e.g. Δ(9) -tetrahydrocannabinol (THC)) in Cannabis sativa are formed from the short-chain fatty acyl-coenzyme A (CoA) precursor hexanoyl-CoA. Cannabinoids are synthesized in glandular trichomes present mainly on female flowers. We quantified hexanoyl-CoA using LC-MS/MS and found levels of 15.5 pmol g(-1) fresh weight in female hemp flowers with lower amounts in leaves, stems and roots. This pattern parallels the accumulation of the end-product cannabinoid, cannabidiolic acid (CBDA). To search for the acyl-activating enzyme (AAE) that synthesizes hexanoyl-CoA from hexanoate, we analyzed the transcriptome of isolated glandular trichomes. We identified 11 unigenes that encoded putative AAEs including CsAAE1, which shows high transcript abundance in glandular trichomes. In vitro assays showed that recombinant CsAAE1 activates hexanoate and other short- and medium-chained fatty acids. This activity and the trichome-specific expression of CsAAE1 suggest that it is the hexanoyl-CoA synthetase that supplies the cannabinoid pathway. CsAAE3 encodes a peroxisomal enzyme that activates a variety of fatty acid substrates including hexanoate. Although phylogenetic analysis showed that CsAAE1 groups with peroxisomal AAEs, it lacked a peroxisome targeting sequence 1 (PTS1) and localized to the cytoplasm. We suggest that CsAAE1 may have been recruited to the cannabinoid pathway through the loss of its PTS1, thereby redirecting it to the cytoplasm. To probe the origin of hexanoate, we analyzed the trichome expressed sequence tag (EST) dataset for enzymes of fatty acid metabolism. The high abundance of transcripts that encode desaturases and a lipoxygenase suggests that hexanoate may be formed through a pathway that involves the oxygenation and breakdown of unsaturated fatty acids.

Density-Dependent Prophylaxis in Freshwater Snails Driven by Oxylipin Chemical Cues.

While animal aggregations can benefit the fitness of group members, the behaviour may also lead to higher risks of parasite infection as group density increases. Some animals are known to moderate their investment in immunity relative to the risk of infection. These animals exhibit density-dependent prophylaxis (DDP) by increasing their immune investment as group density increases. Despite being documented in many taxa, the mechanisms of DDP remain largely unexplored. Snails are known to aggregate and experience large fluctuations in density and serve as required hosts for many parasites. Further, they are known to use chemical cues to aggregate. To test whether freshwater snails exhibit DDP and investigate the role that chemical signaling compounds may play in triggering this phenomenon, we performed four experiments on the freshwater snail Stagnicola elodes, which is a common host for many trematode parasite species. First, we tested if DDP occurred in snails in laboratory-controlled conditions (control vs snail-conditioned water) and whether differences in exposure to chemical cues affected immune function. Second, we used gas chromatography to characterize fatty acids expressed in snail-conditioned water to determine if precursors for particular signaling molecules, such as oxylipins, were being produced by snails. Third, we characterized the oxylipins released by infected and uninfected field-collected snails, to better understand how differences in oxylipin cocktails may play a role in inducing DDP. Finally, we tested the immune response of snails exposed to four oxylipins to test the ability of specific oxylipins to affect DDP. We found that snails exposed to water with higher densities of snails and raised in snail-conditioned water had higher counts of haemocytes. Additionally, lipid analysis demonstrated that fatty acid molecules that are also precursors for oxylipins were present in snail-conditioned water. Trematode-infected snails emitted 50 oxylipins in higher amounts, with 24 of these oxylipins only detected in this group. Finally, oxylipins that were higher in infected snails induced naïve snails to increase their immune responses compared to sham-exposed snails. Our results provide evidence that snails exhibit DDP, and the changes in oxylipins emitted by infected hosts may be one of the molecular mechanisms driving this phenomenon.

Separation of isomeric short-chain acyl-CoAs in plant matrices using ultra-performance liquid chromatography coupled with tandem mass spectrometry.

Acyl coenzyme A (acyl-CoA) thioesters are important intermediates in cellular metabolism and being able to distinguish among them is critical to fully understanding metabolic pathways in plants. Although significant advances have been made in the identification and quantification of acyl-CoAs using liquid chromatography tandem mass spectrometry (LC-MS/MS), separation of isomeric species such as isobutyryl- and n-butyrl-CoA has remained elusive. Here we report an ultra-performance liquid chromatography (UPLC)-MS/MS method for quantifying short-chain acyl-CoAs including isomeric species n-butyryl-CoA and isobutyryl-CoA as well as n-valeryl-CoA and isovaleryl-CoA. The method was applied to the analysis of extracts of hop (Humulus lupulus) and provided strong evidence for the existence of an additional structural isomer of valeryl-CoA, 2-methylbutyryl-CoA, as well as an unexpected isomer of hexanoyl-CoA. The results showed differences in the acyl-CoA composition among varieties of Humulus lupulus, both in glandular trichomes and cone tissues. When compared with the analysis of hemp (Cannabis sativa) extracts, the contribution of isobutyryl-CoAs in hop was greater as would be expected based on the downstream polyketide products. Surprisingly, branched chain valeryl-CoAs (isovaleryl-CoA and 2-methylbutyryl-CoA) were the dominant form of valeryl-CoAs in both hop and hemp. The capability to separate these isomeric forms will help to understand biochemical pathways leading to specialized metabolites in plants.

The Genetic Structure of Marijuana and Hemp.

Despite its cultivation as a source of food, fibre and medicine, and its global status as the most used illicit drug, the genus Cannabis has an inconclusive taxonomic organization and evolutionary history. Drug types of Cannabis (marijuana), which contain high amounts of the psychoactive cannabinoid Δ9-tetrahydrocannabinol (THC), are used for medical purposes and as a recreational drug. Hemp types are grown for the production of seed and fibre, and contain low amounts of THC. Two species or gene pools (C. sativa and C. indica) are widely used in describing the pedigree or appearance of cultivated Cannabis plants. Using 14,031 single-nucleotide polymorphisms (SNPs) genotyped in 81 marijuana and 43 hemp samples, we show that marijuana and hemp are significantly differentiated at a genome-wide level, demonstrating that the distinction between these populations is not limited to genes underlying THC production. We find a moderate correlation between the genetic structure of marijuana strains and their reported C. sativa and C. indica ancestry and show that marijuana strain names often do not reflect a meaningful genetic identity. We also provide evidence that hemp is genetically more similar to C. indica type marijuana than to C. sativa strains.

Transcriptome analysis based on next-generation sequencing of non-model plants producing specialized metabolites of biotechnological interest.

Plants produce a vast array of specialized metabolites, many of which are used as pharmaceuticals, flavors, fragrances, and other high-value fine chemicals. However, most of these compounds occur in non-model plants for which genomic sequence information is not yet available. The production of a large amount of nucleotide sequence data using next-generation technologies is now relatively fast and cost-effective, especially when using the latest Roche-454 and Illumina sequencers with enhanced base-calling accuracy. To investigate specialized metabolite biosynthesis in non-model plants we have established a data-mining framework, employing next-generation sequencing and computational algorithms, to construct and analyze the transcriptomes of 75 non-model plants that produce compounds of interest for biotechnological applications. After sequence assembly an extensive annotation approach was applied to assign functional information to over 800,000 putative transcripts. The annotation is based on direct searches against public databases, including RefSeq and InterPro. Gene Ontology (GO), Enzyme Commission (EC) annotations and associated Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway maps are also collected. As a proof-of-concept, the selection of biosynthetic gene candidates associated with six specialized metabolic pathways is described. A web-based BLAST server has been established to allow public access to assembled transcriptome databases for all 75 plant species of the PhytoMetaSyn Project (www.phytometasyn.ca).

The draft genome and transcriptome of Cannabis sativa.

BACKGROUND: Cannabis sativa has been cultivated throughout human history as a source of fiber, oil and food, and for its medicinal and intoxicating properties. Selective breeding has produced cannabis plants for specific uses, including high-potency marijuana strains and hemp cultivars for fiber and seed production. The molecular biology underlying cannabinoid biosynthesis and other traits of interest is largely unexplored. RESULTS: We sequenced genomic DNA and RNA from the marijuana strain Purple Kush using shortread approaches. We report a draft haploid genome sequence of 534 Mb and a transcriptome of 30,000 genes. Comparison of the transcriptome of Purple Kush with that of the hemp cultivar 'Finola' revealed that many genes encoding proteins involved in cannabinoid and precursor pathways are more highly expressed in Purple Kush than in 'Finola'. The exclusive occurrence of Δ9-tetrahydrocannabinolic acid synthase in the Purple Kush transcriptome, and its replacement by cannabidiolic acid synthase in 'Finola', may explain why the psychoactive cannabinoid Δ9-tetrahydrocannabinol (THC) is produced in marijuana but not in hemp. Resequencing the hemp cultivars 'Finola' and 'USO-31' showed little difference in gene copy numbers of cannabinoid pathway enzymes. However, single nucleotide variant analysis uncovered a relatively high level of variation among four cannabis types, and supported a separation of marijuana and hemp. CONCLUSIONS: The availability of the Cannabis sativa genome enables the study of a multifunctional plant that occupies a unique role in human culture. Its availability will aid the development of therapeutic marijuana strains with tailored cannabinoid profiles and provide a basis for the breeding of hemp with improved agronomic characteristics.