Deep Transcriptome Analysis Research Articles

Rapidly Evolving Genes and Stress Adaptation of Two Desert Poplars, Populus euphratica and P. pruinosa

Understanding which genes have evolved rapidly with the recent tree speciation in arid habitats can provide valuable insights into different adaptation mechanisms. We employed a comparative evolutionary analysis of expressed sequence tags (ESTs) from two desert poplars, Populus pruinosa and P. euphratica, which diverged in the recent past. Following an approach taken previously with P. euphratica, we conducted a deep transcriptomic analysis of P. pruinosa. To maximize representation of conditional transcripts, mRNA was obtained from living tissues of two types of callus and desert-grown trees. De novo assembly generated 114,866 high-quality unique sequences using Solexa sequence data. Following assembly we were able to identify, with high confidence, 2859 orthologous sequence pairs between the two species. Based on the ratio of nonsynonymous (Ka) to synonymous (Ks) substitutions, we identified a total of 84 (2.9%) ortholog pairs exhibiting rapid evolution with signs of strong selection (Ka/Ks>1). Genes homologous to these ortholog pairs in model species are mainly involved in ‘responses to stress’, ‘ubiquitin-dependent protein catabolic processes’, and ‘biological regulation’. Finally, we examined the expression patterns of candidate genes with rapid evolution in response to salt stress. Only one pair of orthologs up-regulated their expression in both species while three and four genes were found to up-regulated in P. pruinosa and in P. euphratica respectively. Our findings together suggest that the genes at the same category or network but with differentiated expressions or functions may have evolved rapidly during adaptive divergence of the two species to differentiated salty desert habitats.

The duck genome and transcriptome provide insight into an avian influenza virus reservoir species

The duck (Anas platyrhynchos) is one of the principal natural hosts of influenza A viruses. We present the duck genome sequence and perform deep transcriptome analyses to investigate immune-related genes. Our data indicate that the duck possesses a contractive immune gene repertoire, as in chicken and zebra finch, and this repertoire has been shaped through lineage-specific duplications. We identify genes that are responsive to influenza A viruses using the lung transcriptomes of control ducks and ones that were infected with either a highly pathogenic (A/duck/Hubei/49/05) or a weakly pathogenic (A/goose/Hubei/65/05) H5N1 virus. Further, we show how the duck’s defense mechanisms against influenza infection have been optimized through the diversification of its β-defensin and butyrophilin-like repertoires. These analyses, in combination with the genomic and transcriptomic data, provide a resource for characterizing the interaction between host and influenza viruses.

Phytochemical Genomics on the Way

This is a really exciting time for phytochemistry research. Plant scientists studying specialized plant metabolism (previously often known as ‘secondary metabolism’) can foresee tremendous advances in their understanding of the mechanisms underlying the production and function of an enormous variety of plant metabolites, and of their regulation and evolution. Such remarkable progress would never have been expected 10 years ago and has taken place in just the last few years (Yonekura-Sakakibara and Saito 2009). Two technological breakthroughs seem to have been the key triggers for this unforeseen progress: metabolomics and mass DNA sequencing. Metabolomics is <15 years old, yet the concept was immediately applied to the plant field as scientists realized its tremendous potential for understanding the capacity of plants for making a huge array of compounds (Bino et al. 2004). Metabolomics now makes possible comprehensive profiling of nearly all the metabolites that accumulate in plant cells (Saito and Matsuda 2010). The second technological breakthrough was ultra-high-throughput DNA sequencing, with so-called ‘next-generation DNA sequencers’ (Wang et al. 2009, Ozsolak and Milos 2011). Until the realization of this technology, no one would have anticipated genomics-based studies of hundreds of medicinal or exotic plant species— genomics was only for model plants or major crops, whose genome sequences could only be revealed by large international research consortia (or by huge well-funded undertakings). However, newly developed mass DNA sequencers have made comprehensive genome and transcriptome sequencing possible for ordinary phytochemical researchers who are investigating the synthesis of particular plant metabolites out of their own curiosity. By combining data sets from metabolomics and genomics/transcriptomics, one can understand relationships between metabolites and genes that allow the generation of testable hypotheses about the functions of genes and metabolites (Higashi and Saito 2013, Hur et al. 2013, Yonekura-Sakakibara et al. 2013). During this integration and hypothesis-generation phase, databases and bioinformatics play indispensable roles in reducing this ‘systems biology’ workload. Such studies will provide us with an understanding of the evolution of genes and pathways and hence in-depth fundamental insights into plant life. This special focus issue (SFI) on ‘Phytochemical Genomics’ is intended to serve as a synthesis of the most up to date information in this field, indicating not only current trends but also its future prospects. For this SFI of Plant and Cell Physiology, we invited one review article and six original articles reporting on the exciting and growing field of phytochemical genomics. The huge numbers of specialized metabolites produced in plants are mainly classified into three basic groups: alkaloids, terpenoids and phenolic compounds. Among alkaloids, benzylisoquinoline alkaloids (BIAs), including the narcotic analgesics codeine and morphine, are one of the groups of compounds most exploited for their medicinal properties. In this SFI, Hagel and Facchini (2013; see pp. 647–672) review recent advances in biosynthetic studies of BIAs including transcriptomics, proteomics and metabolomics; they also discuss the application of synthetic biology to the development of production by microbes as an alternative to plants as a potential commercial source of valuable BIAs. Another three deep transcritome analyses are reported in this SFI. Catharanthus roseus synthesizes numerous terpenoid indole alkaloids, such as the anticancer drugs vinblastine and vincristine. So far, pathway databases and metabolic networks reconstructed from omics sequence data have been mainly built for model plant species such as Arabidopis and rice. Van Moerkercke et al. (2013; see pp. 673–685) conducted a deep transcriptome analysis in the medicinal plant, C. roseus, and constructed a detailed metabolic pathway database ‘CathaCyc’ that contains 390 pathways with more than a thousand assigned enzymes. Yamazaki et al. (2013; see pp. 686–696) coupled the deep transcriptome analysis with untargeted metabolic profiling in Ophiorrhiza pumila, which produces the anti-cancer alkaloids camptothecin and anthraquinones. Yamazaki and colleagues compared both transcriptome and metabolome data sets in the alkaloid-producing hairy root with a cell suspension that does not produce alkaloids. Glycyrrhiza plants produce various phytochemicals, among them many different terpenoids and flavonoids. Among the phytochemicals, glycyrrhizin exhibits several sorts of pharmacological activity and is also used as a natural sweetener. Two cytochrome P450 monooxygenase (P450s) genes for glycyrrhizin biosynthesis have been recently identified (Seki et al. 2008, Seki et al. 2011). Ramilowski et al. (2013; see pp. 697–710) conducted deep transcriptome analysis in G. uralensis, and

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.

Use of Mutagenesis, Genetic Mapping and Next Generation Transcriptomics to Investigate Insecticide Resistance Mechanisms

Insecticide resistance is a worldwide problem with major impact on agriculture and human health. Understanding the underlying molecular mechanisms is crucial for the management of the phenomenon; however, this information often comes late with respect to the implementation of efficient counter-measures, particularly in the case of metabolism-based resistance mechanisms. We employed a genome-wide insertional mutagenesis screen to Drosophila melanogaster, using a Minos-based construct, and retrieved a line (MiT[w−]3R2) resistant to the neonicotinoid insecticide Imidacloprid. Biochemical and bioassay data indicated that resistance was due to increased P450 detoxification. Deep sequencing transcriptomic analysis revealed substantial over- and under-representation of 357 transcripts in the resistant line, including statistically significant changes in mixed function oxidases, peptidases and cuticular proteins. Three P450 genes (Cyp4p2, Cyp6a2 and Cyp6g1) located on the 2R chromosome, are highly up-regulated in mutant flies compared to susceptible Drosophila. One of them (Cyp6g1) has been already described as a major factor for Imidacloprid resistance, which validated the approach. Elevated expression of the Cyp4p2 was not previously documented in Drosophila lines resistant to neonicotinoids. In silico analysis using the Drosophila reference genome failed to detect transcription binding factors or microRNAs associated with the over-expressed Cyp genes. The resistant line did not contain a Minos insertion in its chromosomes, suggesting a hit-and-run event, i.e. an insertion of the transposable element, followed by an excision which caused the mutation. Genetic mapping placed the resistance locus to the right arm of the second chromosome, within a ∼1 Mb region, where the highly up-regulated Cyp6g1 gene is located. The nature of the unknown mutation that causes resistance is discussed on the basis of these results.

The genome of Tetranychus urticae reveals herbivorous pest adaptations

The spider mite Tetranychus urticae is a cosmopolitan agricultural pest with an extensive host plant range and an extreme record of pesticide resistance. Here we present the completely sequenced and annotated spider mite genome, representing the first complete chelicerate genome. At 90 megabases T. urticae has the smallest sequenced arthropod genome. Compared with other arthropods, the spider mite genome shows unique changes in the hormonal environment and organization of the Hox complex, and also reveals evolutionary innovation of silk production. We find strong signatures of polyphagy and detoxification in gene families associated with feeding on different hosts and in new gene families acquired by lateral gene transfer. Deep transcriptome analysis of mites feeding on different plants shows how this pest responds to a changing host environment. The T. urticae genome thus offers new insights into arthropod evolution and plant–herbivore interactions, and provides unique opportunities for developing novel plant protection strategies.

Deep sequencing-based transcriptome analysis of Plutella xylostella larvae parasitized by Diadegma semiclausum

BackgroundParasitoid insects manipulate their hosts' physiology by injecting various factors into their host upon parasitization. Transcriptomic approaches provide a powerful approach to study insect host-parasitoid interactions at the molecular level. In order to investigate the effects of parasitization by an ichneumonid wasp (Diadegma semiclausum) on the host (Plutella xylostella), the larval transcriptome profile was analyzed using a short-read deep sequencing method (Illumina). Symbiotic polydnaviruses (PDVs) associated with ichneumonid parasitoids, known as ichnoviruses, play significant roles in host immune suppression and developmental regulation. In the current study, D. semiclausum ichnovirus (DsIV) genes expressed in P. xylostella were identified and their sequences compared with other reported PDVs. Five of these genes encode proteins of unknown identity, that have not previously been reported.ResultsDe novo assembly of cDNA sequence data generated 172,660 contigs between 100 and 10000 bp in length; with 35% of > 200 bp in length. Parasitization had significant impacts on expression levels of 928 identified insect host transcripts. Gene ontology data illustrated that the majority of the differentially expressed genes are involved in binding, catalytic activity, and metabolic and cellular processes. In addition, the results show that transcription levels of antimicrobial peptides, such as gloverin, cecropin E and lysozyme, were up-regulated after parasitism. Expression of ichnovirus genes were detected in parasitized larvae with 19 unique sequences identified from five PDV gene families including vankyrin, viral innexin, repeat elements, a cysteine-rich motif, and polar residue rich protein. Vankyrin 1 and repeat element 1 genes showed the highest transcription levels among the DsIV genes.ConclusionThis study provides detailed information on differential expression of P. xylostella larval genes following parasitization, DsIV genes expressed in the host and also improves our current understanding of this host-parasitoid interaction.

Genome-scale transcriptome analysis of the desert poplar, Populus euphratica

Populus euphratica is well-adapted to extreme desert environments and is an important model species for studying the effects of abiotic stresses on trees. Here we present the first deep transcriptomic analysis of this species. To maximize representation of conditional transcripts, mRNA was obtained from living tissues of desert-grown trees and two types of callus (salt-stressed and unstressed). De novo assembly generated 86,777 Unigenes using Solexa sequence data. These sequences covered 92% of previously reported P. euphratica expressed sequence tags (ESTs) and 90% of the TIGR poplar ESTs, and a total of 58,499 high-quality unique sequences were annotated by BLAST similarity searches against public databases. We found that 27% of the total Unigenes were differentially expressed (up- or down-regulated) in response to salt stress in P. euphratica callus. These differentially expressed genes are mainly involved in transport, transcription, cellular communication and metabolism. In addition, we found that numerous putative genes involved in ABA regulation and biosynthesis were also differentially regulated. This study represents the deepest transcriptomic and gene-annotation analysis of P. euphratica to date. The genetic knowledge acquired should be very useful for future studies of the molecular adaptation of this tree species to abiotic stress and facilitate genetic manipulation of other poplar species.

Integration of deep transcriptome and proteome analyses reveals the components of alkaloid metabolism in opium poppy cell cultures

BackgroundPapaver somniferum (opium poppy) is the source for several pharmaceutical benzylisoquinoline alkaloids including morphine, the codeine and sanguinarine. In response to treatment with a fungal elicitor, the biosynthesis and accumulation of sanguinarine is induced along with other plant defense responses in opium poppy cell cultures. The transcriptional induction of alkaloid metabolism in cultured cells provides an opportunity to identify components of this process via the integration of deep transcriptome and proteome databases generated using next-generation technologies.ResultsA cDNA library was prepared for opium poppy cell cultures treated with a fungal elicitor for 10 h. Using 454 GS-FLX Titanium pyrosequencing, 427,369 expressed sequence tags (ESTs) with an average length of 462 bp were generated. Assembly of these sequences yielded 93,723 unigenes, of which 23,753 were assigned Gene Ontology annotations. Transcripts encoding all known sanguinarine biosynthetic enzymes were identified in the EST database, 5 of which were represented among the 50 most abundant transcripts. Liquid chromatography-tandem mass spectrometry (LC-MS/MS) of total protein extracts from cell cultures treated with a fungal elicitor for 50 h facilitated the identification of 1,004 proteins. Proteins were fractionated by one-dimensional SDS-PAGE and digested with trypsin prior to LC-MS/MS analysis. Query of an opium poppy-specific EST database substantially enhanced peptide identification. Eight out of 10 known sanguinarine biosynthetic enzymes and many relevant primary metabolic enzymes were represented in the peptide database.ConclusionsThe integration of deep transcriptome and proteome analyses provides an effective platform to catalogue the components of secondary metabolism, and to identify genes encoding uncharacterized enzymes. The establishment of corresponding transcript and protein databases generated by next-generation technologies in a system with a well-defined metabolite profile facilitates an improved linkage between genes, enzymes, and pathway components. The proteome database represents the most relevant alkaloid-producing enzymes, compared with the much deeper and more complete transcriptome library. The transcript database contained full-length mRNAs encoding most alkaloid biosynthetic enzymes, which is a key requirement for the functional characterization of novel gene candidates.

Deep and Comparative Transcriptome Analysis of Rice Plants Infested by the Beet Armyworm (Spodoptera exigua) and Water Weevil (Lissorhoptrus oryzophilus)

Abstract The beet armyworm (Spodoptera exigua) and the rice water weevil (Lissorhoptrus oryzophilus) are two important insect pests in rice production. To identify insect-responsive genes in rice, we performed a deep transcriptome analysis of Nipponbare rice leaves infested with both beet armyworm and water weevil using massively parallel signature sequencing (MPSS). Many antisense, alternative, and novel transcripts were commonly and specifically induced and suppressed in the infested tissue. Key genes involved in the defense metabolic pathways such as salicylic acid and jasmonic acid biosynthesis pathways were up-regulated in the infested leaves. To validate the MPSS results, we analyzed the transcriptome of the rice leaves infested with water weevils using Solexa’s sequencing-by-synthesis (SBS) method. The MPSS and SBS data were highly correlated (Pearson’s correlation coefficient = 0.85), and 83% of genes had similar gene expression in both libraries. Our comprehensive and in-depth survey of the insect-infested libraries provides a rich genomic resource for further analyzing the function of key regulatory genes involved in insect resistance in rice.

Transcriptome response to pollutants and insecticides in the dengue vector Aedes aegypti using next-generation sequencing technology

BackgroundThe control of mosquitoes transmitting infectious diseases relies mainly on the use of chemical insecticides. However, mosquito control programs are now threatened by the emergence of insecticide resistance. Hitherto, most research efforts have been focused on elucidating the molecular basis of inherited resistance. Less attention has been paid to the short-term response of mosquitoes to insecticides and pollutants which could have a significant impact on insecticide efficacy. Here, a combination of LongSAGE and Solexa sequencing was used to perform a deep transcriptome analysis of larvae of the dengue vector Aedes aegypti exposed for 48 h to sub-lethal doses of three chemical insecticides and three anthropogenic pollutants.ResultsThirty millions 20 bp cDNA tags were sequenced, mapped to the mosquito genome and clustered, representing 6850 known genes and 4868 additional clusters not located within predicted genes. Mosquitoes exposed to insecticides or anthropogenic pollutants showed considerable modifications of their transcriptome. Genes encoding cuticular proteins, transporters, and enzymes involved in the mitochondrial respiratory chain and detoxification processes were particularly affected. Genes and molecular mechanisms potentially involved in xenobiotic response and insecticide tolerance were identified.ConclusionsThe method used in the present study appears as a powerful approach for investigating fine transcriptome variations in genome-sequenced organisms and can provide useful informations for the detection of novel transcripts. At the biological level, despite low concentrations and no apparent phenotypic effects, the significant impact of these xenobiotics on mosquito transcriptomes raise important questions about the 'hidden impact' of anthropogenic pollutants on ecosystems and consequences on vector control.

A combination of LongSAGE with Solexa sequencing is well suited to explore the depth and the complexity of transcriptome

Background"Open" transcriptome analysis methods allow to study gene expression without a priori knowledge of the transcript sequences. As of now, SAGE (Serial Analysis of Gene Expression), LongSAGE and MPSS (Massively Parallel Signature Sequencing) are the mostly used methods for "open" transcriptome analysis. Both LongSAGE and MPSS rely on the isolation of 21 pb tag sequences from each transcript. In contrast to LongSAGE, the high throughput sequencing method used in MPSS enables the rapid sequencing of very large libraries containing several millions of tags, allowing deep transcriptome analysis. However, a bias in the complexity of the transcriptome representation obtained by MPSS was recently uncovered.ResultsIn order to make a deep analysis of mouse hypothalamus transcriptome avoiding the limitation introduced by MPSS, we combined LongSAGE with the Solexa sequencing technology and obtained a library of more than 11 millions of tags. We then compared it to a LongSAGE library of mouse hypothalamus sequenced with the Sanger method.ConclusionWe found that Solexa sequencing technology combined with LongSAGE is perfectly suited for deep transcriptome analysis. In contrast to MPSS, it gives a complex representation of transcriptome as reliable as a LongSAGE library sequenced by the Sanger method.

Deep and comparative analysis of the mycelium and appressorium transcriptomes of Magnaporthe grisea using MPSS, RL-SAGE, and oligoarray methods

BackgroundRice blast, caused by the fungal pathogen Magnaporthe grisea, is a devastating disease causing tremendous yield loss in rice production. The public availability of the complete genome sequence of M. grisea provides ample opportunities to understand the molecular mechanism of its pathogenesis on rice plants at the transcriptome level. To identify all the expressed genes encoded in the fungal genome, we have analyzed the mycelium and appressorium transcriptomes using massively parallel signature sequencing (MPSS), robust-long serial analysis of gene expression (RL-SAGE) and oligoarray methods.ResultsThe MPSS analyses identified 12,531 and 12,927 distinct significant tags from mycelia and appressoria, respectively, while the RL-SAGE analysis identified 16,580 distinct significant tags from the mycelial library. When matching these 12,531 mycelial and 12,927 appressorial significant tags to the annotated CDS, 500 bp upstream and 500 bp downstream of CDS, 6,735 unique genes in mycelia and 7,686 unique genes in appressoria were identified. A total of 7,135 mycelium-specific and 7,531 appressorium-specific significant MPSS tags were identified, which correspond to 2,088 and 1,784 annotated genes, respectively, when matching to the same set of reference sequences. Nearly 85% of the significant MPSS tags from mycelia and appressoria and 65% of the significant tags from the RL-SAGE mycelium library matched to the M. grisea genome. MPSS and RL-SAGE methods supported the expression of more than 9,000 genes, representing over 80% of the predicted genes in M. grisea. About 40% of the MPSS tags and 55% of the RL-SAGE tags represent novel transcripts since they had no matches in the existing M. grisea EST collections. Over 19% of the annotated genes were found to produce both sense and antisense tags in the protein-coding region. The oligoarray analysis identified the expression of 3,793 mycelium-specific and 4,652 appressorium-specific genes. A total of 2,430 mycelial genes and 1,886 appressorial genes were identified by both MPSS and oligoarray.ConclusionThe comprehensive and deep transcriptome analysis by MPSS and RL-SAGE methods identified many novel sense and antisense transcripts in the M. grisea genome at two important growth stages. The differentially expressed transcripts that were identified, especially those specifically expressed in appressoria, represent a genomic resource useful for gaining a better understanding of the molecular basis of M. grisea pathogenicity. Further analysis of the novel antisense transcripts will provide new insights into the regulation and function of these genes in fungal growth, development and pathogenesis in the host plants.