Genome-Wide Expression Profiling of Genes Associated with the Lr47-Mediated Wheat Resistance to Leaf Rust (Puccinia triticina)

Puccinia triticina (Pt), the causal agent of wheat leaf rust, is one of the most destructive fungal pathogens threatening global wheat cultivations. The rational utilization of leaf rust resistance (Lr) genes is still the most efficient method for the control of such diseases. The Lr47 gene introgressed from chromosome 7S of Aegilops speltoides still showed high resistance to the majority of Pt races collected in China. However, the Lr47 gene has not been cloned yet, and the regulatory network of the Lr47-mediated resistance has not been explored. In the present investigation, transcriptome analysis was applied on RNA samples from three different wheat lines (“Yecora Rojo”, “UC1037”, and “White Yecora”) carrying the Lr47 gene three days post-inoculation with the epidemic Pt race THTT. A comparison between Pt-inoculated and water-inoculated “Lr47-Yecora Rojo” lines revealed a total number of 863 upregulated (q-value < 0.05 and log2foldchange > 1) and 418 downregulated (q-value < 0.05 and log2foldchange < −1) genes. Specifically, differentially expressed genes (DEGs) located on chromosomes 7AS, 7BS, and 7DS were identified, ten of which encoded receptor-like kinases (RLKs). The expression patterns of these RLK genes were further determined by a time-scale qRT-PCR assay. Moreover, heatmaps for the expression profiles of pathogenesis-related (PR) genes and several transcription factor gene families were generated. Using a transcriptomic approach, we initially profiled the transcriptional changes associated with the Lr47-mediated resistance. The identified DEGs, particularly those genes encoding RLKs, might serve as valuable genetic resources for the improvement of wheat resistance to Pt.


Introduction
Puccinia triticina (Pt), the causal agent of wheat leaf rust, is a wide-spread fungal pathogen affecting global wheat production. Wheat leaf rust causes yield losses ranging from 5% to 20%, which can reach approximately 50% during epidemics [1,2]. Although yield losses from leaf rust are normally less damaging than those from stem rust and stripe rust, it causes high level of annual losses due to its more frequent occurrence in more world-wide regions [1,3]. In China, due to the climate change and

Wheat Lines Carrying the Lr47 Gene Showed High Resistance to the Epidemic Leaf Rust Race THTT
In one of our earlier studies, a wheat line introgressed with the Lr47 gene exerted high resistance to the majority of leaf rust races collected in China, including FHPR, THJS, THTS, FHJQ, PHSL, PHJT, PHSS, KHSS, PHTT, THTT, FHTR, FHHT, TGGT, FHTT, and FGMT [23]. The leaf rust race THTT was predominant in China and showed high virulence to most of the designated Lr genes [5]. In the current investigation, we collected wheat lines introgressed with the Lr47 gene in different spring common wheat cultivars "Yecora Rojo", "UC1037", and "White Yecora" from a previous research [23]. The "Lr47-Yecora Rojo" was a BC 6 F 5 line and thus mainly used in the subsequent experiments. Seedling plants of wheat lines "Lr47-Yecora Rojo-BC 6 F 5 ", "Lr47-UC1037", and "Lr47-White Yecora" were inoculated with water-suspended uredospores of the leaf rust race THTT. A susceptible wheat line "Thatcher" without any Lr genes served as a control. Considerable sporulation of Pt was observed in the susceptible line 8 dpi, whereas wheat lines carrying the Lr47 gene showed necrotic spots phenotype ( Figure 1).

3
(AABBDD, "Chinese Spring" TGACv1 version) to determine the DEGs regulated by the Lr47 gene during the wheat resistance to Pt infection. The genome-wide expression patterns of genes associated with the Lr47-mediated resistance was initially profiled, which would provide valuable genetic resources for the improvement of wheat resistance to Pt.

Wheat Lines Carrying the Lr47 Gene Showed High Resistance to the Epidemic Leaf Rust Race THTT
In one of our earlier studies, a wheat line introgressed with the Lr47 gene exerted high resistance to the majority of leaf rust races collected in China, including FHPR, THJS, THTS, FHJQ, PHSL, PHJT, PHSS, KHSS, PHTT, THTT, FHTR, FHHT, TGGT, FHTT, and FGMT [23]. The leaf rust race THTT was predominant in China and showed high virulence to most of the designated Lr genes [5]. In the current investigation, we collected wheat lines introgressed with the Lr47 gene in different spring common wheat cultivars "Yecora Rojo", "UC1037", and "White Yecora" from a previous research [23]. The "Lr47-Yecora Rojo" was a BC6F5 line and thus mainly used in the subsequent experiments. Seedling plants of wheat lines "Lr47-Yecora Rojo-BC6F5", "Lr47-UC1037", and "Lr47-White Yecora" were inoculated with water-suspended uredospores of the leaf rust race THTT. A susceptible wheat line "Thatcher" without any Lr genes served as a control. Considerable sporulation of Pt was observed in the susceptible line 8 dpi, whereas wheat lines carrying the Lr47 gene showed necrotic spots phenotype ( Figure 1).

Figure 1.
Wheat lines carrying the Lr47 gene showed high resistance to the prevalent leaf rust race THTT. The seedling plants of three wheat lines carrying the Lr47 gene in different background of "Yecora Rojo", "UC1037", and "White Yecora" were spray-inoculated with water-suspended uredospores of the prevalent leaf rust race THTT. A susceptible line "Thatcher" without any Lr genes served as a control. An abundant sporulation of leaf rust was observed in the susceptible line 8 dpi, whereas wheat lines carrying the Lr47 gene showed necrotic spots phenotype.

The Wheat Lines Carrying the Lr47 Gene were Subjected to RNA-Seq Analysis upon Pt Infection
The third leaves of the seedling plants of wheat lines "Lr47-Yecora Rojo-BC6F5", "Lr47-UC1037", and "Lr47-White Yecora" were inoculated with the water-suspended uredospores of leaf rust race THTT. Seedlings of "Lr47-Yecora Rojo-BC6F5" sprayed with water served as a mock control. RNA samples were harvested at 3 dpi from water-inoculated "Lr47-Yecora Rojo-BC6F5", as well as Pt-inoculated "Lr47-Yecora Rojo-BC6F5", "Lr47-UC1037", and "Lr47-White Figure 1. Wheat lines carrying the Lr47 gene showed high resistance to the prevalent leaf rust race THTT. The seedling plants of three wheat lines carrying the Lr47 gene in different background of "Yecora Rojo", "UC1037", and "White Yecora" were spray-inoculated with water-suspended uredospores of the prevalent leaf rust race THTT. A susceptible line "Thatcher" without any Lr genes served as a control. An abundant sporulation of leaf rust was observed in the susceptible line 8 dpi, whereas wheat lines carrying the Lr47 gene showed necrotic spots phenotype.

The Wheat Lines Carrying the Lr47 Gene were Subjected to RNA-Seq Analysis upon Pt Infection
The third leaves of the seedling plants of wheat lines "Lr47-Yecora Rojo-BC 6 F 5 ", "Lr47-UC1037", and "Lr47-White Yecora" were inoculated with the water-suspended uredospores of leaf rust race THTT. Seedlings of "Lr47-Yecora Rojo-BC 6 F 5 " sprayed with water served as a mock control. RNA samples were harvested at 3 dpi from water-inoculated "Lr47-Yecora Rojo-BC 6 F 5 ", as well as Pt-inoculated "Lr47-Yecora Rojo-BC 6 F 5 ", "Lr47-UC1037", and "Lr47-White Yecora". Three biological replicates were collected from each of the inoculated materials. A total number of 12 RNA samples were extracted and subjected to a 12-Gb RNA sequencing approach (Supplementary Table S1). An approximate number of 90-139 million 150-bp pair-end reads were sequenced from each sample and more than 92% of the reads were mapped on the reference genome of Triticum aestivum (AABBDD, "Chinese Spring" TGACv1 version, Supplementary Table S2). A total number of 128,794 genes were annotated in the transcriptome, with 14,366 "Novel" transcripts represented contigs that could not be found in the reference genome. For the overall gene expression abundance, we observed clear correlations (R 2 > 0.92) among biological replicates (Supplementary Figure S1). The expression levels of the transcripts in the database were estimated using the fragments per kilobase of transcript per million mapped reads (FPKM) values. The DEGs were identified by DESeq2. All raw data for the transcriptome assembly were deposited in NCBI under BioProject PRJNA498364.
Based on the GO annotations, the overall information for either upregulated or downregulated DEGs was summarized in three main categories of "molecular function", "cellular components", and "biological process" (Supplementary Figure S2). For the molecular function, the majority of both upregulated and downregulated DEGs were significantly associated with "binding" and "catalytic activity", whereas a portion of genes enriched in "response to stimulus" in the upregulated DEGs was higher than that in the downregulated ones. For the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, several DEGs, including Ca 2+ pathway genes (CNGCs, CDPK, Rboh, and CaM), FLS2, RPM1, PR1, and HSP90, were annotated in the "plant-pathogen interaction" pathway ( Figure 2 and Table 1). 5 Figure 2. The "plant-pathogen interaction" KEGG pathway annotations for the identified DEGs. Nodes of KEGG pathway enriched with DEGs were colored in green. Detailed information for the significant upregulated (frame in red color) or downregulated (frame in bright green color) DEGs was presented in Table 1.

Figure 3.
Distribution map of transcriptional-induced genes located on chromosomes 7AS, 7BS, and 7DS. The race-specific resistance gene Lr47 was introgressed from chromosome 7S of Aegilops speltoides to chromosome 7AS of hexaploid wheat Triticum aestivum. Previous genetic study identified a PCR marker Xabc465 that cosegregated with the Lr47 gene. A total number of 28 upregulated genes located on chromosomes 7AS, 7BS, and 7DS were identified from the assembled transcriptome, and ten of which encoded RLKs (labeled in red color, detailed information in Table 2). The physical position for each of the gene was determined and anchored to the distribution map using MapMaker software. . Distribution map of transcriptional-induced genes located on chromosomes 7AS, 7BS, and 7DS. The race-specific resistance gene Lr47 was introgressed from chromosome 7S of Aegilops speltoides to chromosome 7AS of hexaploid wheat Triticum aestivum. Previous genetic study identified a PCR marker Xabc465 that co-segregated with the Lr47 gene. A total number of 28 upregulated genes located on chromosomes 7AS, 7BS, and 7DS were identified from the assembled transcriptome, and ten of which encoded RLKs (labeled in red color, detailed information in Table 2). The physical position for each of the gene was determined and anchored to the distribution map using MapMaker software. Putative LRR receptor-like serine/threonine-protein kinase R47_CK: water-inoculated "Lr47-Yecora Rojo-BC 6 F 5 ", R47_PT: Pt-inoculated "Lr47-Yecora Rojo-BC 6 F 5 ", R47_1037: Pt-inoculated "Lr47-UC1037", R47_WY: Pt-inoculated "Lr47-White Yecora".
To validate the expression profiles of these chromosome 7AS/7BS/7DS-located upregulated RLK genes, a time-scale qRT-PCR assay was performed. Fully-expanded third leaves of the "Lr47-Yecora Rojo-BC 6 F 5 " line were spray-inoculated with uredospores of the leaf rust race THTT. RNA samples were harvested at 0, 2, 3, 5, and 8 dpi. Seedling plants of the "Lr47-Yecora Rojo-BC 6 F 5 " line sprayed with water served as a control. The expression level of a wheat gene encoding pathogenesis-related protein 1 (TaPR1, GenBank accession FJ815169.1) was significantly (*P < 0.05) upregulated at 2 dpi during the Lr47-mediated resistance (Figure 4), indicating that a regular plant defense response was triggered. The expression levels of genes encoding LRR-RLK_AA1807990, LRR-RLK_AA1807960, RPM1_AA1836940, WAK2_AA1839280, L-type RLK_AA1946750, and ZmPK1_AA1953360, were significantly (* P < 0.05) induced at 3 dpi. On the other hand, we did not observe any significant inductions of genes encoding LRR-RLK_AA1822740 and LRR-RLK_AA2043410.

1
To validate the expression profiles of these chromosome 7AS/7BS/7DS-located upregulated RLK genes, a time-scale qRT-PCR assay was performed. Fully-expanded third leaves of the "Lr47-Yecora Rojo-BC6F5" line were spray-inoculated with uredospores of the leaf rust race THTT. RNA samples were harvested at 0, 2, 3, 5, and 8 dpi. Seedling plants of the "Lr47-Yecora Rojo-BC6F5" line sprayed with water served as a control. The expression level of a wheat gene encoding pathogenesis-related protein 1 (TaPR1, GenBank accession FJ815169.1) was significantly (*P < 0.05) upregulated at 2 dpi during the Lr47-mediated resistance (Figure 4), indicating that a regular plant defense response was triggered. The expression levels of genes encoding LRR-RLK_AA1807990, LRR-RLK_AA1807960, RPM1_AA1836940, WAK2_AA1839280, L-type RLK_AA1946750, and ZmPK1_AA1953360, were significantly (* P < 0.05) induced at 3 dpi. On the other hand, we did not observe any significant inductions of genes encoding LRR-RLK_AA1822740 and LRR-RLK_AA2043410. . The expression profiles of selected genes encoding RLKs during the Lr47-mediated resistance were determined by qRT-PCR assay. Fully-expanded third leaves of "Lr47-Yecora Rojo BC6F5" line were spray-inoculated with leaf rust race THTT. Seedling plants of the same material sprayed with water served as a control. RNA samples were harvested at 0, 2, 3, 5, and 8 dpi. Five independent biological replicates were included. The transcript levels of all genes were expressed as linearized fold-TaActin levels using the 2 −ΔCt method. Wheat pathogenesis-related TaPR1 gene was employed as a positive control. Mean and standard error of relative expressions were calculated, and two-sample t-test (*P < 0.05) was conducted using Microsoft Excel software.

The Expression Patterns of Genes Encoding PR Proteins and Several Transcription Factor Families during the Lr47-Mediated Resistance were Profiled
To determine the possible downstream genes and key regulators of the Lr47-mediated resistance, the expression patterns of the genes encoding PR proteins and several transcription factor families were profiled using their FPKM values in the transcriptome database. Until the time of conducting our investigation, a total number of 18 gene families had been named as PR genes from various plant species in response to pathogen infections, but only a few of them had been cloned or Figure 4. The expression profiles of selected genes encoding RLKs during the Lr47-mediated resistance were determined by qRT-PCR assay. Fully-expanded third leaves of "Lr47-Yecora Rojo BC 6 F 5 " line were spray-inoculated with leaf rust race THTT. Seedling plants of the same material sprayed with water served as a control. RNA samples were harvested at 0, 2, 3, 5, and 8 dpi. Five independent biological replicates were included. The transcript levels of all genes were expressed as linearized fold-TaActin levels using the 2 −∆Ct method. Wheat pathogenesis-related TaPR1 gene was employed as a positive control. Mean and standard error of relative expressions were calculated, and two-sample t-test (*P < 0.05) was conducted using Microsoft Excel software.

The Expression Patterns of Genes Encoding PR Proteins and Several Transcription Factor Families during the Lr47-Mediated Resistance were Profiled
To determine the possible downstream genes and key regulators of the Lr47-mediated resistance, the expression patterns of the genes encoding PR proteins and several transcription factor families were profiled using their FPKM values in the transcriptome database. Until the time of conducting our investigation, a total number of 18 gene families had been named as PR genes from various plant species in response to pathogen infections, but only a few of them had been cloned or characterized in Triticeae crops of wheat and barley [24]. In the present investigation, we found that most of the PR genes excluding TaPR9 were significantly upregulated during the Lr47-mediated resistance at 3 dpi ( Figure 5). Wheat chemical-induced (WCI) genes and barley chemical-induced (BCI) genes have been reported to be sensitive to the treatment of benzothiadiazole (BTH, a salicylic acid homolog) [25,26]. Surprisingly, the majority of the WCI genes and homologs of BCI genes, excluding BCI4, BCI6, and BCI9, also showed significant inductions upon Pt infection in wheat lines carrying the Lr47 gene ( Figure 5).
2 characterized in Triticeae crops of wheat and barley [24]. In the present investigation, we found that most of the PR genes excluding TaPR9 were significantly upregulated during the Lr47-mediated resistance at 3 dpi ( Figure 5). Wheat chemical-induced (WCI) genes and barley chemical-induced (BCI) genes have been reported to be sensitive to the treatment of benzothiadiazole (BTH, a salicylic acid homolog) [25,26]. Surprisingly, the majority of the WCI genes and homologs of BCI genes, excluding BCI4, BCI6, and BCI9, also showed significant inductions upon Pt infection in wheat lines carrying the Lr47 gene ( Figure 5). Figure 5. The expression profiles of pathogenesis-related (PR) and BTH-sensitive genes during the Lr47-mediated resistance. FPKM values for each of the genes in "LR47_BC6F5_Pt", "LR47_1037_Pt", and "LR47_WY_Pt" were relative to that in "LR47_BC6F5_CK". A log2foldchange data transformation was conducted using Microsoft Excel software. A heatmap was generated by MeV software using the relative expression data of PR and BTH-sensitive (WCI and wheat homologs of BCI) genes. Genes with similar expression patterns were clustered using the "hierarchical clustering" function of the MeV software. LR47_BC6F5_CK: water-inoculated "Lr47-Yecora Rojo-BC6F5", LR47_BC6F5_Pt: Pt-inoculated "Lr47-Yecora Rojo-BC6F5", LR47_1037_Pt: Pt-inoculated "Lr47-UC1037", LR47_WY_Pt: Pt-inoculated "Lr47-White Yecora".
The expression profiles of several gene families encoding transcription factors, including WRKY, ERF, MYB, bHLH, bZIP, and NAC, were generated using their FPKM values, respectively. Interestingly, large portions of WRKY and ERF genes were significantly upregulated during the Lr47mediated resistance (Figure 6), whereas only a few members of MYB, bHLH, and bZIP genes showed constant inductions among wheat lines carrying the Lr47 gene in different genetic backgrounds (Supplementary Figure S4). Several polygenetic analyses were conducted using transcription factors encoded by the selected differentially expressed WRKY, ERF, MYB, bHLH, and bZIP genes and their closest homologs in other plant species, respectively (Supplementary Figures S5, S6, and S7). The expression profiles of pathogenesis-related (PR) and BTH-sensitive genes during the Lr47-mediated resistance. FPKM values for each of the genes in "LR47_BC6F5_Pt", "LR47_1037_Pt", and "LR47_WY_Pt" were relative to that in "LR47_BC6F5_CK". A log2foldchange data transformation was conducted using Microsoft Excel software. A heatmap was generated by MeV software using the relative expression data of PR and BTH-sensitive (WCI and wheat homologs of BCI) genes. Genes with similar expression patterns were clustered using the "hierarchical clustering" function of the MeV software. LR47_BC6F5_CK: water-inoculated "Lr47-Yecora Rojo-BC 6 F 5 ", LR47_BC6F5_Pt: Pt-inoculated "Lr47-Yecora Rojo-BC 6 F 5 ", LR47_1037_Pt: Pt-inoculated "Lr47-UC1037", LR47_WY_Pt: The expression profiles of several gene families encoding transcription factors, including WRKY, ERF, MYB, bHLH, bZIP, and NAC, were generated using their FPKM values, respectively. Interestingly, large portions of WRKY and ERF genes were significantly upregulated during the Lr47-mediated resistance (Figure 6), whereas only a few members of MYB, bHLH, and bZIP genes showed constant inductions among wheat lines carrying the Lr47 gene in different genetic backgrounds (Supplementary Figure S4). Several polygenetic analyses were conducted using transcription factors encoded by the selected differentially expressed WRKY, ERF, MYB, bHLH, and bZIP genes and their closest homologs in other plant species, respectively (Supplementary Figures S5-S7). FPKM values for each of the genes in "LR47_BC6F5_Pt", "LR47_1037_Pt", and "LR47_WY_Pt" were relative to that in "LR47_BC6F5_CK". A log2foldchange data transformation was conducted using Microsoft Excel software. Heatmaps were generated by MeV software using the relative expression data of the genes encoding the WRKY and ERF transcription factors, respectively. Genes with similar expression patterns were clustered using the "hierarchical clustering" function of the MeV software. Constantly induced genes were labeled with asterisk (*), and the deduced proteins of which were employed to generate polygenetic trees in Supplementary Figures S5 and S6, respectively. LR47_BC6F5_CK: water-inoculated "Lr47-Yecora Rojo-BC6F5", LR47_BC6F5_Pt: Pt-inoculated "Lr47-Yecora Rojo-BC6F5", LR47_1037_Pt: Pt-inoculated "Lr47-UC1037", LR47_WY_Pt: Pt-inoculated "Lr47-White Yecora".

Discussion
In the present investigation, the transcriptional changes associated with the Lr47-mediated wheat resistance to leaf rust were initially profiled. Approximately, only one-third of the designated Lr resistance loci are effective against the prevalent Pt race THTT in China [23]. The wheat isogenic line introgressed with the Lr47 gene still showed high resistance to the majority of leaf rust races collected in China [27]. The seedling resistance gene Lr47 was introgressed from chromosome 7S of Aegilops speltoides to chromosome 7AS of common wheat [19]. However, this gene has not been widely-used in wheat breeding programs in China. In the present investigation, wheat lines carrying Figure 6. The expression profiles of the genes encoding (A) WRKY and (B) ERF transcription factors during the Lr47-mediated resistance. FPKM values for each of the genes in "LR47_BC6F5_Pt", "LR47_1037_Pt", and "LR47_WY_Pt" were relative to that in "LR47_BC6F5_CK". A log2foldchange data transformation was conducted using Microsoft Excel software. Heatmaps were generated by MeV software using the relative expression data of the genes encoding the WRKY and ERF transcription factors, respectively. Genes with similar expression patterns were clustered using the "hierarchical clustering" function of the MeV software. Constantly induced genes were labeled with asterisk (*), and the deduced proteins of which were employed to generate polygenetic trees in Supplementary Figures S5 and S6, respectively. LR47_BC6F5_CK: water-inoculated "Lr47-Yecora Rojo-BC 6 F 5 ", LR47_BC6F5_Pt: Pt-inoculated "Lr47-Yecora Rojo-BC 6 F 5 ", LR47_1037_Pt: Pt-inoculated "Lr47-UC1037", LR47_WY_Pt: Pt-inoculated "Lr47-White Yecora".

Discussion
In the present investigation, the transcriptional changes associated with the Lr47-mediated wheat resistance to leaf rust were initially profiled. Approximately, only one-third of the designated Lr resistance loci are effective against the prevalent Pt race THTT in China [23]. The wheat isogenic line introgressed with the Lr47 gene still showed high resistance to the majority of leaf rust races collected in China [27]. The seedling resistance gene Lr47 was introgressed from chromosome 7S of Aegilops speltoides to chromosome 7AS of common wheat [19]. However, this gene has not been widely-used in wheat breeding programs in China. In the present investigation, wheat lines carrying the Lr47 gene in different genetic backgrounds ("Yecora Rojo", "UC1037", and "White Yecora") may serve as valuable germplasm resources (Figure 1). Despite their excellent resistance to leaf rust, wheat lines carrying the Lr47 gene had clear drawbacks in terms of grain yield (overall reduction by 3.8%, 220 kg/ha) and flour yield (21.8 g/kg), possibly due to the introgression of other genes co-segregated with Lr47 from Aegilops speltoides [22]. Consequently, gene cloning and mechanism exploration of the Lr47 gene will considerably facilitate the further practical application of this gene.
Due to the large (17 Gb) and polyploid (AABBDD) genome of common wheat, the identification of defense-related genes upon pathogen infection has never been an easy task. A total number of 863 upregulated and 418 downregulated DEGs were annotated in the transcriptome database. Clear activations of genes encoding RLKs were detected, some of which may be involved as key sensors or signal transducers during the HR conferred by the Lr47 gene. In an earlier investigation, a large amount of genes encoding RLKs were transcriptional induced upon infections of wheat stripe rust and powdery mildew [16]. Several genes encoding RLKs were upregulated in rice inoculated with non-host pathogen wheat leaf rust [28]. Interestingly, we also detected higher accumulation of genes encoding Cytochrome P450, metabolic products of which may be critically involved in the plant defense against insects and pathogens [29]. For the downregulated DEGs, we found that several heat-stress-responsive genes were significantly suppressed, including genes encoding heat-shock protein, heat-stress transcription factor, and the abscisic acid (ABA) receptor PYL4. The signaling hormone ABA was reported to be accumulated upon heat stress [30], whereas it negatively affected disease resistance by interfering with other phytohormones, such as salicylic acid, jasmonic acid, and ethylene [31]. Based on the suppression of multiple ABA-responsive genes established in our transcriptome database, we speculated that ABA exerted a negative role during the Lr47-mediated resistance.
For the enrichment of DEGs in the KEGG pathway of "plant-pathogen interaction", we found a clear involvement of the calcium (Ca 2+ ) signal pathway in the Lr47-mediated resistance ( Figure 2 and Table 1). Several genes in this pathway, including those encoding calcineurin B-like proteins (CBLs), CBL-interacting protein kinases (CIPKs), and cyclic nucleotide gated channels (CNGCs), were previously reported to participate in the wheat resistance to stripe rust [32,33].
Currently, all the cloned race-specific Lr genes (Lr1, Lr10, and Lr21) in wheat encoded NBS-LRR RLKs [7][8][9]. Approximately, 3400 full-length NBS-LRR RLK loci have been identified in common wheat, 1540 of which have been further confirmed as complete genes [34]. RLKs play crucial roles in both pathogen recognition and signal transduction during plant-pathogen interactions [35]. They may act through direct detection of the presence of pathogen effector molecules [36], but, more frequently, they exert their functions indirectly through sensing the host target modification by the pathogen effector [37].
In the current study, a total number of ten chromosome 7AS/7BS/7DS-located upregulated RLK genes were identified ( Table 2). The physical positions of such genes, as well as the linkage marker Xabc465 for Lr47 introgression, were predicted based on the blast results in the common wheat reference genome (Figure 3). Most of the identified RLK genes were significantly upregulated during the Lr47-mediated resistance in our subsequent qRT-PCR validation assay (Figure 4), the functions of which require further exploration.
We further extracted the expression data of several gene families that may be functioning during wheat resistance to leaf rust. For instance, a total number of 18 gene families were named as PR genes from various plant species in response to pathogen infections, but only a few of them were cloned or characterized in wheat [24]. Another group of genes, including WCI and BCI genes, was previously reported to be sensitive to the BTH treatment [25,26]. In the present investigation, most of the PR, WCI, and homologs of BCI genes were significantly upregulated upon Pt infection during the Lr47-mediated resistance ( Figure 5), indicating that a wide range of salicylic acid-related plant defenses were recruited by the Lr47 gene.
Transcription factors are key components of the gene regulatory networks in various biological processes. Several transcription factor gene families, including WRKY, ERF, MYB, bHLH, bZIP, and NAC, have been reported to play crucial roles during plant defense and stress responses [38][39][40][41]. In the present investigation, we detected more pronounced inductions in the WRKY and ERF transcription factors during the Lr47-mediated resistance ( Figure 6). The naming system for transcription factors in common wheat or its relative species was confused. Our polygenetic analyses for the selected transcription factors provided initial clues to explore the key nodes of the Lr47-mediated resistance (Supplementary Figures S5-S7).
Wheat lines carrying the Lr47 gene showed great potential in improving wheat resistance to Pt infection in China. In the current study, we initially profiled the regulatory network of the Lr47-mediated resistance. A total of ten chromosome 7AS/7BS/7DS-located upregulated RLK genes were identified, and their expression patterns during the Lr47-mediated resistance were further determined by qRT-PCR assay. Specific PR genes, as well as those differentially expressed transcription factor genes, might be valuable genetic resources for the improvement of wheat resistance to Pt. The transcriptome itself will greatly facilitate the cloning process of the Lr47 gene.

Plants Growth and Leaf Rust Inoculation
Wheat lines introgressed with the Lr47 gene in different genetic backgrounds of "Yecora Rojo", "UC1037", and "White Yecora" were derived from previous study [22] and collected from Prof. Jorge Dubcovsky's lab at University of California, Davis. The epidemic Pt race THTT was isolated from the field as described in one of our former studies [23] and preserved in our lab. Wheat seedlings were kept in the greenhouse and inoculated with leaf rust as described [41]. The third leaves of the seedlings were sprayed with water-suspended uredospores of Pt race THTT upon their full expansion. Seedlings sprayed with water served as a mock control. Phenotype of the leaf rust was photographed at 10 dpi. RNA samples for Illumina sequencing were harvested at 3 dpi. For qRT-PCR assay, inoculated leaves were sampled from the "Lr47-Yecora Rojo-BC 6 F 5 " line at 0, 2, 3, 5, and 8 dpi. We set four independent biological replicates for each time point. All samples were rapidly frozen in liquid nitrogen and then ground into a fine powder.

RNA-seq and Bioinformatics Analysis
Total RNA for each of the samples was extracted using an RNA extraction kit (Qiagen, Hilden, Germany) following the manufacturer's instructions. Library preparation and sequencing procedure were conducted by Novogene Co., Ltd. following the KAPA Library preparation and Illumina HiSeq protocols (Illumina, San Diego, CA, USA). The RNA sequencing was run on a HiSeq 1000 instrument. A quality analysis was conducted to remove adapters and overrepresented sequences using in-house perl scripts. Raw reads with adapters, N ratio higher than 10%, or ratio of poor base (Q phred ≤ 20) higher than 50%, were removed. Clean reads were mapped on the reference genome of Triticum aestivum (AABBDD, "Chinese Spring" TGACv1 version) [42] using TopHat 2.0.8 software with default parameter settings for 150-bp pair-end sequencing reads [43]. Multi-mapped reads and low mapping reads were identified using hisat2 software, and reads with MAPQ = 60 were designated as high qualified uniq reads. Transcripts were extracted using Cufflinks v2.1.1 software with default parameter. All transcripts were compared with gene models in the reference genome using Cuffcompare, and those assembled contigs that could not be found in the reference genome (class_code "u") were annotated as "Novel" transcripts. Then, the expression levels of genes were determined using HTSeq v0.9.1 software [44]. Briefly, raw_readcount from HTSeq was employed to calculate the expected number of fragments per kilobase of transcript sequence per Millions base pairs (FPKM) sequenced, which have considered the effect of sequencing depth and gene length for the reads count at the same time, and is currently the most commonly used method for estimating gene expression levels [44].
DEGs were isolated by filtering the expression levels of genes between Pt-inoculated "Lr47-Yecora Rojo-BC 6 F 5 " and water-inoculated "Lr47-Yecora Rojo-BC 6 F 5 " for an FDR-adjusted p value < 0.05 using the DESeq2 software [45]. The GOseq package software was employed to give each of the genes in the transcriptome database a GO annotation [46]. Summaries of GO annotation categories for the selected genes were generated using function of "GOLevel2 Counter" in TBtools software. Using FPKM values for the selected genes from the transcriptome database, several heatmaps were generated by MeV software. Gene clusters were divided based on the expression patterns using the "hierarchical clustering" function of the MeV software following the "Pearson correlation" distance metric and "average linkage clustering" linkage method. To characterize the complex biological behaviors of the Lr47-mediated resistance, KEGG annotations were utilized to analyze the regulatory pathways enriched with DEGs. The physical positions for DEGs located on chromosomes 7AS, 7BS, and 7DS were collected from the high qualified reference genome of chromosome 7 [47]. The distribution map was generated using MapMaker software version 3.0. Multiple sequence alignments were conducted using the MUSCLE method and neighbor-joining trees were generated through the MEGA software version 7.0. The confidence of nodes in each of the neighbor-joining tree were calculated using 1000 bootstrap cycles.

qRT-PCR Assay
RNAs from leaf samples were isolated using a Qiagen RNA extraction kit. A Takara Reverse Transcription Kit was employed for the synthesis of the first-strand cDNA. Equal amount of the total RNA for each of the sample was used to synthesize cDNA. qRT-PCR primers for each of the tested genes were designed (Supplementary Table S3). The wheat TaActin gene (GenBank accession AB181991.1) was employed as an inner reference gene [48]. The amplification efficiency was determined by a preliminary qRT-PCR assay using six 2-fold diluted cDNA samples (1:1, 1:2, 1:4, 1:8, 1:16, and 1:32). qRT-PCR reactions were performed using TransGen SYBR Green qPCR mix with a Bio-Rad CFX Manager instrument under the following conditions: 3 min at 95 • C; 40 cycles of 10 s at 95 • C, 10 s at 60 • C, and 10 s at 72 • C. Melting curves ranging from 60 • C to 94 • C were recorded by the instrument to evaluate the amplified product. The expression levels of target genes were relative to that of the reference TaActin gene in the same sample using the 2 -∆Ct method [49,50]. The relative expression represents the ratio between the initial number of molecules of the target gene and that of TaActin. Therefore, the Y scales of the qRT-PCR results are comparable across genes and experiments.