AMMI and GGE biplot analysis of yield under terminal heat tolerance in wheat.

BACKGROUND: Wheat is an important cereal crop that helps to meet the food grain needs of people all over the world. Heat stress is one of the most significant abiotic stresses that wheat crops face during terminal growth stages in the wheat growing regions like India. It is very important to identify heat tolerant genotypes to be used as donors for breeding tolerant varieties. METHODS: Thirty-six wheat genotypes were evaluated under different sowing dates viz., Timely sown (TS), Late sown (LS) and very late sown (VLS), and the fourth was sown in the Temperature controlled phenotyping facility (TCPF) across two years. Genotypes were planted following lattice square design with two replications. Data was recorded for yield and yield contributing traits and analysed using selection indices as well AMMI and GGE biplot stability models. RESULTS: Heat stress affected all the traits under different heat environments which ranged from 1.6% (Spikelet number) to 37.2% (grain yield). Regression analysis indicated that the thousand grains weight (R2 = 0.50) contributed significantly towards grain yield under heat stress. Stress susceptibility index (SSI) found genotypes GW322, RAJ3765, Raj4037and MACS6145 as heat tolerant whereas, Stress Tolerance Index (STI) identified C306, HD2967, WH1080, WH730, DBW90, HD2932, DBW17, RAJ3765 as heat tolerant and high yielding. AMMI biplot analysis indicated stable genotypes DBW90, WH730, RAJ4083, CBW38, HD2932, NI5439, WR544, whereas GGE biplot analysis revealed stable genotypes NIAW34, NI5439, RAJ4083, DBW90, PBW590, Raj3765, HUW 510, WH730, HD2967 and UP2382. CONCLUSION: Heat stress affects significantly all yield contributing traits. Thousand grain weight was the most important trait that can be used as a selection criterion for selecting tolerant lines. Based on selection indices and both AMMI and GGE analysis, genotype RAJ3765 was identified to be highly heat tolerant with good grain yield.

Lr80: A new and widely effective source of leaf rust resistance of wheat for enhancing diversity of resistance among modern cultivars.

KEY MESSAGE: A new leaf rust resistance gene Lr80 was identified and closely linked markers were developed for its successful pyramiding with other marker-tagged genes to achieve durable control of leaf rust. Common wheat landrace Hango-2, collected in 2006 from the Himalayan area of Hango, District Kinnaur, in Himachal Pradesh, exhibited a very low infection type (IT;) at the seedling stage to all Indian Puccinia triticina (Pt) pathotypes, except the pathotype 5R9-7 which produced IT 3+. Genetic analysis based on Agra Local/Hango-2-derived F3 families indicated monogenic control of leaf rust resistance, and the underlying locus was temporarily named LrH2. Bulked segregant analysis using 303 simple sequence repeat (SSR) markers located LrH2 in the short arm of chromosome 2D. An additional set of 10 chromosome 2DS-specific markers showed polymorphism between the parents and these were mapped on the entire Agra Local/Hango-2 F3 population. LrH2 was flanked by markers cau96 (distally) and barc124 (proximally). The 90 K Infinium SNP array was used to identify SNP markers linked with LrH2. Markers KASP_17425 and KASP_17148 showed association with LrH2. Comparison of seedling leaf rust response data and marker locations across different maps demonstrated the uniqueness of LrH2 and it was formally named Lr80. The Lr80-linked markers KASP_17425, KASP_17148 and barc124 amplified alleles/products different to Hango-2 in 82 Australian cultivars indicating their robustness for marker-assisted selection of this gene in wheat breeding programs.

Quantitative proteomic analysis reveals novel stress-associated active proteins (SAAPs) and pathways involved in modulating tolerance of wheat under terminal heat.

Terminal heat stress has detrimental effect on the growth and yield of wheat. Very limited information is available on heat stress-associated active proteins (SAAPs) in wheat. Here, we have identified 159 protein groups with 4271 SAAPs in control (22 ± 3 °C) and HS-treated (38 °C, 2 h) wheat cvs. HD2985 and HD2329 using iTRAQ. We identified 3600 proteins to be upregulated and 5825 proteins to be downregulated in both the wheat cvs. under HS. We observed 60.3% of the common SAAPs showing upregulation in HD2985 (thermotolerant) and downregulation in HD2329 (thermosusceptible) under HS. GO analysis showed proton transport (molecular), photosynthesis (biological), and ATP binding (cellular) to be most altered under HS. Most of the SAAPs identified were observed to be chloroplast localized and involved in photosynthesis. Carboxylase enzyme was observed most abundant active enzymes in wheat under HS. An increase in the degradative isoenzymes (α/ß-amylases) was observed, as compared to biosynthesis enzymes (ADP-glucophosphorylase, soluble starch synthase, etc.) under HS. Transcript profiling showed very high relative fold expression of HSP17, CDPK, Cu/Zn SOD, whereas downregulation of AGPase, SSS under HS. The identified SAAPs can be used for targeted protein-based precision wheat-breeding program for the development of 'climate-smart' wheat.

Gamma irradiation protect the developing wheat endosperm from oxidative damage by balancing the trade-off between the defence network and grains quality.

Gamma irradiation has been reported to modulate the biochemical and molecular parameters associated with the tolerance of plant species under biotic/ abiotic stress. Wheat is highly sensitive to heat stress (HS), as evident from the decrease in the quantity and quality of the total grains. Here, we studied the effect of pre-treatment of wheat dry seeds with different doses of gamma irradiation (0.20, 0.25 and 0.30 kGy) on tolerance level and quality of developing wheat endospermic tissue under HS (38 °C, 1 h; continuously for three days). Expression analysis of genes associated with defence and starch metabolism in developing grains showed maximum transcripts of HSP17 (in response to 0.25 kGy + HS) and AGPase (under 0.30 kGy), as compared to control. Gamma irradiation was observed to balance the accumulation of H2O2 by enhancing the activities of SOD and GPx in both the cvs. under HS. Gamma irradiation was observed to stabilize the synthesis of starch and amylose by regulating the activities of AGPase, SSS and α-amylase under HS. The appearance of isoforms of gliadins (α, ß, γ, ω) were observed more in gamma irradiated seeds (0.20 kGy), as compared to control. Gamma irradiation (0.25 kGy in HD3118 & 0.20 kGy in HD3086) was observed to have positive effect on the width, length and test seed weight of the grains under HS. The information generated in present investigation provides easy, cheap and user-friendly technology to mitigate the effect of terminal HS on the grain-development process of wheat along with development of robust seeds with high nutrient density.

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.

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.

X-ray knee as a screening tool for osteoporosis.

Cortical thickness (Cor-Th) of tibia varies considerably on X-ray knees. It was hypothesized that Cor-Th can be used for preliminary prediction of BMD. Ninety nine patients underwent a digital X-ray left knee fixed flexion PA view with an external calibration scale attached to X-ray plate and BMD by DXA using GE lunar machine (Madison, Wisconsin.). Cor-Th was measured at 5 selected levels (A,B,C,D, and E) ranging from 5-7 cm below the tibial plateau on its medial aspect. T-scores were recorded for BMD at AP spine, left forearm and left femur. Cor-Th of tibia at each level significantly correlated with each site of BMD measurement namely AP spine, left femur and left forearm. This correlation varied in the range from 0.241 to 0.426. For AP spine, it was maximum at level C (r=0.347, p<0.001) whereas for left femur and forearm sites, it was maximum at level B (r=0.426 &r=0.373 respectively, p<0.001). The correlation of Cor-Th with BMD varied with age. Above 56 years of age, Cor-Th at each level significantly correlated to BMD at each site. Medial tibial cortical thickness, 6 cm (level C) below tibial plateau can be used as preliminary predictor of patients who need a DXA scan.

Identification of Genomic Regions and Sources for Wheat Blast Resistance through GWAS in Indian Wheat Genotypes.

Wheat blast (WB) is a devastating fungal disease that has recently spread to Bangladesh and poses a threat to the wheat production in India, which is the second-largest wheat producing country in the world. In this study, 350 Indian wheat genotypes were evaluated for WB resistance in 12 field experiments in three different locations, namely Jashore in Bangladesh and Quirusillas and Okinawa in Bolivia. Single nucleotide polymorphisms (SNPs) across the genome were obtained using DArTseq® technology, and 7554 filtered SNP markers were selected for a genome-wide association study (GWAS). All the three GWAS approaches used identified the 2NS translocation as the only major source of resistance, explaining up to 32% of the phenotypic variation. Additional marker-trait associations were located on chromosomes 2B, 3B, 4D, 5A and 7A, and the combined effect of three SNPs (2B_180938790, 7A_752501634 and 5A_618682953) showed better resistance, indicating their additive effects on WB resistance. Among the 298 bread wheat genotypes, 89 (29.9%) carried the 2NS translocation, the majority of which (60 genotypes) were CIMMYT introductions, and 29 were from India. The 2NS carriers with a grand mean WB index of 6.6 showed higher blast resistance compared to the non-2NS genotypes with a mean index of 46.5. Of the 52 durum wheats, only one genotype, HI 8819, had the 2NS translocation and was the most resistant, with a grand mean WB index of 0.93. Our study suggests that the 2NS translocation is the only major resistance source in the Indian wheat panel analysed and emphasizes the urgent need to identify novel non-2NS resistance sources and genomic regions.

Screening and Mapping for Head Blast Resistance in a Panel of CIMMYT and South Asian Bread Wheat Germplasm.

Wheat blast (WB) is a destructive disease in South America and its first outbreak in Bangladesh in 2016 posed a great risk to food security of South Asian countries. A genome wide association study (GWAS) was conducted on a diverse panel of 184 wheat genotypes from South Asia and CIMMYT. Phenotyping was conducted in eight field experiments in Bolivia and Bangladesh and a greenhouse experiment in the United States. Genotypic data included 11,401 SNP markers of the Illumina Infinium 15K BeadChip and four additional STS markers on the 2NS/2AS translocation region. Accessions with stable WB resistance across experiments were identified, which were all 2NS carriers. Nevertheless, a dozen moderately resistant 2AS lines were identified, exhibiting big variation among experiments. Significant marker-trait associations (MTA) were detected on chromosomes 1BS, 2AS, 6BS, and 7BL; but only MTAs on 2AS at the 2NS/2AS translocation region were consistently significant across experiments. The resistant accessions identified in this study could be used in production in South Asian countries as a preemptive strategy to prevent WB outbreak.

New Genotypes and Genomic Regions for Resistance to Wheat Blast in South Asian Germplasm.

Wheat blast (WB) disease, since its first identification in Bangladesh in 2016, is now an established serious threat to wheat production in South Asia. There is a need for sound knowledge about resistance sources and associated genomic regions to assist breeding programs. Hence, a panel of genotypes from India and Bangladesh was evaluated for wheat blast resistance and a genome-wide association study (GWAS) was performed. Disease evaluation was done during five crop seasons-at precision phenotyping platform (PPPs) for wheat blast disease at Jashore (2018-19), Quirusillas (2018-19 and 2019-20) and Okinawa (2019 and 2020). Single nucleotide polymorphisms (SNP) across the genome were obtained using DArTseq genotyping-by-sequencing platform, and in total 5713 filtered markers were used. GWAS revealed 40 significant markers associated with WB resistance, of which 33 (82.5%) were in the 2NS/2AS chromosome segment and one each on seven chromosomes (3B, 3D, 4A, 5A, 5D, 6A and 6B). The 2NS markers contributed significantly in most of the environments, explaining an average of 33.4% of the phenotypic variation. Overall, 22.4% of the germplasm carried 2NS/2AS segment. So far, 2NS translocation is the only effective WB resistance source being used in the breeding programs of South Asia. Nevertheless, the identification of non-2NS/2AS genomic regions for WB resistance provides a hope to broaden and diversify resistance for this disease in years to come.

Genomic Selection for Wheat Blast in a Diversity Panel, Breeding Panel and Full-Sibs Panel.

Wheat blast is an emerging threat to wheat production, due to its recent migration to South Asia and Sub-Saharan Africa. Because genomic selection (GS) has emerged as a promising breeding strategy, the key objective of this study was to evaluate it for wheat blast phenotyped at precision phenotyping platforms in Quirusillas (Bolivia), Okinawa (Bolivia) and Jashore (Bangladesh) using three panels: (i) a diversity panel comprising 172 diverse spring wheat genotypes, (ii) a breeding panel comprising 248 elite breeding lines, and (iii) a full-sibs panel comprising 298 full-sibs. We evaluated two genomic prediction models (the genomic best linear unbiased prediction or GBLUP model and the Bayes B model) and compared the genomic prediction accuracies with accuracies from a fixed effects model (with selected blast-associated markers as fixed effects), a GBLUP + fixed effects model and a pedigree relationships-based model (ABLUP). On average, across all the panels and environments analyzed, the GBLUP + fixed effects model (0.63 ± 0.13) and the fixed effects model (0.62 ± 0.13) gave the highest prediction accuracies, followed by the Bayes B (0.59 ± 0.11), GBLUP (0.55 ± 0.1), and ABLUP (0.48 ± 0.06) models. The high prediction accuracies from the fixed effects model resulted from the markers tagging the 2NS translocation that had a large effect on blast in all the panels. This implies that in environments where the 2NS translocation-based blast resistance is effective, genotyping one to few markers tagging the translocation is sufficient to predict the blast response and genome-wide markers may not be needed. We also observed that marker-assisted selection (MAS) based on a few blast-associated markers outperformed GS as it selected the highest mean percentage (88.5%) of lines also selected by phenotypic selection and discarded the highest mean percentage of lines (91.8%) also discarded by phenotypic selection, across all panels. In conclusion, while this study demonstrates that MAS might be a powerful strategy to select for the 2NS translocation-based blast resistance, we emphasize that further efforts to use genomic tools to identify non-2NS translocation-based blast resistance are critical.

Characterization of the starch synthase under terminal heat stress and its effect on grain quality of wheat.

Terminal heat stress (HS) is a key barrier for wheat grain yield and quality. Various physiochemical and molecular parameters such as photosynthetic rate, expression analysis and activity of starch synthase (SS), total starch, amylose and amylopectin content, total amylolytic activity, and total antioxidant capacity (TAC) were analysed in wheat cvs.HD3059 (thermotolerant) and BT-Schomburgk (thermosusceptible) at grain-filling stage under HS (32 °C and 40 °C, 1 h). The decrease in photosynthetic rate was observed under HS. Expression analysis of the SS gene at transcript level showed downregulation in both the wheat cvs.HD3059 and BT-Schombugk under HS (32 °C and 40 ºC, 1 h) as compared to the control. Although the downregulation of SS gene transcript expression was less in HD3059 than BT-Schombugk. Both the cultivars showed decrease in starch synthase activity and starch content under HS and the overall content was higher in HD3059, compared to BT-Schomburgk. Higher total amylolytic activity and amylose content were observed in BT-Schomburgk. Scanning electron microscopy (SEM) showed un-structured starch granules under HS. Total antioxidant capacity (TAC) was found higher in HD3059 (14.07 mM FeSO4 gm-1 FW) compared to BT-Schomburgk (8.89 mM FeSO4 gm-1 FW) under HS (40 ºC, 1 h). Findings suggest that HS during grain filling stage had more severe impact on the overall physiochemical properties of the wheat grain. Thus the starch bisynthesis pathway associated gene(s) could be exploit to enhance the yield and quality of wheat under heat stress.

Heterologous expression and characterization of novel manganese superoxide dismutase (Mn-SOD) - A potential biochemical marker for heat stress-tolerance in wheat (Triticum aestivum).

Heat stress causes oxidative bursts damaging the organelles and nascent proteins. Plants have inherited antioxidant defense system to neutralize the effect of reactive oxygen species. Superoxide dismutase provides first line of defense against the HS by regulating the accumulation of peroxide radicals inside the cells. Here, we report identification and cloning of putative manganese superoxide dismutase (Mn-SOD) gene of ~733 nt from wheat cv. HD2985 through de novo assembly. The gene was observed to localize on Chr 6D with a mitochondrial targeting peptide sequence and iron/manganese domain. We predicted 147 homologs of Mn-SOD in eukaryotes with diverse speciation nodes. A recombinant Mn-SOD protein of ~25.5 kDa was purified through heterologous expression system. Kinetics assay of recombinant protein showed optimum pH of 8.0, optimum temperature of 35 °C and Km and Vmax values of 1.51 µM and 9.45 U/mg proteins, respectively. Maximum expression and activity of Mn-SOD was observed in leaves from Raj3765, as compared to stem and spike during milky-ripe stage under differential HS. In gel activity assay showed the appearance of all the three isoforms of SOD in thermotolerant cv. under HS. Mn-SOD, being active at pivotal position, can be also used as potential biochemical marker in wheat breeding program.

Harnessing genetic potential of wheat germplasm banks through impact-oriented-prebreeding for future food and nutritional security.

The value of exotic wheat genetic resources for accelerating grain yield gains is largely unproven and unrealized. We used next-generation sequencing, together with multi-environment phenotyping, to study the contribution of exotic genomes to 984 three-way-cross-derived (exotic/elite1//elite2) pre-breeding lines (PBLs). Genomic characterization of these lines with haplotype map-based and SNP marker approaches revealed exotic specific imprints of 16.1 to 25.1%, which compares to theoretical expectation of 25%. A rare and favorable haplotype (GT) with 0.4% frequency in gene bank identified on chromosome 6D minimized grain yield (GY) loss under heat stress without GY penalty under irrigated conditions. More specifically, the 'T' allele of the haplotype GT originated in Aegilops tauschii and was absent in all elite lines used in study. In silico analysis of the SNP showed hits with a candidate gene coding for isoflavone reductase IRL-like protein in Ae. tauschii. Rare haplotypes were also identified on chromosomes 1A, 6A and 2B effective against abiotic/biotic stresses. Results demonstrate positive contributions of exotic germplasm to PBLs derived from crosses of exotics with CIMMYT's best elite lines. This is a major impact-oriented pre-breeding effort at CIMMYT, resulting in large-scale development of PBLs for deployment in breeding programs addressing food security under climate change scenarios.

Correction to: Strategic crossing of biomass and harvest index-source and sink-achieves genetic gains in wheat.

[This corrects the article DOI: 10.1007/s10681-017-2040-z.].

Biochemical Defense Response: Characterizing the Plasticity of Source and Sink in Spring Wheat under Terminal Heat Stress.

Wheat is highly prone to terminal heat stress (HS) under late-sown conditions. Delayed- sowing is one of the preferred methods to screen the genotypes for thermotolerance under open field conditions. We investigated the effect of terminal HS on the thermotolerance of four popular genotypes of wheat i.e. WR544, HD2967, HD2932, and HD2285 under field condition. We observed significant variations in the biochemical parameters like protein content, antioxidant activity, proline and total reducing sugar content in leaf, stem, and spike under normal (26 ± 2°C) and terminal HS (36 ± 2°C) conditions. Maximum protein, sugars and proline was observed in HD2967, as compared to other cultivars under terminal HS. Wheat cv. HD2967 showed more adaptability to the terminal HS. Differential protein-profiling in leaves, stem and spike of HD2967 under normal (26 ± 2°C) and terminal HS (36 ± 2°C) showed expression of some unique protein spots. MALDI-TOF/MS analysis showed the DEPs as RuBisCO (Rub), RuBisCO activase (Rca), oxygen evolving enhancer protein (OEEP), hypothetical proteins, etc. Expression analysis of genes associated with photosynthesis (Rub and Rca) and starch biosynthesis pathway (AGPase, SSS and SBE) showed significant variations in the expression under terminal HS. HD2967 showed better performance, as compared to other cultivars under terminal HS. SSS activity observed in HD2967 showed more stability under terminal HS, as compared with other cultivars. Triggering of different biochemical parameters in response to terminal HS was observed to modulate the plasticity of carbon assimilatory pathway. The identified DEPs will enrich the proteomic resources of wheat and will provide a potential biochemical marker for screening wheat germplasm for thermotolerance. The model hypothesized will help the researchers to work in a more focused way to develop terminal heat tolerant wheat without compromising with the quality and quantity of grains.

Genomic Selection in the Era of Next Generation Sequencing for Complex Traits in Plant Breeding.

Genomic selection (GS) is a promising approach exploiting molecular genetic markers to design novel breeding programs and to develop new markers-based models for genetic evaluation. In plant breeding, it provides opportunities to increase genetic gain of complex traits per unit time and cost. The cost-benefit balance was an important consideration for GS to work in crop plants. Availability of genome-wide high-throughput, cost-effective and flexible markers, having low ascertainment bias, suitable for large population size as well for both model and non-model crop species with or without the reference genome sequence was the most important factor for its successful and effective implementation in crop species. These factors were the major limitations to earlier marker systems viz., SSR and array-based, and was unimaginable before the availability of next-generation sequencing (NGS) technologies which have provided novel SNP genotyping platforms especially the genotyping by sequencing. These marker technologies have changed the entire scenario of marker applications and made the use of GS a routine work for crop improvement in both model and non-model crop species. The NGS-based genotyping have increased genomic-estimated breeding value prediction accuracies over other established marker platform in cereals and other crop species, and made the dream of GS true in crop breeding. But to harness the true benefits from GS, these marker technologies will be combined with high-throughput phenotyping for achieving the valuable genetic gain from complex traits. Moreover, the continuous decline in sequencing cost will make the WGS feasible and cost effective for GS in near future. Till that time matures the targeted sequencing seems to be more cost-effective option for large scale marker discovery and GS, particularly in case of large and un-decoded genomes.

SSH Analysis of Endosperm Transcripts and Characterization of Heat Stress Regulated Expressed Sequence Tags in Bread Wheat.

Heat stress is one of the major problems in agriculturally important cereal crops, especially wheat. Here, we have constructed a subtracted cDNA library from the endosperm of HS-treated (42°C for 2 h) wheat cv. HD2985 by suppression subtractive hybridization (SSH). We identified ~550 recombinant clones ranging from 200 to 500 bp with an average size of 300 bp. Sanger's sequencing was performed with 205 positive clones to generate the differentially expressed sequence tags (ESTs). Most of the ESTs were observed to be localized on the long arm of chromosome 2A and associated with heat stress tolerance and metabolic pathways. Identified ESTs were BLAST search using Ensemble, TriFLD, and TIGR databases and the predicted CDS were translated and aligned with the protein sequences available in pfam and InterProScan 5 databases to predict the differentially expressed proteins (DEPs). We observed eight different types of post-translational modifications (PTMs) in the DEPs corresponds to the cloned ESTs-147 sites with phosphorylation, 21 sites with sumoylation, 237 with palmitoylation, 96 sites with S-nitrosylation, 3066 calpain cleavage sites, and 103 tyrosine nitration sites, predicted to sense the heat stress and regulate the expression of stress genes. Twelve DEPs were observed to have transmembrane helixes (TMH) in their structure, predicted to play the role of sensors of HS. Quantitative Real-Time PCR of randomly selected ESTs showed very high relative expression of HSP17 under HS; up-regulation was observed more in wheat cv. HD2985 (thermotolerant), as compared to HD2329 (thermosusceptible) during grain-filling. The abundance of transcripts was further validated through northern blot analysis. The ESTs and their corresponding DEPs can be used as molecular marker for screening or targeted precision breeding program. PTMs identified in the DEPs can be used to elucidate the thermotolerance mechanism of wheat-a novel step toward the development of "climate-smart" wheat.

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

Wheat is a staple food worldwide and provides 40% of the calories in the diet. Climate change and global warming pose a threat to wheat production, however, and demand a deeper understanding of how heat stress might impact wheat production and wheat biology. However, it is difficult to identify novel heat stress associated genes when the genomic information is not available. Wheat has a very large and complex genome that is about 37 times the size of the rice genome. The present study sequenced the whole transcriptome of the wheat cv. HD2329 at the flowering stage, under control (22°±3°C) and heat stress (42°C, 2 h) conditions using Illumina HiSeq and Roche GS-FLX 454 platforms. We assembled more than 26.3 and 25.6 million high-quality reads from the control and HS-treated tissues transcriptome sequences respectively. About 76,556 (control) and 54,033 (HS-treated) contigs were assembled and annotated de novo using different assemblers and a total of 21,529 unigenes were obtained. Gene expression profile showed significant differential expression of 1525 transcripts under heat stress, of which 27 transcripts showed very high (>10) fold upregulation. Cellular processes such as metabolic processes, protein phosphorylation, oxidations-reductions, among others were highly influenced by heat stress. In summary, these observations significantly enrich the transcript dataset of wheat available on public domain and show a de novo approach to discover the heat-responsive transcripts of wheat, which can accelerate the progress of wheat stress-genomics as well as the course of wheat breeding programs in the era of climate change.