Genome-wide expression profiling of the transcriptomes of four Paulownia tomentosa accessions in response to drought

Genome-wide expression analysis of salt-stressed diploid and autotetraploid Paulownia tomentosa

Paulownia tomentosa is a fast-growing tree species with multiple uses. It is grown worldwide, but is native to China, where it is widely cultivated in saline regions. We previously confirmed that autotetraploid P. tomentosa plants are more stress-tolerant than the diploid plants. However, the molecular mechanism underlying P. tomentosa salinity tolerance has not been fully characterized. Using the complete Paulownia fortunei genome as a reference, we applied next-generation RNA-sequencing technology to analyze the effects of salt stress on diploid and autotetraploid P. tomentosa plants. We generated 175 million clean reads and identified 15,873 differentially expressed genes (DEGs) from four P. tomentosa libraries (two diploid and two autotetraploid). Functional annotations of the differentially expressed genes using the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes databases revealed that plant hormone signal transduction and photosynthetic activities are vital for plant responses to high-salt conditions. We also identified several transcription factors, including members of the AP2/EREBP, bHLH, MYB, and NAC families. Quantitative real-time PCR analysis validated the expression patterns of eight differentially expressed genes. Our findings and the generated transcriptome data may help to accelerate the genetic improvement of cultivated P. tomentosa and other plant species for enhanced growth in saline soils.

Comparative proteomic analysis of autotetraploid and diploid Paulownia tomentosa reveals proteins associated with superior photosynthetic characteristics and stress adaptability in autotetraploid Paulownia

To enlarge the germplasm resource of Paulownia plants, we used colchicine to induce autotetraploid Paulownia tomentosa, as reported previously. Compared with its diploid progenitor, autotetraploid P. tomentosa exhibits better photosynthetic characteristics and higher stress resistance. However, the underlying mechanism for its predominant characteristics has not been determined at the proteome level. In this study, isobaric tag for relative and absolute quantitation coupled with liquid chromatography-tandem mass spectrometry was employed to compare proteomic changes between autotetraploid and diploid P. tomentosa. A total of 1427 proteins were identified in our study, of which 130 proteins were differentially expressed between autotetraploid and diploid P. tomentosa. Functional analysis of differentially expressed proteins revealed that photosynthesis-related proteins and stress-responsive proteins were significantly enriched among the differentially expressed proteins, suggesting they may be responsible for the photosynthetic characteristics and stress adaptability of autotetraploid P. tomentosa. The correlation analysis between transcriptome and proteome data revealed that only 15 (11.5%) of the differentially expressed proteins had corresponding differentially expressed unigenes between diploid and autotetraploid P. tomentosa. These results indicated that there was a limited correlation between the differentially expressed proteins and the previously reported differentially expressed unigenes. This work provides new clues to better understand the superior traits in autotetraploid P. tomentosa and lays a theoretical foundation for developing Paulownia breeding strategies in the future.

Application of Genomic Technologies to the Breeding of Trees

The recent introduction of next generation sequencing (NGS) technologies represents a major revolution in providing new tools for identifying the genes and/or genomic intervals controlling important traits for selection in breeding programs. In perennial fruit trees with long generation times and large sizes of adult plants, the impact of these techniques is even more important. High-throughput DNA sequencing technologies have provided complete annotated sequences in many important tree species. Most of the high-throughput genotyping platforms described are being used for studies of genetic diversity and population structure. Dissection of complex traits became possible through the availability of genome sequences along with phenotypic variation data, which allow to elucidate the causative genetic differences that give rise to observed phenotypic variation. Association mapping facilitates the association between genetic markers and phenotype in unstructured and complex populations, identifying molecular markers for assisted selection and breeding. Also, genomic data provide in silico identification and characterization of genes and gene families related to important traits, enabling new tools for molecular marker assisted selection in tree breeding. Deep sequencing of transcriptomes is also a powerful tool for the analysis of precise expression levels of each gene in a sample. It consists in quantifying short cDNA reads, obtained by NGS technologies, in order to compare the entire transcriptomes between genotypes and environmental conditions. The miRNAs are non-coding short RNAs involved in the regulation of different physiological processes, which can be identified by high-throughput sequencing of RNA libraries obtained by reverse transcription of purified short RNAs, and by in silico comparison with known miRNAs from other species. All together, NGS techniques and their applications have increased the resources for plant breeding in tree species, closing the former gap of genetic tools between trees and annual species.

Drought stress-induced changes of microRNAs in diploid and autotetraploid Paulownia tomentosa

Drought stress adversely affects plant productivity. Growth and timber production of Paulownia trees are limited under drought stress. Changes in gene expression patterns and miRNA in different ploidy of Paulownia tomentosa have been investigated. However, the responses of P. tomentosa to drought stress at the microRNA (miRNA) level have not been reported so far. To identify miRNA candidates and their target genes involved in the drought stress response in diploid and tetraploid P. tomentosa, four small RNA and four degradome libraries from diploid and autotetraploid P. tomentosa under normal and drought stress conditions were constructed and sequenced. A total of 41 conserved and 90 novel miRNAs were identified. Among these miRNAs, 67 (26 conserved and 41 novel) and 53 (six conserved and 47 novel) were significantly differentially expressed in response to drought stress in diploid and autotetraploid P. tomentosa, respectively. Degradome analysis identified 356 candidate miRNA target genes that encoded proteins with functions that included plant defense, transcriptional regulation, and hormone metabolism. In particular, miR4 and miR156 were identified only in autotetraploid P. tomentosa under drought stress. These results will help us build a foundation for future studies of the biological functions of miRNA-mediated gene regulation in P. tomentosa.

Exploring drought stress-regulated genes in senna (Cassia angustifolia Vahl.): a transcriptomic approach

De novo assembly of reads produced by next-generation sequencing (NGS) technologies offers a rapid approach to obtain expressed gene sequences for non-model organisms. Senna (Cassia angustifolia Vahl.) is a drought-tolerant annual undershrub of Caesalpiniaceae, a subfamily of Fabaceae. There are insufficient transcriptomic and genomic data in public databases for understanding the molecular mechanism underlying the drought tolerance of senna. Therefore, the main purpose of this study was to know the transcriptome profile of senna, with special reference to drought stress. RNA from two different stages of leaf development was extracted and sequenced separately using the Illumina technology. A total of 200 million reads were generated, and a de novo assembly of processed reads in the pooled transcriptome using Trinity yielded 43,413 transcripts which were further annotated using NCBI BLAST with "green plant database (txid 33090)," Swiss Prot, Kyoto Encyclopedia of Genes and Genomes (KEGG), Clusters of Orthologous Groups (COG), and Gene Ontology (GO). Out of the total transcripts, 42,280 (95.0 %) were annotated by BLASTX against the green plant database of NCBI. Senna transcriptome showed the highest similarity to Glycine max (41 %), followed by Phaseolus vulgaris (16 %), Cicer arietinum (15 %), and Medicago trancatula (5 %). The highest number of GO terms were enriched for the molecular functions category; of these "catalytic activity" (GO: 0003824) (25.10 %) and "binding activity" (GO: 0005488) (20.10 %) were most abundantly represented. We used InterProscan to see protein similarity at domain level; a total of 33,256 transcripts were annotated against the Pfam domains. The transcripts were assigned with various KEGG pathways. Coding DNA sequences (CDS) encoding various drought stress-regulated pathways such as signaling factors, protein-modifying/degrading enzymes, biosynthesis of phytohormone, phytohormone signaling, osmotically active compounds, free radical scavengers, chlorophyll metabolism, leaf cuticular wax, polyamines, and protective proteins were identified through BLASTX search. The lucine-rich repeat kinase family was the most abundantly found group of protein kinases. Orphan, bHLH, and bZIP family TFs were the most abundantly found in senna. Six genes encoding MYC2 transcription factor, 9-cis-epoxycarotenoid dioxygenase (NCED), l -ascorbate peroxidase (APX), aminocyclopropane carboxylate oxidase (ACO), abscisic acid 8'-hydroxylase (ABA), and WRKY transcription factor were confirmed through reverse transcriptase-PCR (RT-PCR) and Sanger sequencing for the first time in senna. The potential drought stress-related transcripts identified in this study provide a good start for further investigation into the drought adaptation in senna. Additionally, our transcriptome sequences are the valuable resource for accelerated genomics-assisted genetic improvement programs and facilitate manipulation of biochemical pathways for developing drought-tolerant genotypes of crop plants.

Quantitative Proteomic and Transcriptomic Study on Autotetraploid Paulownia and Its Diploid Parent Reveal Key Metabolic Processes Associated with Paulownia Autotetraploidization

Polyploidy plays a very important role in speciation and plant evolution by way of genomic merging and doubling. In the process of polyploidy, rapid genomic, and transcriptomic changes have been observed and researched. However, proteomic divergence caused by the effects of polyploidization is still poorly understood. In the present study, we used iTRAQ coupled with mass spectrometry to quantitatively analyze proteomic changes in the leaves of autotetraploid Paulownia and its diploid parent. A total of 2963 proteins were identified and quantified. Among them, 463 differentially abundant proteins were detected between autotetraploid Paulownia and its diploid parent, and 198 proteins were found to be non-additively abundant in autotetraploid Paulownia, suggesting the presence of non-additive protein regulation during genomic merger and doubling. We also detected 1808 protein-encoding genes in previously published RNA sequencing data. We found that 59 of the genes that showed remarkable changes at mRNA level encoded proteins with consistant changes in their abundance levels, while a further 48 genes that showed noteworthy changes in their expression levels encoded proteins with opposite changes in their abundance levels. Proteins involved in posttranslational modification, protein turnover, and response to stimulus, were significantly enriched among the non-additive proteins, which may provide some of the driving power for variation and adaptation in autopolyploids. Quantitative real-time PCR analysis verified the expression patterns of related protein-coding genes. In addition, we found that the percentage of differentially abundant proteins that matched previously reported differentially expressed genes was relatively low.

Changes in Transcript Related to Osmosis and Intracellular Ion Homeostasis in Paulownia tomentosa under Salt Stress

Paulownia tomentosa is an important economic and greening tree species that is cultivated widely, including salt environment. Our previous studies indicated its autotetraploid induced by colchicine showed better stress tolerance, but the underlying molecular mechanism related to ploidy and salt stress is still unclear. To investigate this issue, physiological measurements and transcriptome profiling of diploid and autotetraploid plants untreated and treated with NaCl were performed. Through the comparisons among four accessions, for one thing, we found different physiological changes between diploid and autotetraploid P. tomentosa; for another, and we detected many differentially expressed unigenes involved in salt stress response. These differentially expressed unigenes were assigned to several metabolic pathways, including "plant hormone signal transduction," "RNA transporter," "protein processing in endoplasmic reticulum," and "plant-pathogen interaction," which constructed the complex regulatory network to maintain osmotic and intracellular ion homeostasis. Quantitative real-time polymerase chain reaction was used to confirm the expression patterns of 20 unigenes. The results establish the foundation for the genetic basis of salt tolerance in P. tomentosa, which in turn accelerates Paulownia breeding and expands available arable land.

Comparative Analysis and Identification of miRNAs and Their Target Genes Responsive to Salt Stress in Diploid and Tetraploid Paulownia fortunei Seedlings

Salt stress is a global environmental problem that affects plant growth and development. Paulownia fortunei is an adaptable and fast-growing deciduous tree native to China that is environmentally and economically important. MicroRNAs (miRNAs) play important regulatory roles in growth, development, and stress responses in plants. MiRNAs that respond to biotic stresses have been identified; however, how miRNAs in P. fortunei respond to salt stress has not yet been reported. To identify salt-stress-responsive miRNAs and predict their target genes, four small RNA and four degradome libraries were constructed from NaCl-treated and NaCl-free leaves of P. fortunei seedlings. The results indicated that salt stress had different physiological effects on diploid and tetraploid P. fortunei. We detected 53 conserved miRNAs belonging to 17 miRNA families and 134 novel miRNAs in P. fortunei. Comparing their expression levels in diploid and tetraploid P. fortunei, we found 10 conserved and 10 novel miRNAs that were significantly differentially expressed under salt treatment, among them eight were identified as miRNAs probably associated with higher salt tolerance in tetraploid P. fortunei than in diploid P. fortunei. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses were performed to predict the functions of the target genes of the conserved and novel miRNAs. The expressions of 10 differentially expressed miRNAs were validated by quantitative real-time polymerase chain reaction (qRT-PCR). This is the first report on P. fortunei miRNAs and their target genes under salt stress. The results provided information at the physiological and molecular levels for further research into the response mechanisms of P. fortunei to salt stress.

Transcriptional Responses in Root and Leaf of Prunus persica under Drought Stress Using RNA Sequencing

Prunus persica L. Batsch, or peach, is one of the most important crops and it is widely established in irrigated arid and semi-arid regions. However, due to variations in the climate and the increased aridity, drought has become a major constraint, causing crop losses worldwide. The use of drought-tolerant rootstocks in modern fruit production appears to be a useful method of alleviating water deficit problems. However, the transcriptomic variation and the major molecular mechanisms that underlie the adaptation of drought-tolerant rootstocks to water shortage remain unclear. Hence, in this study, high-throughput sequencing (RNA-seq) was performed to assess the transcriptomic changes and the key genes involved in the response to drought in root tissues (GF677 rootstock) and leaf tissues (graft, var. Catherina) subjected to 16 days of drought stress. In total, 12 RNA libraries were constructed and sequenced. This generated a total of 315 M raw reads from both tissues, which allowed the assembly of 22,079 and 17,854 genes associated with the root and leaf tissues, respectively. Subsets of 500 differentially expressed genes (DEGs) in roots and 236 in leaves were identified and functionally annotated with 56 gene ontology (GO) terms and 99 metabolic pathways, which were mostly associated with aminobenzoate degradation and phenylpropanoid biosynthesis. The GO analysis highlighted the biological functions that were exclusive to the root tissue, such as "locomotion," "hormone metabolic process," and "detection of stimulus," indicating the stress-buffering role of the GF677 rootstock. Furthermore, the complex regulatory network involved in the drought response was revealed, involving proteins that are associated with signaling transduction, transcription and hormone regulation, redox homeostasis, and frontline barriers. We identified two poorly characterized genes in P. persica: growth-regulating factor 5 (GRF5), which may be involved in cellular expansion, and AtHB12, which may be involved in root elongation. The reliability of the RNA-seq experiment was validated by analyzing the expression patterns of 34 DEGs potentially involved in drought tolerance using quantitative reverse transcription polymerase chain reaction. The transcriptomic resources generated in this study provide a broad characterization of the acclimation of P. persica to drought, shedding light on the major molecular responses to the most important environmental stressor.