Elsa Góngora-Castillo

Architecture and evolution of a minute plant genome

It has been argued that the evolution of plant genome size is principally unidirectional and increasing owing to the varied action of whole-genome duplications (WGDs) and mobile element proliferation. However, extreme genome size reductions have been reported in the angiosperm family tree. Here we report the sequence of the 82-megabase genome of the carnivorous bladderwort plant Utricularia gibba. Despite its tiny size, the U. gibba genome accommodates a typical number of genes for a plant, with the main difference from other plant genomes arising from a drastic reduction in non-genic DNA. Unexpectedly, we identified at least three rounds of WGD in U. gibba since common ancestry with tomato (Solanum) and grape (Vitis). The compressed architecture of the U. gibba genome indicates that a small fraction of intergenic DNA, with few or no active retrotransposons, is sufficient to regulate and integrate all the processes required for the development and reproduction of a complex organism.

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.

Comparative Transcriptome Atlases Reveal Altered Gene Expression Modules between Two Cleomaceae C3 and C4 Plant Species

Transcriptome atlases of Tarenaya hassleriana, a C3 plant, and the C4-species Gynandropsis gynandra were developed. Integrated analysis of RNA-seq and leaf anatomical data indicates altered regulation of the cell cycle and endoreduplication during the establishment of C4 leaf anatomy. C4 genes have been recruited to photosynthetic gene expression regulation from ancestral nonphotosynthetic functions in C3. C4 photosynthesis outperforms the ancestral C3 state in a wide range of natural and agro-ecosystems by affording higher water-use and nitrogen-use efficiencies. It therefore represents a prime target for engineering novel, high-yielding crops by introducing the trait into C3 backgrounds. However, the genetic architecture of C4 photosynthesis remains largely unknown. To define the divergence in gene expression modules between C3 and C4 photosynthesis during leaf ontogeny, we generated comprehensive transcriptome atlases of two Cleomaceae species, Gynandropsis gynandra (C4) and Tarenaya hassleriana (C3), by RNA sequencing. Overall, the gene expression profiles appear remarkably similar between the C3 and C4 species. We found that known C4 genes were recruited to photosynthesis from different expression domains in C3, including typical housekeeping gene expression patterns in various tissues as well as individual heterotrophic tissues. Furthermore, we identified a structure-related module recruited from the C3 root. Comparison of gene expression patterns with anatomy during leaf ontogeny provided insight into genetic features of Kranz anatomy. Altered expression of developmental factors and cell cycle genes is associated with a higher degree of endoreduplication in enlarged C4 bundle sheath cells. A delay in mesophyll differentiation apparent both in the leaf anatomy and the transcriptome allows for extended vein formation in the C4 leaf.

Coupling Deep Transcriptome Analysis with Untargeted Metabolic Profiling in Ophiorrhiza pumila to Further the Understanding of the Biosynthesis of the Anti-cancer Alkaloid Camptothecin and Anthraquinones

The Rubiaceae species, Ophiorrhiza pumila, accumulates camptothecin, an anti-cancer alkaloid with a potent DNA topoisomerase I inhibitory activity, as well as anthraquinones that are derived from the combination of the isochorismate and hemiterpenoid pathways. The biosynthesis of these secondary products is active in O. pumila hairy roots yet very low in cell suspension culture. Deep transcriptome analysis was conducted in O. pumila hairy roots and cell suspension cultures using the Illumina platform, yielding a total of 2 Gb of sequence for each sample. We generated a hybrid transcriptome assembly of O. pumila using the Illumina-derived short read sequences and conventional Sanger-derived expressed sequence tag clones derived from a full-length cDNA library constructed using RNA from hairy roots. Among 35,608 non-redundant unigenes, 3,649 were preferentially expressed in hairy roots compared with cell suspension culture. Candidate genes involved in the biosynthetic pathway for the monoterpenoid indole alkaloid camptothecin were identified; specifically, genes involved in post-strictosamide biosynthetic events and genes involved in the biosynthesis of anthraquinones and chlorogenic acid. Untargeted metabolomic analysis by Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR-MS) indicated that most of the proposed intermediates in the camptothecin biosynthetic pathway accumulated in hairy roots in a preferential manner compared with cell suspension culture. In addition, a number of anthraquinones and chlorogenic acid preferentially accumulated in hairy roots compared with cell suspension culture. These results suggest that deep transcriptome and metabolome data sets can facilitate the identification of genes and intermediates involved in the biosynthesis of secondary products including camptothecin in O. pumila.

Bioinformatics challenges in de novo transcriptome assembly using short read sequences in the absence of a reference genome sequence

Plant natural product research can be facilitated through genome and transcriptome sequencing approaches that generate informative sequence and expression datasets that enable characterization of biochemical pathways of interest. As the overwhelming majority of plant-derived natural products are derived from species with little, if any, sequence and/or genomic resources, the ability to perform whole genome shotgun sequencing and assembly has been and will continue to be transformative as access to a genome sequence provides molecular resources and a context for discovery and characterization of biosynthetic pathways. Due to the reduced size and complexity of the transcriptome relative to the genome, transcriptome sequencing provides a rapid, inexpensive approach to access gene sequences, gene expression abundances, and gene expression patterns in any species, including those that lack a reference genome sequence. To date, successful applications of RNA sequencing in conjunction with de novo transcriptome assembly has enabled identification of new genes in an array of biochemical pathways in plants. While sequencing technologies are well developed, challenges remain in the handling and analysis of transcriptome sequences. In this Highlight article, we provide an overview of the bioinformatics challenges associated with transcriptome analyses using short read sequences and how to address these issues in plant species that lack a reference genome.

Genomic approaches for interrogating the biochemistry of medicinal plant species.

Development of next-generation sequencing, coupled with the advancement of computational methods, has allowed researchers to access the transcriptomes of recalcitrant genomes such as those of medicinal plant species. Through the sequencing of even a few cDNA libraries, a broad representation of the transcriptome of any medicinal plant species can be obtained, providing a robust resource for gene discovery and downstream biochemical pathway discovery. When coupled to estimation of expression abundances in specific tissues from a developmental series, biotic stress, abiotic stress, or elicitor challenge, informative coexpression and differential expression estimates on a whole transcriptome level can be obtained to identify candidates for function discovery.

Transcriptome analysis of symptomatic and recovered leaves of geminivirus-infected pepper ( Capsicum annuum )

BackgroundGeminiviruses are a large and important family of plant viruses that infect a wide range of crops throughout the world. The Begomovirus genus contains species that are transmitted by whiteflies and are distributed worldwide causing disease on an array of horticultural crops. Symptom remission, in which newly developed leaves of systemically infected plants exhibit a reduction in symptom severity (recovery), has been observed on pepper (Capsicum annuum) plants infected with Pepper golden mosaic virus (PepGMV). Previous studies have shown that transcriptional and post-transcriptional gene silencing mechanisms are involved in the reduction of viral nucleic acid concentration in recovered tissue. In this study, we employed deep transcriptome sequencing methods to assess transcriptional variation in healthy (mock), symptomatic, and recovered pepper leaves following PepGMV infection.ResultsDifferential expression analyses of the pepper leaf transcriptome from symptomatic and recovered stages revealed a total of 309 differentially expressed genes between healthy (mock) and symptomatic or recovered tissues. Computational prediction of differential expression was validated using quantitative reverse-transcription PCR confirming the robustness of our bioinformatic methods. Within the set of differentially expressed genes associated with the recovery process were genes involved in defense responses including pathogenesis-related proteins, reactive oxygen species, systemic acquired resistance, jasmonic acid biosynthesis, and ethylene signaling. No major differences were found when compared the differentially expressed genes in symptomatic and recovered tissues. On the other hand, a set of genes with novel roles in defense responses was identified including genes involved in histone modification. This latter result suggested that post-transcriptional and transcriptional gene silencing may be one of the major mechanisms involved in the recovery process. Genes orthologous to the C. annuum proteins involved in the pepper-PepGMV recovery response were identified in both Solanum lycopersicum and Solanum tuberosum suggesting conservation of components of the viral recovery response in the Solanaceae.ConclusionThese data provide a valuable source of information for improving our understanding of the underlying molecular mechanisms by which pepper leaves become symptomless following infection with geminiviruses. The identification of orthologs for the majority of genes differentially expressed in recovered tissues in two major solanaceous crop species provides the basis for future comparative analyses of the viral recovery process across related taxa.

Functional identification of valerena-1,10-diene synthase, a terpene synthase catalyzing a unique chemical cascade in the biosynthesis of biologically active sesquiterpenes in Valeriana officinalis

Background: Therapeutic values of Valeriana officinalis have been associated with sesquiterpenes whose biosynthetic origins have remained enigmatic. Results: A cyclobutenyl intermediate in the catalytic cascade of valerena-1,10-diene synthase is reported. Conclusion: A new class of sesquiterpene synthases for the biosynthesis of sesquiterpenes harboring isobutenyl functional groups is proposed. Significance: Similar catalytic mechanisms from evolutionarily diverse organisms are proposed and portend sources for sesquiterpene diversity. Valerian is an herbal preparation from the roots of Valeriana officinalis used as an anxiolytic and sedative and in the treatment of insomnia. The biological activities of valerian are attributed to valerenic acid and its putative biosynthetic precursor valerenadiene, sesquiterpenes, found in V. officinalis roots. These sesquiterpenes retain an isobutenyl side chain whose origin has been long recognized as enigmatic because a chemical rationalization for their biosynthesis has not been obvious. Using recently developed metabolomic and transcriptomic resources, we identified seven V. officinalis terpene synthase genes (VoTPSs), two that were functionally characterized as monoterpene synthases and three that preferred farnesyl diphosphate, the substrate for sesquiterpene synthases. The reaction products for two of the sesquiterpene synthases exhibiting root-specific expression were characterized by a combination of GC-MS and NMR in comparison to the terpenes accumulating in planta. VoTPS7 encodes for a synthase that biosynthesizes predominately germacrene C, whereas VoTPS1 catalyzes the conversion of farnesyl diphosphate to valerena-1,10-diene. Using a yeast expression system, specific labeled [13C]acetate, and NMR, we investigated the catalytic mechanism for VoTPS1 and provide evidence for the involvement of a caryophyllenyl carbocation, a cyclobutyl intermediate, in the biosynthesis of valerena-1,10-diene. We suggest a similar mechanism for the biosynthesis of several other biologically related isobutenyl-containing sesquiterpenes.

A Root-Expressed l-Phenylalanine:4-Hydroxyphenylpyruvate Aminotransferase Is Required for Tropane Alkaloid Biosynthesis in Atropa belladonna

This work helps uncover the biosynthetic pathway of the medicinally important tropane alkaloids in Deadly Nightshade. Transcriptomics, silencing, and biochemical analysis reveal the phylogenetically distinct, root-specific aromatic amino acid aminotransferase Ab-ArAT4, which preferentially catalyzes the first step in the formation of littorine, a key intermediate in tropane alkaloid biosynthesis. The tropane alkaloids, hyoscyamine and scopolamine, are medicinal compounds that are the active components of several therapeutics. Hyoscyamine and scopolamine are synthesized in the roots of specific genera of the Solanaceae in a multistep pathway that is only partially elucidated. To facilitate greater understanding of tropane alkaloid biosynthesis, a de novo transcriptome assembly was developed for Deadly Nightshade (Atropa belladonna). Littorine is a key intermediate in hyoscyamine and scopolamine biosynthesis that is produced by the condensation of tropine and phenyllactic acid. Phenyllactic acid is derived from phenylalanine via its transamination to phenylpyruvate, and mining of the transcriptome identified a phylogenetically distinct aromatic amino acid aminotransferase (ArAT), designated Ab-ArAT4, that is coexpressed with known tropane alkaloid biosynthesis genes in the roots of A. belladonna. Silencing of Ab-ArAT4 disrupted synthesis of hyoscyamine and scopolamine through reduction of phenyllactic acid levels. Recombinant Ab-ArAT4 preferentially catalyzes the first step in phenyllactic acid synthesis, the transamination of phenylalanine to phenylpyruvate. However, rather than utilizing the typical keto-acid cosubstrates, 2-oxoglutarate, pyruvate, and oxaloacetate, Ab-ArAT4 possesses strong substrate preference and highest activity with the aromatic keto-acid, 4-hydroxyphenylpyruvate. Thus, Ab-ArAT4 operates at the interface between primary and specialized metabolism, contributing to both tropane alkaloid biosynthesis and the direct conversion of phenylalanine to tyrosine.

The capsicum transcriptome DB: a "hot" tool for genomic research

Chili pepper (Capsicum annuum) is an economically important crop with no available public genome sequence. We describe a genomic resource to facilitate Capsicum annuum research. A collection of Expressed Sequence Tags (ESTs) derived from five C. annuum organs (root, stem, leaf, flower and fruit) were sequenced using the Sanger method and multiple leaf transcriptomes were deeply sampled using with GS-pyrosequencing. A hybrid assembly of 1,324,516 raw reads yielded 32,314 high quality contigs as validated by coverage and identity analysis with existing pepper sequences. Overall, 75.5% of the contigs had significant sequence similarity to entries in nucleic acid and protein databases; 23% of the sequences have not been previously reported for C. annuum and expand sequence resources for this species. A MySQL database and a user-friendly Web interface were constructed with search-tools that permit queries of the ESTs including sequence, functional annotation, Gene Ontology classification, metabolic pathways, and assembly information. The Capsicum Transcriptome DB is free available from http://www.bioingenios.ira.cinvestav.mx:81/Joomla/