Harnessing Next Generation Sequencing in Climate Change: RNA-Seq Analysis of Heat Stress-Responsive Genes in Wheat (Triticum aestivum L.)

Physical map of QTL for eleven agronomic traits across fifteen environments, identification of related candidate genes, and development of KASP markers with emphasis on terminal heat stress tolerance in common wheat

KEY MESSAGE

Key message This study identified stable QTL, promising candidate genes and developed novel KASP markers for heat tolerance, providing genomic resources to assist breeding for the development of high-yielding and heat-tolerant wheat germplasm and varieties. To understand the genetic architecture of eleven agronomic traits under heat stress, we used a doubled-haploid population (177 lines) derived from a heat-sensitive cultivar (PBW343) and a heat-tolerant genotype (KSG1203). This population was evaluated under timely, late and very late sown conditions over locations and years comprising fifteen environments. Best linear unbiased estimates and a genetic map (5,710 SNPs) developed using sequencing-based genotyping were used for QTL mapping. The identified 66 QTL (20 novel) were integrated into wheat physical map (14,263.4 Mb). These QTL explained 5.3% (QDth.ccsu-4A for days to heading and QDtm.ccsu-5B for days to maturity) to 24.9% (QGfd.ccsu-7D for grain filling duration) phenotypic variation. Thirteen stable QTL explaining high phenotypic variation were recommended for marker-assisted recurrent selection (MARS) for optimum/heat stress environments. Selected QTL were validated by their presence in high-yielding doubled-haploid lines. Some QTL for 1000-grain weight (TaERF3-3B, TaFER-5B, and TaZIM-A1), grain yield (TaCol-B5), and developmental traits (TaVRT-2) were co-localized with known genes. Specific known genes for traits like abiotic/biotic stress, grain quality and yield were co-located with 26 other QTL. Furthermore, 209 differentially expressed candidate genes for heat tolerance in plants that encode 28 different proteins were identified. KASP markers for three major/stable QTL, namely QGfd.ccsu-7A for grain filling duration on chromosome 7A (timely sown), QNgs.ccsu-3A for number of grains per spike on 3A, and QDth.ccsu-7A for days to heading on 7A (late and very late sown) environments were developed for MARS focusing on the development of heat-tolerant wheat varieties/germplasm.

Characterization of putative calcium-dependent protein kinase-1 (TaCPK-1) gene: hubs in signalling and tolerance network of wheat under terminal heat

Calcium-dependent protein kinase (CDPK) is member of one of the most important signalling cascades operating inside the plant system due to its peculiar role as thermo-sensor. Here, we identified 28 full length putative CDPKs from wheat designated as TaCDPK (1-28). Based on digital gene expression, we cloned full length TaCPK-1 gene of 1691 nucleotides with open reading frame (ORF) of 548 amino acids (accession number OP125853). The expression of TaCPK-1 was observed maximum (3.1-fold) in leaf of wheat cv. HD2985 (thermotolerant) under T2 (38 ± 3 °C, 2 h), as compared to control. A positive correlation was observed between the expression of TaCPK-1 and other stress-associated genes (MAPK6, CDPK4, HSFA6e, HSF3, HSP17, HSP70, SOD and CAT) involved in thermotolerance. Global protein kinase assay showed maximum activity in leaves, as compared to root, stem and spike under heat stress. Immunoblot analysis showed abundance of CDPK protein in wheat cv. HD2985 (thermotolerant) in response to T2 (38 ± 3 °C, 2 h), as compared to HD2329 (thermosusceptible). Calcium ion (Ca2+), being inducer of CDPK, showed strong Ca-signature in the leaf tissue (Ca-622 ppm) of thermotolerant wheat cv. under heat stress, whereas it was minimum (Ca-201 ppm) in spike tissue. We observed significant variations in the ionome of wheat under HS. To conclude, TaCPK-1 plays important role in triggering signaling network and in modulation of HS-tolerance in wheat.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-03989-6.

Discovering the regulators of heat stress tolerance in Ziziphus nummularia (Burm.f) wight and walk.-arn

Ziziphus nummularia an elite heat-stress tolerant shrub, grows in arid regions of desert. However, its molecular mechanism responsible for heat stress tolerance is unexplored. Therefore, we analysed whole transcriptome of Jaisalmer (heat tolerant) and Godhra (heat sensitive) genotypes of Z. nummularia to understand its molecular mechanism responsible for heat stress tolerance. De novo assembly of 16,22,25,052 clean reads yielded 276,029 transcripts. A total of 208,506 unigenes were identified which contains 4290 and 1043 differentially expressed genes (DEG) in TGO (treated Godhra at 42 °C) vs. CGO (control Godhra) and TJR (treated Jaisalmer at 42 °C) vs. CJR (control Jaisalmer), respectively. A total of 987 (67 highly enriched) and 754 (34 highly enriched) pathways were obsorved in CGO vs. TGO and CJR vs. TJR, respectively. Antioxidant pathways and TFs like Homeobox, HBP, ARR, PHD, GRAS, CPP, and E2FA were uniquely observed in Godhra genotype and SET domains were uniquely observed in Jaisalmer genotype. Further transposable elements were highly up-regulated in Godhra genotype but no activation in Jaisalmer genotype. A total of 43,093 and 39,278 simple sequence repeats were identified in the Godhra and Jaisalmer genotypes, respectively. A total of 10 DEGs linked to heat stress were validated in both genotypes for their expression under different heat stresses using quantitative real-time PCR. Comparing expression patterns of the selected DEGs identified ClpB1 as a potential candidate gene for heat tolerance in Z. nummularia. Here we present first characterized transcriptome of Z. nummularia in response to heat stress for the identification and characterization of heat stress-responsive genes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-024-01431-y.

Chemical-tag-based semi-annotated metabolomics facilitates gene identification and specialized metabolic pathway elucidation in wheat

The importance of metabolite modification and species-specific metabolic pathways has long been recognized. However, linking the chemical structure of metabolites to gene function in order to explore the genetic and biochemical basis of metabolism has not yet been reported in wheat (Triticum aestivum). Here, we profiled metabolic fragment enrichment in wheat leaves and consequently applied chemical-tag-based semi-annotated metabolomics in a genome-wide association study in accessions of wheat. The studies revealed that all 1,483 quantified metabolites have at least one known functional group whose modification is tailored in an enzyme-catalyzed manner and eventually allows efficient candidate gene mining. A Triticeae crop-specific flavonoid pathway and its underlying metabolic gene cluster were elucidated in further functional studies. Additionally, upon overexpressing the major effect gene of the cluster TraesCS2B01G460000 (TaOMT24), the pathway was reconstructed in rice (Oryza sativa), which lacks this pathway. The reported workflow represents an efficient and unbiased approach for gene mining using forward genetics in hexaploid wheat. The resultant candidate gene list contains vast molecular resources for decoding the genetic architecture of complex traits and identifying valuable breeding targets and will ultimately aid in achieving wheat crop improvement.

Identifying long non-coding RNAs involved in heat stress response during wheat pollen development

Introduction

Wheat is a staple food crop for over one-third of the global population. However, the stability of wheat productivity is threatened by heat waves associated with climate change. Heat stress at the reproductive stage can result in pollen sterility and failure of grain development.

Methods

This study used transcriptome data analysis to explore the specific expression of long non-coding RNAs (lncRNAs) in response to heat stress during pollen development in four wheat cultivars.

Results and discussion

We identified 11,054 lncRNA-producing loci, of which 5,482 lncRNAs showed differential expression in response to heat stress. Heat-responsive lncRNAs could target protein-coding genes in cis and trans and in lncRNA-miRNA-mRNA regulatory networks. Gene ontology analysis predicted that target protein-coding genes of lncRNAs regulate various biological processes such as hormonal responses, protein modification and folding, response to stress, and biosynthetic and metabolic processes. We also noted some paired lncRNA/protein-coding gene modules and some lncRNA-miRNA-mRNA regulatory modules shared in two or more wheat cultivars. These modules were related to regulating plant responses to heat stress, such as heat-shock proteins and transcription factors, and protein domains, such as MADS-box, Myc-type, and Alpha crystallin/Hsp20 domain.

Conclusion

Our results provide the basic knowledge and molecular resources for future functional studies investigating wheat reproductive development under heat stress.

QTL cluster analysis and marker development for kernel traits based on DArT markers in spring bread wheat (Triticum aestivum L.)

Genetic dissection of yield component traits including kernel characteristics is essential for the continuous improvement in wheat yield. In the present study, one recombinant inbred line (RIL) F₆ population derived from a cross between Avocet and Chilero was used to evaluate the phenotypes of kernel traits of thousand-kernel weight (TKW), kernel length (KL), and kernel width (KW) in four environments at three experimental stations during the 2018-2020 wheat growing seasons. The high-density genetic linkage map was constructed with the diversity arrays technology (DArT) markers and the inclusive composite interval mapping (ICIM) method to identify the quantitative trait loci (QTLs) for TKW, KL, and KW. A total of 48 QTLs for three traits were identified in the RIL population on the 21 chromosomes besides 2A, 4D, and 5B, accounting for 3.00%-33.85% of the phenotypic variances. Based on the physical positions of each QTL, nine stable QTL clusters were identified in the RILs, and among these QTL clusters, TaTKW-1A was tightly linked to the DArT marker interval 3950546-1213099, explaining 10.31%-33.85% of the phenotypic variances. A total of 347 high-confidence genes were identified in a 34.74-Mb physical interval. TraesCS1A02G045300 and TraesCS1A02G058400 were among the putative candidate genes associated with kernel traits, and they were expressed during grain development. Moreover, we also developed high-throughput kompetitive allele-specific PCR (KASP) markers of TaTKW-1A, validated in a natural population of 114 wheat varieties. The study provides a basis for cloning the functional genes underlying the QTL for kernel traits and a practical and accurate marker for molecular breeding.

ZmDWF1 regulates leaf angle in maize

Leaf angle (LA) is a critical agronomic trait enhancing grain yield under high-density planting in maize. A number of researches have been conducted in recent years to investigate the quantitative trait loci/genes responsible for LA variation, while only a few genes were identified through map-based cloning. Here we cloned the ZmDWF1 gene, which was previously reported to encode Δ24-sterol reductase in the brassinosteroids (BRs) biosynthesis pathway. Overexpression of ZmDWF1 resulted in enlarged LA, indicating that ZmDWF1 is a positive regulator of LA in maize. To reveal the regulatory framework of ZmDWF1, we conducted RNA-Sequencing and yeast-two hybrid (Y2H) screening analysis. RNA-Sequencing analyzing results indicate ZmDWF1 mainly affected expression level of genes involved in cell wall associated metabolism and hormone metabolism including BR, gibberellin, and auxin. Y2H screening with Bi-FC assay confirmed three proteins (ZmPP2C-1, ZmROF1, and ZmTWD1) interacting with ZmDWF1. We revealed a new regulatory network of ZmDWF1 gene in controlling plant architecture in maize.

Transcriptome Profiling in Leaves of Wheat Genotype under Heat Stress

Hexaploid wheat is the main cereal food crop for most people but it is highly influenced by climatic variations. The influence of these climatic variations was studies in wheat genotype WH -1184 in field conditions under two different environments (normal and late sown) and it was found that the genotype is less yielding under late sown conditions. To study the effects of heat stress at transcript level, it was grown under two different conditions (WH-1184 control and heat treated) in pots and transcriptome analysis based on Illumina Novoseq 6000 was carried out for the identification of the differentially expressed genes (DEGs) and metabolic processes or gene regulations influenced by heat stress which lead to a reduction in both quality and quantity of wheat production. These DEGs were utilized to set up a subsequent unigene assembly and GO analysis was performed using unigenes to analyze functions of DEGs which were classified into three main domains, i.e., biological process, cellular component, and molecular function. KEGG (Kyoto Encyclopedia of Genes and Genomes) ontology was used to visualize the physiological processes or to identify KEGG pathways that provide plants their ability to shield in adverse conditions of heat stress. From KEGG ontology, it was reported that genes which encoded protein detoxification and ABC1 domain-containing protein were upregulated while genes thatencoded glutathione transferase (GST), peroxidase, and chitinase enzymes were downregulated. Downregulation of these enzymes during heat stress causes oxidative damages in plants while upregulated proteins play a main role in detoxification to protect plants from heat stress. It was hypothesized that the yield of WH-1184 decreased 44% under heat stress due to the downregulation of genes that encoded GST, peroxidase, and chitinase enzymes which can protect plants from oxidative damage. Hence, upregulation of these genes might be helpful for the adaptation of this genotype under heat stress condition.

Identification of long non-coding RNA-microRNA-mRNA regulatory modules and their potential roles in drought stress response in wheat (Triticum aestivum L.)

Drought is one of the most severe abiotic stresses that influence wheat production across the globe. Understanding the molecular regulatory network of wheat in response to drought is of great importance in molecular breeding. Noncoding RNAs influence plant development and resistance to abiotic stresses by regulating gene expression. In this study, whole-transcriptome sequencing was performed on the seedlings of two wheat varieties with contrasting levels of drought tolerance under drought and control conditions to identify long noncoding RNAs (lncRNAs), micro RNAs (miRNAs), and mRNAs related to drought stress and explore the potential lncRNA-miRNA-mRNA regulatory modules in controlling wheat drought stress response. A total of 1515 differentially expressed lncRNAs (DELs), 209 differentially expressed miRNAs (DEMs), and 20462 differentially expressed genes (DEGs) were identified. Of the 20462 DEGs, 1025 were identified as potential wheat drought resistance-related DEGs. Based on the regulatory relationship and expression patterns of DELs, DEMs, and DEGs, 10 DEL-DEM-DEG regulatory modules related to wheat drought stress response were screened, and preliminary expression verification of two important candidate modules was performed. Our results revealed the possible roles of lncRNA-miRNA-mRNA modules in regulatory networks related to drought tolerance and provided useful information as valuable genomic resources in molecular breeding of wheat.

Comparison of Gene Expression Changes in Three Wheat Varieties with Different Susceptibilities to Heat Stress Using RNA-Seq Analysis

Wheat is highly susceptible to heat stress, which significantly reduces grain yield. In this study, we used RNA-seq technology to analyze the transcript expression at three different time-points after heat treatment in three cultivars differing in their susceptibility to heat stress: Jopum, Keumkang, and Olgeuru. A total of 11,751, 8850, and 14,711; 10,959, 7946, and 14,205; and 22,895, 13,060, and 19,408 differentially-expressed genes (log2 fold-change > 1 and FDR (padj) < 0.05) were identified in Jopum, Keumkang, and Olgeuru in the control vs. 6-h, in the control vs. 12-h, and in the 6-h vs. 12-h heat treatment, respectively. Functional enrichment analysis showed that the biological processes for DEGs, such as the cellular response to heat and oxidative stress—and including the removal of superoxide radicals and the positive regulation of superoxide dismutase activity—were significantly enriched among the three comparisons in all three cultivars. Furthermore, we investigated the differential expression patterns of reactive oxygen species (ROS)-scavenging enzymes, heat shock proteins, and heat-stress transcription factors using qRT-PCR to confirm the differences in gene expression among the three varieties under heat stress. This study contributes to a better understanding of the wheat heat-stress response at the early growth stage and the varietal differences in heat tolerance.

Transcriptome profiling reveals the genes and pathways involved in thermo-tolerance in wheat (Triticum aestivum L.) genotype Raj 3765

Wheat, one of the most widely consumed staple food crops globally, is relatively vulnerable to high temperature-induced heat stress. It is therefore essential to gain more insight into the comprehensive mechanism of thermotolerance of wheat in order to safeguard its production. In view of this, we analysed heat stress responsive transcriptome data of wheat to determine its gene expression level under heat stress. A total of 7990 DEGs, including 4483 up-regulated and 3507 down regulated genes were identified. Gene Ontology (GO) analysis categorized 3910 DEGs into different ontology families. 146 pathways involving 814 DEGs were enriched during KEGG analysis. Metabolic pathways and biosynthesis of secondary metabolites were the major pathways enriched. MYB (myeloblastosis) transcription factors (TFs) and many other TFs as bHLH, WRKY, NAC, ERF, were determined to be quite abundant in the DEGs. Since various reports indicate that these TFs play important role in plants abiotic stress, it is an indication that our DEGs are functional in heat stress tolerance. Verification of few selected DEGs using RT-qPCR produced expression levels similar to the transcriptome data. This indicates that the transcriptome data is reliable. These results could be helpful in enhancing our understanding of the mechanism underlying thermotolerance in wheat.

Identification of Novel Genes Associated with Partial Resistance to Aphanomyces Root Rot in Field Pea by BSR-Seq Analysis

Aphanomyces root rot, caused by Aphanomyces euteiches, causes severe yield loss in field pea (Pisum sativum). The identification of a pea germplasm resistant to this disease is an important breeding objective. Polygenetic resistance has been reported in the field pea cultivar '00-2067'. To facilitate marker-assisted selection (MAS), bulked segregant RNA-seq (BSR-seq) analysis was conducted using an F8 RIL population derived from the cross of 'Carman' × '00-2067'. Root rot development was assessed under controlled conditions in replicated experiments. Resistant (R) and susceptible (S) bulks were constructed based on the root rot severity in a greenhouse study. The BSR-seq analysis of the R bulks generated 44,595,510~51,658,688 reads, of which the aligned sequences were linked to 44,757 genes in a reference genome. In total, 2356 differentially expressed genes were identified, of which 44 were used for gene annotation, including defense-related pathways (jasmonate, ethylene and salicylate) and the GO biological process. A total of 344.1 K SNPs were identified between the R and S bulks, of which 395 variants were located in 31 candidate genes. The identification of novel genes associated with partial resistance to Aphanomyces root rot in field pea by BSR-seq may facilitate efforts to improve management of this important disease.

Transcriptomic and physiological responses of contrasting maize genotypes to drought stress

Drought is a significant environmental stress factor that adversely affects maize productivity. However, many details regarding the molecular mechanisms of maize against drought are still unclear. In this study, leaf transcriptomics and physiological traits of two maize genotypes with differing drought resistance were analyzed. Transcriptome sequencing identified 8985 and 7305 differentially expressed genes (DEGs) in SD902 and SD609, respectively. Functional analysis suggested that numerous genes are highly involved in oxidative defense, protein modification, photosynthesis, phytohormone response, MAPK signaling, and transcription factors (TFs). Compared to SD902, SD609 had a higher expression of DEGs related to antioxidant enzymes, photosynthetic electron transport, heat shock proteins, and indole-3-acetic acid (IAA) signaling under drought conditions, which might contribute to its tolerance mechanisms to drought. Stress-induced TFs may play a crucial regulatory role in genotypic differences. Moreover, the physiological changes and gene expression abundance determined using quantitative reverse transcription polymerase chain reaction were consistent with the RNA sequencing data. The study results suggest that the higher drought tolerance of SD609 than SD902 can be attributed to stronger stress defense capabilities, IAA signal transduction, and more stable photosynthesis. Our findings provide new insights into the molecular mechanisms of maize against drought stress, and the candidate genes identified may be used in breeding drought-tolerant maize cultivars.

RNA-Seq Analysis of Developing Grains of Wheat to Intrigue Into the Complex Molecular Mechanism of the Heat Stress Response

Heat stress is one of the significant constraints affecting wheat production worldwide. To ensure food security for ever-increasing world population, improving wheat for heat stress tolerance is needed in the presently drifting climatic conditions. At the molecular level, heat stress tolerance in wheat is governed by a complex interplay of various heat stress-associated genes. We used a comparative transcriptome sequencing approach to study the effect of heat stress (5°C above ambient threshold temperature of 20°C) during grain filling stages in wheat genotype K7903 (Halna). At 7 DPA (days post-anthesis), heat stress treatment was given at four stages: 0, 24, 48, and 120 h. In total, 115,656 wheat genes were identified, including 309 differentially expressed genes (DEGs) involved in many critical processes, such as signal transduction, starch synthetic pathway, antioxidant pathway, and heat stress-responsive conserved and uncharacterized putative genes that play an essential role in maintaining the grain filling rate at the high temperature. A total of 98,412 Simple Sequences Repeats (SSR) were identified from de novo transcriptome assembly of wheat and validated. The miRNA target prediction from differential expressed genes was performed by psRNATarget server against 119 mature miRNA. Further, 107,107 variants including 80,936 Single nucleotide polymorphism (SNPs) and 26,171 insertion/deletion (Indels) were also identified in de novo transcriptome assembly of wheat and wheat genome Ensembl version 31. The present study enriches our understanding of known heat response mechanisms during the grain filling stage supported by discovery of novel transcripts, microsatellite markers, putative miRNA targets, and genetic variant. This enhances gene functions and regulators, paving the way for improved heat tolerance in wheat varieties, making them more suitable for production in the current climate change scenario.

The genetic and molecular basis for improving heat stress tolerance in wheat

Wheat production requires at least ~ 2.4% increase per year rate by 2050 globally to meet food demands. However, heat stress results in serious yield loss of wheat worldwide. Correspondingly, wheat has evolved genetic basis and molecular mechanisms to protect themselves from heat-induced damage. Thus, it is very urgent to understand the underlying genetic basis and molecular mechanisms responsive to elevated temperatures to provide important strategies for heat-tolerant varieties breeding. In this review, we focused on the impact of heat stress on morphology variation at adult stage in wheat breeding programs. We also summarize the recent studies of genetic and molecular factors regulating heat tolerance, including identification of heat stress tolerance related QTLs/genes, and the regulation pathway in response to heat stress. In addition, we discuss the potential ways to improve heat tolerance by developing new technologies such as genome editing. This review of wheat responses to heat stress may shed light on the understanding heat-responsive mechanisms, although the regulatory network of heat tolerance is still ambiguous in wheat.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42994-021-00064-z.

Identification of Potential Cytokinin Responsive Key Genes in Rice Treated With Trans-Zeatin Through Systems Biology Approach

Rice is an important staple food grain consumed by most of the population around the world. With climate and environmental changes, rice has undergone a tremendous stress state which has impacted crop production and productivity. Plant growth hormones are essential component that controls the overall outcome of the growth and development of the plant. Cytokinin is a hormone that plays an important role in plant immunity and defense systems. Trans-zeatin is an active form of cytokinin that can affect plant growth which is mediated by a multi-step two-component phosphorelay system that has different roles in various developmental stages. Systems biology is an approach for pathway analysis to trans-zeatin treated rice that could provide a deep understanding of different molecules associated with them. In this study, we have used a weighted gene co-expression network analysis method to identify the functional modules and hub genes involved in the cytokinin pathway. We have identified nine functional modules comprising of different hub genes which contribute to the cytokinin signaling route. The biological significance of these identified hub genes has been tested by applying well-proven statistical techniques to establish the association with the experimentally validated QTLs and annotated by the DAVID server. The establishment of key genes in different pathways has been confirmed. These results will be useful to design new stress-resistant cultivars which can provide sustainable yield in stress-specific conditions.

Transcriptome based identification and validation of heat stress transcription factors in wheat progenitor species Aegilops speltoides

Wheat, one of the major cereal crops worldwide, get adversely affected by rising global temperature. We have identified the diploid B genome progenitor of wheat, Aegilops speltoides (SS), as a potential donor for heat stress tolerance. Therefore, the present work was planned to study the total transcriptome profile of heat stress-tolerant Ae. speltoides accession pau3809 (AS3809) and compare with that of tetraploid and hexaploid wheat cultivars PDW274 and PBW725, respectively. The comparative transcriptome was utilized to identify and validate heat stress transcription factors (HSFs), the key genes involved in imparting heat stress tolerance. Transcriptome analysis led to the identification of a total of 74 K, 68 K, and 76 K genes in AS3809, PDW274, and PBW725, respectively. There was a high uniformity of GO profiles under the biological, molecular, and cellular functions across the three wheat transcriptomes, suggesting the conservation of gene function. Twelve HSFs having the highest FPKM value were identified in the AS3809 transcriptome data, while six of these HSFs namely HSFA3, HSFA5, HSFA9, HSFB2a, HSFB2b, and HSFC1b, were validated with qRT PCR. These six HSFs were identified as an important component of thermotolerance in AS3809 as evident from their comparative higher expression under heat stress.

Reproductive developmental transcriptome analysis of Tripidium ravennae (Poaceae)

BACKGROUND

Tripidium ravennae is a cold-hardy, diploid species in the sugarcane complex (Poaceae subtribe Saccharinae) with considerable potential as a genetic resource for developing improved bioenergy and ornamental grasses. An improved understanding of the genetic regulation of reproductive processes (e.g., floral induction, inflorescence development, and seed development) will enable future applications of precision breeding and gene editing of floral and seed development. In particular, the ability to silence reproductive processes would allow for developing seedless forms of valuable but potentially invasive plants. The objective of this research was to characterize the gene expression environment of reproductive development in T. ravennae.

RESULTS

During the early phases of inflorescence development, multiple key canonical floral integrators and pathways were identified. Annotations of type II subfamily of MADS-box transcription factors, in particular, were over-represented in the GO enrichment analyses and tests for differential expression (FDR p-value < 0.05). The differential expression of floral integrators observed in the early phases of inflorescence development diminished prior to inflorescence determinacy regulation. Differential expression analysis did not identify many unique genes at mid-inflorescence development stages, though typical biological processes involved in plant growth and development expressed abundantly. The increase in inflorescence determinacy regulatory elements and putative homeotic floral development unigenes at mid-inflorescence development coincided with the expression of multiple meiosis annotations and multicellular organism developmental processes. Analysis of seed development identified multiple unigenes involved in oxidative-reductive processes.

CONCLUSION

Reproduction in grasses is a dynamic system involving the sequential coordination of complex gene regulatory networks and developmental processes. This research identified differentially expressed transcripts associated with floral induction, inflorescence development, and seed development in T. ravennae. These results provide insights into the molecular regulation of reproductive development and provide a foundation for future investigations and analyses, including genome annotation, functional genomics characterization, gene family evolutionary studies, comparative genomics, and precision breeding.

Transcriptome analysis of Sonneratia caseolaris seedlings under chilling stress

Sonneratia caseolaris is a native mangrove species found in China. It is fast growing and highly adaptable for mangrove afforestation, but suffered great damage by chilling event once introduced to high latitude area. To understand the response mechanisms under chilling stress, physiological and transcriptomic analyses were conducted. The relative electrolyte conductivity, malondialdehyde (MDA) content, soluble sugar content and soluble protein content increased significantly under chilling stress. This indicated that S. caseolaris suffered great damage and increased the levels of osmoprotectants in response to the chilling stress. Gene expression comparison analysis of S. caseolaris leaves after 6 h of chilling stress was performed at the transcriptional scale using RNA-Seq. A total of 168,473 unigenes and 3,706 differentially expressed genes (DEGs) were identified. GO and KEGG enrichment analyses showed that the DEGs were mainly involved in carbohydrate metabolism, antioxidant enzyme, plant hormone signal transduction, and transcription factors (TFs). Sixteen genes associated with carbohydrate metabolism, antioxidant enzyme, phytohormones and TFs were selected for qRT-PCR verification, and they indicated that the transcriptome data were reliable. Our work provided a comprehensive review of the chilling response of S. caseolaris at both physiological and transcriptomic levels, which will prove useful for further studies on stress-responses in mangrove plants.

Comparative transcriptome analyses revealed different heat stress responses in pigeonpea (Cajanus cajan) and its crop wild relatives

Comparative transcriptome analyses accompanied by biochemical assays revealed high variability in heat stress response in Cajanus species. Among the studied species, C. scarabaeoides was the most thermotolerant followed by C. cajanifolius, C. cajan, and C. acutifolius. Pigeonpea is one of the climate-resilient grain legumes. Though the optimum temperature for cultivated pigeonpea is ~ 25-35 °C, its wild relatives grow in temperatures ranging between 18 and 45 °C. To gain insight into molecular mechanisms responsible for the heat stress tolerance in pigeonpea, we conducted time-series transcriptome analysis of one pigeonpea cultivar (Cajanus cajan) and two wild relatives, Cajanus acutifolius, and Cajanus scarabaeoides subjected to heat stress at 42 ± 2 ºC for 30 min and 3 h. A total of 9521, 12,447, and 5282 identified transcripts were differentially expressed in C. cajan, C. acutifolius, and C. scarabaeoides, respectively. In this study, we observed that a significant number of genes undergo alternative splicing in a species-specific pattern during heat stress. Gene expression profiling analysis, histochemical assay, chlorophyll content, and electrolyte leakage assay showed that C. scarabaeoides has adaptive features for heat stress tolerance. The gene set enrichment analyses of differentially expressed genes in these Cajanus species during heat stress revealed that oxidoreductase activity, transcription factor activity, oxygen-evolving complex, photosystem-II, thylakoid, phenylpropanoid biosynthetic process, secondary metabolic process, and flavonoid biosynthetic process were highly affected. The histochemical assay showed more lipid peroxidation in C. acutifolius compared to other Cajanus species inferring the presence of higher quantities of polyunsaturated fatty acids in the plasma membrane which might have led to severe damage of membrane-bound organelles like chloroplast, and high electrolyte leakage during heat stress. This study paves the way for the identification of candidate genes, which can be useful for the development of thermo-tolerant pigeonpea cultivars.

Autophagy in Plant Abiotic Stress Management

Plants can be considered an open system. Throughout their life cycle, plants need to exchange material, energy and information with the outside world. To improve their survival and complete their life cycle, plants have developed sophisticated mechanisms to maintain cellular homeostasis during development and in response to environmental changes. Autophagy is an evolutionarily conserved self-degradative process that occurs ubiquitously in all eukaryotic cells and plays many physiological roles in maintaining cellular homeostasis. In recent years, an increasing number of studies have shown that autophagy can be induced not only by starvation but also as a cellular response to various abiotic stresses, including oxidative, salt, drought, cold and heat stresses. This review focuses mainly on the role of autophagy in plant abiotic stress management.

Role of ATP-binding cassette transporters in maintaining plant homeostasis under abiotic and biotic stresses

The ATP-binding cassette (ABC) transporters belong to a large protein family predominantly present in diverse species. ABC transporters are driven by ATP hydrolysis and can act as exporters as well as importers. These proteins are localized in the membranes of chloroplasts, mitochondria, peroxisomes and vacuoles. ABC proteins are involved in regulating diverse biological processes in plants, such as growth, development, uptake of nutrients, tolerance to biotic and abiotic stresses, tolerance to metal toxicity, stomatal closure, shape and size of grains, protection of pollens, transport of phytohormones, etc. In mitochondria and chloroplast, the iron metabolism and its transport across the membrane are mediated by ABC transporters. Tonoplast-localized ABC transporters are involved in internal detoxification of metal ion; thus protecting against the DNA impairment and maintaining cell growth. ABC transporters are involved in the transport of secondary metabolites inside the cells. Microorganisms also engage a large number of ABC transporters to import and expel substrates decisive for their pathogenesis. ABC transporters also suppress the seed embryonic growth until favorable conditions come. This review aims at giving insights on ABC transporters, their evolution, structure, functions and roles in different biological processes for helping the terrestrial plants to survive under adverse environmental conditions. These specialized plant membrane transporters ensure a sustainable economic yield and high-quality products, especially under unfavorable conditions of growth. These transporters can be suitably manipulated to develop 'Plants for the Future'.

Roles of Brassinosteroids in Mitigating Heat Stress Damage in Cereal Crops

Heat stress causes huge losses in the yield of cereal crops. Temperature influences the rate of plant metabolic and developmental processes that ultimately determine the production of grains, with high temperatures causing a reduction in grain yield and quality. To ensure continued food security, the tolerance of high temperature is rapidly becoming necessary. Brassinosteroids (BR) are a class of plant hormones that impact tolerance to various biotic and abiotic stresses and regulate cereal growth and fertility. Fine-tuning the action of BR has the potential to increase cereals' tolerance and acclimation to heat stress and maintain yields. Mechanistically, exogenous applications of BR protect yields through amplifying responses to heat stress and rescuing the expression of growth promoters. Varied BR compounds and differential signaling mechanisms across cereals point to a diversity of mechanisms that can be leveraged to mitigate heat stress. Further, hormone transport and BR interaction with other molecules in plants may be critical to utilizing BR as protective agrochemicals against heat stress. Understanding the interplay between heat stress responses, growth processes and hormone signaling may lead us to a comprehensive dogma of how to tune BR application for optimizing cereal growth under challenging environments in the field.

RNA-seq analysis reveals different drought tolerance mechanisms in two broadly adapted wheat cultivars 'TAM 111' and 'TAM 112'

Wheat cultivars 'TAM 111' and 'TAM 112' have been dominantly grown in the Southern U.S. Great Plains for many years due to their high yield and drought tolerance. To identify the molecular basis and genetic control of drought tolerance in these two landmark cultivars, RNA-seq analysis was conducted to compare gene expression difference in flag leaves under fully irrigated (wet) and water deficient (dry) conditions. A total of 2254 genes showed significantly altered expression patterns under dry and wet conditions in the two cultivars. TAM 111 had 593 and 1532 dry-wet differentially expressed genes (DEGs), and TAM 112 had 777 and 1670 at heading and grain-filling stages, respectively. The two cultivars have 1214 (53.9%) dry-wet DEGs in common, which agreed with their excellent adaption to drought, but 438 and 602 dry-wet DEGs were respectively shown only in TAM 111 and TAM 112 suggested that each has a specific mechanism to cope with drought. Annotations of all 2254 genes showed 1855 have functions related to biosynthesis, stress responses, defense responses, transcription factors and cellular components related to ion or protein transportation and signal transduction. Comparing hierarchical structure of biological processes, molecule functions and cellular components revealed the significant regulation differences between TAM 111 and TAM 112, particularly for genes of phosphorylation and adenyl ribonucleotide binding, and proteins located in nucleus and plasma membrane. TAM 112 showed more active than TAM 111 in response to drought and carried more specific genes with most of them were up-regulated in responses to stresses of water deprivation, heat and oxidative, ABA-induced signal pathway and transcription regulation. In addition, 258 genes encoding predicted uncharacterized proteins and 141 unannotated genes with no similar sequences identified in the databases may represent novel genes related to drought response in TAM 111 or TAM 112. This research thus revealed different drought-tolerance mechanisms in TAM 111 and TAM 112 and identified useful drought tolerance genes for wheat adaption. Data of gene sequence and expression regulation from this study also provided useful information of annotating novel genes associated with drought tolerance in the wheat genome.

De Novo Transcriptome Analysis of Durum Wheat Flag Leaves Provides New Insights Into the Regulatory Response to Elevated CO2 and High Temperature

Global warming is becoming a significant problem for food security, particularly in the Mediterranean basin. The use of molecular techniques to study gene-level responses to environmental changes in non-model organisms is increasing and may help to improve the mechanistic understanding of durum wheat response to elevated CO2 and high temperature. With this purpose, we performed transcriptome RNA sequencing (RNA-Seq) analyses combined with physiological and biochemical studies in the flag leaf of plants grown in field chambers at ear emergence. Enhanced photosynthesis by elevated CO2 was accompanied by an increase in biomass and starch and fructan content, and a decrease in N compounds, as chlorophyll, soluble proteins, and Rubisco content, in association with a decline of nitrate reductase and initial and total Rubisco activities. While high temperature led to a decline of chlorophyll, Rubisco activity, and protein content, the glucose content increased and starch decreased. Furthermore, elevated CO2 induced several genes involved in mitochondrial electron transport, a few genes for photosynthesis and fructan synthesis, and most of the genes involved in secondary metabolism and gibberellin and jasmonate metabolism, whereas those related to light harvesting, N assimilation, and other hormone pathways were repressed. High temperature repressed genes for C, energy, N, lipid, secondary, and hormone metabolisms. Under the combined increases in atmospheric CO2 and temperature, the transcript profile resembled that previously reported for high temperature, although elevated CO2 partly alleviated the downregulation of primary and secondary metabolism genes. The results suggest that there was a reprogramming of primary and secondary metabolism under the future climatic scenario, leading to coordinated regulation of C-N metabolism towards C-rich metabolites at elevated CO2 and a shift away from C-rich secondary metabolites at high temperature. Several candidate genes differentially expressed were identified, including protein kinases, receptor kinases, and transcription factors.

Identification of putative key genes for coastal environments and cold adaptation in mangrove Kandelia obovata through transcriptome analysis

Mangrove forests are an important contributor to the coastal marine environment. They have developed unique adaptations to the harsh coastal wetland, yet their geographic distribution is limited by environmental temperature. The adaptive strategies of mangrove at the molecular level, however, have not been addressed. In the present work, transcriptome analyses were performed on different cold damaged plants of a mangrove species, Kandelia obovata. From the samples collected in the field after a cold stress, we found that distinct expression profiles of many key genes are related to extreme temperature responses. These include transcription factors such as WRKY and bHLH, and other genes encoding proteins like SnRK2, PR-1, KCS, involving in the pathways of plant hormones, plant-pathogen interactions, and long chain fatty acid synthesis. We also examined the transcriptomes of eight tissues of K. obovata to identify candidate genes involved in adaptation and development. While stress-responsive genes were globally expressed, tissue-specific genes with diverse functions might be involved in tissue development and adaptability. For examples, genes encoding CYP724B1 and ABCB1 were specifically expressed in the fruit and root, respectively. Additionally, 26 genes were identified as positively selected genes in K. obovata, six of them were found to be involved in chilling stress response, seed germination and oxidation-reduction processes, suggesting their roles in stressful environment adaptation. Together, these results shed light into the K. obovata's natural responses to cold snaps at the molecular level, and reveal a global gene expression portrait across different tissues. It also provides a transcriptome resource for further molecular ecology studies and conservation planning of this and other mangrove plants in their native and adopted environments.

iTRAQ-based quantitative proteome analysis reveals metabolic changes between a cleistogamous wheat mutant and its wild-type wheat counterpart

Background

Wheat is one of the most important staple crops worldwide. Fusarium head blight (FHB) severely affects wheat yield and quality. A novel bread wheat mutant, ZK001, characterized as cleistogamic was isolated from a non-cleistogamous variety Yumai 18 (YM18) through static magnetic field mutagenesis. Cleistogamy is a promising strategy for controlling FHB. However, little is known about the mechanism of cleistogamy in wheat.

Methods

We performed a FHB resistance test to identify the FHB infection rate of ZK001. We also measured the agronomic traits of ZK001 and the starch and total soluble sugar contents of lodicules in YM18 and ZK001. Finally, we performed comparative studies at the proteome level between YM18 and ZK001 based on the proteomic technique of isobaric tags for relative and absolute quantification.

Results

The infection rate of ZK001 was lower than that of its wild-type and Aikang 58. The abnormal lodicules of ZK001 lost the ability to push the lemma and palea apart during the flowering stage. Proteome analysis showed that the main differentially abundant proteins (DAPs) were related to carbohydrate metabolism, protein transport, and calcium ion binding. These DAPs may work together to regulate cellular homeostasis, osmotic pressure and the development of lodicules. This hypothesis is supported by the analysis of starch, soluble sugar content in the lodicules as well as the results of Quantitative reverse transcription polymerase chain reaction.

Conclusions

Proteomic analysis has provided comprehensive information that should be useful for further research on the lodicule development mechanism in wheat. The ZK001 mutant is optimal for studying flower development in wheat and could be very important for FHB resistant projects via conventional crossing.

Contribution of time of day and the circadian clock to the heat stress responsive transcriptome in Arabidopsis

In Arabidopsis, a large subset of heat responsive genes exhibits diurnal or circadian oscillations. However, to what extent the dimension of time and/or the circadian clock contribute to heat stress responses remains largely unknown. To determine the direct contribution of time of day and/or the clock to differential heat stress responses, we probed wild-type and mutants of the circadian clock genes CCA1, LHY, PRR7, and PRR9 following exposure to heat (37 °C) and moderate cold (10 °C) in the early morning (ZT1) and afternoon (ZT6). Thousands of genes were differentially expressed in response to temperature, time of day, and/or the clock mutation. Approximately 30% more genes were differentially expressed in the afternoon compared to the morning, and heat stress significantly perturbed the transcriptome. Of the DEGs (~3000) specifically responsive to heat stress, ~70% showed time of day (ZT1 or ZT6) occurrence of the transcriptional response. For the DEGs (~1400) that are shared between ZT1 and ZT6, we observed changes to the magnitude of the transcriptional response. In addition, ~2% of all DEGs showed differential responses to temperature stress in the clock mutants. The findings in this study highlight a significant role for time of day in the heat stress responsive transcriptome, and the clock through CCA1 and LHY, appears to have a more profound role than PRR7 and PRR9 in modulating heat stress responses during the day. Our results emphasize the importance of considering the dimension of time in studies on abiotic stress responses in Arabidopsis.

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.

Wheat seed transcriptome reveals genes controlling key traits for human preference and crop adaptation

Analysis of the transcriptome of the developing wheat grain has associated expression of genes with traits involving production (e.g. yield) and quality (e.g. bread quality). Photosynthesis in the grain may be important in retaining carbon that would be lost in respiration during grain filling and may contribute to yield in the late stages of seed formation under warm and dry environments. A small number of genes have been identified as having been selected by humans to optimize the performance of wheat for foods such as bread. Genes determining flour yield in milling have been discovered. Hardness is explained by variations in expression of pin genes. Knowledge of these genes should dramatically improve the efficiency of breeding better climate adapted wheat genotypes.

Transcriptional profiling and genes involved in acquired thermotolerance in Banana: a non-model crop

Banana is a non- model crop plant, and one of the most important crops in the tropics and sub tropics. Heat stress is the major abiotic stress affecting banana crop production because of its long growth period and is likely to become a threat due to global warming. To understand an acquired thermotolerance phenomenon at the molecular level, the RNA-seq approach was employed by adapting TIR method. A total of 136.38 million high quality reads were assembled. Differentially expressed genes under induction (I) was 3936, I + L was 2268 and lethal stress was 907 compared to control. Gene ontology and DGE analysis showed that genes related to heat shock factors, heat shock proteins, stress associated proteins, ROS scavenging, fatty acid metabolism, protein modification were significantly up regulated during induction, thus preparing the organism or tissue at molecular and cellular level for acquired thermotolerance. KEGG pathway analysis revealed the significant enrichment of pathways involved in protein processing, MAPK signaling and HSPs which indicates that these processes are conserved and involved in thermo tolerance. Thus, this study provides insights into the acquired thermotolerance phenomena in plants especially banana.

Next-Generation Sequencing Reveals Novel Mutations in X-linked Intellectual Disability

Robust diagnostics for many human genetic disorders are much needed in the pursuit of global personalized medicine. Next-generation sequencing now offers new promise for biomarker and diagnostic discovery, in developed as well as resource-limited countries. In this broader global health context, X-linked intellectual disability (XLID) is an inherited genetic disorder that is associated with a range of phenotypes impacting societies in both developed and developing countries. Although intellectual disability arises due to diverse causes, a substantial proportion is caused by genomic alterations. Studies have identified causal XLID genomic alterations in more than 100 protein-coding genes located on the X-chromosome. However, the causes for a substantial number of intellectual disability and associated phenotypes still remain unknown. Identification of causative genes and novel mutations will help in early diagnosis as well as genetic counseling of families. Advent of next-generation sequencing methods has accelerated the discovery of new genes involved in mental health disorders. In this study, we analyzed the exomes of three families from India with nonsyndromic XLID comprising seven affected individuals. The affected individuals had varying degrees of intellectual disability, microcephaly, and delayed motor and language milestones. We identified potential causal variants in three XLID genes, including PAK3 (V294M), CASK (complex structural variant), and MECP2 (P354T). Our findings reported in this study extend the spectrum of mutations and phenotypes associated with XLID, and calls for further studies of intellectual disability and mental health disorders with use of next-generation sequencing technologies.

Exploring the heat-responsive chaperones and microsatellite markers associated with terminal heat stress tolerance in developing wheat

Global warming is a major threat for agriculture and food security, and in many cases the negative impacts are already apparent. Wheat is one of the most important staple food crops and is highly sensitive to the heat stress (HS) during reproductive and grain-filling stages. Here, whole transcriptome analysis of thermotolerant wheat cv. HD2985 was carried out at the post-anthesis stage under control (22 ± 3 °C) and HS-treated (42 °C, 2 h) conditions using Illumina Hiseq and Roche GS-FLX 454 platforms. We assembled ~24 million (control) and ~23 million (HS-treated) high-quality trimmed reads using different assemblers with optimal parameters. De novo assembly yielded 52,567 (control) and 59,658 (HS-treated) unigenes. We observed 785 transcripts to be upregulated and 431 transcripts to be downregulated under HS; 78 transcripts showed >10-fold upregulation such as HSPs, metabolic pathway-related genes, etc. Maximum number of upregulated genes was observed to be associated with processes such as HS-response, protein-folding, oxidation-reduction and photosynthesis. We identified 2008 and 2483 simple sequence repeats (SSRs) markers from control and HS-treated samples; 243 SSRs were observed to be overlying on stress-associated genes. Polymorphic study validated four SSRs to be heat-responsive in nature. Expression analysis of identified differentially expressed transcripts (DETs) showed very high fold increase in the expression of catalytic chaperones (HSP26, HSP17, and Rca) in contrasting wheat cvs. HD2985 and HD2329 under HS. We observed positive correlation between RNA-seq and qRT-PCR expression data. The present study culminated in greater understanding of the heat-response of tolerant genotype and has provided good candidate genes for the marker development and screening of wheat germplasm for thermotolerance.

Identification of Putative RuBisCo Activase (TaRca1)-The Catalytic Chaperone Regulating Carbon Assimilatory Pathway in Wheat (Triticum aestivum) under the Heat Stress

RuBisCo activase (Rca) is a catalytic chaperone involved in modulating the activity of RuBisCo (key enzyme of photosynthetic pathway). Here, we identified eight novel transcripts from wheat through data mining predicted to be Rca and cloned a transcript of 1.4 kb from cv. HD2985, named as TaRca1 (GenBank acc. no. KC776912). Single copy number of TaRca1 was observed in wheat genome. Expression analysis in diverse wheat genotypes (HD2985, Halna, PBW621, and HD2329) showed very high relative expression of TaRca1 in Halna under control and HS-treated, as compared to other cultivars at different stages of growth. TaRca1 protein was predicted to be chloroplast-localized with numerous potential phosphorylation sites. Northern blot analysis showed maximum accumulation of TaRca1 transcript in thermotolerant cv. during mealy-ripe stage, as compared to thermosusceptible. Decrease in the photosynthetic parameters was observed in all the cultivars, except PBW621 in response to HS. We observed significant increase in the Rca activity in all the cultivars under HS at different stages of growth. HS causes decrease in the RuBisCo activity; maximum reduction was observed during pollination stage in thermosusceptible cvs. as validated through immunoblotting. We observed uniform carbon distribution in different tissues of thermotolerant cvs., as compared to thermosusceptible. Similarly, tolerance level of leaf was observed maximum in Halna having high Rca activity under HS. A positive correlation was observed between the transcript and activity of TaRca1 in HS-treated Halna. Similarly, TaRca1 enzyme showed positive correlation with the activity of RuBisCo. There is, however, need to manipulate the thermal stability of TaRca1 enzyme through protein engineering for sustaining the photosynthetic rate under HS-a novel approach toward development of "climate-smart" crop.

Genetic Approaches to Study Plant Responses to Environmental Stresses: An Overview

The assessment of gene expression levels is an important step toward elucidating gene functions temporally and spatially. Decades ago, typical studies were focusing on a few genes individually, whereas now researchers are able to examine whole genomes at once. The upgrade of throughput levels aided the introduction of systems biology approaches whereby cell functional networks can be scrutinized in their entireties to unravel potential functional interacting components. The birth of systems biology goes hand-in-hand with huge technological advancements and enables a fairly rapid detection of all transcripts in studied biological samples. Even so, earlier technologies that were restricted to probing single genes or a subset of genes still have their place in research laboratories. The objective here is to highlight key approaches used in gene expression analysis in plant responses to environmental stresses, or, more generally, any other condition of interest. Northern blots, RNase protection assays, and qPCR are described for their targeted detection of one or a few transcripts at a once. Differential display and serial analysis of gene expression represent non-targeted methods to evaluate expression changes of a significant number of gene transcripts. Finally, microarrays and RNA-seq (next-generation sequencing) contribute to the ultimate goal of identifying and quantifying all transcripts in a cell under conditions or stages of study. Recent examples of applications as well as principles, advantages, and drawbacks of each method are contrasted. We also suggest replacing the term “Next-Generation Sequencing (NGS)” with another less confusing synonym such as “RNA-seq”, “high throughput sequencing”, or “massively parallel sequencing” to avoid confusion with any future sequencing technologies.

Heat Stress Regulates the Expression of Genes at Transcriptional and Post-Transcriptional Levels, Revealed by RNA-seq in Brachypodium distachyon

Heat stress greatly affects plant growth/development and influences the output of crops. With the increased occurrence of extreme high temperature, the negative influence on cereal products from heat stress becomes severer and severer. It is urgent to reveal the molecular mechanism in response to heat stress in plants. In this research, we used RNA-seq technology to identify differentially expressed genes (DEGs) in leaves of seedlings, leaves and inflorescences at heading stage of Brachypodium distachyon, one model plant of grasses. Results showed many genes in responding to heat stress. Of them, the expression level of 656 DEGs were altered in three groups of samples treated with high temperature. Gene ontology (GO) analysis showed that the highly enriched DEGs were responsible for heat stress and protein folding. According to KEGG pathway analysis, the DEGs were related mainly to photosynthesis-antenna proteins, the endoplasmic reticulum, and the spliceosome. Additionally, the expression level of 454 transcription factors belonging to 49 gene families was altered, as well as 1,973 splicing events occurred after treatment with high temperature. This research lays a foundation for characterizing the molecular mechanism of heat stress response and identifying key genes for those responses in plants. These findings also clearly show that heat stress regulates the expression of genes not only at transcriptional level, but also at post-transcriptional level.