Transcriptome characterization and differential expression analysis of cold-responsive genes in young spikes of common wheat

Integration of the transcriptome and metabolome reveals the mechanism of resistance to low phosphorus in wild soybean seedling leaves

Plant growth, development, yield and quality are limited by barren soil. Soil phosphorus deficiency is one of the common factors causing soil barrenness. Plants have evolved morphological, physiological and molecular adaptations to resist to phosphorus deficiency. Wild soybean, a wild relative of cultivated soybean, has an obvious genetic relationship with cultivated soybean and has many beneficial characteristics such as strong low phosphorus resistance. Therefore, in this study, the integration analysis of transcriptome and metabolome of wild and cultivated soybean seedlings leaves were applied under phosphorus deficiency to reveal the mechanism of resistance to low phosphorus stress in wild soybean leaves, especially the key role of membrane phospholipid reuse and protection. Under phosphorus deficiency, wild soybean resisted low phosphorus stress by enhancing phosphorus reuse and strengthening membrane protection mechanisms, that is, by enhancing phospholipid metabolism, degrading membrane phospholipids, releasing phosphorus, increasing phosphorus reuse, and enhancing galactolipid biosynthesis. This, in turn, produced digalactosyl diacylglycerol to replace missing phospholipids for membrane maintenance and enhanced glutathione metabolism to protect the membrane system from damage. At the same time, phosphorus deficiency increased the levels of the intermediate metabolites glycine and ornithine, while significantly regulating the expression of transcription factors WRKY75 and MYB86. The enhancement of these metabolic pathways and the significant regulation of gene expression play an important role in improving the low phosphorus tolerance of wild soybean. This study will provide a useful theoretical basis for breeding soybean with low phosphorus tolerance.

Physiology and Molecular Breeding in Sustaining Wheat Grain Setting and Quality under Spring Cold Stress

Spring cold stress (SCS) compromises the reproductive growth of wheat, being a major constraint in achieving high grain yield and quality in winter wheat. To sustain wheat productivity in SCS conditions, breeding cultivars conferring cold tolerance is key. In this review, we examine how grain setting and quality traits are affected by SCS, which may occur at the pre-anthesis stage. We have investigated the physiological and molecular mechanisms involved in floret and spikelet SCS tolerance. It includes the protective enzymes scavenging reactive oxygen species (ROS), hormonal adjustment, and carbohydrate metabolism. Lastly, we explored quantitative trait loci (QTLs) that regulate SCS for identifying candidate genes for breeding. The existing cultivars for SCS tolerance were primarily bred on agronomic and morphophysiological traits and lacked in molecular investigations. Therefore, breeding novel wheat cultivars based on QTLs and associated genes underlying the fundamental resistance mechanism is urgently needed to sustain grain setting and quality under SCS.

Comparative Physiology and Transcriptome Analysis of Young Spikes in Response to Late Spring Coldness in Wheat (Triticum aestivum L.)

Late spring coldness (LSC) is critical for wheat growth and development in the Huang-Huai valleys of China. However, little is known about the molecular mechanisms for young spikes responding to low temperature (LT) stress during anther connective tissue formation phase (ACFP). To elucidate the molecular mechanisms associated with low temperature, we performed a comparative transcriptome analysis of wheat cultivars Xinmai26 (XM26: cold-sensitive) and Yannong19 (YN19: cold-tolerant) using RNA-seq data. Over 4000 differently expressed genes (DEGs) were identified under low temperature conditions (T1: 4°C) and freezing conditions (T2: -4°C) compared with control (CK: 16°C). The number of DEGs associated with two cultivars at two low temperature treatments (T1: 4°C and T2: -4°C) were 834, 1,353, 231, and 1,882 in four comparison groups (Xinmai26-CK vs. Xinmai26-T1, Xinmai26-CK vs. Xinmai26-T2, Yannong19-CK vs. Yannong19-T1, and Yannong19-CK vs. Yannong19-T2), respectively. Furthermore, to validate the accuracy of RNA-seq, 16 DEGs were analyzed using quantitative real-time RT-PCR. Several transcriptome changes were observed through Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway functional enrichment analysis in plant hormone signal transduction, circadian rhythm-plant, and starch and sucrose metabolism under low temperature. In addition, 126 transcription factors (TFs), including AP2-ERF, bHLH, WRKY, MYB, HSF, and members of the bZIP family, were considered as cold-responsive. It is the first study to investigate DEGs associated with low temperature stress at the transcriptome level in two wheat cultivars with different cold resistance capacities. Most likely, the variations in transcription factors (TFs) regulation, and starch and sucrose metabolism contribute to different cold resistance capacities in the two cultivars. Further, physiological activities of superoxide dismutase (SOD), peroxidase (POD), catalase (CAT) enzymes, malondialdehyde (MDA), soluble sugar (SS), and sucrose contents were evaluated to investigate the negative impacts of low temperature in both cultivars. These findings provide new insight into the molecular mechanisms of plant responses to low temperature and potential candidate genes that required for improving wheat's capacity to withstand low temperature stress.

Comparative Transcriptome Analyses Revealed Conserved and Novel Responses to Cold and Freezing Stress in Brassica napus L

Oil rapeseed (Brassica napus L.) is a typical winter biennial plant, with high cold tolerance during vegetative stage. In recent years, more and more early-maturing rapeseed varieties were planted across China. Unfortunately, the early-maturing rapeseed varieties with low cold tolerance have higher risk of freeze injury in cold winter and spring. Little is known about the molecular mechanisms for coping with different low-temperature stress conditions in rapeseed. In this study, we investigated 47,328 differentially expressed genes (DEGs) of two early-maturing rapeseed varieties with different cold tolerance treated with cold shock at chilling (4°) and freezing (−4°) temperatures, as well as chilling and freezing stress following cold acclimation or control conditions. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated that two conserved (the primary metabolism and plant hormone signal transduction) and two novel (plant-pathogen interaction pathway and circadian rhythms pathway) signaling pathways were significantly enriched with differentially-expressed transcripts. Our results provided a foundation for understanding the low-temperature stress response mechanisms of rapeseed. We also propose new ideas and candidate genes for genetic improvement of rapeseed tolerance to cold stresses.

Transcriptomic response of durum wheat to cold stress at reproductive stage

Understanding the genetic basis of cold tolerance is a key step towards obtaining new and improved crop varieties. Current geographical distribution of durum wheat in Argentina exposes the plants to frost damage when spikes have already emerged. Biochemical pathways involved in cold tolerance are known to be early activated at above freezing temperatures. In this study we reported the transcriptome of CBW0101 spring durum wheat by merging data from untreated control and cold (5 °C) treated plant samples at reproductive stage. A total of 128,804 unigenes were predicted. Near 62% of the unigenes were annotated in at least one database. In total 876 unigenes were differentially expressed (DEGs), 562 were up-regulated and 314 down-regulated in treated samples. DEGs are involved in many critical processes including, photosynthetic activity, lipid and carbohydrate synthesis and accumulation of amino acids and seed proteins. Twenty-eight transcription factors (TFs) belonging to 14 families resulted differentially expressed from which eight families comprised of only TFs induced by cold. We also found 31 differentially expressed Long non-coding RNAs (lncRNAs), most of them up-regulated in treated plants. Two of these lncRNAs could operate via microRNAs (miRNAs) target mimic. Our results suggest a reprogramming of expression patterns in CBW0101 that affects a number of genes that is closer to the number reported in winter genotypes. These observations could partially explain its moderate tolerance (low proportion of frost-damaged spikes) when exposed to freezing days in the field.

De novo assembly of wheat root transcriptomes and transcriptional signature of longitudinal differentiation

Hidden underground, root systems constitute an important part of the plant for its development, nourishment and sensing the soil environment around it, but we know very little about its genetic regulation in crop plants like wheat. In the present study, we de novo assembled the root transcriptomes in reference cultivar Chinese Spring from RNA-seq reads generated by the 454-GS-FLX and HiSeq platforms. The FLX reads were assembled into 24,986 transcripts with completeness of 54.84%, and the HiSeq reads were assembled into 91,543 high-confidence protein-coding transcripts, 2,404 low-confidence protein-coding transcripts, and 13,181 non-coding transcripts with the completeness of >90%. Combining the FLX and HiSeq assemblies, we assembled a root transcriptome of 92,335 ORF-containing transcripts. Approximately 7% of the coding transcripts and ~2% non-coding transcripts are not present in the current wheat genome assembly. Functional annotation of both assemblies showed similar gene ontology patterns and that ~7% coding and >5% non-coding transcripts are root-specific. Transcription quantification identified 1,728 differentially expressed transcripts between root tips and maturation zone, and functional annotation of these transcripts captured a transcriptional signature of longitudinal development of wheat root. With the transcriptomic resources developed, this study provided the first view of wheat root transcriptome under different developmental zones and laid a foundation for molecular studies of wheat root development and growth using a reverse genetic approach.

Response of microRNAs to cold treatment in the young spikes of common wheat

Background

MicroRNAs (miRNAs) are a class of small non-coding RNAs that play important roles in biotic and abiotic stresses by regulating their target genes. For common wheat, spring frost damage frequently occurs, especially when low temperature coincides with plants at early floral organ differentiation, which may result in significant yield loss. Up to date, the role of miRNAs in wheat response to frost stress is not well understood.

Results

We report here the sequencing of small RNA transcriptomes from the young spikes that were treated with cold stress and the comparative analysis with those of the control. A total of 192 conserved miRNAs from 105 families and nine novel miRNAs were identified. Among them, 34 conserved and five novel miRNAs were differentially expressed between the cold-stressed samples and the controls. The expression patterns of 18 miRNAs were further validated by quantitative real time polymerase chain reaction (qRT-PCR). Moreover, nearly half of the miRNAs were cross inducible by biotic and abiotic stresses when compared with previously published work. Target genes were predicted and validated by degradome sequencing. Gene Ontology (GO) enrichment analysis showed that the target genes of differentially expressed miRNAs were enriched for response to the stimulus, regulation of transcription, and ion transport functions. Since many targets of differentially expressed miRNAs were transcription factors that are associated with floral development such as ARF, SPB (Squamosa Promoter Binding like protein), MADS-box (MCM1, AG, DEFA and SRF), MYB, SPX (SYG1, Pho81 and XPR1), TCP (TEOSINTE BRANCHED, Cycloidea and PCF), and PPR (PentatricoPeptide Repeat) genes, cold-altered miRNA expression may cause abnormal reproductive organ development.

Conclusion

Analysis of small RNA transcriptomes and their target genes provide new insight into miRNA regulation in developing wheat inflorescences under cold stress. MiRNAs provide another layer of gene regulation in cold stress response that can be genetically manipulated to reduce yield loss in wheat.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3556-2) contains supplementary material, which is available to authorized users.

Identification of Winter-Responsive Proteins in Bread Wheat Using Proteomics Analysis and Virus-Induced Gene Silencing (VIGS)

Proteomic approaches were applied to identify protein spots involved in cold responses in wheat. By comparing the differentially accumulated proteins from two cultivars (UC1110 and PI 610750) and their derivatives, as well as the F10 recombinant inbred line population differing in cold-tolerance, a total of 20 common protein spots representing 16 unique proteins were successfully identified using 2-DE method. Of these, 14 spots had significantly enhanced abundance in the cold-sensitive parental cultivar UC1110 and its 20 descendant lines when compared with the cold-tolerant parental cultivar PI 610750 and its 20 descendant lines. Six protein spots with reduced abundance were also detected. The identified protein spots are involved in stress/defense, carbohydrate metabolism, protein metabolism, nitrogen metabolism, energy metabolism, and photosynthesis. The 20 differentially expressed protein spots were chosen for quantitative real-time polymerase chain reaction (qRT-PCR) to investigate expression changes at the RNA level. The results indicated that the transcriptional expression patterns of 11 genes were consistent with their protein expression models. Among the three unknown proteins, Spot 20 (PAP6-like) showed high sequence similarities with PAP6. qRT-PCR results implied that cold and salt stresses increased the expression of PAP6-like in wheat leaves. Furthermore, VIGS (virus-induced gene silencing)-treated plants generated for PAP6-like were subjected to freezing stress, these plants had more serious droop and wilt, an increased rate of relative electrolyte leakage, reduced relative water content (RWC) and decreased tocopherol levels when compared with viral control plants. However, the plants that were silenced for the other two unknown proteins had no significant differences in comparison to the BSMV0-inoculated plants under freezing conditions. These results indicate that PAP6-like possibly plays an important role in conferring cold tolerance in wheat.

Transcriptome Analysis Provides Insights into the Mechanisms Underlying Wheat Plant Resistance to Stripe Rust at the Adult Plant Stage

Stripe rust (or yellow rust), which is caused by Puccinia striiformis f. sp. tritici (Pst), is one of the most devastating wheat diseases worldwide. The wheat cultivar Xingzi 9104 (XZ) is an elite wheat germplasm that possesses adult plant resistance (APR), which is non–race-specific and durable. Thus, to better understand the mechanism underlying APR, we performed transcriptome sequencing of wheat seedlings and adult plants without Pst infection, and a total of 157,689 unigenes were obtained as a reference. In total, 2,666, 783 and 2,587 differentially expressed genes (DEGs) were found to be up- or down-regulated after Pst infection at 24, 48 and 120 hours post-inoculation (hpi), respectively, based on a comparison of Pst- and mock-infected plants. Among these unigenes, the temporal pattern of the up-regulated unigenes exhibited transient expression patterns during Pst infection, as determined through a Gene Ontology (GO) enrichment analysis. In addition, a Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that many biological processes, including phenylpropanoid biosynthesis, reactive oxygen species, photosynthesis and thiamine metabolism, which mainly control the mechanisms of lignification, reactive oxygen species and sugar, respectively, are involved in APR. In particular, the continuous accumulation of reactive oxygen species may potentially contribute to the ability of the adult plant to inhibit fungal growth and development. To validate the bioinformatics results, 6 candidate genes were selected for further functional identification using the virus-induced gene silencing (VIGS) system, and 4 candidate genes likely contribute to plant resistance against Pst infection. Our study provides new information concerning the transcriptional changes that occur during the Pst-wheat interaction at the adult stage and will help further our understanding of the detailed mechanisms underlying APR to Pst.

Transcriptome-based gene expression profiling identifies differentially expressed genes critical for salt stress response in radish (Raphanus sativus L.)

Transcriptome-based gene expression analysis identifies many critical salt-responsive genes in radish and facilitates further dissecting the molecular mechanism underlying salt stress response. Salt stress severely impacts plant growth and development. Radish, a moderately salt-sensitive vegetable crop, has been studied for decades towards the physiological and biochemical performances under salt stress. However, no systematic study on isolation and identification of genes involved in salt stress response has been performed in radish, and the molecular mechanism governing this process is still indistinct. Here, the RNA-Seq technique was applied to analyze the transcriptomic changes on radish roots treated with salt (200 mM NaCl) for 48 h in comparison with those cultured in normal condition. Totally 8709 differentially expressed genes (DEGs) including 3931 up- and 4778 down-regulated genes were identified. Functional annotation analysis indicated that many genes could be involved in several aspects of salt stress response including stress sensing and signal transduction, osmoregulation, ion homeostasis and ROS scavenging. The association analysis of salt-responsive genes and miRNAs exhibited that 36 miRNA-mRNA pairs had negative correlationship in expression trends. Reverse-transcription quantitative PCR (RT-qPCR) analysis revealed that the expression profiles of DEGs were in line with results from the RNA-Seq analysis. Furthermore, the putative model of DEGs and miRNA-mediated gene regulation was proposed to elucidate how radish sensed and responded to salt stress. This study represents the first comprehensive transcriptome-based gene expression profiling under salt stress in radish. The outcomes of this study could facilitate further dissecting the molecular mechanism underlying salt stress response and provide a valuable platform for further genetic improvement of salt tolerance in radish breeding programs.

Comparative transcriptome analysis to reveal genes involved in wheat hybrid necrosis

Wheat hybrid necrosis is an interesting genetic phenomenon that is found frequently and results in gradual death or loss of productivity of wheat. However, the molecular basis and mechanisms of this genetic phenomenon are still not well understood. In this study, the transcriptomes of wheat hybrid necrosis F₁ and its parents (Neimai 8 and II469) were investigated using digital gene expression (DGE). A total of 1300 differentially expressed genes were identified, indicating that the response to hybrid necrosis in wheat is complicated. The assignments of the annotated genes based on Gene Ontology (GO) revealed that most of the up-regulated genes belong to “universal stress related”, “DNA/RNA binding”, “protein degradation” functional groups, while the down-regulated genes belong to “carbohydrate metabolism” and “translation regulation” functional groups. These findings suggest that these pathways were affected by hybrid necrosis. Our results provide preliminarily new insight into the underlying molecular mechanisms of hybrid necrosis and will help to identify important candidate genes involved in wheat hybrid necrosis.