Gene expression profiling of soybean leaves and roots under salt, saline-alkali and drought stress by high-throughput Illumina sequencing

Wild soybean roots depend on specific transcription factors and oxidation reduction related genesin response to alkaline stress

Soil alkalinity is an important environmental problem limiting agricultural productivity. Wild soybean (Glycine soja) shows strong alkaline stress tolerance, so it is an ideal plant candidate for studying the molecular mechanisms of alkaline tolerance and identifying alkaline stress-responsive genes. However, limited information is available about G. soja responses to alkaline stress on a genomic scale. Therefore, in the present study, we used RNA sequencing to compare transcript profiles of G. soja root responses to sodium bicarbonate (NaHCO3) at six time points, and a total of 68,138,478 pairs of clean reads were obtained using the Illumina GAIIX. Expression patterns of 46,404 G. soja genes were profiled in all six samples based on RNA-seq data using Cufflinks software. Then, t12 transcription factors from MYB, WRKY, NAC, bZIP, C2H2, HB, and TIFY families and 12 oxidation reduction related genes were chosen and verified to be induced in response to alkaline stress by using quantitative real-time polymerase chain reaction (qRT-PCR). The GO functional annotation analysis showed that besides "transcriptional regulation" and "oxidation reduction," these genes were involved in a variety of processes, such as "binding" and "response to stress." This is the first comprehensive transcriptome profiling analysis of wild soybean root under alkaline stress by RNA sequencing. Our results highlight changes in the gene expression patterns and identify a set of genes induced by NaHCO3 stress. These findings provide a base for the global analyses of G. soja alkaline stress tolerance mechanisms.

Mongolian Almond (Prunus mongolica Maxim): The Morpho-Physiological, Biochemical and Transcriptomic Response to Drought Stress

Prunus mongolica Maxim, which is widely established in the Gobi Desert, shows extreme tolerance to drought. However, there is a lack of available transcriptomic resources for this species related to its response to water deficiency. To investigate the mechanisms that allow P. mongolica to maintain growth in extremely arid environments, the response of P. mongolica seedlings to drought stress was analyzed using morphological, physiological, biochemical and high-throughput sequencing approaches. We generated 28,713,735 and 26,650,133 raw reads from no-stress control and drought-stressed P. mongolica seedlings, respectively. In total, we obtained 67,352 transcripts with an average length of 874.44 bp. Compared with the no-stress control, 3,365 transcripts were differentially expressed in the drought-stressed seedlings, including 55.75% (1,876 transcripts) up-regulated and 44.25% (1,489 transcripts) down-regulated transcripts. The photosynthesis response showed a decreasing tendency under drought stress, but the changes in the levels of hormones (auxins, cytokinins and abscisic acid) resulted in the closing of stomata and decreased cell enlargement and division; these changes were effective for promoting P. mongolica survival in Gobi Desert. Next, we analyzed the aquaporin and superoxide dismutase gene families due to their importance in plant resistance to drought stress. We found that all of the plasma membrane intrinsic protein transcripts were down-regulated in the drought-stressed treatment, whereas drought did not affect the expression of nodulin intrinsic protein or small basic intrinsic protein transcripts in P. mongolica seedlings. In addition, activation of iron superoxide dismutase transcription and enhanced transcription of manganese superoxide dismutase were observed in P. mongolica to promote tolerance of drought stress. This study identified drought response genes in P. mongolica seedlings. Our results provide a significant contribution to the understanding of how P. mongolica responds to drought stress at the transcriptome level, which may help to elucidate molecular mechanisms associated with the drought response of almond plants.

Transcriptome-based gene profiling provides novel insights into the characteristics of radish root response to Cr stress with next-generation sequencing

Radish (Raphanus sativus L.) is an important worldwide root vegetable crop with high nutrient values and is adversely affected by non-essential heavy metals including chromium (Cr). Little is known about the molecular mechanism underlying Cr stress response in radish. In this study, RNA-Seq technique was employed to identify differentially expressed genes (DEGs) under Cr stress. Based on de novo transcriptome assembly, there were 30,676 unigenes representing 60,881 transcripts isolated from radish root under Cr stress. Differential gene analysis revealed that 2985 uingenes were significantly differentially expressed between Cr-free (CK) and Cr-treated (Cr600) libraries, among which 1424 were up-regulated and 1561 down-regulated. Gene ontology (GO) analysis revealed that these DEGs were mainly involved in primary metabolic process, response to abiotic stimulus, cellular metabolic process and small molecule metabolic process. Kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis showed that the DEGs were mainly involved in protein processing in endoplasmic reticulum, starch and sucrose metabolism, amino acid metabolism, glutathione metabolism, drug and xenobiotics by cytochrome P450 metabolism. RT-qPCR analysis showed that the expression patterns of 12 randomly selected DEGs were highly accordant with the results from RNA-seq. Furthermore, many candidate genes including signaling protein kinases, transcription factors and metal transporters, chelate compound biosynthesis and antioxidant system, were involved in defense and detoxification mechanisms of Cr stress response regulatory networks. These results would provide novel insight into molecular mechanism underlying plant responsiveness to Cr stress and facilitate further genetic manipulation on Cr uptake and accumulation in radish.

Response to drought and salt stress in leaves of poplar (Populus alba × Populus glandulosa): expression profiling by oligonucleotide microarray analysis

Drought and salt stresses are major environmental constraints on forest productivity. To identify genes responsible for stress tolerance, we conducted a genome-wide analysis in poplar (Populus alba × Populus glandulosa) leaves exposed to drought and salt (NaCl) stresses. We investigated gene expression at the mRNA level using oligonucleotide microarrays containing 44,718 genes from Populus trichocarpa. A total of 1604 and 1042 genes were up-regulated (≥2-fold; P value < 0.05) by drought and salt stresses, respectively, and 765 genes were up-regulated by both stresses. In addition, 2742 and 1685 genes were down-regulated by drought and salt stresses, respectively, and 1564 genes were down-regulated by both stresses. The large number of genes regulated by both stresses suggests that crosstalk occurs between the drought and salt stress responses. Most up-regulated genes were involved in functions such as subcellular localization, signal transduction, metabolism, and transcription. Among the up-regulated genes, we identified 47 signaling proteins, 65 transcription factors, and 43 abiotic stress-related genes. Several genes were modulated by only one of the two stresses. About 25% of the genes significantly regulated by these stresses are of unknown function, suggesting that poplar may provide an opportunity to discover novel stress-related genes.

Transcriptomic complexity in young maize primary roots in response to low water potentials

BACKGROUND

Widespread and more frequently occurring drought conditions are a consequence of global warming and increase the demand for tolerant crop varieties to feed the growing world population. A better understanding of the molecular mechanisms underlying the water deficit response of crops will enable targeted breeding strategies to develop robust cultivars.

RESULTS

In the present study, the transcriptional response of maize (Zea mays L.) primary roots to low water potentials was monitored by RNA sequencing (RNA-Seq) experiments. After 6 h and 24 h of mild (-0.2 MPa) and severe (-0.8 MPa) water deficit conditions, the primary root transcriptomes of seedlings grown under water deficit and control conditions were compared. The number of responsive genes was dependent on and increased with intensification of water deficit treatment. After short-term mild and severe water deficit 249 and 3,000 genes were differentially expressed, respectively. After a 24 h treatment the number of affected genes increased to 7,267 and 12,838 for mild and severe water deficit, respectively, including more than 80% of the short-term responsive genes. About half of the differentially expressed genes were up-regulated and maximal fold-changes increased with treatment intensity to more than 300-fold. A consensus set of 53 genes was differentially regulated independently of the nature of deficit treatment. Characterization revealed an overrepresentation of the Gene Ontology (GO) categories "oxidoreductase activity" and "heme binding" among regulated genes connecting the water deficit response to ROS metabolism.

CONCLUSION

This study gives a comprehensive insight in water deficit responsive genes in young maize primary roots and provides a set of candidate genes that merit further genetic analyses in the future.

Genome-wide transcriptome analysis of genes involved in flavonoid biosynthesis between red and white strains of Magnolia sprengeri pamp

BACKGROUND

Magnolia sprengeri Pamp is one of the most highly valuable medicinal and ornamental plants of the Magnolia Family. The natural color of M. sprengeri is variable. The complete genome sequence of M. sprengeri is not available; therefore we sequenced the transcriptome of white and red petals of M. sprengeri using Illumina technology. We focused on the identity of structural and regulatory genes encoding the enzymes involved in the determination of flower color.

RESULTS

We sequenced and annotated a reference transcriptome for M. sprengeri, and aimed to capture the transcriptional determinanats of flower color. We sequenced a normalized cDNA library of white and red petals using Illumina technology. The resulting reads were assembled into 77,048 unique sequences, of which 28,243 could be annotated by Gene Ontology (GO) analysis, while 48,805 transcripts lacked GO annotation. The main enzymes involved in the flavonoid biosynthesis, such as phenylalanine ammonia-Lyase, cinnamat-4-Hydroxylase, dihydroflavonol-4-reductase, flavanone 3-hydroxylase, flavonoid-3'-hydroxylase, flavonol synthase, chalcone synthase and anthocyanidin synthase, were identified in the transcriptome. A total of 270 transcription factors were sorted into three families, including MYB, bHLH and WD40 types. Among these transcription factors, eight showed 4-fold or greater changes in transcript abundance in red petals compared with white petals. High-performance liquid chromatography analysis of anthocyanin compositions showed that the main anthocyanin in the petals of M. sprengeri is cyanidin-3-O-glucoside chloride and its content in red petals was 26-fold higher than that in white petals.

CONCLUSION

This study presents the first next-generation sequencing effort and transcriptome analysis of a non-model plant from the Family Magnoliaceae. Genes encoding key enzymes were identified and the metabolic pathways involved in biosynthesis and catabolism of M. sprengeri flavonoids were reconstructed. Identification of these genes and pathways adds to the current knowledge of the molecular biology and biochemistry of their production in plant. Such insights into the mechanisms supporting metabolic processes could be used to genetically to enhance flower color among members of the Magnoliaceae.

De novo sequencing, assembly, and analysis of the Taxodium 'Zhongshansa' roots and shoots transcriptome in response to short-term waterlogging

BACKGROUND

Taxodium is renowned for its strong tolerance to waterlogging stress, thus it has great ecological and economic potential. However, the scant genomic resources in genus Taxodium have greatly hindered further exploration of its underlying flood-tolerance mechanism. Taxodium 'Zhongshansa' is an interspecies hybrid of T. distichum and T. mucronatum, and has been widely planted in southeastern China. To understand the genetic basis of its flood tolerance, we analyzed the transcriptomes of Taxodium 'Zhongshansa' roots and shoots in response to short-term waterlogging.

RESULTS

RNA-seq was used to analyze genome-wide transcriptome changes of Taxodium 'Zhongshansa 406' clone root and shoot treated with 1 h of soil-waterlogging stress. After de novo assembly, 108,692 unigenes were achieved, and 70,260 (64.64%) of them were annotated. There were 2090 differentially expressed genes (DEGs) found in roots and 394 in shoots, with 174 shared by both of them, indicating that the aerial parts were also affected. Under waterlogging stress, the primary reaction of hypoxic-treated root was to activate the antioxidative defense system to prevent cells experiencing reactive oxygen species (ROS) poisoning. As respiration was inhibited and ATP decreased, another quick coping mechanism was repressing the energy-consuming biosynthetic processes through the whole plant. The glycolysis and fermentation pathway was activated to maintain ATP production in the hypoxic root. Constantly, the demand for carbohydrates increased, and carbohydrate metabolism were accumulated in the root as well as the shoot, possibly indicating that systemic communications between waterlogged and non-waterlogged tissues facilated survival. Amino acid metabolism was also greatly influenced, with down-regulation of genes involvedin serine degradation and up-regulation of aspartic acid degradation. Additionally, a non-symbiotic hemoglobin class 1 gene was up-regulated, which may also help the ATP production. Moreover, the gene expression pattern of 5 unigenes involving in the glycolysis pathway revealed by qRT-PCR confirmed the RNA-Seq data.

CONCLUSIONS

We conclude that ROS detoxification and energy maintenance were the primary coping mechanisms of 'Zhongshansa' in surviving oxygen deficiency, which may be responsible for its remarkable waterlogging tolerance. Our study not only provided the first large-scale assessment of genomic resources of Taxodium but also guidelines for probing the molecular mechanism underlying 'Zhongshansa' waterlogging tolerance.

Comparative ecophysiological study of salt stress for wild and cultivated soybean species from the Yellow River Delta, China

Osmotic and ionic stresses were the primary and instant damage produced by salt stress. They can also bring about other secondary stresses. Soybean is an important economic crop and the wild soybean aroused increasing attention for its excellent performance in salt resistance. For this reason, we compared the different performances of Glycine max L. (ZH13) and Glycine soja L. (BB52) in both young and mature seedlings, hoping to clarify the specific reasons. Our research revealed that, compared to the cultivated soybean, the wild soybean was able to maintain higher water potential and relative water content (RWC), accumulate more amount of proline and glycine betaine, reduce the contents of Na(+) and Cl(-) by faster efflux, and cut down the efflux of the K(+) as well as keep higher K(+)/Na(+) ratio. And what is more is that, almost all the excel behaviors became particularly obvious under higher NaCl concentration (300 mM). Therefore, according to all the detections and comparisons, we concluded that the wild soybean had different tolerance mechanisms and better salt resistance. It should be used as eminent germplasm resource to enhance the resistant ability of cultivated soybean or even other crops.

Identification of differentially expressed genes in leaf of Reaumuria soongorica under PEG-induced drought stress by digital gene expression profiling

Reaumuria soongorica (Pall.) Maxim., a resurrection semi-shrub, is a typical constructive and dominant species in desert ecosystems in northwestern China. However, the gene expression characteristics of R. soongorica under drought stress have not been elucidated. Digital gene expression analysis was performed using Illumina technique to investigate differentially expressed genes (DEGs) between control and PEG-treated samples of R. soongorica. A total of 212,338 and 211,052 distinct tags were detected in the control and PEG-treated libraries, respectively. A total of 1,325 genes were identified as DEGs, 379 (28.6%) of which were up-regulated and 946 (71.4%) were down-regulated in response to drought stress. Functional annotation analysis identified numerous drought-inducible genes with various functions in response to drought stress. A number of regulatory proteins, functional proteins, and proteins induced by other stress factors in R. soongorica were identified. Alteration in the regulatory proteins (transcription factors and protein kinase) may be involved in signal transduction. Functional proteins, including flavonoid biosynthetic proteins, late embryogenesis abundant (LEA) proteins, small heat shock proteins (sHSP), and aquaporin and proline transporter may play protective roles in response to drought stress. Flavonoids, LEA proteins and sHSP function as reactive oxygen species scavenger or molecular chaperone. Aquaporin and proline transporters regulate the distribution of water and proline throughout the whole plant. The tolerance ability of R. soongorica may be gained through effective signal transduction and enhanced protection of functional proteins to reestablish cellular homeostasis. DEGs obtained in this study may provide useful insights to help further understand the drought-tolerant mechanism of R. soongorica.

Transcriptome analysis of salt tolerant common bean (Phaseolus vulgaris L.) under saline conditions

Salinity is one of the important abiotic stress factors that limit crop production. Common bean, Phaseolus vulgaris L., a major protein source in developing countries, is highly affected by soil salinity and the information on genes that play a role in salt tolerance is scarce. We aimed to identify differentially expressed genes (DEGs) and related pathways by comprehensive analysis of transcriptomes of both root and leaf tissues of the tolerant genotype grown under saline and control conditions in hydroponic system. We have generated a total of 158 million high-quality reads which were assembled into 83,774 all-unigenes with a mean length of 813 bp and N50 of 1,449 bp. Among the all-unigenes, 58,171 were assigned with Nr annotations after homology analyses. It was revealed that 6,422 and 4,555 all-unigenes were differentially expressed upon salt stress in leaf and root tissues respectively. Validation of the RNA-seq quantifications (RPKM values) was performed by qRT-PCR (Quantitative Reverse Transcription PCR) analyses. Enrichment analyses of DEGs based on GO and KEGG databases have shown that both leaf and root tissues regulate energy metabolism, transmembrane transport activity, and secondary metabolites to cope with salinity. A total of 2,678 putative common bean transcription factors were identified and classified under 59 transcription factor families; among them 441 were salt responsive. The data generated in this study will help in understanding the fundamentals of salt tolerance in common bean and will provide resources for functional genomic studies.

Early transcriptional response of soybean contrasting accessions to root dehydration

Drought is a significant constraint to yield increase in soybean. The early perception of water deprivation is critical for recruitment of genes that promote plant tolerance. DeepSuperSAGE libraries, including one control and a bulk of six stress times imposed (from 25 to 150 min of root dehydration) for drought-tolerant and sensitive soybean accessions, allowed to identify new molecular targets for drought tolerance. The survey uncovered 120,770 unique transcripts expressed by the contrasting accessions. Of these, 57,610 aligned with known cDNA sequences, allowing the annotation of 32,373 unitags. A total of 1,127 unitags were up-regulated only in the tolerant accession, whereas 1,557 were up-regulated in both as compared to their controls. An expression profile concerning the most representative Gene Ontology (GO) categories for the tolerant accession revealed the expression "protein binding" as the most represented for "Molecular Function", whereas CDPK and CBL were the most up-regulated protein families in this category. Furthermore, particular genes expressed different isoforms according to the accession, showing the potential to operate in the distinction of physiological behaviors. Besides, heat maps comprising GO categories related to abiotic stress response and the unitags regulation observed in the expression contrasts covering tolerant and sensitive accessions, revealed the unitags potential for plant breeding. Candidate genes related to "hormone response" (LOX, ERF1b, XET), "water response" (PUB, BMY), "salt stress response" (WRKY, MYB) and "oxidative stress response" (PER) figured among the most promising molecular targets. Additionally, nine transcripts (HMGR, XET, WRKY20, RAP2-4, EREBP, NAC3, PER, GPX5 and BMY) validated by RT-qPCR (four different time points) confirmed their differential expression and pointed that already after 25 minutes a transcriptional reorganization started in response to the new condition, with important differences between both accessions.

Transcriptome profiling of radish (Raphanus sativus L.) root and identification of genes involved in response to Lead (Pb) stress with next generation sequencing

Lead (Pb), one of the most toxic heavy metals, can be absorbed and accumulated by plant roots and then enter the food chain resulting in potential health risks for human beings. The radish (Raphanus sativus L.) is an important root vegetable crop with fleshy taproots as the edible parts. Little is known about the mechanism by which radishes respond to Pb stress at the molecular level. In this study, Next Generation Sequencing (NGS)-based RNA-seq technology was employed to characterize the de novo transcriptome of radish roots and identify differentially expressed genes (DEGs) during Pb stress. A total of 68,940 assembled unique transcripts including 33,337 unigenes were obtained from radish root cDNA samples. Based on the assembled de novo transcriptome, 4,614 DEGs were detected between the two libraries of untreated (CK) and Pb-treated (Pb1000) roots. Gene Ontology (GO) and pathway enrichment analysis revealed that upregulated DEGs under Pb stress are predominately involved in defense responses in cell walls and glutathione metabolism-related processes, while downregulated DEGs were mainly involved in carbohydrate metabolism-related pathways. The expression patterns of 22 selected genes were validated by quantitative real-time PCR, and the results were highly accordant with the Solexa analysis. Furthermore, many candidate genes, which were involved in defense and detoxification mechanisms including signaling protein kinases, transcription factors, metal transporters and chelate compound biosynthesis related enzymes, were successfully identified in response to heavy metal Pb. Identification of potential DEGs involved in responses to Pb stress significantly reflected alterations in major biological processes and metabolic pathways. The molecular basis of the response to Pb stress in radishes was comprehensively characterized. Useful information and new insights were provided for investigating the molecular regulation mechanism of heavy metal Pb accumulation and tolerance in root vegetable crops.