A study of transcriptome in leaf rust infected bread wheat involving seedling resistance gene Lr28

Meta-QTL analysis and identification of candidate genes for multiple-traits associated with spot blotch resistance in bread wheat

In bread wheat, a literature search gave 228 QTLs for six traits, including resistance against spot blotch and the following five other related traits: (i) stay green; (ii) flag leaf senescence; (iii) green leaf area duration; (iv) green leaf area of the main stem; and (v) black point resistance. These QTLs were used for metaQTL (MQTL) analysis. For this purpose, a consensus map with 72,788 markers was prepared; 69 of the above 228 QTLs, which were suitable for MQTL analysis, were projected on the consensus map. This exercise resulted in the identification of 16 meta-QTLs (MQTLs) located on 11 chromosomes, with the PVE ranging from 5.4% (MQTL7) to 21.8% (MQTL5), and the confidence intervals ranging from 1.5 to 20.7 cM (except five MQTLs with a range of 36.1-57.8 cM). The number of QTLs associated with individual MQTLs ranged from a maximum of 17 in MQTL3 to 8 each in MQTL5 and MQTL8 and 5 each in MQTL7 and MQTL14. The 16 MQTLs, included 12 multi-trait MQTLs; one of the MQTL also overlapped a genomic region carrying the major spot blotch resistance gene Sb1. Of the total 16 MQTLs, 12 MQTLs were also validated through marker-trait associations that were available from earlier genome-wide association studies. The genomic regions associated with MQTLs were also used for the identification of candidate genes (CGs) and led to the identification of 516 CGs encoding 508 proteins; 411 of these proteins are known to be associated with resistance against several biotic stresses. In silico expression analysis of CGs using transcriptome data allowed the identification of 71 differentially expressed CGs, which were examined for further possible studies. The findings of the present study should facilitate fine-mapping and cloning of genes, enabling Marker Assisted Selection.

Comparative transcriptomics of stem rust resistance in wheat NILs mediated by Sr24 rust resistance gene

Stem rust of wheat is a deleterious fungal disease across the globe causing severe yield losses. Although, many stem rust resistance genes (Sr) are being used in wheat breeding programs, new emerging stem rust pathotypes are a challenge to important Sr genes. In recent years, multiple studies on leaf and yellow rust molecular mechanism have been done, however, for stem rust such studies are lacking. Current study investigated stem rust induced response in the susceptible wheat genotype C306 and its Near Isogenic Line (NIL) for Sr24 gene, HW2004, using microarray analysis to understand the transcriptomic differences at different stages of infection. Results showed that HW2004 has higher basal levels of several important genes involved in pathogen detection, defence, and display early activation of multiple defence mechanisms. Further Gene Ontology (GO) and pathway analysis identified important genes responsible for pathogen detection, downstream signalling cascades and transcription factors (TFs) involved in activation and mediation of defence responses. Results suggest that generation of Reactive Oxygen Species (ROS), cytoskeletal rearrangement, activation of multiple hydrolases, and lipid metabolism mediated biosynthesis of certain secondary metabolites are collectively involved in Sr24-mediated defence in HW2004, in response to stem rust infection. Novel and unannotated, but highly responsive genes were also identified, which may also contribute towards resistance phenotype. Furthermore, certain DEGs also mapped close to the Sr24-linked marker on Thinopyrum elongatum translocated fragment on wheat 3E chromosome, which advocate further investigations for better insights of the Sr24-mediated stem rust resistance.

Candidate effectors for leaf rust resistance gene Lr28 identified through transcriptome and in-silico analysis

Puccinia spp. causing rust diseases in wheat and other cereals secrete several specialized effector proteins into host cells. Characterization of these proteins and their interaction with host's R proteins could greatly help to limit crop losses due to diseases. Prediction of effector proteins by combining the transcriptome analysis and multiple in-silico approaches is gaining importance in revealing the pathogenic mechanism. The present study involved identification of 13 Puccinia triticina (Pt) coding sequences (CDSs), through transcriptome analysis, that were differentially expressed during wheat-leaf rust interaction; and prediction of their effector like features using different in-silico tools. NCBI-BLAST and pathogen-host interaction BLAST (PHI-BLAST) tools were used to annotate and classify these sequences based on their most closely matched counterpart in both the databases. Homology between CDSs and the annotated sequences in the NCBI database ranged from 79 to 94% and with putative effectors of other plant pathogens in PHI-BLAST from 24.46 to 54.35%. Nine of the 13 CDSs had effector-like features according to EffectorP 3.0 (≥0.546 probability of these sequences to be effector). The qRT-PCR expression analysis revealed that the relative expression of all CDSs in compatible interaction (HD2329) was maximum at 11 days post inoculation (dpi) and that in incompatible interactions (HD2329 + Lr28) was maximum at 3 dpi in seven and 9 dpi in five CDSs. These results suggest that six CDSs (>0.8 effector probability as per EffectorP 3.0) could be considered as putative Pt effectors. The molecular docking and MD simulation analysis of these six CDSs suggested that candidate Lr28 protein binds more strongly to candidate effector c14094_g1_i1 to form more stable complex than the remaining five. Further functional characterization of these six candidate effectors should prove useful for a better understanding of wheat-leaf rust interaction. In turn, this should facilitate effector-based leaf rust resistance breeding in wheat.

Wheat Omics: Advancements and Opportunities

Plant omics, which includes genomics, transcriptomics, metabolomics and proteomics, has played a remarkable role in the discovery of new genes and biomolecules that can be deployed for crop improvement. In wheat, great insights have been gleaned from the utilization of diverse omics approaches for both qualitative and quantitative traits. Especially, a combination of omics approaches has led to significant advances in gene discovery and pathway investigations and in deciphering the essential components of stress responses and yields. Recently, a Wheat Omics database has been developed for wheat which could be used by scientists for further accelerating functional genomics studies. In this review, we have discussed various omics technologies and platforms that have been used in wheat to enhance the understanding of the stress biology of the crop and the molecular mechanisms underlying stress tolerance.

Leaf rust responsive miRNA and their target genes in wheat

Small RNA sequencing (sRNA-seq) and degradome analysis were used for the identification of miRNAs and their target host genes in a pair of near-isogenic lines (NILs), which differed for the presence of leaf rust resistance gene Lr28. The study led to identification of (i) 506 known and 346 novel miRNAs; and (ii) 5054 target genes including 4557 in silico predicted and 497 degradome-based genes using 105 differentially expressed (DE) miRNAs. A subset of 128 targets (67 in silico + 61 degradome-based) was differentially expressed in RNA-seq data that was generated by us earlier using the same pair of NILs; among these 128 targets, 58 target genes exhibited an inverse relationship with the DE miRNAs (expression of miRNAs and activation/suppression of target genes). Eight miRNAs which belonged to the conserved miRNA families and were known to be induced in response to fungal diseases in plants included the following: miR156, miR158, miR159, miR168, miR169, miR172, miR319, miR396. The target genes belonged to the following classes of genes known to be involved in downstream disease resistance pathways; peroxidases, sugar transporters, auxin response signaling, oxidation-reduction, etc. It was also noticed that although a majority of miRNAs and target genes followed the above classical inverse relationship, there were also examples, where no such relationship was observed. Among the target genes, there were also 51 genes that were not only regulated by miRNAs, but were also differentially methylated at sequences including the following segments: promotors, introns, TSS, exons. The results of the present study suggest a complex interplay among miRNA genes, target genes, and various epigenetic controls, which regulate the expression of genes involved in downstream pathways for disease resistance.

BS-Seq reveals major role of differential CHH methylation during leaf rust resistance in wheat (Triticum aestivum L.)

Epigenetic regulation of the activity of defense genes during onset of diseases or resistance against diseases in plants is an active area of research. In the present study, a pair of wheat NILs for leaf rust resistance gene Lr28 (R) in the background of an Indian cultivar HD2329 (S) was used for a study of DNA methylation mediated regulation of gene expression. Leaf samples were collected at 0 h before (S0 and R0) and 96 h after inoculation (S96 and R96). The DNA samples were subjected to BS-Seq and sequencing data were used for identification of differentially methylated/demethylated regions/genes (DMRs and DMGs). Following four pairs of comparisons were used for this purpose: S0 vs S96; S0 vs R0; R0 vs R96; S96 vs R96. Major role of CHH methylation relative to that of CG and CHG methylation was observed. Some important observations include the following: (i) abundance of CHH methylation among DMRs; (ii) predominance of DMRs in intergenic region, relative to other genomic regions (promoters, exons, introns, TSS and TTS); (iii) abundance of transposable elements (TEs) in DMRs with CHH context; (iv) demethylation mediated high expression of genes during susceptible reaction (S0 vs S96) and methylation mediated low expression of genes during resistant reaction (R0 vs R96 and S96 vs R96); (v) major genes under regulation encode proteins, which differ from those encoded by genes regulated during susceptible reaction and (vi) ~ 500 DMGs carried differential binding sites for H3K4/K27me3 marks suggesting joint involvement of DNA and H3 methylation. Thus, CHH methylation either alone or in combination with histone methylation plays a major role in regulating the expression of genes involved in wheat-leaf rust interaction.

Omics for the Improvement of Abiotic, Biotic, and Agronomic Traits in Major Cereal Crops: Applications, Challenges, and Prospects

Omics technologies, namely genomics, transcriptomics, proteomics, metabolomics, and phenomics, are becoming an integral part of virtually every commercial cereal crop breeding program, as they provide substantial dividends per unit time in both pre-breeding and breeding phases. Continuous advances in omics assure time efficiency and cost benefits to improve cereal crops. This review provides a comprehensive overview of the established omics methods in five major cereals, namely rice, sorghum, maize, barley, and bread wheat. We cover the evolution of technologies in each omics section independently and concentrate on their use to improve economically important agronomic as well as biotic and abiotic stress-related traits. Advancements in the (1) identification, mapping, and sequencing of molecular/structural variants; (2) high-density transcriptomics data to study gene expression patterns; (3) global and targeted proteome profiling to study protein structure and interaction; (4) metabolomic profiling to quantify organ-level, small-density metabolites, and their composition; and (5) high-resolution, high-throughput, image-based phenomics approaches are surveyed in this review.

Transcriptome-based analysis of resistance mechanism to black point caused by Bipolaris sorokiniana in wheat

Black point is a cereal disease caused by complex pathogens, of which the pathogenicity of Bipolaris sorokiniana is the most serious in wheat. Resistance to black point is quantitative in nature, and thus the mechanism is poorly understood. We conducted a comparative transcriptome analysis to identify differentially expressed genes (DEGs) in black point-slightly susceptible and -highly susceptible wheat lines at different timepoints following B. sorokiniana inoculation. DEGs associated with photosynthesis were upregulated in black point-slightly susceptible lines. The top Gene Ontology enrichment terms for biological processes were oxidation–reduction, response to cold, salt stress, oxidative stress, and cadmium ion; terms for cellular component genes were mainly involved in plasma membrane and cytoplasmic membrane-bounded vesicle, whereas those for molecular function were heme binding and peroxidase activity. Moreover, activities of antioxidant enzymes superoxide dismutase, catalase, and peroxidase were higher in slightly susceptible lines than those in highly susceptible lines (except peroxidase 12–24 days post-inoculation). Thus, resistance to B. sorokiniana-caused black point in wheat was mainly related to counteracting oxidative stress, although the specific metabolic pathways require further study. This study presents new insights for understanding resistance mechanisms of selected wheat lines to black point.

QTL mapping of adult plant and seedling resistance to leaf rust (Puccinia triticina Eriks.) in a multiparent advanced generation intercross (MAGIC) wheat population

The Bavarian MAGIC Wheat population, comprising 394 F6:8 recombinant inbred lines was phenotyped for Puccinia triticina resistance in multi-years' field trials at three locations and in a controlled environment seedling test. Simple intervall mapping revealed 19 QTL, corresponding to 11 distinct chromosomal regions. The biotrophic rust fungus Puccinia triticina is one of the most important wheat pathogens with the potential to cause yield losses up to 70%. Growing resistant cultivars is the most cost-effective and environmentally friendly way to encounter this problem. The emergence of leaf rust races being virulent against common resistance genes increases the demand for wheat varieties with novel resistances. In the past decade, the use of complex experimental populations, like multiparent advanced generation intercross (MAGIC) populations, has risen and offers great advantages for mapping resistances. The genetic diversity of multiple parents, which has been recombined over several generations, leads to a broad phenotypic diversity, suitable for high-resolution mapping of quantitative traits. In this study, interval mapping was performed to map quantitative trait loci (QTL) for leaf rust resistance in the Bavarian MAGIC Wheat population, comprising 394 F_6:8 recombinant inbred lines (RILs). Phenotypic evaluation of the RILs for adult plant resistance was carried out in field trials at three locations and two years, as well as in a controlled-environment seedling inoculation test. In total, interval mapping revealed 19 QTL, which corresponded to 11 distinct chromosomal regions controlling leaf rust resistance. Six of these regions may represent putative new QTL. Due to the elite parental material, RILs identified to be resistant to leaf rust can be easily introduced in breeding programs.

Genome-wide analysis of H3K4me3 and H3K27me3 modifications due to Lr28 for leaf rust resistance in bread wheat (Triticum aestivum)

Present study revealed a complex relationship among histone H3 methylation (examined using H3K4/K27me3 marks), cytosine DNA methylation and differential gene expression during Lr28 mediated leaf rust resistance in wheat. During the present study, genome-wide histone modifications were examined in a pair of near isogenic lines (NILs) (with and without Lr28 in the background of cv. HD2329). The two histone marks used included H3K4me3 (an activation mark) and H3K27me3 (a repression mark). The results were compared with levels of expression (using RNA-seq) and DNA methylation (MeDIP) data obtained using the same pair of NILs. Some of the salient features of the present study include the following: (i) large scale differential binding sites (DBS) were available for only H3K4me3 in the susceptible cultivar, but for both H3K4me3 and H3K27me3 in its resistant NIL; (ii) DBSs for H3K27me3 mark were more abundant (> 80%) in intergenic regions, whereas DBSs for H3K4me3 were distributed in all genomic regions including exons, introns, intergenic, TTS (transcription termination sites) and promoters; (iii) fourteen (14) genes associated with DBSs showed co-localization for both the marks; (iv) only a small fraction (7% for H3K4me3 and 12% for H3K27me3) of genes associated with DBSs matched with the levels of gene expression inferred from RNA-seq data; (v) validation studies using qRT-PCR were conducted on 26 selected representative genes; results for only 11 genes could be validated. The proteins encoded by important genes involved in promoting infection included domains generally carried by R gene proteins such as Mlo like protein, protein kinases and purple acid phosphatase. Similarly, proteins encoded by genes involved in resistance included those carrying domains for lectin kinase, R gene, aspartyl protease, etc. Overall, the results suggest a very complex network of downstream genes that are expressed during compatible and incompatible interactions; some of the genes identified during the present study may be used in future validation studies involving RNAi/overexpression approaches.

The progress of leaf rust research in wheat

Leaf rust (also called brown rust) in wheat, caused by fungal pathogen Puccinia triticina Erikss. (Pt) is one of the major constraints in wheat production worldwide. Pt is widespread with diverse population structure and undergoes rapid evolution to produce new virulent races against resistant cultivars that are regularly developed to provide resistance against the prevailing races of the pathogen. Occasionally, the disease may also take the shape of an epidemic in some wheat-growing areas causing major economic losses. In the recent past, substantial progress has been made in characterizing the sources of leaf rust resistance including non-host resistance (NHR). Progress has also been made in elucidating the population biology of Pt and the mechanisms of wheat-Pt interaction. So far, ∼80 leaf rust resistance genes (Lr genes) have been identified and characterized; some of them have also been used for the development of resistant wheat cultivars. It has also been shown that a gene-for-gene relationship exists between individual wheat Lr genes and the corresponding Pt Avr genes so that no Lr gene can provide resistance unless the prevailing race of the pathogen carries the corresponding Avr gene. Several Lr genes have also been cloned and their products characterized, although no Avr gene corresponding a specific Lr gene has so far been identified. However, several candidate effectors for Pt have been identified and functionally characterized using genome-wide analyses, transcriptomics, RNA sequencing, bimolecular fluorescence complementation (BiFC), virus-induced gene silencing (VIGS), transient expression and other approaches. This review summarizes available information on different aspects of the pathogen Pt, genetics/genomics of leaf rust resistance in wheat including cloning and characterization of Lr genes and epigenetic regulation of disease resistance.

Large-scale stage-specific regulation of gene expression during host-pathogen interactions in CSP44 bread wheat carrying APR gene Lr48

Genome-wide transcriptome analysis was undertaken in a leaf-rust resistant bread wheat line CSP44 (selected from Australian cv. Condor) carrying the adult plant resistance (APR) gene Lr48. Two pre-adult plant (P-AP) susceptible stages (S48 and S96) and two adult plant (AP) resistant stages (R48 and R96) were used for RNA-seq. At the susceptible P-AP stage (during S48 to S96), expression increased in 2062 genes, and declined in 130 genes; 1775 of 2062 differentially expressed genes (DEGs) also exhibited high expression during early incompatible stage R48. Comparison of S96 with R96 showed that the expression of 80 genes was enhanced and that of 208 genes declined at the AP stage. At the resistant AP stage (during R48 to R96), expression of mere 25 genes increased and that of 126 genes declined. Apparently, the resistance during late adult stage (R96) is caused by regulation of the expression of relatively fewer genes, although at pre-adult stage (S48 to S96), expression of large number of genes increased; expression of majority of these genes kept on increasing during adult stage at R48 also. These and other results of the present study suggest that APR may mimic some kind of systemic acquired resistance (SAR). The host-specific DEGs belonged to 10 different classes including genes involved in defence, transport, epigenetics, photosynthesis, genes encoding some transcription factors etc. The pathogen (Puccinia triticina) specific DEGs (including three genes encoding known biotrophic effectors) seem to help the pathogen in infection/growth through large-scale stage-specific enhanced expression of host's genes. A putative candidate gene for Lr48 containing protein kinase domain (its ortholog in rice encoding OsWAK8) was also identified.

Further studies on sugar transporter (SWEET) genes in wheat (Triticum aestivum L.)

SWEET proteins represent one of the largest sugar transporter family in the plant kingdom and play crucial roles in plant development and stress responses. In the present study, a total of 108 TaSWEET genes distributed on all the 21 wheat chromosomes were identified using the latest whole genome sequence (as against 59 genes reported in an earlier report). These 108 genes included 14 of the 17 types reported in Arabidopsis and also included three novel types. Tandem duplications (22) and segmental duplications (5) played a significant role in the expansion of TaSWEET family. A number of cis-elements were also identified in the promoter regions of TaSWEET genes, indicating response of TaSWEET genes during development and also during biotic/abiotic stresses. The TaSWEET proteins carried 4-7 trans-membrane helices (TMHs) showing diversity in structure. Phylogenetic analysis using SWEET proteins of wheat and 8 other species gave four well-known clusters. Expression analysis involving both in silico and in planta indicated relatively higher expression of TaSWEET genes in water/heat sensitive and leaf rust resistant genotypes. The results provided insights into the functional role of TaSWEETs in biotic and abiotic stresses, which may further help in planning strategies to develop high yielding wheat varieties tolerant to environmental stresses.

H3K4/K9 acetylation and Lr28-mediated expression of six leaf rust responsive genes in wheat (Triticum aestivum)

Development of leaf rust-resistant cultivars is a priority during wheat breeding, since leaf rust causes major losses in yield. Resistance against leaf rust due to Lr genes is partly controlled by epigenetic modifications including histone acetylation that is known to respond to biotic/abiotic stresses. In the present study, enrichment of H3K4ac and H3K9ac in promoters of six defense responsive genes (N-acetyltransferase, WRKY 40, WRKY 70, ASR1, Peroxidase 12 and Sarcosine oxidase) was compared with their expression in a pair of near-isogenic lines (NILs) for the gene Lr28 following inoculation with leaf rust pathotype '77-5'; ChIP-qPCR was used for this purpose. The proximal and distal promoters of these genes contained a number of motifs that are known to respond to biotic stresses. The enrichment of two acetylation marks changed with passage of time; changes in expression of two of the six genes (N-acetyltransferase and peroxidase12), largely matched with changes in H3K4/H3K9 acetylation patterns of the two promoter regions. For example, enrichment of both the marks matched with higher expression of N-acetyltransferase gene in susceptible NIL and the deacetylation (H3K4ac) largely matched with reduced gene expression in resistant NIL. In peroxidase12, enrichment of H3K4ac and H3K9ac largely matched with higher expression in both the NILs. In the remaining four genes, changes in H3 acetylation did not always match with gene expression levels. This indicated complexity in the regulation of the expression of these remaining four genes, which may be controlled by other epigenetic/genetic regulatory mechanisms that need further analysis.