Development of Molecular Markers Linked to Powdery Mildew Resistance Gene Pm4b by Combining SNP Discovery from Transcriptome Sequencing Data with Bulked Segregant Analysis (BSR-Seq) in Wheat

Genome-wide association studies on resistance to powdery mildew in cultivated emmer wheat

Powdery mildew, caused by the fungal pathogen Blumeria graminis (DC.) E. O. Speer f. sp. tritici Em. Marchal (Bgt), is a constant threat to global wheat (Triticum aestivum L.) production. Although ∼100 powdery mildew (Pm) resistance genes and alleles have been identified in wheat and its relatives, more is needed to minimize Bgt's fast evolving virulence. In tetraploid wheat (Triticum turgidum L.), wild emmer wheat [T. turgidum ssp. dicoccoides (Körn. ex Asch. & Graebn.) Thell.] accessions from Israel have contributed many Pm resistance genes. However, the diverse genetic reservoirs of cultivated emmer wheat [T. turgidum ssp. dicoccum (Schrank ex Schübl.) Thell.] have not been fully exploited. In the present study, we evaluated a diverse panel of 174 cultivated emmer accessions for their reaction to Bgt isolate OKS(14)-B-3-1 and found that 66% of accessions, particularly those of Ethiopian (30.5%) and Indian (6.3%) origins, exhibited high resistance. To determine the genetic basis of Bgt resistance in the panel, genome-wide association studies were performed using 46,383 single nucleotide polymorphisms (SNPs) from genotype-by-sequencing and 4331 SNPs from the 9K SNP Infinium array. Twenty-five significant SNP markers were identified to be associated with Bgt resistance, of which 21 SNPs are likely novel loci, whereas four possibly represent emmer derived Pm4a, Pm5a, PmG16, and Pm64. Most novel loci exhibited minor effects, whereas three novel loci on chromosome arms 2AS, 3BS, and 5AL had major effect on the phenotypic variance. This study demonstrates cultivated emmer as a rich source of powdery mildew resistance, and the resistant accessions and novel loci found herein can be utilized in wheat breeding programs to enhance Bgt resistance in wheat.

Wheat powdery mildew resistance gene Pm13 encodes a mixed lineage kinase domain-like protein

Wheat powdery mildew is one of the most destructive diseases threatening global wheat production. The wild relatives of wheat constitute rich sources of diversity for powdery mildew resistance. Here, we report the map-based cloning of the powdery mildew resistance gene Pm13 from the wild wheat species Aegilops longissima. Pm13 encodes a mixed lineage kinase domain-like (MLKL) protein that contains an N-terminal-domain of MLKL (MLKL_NTD) domain in its N-terminus and a C-terminal serine/threonine kinase (STK) domain. The resistance function of Pm13 is validated by mutagenesis, gene silencing, transgenic assay, and allelic association analyses. The development of introgression lines with significantly reduced chromosome segments of Ae. longissima encompassing Pm13 enables widespread deployment of this gene into wheat cultivars. The cloning of Pm13 may provide valuable insights into the molecular mechanisms underlying Pm13-mediated powdery mildew resistance and highlight the important roles of kinase fusion proteins (KFPs) in wheat immunity.

Genomics empowering conservation action and improvement of celery in the face of climate change

MAIN CONCLUSION

Integration of genomic approaches like whole genome sequencing, functional genomics, evolutionary genomics, and CRISPR/Cas9-based genome editing has accelerated the improvement of crop plants including leafy vegetables like celery in the face of climate change. The anthropogenic climate change is a real peril to the existence of life forms on our planet, including human and plant life. Climate change is predicted to be a significant threat to biodiversity and food security in the coming decades and is rapidly transforming global farming systems. To avoid the ghastly future in the face of climate change, the elucidation of shifts in the geographical range of plant species, species adaptation, and evolution is necessary for plant scientists to develop climate-resilient strategies. In the post-genomics era, the increasing availability of genomic resources and integration of multifaceted genomics elements is empowering biodiversity conservation action, restoration efforts, and identification of genomic regions adaptive to climate change. Genomics has accelerated the true characterization of crop wild relatives, genomic variations, and the development of climate-resilient varieties to ensure food security for 10 billion people by 2050. In this review, we have summarized the applications of multifaceted genomic tools, like conservation genomics, whole genome sequencing, functional genomics, genome editing, pangenomics, in the conservation and adaptation of plant species with a focus on celery, an aromatic and medicinal Apiaceae vegetable. We focus on how conservation scientists can utilize genomics and genomic data in conservation and improvement.

Genetic Mapping of a Novel Gene PmAege7M from Aegilops geniculata Conferring Resistance to Wheat Powdery Mildew

Powdery mildew, caused by Blumeria graminis f. sp. tritici (Bgt), is one of the most damaging foliage diseases of wheat across the world. Aegilops geniculata Roth is a valuable gene resource for enhancing wheat resistance to powdery mildew. This study identified Ae. geniculata accession PI 487224 as immune and PI 487228 as susceptible to powdery mildew. Genetic analysis of the F1, F2, and F2:3 progeny derived from PI 487224 × PI 487228 showed that powdery mildew resistance in PI 487224 was controlled by two independent dominant genes located on two different nonhomologous chromosomes. By combing bulked segregant RNA-Seq, genetic linkage analysis of a single resistance gene segregation population, and marker analysis of a set of 14 wheat-Ae. geniculata chromosome addition lines, one of the resistance genes, temperately designated PmAege7M, was mapped to a 4.9-cM interval flanked by markers STS7-55926 and SNP7-45792/STS7-65911 on the long arm of chromosome 7 Mg of PI 487224, spanning 604.73 to 622.82 Mb on the 7D long arm based on the Ae. tauschii reference genome (Aet_v4.0). The map and closely linked markers of PmAege7M from Ae. geniculata in this study will facilitate the transfer of PmAege7M into common wheat and fine mapping of the gene.

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.

Genetic basis of resistance against powdery mildew in the wheat cultivar "Tabasco"

European winter wheat cultivar "Tabasco" was reported to have resistance to powdery mildew disease caused by Blumeria graminis f. sp. tritici (Bgt) in China. In previous studies, Tabasco was reported to have the resistance gene designated as Pm48 on the short arm of chromosome 5D when a mapping population was phenotyped with pathogen isolate Bgt19 collected in China and was genotyped with simple sequence repeat (SSR) markers. In this study, single-nucleotide polymorphism (SNP) chips were used to rapidly determine the resistance gene by mapping a new F2 population that was developed from Tabasco and a susceptible cultivar "Ningmaizi119" and inoculated with pathogen isolate NCF-D-1-1 that was collected in the USA. The segregation of resistance in the population was found to link with Pm2 which was identified in Tabasco. Therefore, it was concluded that the previously reported Pm48 on chromosome arm 5DS in Tabasco should be the Pm2 gene on the same chromosome. The Pm2 was also found in European cultivars "Mattis" and "Claire" but not in any of the accessions from diploid wheat Aegilops tauschii or modern cultivars such as "Gallagher," "Smith's Gold," and "OK Corral" being used in the Great Plains in the USA. A KASP marker was developed to track the resistance allele Pm2 in wheat breeding.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01402-3.

Resistance to Powdery Mildew Is Conferred by Different Genetic Loci at the Adult-Plant and Seedling Stages in Winter Wheat Line Tianmin 668

Winter wheat line Tianmin 668 was crossed with susceptible cultivar Jingshuang 16 to develop 216 recombinant inbred lines (RILs) for dissecting its adult-plant resistance (APR) and all-stage resistance (ASR) against powdery mildew. The RIL population was genotyped on a 16K genotyping by target sequencing single-nucleotide polymorphism array and phenotyped in six field trials and in the greenhouse. Three loci-QPmtj.caas-2BL, QPmtj.caas-2AS, and QPmtj.caas-5AL-conferring APR to powdery mildew were detected on chromosomes 2BL, 2AS, and 5AL, respectively, of Tianmin 668. The effect of resistance to powdery mildew for QPmtj.caas-2BL was greater than that of the other two loci. A Kompetitive allele-specific PCR marker specific for QPmtj.caas-2BL was developed and verified on 402 wheat cultivars or breeding lines. Results of virulence and avirulence patterns to 17 Blumeria graminis f. sp. tritici isolates, bulked segregant analysis-RNA-sequencing, and a genetic linkage mapping identified a resistance allele at locus Pm4 in Tianmin 668 based on the seedling phenotypes of the RIL population. The PCR-based DNA sequence alignment and cosegregation of the functional marker with the phenotypes of the RIL population demonstrated that Pm4d was responsible for the ASR to isolate Bgt1 in Tianmin 668. The dissection of genetic loci for APR and ASR may facilitate the application of Tianmin 668 in developing powdery mildew-resistant wheat cultivars.

Transfer of the high-temperature adult-plant stripe rust resistance gene Yr62 in four Chinese wheat cultivars

Wheat stripe rust is one of the diseases that seriously affect wheat production worldwide. Breeding resistant cultivars is an effective way to control this disease. The wheat stripe rust resistance gene Yr62 has high-temperature adult-plant resistance (HTAP). In this study, PI 660,060, a single Yr62 gene line, was crossed with four Chinese wheat cultivars, LunXuan987 (LX987), Bainongaikang58 (AK58), ZhengMai9023 (ZM9023), and HanMai6172 (H6172). F1 seeds of four cross combinations were planted and self-crossed to develop the advance generations in the field. The seeds of each cross were mixed harvested and about 2400 to 3000 seeds were sown in each generation for F1 to F4 to maintain the maximum possible genotypes. Forty-five lines were selected and evaluated for resistance to stripe rust and agronomic traits, including plant height, number of grains per spike, and tiller number, in F5 and F6. Then, 33 lines with good agronomic traits and high disease resistance were developed to F9 generation. SSR markers Xgwm251 and Xgwm192 flank linked with the Yr62 were used to detect the presence of Yr62 in these 33 F9 lines. Of these, 22 lines were confirmed with the resistance gene Yr62. Finally, nine lines with good agronomic traits and disease resistance were successfully selected. The selected wheat lines in this study provide material support for the future breeding of wheat for stripe rust resistance.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01393-1.

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.

Pm SN15218 : A Potential New Powdery Mildew Resistance Gene on Wheat Chromosome 2AL

Powdery mildew, caused by Blumeria graminis f. sp. tritici (Bgt), is a devastating fungal disease that seriously damages the yield and quality of wheat in many regions of the world. Identifying new resistance genes and breeding new resistant varieties are effective methods to control this disease. The breeding line SN15218 shows good resistance against powdery mildew. We, therefore, developed an F2 population and 287 F2:3 families crossed between SN15218 and the powdery mildew susceptible cultivar Huixianhong (HXH). Genetic analysis indicated that a single dominant gene, designated herein Pm SN15218 , conferred resistance to the Bgt isolate E09 in SN15218. Bulked segregant RNA-Seq (BSR-Seq) analysis revealed that Pm SN15218 is located in a ∼25-Mb interval on chromosome 2AL. Using the polymorphism information between SN15218 and HXH, we developed 13 polymerase chain reaction (PCR) markers and mapped this gene to a 0.5-cM genetic interval between the two flanking markers PmM12 and PmM14, corresponding to a 6.01-Mb physical region in the Chinese Spring reference genome. The results of molecular marker analysis, allelic tests of resistance spectrum, and DNA resequencing indicated that Pm SN15218 is distinct from the known resistance gene Pm4b on 2AL.

Virulence Structure and Genetic Diversity of Blumeria graminis f. sp. avenae from Different Regions of Europe

The structure and dynamics of changes in pathogen populations can be analysed by assessing the level of virulence and genetic diversity. The aim of the present study was to determine the diversity of Blumeria graminis f. sp. avenae populations. Diversity and virulence of B. graminis f. sp. avenae was assessed based on 80 single-spore isolates collected in different European countries such as Poland (40 isolates), Germany (10), Finland (10), Czech Republic (10) and Ireland (10) using ISSR (Inter-Simple Sequence Repeats) and SCoT (Start Codon Targeted) markers. This work demonstrated differences in virulence of B. graminis f. sp. avenae isolates sampled from different countries. Molecular analysis showed that both systems were useful for assessing genetic diversity, but ISSR markers were superior and generated more polymorphic products, as well as higher PIC and RP values. UPMGA and PCoA divided the isolates into groups corresponding with their geographical origin. In conclusion, the low level of genetic differentiation of the analysed isolates has suggested that the evolution of B. graminis f. sp. Avenae population is slow, and thus the evolutionary potential of the pathogen is low. This work paves the way for future studies on B. graminis f. sp. Avenae population structure and dynamics based on genetic variability.

Pto Interaction Proteins: Critical Regulators in Plant Development and Stress Response

Pto interaction (Pti) proteins are a group of proteins that can be phosphorylated by serine/threonine protein kinase Pto, which have diverse functions in plant development and stress response. In this study, we analyzed the phylogenetic relationship, gene structure, and conserved motifs of Pti1s and predicted the potential cis-elements in the promoters of Pti1 genes using bioinformatics methods. Importantly, we systematically summarized the diverse functions of Pti1s in tomato, rice, Arabidopsis, potato, apple, and cucumber. The potential cis-elements in promoters of Pti1s decide their functional diversity in response to various biotic and abiotic stresses. The protein kinase Pti1 was phosphorylated by Pto and then modulated the downstream signaling pathways for PTI and ETI in the disease insistence process. In addition, some transcription factors have been defined as Ptis (e.g., Pti4, Pti5, and Pti6) originally, which actually were ethylene-response factors (ERFs). Pti4, Pti5, and Pti6 were modulated by salicylic acid (SA), jasmonate (JA), and ethylene signaling pathways and regulated diverse defense-related gene expression to cope with Pst infection and insect wounding.

Identification of a Pm4 Allele as a Powdery Mildew Resistance Gene in Wheat Line Xiaomaomai

Powdery mildew, caused by Blumeria graminis f. sp. tritici (Bgt), is one of the most destructive foliar diseases of wheat. In this study, we combined the bulked segregant RNA sequencing (BSR-seq) and comparative genomics analysis to localize the powdery mildew resistance gene in Chinese landrace Xiaomaomai. Genetic analysis of F₁ plants from a crossing of Xiaomaomai × Lumai23 and the derived F₂ population suggests that a single recessive gene, designated as pmXMM, confers the resistance in this germplasm. A genetic linkage map was constructed using the newly developed SNP markers and pmXMM was mapped to the distal end of chromosome 2AL. The two flanking markers 2AL15 and 2AL34 were closely linked to pmXMM at the genetic distance of 3.9 cM and 1.4 cM, respectively. Using the diagnostic primers of Pm4, we confirmed that Xiaomaomai carries a Pm4 allele and the gene function was further validated by the virus-induced gene silencing (VIGS). In addition, we systematically analyzed pmXMM in comparison with the other Pm4 alleles. The results suggest that pmXMM is identical to Pm4d and Pm4e at sequence level. Pm4b is also not different from Pm4c according to their genome/amino acid sequences. Only a few nucleotide variances were detected between pmXMM and Pm4a/b, which indicate the haplotype variation of the Pm4 gene.

Resistance of the Wheat Cultivar ‘Renan’ to Septoria Leaf Blotch Explained by a Combination of Strain Specific and Strain Non-Specific QTL Mapped on an Ultra-Dense Genetic Map

Quantitative resistance is considered more durable than qualitative resistance as it does not involve major resistance genes that can be easily overcome by pathogen populations, but rather a combination of genes with a lower individual effect. This durability means that quantitative resistance could be an interesting tool for breeding crops that would not systematically require phytosanitary products. Quantitative resistance has yet to reveal all of its intricacies. Here, we delve into the case of the wheat/Septoria tritici blotch (STB) pathosystem. Using a population resulting from a cross between French cultivar Renan, generally resistant to STB, and Chinese Spring, a cultivar susceptible to the disease, we built an ultra-dense genetic map that carries 148,820 single nucleotide polymorphism (SNP) markers. Phenotyping the interaction was done with two different Zymoseptoria tritici strains with contrasted pathogenicities on Renan. A linkage analysis led to the detection of three quantitative trait loci (QTL) related to resistance in Renan. These QTL, on chromosomes 7B, 1D, and 5D, present with an interesting diversity as that on 7B was detected with both fungal strains, while those on 1D and 5D were strain-specific. The resistance on 7B was located in the region of Stb8 and the resistance on 1D colocalized with Stb19. However, the resistance on 5D was new, so further designated Stb20q. Several wall-associated kinases (WAK), nucleotide-binding and leucine-rich repeats (NB-LRR) type, and kinase domain carrying genes were present in the QTL regions, and some of them were expressed during the infection. These results advocate for a role of Stb genes in quantitative resistance and for resistance in the wheat/STB pathosystem being as a whole quantitative and polygenic.

Identification of Resistant Germplasm and Detection of Genes for Resistance to Powdery Mildew and Leaf Rust from 2,978 Wheat Accessions

Powdery mildew and leaf rust, caused by Blumeria graminis f. sp. tritici and Puccinia triticina, respectively, are widespread diseases of wheat worldwide. The use of resistant cultivars is considered the most economical, environment-friendly, and effective method to control these diseases. In the present study, a collection of 2,978 wheat accessions consisting of 1,394 advanced breeding lines, 1,078 Chinese cultivars, 291 introduced cultivars, 132 lines containing alien chromosomes, and 83 landraces was tested for reactions to powdery mildew and leaf rust. The results indicated that 659 wheat accessions (22.1%) were highly resistant to a widely prevalent B. graminis f. sp. tritici isolate, E09, at the seedling stage, and 390 were consistently resistant to the mixture of B. graminis f. sp. tritici isolates at the adult plant stage. Meanwhile, 63 accessions (2.1%) were highly resistant to leaf rust at the adult plant stage, of which 54 were resistant to a predominant and highly virulent P. triticina race, THTT, at the seedling stage. Notably, 17 accessions were resistant to both powdery mildew and leaf rust. To detect known genes for resistance to powdery mildew and leaf rust, these accessions were tested with gene-specific or tightly linked markers for seven powdery mildew genes (Pm genes; Pm2, Pm4, Pm5, Pm6, Pm8, Pm21, and Pm24) and 10 Lr genes (Lr1, Lr9, Lr10, Lr19, Lr20, Lr24, Lr26, Lr34, Lr37, and Lr46). Of the 659 powdery mildew-resistant accessions, 328 might carry single Pm genes and 191 carry combined Pm genes. Pm2 was detected at the highest frequency of 59.6%, followed by Pm8, Pm6, Pm21, Pm4, and Pm5, whereas Pm24 was not detected. In addition, 139 accessions might contain unknown Pm genes different from those tested in this study. In the 63 accessions resistant to leaf rust, four leaf rust genes (Lr genes; Lr1, Lr10, Lr26, and Lr34) were detected in 41 accessions singly or in combination, whereas six genes (Lr9, Lr19, Lr20, Lr24, Lr37, and Lr46) were not detected. Twenty-two accessions might contain unknown Lr genes different from those tested in this study. This study not only provided important information for rationally distributing resistance genes in wheat breeding programs, but also identified resistant germplasm that might have novel genes to enrich the diversity of resistance sources.

Characterization of a new gene for resistance to wheat powdery mildew on chromosome 1RL of wild rye Secale sylvestre

PmSESY, a new wheat powdery mildew resistance gene was characterized and genetically mapped to the terminal region of chromosome 1RL of wild rye Secale sylvestre. The genus Secale is an important resource for wheat improvement. The Secale species are usually considered as non-adapted hosts of Blumeria graminis f. sp. tritici (Bgt) that causes wheat powdery mildew. However, as a wild species of cultivated rye, S. sylvestre is rarely studied. Here, we reported that 25 S. sylvestre accessions were susceptible to isolate BgtYZ01, whereas the other five confer effective resistance to all the tested isolates of Bgt. A population was then constructed by crossing the resistant accession SESY-01 with the susceptible accession SESY-11. Genetic analysis showed that the resistance in SESY-01 was controlled by a single dominant gene, temporarily designated as PmSESY. Subsequently, combining bulked segregant RNA-Seq (BSR-Seq) analysis with molecular analysis, PmSESY was mapped into a 1.88 cM genetic interval in the terminus of the long arm of 1R, which was closely flanked by markers Xss06 and Xss09 with genetic distances of 0.87 cM and 1.01 cM, respectively. Comparative mapping demonstrated that the corresponding physical region of the PmSESY locus was about 3.81 Mb in rye cv. Lo7 genome, where 30 disease resistance-related genes were annotated, including five NLR-type disease resistance genes, three kinase family protein genes, three leucine-rich repeat receptor-like protein kinase genes and so on. This study gives a new insight into S. sylvestre that shows divergence in response to Bgt and reports a new powdery mildew resistance gene that has potential to be used for resistance improvement in wheat.

Bulked Segregant RNA-Seq Provides Distinctive Expression Profile Against Powdery Mildew in the Wheat Genotype YD588

Wheat powdery mildew, caused by the fungal pathogen Blumeria graminis f. sp. tritici (Bgt), is a destructive disease leading to huge yield losses in production. Host resistance can greatly contribute to the control of the disease. To explore potential genes related to the powdery mildew (Pm) resistance, in this study, we used a resistant genotype YD588 to investigate the potential resistance components and profiled its expression in response to powdery mildew infection. Genetic analysis showed that a single dominant gene, tentatively designated PmYD588, conferred resistance to powdery mildew in YD588. Using bulked segregant RNA-Seq (BSR-Seq) and single nucleotide polymorphism (SNP) association analysis, two high-confidence candidate regions were detected in the chromosome arm 2B, spanning 453,752,054-506,356,791 and 584,117,809-664,221,850 bp, respectively. To confirm the candidate region, molecular markers were developed using the BSR-Seq data and mapped PmYD588 to an interval of 4.2 cM by using the markers YTU588-004 and YTU588-008. The physical position was subsequently locked into the interval of 647.1-656.0 Mb, which was different from those of Pm6, Pm33, Pm51, Pm52, Pm63, Pm64, PmQ, PmKN0816, MlZec1, and MlAB10 on the same chromosome arm in its position, suggesting that it is most likely a new Pm gene. To explore the potential regulatory genes of the R gene, 2,973 differentially expressed genes (DEGs) between the parents and bulks were analyzed using gene ontology (GO), clusters of orthologous group (COG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. Based on the data, we selected 23 potential regulated genes in the enriched pathway of plant-pathogen interaction and detected their temporal expression patterns using an additional set of wheat samples and time-course analysis postinoculation with Bgt. As a result, six disease-related genes showed distinctive expression profiles after Bgt invasion and can serve as key candidates for the dissection of resistance mechanisms and improvement of durable resistance to wheat powdery mildew.

Mapping Powdery Mildew Resistance Gene pmYBL on Chromosome 7B of Chinese Wheat (Triticum aestivum L.) Landrace Youbailan

Chinese wheat landrace Youbailan has excellent resistance to powdery mildew caused by Blumeria graminis f. sp. tritici. In the present study, genetic analysis indicated that a recessive gene, tentatively designated pmYBL, was responsible for the powdery mildew resistance of Youbailan. pmYBL was located in the 695-to-715-Mb genomic region of chromosome 7BL, with 19 gene-linked single-nucleotide polymorphism (SNP) markers. It was flanked by SNP1-12 and SNP1-2 with genetic distances of 0.6 and 1.8 centimorgans, respectively. The disease reaction patterns of Youbailan and four cultivars (lines) carrying the powdery mildew (Pm) genes located on chromosome arm 7BL indicated that pmYBL may be allelic or closely linked to these genes. All of the SNP markers linked to pmYBL were diagnostic, indicating that these markers will be useful for pyramiding pmYBL using marker-assisted selection.

Repressed Gene Expression of Photosynthetic Antenna Proteins Associated with Yellow Leaf Variation as Revealed by Bulked Segregant RNA-seq in Tea Plant Camellia sinensis

The young leaves and shoots of albino tea cultivars are usually characterized as having a yellow or pale color, high amino acid, and low catechin. Increasing attention has been paid to albino tea cultivars in recent years because their tea generally shows high umami and reduced astringency. However, the genetic mechanism of yellow-leaf variation in albino tea cultivar has not been elucidated clearly. In this study, bulked segregant RNA-seq (BSR-seq) was performed on bulked yellow- and green-leaf hybrid progenies from a leaf color variation population. A total of 359 and 1134 differentially expressed genes (DEGs) were identified in the yellow and green hybrid bulked groups (Y_f vs G_f) and parent plants (Y_p vs G_p), respectively. The significantly smaller number of DEGs in Y_f versus G_f than in Y_p versus G_p indicated that individual differences could be reduced within the same hybrid progeny. Analysis of Gene Ontology and Kyoto Encyclopedia of Genes and Genomes revealed that the photosynthetic antenna protein was most significantly enriched in either the bulked groups or their parents. Interaction was found among light-harvesting chlorophyll a/b -binding proteins (LHC), heat shock proteins (HSPs), and enzymes involved in cuticle formation. Combined with the transcriptomic expression profile, results showed that the repressed genes encoding LHC were closely linked to aberrant chloroplast development in yellow-leaf tea plants. Furthermore, the photoprotection and light stress response possessed by genes involved in HSP protein interaction and cuticle formation were discussed. The expression profile of DEGs was verified via quantitative real-time PCR analysis of the bulked samples and other F₁ individuals. In summary, using BSR-seq on a hybrid population eliminated certain disturbing effects of genetic background and individual discrepancy, thereby helping this study to intensively focus on the key genes controlling leaf color variation in yellow-leaf tea plants.

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.

Fine-Mapping of the Powdery Mildew Resistance Gene mlxbd in the Common Wheat Landrace Xiaobaidong

Powdery mildew, which is caused by Blumeria graminis f. sp. tritici (Bgt), is a disease of wheat worldwide. Xiaobaidong is a Chinese wheat landrace, which still maintains good resistance against powdery mildew. To obtain more genetic markers closely linked to the powdery mildew resistance gene mlxbd and narrow the candidate region for its isolation, new simple sequence repeats and cross intron-spanning markers were designed based on the genome sequence of Triticum aestivum cultivar Chinese Spring chromosome 7BL. The flanking markers 7BLSSR49 and WGGC5746 were found to be tightly linked to mlxbd at genetic distances of 0.4 cM and 0.3 cM, respectively. The resistance locus was mapped to a 63.40 kb and 0.29 Mb region of the Chinese Spring genome and Zavitan genome, respectively. The linked markers of mlxbd could be used as diagnostic markers for mlxbd. The linked molecular markers and delineated genomic region in the sequenced Chinese Spring genome will assist the future map-based cloning of mlxbd.

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.

Alternative splicing of a barley gene results in an excess-tillering and semi-dwarf mutant

An excess-tillering semi-dwarf gene Hvhtd was identified from an EMS-induced mutant in barley and alternative splicing results in excess-tillering semi-dwarf traits. Tillering and plant height are important traits determining plant architecture and grain production in cereal crops. This study identified an excess-tillering semi-dwarf mutant (htd) from an EMS-treated barley population. Genetic analysis of the F₁, F₂, and F_2:3 populations showed that a single recessive gene controlled the excess-tillering semi-dwarf in htd. Using BSR-Seq and gene mapping, the Hvhtd gene was delimited within a 1.8 Mb interval on chromosome 2HL. Alignment of the RNA-Seq data with the functional genes in the interval identified a gene HORVU2Hr1G098820 with alternative splicing between exon2 and exon3 in the mutant, due to a G to A single-nucleotide substitution at the exon and intron junction. An independent mutant with a similar phenotype confirmed the result, with alternative splicing between exon3 and exon4. In both cases, the alternative splicing resulted in a non-functional protein. And the gene HORVU2Hr1G098820 encodes a trypsin family protein and may be involved in the IAA signaling pathway and differs from the mechanism of Green Revolution genes in the gibberellic acid metabolic pathway.

A predicted NEDD8 conjugating enzyme gene identified as a Capsicum candidate Rf gene using bulk segregant RNA sequencing

Cytoplasmic male sterility (CMS) is an important tool for producing F₁ hybrids, which can exhibit heterosis. The companion system, restorer-of-fertility (Rf), is poorly understood at the molecular level and would be valuable in producing restorer lines for hybrid seed production. The identity of the Rf gene in Capsicum (pepper) is currently unclear. In this study, using bulked segregant RNA sequencing (BSR-seq), a strong candidate Rf gene, Capana06g002866, which is annotated as a NEDD8 conjugating enzyme E2, was identified. Capana06g002866 has an ORF of 555 bp in length encoding 184 amino acids; it can be cloned from F₁ plants from the hybridization of the CMS line 8A and restorer line R1 but is not found in CMS line 8A. With qRT-PCR validation, Capana06g002866 was found to be upregulated in restorer accessions compared to sterile accessions. The relative expression in flower buds increased with the developmental stage in F₁ plants, while the expression was very low in all flower bud stages of the CMS lines. These results provide new insights into the Rf gene in pepper and will be useful for other crops utilizing the CMS system.

Current achievements and future prospects in the genetic breeding of chrysanthemum: a review

Chrysanthemum (Chrysanthemum morifolium Ramat.) is a leading flower with applied value worldwide. Developing new chrysanthemum cultivars with novel characteristics such as new flower colors and shapes, plant architectures, flowering times, postharvest quality, and biotic and abiotic stress tolerance in a time- and cost-efficient manner is the ultimate goal for breeders. Various breeding strategies have been employed to improve the aforementioned traits, ranging from conventional techniques, including crossbreeding and mutation breeding, to a series of molecular breeding methods, including transgenic technology, genome editing, and marker-assisted selection (MAS). In addition, the recent extensive advances in high-throughput technologies, especially genomics, transcriptomics, proteomics, metabolomics, and microbiomics, which are collectively referred to as omics platforms, have led to the collection of substantial amounts of data. Integration of these omics data with phenotypic information will enable the identification of genes/pathways responsible for important traits. Several attempts have been made to use emerging molecular and omics methods with the aim of accelerating the breeding of chrysanthemum. However, applying the findings of such studies to practical chrysanthemum breeding remains a considerable challenge, primarily due to the high heterozygosity and polyploidy of the species. This review summarizes the recent achievements in conventional and modern molecular breeding methods and emerging omics technologies and discusses their future applications for improving the agronomic and horticultural characteristics of chrysanthemum.

Fine mapping of the wheat powdery mildew resistance gene Pm52 using comparative genomics analysis and the Chinese Spring reference genomic sequence

A high-resolution genetic linkage map was constructed using the comparative genomics analysis approach and the wheat reference genome, which placed wheat powdery mildew resistance gene Pm52 in a 0.21-cM genetic interval on chromosome arm 2BL. The gene Pm52 confers resistance to powdery mildew and has been previously mapped on chromosome arm 2BL in winter wheat cultivar Liangxing 99. Because of its effectiveness against the disease, this study was initiated to finely map Pm52 using the comparative genomics analysis approach and the wheat reference genomic sequence. Based on the EST sequences that were located in the chromosome region flanking Pm52, four EST-SSR markers were developed, and another nine SSR markers were developed using the comparative genomics technology. These thirteen markers were integrated into a genetic linkage map using an F_2:3 subpopulation of the Liangxing 99 × Zhongzuo 9504 cross. Pm52 was mapped within a 3.2-cM genetic interval in the subpopulation that corresponded to a ~40-Mb genomic interval on chromosome arm 2BL of the Chinese Spring reference genome. The Pm52-flanking markers Xicsl163 and Xicsl62 identified 344 recombinant individuals from 8820 F₂ plants. Nine SSR markers generated from the Chinese Spring genomic interval were incorporated into a high-resolution genetic linkage map, which placed Pm52 in a 0.21-cM genetic interval corresponding to 5.6-Mb genomic region. The constructed high-resolution genetic linkage map will facilitate the map-based cloning of Pm52 and its marker-assisted selection.

Rapid identification of a candidate nicosulfuron sensitivity gene (Nss) in maize (Zea mays L.) via combining bulked segregant analysis and RNA-seq

A candidate nicosulfuron sensitivity gene Nss was identified by combining bulked segregant analysis and RNA-seq. Multiple mutations of this gene were discovered in nicosulfuron-sensitive maize compared with the tolerant. It has been demonstrated that variabilities exist in maize response to nicosulfuron. Two nicosulfuron-sensitive inbred lines (HB39, HB41) and two tolerant inbred lines (HB05, HB09) were identified via greenhouse and field trials. Genetic analysis indicated that the sensitivity to nicosulfuron in maize was controlled by a single, recessive gene. To precisely and rapidly map the nicosulfuron sensitivity gene (Nss), two independent F₂ segregating populations, Population A (HB41 × HB09) and Population B (HB39 × HB05), were constructed. By applying bulked segregant RNA-Seq (BSR-Seq), the Nss gene was, respectively, mapped on the short arm of chromosome 5 (chr5: 1.1-15.3 Mb) and (chr5: 0.5-18.2 Mb) using two populations, with 14.2 Mb region in common. Further analysis revealed that there were 43 and 119 differentially expressed genes in the mapping intervals, with 18 genes in common. Gene annotation results showed that a cytochrome P450 gene (CYP81A9) appeared to be the candidate gene of Nss associated with nicosulfuron sensitivity in maize. Sequence analysis demonstrated that two common deletion mutations existed in the sensitive maize, which might lead to the nicosulfuron sensitivity in maize. The results will make valuable contributions to the understanding of molecular mechanism of herbicide sensitivity in maize.

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

The gene Pm61 that confers powdery mildew resistance has been previously identified on chromosome arm 4AL in Chinese wheat landrace Xuxusanyuehuang (XXSYH). To facilitate the use of Pm61 in breeding practices, the bulked segregant analysis-RNA-Seq (BSR-Seq) analysis, in combination with the information on the Chinese Spring reference genome sequence, was performed in the F_2:3 mapping population of XXSYH × Zhongzuo 9504. Two single nucleotide polymorphism (SNP), two Kompetitive Allele Specific PCR (KASP), and six simple sequence repeat (SSR) markers, together with previously identified polymorphic markers, saturated the genetic linkage map for Pm61, especially in the proximal side of the target gene that was short of gene-linked markers. In the newly established genetic linkage map, Pm61 was located in a 0.71 cM genetic interval and can be detected in a high throughput scale by the KASP markers Xicsk8 and Xicsk13 or by the standard PCR-based markers Xicscx497 and Xicsx538. The newly saturated genetic linkage map will be useful in molecular marker assisted-selection of Pm61 in breeding for disease resistant cultivar and in its map-based cloning.

Mapping of the male sterile mutant gene ftms in Brassica rapa L. ssp. pekinensis via BSR-Seq combined with whole-genome resequencing

A male sterile mutant was created by ⁶⁰Co γ-rays of microspores isolated from Chinese cabbage DH line 'FT'. A candidate gene for the male sterile trait was identified as Bra010198. Male sterility is used for hybrid seed production in Chinese cabbage. In this study, we derived a male sterile mutant (ftms) from Chinese cabbage DH line 'FT' by irradiating microspores with ⁶⁰Co γ-rays and realized the rapid trait transformation from male fertility to sterility for creating valuable breeding materials. Genetic analysis indicated that the male sterile trait is controlled by a single recessive nuclear gene, ftms. Microspore development in mutant ftms was aborted at the tetrad stage and associated with severely retarded degeneration and vacuolation of tapetum. Using BSR-seq analysis, the candidate region for ftms was mapped on chromosome A05. A large F₂ population was created, and the region was narrowed to approximately 1.7-Mb between markers Indel20 and Indel14 via linkage analysis. The recombination frequency was extremely suppressed because the region was located on the chromosome A05 centromere. Whole-genome resequencing of mutant ftms and wild-type 'FT' aligned only one nonsynonymous SNP to Bra010198; this gene is a homolog of Arabidopsis KNS4/UPEX1, which encodes a putative β-(1,3)-galactosyltransferase that controls pollen exine development. Comparative sequencing verified the SNP position on the fifth exon of Bra010198 in mutant ftms. Further genotyping revealed that the male sterile phenotype was fully co-segregated with this SNP. Quantitative real-time PCR indicated that Bra0101918 specifically expressed in stamen. The data presented herein suggested that Bra010198 is a strong candidate gene for ftms. Hence, we developed a male sterile line for potential application in breeding and expanded the knowledge about the molecular mechanism underlying male sterility in Chinese cabbage.

Molecular mapping of the Cf-10 gene by combining SNP/InDel-index and linkage analysis in tomato (Solanum lycopersicum)

BACKGROUND

Leaf mold, one of the major diseases of tomato caused by Cladosporium fulvum (C. fulvum), can dramatically reduce the yield and cause multimillion dollar losses annually worldwide. Mapping the resistance genes (R genes) of C. fulvum and devising MAS based strategies for breeding new cultivars is an effective approach to improve the resistance in tomato. Up to now, many C. fulvum genes or QTLs have been mapped using different genetic materials, but few studies focused on Cf-10 gene positioning.

RESULTS

In this study, we investigated the genetic rules for Cf-10 and used a novel combinatorial strategy to rapidly map the Cf-10 gene. Initially, the performance of F1, F2 and BC1F1 individuals after infection, demonstrated that the resistance against C. fulvum was controlled by a single dominant gene. Two pools of resistant and susceptible individuals from F2 population were investigated, using mapping by sequencing approach and Cf-10 was found to be localized to 3.35 Mb and 3.74 Mb on chromosome 1, employing SNP/InDel index methods, respectively. After accounting for overlapping regions, these two algorithms yielded a total length of 3.29 Mb, narrowing down the target region. We further developed five serviceable KASP markers for this region based on sequencing data and conducted local QTL mapping using individuals from the F2 population, except for mapping by sequencing as mentioned above. Finally Cf-10 gene was mapped spanning a region of 790 kb, where only one gene (Solyc01g007130.3) was annotated as probable receptor protein kinase TMK1 with a LRR motif, a common R gene characteristic. The RT-qPCR analysis further confirmed the localization and the relative expression of Solyc01g007130.3 in Ontario 792 and was found to be significantly higher than that in Moneymaker at 9 dpi and 12 dpi, respectively.

CONCLUSION

This study proposed a novel combinatorial strategy by combining SNP-index, InDel-index analyses and local QTL mapping using KASP genotyping approach to rapidly map genes responsible for specific traits and provided a robust base for cloning the Cf-10 gene. Furthermore, these analyses suggest that Solyc01g007130.3 is a potential candidate to be regarded as Cf-10 gene.

RNA Sequencing-Based Bulked Segregant Analysis Facilitates Efficient D-genome Marker Development for a Specific Chromosomal Region of Synthetic Hexaploid Wheat

Common wheat originated from interspecific hybridization between cultivated tetraploid wheat and its wild diploid relative Aegilops tauschii followed by amphidiploidization. This evolutionary process can be reproduced artificially, resulting in synthetic hexaploid wheat lines. Here we performed RNA sequencing (RNA-seq)-based bulked segregant analysis (BSA) using a bi-parental mapping population of two synthetic hexaploid wheat lines that shared identical A and B genomes but included with D-genomes of distinct origins. This analysis permitted identification of D-genome-specific polymorphisms around the Net2 gene, a causative locus to hybrid necrosis. The resulting single nucleotide polymorphisms (SNPs) were classified into homoeologous polymorphisms and D-genome allelic variations, based on the RNA-seq results of a parental tetraploid and two Ae. tauschii accessions. The difference in allele frequency at the D-genome-specific SNP sites between the contrasting bulks (ΔSNP-index) was higher on the target chromosome than on the other chromosomes. Several SNPs with the highest ΔSNP-indices were converted into molecular markers and assigned to the Net2 chromosomal region. These results indicated that RNA-seq-based BSA can be applied efficiently to a synthetic hexaploid wheat population to permit molecular marker development in a specific chromosomal region of the D genome.

Fine mapping of powdery mildew resistance gene Pm4e in bread wheat (Triticum aestivum L.)

Fine mapping of wheat powdery mildew-resistance gene Pm4e to a 0.19 cM interval with sequence-based markers provides the foundation for map-based cloning and marker-assisted selection with breeder-friendly markers. Powdery mildew caused by Blumeria graminis f. sp. tritici is a wheat foliar disease that poses a serious threat to global wheat production. Pm4 is a resistance gene locus that has played a key role in controlling this disease in wheat production and a few resistance alleles of this locus have been identified. We have previously mapped the Pm4e allele to a 6.7 cM interval on chromosome 2AL. In this study, Pm4e was delimited to a 0.19 cM interval flanked by Xwgrc763 and Xwgrc865, through employment of a larger segregating population, derived from the cross of resistant parent D29 with susceptible parent Yangmai 158 (Y158), and enrichment of the genetic interval with markers developed on Chinese Spring (C.S.) survey sequence. In this interval, Pm4e co-segregated with a few markers, some of which were either D29-dominant or Y158-dominant, implying great sequence variation in the interval between D29 and Y158. Most of these co-segregation markers could not differentiate the Pm4 alleles from each other. Survey of 55 wheat cultivars with four co-dominant markers showed that the Pm4e-co-segregating loci always co-exist. Annotation of the Pm4e interval-corresponding C.S. sequence revealed more than a dozen resistance gene analogs clustered in a 2.4 Mb region, although C.S. is susceptible to the Pm4e-avirulent isolate Bgt2. This study has established the foundation for map-based cloning of Pm4e. Moreover, some of the co-dominant markers developed in this study could help in marker-assisted transfer of Pm4e into elite cultivars.