An EST dataset for Metasequoia glyptostroboides buds: the first EST resource for molecular genomics studies in Metasequoia

De novo assembly of the Platycladus orientalis (L.) Franco transcriptome provides insight into the development and pollination mechanism of female cone based on RNA-Seq data

For seed-bearing plants, the basis of seed and fruit formation is pollination. The normal progression of pollination is through advances in continuous signal exchange and material transfer, which occur mainly in female reproductive organs; thus, the molecular mechanism of development in female reproductive organs is vital for understanding the principle of pollination. However, molecular biology studies on the development of female cones related to pollination are rare and unclear in gymnosperms, especially in Cupressaceae. In this study, Platycladus orientalis, a monotypic genus within Cupressaceae, was chosen to examine female cone transcriptomes at pre-pollination and pollination stages by Illumina paired-end sequencing technology to de novo sequence six libraries with 3 biological replicates. These libraries were used to construct a P. orientalis transcriptome database containing 71,669 unigenes (4,963 upregulated unigenes and 11,747 downregulated unigenes at the pollination stage) for subsequent analysis. Based on the annotations and expression levels, the functions of differentially expressed unigenes and enriched pathways between the developmental processes of female cones were analysed to detail the preliminary development and pollination mechanism of the female cone. Targeted investigations were specifically performed to determine the elementary mechanism of secretion and functioning of the pollination drop, a vital ovule secretion at the pollination stage. Ultimately, the expression of 15 unigenes selected between two stages were further assessed and confirmed using qRT-PCR, which demonstrated reliable data and significant differences in the expression profiles of key genes. As one of the largest available transcriptomic resources of this species, the database is constructed to prospectively adapt to the physiological and genomic data of woody plants. This work provided the first transcriptome profile of P. orientalis female cones at different developmental stages, and will promote the illumination of the pollination mechanism of P. orientalis, and will serve as the basis for in-depth genomic study in the Cupressaceae family. This initiative will arouse the interest and attention of scholars and pave the way for future studies.

Comparison of Reliable Reference Genes Following Different Hormone Treatments by Various Algorithms for qRT-PCR Analysis of Metasequoia

Quantitative reverse transcription polymerase chain reaction (qRT-PCR) is the most sensitive technique for evaluating gene expression levels. Choosing appropriate reference genes for normalizing target gene expression is important for verifying expression changes. Metasequoia is a high-quality and economically important wood species. However, few systematic studies have examined reference genes in Metasequoia. Here, the expression stability of 14 candidate reference genes in different tissues and following different hormone treatments were analyzed using six algorithms. Candidate reference genes were used to normalize the expression pattern of FLOWERING LOCUS T and pyrabactin resistance-like 8. Analysis using the GrayNorm algorithm showed that ACT2 (Actin 2), HIS (histone superfamily protein H3) and TATA (TATA binding protein) were stably expressed in different tissues. ACT2, EF1α (elongation factor-1 alpha) and HIS were optimal for leaves treated with the flowering induction hormone solution, while Cpn60β (60-kDa chaperonin β-subunit), GAPDH (glyceraldehyde-3-phosphate dehydrogenase) and HIS were the best reference genes for treated buds. EF1α, HIS and TATA were useful reference genes for accurate normalization in abscisic acid-response signaling. Our results emphasize the importance of validating reference genes for qRT-PCR analysis in Metasequoia. To avoid errors, suitable reference genes should be used for different tissues and hormone treatments to increase normalization accuracy. Our study provides a foundation for reference gene normalization when analyzing gene expression in Metasequoia.

Development of 23 novel polymorphic EST-SSR markers for the endangered relict conifer Metasequoia glyptostroboides

PREMISE OF THE STUDY

Metasequoia glyptostroboides is an endangered relict conifer species endemic to China. In this study, expressed sequence tag-simple sequence repeat (EST-SSR) markers were developed using transcriptome mining for future genetic and functional studies.

METHODS AND RESULTS

We collected 97,565 unigene sequences generated by 454 pyrosequencing. A bioinformatics analysis identified 2087 unique and putative microsatellites, from which 96 novel microsatellite markers were developed. Fifty-three of the 96 primer sets successfully amplified clear fragments of the expected sizes; 23 of those loci were polymorphic. The number of alleles per locus ranged from two to eight, with an average of three, and the observed and expected heterozygosity values ranged from 0 to 1.0 and 0.117 to 0.813, respectively.

CONCLUSIONS

These microsatellite loci will enrich the genetic resources to develop functional studies and conservation strategies for this endangered relict species.

Characterization of the Transcriptome of the Xerophyte Ammopiptanthus mongolicus Leaves under Drought Stress by 454 Pyrosequencing

BACKGROUND

Ammopiptanthus mongolicus (Maxim. Ex Kom.) Cheng f., an endangered ancient legume species, endemic to the Gobi desert in north-western China. As the only evergreen broadleaf shrub in this area, A. mongolicus plays an important role in the region's ecological-environmental stability. Despite the strong potential of A. mongolicus in providing new insights on drought tolerance, sequence information on the species in public databases remains scarce. To both learn about the role of gene expression in drought stress tolerance in A. mongolicus and to expand genomic resources for the species, transcriptome sequencing of stress-treated A. mongolicus plants was performed.

RESULTS

Using 454 pyrosequencing technology, 8,480 and 7,474 contigs were generated after de novo assembly of RNA sequences from leaves of untreated and drought-treated plants, respectively. After clustering using TGICL and CAP3 programs, a combined assembly of all reads produced a total of 11,357 putative unique transcripts (PUTs). Functional annotation and classification of the transcripts were conducted by aligning the 11,357 PUTs against the public protein databases and nucleotide database (Nt). Between control and drought-treated plants, 1,620 differentially expressed genes (DEGs) were identified, of which 1,106 were up-regulated and 514 were down-regulated. The differential expression of twenty candidate genes in metabolic pathways and transcription factors families related to stress-response were confirmed by quantitative real-time PCR. Representatives of several large gene families, such as WRKY and P5CS, were identified and verified in A. mongolicus for the first time.

CONCLUSIONS

The additional transcriptome resources, gene expression profiles, functional annotations, and candidate genes provide a more comprehensive understanding of the stress response pathways in xeric-adapted plant species such as A. mongolicus.

Digital gene expression analysis of male and female bud transition in Metasequoia reveals high activity of MADS-box transcription factors and hormone-mediated sugar pathways

Metasequoia glyptostroboides is a famous redwood tree of ecological and economic importance, and requires more than 20 years of juvenile-to-adult transition before producing female and male cones. Previously, we induced reproductive buds using a hormone solution in juvenile Metasequoia trees as young as 5-to-7 years old. In the current study, hormone-treated shoots found in female and male buds were used to identify candidate genes involved in reproductive bud transition in Metasequoia. Samples from hormone-treated cone reproductive shoots and naturally occurring non-cone setting shoots were analyzed using 24 digital gene expression (DGE) tag profiles using Illumina, generating a total of 69,520 putative transcripts. Next, 32 differentially and specifically expressed transcripts were determined using quantitative real-time polymerase chain reaction, including the upregulation of MADS-box transcription factors involved in male bud transition and flowering time control proteins involved in female bud transition. These differentially expressed transcripts were associated with 243 KEGG pathways. Among the significantly changed pathways, sugar pathways were mediated by hormone signals during the vegetative-to-reproductive phase transition, including glycolysis/gluconeogenesis and sucrose and starch metabolism pathways. Key enzymes were identified in these pathways, including alcohol dehydrogenase (NAD) and glutathione dehydrogenase for the glycolysis/gluconeogenesis pathway, and glucanphosphorylase for sucrose and starch metabolism pathways. Our results increase our understanding of the reproductive bud transition in gymnosperms. In addition, these studies on hormone-mediated sugar pathways increase our understanding of the relationship between sugar and hormone signaling during female and male bud initiation in Metasequoia.

Populus euphratica: the transcriptomic response to drought stress

Populus euphratica Olivier is widely established in arid and semiarid regions but lags in the availability of transcriptomic resources in response to water deficiency. To investigate the mechanisms that allow P. euphratica to maintain growth in arid regions, the responses of the plant to soil water deficit were analyzed at a systems level using physiological and pyrosequencing approaches. We generated 218,601 and 287,120 reads from non-stressed control and drought-stressed P. euphratica leaves respectively, totaling over 200 million base pairs. After assembly, 24,013 transcripts were yielded with an average length of 1,128 bp. We determined 2,279 simple sequence repeats, which may have possible information for understanding drought adaption of woody plants. Stomatal closure was inhibited under moderate drought to maintain a relatively high rate of CO2 assimilation and water transportation, which was supposed to be important for P. euphratica to maintain normal growth and develop vigorous root systems in an adverse environment. This was accompanied by strong transcriptional remodeling of stress-perception, signaling and transcription regulation, photoprotective system, oxidative stress detoxification, and other stress responsive genes. In addition, genes involved in stomatal closure inhibition, ascorbate-glutathione pathway and ubiquitin-proteasome system that may specially modulate the drought stress responses of P. euphratica are highlighted. Our analysis provides a comprehensive picture of how P. euphratica responds to drought stress at physiological and transcriptome levels which may help to understand molecular mechanisms associated with drought response and could be useful for genetic engineering of woody plants.