Xin Du1,2,3Xianbo Feng1,2,3Ruoxuan Li1,2,3Yanlong Jin1,4Lihui Shang1,2,3
Jixin Zhao1,2,3Changyou Wang1,2,3
Tingdong Li1,2,3Chunhuan Chen1,3Zengrong Tian1,2,3Pingchuan Deng1,2,3*
Wanquan Ji1,2,3*- 1College of Agronomy, Northwest A & F University, Yangling, Shaanxi, China
- 2State Key Laboratory of Crop Stress Biology for Arid Areas, Northwest A & F University, Yangling, Shaanxi, China
- 3Shaanxi Research Station of Crop Gene Resources and Germplasm Enhancement, Ministry of Agriculture, Yangling, Shaanxi, China
- 4Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing, China
Leymus mollis (Trin.) Pilg. (2n = 4x = 28, NsNsXmXm) potentially harbours useful genes that might contribute to the improvement of wheat. We describe M862 as a novel wheat-L. mollis alien disomic substitution line from a cross between wheat cv. 7182 and octoploid Tritileymus M47. Cytological observations indicate that M862 has a chromosome constitution of 2n = 42 = 21II. Two 4D chromosomes of wheat substituted by two L. mollis Ns chromosomes were observed, using the GISH and ND-FISH analyses. Molecular marker, 55K SNP array and wheat-P. huashanica liquid array (GenoBaits®WheatplusPh) analyses further indicate that the alien chromosomes are L. mollis 4Ns. Therefore, it was deduced that M862 was a wheat-L. mollis 4Ns(4D) alien disomic substitution line. There were also changes in chromosomes 1A, 1D, 2B and 5A detected by ND-FISH analysis. Transcriptome sequencing showed that the structural variation of 1D, 1A and 5A may have smaller impact on gene expression than that for 2B. In addition, a total of 16 markers derived from Lm#4Ns were developed from transcriptome sequences, and these proved to be highly effective for tracking the introduced chromosome. M862 showed reduced height, larger grains (weight and width), and was highly resistance to CYR32 and CYR34 stripe rust races at the seedling stage and mixed stripe rust races (CYR32, CYR33 and CYR34) at the adult stage. It was also resistance to Fusarium head blight (FHB). This alien disomic substitution line M862 may be exploited as an important genetic material in the domestication of stipe rust and FHB resistance wheat varieties.
1 Introduction
Wheat is the most extensive and largest cultivated cereal worldwide, providing a variety of nutrients and about 20% of the calories consumed by humans (Shewry and Hey, 2015). Wheat production is threatened by various biotic and abiotic stresses, whereas beneficial genetic materials are increasingly scarce. Wild relatives have numerous agronomically valuable genes that can be introduced into common wheat by distant hybridization (Du et al., 2013; Liu et al., 2019), serving as an important way to broaden the genetic base of modern wheat varieties (Sepsi et al., 2008; Wang et al., 2020c).
Stripe rust (YR) caused by Puccinia striiformis f. sp. tritici (Pst) is a severe economic disease resulting in several pandemics and significant production losses worldwide (Yaniv et al., 2015; Sha et al., 2016a). At present, 88% of wheat yield around the world is threatened by stripe rust, and surpasses one billion dollars per year (Beddow et al., 2015). Fusarium head blight (FHB), caused by Fusarium graminearum, is emerging as a widely distributed and devastating disease to global wheat production (Wang et al., 2020). Deoxynivalenol (DON) is the most common and toxic fusarium mycotoxin, as well as the crucial virulence factor of F. graminosa infection (Chen et al., 2019). The rapid evolution and spread of new virulent wheat stripe rust races (Zhan et al., 2015), as well as limited FHB resistance sources, promote characterization and deployment of wild relatives to overcome these threats (Wang et al., 2020b).
Leymus mollis (Trin.) Pilger (2n = 4 x = 28, NsNsXmXm), a perennial tetraploid species in Leymus mollis (Gramineae: Myaceae), is mainly distributed in the northern coastal areas of China, Russia (Far East), North Korea, Japan and North America (Fan et al., 2009). It possesses strong and short stems, large spikes, multiple spikelets, tolerance to drought, cold, alkaline and saline conditions, in addition to resistance to many fungal and bacterial diseases (stripe rust, powdery mildew, leaf rust and Fusarium head blight) (Fatih, 1983; Mujeeb-Kazi et al., 1983; Merker, 1992), serving as an important tertiary gene pool for wheat breeding (Kishii et al., 2003). Leymus genome composition is generally regarded as NsNsXmXm. It has been proved that Ns subgenome in Leymus originated from Psathrostachys (Zhang and Dvôrák, 1991; Zhang et al., 2006; Fan et al., 2009; Culumber et al., 2011; Sha et al., 2016a; Sha et al., 2017a). Meanwhile, the origin of the Xm subgenome is still controversial, which had been proposed to be donated by Agropyron (P), Eremopyrum (F) (Fan et al., 2009; Sha et al., 2017a) or Psathyrostachys (Ns) (Sha et al., 2017a) or vary among species (Sha et al., 2016a).
Since the 1960s, wheat and L. mollis hybridization had been applied successfully (Tsitsin, 1965). The first hybrid seeds were obtained by embryo rescue and colchicine treatment of the hybrid between common wheat 7182 and L. mollis in China (Chen et al., 1985). The derivatives of wheat and L. mollis have appeared successively. Octoploid Tritileymus (2n = 8x = 56, AABBDDNsNs or AABBDDXmXm) displayed excellent agronomic characteristics, such as large grain and panicle, many flowers, and resistance to various fungal diseases (Wang et al., 2000; Wang et al., 2012). A series of addition, substitution, and translocation lines were obtained by crossing durum wheat or common wheat with octoploid Tritileymus (2n = 56, AABBDDNsNs) (Bao et al., 2012; Zhao et al., 2013; Pang et al., 2014a; Pang et al., 2014b; Zhao et al., 2016; Zhao et al., 2019). Several wheat-L. mollis addition lines and substitution lines with resistance to stripe rust and powdery mildew using a cross between common wheat 7182 and partial amphiploid M47 (2n = 8x = 56, AABBDDNsNs) were created previously (Yang et al., 2014a; Yang et al., 2015; Yang et al., 2017; Zhang et al., 2017a; Yang et al., 2020). Furthermore, a few studies have reported that the useful trait on the resistance to FHB for L. mollis have been introgressed into cultivated wheat varieties (Chen et al., 2005; Qi et al., 2008; Zhao et al., 2019).
Molecular marker analysis is the most convenient method to identify and trace the alien genetic material in the common wheat background (Wang et al., 2020c). The main markers include simple sequence repeats (SSRs), expressed sequence tag (EST), and PCR-based Landmark Unique Gene (PLUG) (Du et al., 2013; Pang et al., 2014a; Yang et al., 2020). However, these markers have some disadvantages such as poor stability and repeatability, as well as low specificity. Currently, high-throughput sequencing technology has been widely applied in the development of wild relatives specific markers as their high-efficiency and low-cost, including SLAF-seq (Chen et al., 2013; Zhang et al., 2015; Liu et al., 2018; Song et al., 2020; Wang et al., 2020c), GBS (Li et al., 2019a) and transcriptome sequencing (Li et al., 2017). So far, only two studies have reported the development of markers tightly linked to L. mollis, including 14,530 SNP markers developed based on genome DArTseq (Edet et al., 2018) and 13 DNA markers developed for 2Ns chromosome by SLAF-seq (Feng et al., 2022).
It is far from enough to introduce the exogenous genetic material from L. mollis into the background of wheat. Therefore, researchers need to continuously create distant hybridization materials of wheat and L. mollis with excellent disease resistance or other beneficial traits. In this study, a novel wheat-L. mollis 4Ns(4D) alien disomic substitution line M862 was obtained from a cross between common wheat 7182 and octoploid Tritileymus M47. M862 was analyzed by cytological identification, genome in situ hybridization (GISH), fluorescence in situ hybridization (FISH), wheat 55K SNP array, wheat-P. huashanica liquid array and molecular marker (SSR, EST, PLUG) identification, as well as disease resistance and agronomic traits evaluation. The results demonstrated that the wheat-L. mollis 4Ns(4D) alien disomic substitution line showed high resistance to stripe rust and FHB. Meanwhile, we develop specific molecular markers of Lm#4Ns, which will improve the efficiency of wheat breeding.
2 Materials and methods
2.1 Plant materials
Common wheat line 7182 (2n = 6x = 42, AABBDD), Huixianghong (HXH), Sumai3 (SM3), Xiaoyan22 (XY22), Chinese Spring (CS), Aikang58 (AK58), wheat-related species of Leymus mollis (Trin.) Pilger (2n = 4x = 28, NsNsXmXm), Psathyrostachys huashanica Keng ex P.C.Kuo (2n = 2x = 14, NsNs), Thinopyrum ponticum (Podp.) Barkworth & DR Dewey (2n = 10x = 70, EeEeEbEbExExStStStSt or JJJJJJJsJsJsJs), Secale cereale L. (2n = 2x = 14, RR) and Hordeum vulgare L. (2n=2x=14, HH), the progeny of wheat-L. mollis M47 (2n = 8x = 56, AABBDDNsNs), alien disomic substitution lines of Lm#1Ns(1D) and Lm#4Ns(4D), disomic addition lines of Lm#2Ns, Lm#3Ns, Lm#5Ns, Lm#6Ns and Lm#7Ns were used in this study (Supplementary Table 1). Among them, Lm#4Ns(4D) line M862 was derived from the BC1F4 progeny of a 7182/M47//M47 cross. M47 was a partial amphidiploid derived from common wheat 7182 and L. mollis (Yang et al., 2014a).
2.2 Cytological identification
For mitosis analyses, root tips were collected from ten independent healthy plants grown in the field during late March over two growing seasons (2021 and 2022, five plants per year). The excised roots were immediately placed on ice for 24 h. Subsequently, the roots were rinsed thoroughly with water to remove soil, then fixed in 3:1 Carnoy’s fluid (3:1 ethanol: glacial acetic acid) at room temperature. After fixation for 48 h, the root tips were stained with 1% acetic acid magenta for 12 h and stored at 4°C until further use. After removing the root caps from well-stained root tips, 1 mm of the meristematic zones were immersed in a drop of 45% acetic acid on a clean slide for cytogenetic analysis.
For meiosis analyses, young spikes were excised from ten healthy plants grown in the field during early April over two growing seasons (2021 and 2022, five plants per year), and fixed in 6:3:1 Carnoy’s fluid (6:3:1 ethanol: chloroform: glacial acetic acid) for 48-72 h (room temperature). Then, the anthers were stored at 4°C) until further use. The fixed anthers were squashed on slides and stained with 1% aceto-carmine before examination. Homologous chromosomes pairing in pollen mother cells were observed by Olympus BX43 microscope (Tokyo, Japan).
2.3 GISH, sequential FISH-GISH and ND-FISH analysis
Root tip processing procedures were performed as previously described (Han et al., 2006). The genomic DNA of L. mollis and P. huashanica were purified and used as probes by labeling with fluorescein-12-dUTP and Texas Red-5-dCTP. GISH was used to detect the alien chromosome composition of M862 (Zhang et al., 2017a). Sequential FISH-GISH was performed to analyze the karyotypes of M862 and foreign chromosomes on the same slide. The probes of Oligo-pSc119.2 (6-FAM-5′, green), Oligo-pTa535 (Tamra-5′, red), and Oligo-D (6-FAM-5′, green) were synthesized by Invitrogen Biotechnology Co., Ltd (Shanghai) (Tang et al., 2018), which were further used to detect the chromosome composition of D subgenome in M862 as previously described (Tang et al., 2014). 4, 6-diamidino-2-phenylindole (DAPI) was used to counterstain chromosomes. The slides were examined with the Olympus BX63 fluorescence microscope equipped with a DP80 camera (Tokyo, Japan).
2.4 Homoeologous group identification
2.4.1 Wheat 55K SNP Array and Wheat-P. huashanica Liquid Array Analysis
Genomic DNA of M862, 7182, and L. mollis were hybridized to Wheat 55K SNP genotyping array by China Golden Marker Biotechnology Company (Beijing). There were 53007 microarray probe sites per chip, including diploid markers. Based on IWGSC-RefSeq-v1.0 (https://wheat-urgi.versailles.inra.fr/SeqRepository/), a total of 49,060 markers had specific physical location information, evenly covering the whole wheat genome. The average distance between the markers was 0.2~0.3Mb. A dish quality control (DQC) value ≥ 0.82 and call rate (CR) ≥ 90 were the thresholds used as the criteria for SNP filtering. Collected data were analyzed using an Excel V2019 and Origin V2019 (OriginLab, USA).
Additionally, a wheat-P. huashanica liquid array (GenoBaits®WheatplusPh) was further used to clarify the alien chromosomes homoeologous group of M862. Genomic DNA of M862 was sent to MolBreeding Biotechnology Company (Shijiazhuang, China) for DNA library preparation and probe hybridization, which was further sequenced on the Illumina HiSeq X. The sequencing data were first checked for quality control, then aligned to the reference genome using BWA (Li, 2012). The sequencing depth for each probe was calculated by bamdst v.1.0.6 (https://github.com/shiquan/bamdst), then plotted using R package.
2.4.2 Molecular makers analysis
122 SSR makers, 34 EST-STS markers (http://wheat.pw.usda.gov/SNP/new/pcr_primers.shtml) and 13 PLUG (Polymerase Chain Reaction-based landmark unique gene) markers (Hu et al., 2011) were used to further determine the homoeologous group of the alien chromosomes. All primers were synthesized by Beijing AuGCT DNA-SYN Biotechnology Co., Ltd. As previously described (Yang et al., 2014a), 8% non-denaturing polyacrylamide gel was used to separate the PCR products of EST-STS and PLUG markers. PCR products of PLUG makers were digested with HaeIII at 37°C for 3 h or TaqI at 65°C for 2 h (Zhang et al., 2017a).
2.5 Disease resistance assessment
The Pst races CYR32 and CYR34 were used to assess stripe rust resistance at the seedling stage, while the mixed races (including CYR32, CYR33, and CYR34) were used to identify stripe rust resistance at the adult stage. At the seedling stage, 10 plants per pot were inoculated with CYR32 and CYR34 Pst races using spraying method in three independent biological repeats. After inoculation, the seedlings were firstly placed in a humidor that was kept in the dark for 24 hours at 8°C. Then, the incubator parameters were adjusted to 14h in the daytime, 60% humidity, 6000Lx, and 8 h at night with 70% humidity. When the control material HXH was fully infected, we began to assess the disease resistance every 4 days until 21 days after inoculation (Wang et al., 2020c).
Adult stage stripe rust resistance trials were performed under field conditions over two growing seasons (2021 and 2022) using mixtures of Pst races (CYR32, CYR33 and CYR34). All materials were planted in two rows with 20 plants. At the jointing stage, the mixed Pst races and talcum powder were inoculated by shaking powder method in a ratio of 1:200 after spaying water onto plants (Wang et al., 2020c). When the disease of HXH was completely developed, the symptoms of the infected plants were recorded, again 7 days later. Immunity type (IT) was divided into 6 grades from 0-4 (IT 0, IT 0;, IT 1, IT 2, IT 3, and IT 4) as previously reported (Ma et al., 1995). In brief, IT 0 = immune, no visible uredia and necrosis on leaves; IT 0; = nearly immune, no uredia with hypersensitive flecks on leaves; IT 1 = very resistant, few small uredia clearly visible embedded in well-defined necrotic areas; IT 2 = moderately resistant, few small- to medium-sized uredia surrounded by clearly-defined but less extensive necrotic areas; IT 3 = moderately susceptible, a lot of medium-sized uredia, no necrosis, but with chlorosis on leaves; IT 4 = highest susceptible, a large number of large-sized uredia without necrosis on leaves.
M862 and its parents were used to assess FHB resistance over two growing seasons (2021 and 2022) in the field. F. graminearum strain PH1 was used to evaluate FHB resistance. SM3 and XY22 were used as susceptible and resistant controls for FHB, respectively. At the flowering stage, 10 spikes were selected from each material in the field. The 5th floret from the top of each selected spike was inoculated by injecting 10 μL micro-conidial suspension (250,000 spores ml-1). Afterward, transparent plastic bags sprayed with mist water were placed on the injection spikes and kept at a high temperature until 72 h after inoculation. Disease severity in the inoculated wheat spikes was assessed 21 days after inoculation (Guo et al., 2015). The infected spikelet rate (ISR, %), the severity level (SL, 0 to 4) and resistance evaluation (RE, I, R, MR, MS, S) of FHB were recorded based on the Agricultural Industry Standard of the People’s Republic of China (PRC, 2007). FHB Type II resistance was defined as suppressing the spread of the disease from the point of infection.
2.6 Agronomic traits assessment
M862 and its parents were evaluated for the following agronomic traits: plant height, tiller number, spike length, spikelets number, florets number per spikelet, kernel number per spike, 1000-kernel weight, kernel length, kernel width and awnedness. Ten plants were measured for each material over two growing seasons (2021 and 2022). The length and width of grain were performed on 200 to 250 grains per line using the SC-G automated seed testing system (Wseen Detection Technology Co., Ltd., Hangzhou, China) (repeated five times). Data were tested for statistical significance by using the IBM SPSS Statistics 23 software (IBM Corp., Armonk, NY, USA).
2.7 Transcriptome sequencing and data analysis
RNA was extracted from the fresh leaves (third leaf stage) of M862 and 7182 using TRIZOL reagent (Thermo Fisher Scientific Inc., Shanghai, China) according to the manufacturer’s instructions. RNA quality was assessed using Agient 2100 Bioanalyzer. Transcriptome sequencing was performed using an Illumina NovaSeq 6000 platform at Beijing Biomarker Technologies Corporation (Beijing, China), with three biological replicates for each material. Clean data were obtained by removing the reads containing adaptor, poly-N and low quality using FASTX-Toolkit (Schmieder and Edwards, 2011). High-quality clean data from each sample were separately aligned to the wheat reference genome (IWGSC RefSeq v1.1) by Hisat2 (Pertea et al., 2016) (Supplementary Figure 1). To construct the candidate unigene dataset for Lm#4Ns, the unmapped reads for all samples were used in de novo transcriptome assembly using Trinity software package (Grabherr et al., 2011) with default parameters. The high-quality clean data from each sample were aligned to wheat reference sequences and de novo assembly by Hisat2 (Pertea et al., 2016). Only uniquely mapped reads were used for gene-level quantification using StringTie (Pertea et al., 2016). Differentially expressed genes (DEGs) were calculated by DESeq2 (FPKM≥1, Log2|fold change|≥1, adjusted p ≤ 0.05). Unigenes expressed in M862 but not expressed in 7182 were retained as Lm#4Ns candidate genes. To further evaluate the Lm#4Ns gene, the sequences were used to search against the existing protein database of Swiss-Prot, Triticum urartu (A), Aegilops speltoides (B), Aegilops tauschii (D), Hordeum vulgare (H), Secale cereale (R), Thinopyrum elongatum (E), Aegilops longissima (S). GO enrichment analysis was performed by AgriGO v2.0 analysis tools (http://bioinfo.cau.edu.cn/agriGO/). A corrected p ≤ 0.01 was set as a threshold to determine the significant enrichments.
2.8 Lm#4Ns-specific marker development based on the RNA−seq data
The Lm#4Ns candidate unigenes that meet the following criteria were used to develop molecular marker: (I) unigene sequences that had the highest identity to the fourth group chromosome in at least one related species; (II) a sequence identity less than 50% by comparison with the wheat reference genome sequence (IWGSC-RefSeq-v1.1). (III) more than 50% by comparison with the P. huashanica Keng reference genome sequence (unpublished). Primers were designed in batches using the primer3_core program (Deng et al., 2016). The major parameters for primer designing were set as follows: 18-27 bp in primer length, 57-63°C in melting temperature, 40-60% in GC content and 200-700 bp in product size. The designed primers were searched against the genomes of wheat, barley and rye on WheatOmics 1.0 (http://wheatomics.sdau.edu.cn/), and sequence-free amplification was finally used to synthesize the primers at AuGCT (Beijing AuGCT Biotechnology Co. Ltd., Beijing, China). Markers were firstly screened in M862 and its parents, then verified in the common wheat (CS and AK58), wheat-related species (P. huashanica, Th. ponticum, S. cereale and H. vulgare), and wheat-L. mollis derived progenies (Lm#1Ns(1D), Lm#2Ns, Lm#3Ns, Lm#5Ns, Lm#6Ns and Lm#7Ns).
3 Results
3.1 Cytological characterization of M862
M862 was derived from the BC1F4 progeny of the cross between common wheat cultivar 7182 and wheat-L. mollis partial amphidiploid M47 (2n = 56). The mitotic observation of M862 showed that 156 out of 161 cells (96.89%) in 28 root tips from ten plants (2021 and 2022, five plants per year) had a chromosome number of 2n = 42 (Figure 1A). There were 2 and 3 cells with 41 and 40 chromosomes, respectively. The meiotic analysis of 117 pollen mother cells (PMCs) at metaphase I (MI) in two years showed that 110 (94.02%) PMCs had a chromosome configuration of 2n = 21 II (Figure 1B), and 7 PMCs had 2n = 20 II + 2 I. While no trivalents or quadrivalents were observed at MI, and no chromosomes lagged at anaphase I. 21’+21’ form was produced by equal separation of chromosomes at anaphase of meiosis II (Figure 1C). Cytological observations indicated that M862 was a highly stable line.
Figure 1 Cytological observations of wheat-Leymus mollis substitution line M862. (A) Chromosomes of root tip cell at mitotic metaphase, 2n = 42. (B) Chromosome characteristics of pollen mother cell at meiotic metaphase I, 2n = 21II. (C) Chromosome characteristics of pollen mother cell at meiotic anaphase I, 2n = 21’ + 21’. Scale bar = 10 μm.
3.2 GISH, sequential FISH-GISH and ND-FISH analyses of M862
GISH was employed to analyze the source of alien chromosomes for M862 at mitotic metaphase by labeled genomic DNA of L. mollis and P. huashanica. It was observed that two chromosomes had the green signal as evidence for L. mollis (Figure 2A). On the same slide, the two L. mollis chromosomes were labeled by the red signal of P. huashanica (Figure 2B). The remaining 40 chromosomes labeled with a blue signal were derived from wheat. These results indicated that the two alien chromosomes were derived from L. mollis Ns genome.
Figure 2 In situ hybridization identification of wheat-Leymus mollis substitution line M862. (A) GISH analysis of M862 using L. mollis genomic DNA as the probe (green). (B) Sequential GISH analysis of M862 using P. hushanica genomic DNA as the probe (red). (C) mc-FISH analysis of M862 using probes Oligo-pSc119.2 (green) and Oligo-pTa535 (red). (D) GISH analysis of M862 using L. mollis genomic DNA as the probe (green). (E) FISH analysis of M862 using Oligo-D probe (green), the white arrow referred 12 D chromosomes. (F) Sequential FISH analysis of M862 using probes Oligo-pSc119.2 (green) and Oligo-pTa535 (red). (G) FISH karyotype of parent 7182 using probes Oligo-pSc119.2 (green) and Oligo-pTa535 (red). (H) FISH karyotype of M862 using probes Oligo-pSc119.2 (green) and Oligo-pTa535 (red). (I) FISH signal comparison between CS, 7182, and M862 on chromosomes 1D, 2B and 6B. The alien chromosomes are indicated by yellow arrows. Chromosomes were counterstained with DAPI (blue). Scale bar = 10 μm.
Sequential FISH-GISH results showed that the 4D chromosomes were absent and two chromosomes of L. mollis had weak red signal labeled with Oligo-pTa535, while there were no hybridization signals on alien chromosomes when labeled with Oligo-pSc119.2 (Figures 2C, D). 12 D-genome chromosomes were labeled with the green signal in M862 using Oligo-D (Figure 2E). ND-FISH analysis further indicated that two D chromosomes were missed (Figure 2F). Therefore, two 4D chromosomes were replaced by two Ns chromosomes of L. mollis in M862.
Oligo-pSc119.2 (green) and Oligo-pTa535 (red) were used to construct the FISH karyotype of 7182 (Figure 2G) and M862 (Figure 2H). The chromosome structural variation of M862 was analyzed based on the 7182 FISH karyotype (Figures 2H, I). It was found that the end of 1DL and 2BL chromosomes were labeled with the green signal in M862, whereas the 2BS chromosome of both 7182 and M862 was labeled with the red signal. In addition, Oligo-pTa535 probe signal intensity differences could be seen in 1AS (much stronger) and 5AL (much weaker) in M862. These results suggest that the alien substitution of wheat 4D chromosomes may be accompanied by chromosome structures mutations in chromosomes 1A, 1D, 2B and 5A.
3.3 Analysis of foreign chromosome homoeologous group in M862
The chromosome composition of M862 was assessed by wheat 55K SNP array. A total of 51256, 51531 and 30076 SNP loci were detected in M862, 7182 and L. mollis, respectively (Supplementary Table 2). The miss rate (52.5%), homozygous rate (31.1%) and heterozygous rate (16.4%) of SNP in 4D chromosome were substantially obviously higher than those in other chromosomes (Figure 3A). Correspondingly, the maximum percentage of SNP loci shared between M862 and L. mollis was on 4D chromosome (36.8%), which also displayed the minimum percentage of SNP loci shared between M862 and 7182 (24.1%) (Figure 3B).
Figure 3 Molecular marker analysis of wheat-Leymus mollis substitution line M862. (A, B) Wheat 55K array analysis. (A) Obvious heterozygous and deletions on the 4D chromosome. (B) Obvious crossing point in terms of the position of the 4D chromosome. (C, D) Wheat-P. huashanica liquid array (GenoBaits®WheatplusPh) analysis. (C) The genomic sequence of the wheat parent line 7182 is the main distribution on chromosomes 1A-7D. (D) The sequences of M862 are rare distribution in 4D, but intensively distributed in 4Ns. (E) SSR, EST and PLUG markers analysis. M: D2000, 1: Wheat parent line 7182; 2. P. huashanica; 3: L. mollis; 4: Wheat-L. mollis partial amphidiploid line M47; 5: M862. The x-axis refers to the position of markers on wheat in (A) and (B). The red arrow represents the specific band of the 4Ns chromosome of L. mollis or P. huashanica.
GenoBaits®WheatplusPh, a liquid array contained 90000 capture probes derived from both wheat (10000 capture probes) and P. huashanica (80000 capture probes), displayed in a way that could be understood more intuitively. For 7182, the detecting signal was evenly distributed along 21 wheat chromosomes and showed low signal intensity along seven P. huashanica chromosomes (1Ns-7Ns) (Figure 3C). In contrast, detecting signals in M862 occurred quite rarely on wheat 4D chromosome, but were significantly enriched on the P. huashanica 4Ns chromosome (Figure 3D).
A total of 169 molecular markers in the fourth homoeologous group were used to verify the accuracy of chip results. We observed that 36 (21.30%) markers amplified the same special bands in M862 and L. mollis, including 15 (12.30%) SSR makers, 10 (29.41%) EST-STS makers and 11 (84.62%) PLUG makers. Among them, 14 (38.89%, 14/36) amplified the same bands among M862, L. mollis and P. huashanica (Figure 3E; Supplementary Table 3). The remaining 22 markers displayed different amplification products between L. mollis and P. huashanica, indicating that there were some differences between the Ns chromosomes of L. mollis and P. huashanica. Overall, the molecular marker and high-throughput assay further confirmed that M862 was a wheat-L. mollis 4Ns(4D) alien disomic substitution line.
3.4 Evaluation of disease resistance and agronomy traits for M862
M862 and its parents were inoculated with CYR32 and CYR34 Pst races to determine stripe rust resistance at seedling stage. Half a month after inoculation, the susceptible control HXH had become fully diseased. Three successive investigations of disease resistance showed that 7182 was susceptible to CYR32 and CYR34. Nevertheless, the immunity type of M862 to CYR32 and CYR34 was corresponding to 0 and 2, indicating immunity and middle resistance, respectively. Correspondingly, its parents (M47) displayed high-level resistance to CYR32 and immunity to CYR34 (Figure 4A; Supplementary Table 4).
Figure 4 Evaluation of disease resistance and agronomic traits. (A, B) Identification of Stripe rust resistance of M862 and its parents. (A) Seedling stage reactions to Pst race CYR32 (left) and CYR34 (right). (B) Adult stage reactions with mixture Pst race CYR32 and CYR33. (C) FHB reactions to F.g PH1 of M862 and its parents. (D-O) Analysis of agronomic traits of M862 and its parents (p < 0.05). (D) Adult plant. (E) Spikes. (F) Spikelets. (G) 10-seeds length and width. 1: HXH; 2: 7182 wheat; 3: M47 part. amp; 4: M862 4Ns(4D); 5: L. mollis. (H) Plant height. (I) Tillers. (J) Spike length. (K) Spikelets per spike. (L) Kernels per spike. (M) Thousand kernel weight. (N) Kernel length. (O) Kernel width.
At the adult stage, the mixed Pst races CYR32, CYR33 and CYR34 were used to assess M862 and its parents over two growing seasons (2021, 2022). The control HXH was highly susceptible to the mixed Pst races CYR32 and CYR33 (Figure 4B; Supplementary Table 4), displaying a large number of visible uredia and a lack of visible cell death. That of 7182 was moderate sensitivity. However, L. mollis and its derived materials (M47 and M862) showed high resistance, with very few or no uredia and prominent cell death. The immunity types of M47 and M862 were 0 and 1, respectively. These results demonstrated that M862 carried new stripe rust resistance genes from L. mollis and was effective at the seedling and the adult stages, indicating 4Ns chromosome of L. mollis may carry a broad spectrum of stripe rust resistance genes.
After two consecutive years (2021 and 2022) of FHB susceptibility identification, the results of all materials were relatively consistent. The control XY22 and SM3 were susceptible and resistant to FHB, respectively. The infected spikelet rate (ISR = 3.45%, 10.50%) and severity level (SL =1, 1.4) of M862 were significantly lower than the parent 7182, but it was similar to the parent M47 (Figure 4C; Supplementary Table 5). The results showed that M862 inherited resistance gene(s) of FHB from M47, and finally traced back to the resistance source of L. mollis.
The agronomic traits were evaluated for M862 and its parents 7182 and M47 (Figures 4D–O). For M862, the plant height and spike length were significantly shorter than its parents, whereas the tillers and number of kernels per spike were similar to its parents. M862 had more 1000-kernel weight and wider kernel width than its parents. The kernel length of M862 was similar to M47 but longer than 7182. M862 had ultrashort awn, which was different from M47 and 7182. In conclusion, M862 is a potential material for increasing yield, because of its excellent properties including increased 1000-grain weight and kernel width, and keeping other traits constant.
3.5 Impact of Lm#4Ns on wheat gene expression
To investigate the impact of the Lm#4Ns on gene expression, we compared wheat gene expression profiles between leaf (3-leaf seedling stage) of 7182 and M862 using transcriptome sequencing. After removing low-quality reads, a total of 1.94 × 108 paired reads were obtained, with an average of 10.19 G clean bases for each sample (Supplementary Table 6). On average, approximately 89.6% of reads were aligned to the Chinese Spring genome (IWGSC-RefSeq-v1.1), supporting the expression evidence for 37493 high-confidence wheat genes (FPKM ≥ 1 in at least one sample). Expression analysis revealed a total number of 6894 significant differentially expressed genes (DEGs) (adjusted p≤ 0.05; Log2|fold change| ≥ 1) (Supplementary Table 7). Among them, 3644 genes (52.86%) were down-regulated in M862, of which 1152 genes (31.61%) were distributed on chromosome 4D (Figure 5A; Supplementary Table 8). Compared to 7182, the substitution of 4D chromosome by Lm#4Ns in M862 would up-regulated the expression of 3250 genes, which was more highly distributed on chromosome 2D (635 up-regulated genes, 19.54%).
Figure 5 Characterization of Lm#4Ns candidate genes using transcriptome sequencing. (A) The chromosome distribution of DEGs between M862 and 7182. (B) Enrichment of up-regulated genes using Gene Ontology (GO). (C) The top 10 enriched GO terms for down-regulated genes. (D) Volcano plot indicating candidate genes for Lm#4Ns. (E) The distribution of Lm#4Ns candidate genes among each homologous chromosome group. (F) Venn-diagram of Lm#4Ns candidate genes have orthologs located on chromosomes of homoeologous group 4 among related species. T. urartu (A), A speltoides (B), A tauschii (D), H vulgare (H), S. cereale (R), Th. elongatum (E) and (A) longissima (S) genomes.
Gene ontology (GO) analysis showed that the 1531 up-regulated genes were enriched in molecular function mainly related to aspartic-type endopeptidase activity, and cation transmembrane transporter activity (Figure 5B; Supplementary Table 9). Meanwhile, the 3363 down-regulated genes were significantly enriched in nucleic acid binding, heterocyclic compound binding, organic cyclic compound binding, transition metal ion binding, cellular nitrogen compound biosynthetic process, and biosynthetic process (Figure 5C, Supplementary Table 10). Overall, these results indicated that the substitution of 4D chromosome by Lm#4Ns had a certain effect on wheat gene expression and was involved in a subset of biological processes.
3.6 Development and evaluation of Lm#4Ns chromosomes specific molecular markers
To create a reference transcriptome for Lm#4Ns, the unmapped reads from all six samples were then assembled using Trinity software (v2.8.4). We obtained a total of 208,181 trinity unigenes, with an average length of 817 bp (Supplementary Table 6). A total of 12,895 differentially expressed unigenes were found between M862 and 7182, including 7249 up-regulated genes and 5647 down-regulated genes (Figure 5D; Supplementary Table 11). Among them, 4117 unigenes (31.9%) could be annotated with the Swiss-Prot database (Supplementary Table 12). To investigate the homology between Lm#4Ns and wheat related species, all the 12,895 unigenes were searched against the reference genome sequence of T. urartu (A), A. speltoides (B), A. tauschii (D), H. vulgare (H), S. cereale (R), Th. elongatum (E) and A. longissima (S) using local BLAST. The distribution of mapping numbers ranged from 5948 to 7600 for all the comparisons (Figure 5E; Supplementary Table 13). Approximately 26.7% (3,443/12,895) of these unigenes had the highest identity to the ones on the fourth chromosome group in at least one related species (Figure 5F; Supplementary Table 13). A. longissima (S) (13.53%) had the highest number of hits, followed by E (12.77%), D (12.77%), B (12.12%), H (11.37%), R (10.14%) and A (9.55%).
Additionally, a total of 458 Lm#4Ns unigenes were selected following the criteria as described in the methods, which were used to develop potential suitable markers. 449 Lm#4Ns candidate markers were designed, with a successful primer designing potential of 98.03% (Supplementary Table 14). The validation of newly developed Lm#4Ns markers was performed using 50 randomly selected markers. Finally, sixteen (32%) markers had similar amplification products in both M862 and L. mollis, but were absent in 7182 (Figure 6; Supplementary Table 15). These markers were considered as candidate-specific markers for tracing Lm#4Ns chromosomes in wheat background.
Figure 6 Candidate special primers for chromosomes of Lm#4Ns were designed by RNA-seq DEGs. (A) Lm4-2. (B) Lm4-3. (C) Lm4-4. (D) Lm4-5. (E) Lm4-6. (F) Lm4-7. (G) Lm4-10. (H) Lm4-15. (I) Lm4-17. (J) Lm4-22. (K) Lm4-26. (L) Lm4-29. (M) Lm4-36. (N) Lm4-37.(O) Lm4-42. (P) Lm4-43. M: D2000, 1: 7182, 2: L. mollis, 3: M862.
To evaluate the specificity of obtained markers for Lm#4Ns chromosomes, these 16 markers were further amplified in common wheat (CS and AK58), wheat related species (P. huashanica, Th. ponticum, S. cereale, H. vulgare) and wheat-L. mollis derived lines (Lm#1Ns(1D), Lm#2Ns, Lm#3Ns, Lm#5Ns, Lm#6Ns and Lm#7Ns). We observed that seven markers (Lm4-2, Lm4-5, Lm4-6, Lm4-10, Lm4-15, Lm4-22, Lm4-43) were only amplified in L. mollis and its derivative lines (M862 and M47) (Figure 7A; Supplementary Table 15), which could be used to trace Lm#4Ns chromosome. The remaining 9 markers showed polymorphism between Lm#4Ns and other materials (Figure 7B; Supplementary Table 15). Among them, five markers displayed transferable to P. huashanica and Th. ponticum, of which 3 markers (Lm4-4, Lm4-7, Lm4-37) showed similar product size. The markers Lm4-29 and Lm4-42 amplified specific bands in P. huashanica. The other two markers (Lm4-3, Lm4-36) amplified the same bands among Th. ponticum, S. cereale and H. vulgare. In addition, two markers (Lm4-3 and Lm4-26) and marker Lm4-36 amplified specific bands for S. cereale and H. vulgare, respectively. Overall, these 16 newly developed markers were useful resources for tracing Lm#4Ns chromosomes and other Triticeae species during molecular breeding.
Figure 7 The amplification of Lm#4Ns candidate specific primers in wheat, wheat-related species, and different chromosomes of L. mollis. (A) The target bands were amplified on different species or on other chromosomes of L. mollis. (B) The target bands were amplified only on the Lm#4Ns. M: D2000; 1: CS; 2: AK58; 3: 7182; 4: L. mollis; 5: P. huashanica; 6: Th. ponticum; 7: S. cereale; 8: H vulgare; 9: M47; 10: Lm#1Ns(1D); 11: Lm#2Ns; 12: Lm#3Ns; 13: Lm#4Ns(4D); 14: Lm#5Ns; 15: Lm#6Ns; 16: Lm#7Ns.
4 Discussion
In situ hybridization is a visual and accurate technique that is widely used in wheat distant hybridization (Metzlaff et al., 1986; Rayburn and Gill, 1987; Yang et al., 2015; Zhang et al., 2017a; Li et al., 2021). Genomic DNA of L. mollis and P. huashanica were usually used as GISH probe to identify the introduced L. mollis chromosome in wheat background (Yang et al., 2014a; Zhang et al., 2017a; Feng et al., 2022). In this study, the combination of GISH and FISH analysis demonstrated that the alien chromosomes in M862 were obtained from the L. mollis Ns genome, supporting that M47 parental amphidiploid line contains chromosomes from only Ns genome (Yang et al., 2014a). Since no FISH karyotype has been published in L. mollis, we performed wheat 55K SNP arrays to further confirm the composition of M862 chromosomes. Meanwhile, it was the first time that the wheat-P. huashanica liquid array (GenoBaits®WheatplusPh) was used to trace the foreign chromosome of wheat-L. mollis derivates. GenoBaits®WheatplusPh, a genotyping by target sequencing (GBTS) platform through in-solution capture, is constituted by 90000 capture probes derived from both wheat (10000 capture probes) and P. huashanica (80000 capture probes). Finally, the accuracy of the liquid array was further verified by combining the molecular markers of the fourth homoeologous group of wheat. Results showed that the GenoBaits®WheatplusPh displayed in a way that could be understood more intuitively to characterize the chromosome composition of wheat-L. mollis derivates. By taking advantage of these high-throughput platforms, their utilization will make a great contribution to chromosome engineering breeding.
The introduction of alien chromosomes into wheat by distant hybridization has occurred homoeologous chromosomal recombination between the wheat chromosomes and those of alien chromosomes (Han et al., 2004; Wang et al., 2020c). For example, Zhong 1 had a chromosome number of 52 that included 10 alien chromosomes and 42 wheat chromosomes, of which one pair of its A subgenome chromosomes contained a terminal Th. intermedium segment on the short arm (Han et al., 2004). Our earlier study showed that the substitution of chromosome 1D by Th. ponticum 1Js in common wheat 7182 was accompanied by chromosome structure mutations in wheat chromosomes 5A and 6B (Wang et al., 2020c). In our study, a comparison of FISH karyotype between M862 and its wheat parent line 7182 indicates that chromosome structure mutations may have occurred in chromosomes 1A, 1D, 2B, and 5A. These observations suggest that the presence of L. mollis 4Ns chromosome may suppress the wheat homoeologous pairing control genes, and hence enhance limited homoeologous chromosomal recombination between wheat chromosomes and those of L. mollis. Further transcriptome analysis was performed to assess the effects of chromosomal recombination on gene expression. Ours results revealed that the expression of 18.39% wheat expressed genes were significantly altered in M862. As expected, wheat chromosome 4D held the highest number of DEGs and majority of them were down-regulated. However, we did not observe obvious DEGs distributed on 4A and 4B chromosomes, suggesting that Lm#4Ns may achieve the best compensatory effect for the absence 4D chromosome after multiple generations of continuous self-crossing. Additionally, the number of DEGs on chromosome 2B ranks the third, but few DEGs are distributed on chromosome 1A, 1D and 5A. By combining the altered karyotypes of 1A, 1D, 2B and 5A in M862, we infer that the structural variation of 1D, 1A and 5A may have smaller impact on gene expression than that for 2B.
Stripe rust and FHB are two major diseases threatening wheat production worldwide. The limited availability of resistance genes in cultivated wheat and the continuous evolution of new pathogen races make it necessary to constantly search for new resistance sources (Li et al., 2020). To date, a series of germplasms with resistance to stripe rust at the adult stage have been developed from the distant hybridization between wheat and L. mollis, including two partial amphidiploids (M47 and M42) (Yang et al., 2014a), two translocation lines (shannong0096 and WL24-4) (Bao et al., 2012; Li et al., 2015), four disomic addition liens (2Ns, 3Ns, 6Ns and 7Ns) (Li et al., 2016; Yang et al., 2017; Du et al., 2020; Yang et al., 2020), and a 7Ns(7D) alien disomic substitution line (Yang et al., 2015). For FHB, a few wheat-Leymus derivatives have been identified as resistance resource (Chen et al., 2005). For example, three alien introgressions sharing a common distal segment of 7Lr#1S designated a novel FHB resistance gene Fhb3 (Qi et al., 2008), a translocation line T3AS-Lr7S (Wang et al., 2009) and a double substitution line DM96 (Lm#2Ns+3Ns) (Zhao et al., 2019) were highly resistant to FHB. According to previous studies, 2Ns, 3Ns, 6Ns and 7Ns chromosomes of L. mollis have stripe rust resistance, while the terminal of 7Lr of Leymus with Fhb3 have FHB resistance. In this study, M862 showed high resistance to stripe rust at the seedling and adult stages, as well as high resistance to FHB. Further analysis demonstrated that the resistance genes originated from Lm#4Ns, indicating that M862 is a novel resistance source of stripe rust and FHB. Particularly, M862 displayed lower plant height, enhanced 1000-grain weight, wider grain width, and keeping other traits constant. These results suggest that M862 could be exploited as a useful intermediate material in wheat improvement and breeding.
It remains challenging to de novo assemble the transcriptome of most Triticeae species and further utilization of transcriptome sequencing in developing species-specific markers in wheat background (Li et al., 2019b). For example, 33 EST markers on the 6P chromosome of Agropyron cristatum (Dai et al., 2012), 75 EST-SSR genome-specific molecular markers of Thinopyrum intermedium (Cui et al., 2016), 25 Dasypyrum villosum 6V#4S−specific markers (Li et al., 2017), 134 Aegilops longissima specific markers (Wang et al., 2018) and 76 markers specific to each chromosome arm of D. villosum#4 (Li et al., 2019b) were developed using transcriptome data. In this study, we successfully screened out 16 markers that could be used to track Lm#4Ns chromosome. Among them, seven were Lm#4Ns specific markers, and the other nine could trace the genes of P. huashanica, Th. ponticum, S. cereale, and H. vulgare to varying degrees. These markers will provide an efficient way to track the 4Ns chromosome of L. mollis and accelerate the process of wheat breeding.
Data availability statement
The sequencing data from transcriptome and Wheat-P. huashanica Liquid Array analysis have been submitted to the National Genomics Data Center (NGDC) of the China National Center for Bioinformation (CNCB) study accession CRA007625. All data supporting the findings of this study are available within the paper and within its supplementary materials published online. The datasets and materials used during current study are also available by emailing the corresponding authors on reasonable request.
Author contributions
The experiment was designed by WJ and PD. The experimental materials were developed by JZ, CW, CC and ZT. PD, XD and YJ analyzed the transcriptome data. PD and XD developed the molecular markers. XD and RL analyzed foreign chromosome using molecular markers. XD, RL, XF, LS evaluated materials disease resistance and agronomy traits. XD and RL validated Lm#4Ns chromosomes specific molecular marker. The manuscript was drafted by XD, WJ, PD and TL. All authors read and approved the final manuscript.
Funding
This work was supported by the National Key Research and Development Program of China (2021YFD1200601-04) and Natural Science Basic Research Plan in Shaanxi Province of China (2022JQ-172).
Acknowledgments
We are grateful to Prof. Xiangqi Zhang for providing materials and Prof. Gourong Wei for providing stripe rust pathogens.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.1012939/full#supplementary-material
References
Bao, Y., Wang, J., He, F., Ma, H., Wang, H. (2012). Molecular cytogenetic identification of a wheat (Triticum aestivum)- American dune grass (Leymus mollis) translocation line resistant to stripe rust. Genet. Mol. Res. 11, 3198–3206. doi: 10.4238/2012.May.22.2
Beddow, J. M., Pardey, P. G., Chai, Y., Hurley, T. M., Kriticos, D. J., Braun, H. J., et al. (2015). Research investment implications of shifts in the global geography of wheat stripe rust. Nat. Plants 1, 15132. doi: 10.1038/nplants.2015.132
Chen, S. Y., Fu, J., Gao, L. Z. (1985). The hybridization between Triticum aestivum and Leymus mollis. Acta Botanica Boreali-Occidentalia Sin. 4, 260–266. doi: CNKI:SUN:DNYX.0.1985-04-002
Chen, S., Huang, Z., Dai, Y., Qin, S., Gao, Y., Zhang, L., et al. (2013). The development of 7E chromosome-specific molecular markers for Thinopyrum elongatum based on SLAF-seq technology. PLoS One 8, e65122. doi: 10.1371/journal.pone.0065122
Chen, Y., Kistler, H. C., Ma, Z. (2019). Fusarium graminearum trichothecene mycotoxins: Biosynthesis, regulation, and management. Annu. Rev. Phytopathol. 57, 15–39. doi: 10.1146/annurev-phyto-082718-100318
Chen, P., Liu, W., Yuan, J., Wang, X., Zhou, B., Wang, S., et al. (2005). Development and characterization of wheat- Leymus racemosus translocation lines with resistance to fusarium head blight. Theor. Appl. Genet. 111, 941–948. doi: 10.1007/s00122-005-0026-z
Cui, Y., Bao, Y. G., Wang, H. G., Li, X. F. (2016). Development of specific molecular markers for Thinopyrum intermedium using RNA-seq data. J. Triticeae Crops 36 (6), 699–707. doi: 10.7606/j.issn.1009-1041.2016.06.03
Culumber, C. M., Larson, S. R., Jensen, K. B., Jones, T. A. (2011). Genetic structure of Eurasian and north American Leymus (Triticeae) wildryes assessed by chloroplast DNA sequences and AFLP profiles. Plant Syst. Evol. 294, 207–225. doi: 10.1007/s00606-011-0455-x
Dai, C., Zhang, J. P., Wu, X. Y., Yang, X. M., Li, X. Q., Liu, W. H., et al. (2012). Development of EST markers specific to Agropyron cristatum chromosome 6P in common wheat background. Acta Agronomica Sin. 38, 1791–1801. doi: 10.3724/SP.J.1006.2012.01791
Deng, P., Wang, M., Feng, K., Cui, L., Tong, W., Song, W., et al. (2016). Genome-wide characterization of microsatellites in triticeae species: abundance, distribution and evolution. Sci. Rep. 6, 32224. doi: 10.1038/srep32224
Du, X., Du, S. S., Ma, B., Wang, Y. Z., Chen, C. H., Ji, W. Q., et al. (2020). Molecular cytogenetics identification of a wheat-Leymus mollis Lm#3Ns disomic addition line with resistance to stripe rust. J. Triticeae Crops 40, 1–11. doi: 10.7606/j.issn.1009-1041.2016.05.03
Du, W., Wang, J., Wang, L., Zhang, J., Chen, X., Zhao, J., et al. (2013). Development and characterization of a Psathyrostachys huashanica keng 7Ns chromosome addition line with leaf rust resistance. PloS One 8, e70879. doi: 10.1371/journal.pone.0070879
Edet, O. U., Gorafi, Y. S. A., Cho, S. W., Kishii, M., Tsujimoto, H. (2018). Novel molecular marker-assisted strategy for production of wheat-Leymus mollis chromosome addition lines. Sci. Rep. 8, 16117. doi: 10.1038/s41598-018-34545-x
Fan, X., Sha, L. N., Yang, R. W., Zhang, H. Q., Kang, H. Y., Ding, C. B., et al. (2009). Phylogeny and evolutionary history of Leymus (Triticeae; poaceae) based on a single-copy nuclear gene encoding plastid acetyl-CoA carboxylase. BMC Evol. Biol. 9, 247. doi: 10.1186/1471-2148-9-247
Fatih, A. (1983). Analysis of the breeding potential of wheat-Agropyron and wheat-Elymus derivatives, 1 agronomic and quality characteristics. Hereditas 98, 287–295. doi: 10.1111/j.1601-5223.1983.tb00607.x
Feng, X., Du, X., Wang, S., Deng, P., Wang, Y., Shang, L., et al. (2022). Identification and DNA marker development for a wheat-Leymus mollis 2Ns (2D) disomic chromosome substitution. Int. J. Mol. Sci. 23 (5), 2676. doi: 10.3390/ijms23052676
Grabherr, M. G., Haas, B. J., Yassour, M., Levin, J. Z., Thompson, D. A., Amit, I., et al. (2011). Full-length transcriptome assembly from RNA-seq data without a reference genome. Nat. Biotechnol. 29, 644–652. doi: 10.1038/nbt.1883
Guo, J., Zhang, X., Hou, Y., Cai, J., Shen, X., Zhou, T., et al. (2015). High-density mapping of the major FHB resistance gene Fhb7 derived from Thinopyrum ponticum and its pyramiding with Fhb1 by marker-assisted selection. Theor. Appl. Genet. 128, 2301–2316. doi: 10.1007/s00122-015-2586-x
Han, F. P., Lamb, J., Birchler, J. (2006). High frequency of centromere inactivation resulting in stable dicentric chromosomes of maize. Proc. Natl. Acad. Sci. 103 (2006), 3238–3243. doi: 10.1073/pnas.0509650103
Han, F., Liu, B., Fedak, G., Liu, Z. (2004). Genomic constitution and variation in five partial amphiploids of wheat–Thinopyrum intermedium as revealed by GISH, multicolor GISH and seed storage protein analysis. Theor. Appl. Genet. 109, 1070–1076. doi: 10.1007/s00122-004-1720-y
Hu, L. J., Li, G. R., Zeng, Z. X., Chang, Z. J., Liu, C., Yang, Z. J. (2011). Molecular characterization of a wheat-Thinopyrum ponticum partial amphiploid and its derived substitution line for resistance to stripe rust. J. Appl. Genet. 52, 279–285. doi: 10.1007/s13353-011-0038-0
Kishii, M., Wang, C., Tsujimoto, H. (2003). Characteristics and behaviour of the chromosomes of Leymus mollis and L. racemosus (Triticeae, poaceae) during mitosis and meiosis. Chromosome Res. 11, 741–748. doi: 10.1023/B:CHRO.0000005774.00726.71
Li, H. (2012). Exploring single-sample SNP and INDEL calling with whole-genome de novo assembly. Bioinformatics 28, 1838–1844. doi: 10.1093/bioinformatics/bts280
Li, J., Dundas, I., Dong, C., Li, G., Trethowan, R., Yang, Z., et al. (2020). Identification and characterization of a new stripe rust resistance gene Yr83 on rye chromosome 6R in wheat. Theor. Appl. Genet. 133, 1095–1107. doi: 10.1007/s00122-020-03534-y
Li, H., Fan, R., Fu, S., Wei, B., Xu, S., Feng, J., et al. (2015). Development of Triticum aestivum-Leymus mollis translocation lines and identification of resistance to stripe rust. J. Genet. Genomics 42, 129–132. doi: 10.1016/j.jgg.2014.11.008
Li, S., Lin, Z., Liu, C., Wang, K., Du, L., Ye, X. (2017). Development and comparative genomic mapping of Dasypyrum villosum 6V#4S-specific PCR markers using transcriptome data. Theor. Appl. Genet. 130, 2057–2068. doi: 10.1007/s00122-017-2942-0
Liu, L., Luo, Q., Teng, W., Li, B., Li, H., Li, Y., et al. (2018). Development of Thinopyrum ponticum-specific molecular markers and FISH probes based on SLAF-seq technology. Planta 247, 1099–1108. doi: 10.1007/s00425-018-2845-6
Liu, L., Yuan, C. Y., Wang, M. N., See, D. R., Zemetra, R. S., Chen, X. M. (2019). QTL analysis of durable stripe rust resistance in the north American winter wheat cultivar skiles. Theor. Appl. Genet. 132, 1677–1691. doi: 10.1007/s00122-019-03307-2
Li, M., Wang, Y., Liu, X., Li, X., Wang, H., Bao, Y. (2021). Molecular cytogenetic identification of a novel wheat-Thinopyrum ponticum 1JS (1B) substitution line resistant to powdery mildew and leaf rust. Front. Plant Sci. 12, 727734. doi: 10.3389/fpls.2021.727734
Li, S., Wang, J., Wang, K., Chen, J., Wang, K., Du, L., et al. (2019b). Development of PCR markers specific to Dasypyrum villosum genome based on transcriptome data and their application in breeding Triticum aestivum-D. villosum4 alien chromosome lines. BMC Genomics 20, 289. doi: 10.1186/s12864-019-5630-4
Li, X., Yang, X. F., Ma, B., Chen, Z. G., Wang, C. Y., Chen, C. H., et al. (2016). Molecular cytogenetics identification of a wheat-Leymus mollis derivative with resistance to stripe rust. J. Triticeae Crops 36, 556–563. doi: 10.1139/gen-2016-0151
Li, D., Zhang, J., Liu, H., Tan, B., Zhu, W., Xu, L., et al. (2019a). Characterization of a wheat-tetraploid Thinopyrum elongatum 1E(1D) substitution line K17-841-1 by cytological and phenotypic analysis and developed molecular markers. BMC Genomics 20, 963. doi: 10.1186/s12864-019-6359-9
Ma, H., Singh, R. P., Mujeeb-Kazi, A. (1995). Suppression/expression of resistance to stripe rust in synthetic hexaploid wheat (Triticum turgidum × T. tauschii). Euphytica 83, 87–93. doi: 10.1007/BF01678034
Merker, A. (1992). The triticeae in cereal breeding. Hereditas 116, 277–280. doi: 10.1111/j.1601-5223.1992.tb00836.x
Metzlaff, M., Troebner, W., Baldauf, F., Schlegel, R., Cullum, J. (1986). Wheat specific repetitive DNA sequences- construction and characterization of four different genomic clones. Theor. Appl. Genet. 72, 207–210. doi: 10.1007/BF00266993
Mujeeb-Kazi, A., Bernard, M., Bekele, G. T., Mirand, J. L. (1983). “Incorporation of alien genetic information from elymus giganteus into triticum aestivum,” in Proceedings of the 6th international wheat genetics symposium plant germplasm institute. Ed. Sakamoto, S. (Kyoto, Japan: Plant Germplasm Institute), 223–231.
Pang, Y., Chen, X., Zhao, J., Du, W., Cheng, X., Wu, J., et al. (2014a). Molecular cytogenetic characterization of a wheat-Leymus mollis 3D(3Ns) substitution line with resistance to leaf rust. J. Genet. Genomics 41, 205–214. doi: 10.1016/j.jgg.2013.11.008
Pang, Y. H., Zhao, J. X., Du, W. L., Li, Y. L., Wang, J., Wang, L. M., et al. (2014b). Cytogenetic and molecular identification of a wheat-Leymus mollis alien multiple substitution line from octoploid Tritileymus × Triticum durum. Genet. Mol. Res. 13, 3903–3913. doi: 10.4238/2014.May.23.1
Pertea, M., Kim, D., Pertea, G. M., Leek, J. T., Salzberg, S. L. (2016). Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and ballgown. Nat. Protoc. 11, 1650–1667. doi: 10.1038/nprot.2016.095
Qi, L. L., Pumphrey, M. O., Friebe, B., Chen, P. D., Gill, B. S. (2008). Molecular cytogenetic characterization of alien introgressions with gene Fhb3 for resistance to fusarium head blight disease of wheat. Theor. Appl. Genet. 117, 1155–1166. doi: 10.1007/s00122-008-0853-9
Rayburn, A., Gill, B. S. (1987). Use of repeated DNA sequences as cytological markers. Am. J. Bot. 74 (4), 574–580. doi: 10.1002/j.1537-2197.1987.tb08678.x
Schmieder, R., Edwards, R. (2011). Quality control and preprocessing of metagenomic datasets. Bioinformatics 27, 863–864. doi: 10.1093/bioinformatics/btr026
Sepsi, A., Molnar, I., Szalay, D., Molnar-Lang, M. (2008). Characterization of a leaf rust-resistant wheat-Thinopyrum ponticum partial amphiploid BE-1, using sequential multicolor GISH and FISH. Theor. Appl. Genet. 116, 825–834. doi: 10.1007/s00122-008-0716-4
Sha, L., Fan, X., Li, J., Liao, J., Zeng, J., Wang, Y., et al. (2017a). Contrasting evolutionary patterns of multiple loci uncover new aspects in the genome origin and evolutionary history of Leymus (Triticeae; poaceae). Mol. Phylogenet. Evol. 114, 175–188. doi: 10.1016/j.ympev.2017.05.015
Sha, L.-N., Fan, X., Zhang, H.-Q., Kang, H.-Y., Wang, Y., Wang, X.-L., et al. (2016a). Phylogeny and molecular evolution of the DMC1 gene in the polyploid genus Leymus (Triticeae: Poaceae) and its diploid relatives. J. Systematics Evol. 54, 250–263. doi: 10.1111/jse.12188
Shewry, P. R., Hey, S. J. (2015). The contribution of wheat to human diet and health. Food Energy Secur 4, 178–202. doi: 10.1002/fes3.64
Song, L., Zhao, H., Zhang, Z., Zhang, S., Liu, J., Zhang, W., et al. (2020). Molecular cytogenetic identification of wheat-Aegilops biuncialis 5M(b) disomic addition line with tenacious and black glumes. Int. J. Mol. Sci. 21 (11), 4053. doi: 10.3390/ijms21114053
Tang, S., Tang, Z., Qiu, L., Yang, Z., Li, G., Lang, T., et al. (2018). Developing new oligo probes to distinguish specific chromosomal segments and the a, b, d genomes of wheat (Triticum aestivum l.) using ND-FISH. Front. Plant Sci. 9, 1104. doi: 10.3389/fpls.2018.01104
Tang, Z., Yang, Z., Fu, S. (2014). Oligonucleotides replacing the roles of repetitive sequences pAs1, pSc119.2, pTa-535, pTa71, CCS1, and pAWRC.1 for FISH analysis. J. Appl. Genet. 55, 313–318. doi: 10.1007/s13353-014-0215-z
Tsitsin, N. (1965). Remote hybridization as a method of creating new species and varieties of plants. Euphytica 14, 326–330. doi: 10.1007/BF00149519
Wang, J., Chen, X.-H., Du, W.-L., Zhao, J.-X., Wu, J., Cheng, X.-N., et al. (2012). Morphological and molecular cytogenetic characterization of partial octoploid Tritileymus. Genet. Resour. Crop Evol. 60, 1453–1462. doi: 10.1007/s10722-012-9933-4
Wang, X. P., Fu, J., Zhang, X. Q., Jing, J. K., Wen, Y. X. (2000). Molecular cytogenetic study on genome constitutions of octoploid Tritileymus. Acta Botanica Boreali-Occidentalia Sin. 42 (6), 582–586.
Wang, K., Lin, Z., Wang, L., Wang, K., Shi, Q., Du, L., et al. (2018). Development of a set of PCR markers specific to Aegilops longissima chromosome arms and application in breeding a translocation line. Theor. Appl. Genet. 131, 13–25. doi: 10.1007/s00122-017-2982-5
Wang, H. W., Sun, S., Ge, W., Zhao, L., Hou, B., Wang, K. (2020b). Horizontal gene transfer of Fhb7 from fungus underlies fusarium head blight resistance in wheat. Science 368 (6493), eaba5435. doi: 10.1126/science.aba5435
Wang, M., Wu, L., Mei, Y., Zhao, Y., Ma, Z., Zhang, X., et al. (2020c). Host-induced gene silencing of multiple genes of Fusarium graminearum enhances resistance to fusarium head blight in wheat. Plant Biotechnol. J. 18, 2373–2375. doi: 10.1111/pbi.13401
Wang, L., Yuan, J., Bie, T., Zhou, B., Chen, P. (2009). Cytogenetic and molecular identification of three Triticum aestivum-Leymus racemosus translocation addition lines. J. Genet. Genomics 36, 379–385. doi: 10.1016/S1673-8527(08)60127-7
Yang, X., Li, X., Wang, C., Chen, C., Tian, Z., Ji, W. (2017). Isolation and molecular cytogenetic characterization of a wheat-Leymus mollis double monosomic addition line and its progenies with resistance to stripe rust. Genome 60, 1029–1036. doi: 10.1139/gen-2017-0078
Yang, X. F., Wang, C. Y., Chen, C. H., Tian, Z. R., Ji, W. Q. (2020). Development and characterization of a wheat-Leymus mollis Lm#7Ns disomic addition line with resistance to stripe rust. Cereal Res. Commun. 48, 467–476. doi: 10.1007/s42976-020-00058-5
Yang, X., Wang, C., Chen, C., Zhang, H., Tian, Z., Li, X., et al. (2014a). Chromosome constitution and origin analysis in three derivatives of Triticum aestivum-Leymus mollis by molecular cytogenetic identification. Genome 57, 583–591. doi: 10.1139/gen-2014-0161
Yang, X., Wang, C., Li, X., Chen, C., Tian, Z., Wang, Y., et al. (2015). Development and molecular cytogenetic identification of a novel wheat-Leymus mollis Lm#7Ns (7D) disomic substitution line with stripe rust resistance. PloS One 10, e0140227. doi: 10.1371/journal.pone.0140227
Yaniv, E., Raats, D., Ronin, Y., Korol, A. B., Grama, A., Bariana, H., et al. (2015). Evaluation of marker-assisted selection for the stripe rust resistance gene Yr15, introgressed from wild emmer wheat. Mol. Breed 35 (1), 43. doi: 10.1007/s11032-015-0238-0
Zhang, H. B., Dvôrák, J. (1991). The genome origin of tetraploid species of Leymus (Poaceae: Triticeae) inferred from variation in repeated nucleotide sequences. Am. J. Bot. 78 (7), 871–884. doi: 10.1002/j.1537-2197.1991.tb14491.x
Zhang, A., Li, W., Wang, C., Yang, X., Chen, C., Zhu, C., et al. (2017a). Molecular cytogenetics identification of a wheat-Leymus mollis double disomic addition line with stripe rust resistance. Genome 60, 375–383. doi: 10.1139/gen-2016-0151
Zhang, H. Q., Yang, R. W., Dou, Q. W., Tsujimoto, H., Zhou, Y. H. (2006). Genome constitutions of hystrix patula, h-duthiei ssp duthiei and h-duthiei ssp longearistata (Poaceae: Triticeae) revealed by meiotic pairing behavior and genomic in-situ hybridization. Chromosome Res. 14, 595–604. doi: 10.1007/s10577-006-1057-2
Zhang, Y., Zhang, J., Huang, L., Gao, A., Zhang, J., Yang, X., et al. (2015). A high-density genetic map for p genome of Agropyron gaertn. based on specific-locus amplified fragment sequencing (SLAF-seq). Planta 242, 1335–1347. doi: 10.1007/s00425-015-2372-7
Zhan, H., Zhang, X., Li, G., Pan, Z., Hu, J., Li, X., et al. (2015). Molecular characterization of a new wheat-Thinopyrum intermedium translocation line with resistance to powdery mildew and stripe rust. Int. J. Mol. Sci. 16, 2162–2173. doi: 10.3390/ijms16012162
Zhao, J. X., Du, W.-L., Wu, J., Cheng, X. N., Gao, Y., Pang, Y.-H., et al. (2013). Development and identification of a wheat-Leymus mollis multiple alien substitution line. Euphytica 190, 45–52. doi: 10.1007/s10681-012-0772-3
Zhao, J. X., Liu, Y., Cheng, X., Pang, Y., Li, J., Su, Z., et al. (2019). Development and identification of a dwarf wheat-Leymus mollis double substitution line with resistance to yellow rust and fusarium head blight. Crop J. 7, 516–526. doi: 10.1016/j.cj.2018.11.012
Keywords: Leymus mollis, alien disomic substitution line, stripe rust, FHB, transcriptome sequencing, liquid array
Citation: Du X, Feng X, Li R, Jin Y, Shang L, Zhao J, Wang C, Li T, Chen C, Tian Z, Deng P and Ji W (2022) Cytogenetic identification and molecular marker development of a novel wheat-Leymus mollis 4Ns(4D) alien disomic substitution line with resistance to stripe rust and Fusarium head blight. Front. Plant Sci. 13:1012939. doi: 10.3389/fpls.2022.1012939
Received: 06 August 2022; Accepted: 03 October 2022;
Published: 01 November 2022.
Edited by:
Marianna Rakszegi, Hungarian Academy of Sciences (MTA), HungaryReviewed by:
Mian Abdur Rehman Arif, Nuclear Institute for Agriculture and Biology, PakistanPeter Sharp, The University of Sydney, Australia
Klaudia Kruppa, ELKH Agricultural Institute, Centre for Agricultural Research, Hungary
Copyright © 2022 Du, Feng, Li, Jin, Shang, Zhao, Wang, Li, Chen, Tian, Deng and Ji. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wanquan Ji, aml3YW5xdWFuMjAwOEAxMjYuY29t; Pingchuan Deng, ZGVuZ3BpbmdjaHVhbkBud3N1YWYuZWR1LmNu