Characterization of the ZCTs, a subgroup of Cys2-His2 zinc finger transcription factors regulating alkaloid biosynthesis in Catharanthus roseus.

KEY MESSAGE: The C. roseus ZCTs are jasmonate-responsive, can be induced by CrMYC2a, and can act as significant regulators of the terpenoid indole alkaloid pathway when highly expressed. Catharanthus roseus is the sole known producer of the anti-cancer terpenoid indole alkaloids (TIAs), vinblastine and vincristine. While the enzymatic steps of the pathway have been elucidated, an understanding of its regulation is still emerging. The present study characterizes an important subgroup of Cys2-His2 zinc finger transcription factors known as Zinc finger Catharanthus Transcription factors (ZCTs). We identified three new ZCT members (named ZCT4, ZCT5, and ZCT6) that clustered with the putative repressors of the TIA pathway, ZCT1, ZCT2, and ZCT3. We characterized the role of these six ZCTs as potential redundant regulators of the TIA pathway, and their tissue-specific and jasmonate-responsive expression. These ZCTs share high sequence conservation in their two Cys2-His2 zinc finger domains but differ in the spacer length and sequence between these zinc fingers. The transient overexpression of ZCTs in seedlings significantly repressed the promoters of the terpenoid (pLAMT) and condensation branch (pSTR1) of the TIA pathway, consistent with that previously reported for ZCT1, ZCT2, and ZCT3. In addition, ZCTs significantly repressed and indirectly activated several promoters of the vindoline pathway (not previously studied). The ZCTs differed in their tissue-specific expression but similarly increased with jasmonate in a dosage-dependent manner (except for ZCT5). We showed significant activation of the pZCT1 and pZCT3 promoters by the de-repressed CrMYC2a, suggesting that the jasmonate-responsive expression of the ZCTs can be mediated by CrMYC2a. In summary, the C. roseus ZCTs are jasmonate-responsive, can be induced by CrMYC2a, and can act as significant regulators of the TIA pathway when highly expressed.

Transcriptomics Identifies Modules of Differentially Expressed Genes and Novel Cyclotides in Viola pubescens.

Viola is a large genus with worldwide distribution and many traits not currently exemplified in model plants including unique breeding systems and the production of cyclotides. Here we report de novo genome assembly and transcriptomic analyses of the non-model species Viola pubescens using short-read DNA sequencing data and RNA-Seq from eight diverse tissues. First, V. pubescens genome size was estimated through flow cytometry, resulting in an approximate haploid genome of 455 Mbp. Next, the draft V. pubescens genome was sequenced and assembled resulting in 264,035,065 read pairs and 161,038 contigs with an N50 length of 3,455 base pairs (bp). RNA-Seq data were then assembled into tissue-specific transcripts. Together, the DNA and transcript data generated 38,081 ab initio gene models which were functionally annotated based on homology to Arabidopsis thaliana genes and Pfam domains. Gene expression was visualized for each tissue via principal component analysis and hierarchical clustering, and gene co-expression analysis identified 20 modules of tissue-specific transcriptional networks. Some of these modules highlight genetic differences between chasmogamous and cleistogamous flowers and may provide insight into V. pubescens' mixed breeding system. Orthologous clustering with the proteomes of A. thaliana and Populus trichocarpa revealed 8,531 sequences unique to V. pubescens, including 81 novel cyclotide precursor sequences. Cyclotides are plant peptides characterized by a stable, cyclic cystine knot motif, making them strong candidates for drug scaffolding and protein engineering. Analysis of the RNA-Seq data for these cyclotide transcripts revealed diverse expression patterns both between transcripts and tissues. The diversity of these cyclotides was also highlighted in a maximum likelihood protein cladogram containing V. pubescens cyclotides and published cyclotide sequences from other Violaceae and Rubiaceae species. Collectively, this work provides the most comprehensive sequence resource for Viola, offers valuable transcriptomic insight into V. pubescens, and will facilitate future functional genomics research in Viola and other diverse plant groups.

Origin and evolution of the octoploid strawberry genome.

Cultivated strawberry emerged from the hybridization of two wild octoploid species, both descendants from the merger of four diploid progenitor species into a single nucleus more than 1 million years ago. Here we report a near-complete chromosome-scale assembly for cultivated octoploid strawberry (Fragaria × ananassa) and uncovered the origin and evolutionary processes that shaped this complex allopolyploid. We identified the extant relatives of each diploid progenitor species and provide support for the North American origin of octoploid strawberry. We examined the dynamics among the four subgenomes in octoploid strawberry and uncovered the presence of a single dominant subgenome with significantly greater gene content, gene expression abundance, and biased exchanges between homoeologous chromosomes, as compared with the other subgenomes. Pathway analysis showed that certain metabolomic and disease-resistance traits are largely controlled by the dominant subgenome. These findings and the reference genome should serve as a powerful platform for future evolutionary studies and enable molecular breeding in strawberry.

Haplotype-phased genome and evolution of phytonutrient pathways of tetraploid blueberry.

BACKGROUND: Highbush blueberry (Vaccinium corymbosum) has long been consumed for its unique flavor and composition of health-promoting phytonutrients. However, breeding efforts to improve fruit quality in blueberry have been greatly hampered by the lack of adequate genomic resources and a limited understanding of the underlying genetics encoding key traits. The genome of highbush blueberry has been particularly challenging to assemble due, in large part, to its polyploid nature and genome size. FINDINGS: Here, we present a chromosome-scale and haplotype-phased genome assembly of the cultivar "Draper," which has the highest antioxidant levels among a diversity panel of 71 cultivars and 13 wild Vaccinium species. We leveraged this genome, combined with gene expression and metabolite data measured across fruit development, to identify candidate genes involved in the biosynthesis of important phytonutrients among other metabolites associated with superior fruit quality. Genome-wide analyses revealed that both polyploidy and tandem gene duplications modified various pathways involved in the biosynthesis of key phytonutrients. Furthermore, gene expression analyses hint at the presence of a spatial-temporal specific dominantly expressed subgenome including during fruit development. CONCLUSIONS: These findings and the reference genome will serve as a valuable resource to guide future genome-enabled breeding of important agronomic traits in highbush blueberry.

Isoprene Acts as a Signaling Molecule in Gene Networks Important for Stress Responses and Plant Growth.

Isoprene synthase converts dimethylallyl diphosphate to isoprene and appears to be necessary and sufficient to allow plants to emit isoprene at significant rates. Isoprene can protect plants from abiotic stress but is not produced naturally by all plants; for example, Arabidopsis (Arabidopsis thaliana) and tobacco (Nicotiana tabacum) do not produce isoprene. It is typically present at very low concentrations, suggesting a role as a signaling molecule; however, its exact physiological role and mechanism of action are not fully understood. We transformed Arabidopsis with a Eucalyptus globulus isoprene synthase The regulatory mechanisms of photosynthesis and isoprene emission were similar to those of native emitters, indicating that regulation of isoprene emission is not specific to isoprene-emitting species. Leaf chlorophyll and carotenoid contents were enhanced by isoprene, which also had a marked positive effect on hypocotyl, cotyledon, leaf, and inflorescence growth in Arabidopsis. By contrast, leaf and stem growth was reduced in tobacco engineered to emit isoprene. Expression of genes belonging to signaling networks or associated with specific growth regulators (e.g. gibberellic acid that promotes growth and jasmonic acid that promotes defense) and genes that lead to stress tolerance was altered by isoprene emission. Isoprene likely executes its effects on growth and stress tolerance through direct regulation of gene expression. Enhancement of jasmonic acid-mediated defense signaling by isoprene may trigger a growth-defense tradeoff leading to variations in the growth response. Our data support a role for isoprene as a signaling molecule.

Subgenome Dominance in an Interspecific Hybrid, Synthetic Allopolyploid, and a 140-Year-Old Naturally Established Neo-Allopolyploid Monkeyflower.

Recent studies have shown that one of the parental subgenomes in ancient polyploids is generally more dominant, having retained more genes and being more highly expressed, a phenomenon termed subgenome dominance. The genomic features that determine how quickly and which subgenome dominates within a newly formed polyploid remain poorly understood. To investigate the rate of emergence of subgenome dominance, we examined gene expression, gene methylation, and transposable element (TE) methylation in a natural, <140-year-old allopolyploid (Mimulus peregrinus), a resynthesized interspecies triploid hybrid (M. robertsii), a resynthesized allopolyploid (M. peregrinus), and progenitor species (M. guttatus and M. luteus). We show that subgenome expression dominance occurs instantly following the hybridization of divergent genomes and significantly increases over generations. Additionally, CHH methylation levels are reduced in regions near genes and within TEs in the first-generation hybrid, intermediate in the resynthesized allopolyploid, and are repatterned differently between the dominant and recessive subgenomes in the natural allopolyploid. Subgenome differences in levels of TE methylation mirror the increase in expression bias observed over the generations following hybridization. These findings provide important insights into genomic and epigenomic shock that occurs following hybridization and polyploid events and may also contribute to uncovering the mechanistic basis of heterosis and subgenome dominance.

Genome-guided investigation of plant natural product biosynthesis.

The medicinal plant Madagascar periwinkle, Catharanthus roseus (L.) G. Don, produces hundreds of biologically active monoterpene-derived indole alkaloid (MIA) metabolites and is the sole source of the potent, expensive anti-cancer compounds vinblastine and vincristine. Access to a genome sequence would enable insights into the biochemistry, control, and evolution of genes responsible for MIA biosynthesis. However, generation of a near-complete, scaffolded genome is prohibitive to small research communities due to the expense, time, and expertise required. In this study, we generated a genome assembly for C. roseus that provides a near-comprehensive representation of the genic space that revealed the genomic context of key points within the MIA biosynthetic pathway including physically clustered genes, tandem gene duplication, expression sub-functionalization, and putative neo-functionalization. The genome sequence also facilitated high resolution co-expression analyses that revealed three distinct clusters of co-expression within the components of the MIA pathway. Coordinated biosynthesis of precursors and intermediates throughout the pathway appear to be a feature of vinblastine/vincristine biosynthesis. The C. roseus genome also revealed localization of enzyme-rich genic regions and transporters near known biosynthetic enzymes, highlighting how even a draft genome sequence can empower the study of high-value specialized metabolites.

Nucleotide polymorphism and copy number variant detection using exome capture and next-generation sequencing in the polyploid grass Panicum virgatum.

Switchgrass (Panicum virgatum) is a polyploid, outcrossing grass species native to North America and has recently been recognized as a potential biofuel feedstock crop. Significant phenotypic variation including ploidy is present across the two primary ecotypes of switchgrass, referred to as upland and lowland switchgrass. The tetraploid switchgrass genome is approximately 1400 Mbp, split between two subgenomes, with significant repetitive sequence content limiting the efficiency of re-sequencing approaches for determining genome diversity. To characterize genetic diversity in upland and lowland switchgrass as a first step in linking genotype to phenotype, we designed an exome capture probe set based on transcript assemblies that represent approximately 50 Mb of annotated switchgrass exome sequences. We then evaluated and optimized the probe set using solid phase comparative genome hybridization and liquid phase exome capture followed by next-generation sequencing. Using the optimized probe set, we assessed variation in the exomes of eight switchgrass genotypes representing tetraploid lowland and octoploid upland cultivars to benchmark our exome capture probe set design. We identified ample variation in the switchgrass genome including 1,395,501 single nucleotide polymorphisms (SNPs), 8173 putative copy number variants and 3336 presence/absence variants. While the majority of the SNPs (84%) detected was bi-allelic, a substantial number was tri-allelic with limited occurrence of tetra-allelic polymorphisms consistent with the heterozygous and polyploid nature of the switchgrass genome. Collectively, these data demonstrate the efficacy of exome capture for discovery of genome variation in a polyploid species with a large, repetitive and heterozygous genome.

Transcriptional and metabolic analysis of senescence induced by preventing pollination in maize.

Transcriptional and metabolic changes were evaluated during senescence induced by preventing pollination in the B73 genotype of maize (Zea mays). Accumulation of free glucose and starch and loss of chlorophyll in leaf was manifested early at 12 d after anthesis (DAA), while global transcriptional and phenotypic changes were evident only at 24 DAA. Internodes exhibited major transcriptomic changes only at 30 DAA. Overlaying expression data onto metabolic pathways revealed involvement of many novel pathways, including those involved in cell wall biosynthesis. To investigate the overlap between induced and natural senescence, transcriptional data from induced senescence in maize was compared with that reported for Arabidopsis (Arabidopsis thaliana) undergoing natural and sugar-induced senescence. Notable similarities with natural senescence in Arabidopsis included up-regulation of senescence-associated genes (SAGs), ethylene and jasmonic acid biosynthetic genes, APETALA2, ethylene-responsive element binding protein, and no apical meristem transcription factors. However, differences from natural senescence were highlighted by unaltered expression of a subset of the SAGs, and cytokinin, abscisic acid, and salicylic acid biosynthesis genes. Key genes up-regulated during sugar-induced senescence in Arabidopsis, including a cysteine protease (SAG12) and three flavonoid biosynthesis genes (PRODUCTION OF ANTHOCYANIN PIGMENT1 (PAP1), PAP2, and LEUCOANTHOCYANIDIN DIOXYGENASE), were also induced, suggesting similarities in senescence induced by pollination prevention and sugar application. Coexpression analysis revealed networks involving known senescence-related genes and novel candidates; 82 of these were shared between leaf and internode networks, highlighting similarities in induced senescence in these tissues. Insights from this study will be valuable in systems biology of senescence in maize and other grasses.

Development of transcriptomic resources for interrogating the biosynthesis of monoterpene indole alkaloids in medicinal plant species.

The natural diversity of plant metabolism has long been a source for human medicines. One group of plant-derived compounds, the monoterpene indole alkaloids (MIAs), includes well-documented therapeutic agents used in the treatment of cancer (vinblastine, vincristine, camptothecin), hypertension (reserpine, ajmalicine), malaria (quinine), and as analgesics (7-hydroxymitragynine). Our understanding of the biochemical pathways that synthesize these commercially relevant compounds is incomplete due in part to a lack of molecular, genetic, and genomic resources for the identification of the genes involved in these specialized metabolic pathways. To address these limitations, we generated large-scale transcriptome sequence and expression profiles for three species of Asterids that produce medicinally important MIAs: Camptotheca acuminata, Catharanthus roseus, and Rauvolfia serpentina. Using next generation sequencing technology, we sampled the transcriptomes of these species across a diverse set of developmental tissues, and in the case of C. roseus, in cultured cells and roots following elicitor treatment. Through an iterative assembly process, we generated robust transcriptome assemblies for all three species with a substantial number of the assembled transcripts being full or near-full length. The majority of transcripts had a related sequence in either UniRef100, the Arabidopsis thaliana predicted proteome, or the Pfam protein domain database; however, we also identified transcripts that lacked similarity with entries in either database and thereby lack a known function. Representation of known genes within the MIA biosynthetic pathway was robust. As a diverse set of tissues and treatments were surveyed, expression abundances of transcripts in the three species could be estimated to reveal transcripts associated with development and response to elicitor treatment. Together, these transcriptomes and expression abundance matrices provide a rich resource for understanding plant specialized metabolism, and promotes realization of innovative production systems for plant-derived pharmaceuticals.