Abiotic Stresses Modulate Landscape of Poplar Transcriptome via Alternative Splicing, Differential Intron Retention, and Isoform Ratio Switching.

Abiotic stresses affect plant physiology, development, growth, and alter pre-mRNA splicing. Western poplar is a model woody tree and a potential bioenergy feedstock. To investigate the extent of stress-regulated alternative splicing (AS), we conducted an in-depth survey of leaf, root, and stem xylem transcriptomes under drought, salt, or temperature stress. Analysis of approximately one billion of genome-aligned RNA-Seq reads from tissue- or stress-specific libraries revealed over fifteen millions of novel splice junctions. Transcript models supported by both RNA-Seq and single molecule isoform sequencing (Iso-Seq) data revealed a broad array of novel stress- and/or tissue-specific isoforms. Analysis of Iso-Seq data also resulted in the discovery of 15,087 novel transcribed regions of which 164 show AS. Our findings demonstrate that abiotic stresses profoundly perturb transcript isoform profiles and trigger widespread intron retention (IR) events. Stress treatments often increased or decreased retention of specific introns - a phenomenon described here as differential intron retention (DIR). Many differentially retained introns were regulated in a stress- and/or tissue-specific manner. A subset of transcripts harboring super stress-responsive DIR events showed persisting fluctuations in the degree of IR across all treatments and tissue types. To investigate coordinated dynamics of intron-containing transcripts in the study we quantified absolute copy number of isoforms of two conserved transcription factors (TFs) using Droplet Digital PCR. This case study suggests that stress treatments can be associated with coordinated switches in relative ratios between fully spliced and intron-retaining isoforms and may play a role in adjusting transcriptome to abiotic stresses.

Plant Reactome: a resource for plant pathways and comparative analysis.

Plant Reactome (http://plantreactome.gramene.org/) is a free, open-source, curated plant pathway database portal, provided as part of the Gramene project. The database provides intuitive bioinformatics tools for the visualization, analysis and interpretation of pathway knowledge to support genome annotation, genome analysis, modeling, systems biology, basic research and education. Plant Reactome employs the structural framework of a plant cell to show metabolic, transport, genetic, developmental and signaling pathways. We manually curate molecular details of pathways in these domains for reference species Oryza sativa (rice) supported by published literature and annotation of well-characterized genes. Two hundred twenty-two rice pathways, 1025 reactions associated with 1173 proteins, 907 small molecules and 256 literature references have been curated to date. These reference annotations were used to project pathways for 62 model, crop and evolutionarily significant plant species based on gene homology. Database users can search and browse various components of the database, visualize curated baseline expression of pathway-associated genes provided by the Expression Atlas and upload and analyze their Omics datasets. The database also offers data access via Application Programming Interfaces (APIs) and in various standardized pathway formats, such as SBML and BioPAX.

Gramene 2016: comparative plant genomics and pathway resources.

Gramene (http://www.gramene.org) is an online resource for comparative functional genomics in crops and model plant species. Its two main frameworks are genomes (collaboration with Ensembl Plants) and pathways (The Plant Reactome and archival BioCyc databases). Since our last NAR update, the database website adopted a new Drupal management platform. The genomes section features 39 fully assembled reference genomes that are integrated using ontology-based annotation and comparative analyses, and accessed through both visual and programmatic interfaces. Additional community data, such as genetic variation, expression and methylation, are also mapped for a subset of genomes. The Plant Reactome pathway portal (http://plantreactome.gramene.org) provides a reference resource for analyzing plant metabolic and regulatory pathways. In addition to ∼ 200 curated rice reference pathways, the portal hosts gene homology-based pathway projections for 33 plant species. Both the genome and pathway browsers interface with the EMBL-EBI's Expression Atlas to enable the projection of baseline and differential expression data from curated expression studies in plants. Gramene's archive website (http://archive.gramene.org) continues to provide previously reported resources on comparative maps, markers and QTL. To further aid our users, we have also introduced a live monthly educational webinar series and a Gramene YouTube channel carrying video tutorials.

Environmental stresses modulate abundance and timing of alternatively spliced circadian transcripts in Arabidopsis.

Environmental stresses profoundly altered accumulation of nonsense mRNAs including intron-retaining (IR) transcripts in Arabidopsis. Temporal patterns of stress-induced IR mRNAs were dissected using both oscillating and non-oscillating transcripts. Broad-range thermal cycles triggered a sharp increase in the long IR CCA1 isoforms and altered their phasing to different times of day. Both abiotic and biotic stresses such as drought or Pseudomonas syringae infection induced a similar increase. Thermal stress induced a time delay in accumulation of CCA1 I4Rb transcripts, whereas functional mRNA showed steady oscillations. Our data favor a hypothesis that stress-induced instabilities of the central oscillator can be in part compensated through fluctuations in abundance and out-of-phase oscillations of CCA1 IR transcripts. Taken together, our results support a concept that mRNA abundance can be modulated through altering ratios between functional and nonsense/IR transcripts. SR45 protein specifically bound to the retained CCA1 intron in vitro, suggesting that this splicing factor could be involved in regulation of intron retention. Transcriptomes of nonsense-mediated mRNA decay (NMD)-impaired and heat-stressed plants shared a set of retained introns associated with stress- and defense-inducible transcripts. Constitutive activation of certain stress response networks in an NMD mutant could be linked to disequilibrium between functional and nonsense mRNAs.

The floral transcriptome of Eucalyptus grandis.

As a step toward functional annotation of genes required for floral initiation and development within the Eucalyptus genome, we used short read sequencing to analyze transcriptomes of floral buds from early and late developmental stages, and compared these with transcriptomes of diverse vegetative tissues, including leaves, roots, and stems. A subset of 4807 genes (13% of protein-coding genes) were differentially expressed between floral buds of either stage and vegetative tissues. A similar proportion of genes were differentially expressed among all tissues. A total of 479 genes were differentially expressed between early and late stages of floral development. Gene function enrichment identified 158 gene ontology classes that were overrepresented in floral tissues, including 'pollen development' and 'aromatic compound biosynthetic process'. At least 40 floral-dominant genes lacked functional annotations and thus may be novel floral transcripts. We analyzed several genes and gene families in depth, including 49 putative biomarkers of floral development, the MADS-box transcription factors, 'S-domain'-receptor-like kinases, and selected gene family members with phosphatidylethanolamine-binding protein domains. Expanded MADS-box gene subfamilies in Eucalyptus grandis included SUPPRESSOR OF OVEREXPRESSION OF CO 1 (SOC1), SEPALLATA (SEP) and SHORT VEGETATIVE PHASE (SVP) Arabidopsis thaliana homologs. These data provide a rich resource for functional and evolutionary analysis of genes controlling eucalypt floral development, and new tools for breeding and biotechnology.

Extensive Transcriptome Changes During Natural Onset and Release of Vegetative Bud Dormancy in Populus.

To survive winter, many perennial plants become endodormant, a state of suspended growth maintained even in favorable growing environments. To understand vegetative bud endodormancy, we collected paradormant, endodormant, and ecodormant axillary buds from Populus trees growing under natural conditions. Of 44,441 Populus gene models analyzed using NimbleGen microarrays, we found that 1,362 (3.1%) were differentially expressed among the three dormancy states, and 429 (1.0%) were differentially expressed during only one of the two dormancy transitions (FDR p-value < 0.05). Of all differentially expressed genes, 69% were down-regulated from paradormancy to endodormancy, which was expected given the lower metabolic activity associated with endodormancy. Dormancy transitions were accompanied by changes in genes associated with DNA methylation (via RNA-directed DNA methylation) and histone modifications (via Polycomb Repressive Complex 2), confirming and extending knowledge of chromatin modifications as major features of dormancy transitions. Among the chromatin-associated genes, two genes similar to SPT (SUPPRESSOR OF TY) were strongly up-regulated during endodormancy. Transcription factor genes and gene sets that were atypically up-regulated during endodormancy include a gene that seems to encode a trihelix transcription factor and genes associated with proteins involved in responses to ethylene, cold, and other abiotic stresses. These latter transcription factors include ETHYLENE INSENSITIVE 3 (EIN3), ETHYLENE-RESPONSIVE ELEMENT BINDING PROTEIN (EBP), ETHYLENE RESPONSE FACTOR (ERF), ZINC FINGER PROTEIN 10 (ZAT10), ZAT12, and WRKY DNA-binding domain proteins. Analyses of phytohormone-associated genes suggest important changes in responses to ethylene, auxin, and brassinosteroids occur during endodormancy. We found weaker evidence for changes in genes associated with salicylic acid and jasmonic acid, and little evidence for important changes in genes associated with gibberellins, abscisic acid, and cytokinin. We identified 315 upstream sequence motifs associated with eight patterns of gene expression, including novel motifs and motifs associated with the circadian clock and responses to photoperiod, cold, dehydration, and ABA. Analogies between flowering and endodormancy suggest important roles for genes similar to SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL), DORMANCY ASSOCIATED MADS-BOX (DAM), and SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1).

Gramene 2013: comparative plant genomics resources.

Gramene (http://www.gramene.org) is a curated online resource for comparative functional genomics in crops and model plant species, currently hosting 27 fully and 10 partially sequenced reference genomes in its build number 38. Its strength derives from the application of a phylogenetic framework for genome comparison and the use of ontologies to integrate structural and functional annotation data. Whole-genome alignments complemented by phylogenetic gene family trees help infer syntenic and orthologous relationships. Genetic variation data, sequences and genome mappings available for 10 species, including Arabidopsis, rice and maize, help infer putative variant effects on genes and transcripts. The pathways section also hosts 10 species-specific metabolic pathways databases developed in-house or by our collaborators using Pathway Tools software, which facilitates searches for pathway, reaction and metabolite annotations, and allows analyses of user-defined expression datasets. Recently, we released a Plant Reactome portal featuring 133 curated rice pathways. This portal will be expanded for Arabidopsis, maize and other plant species. We continue to provide genetic and QTL maps and marker datasets developed by crop researchers. The project provides a unique community platform to support scientific research in plant genomics including studies in evolution, genetics, plant breeding, molecular biology, biochemistry and systems biology.

A genome scale metabolic network for rice and accompanying analysis of tryptophan, auxin and serotonin biosynthesis regulation under biotic stress.

BACKGROUND: Functional annotations of large plant genome projects mostly provide information on gene function and gene families based on the presence of protein domains and gene homology, but not necessarily in association with gene expression or metabolic and regulatory networks. These additional annotations are necessary to understand the physiology, development and adaptation of a plant and its interaction with the environment. RESULTS: RiceCyc is a metabolic pathway networks database for rice. It is a snapshot of the substrates, metabolites, enzymes, reactions and pathways of primary and intermediary metabolism in rice. RiceCyc version 3.3 features 316 pathways and 6,643 peptide-coding genes mapped to 2,103 enzyme-catalyzed and 87 protein-mediated transport reactions. The initial functional annotations of rice genes with InterPro, Gene Ontology, MetaCyc, and Enzyme Commission (EC) numbers were enriched with annotations provided by KEGG and Gramene databases. The pathway inferences and the network diagrams were first predicted based on MetaCyc reference networks and plant pathways from the Plant Metabolic Network, using the Pathologic module of Pathway Tools. This was enriched by manually adding metabolic pathways and gene functions specifically reported for rice. The RiceCyc database is hierarchically browsable from pathway diagrams to the associated genes, metabolites and chemical structures. Through the integrated tool OMICs Viewer, users can upload transcriptomic, proteomic and metabolomic data to visualize expression patterns in a virtual cell. RiceCyc, along with additional species-specific pathway databases hosted in the Gramene project, facilitates comparative pathway analysis. CONCLUSIONS: Here we describe the RiceCyc network development and discuss its contribution to rice genome annotations. As a case study to demonstrate the use of RiceCyc network as a discovery environment we carried out an integrated bioinformatic analysis of rice metabolic genes that are differentially regulated under diurnal photoperiod and biotic stress treatments. The analysis of publicly available rice transcriptome datasets led to the hypothesis that the complete tryptophan biosynthesis and its dependent metabolic pathways including serotonin biosynthesis are induced by taxonomically diverse pathogens while also being under diurnal regulation. The RiceCyc database is available online for free access at http://www.gramene.org/pathway/.

Contrasting patterns of evolution following whole genome versus tandem duplication events in Populus.

Comparative analysis of multiple angiosperm genomes has implicated gene duplication in the expansion and diversification of many gene families. However, empirical data and theory suggest that whole-genome and small-scale duplication events differ with respect to the types of genes preserved as duplicate pairs. We compared gene duplicates resulting from a recent whole genome duplication to a set of tandemly duplicated genes in the model forest tree Populus trichocarpa. We used a combination of microarray expression analyses of a diverse set of tissues and functional annotation to assess factors related to the preservation of duplicate genes of both types. Whole genome duplicates are 700 bp longer and are expressed in 20% more tissues than tandem duplicates. Furthermore, certain functional categories are over-represented in each class of duplicates. In particular, disease resistance genes and receptor-like kinases commonly occur in tandem but are significantly under-retained following whole genome duplication, while whole genome duplicate pairs are enriched for members of signal transduction cascades and transcription factors. The shape of the distribution of expression divergence for duplicated pairs suggests that nearly half of the whole genome duplicates have diverged in expression by a random degeneration process. The remaining pairs have more conserved gene expression than expected by chance, consistent with a role for selection under the constraints of gene balance. We hypothesize that duplicate gene preservation in Populus is driven by a combination of subfunctionalization of duplicate pairs and purifying selection favoring retention of genes encoding proteins with large numbers of interactions.

Global profiling of rice and poplar transcriptomes highlights key conserved circadian-controlled pathways and cis-regulatory modules.

BACKGROUND: Circadian clocks provide an adaptive advantage through anticipation of daily and seasonal environmental changes. In plants, the central clock oscillator is regulated by several interlocking feedback loops. It was shown that a substantial proportion of the Arabidopsis genome cycles with phases of peak expression covering the entire day. Synchronized transcriptome cycling is driven through an extensive network of diurnal and clock-regulated transcription factors and their target cis-regulatory elements. Study of the cycling transcriptome in other plant species could thus help elucidate the similarities and differences and identify hubs of regulation common to monocot and dicot plants. METHODOLOGY/PRINCIPAL FINDINGS: Using a combination of oligonucleotide microarrays and data mining pipelines, we examined daily rhythms in gene expression in one monocotyledonous and one dicotyledonous plant, rice and poplar, respectively. Cycling transcriptomes were interrogated under different diurnal (driven) and circadian (free running) light and temperature conditions. Collectively, photocycles and thermocycles regulated about 60% of the expressed nuclear genes in rice and poplar. Depending on the condition tested, up to one third of oscillating Arabidopsis-poplar-rice orthologs were phased within three hours of each other suggesting a high degree of conservation in terms of rhythmic gene expression. We identified clusters of rhythmically co-expressed genes and searched their promoter sequences to identify phase-specific cis-elements, including elements that were conserved in the promoters of Arabidopsis, poplar, and rice. CONCLUSIONS/SIGNIFICANCE: Our results show that the cycling patterns of many circadian clock genes are highly conserved across poplar, rice, and Arabidopsis. The expression of many orthologous genes in key metabolic and regulatory pathways is diurnal and/or circadian regulated and phased to similar times of day. Our results confirm previous findings in Arabidopsis of three major classes of cis-regulatory modules within the plant circadian network: the morning (ME, GBOX), evening (EE, GATA), and midnight (PBX/TBX/SBX) modules. Identification of identical overrepresented motifs in the promoters of cycling genes from different species suggests that the core diurnal/circadian cis-regulatory network is deeply conserved between mono- and dicotyledonous species.

Gibberellin-associated cisgenes modify growth, stature and wood properties in Populus.

We studied the effects on plant growth from insertion of five cisgenes that encode proteins involved in gibberellin metabolism or signalling. Intact genomic copies of PtGA20ox7, PtGA2ox2,Pt RGL1_1, PtRGL1_2 and PtGAI1 genes from the genome-sequenced Populus trichocarpa clone Nisqually-1 were transformed into Populus tremula × alba (clone INRA 717-1B4), and growth, morphology and xylem cell size characterized in the greenhouse. Each cisgene encompassed 1-2 kb of 5' and 1 kb of 3' flanking DNA, as well as all native exons and introns. Large numbers of independent insertion events per cisgene (19-38), including empty vector controls, were studied. Three of the cisgenic modifications had significant effects on plant growth rate, morphology or wood properties. The PtGA20ox7 cisgene increased rate of shoot regeneration in vitro, accelerated early growth, and variation in growth rate was correlated with PtGA20ox7 gene expression. PtRGL1_1 and PtGA2ox2 caused reduced growth, while PtRGL1_2 gave rise to plants that grew normally but had significantly longer xylem fibres. RT-PCR studies suggested that the lack of growth inhibition observed in PtRGL1_2 cisgenic plants was a result of co-suppression. PtGAI1 slowed regeneration rate and both PtGAI1 and PtGA20ox7 gave rise to increased variance among events for early diameter and volume index, respectively. Our work suggests that cisgenic insertion of additional copies of native genes involved in growth regulation may provide tools to help modify plant architecture, expand the genetic variance in plant architecture available to breeders and accelerate transfer of alleles between difficult-to-cross species.

The genome of woodland strawberry (Fragaria vesca).

The woodland strawberry, Fragaria vesca (2n = 2x = 14), is a versatile experimental plant system. This diminutive herbaceous perennial has a small genome (240 Mb), is amenable to genetic transformation and shares substantial sequence identity with the cultivated strawberry (Fragaria × ananassa) and other economically important rosaceous plants. Here we report the draft F. vesca genome, which was sequenced to ×39 coverage using second-generation technology, assembled de novo and then anchored to the genetic linkage map into seven pseudochromosomes. This diploid strawberry sequence lacks the large genome duplications seen in other rosids. Gene prediction modeling identified 34,809 genes, with most being supported by transcriptome mapping. Genes critical to valuable horticultural traits including flavor, nutritional value and flowering time were identified. Macrosyntenic relationships between Fragaria and Prunus predict a hypothetical ancestral Rosaceae genome that had nine chromosomes. New phylogenetic analysis of 154 protein-coding genes suggests that assignment of Populus to Malvidae, rather than Fabidae, is warranted.

Gramene database in 2010: updates and extensions.

Now in its 10th year, the Gramene database (http://www.gramene.org) has grown from its primary focus on rice, the first fully-sequenced grass genome, to become a resource for major model and crop plants including Arabidopsis, Brachypodium, maize, sorghum, poplar and grape in addition to several species of rice. Gramene began with the addition of an Ensembl genome browser and has expanded in the last decade to become a robust resource for plant genomics hosting a wide array of data sets including quantitative trait loci (QTL), metabolic pathways, genetic diversity, genes, proteins, germplasm, literature, ontologies and a fully-structured markers and sequences database integrated with genome browsers and maps from various published studies (genetic, physical, bin, etc.). In addition, Gramene now hosts a variety of web services including a Distributed Annotation Server (DAS), BLAST and a public MySQL database. Twice a year, Gramene releases a major build of the database and makes interim releases to correct errors or to make important updates to software and/or data.

Genome-wide transcriptome analysis of the transition from primary to secondary stem development in Populus trichocarpa.

BACKGROUND: With its genome sequence and other experimental attributes, Populus trichocarpa has become the model species for genomic studies of wood development. Wood is derived from secondary growth of tree stems, and begins with the development of a ring of vascular cambium in the young developing stem. The terminal region of the developing shoot provides a steep developmental gradient from primary to secondary growth that facilitates identification of genes that play specialized functions during each of these phases of growth. RESULTS: Using a genomic microarray representing the majority of the transcriptome, we profiled gene expression in stem segments that spanned primary to secondary growth. We found 3,016 genes that were differentially expressed during stem development (Q-value 2-fold expression variation), and 15% of these genes encode proteins with no significant identities to known genes. We identified all gene family members putatively involved in secondary growth for carbohydrate active enzymes, tubulins, actins, actin depolymerizing factors, fasciclin-like AGPs, and vascular development-associated transcription factors. Almost 70% of expressed transcription factors were upregulated during the transition to secondary growth. The primary shoot elongation region of the stem contained specific carbohydrate active enzyme and expansin family members that are likely to function in primary cell wall synthesis and modification. Genes involved in plant defense and protective functions were also dominant in the primary growth region. CONCLUSION: Our results describe the global patterns of gene expression that occur during the transition from primary to secondary stem growth. We were able to identify three major patterns of gene expression and over-represented gene ontology categories during stem development. The new regulatory factors and cell wall biogenesis genes that we identified provide candidate genes for further functional characterization, as well as new tools for molecular breeding and biotechnology aimed at improvement of tree growth rate, crown form, and wood quality.

Genome scale transcriptome analysis of shoot organogenesis in Populus.

BACKGROUND: Our aim is to improve knowledge of gene regulatory circuits important to dedifferentiation, redifferentiation, and adventitious meristem organization during in vitro regeneration of plants. Regeneration of transgenic cells remains a major obstacle to research and commercial deployment of most taxa of transgenic plants, and woody species are particularly recalcitrant. The model woody species Populus, due to its genome sequence and amenability to in vitro manipulation, is an excellent species for study in this area. The genes recognized may help to guide the development of new tools for improving the efficiency of plant regeneration and transformation. RESULTS: We analyzed gene expression during poplar in vitro dedifferentiation and shoot regeneration using an Affymetrix array representing over 56,000 poplar transcripts. We focused on callus induction and shoot formation, thus we sampled RNAs from tissues: prior to callus induction, 3 days and 15 days after callus induction, and 3 days and 8 days after the start of shoot induction. We used a female hybrid white poplar clone (INRA 717-1 B4, Populus tremula x P. alba) that is used widely as a model transgenic genotype. Approximately 15% of the monitored genes were significantly up-or down-regulated when controlling the false discovery rate (FDR) at 0.01; over 3,000 genes had a 5-fold or greater change in expression. We found a large initial change in expression after the beginning of hormone treatment (at the earliest stage of callus induction), and then a much smaller number of additional differentially expressed genes at subsequent regeneration stages. A total of 588 transcription factors that were distributed in 45 gene families were differentially regulated. Genes that showed strong differential expression included components of auxin and cytokinin signaling, selected cell division genes, and genes related to plastid development and photosynthesis. When compared with data on in vitro callogenesis in Arabidopsis, 25% (1,260) of up-regulated and 22% (748) of down-regulated genes were in common with the genes regulated in poplar during callus induction. CONCLUSION: The major regulatory events during plant cell organogenesis occur at early stages of dedifferentiation. The regulatory circuits reflect the combinational effects of transcriptional control and hormone signaling, and associated changes in light environment imposed during dedifferentiation.