Development of SNP, KASP, and SSR Markers by BSR-Seq Technology for Saturation of Genetic Linkage Map and Efficient Detection of Wheat Powdery Mildew Resistance Gene Pm61

Fine mapping of powdery mildew resistance gene PmXNM in a Chinese wheat landrace Xiaonanmai

KEY MESSAGE

Powdery mildew resistance gene PmXNM, originated from the Chinese wheat landrace Xiaonanmai, was delimited to a 300.7-kb interval enriched with resistance genes. Powdery mildew, caused by Blumeria graminis f. sp. tritici (Bgt), is a globally devastating disease threatening the yield and quality of wheat worldwide. The use of broad-spectrum disease resistance genes from wheat landraces is an effective strategy to prevent this pathogen. Chinese wheat landrace Xiaonanmai (XNM) was immune to 23 tested Bgt isolates at the seedling stage. The F1, F2, and F2:4 progenies derived from the cross between XNM and Chinese Spring (CS) were used in this study. Genetic analysis revealed that powdery mildew resistance in XNM was controlled by a single dominant gene, temporarily designated PmXNM. Bulked segregant analysis and molecular mapping delimited PmXNM to the distal terminal region of chromosome 4AL flanked by markers caps213923 and kasp511718. The region carrying the PmXNM locus was approximately 300.7 kb and contained nine high-confidence genes according to the reference genome sequence of CS. Five of these genes, annotated as disease resistance RPP13-like proteins 1, were clustered in the target region. Haplotype analysis using the candidate gene-specific markers indicated that the majority of 267 common wheat accessions (75.3%) exhibited extensive gene losses at the PmXNM locus, as confirmed by aligning the targeted genome sequences of CS with those of other sequenced wheat cultivars. Seven candidate gene-specific markers have proven effective for marker-assisted introgression of PmXNM into modern elite cultivars.

Fighting wheat powdery mildew: from genes to fields

KEY MESSAGE

Host resistance conferred by Pm genes provides an effective strategy to control powdery mildew. The study of Pm genes helps modern breeding develop toward more intelligent and customized. Powdery mildew of wheat is one of the most destructive diseases seriously threatening the crop yield and quality worldwide. The genetic research on powdery mildew (Pm) resistance has entered a new era. Many Pm genes from wheat and its wild and domesticated relatives have been mined and cloned. Meanwhile, modern breeding strategies based on high-throughput sequencing and genome editing are emerging and developing toward more intelligent and customized. This review highlights mining and cloning of Pm genes, molecular mechanism studies on the resistance and avirulence genes, and prospects for genomic-assisted breeding for powdery mildew resistance in wheat.

First Report on Development of Genome-Wide Microsatellite Markers for Stock (Matthiola incana L.)

Stock (Matthiola incana (L.) R. Br.) is a famous annual ornamental plant with important ornamental and economic value. The lack of DNA molecular markers has limited genetic analysis, genome evolution, and marker-assisted selective breeding studies of M. incana. Therefore, more DNA markers are needed to support the further elucidation of the biology and genetics of M. incana. In this study, a high-quality genome of M. incana was initially assembled and a set of effective SSR primers was developed at the whole-genome level using genome data. A total of 45,612 loci of SSRs were identified; the di-nucleotide motifs were the most abundant (77.35%). In total, 43,540 primer pairs were designed, of which 300 were randomly selected for PCR validation, and as the success rate for amplification. In addition, 22 polymorphic SSR markers were used to analyze the genetic diversity of 40 stock varieties. Clustering analysis showed that all varieties could be divided into two clusters with a genetic distance of 0.68, which were highly consistent with their flower shape (potted or cut type). Moreover, we have verified that these SSR markers are effective and transferable within the Brassicaceae family. In this study, potential SSR molecular markers were successfully developed for 40 M. incana varieties using whole genome analysis, providing an important genetic tool for theoretical and applied research on M. incana.

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

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

Cis-element amplified polymorphism (CEAP), a novel promoter- and gene-targeted molecular marker of plants

In this study, we selected eight cis-elements: AAAG, ACGTG, CCGA, ACTCAT, GGTCA, TATCC, TGAC and GATAA, which are closely related to plant growth and development, signal transduction and stress response. The CEAP primers were 18 nucleotides long and consisted of a central cis-element nucleotide core flanked by a filler sequence at the 5' end and di- or tri-nucleotides at the 3' end. A total of two hundred and twenty-four primers were developed, and the PCR procedure consisted of 5 cycles of low-temperature annealing and 35 subsequent cycles of annealing at 50°C. The PCR products are electrophoretically separated by 1.8-2.3% agarose. The polymorphism of the CEAP marker was amplified in eight mango (Mangifera indica L.) species. The results showed that the CEAP primers could amplify clear, repeatable bands in mango and combine at least four cis-elements from which a large number of bands were amplified and six highly polymorphic primers for each cis-element can reach an accurate clustering result. The results of CEAP marker assays compared with ISSR, CBDP and iPBS marker assays showed that CEAP marker was better than the other three markers in the number of fragment bands, H and I indexes. In addition, we also tested the CEAP markers in rice, tomato, potato, wax gourd, citrus and longan and the results showed that the CEAP marker assay could amplify clear polymorphic bands in different species. Our results indicate that the CEAP markers could be universally used in different species for genetic diversity analysis, relationship analysis, and marker-assisted selection for breeding.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-022-01212-5.

Mapping QTL associated with partial resistance to Aphanomyces root rot in pea (Pisum sativum L.) using a 13.2 K SNP array and SSR markers

KEY MESSAGE

A stable and major QTL, which mapped to an approximately 20.0 cM region on pea chromosome 4, was identified as the most consistent region conferring partial resistance to Aphanomyces euteiches. Aphanomyces root rot (ARR), caused by Aphanomyces euteiches Drechs., is a destructive soilborne disease of field pea (Pisum Sativum L.). No completely resistant pea germplasm is available, and current ARR management strategies rely on partial resistance and fungicidal seed treatments. In this study, an F8 recombinant inbred line population of 135 individuals from the cross 'Reward' (susceptible) × '00-2067' (tolerant) was evaluated for reaction to ARR under greenhouse conditions with the A. euteiches isolate Ae-MDCR1 and over 2 years in a field nursery in Morden, Manitoba. Root rot severity, foliar weight, plant vigor and height were used as estimates of tolerance to ARR. Genotyping was conducted with a 13.2 K single-nucleotide polymorphism (SNP) array and 222 simple sequence repeat (SSR) markers. Statistical analyses of the phenotypic data indicated significant (P < 0.001) genotypic effects and significant G × E interactions (P < 0.05) in all experiments. After filtering, 3050 (23.1%) of the SNP and 30 (13.5%) of the SSR markers were retained for linkage analysis, which distributed 2999 (2978 SNP + 21 SSR) of the markers onto nine linkage groups representing the seven chromosomes of pea. Mapping of quantitative trait loci (QTL) identified 8 major-effect (R2 > 20%), 13 moderate-effect (10% < R2 < 20%) effect and 6 minor-effect (R2 < 10%) QTL. A genomic region on chromosome 4, delimited by the SNP markers PsCam037549_22628_1642 and PsCam026054_14999_2864, was identified as the most consistent region responsible for partial resistance to A. euteiches isolate Ae-MDCR1. Other genomic regions important for resistance were of the order chromosome 5, 6 and 7.

Characterization of PmHHXM, a New Broad-Spectrum Powdery Mildew Resistance Gene in Chinese Wheat Landrace Honghuaxiaomai

Powdery mildew, caused by fungal pathogen Blumeria graminis f. sp. tritici, is an agronomically important and widespread wheat disease causing severe yield losses. Deployment of broad-spectrum disease resistance genes is the preferred strategy to prevent this pathogen. Chinese wheat landrace Honghuaxiaomai (HHXM) was resistant to all 23 tested B. graminis f. sp. tritici isolates at the seedling stage. The F₁, F₂, and F_2:3 progenies derived from the cross HHXM × Yangmai 158 were used in this study, and genetic analysis revealed that a single dominant gene, designated PmHHXM, conferred resistance to B. graminis f. sp. tritici isolate E09. Bulked segregant analysis and molecular mapping initially located PmHHXM to the distal region of chromosome 4AL. To fine map PmHHXM, we identified two critical recombinants from 592 F₂ plants and delimited PmHHXM to a 0.18-cM Xkasp475200 to Xhnu552 interval covering 1.77 Mb, in which a number of disease resistance-related gene clusters were annotated. Comparative mapping of this interval revealed a perturbed synteny among Triticeae species. This study reports the new powdery mildew resistance gene PmHHXM, which seems different from three known quantitative trait loci/genes identified on chromosome 4AL and has significant values for further genetic improvement. Analysis of the polymorphisms of 13 cosegregating markers between HHXM and 170 modern wheat cultivars indicates that Xhnu227 and Xsts478700 developed here are ideal for marker-assisted introgression of this locus in wheat breeding.

Genomic Approaches to Identify Molecular Bases of Crop Resistance to Diseases and to Develop Future Breeding Strategies

Plant diseases are responsible for substantial crop losses each year and affect food security and agricultural sustainability. The improvement of crop resistance to pathogens through breeding represents an environmentally sound method for managing disease and minimizing these losses. The challenge is to breed varieties with a stable and broad-spectrum resistance. Different approaches, from markers to recent genomic and 'post-genomic era' technologies, will be reviewed in order to contribute to a better understanding of the complexity of host-pathogen interactions and genes, including those with small phenotypic effects and mechanisms that underlie resistance. An efficient combination of these approaches is herein proposed as the basis to develop a successful breeding strategy to obtain resistant crop varieties that yield higher in increasing disease scenarios.

QTL Mapping and Validation for Kernel Area and Circumference in Common Wheat via High-Density SNP-Based Genotyping

As an important component, 1,000 kernel weight (TKW) plays a significant role in the formation of yield traits of wheat. Kernel size is significantly positively correlated to TKW. Although numerous loci for kernel size in wheat have been reported, our knowledge on loci for kernel area (KA) and kernel circumference (KC) remains limited. In the present study, a recombinant inbred lines (RIL) population containing 371 lines genotyped using the Wheat55K SNP array was used to map quantitative trait loci (QTLs) controlling the KA and KC in multiple environments. A total of 54 and 44 QTLs were mapped by using the biparental population or multienvironment trial module of the inclusive composite interval mapping method, respectively. Twenty-two QTLs were considered major QTLs. BLAST analysis showed that major and stable QTLs QKc.sau-6A.1 (23.12-31.64 cM on 6A) for KC and QKa.sau-6A.2 (66.00-66.57 cM on 6A) for KA were likely novel QTLs, which explained 22.25 and 20.34% of the phenotypic variation on average in the 3 year experiments, respectively. Two Kompetitive allele-specific PCR (KASP) markers, KASP-AX-109894590 and KASP-AX-109380327, were developed and tightly linked to QKc.sau-6A.1 and QKa.sau-6A.2, respectively, and the genetic effects of the different genotypes in the RIL population were successfully confirmed. Furthermore, in the interval where QKa.sau-6A.2 was located on Chinese Spring and T. Turgidum ssp. dicoccoides reference genomes, only 11 genes were found. In addition, digenic epistatic QTLs also showed a significant influence on KC and KA. Altogether, the results revealed the genetic basis of KA and KC and will be useful for the marker-assisted selection of lines with different kernel sizes, laying the foundation for the fine mapping and cloning of the gene(s) underlying the stable QTLs detected in this study.

High-Resolution Mapping of the Novel Early Leaf Senescence Gene Els2 in Common Wheat

Early leaf senescence negatively impacts the grain yield in wheat (Triticum aestivum L.). Induced mutants provide an important resource for mapping and cloning of genes for early leaf senescence. In our previous study, Els2, a single incomplete dominance gene, that caused early leaf senescence phenotype in the wheat mutant LF2099, had been mapped on the long arm of chromosome 2B. The objective of this study was to develop molecular markers tightly linked to the Els2 gene and construct a high-resolution map surrounding the Els2 gene. Three tightly linked single-nucleotide polymorphism (SNP) markers were obtained from the Illumina Wheat 90K iSelect SNP genotyping array and converted to Kompetitive allele-specific polymerase chain reaction (KASP) markers. To saturate the Els2 region, the Axiom® Wheat 660K SNP array was used to screen bulked extreme phenotype DNA pools, and 9 KASP markers were developed. For fine mapping of the Els2 gene, these KASP markers and previously identified polymorphic markers were analyzed in a large F₂ population of the LF2099 × Chinese Spring cross. The Els2 gene was located in a 0.24-cM genetic region flanked by the KASP markers AX-111643885 and AX-111128667, which corresponded to a physical interval of 1.61 Mb in the Chinese Spring chromosome 2BL containing 27 predicted genes with high confidence. The study laid a foundation for a map-based clone of the Els2 gene controlling the mutation phenotype and revealing the molecular regulatory mechanism of wheat leaf senescence.

Appraisal of wheat genomics for gene discovery and breeding applications: a special emphasis on advances in Asia

We discussed the most recent efforts in wheat functional genomics to discover new genes and their deployment in breeding with special emphasis on advances in Asian countries. Wheat research community is making significant progress to bridge genotype-to-phenotype gap and then applying this knowledge in genetic improvement. The advances in genomics and phenomics have intrigued wheat researchers in Asia to make best use of this knowledge in gene and trait discovery. These advancements include, but not limited to, map-based gene cloning, translational genomics, gene mapping, association genetics, gene editing and genomic selection. We reviewed more than 57 homeologous genes discovered underpinning important traits and multiple strategies used for their discovery. Further, the complementary advancements in wheat phenomics and analytical approaches to understand the genetics of wheat adaptability, resilience to climate extremes and resistance to pest and diseases were discussed. The challenge to build a gold standard reference genome sequence of bread wheat is now achieved and several de novo reference sequences from the cultivars representing different gene pools will be available soon. New pan-genome sequencing resources of wheat will strengthen the foundation required for accelerated gene discovery and provide more opportunities to practice the knowledge-based breeding.

Identification of a Recessive Gene PmQ Conferring Resistance to Powdery Mildew in Wheat Landrace Qingxinmai Using BSR-Seq Analysis

Wheat powdery mildew is caused by Blumeria graminis f. sp. tritici (Bgt), a biotrophic fungal species. It is very important to mine new powdery mildew (Pm) resistance genes for developing resistant wheat cultivars to reduce the deleterious effects of the disease. This study was carried out to characterize the Pm gene in Qingxinmai, a winter wheat landrace from Xinjiang, China. Qingxinmai is resistant to many Bgt isolates collected from different wheat fields in China. F₁, F₂, and F_2:3 generations of the cross between Qingxinmai and powdery mildew susceptible line 041133 were developed. It was confirmed that a single recessive gene, PmQ, conferred the seedling resistance to a Bgt isolate in Qingxinmai. Bulked segregant analysis-RNA-Seq (BSR-Seq) was performed on the bulked homozygous resistant and susceptible F_2:3 families, which detected 57 single nucleotide polymorphism (SNP) variants that were enriched in a 40 Mb genomic interval on chromosome arm 2BL. Based on the flanking sequences of the candidate SNPs extracted from the Chinese Spring reference genome, 485 simple sequence repeat (SSR) markers were designed. Six polymorphic SSR markers, together with nine markers that were anchored on chromosome arm 2BL, were used to construct a genetic linkage map for PmQ. This gene was placed in a 1.4 cM genetic interval between markers Xicsq405 and WGGBH913 corresponding to 4.9 Mb physical region in the Chinese Spring reference genome. PmQ differed from most of the other Pm genes identified on chromosome arm 2BL based on its position and/or origin. However, this gene and Pm63 from an Iranian common wheat landrace were located in a similar genomic region, so they may be allelic.