RNA-Seq bulked segregant analysis enables the identification of high-resolution genetic markers for breeding in hexaploid wheat

Autoactive CNGC15 enhances root endosymbiosis in legume and wheat

Nutrient acquisition is crucial for sustaining life. Plants develop beneficial intracellular partnerships with arbuscular mycorrhiza (AM) and nitrogen-fixing bacteria to surmount the scarcity of soil nutrients and tap into atmospheric dinitrogen, respectively1,2. Initiation of these root endosymbioses requires symbiont-induced oscillations in nuclear calcium (Ca2+) concentrations in root cells3. How the nuclear-localized ion channels, cyclic nucleotide-gated channel (CNGC) 15 and DOESN'T MAKE INFECTIONS1 (DMI1)4 are coordinated to specify symbiotic-induced nuclear Ca2+ oscillations remains unknown. Here we discovered an autoactive CNGC15 mutant that generates spontaneous low-frequency Ca2+ oscillations. While CNGC15 produces nuclear Ca2+ oscillations via a gating mechanism involving its helix 1, DMI1 acts as a pacemaker to specify the frequency of the oscillations. We demonstrate that the specificity of symbiotic-induced nuclear Ca2+ oscillations is encoded in its frequency. A high frequency activates endosymbiosis programmes, whereas a low frequency modulates phenylpropanoid pathways. Consequently, the autoactive cngc15 mutant, which is capable of generating both frequencies, has increased flavonoids that enhance AM, root nodule symbiosis and nutrient acquisition. We transferred this trait to wheat, resulting in field-grown wheat with increased AM colonization and nutrient acquisition. Our findings reveal a new strategy to boost endosymbiosis in the field and reduce inorganic fertilizer use while sustaining plant growth.

QTL mapping and BSR-seq revealed loci and candidate genes associated with the sporadic multifoliolate phenotype in soybean (Glycine max)

KEY MESSAGE

The QTLs and candidate genes governing the multifoliolate phenotype were identified by combining linkage mapping with BSR-seq, revealing a possible interplay between genetics and the environment in soybean leaf development. Soybean, as a legume, is typified by trifoliolate leaves. Although multifoliolate leaves (compound leaves with more than three leaflets each) have been reported in soybean, including sporadic appearances in the first compound leaves in a recombinant inbred line (RIL) population from a cross between cultivated soybean C08 and wild soybean W05 from this study, the genetic basis of this phenomenon is still unclear. Here, we integrated quantitative trait locus (QTL) mapping with bulked segregant RNA sequencing (BSR-seq) to identify the genetic loci associated with the multifoliolate phenotype in soybean. Using linkage mapping, ten QTLs related to the multifoliolate trait were identified. Among these, a significant and major QTL, qMF-2-1 on chromosome 2 and consistently detected across biological replicates, explained more than 10% of the phenotypic variation. Together with BSR-seq analyses, which analyzed the RILs with the highest multifoliolate frequencies and those with the lowest frequencies as two distinct bulks, two candidate genes were identified: Glyma.06G204300 encoding the transcription factor TCP5, and Glyma.06G204400 encoding LONGIFOLIA 2 (LNG2). Transcriptome analyses revealed that stress-responsive genes were significantly differentially expressed between high-multifoliolate occurrence lines and low occurrence ones, indicating environmental factors probably influence the appearance of multifoliolate leaves in soybean through stress-responsive genes. Hence, this study offers new insights into the genetic mechanism behind the multifoliolate phenotype in soybean.

Characterization of a Powdery Mildew Resistance Gene in Wheat Breeding Line Jingzi 102 Using Bulk Segregant RNA Sequencing

Wheat (Triticum aestivum L.) is one of the most important crops worldwide. Powdery mildew caused by Blumeria graminis f. sp. tritici is a destructive disease threatening wheat yield and quality. The utilization of resistant genes and cultivars is considered the most economical, environmentally friendly, and effective method to control powdery mildew. Wheat breeding line Jingzi 102 was highly resistant to powdery mildew at both seedling and adult plant stages. Genetic analysis of F1, F2, and F2:3 populations of "Jingzi 102 × Shixin 828" showed that the resistance of Jingzi 102 against powdery mildew isolate E09 at the seedling stage was controlled by a single dominant gene, temporarily designated PmJZ. Using bulked segregant RNA sequencing combined with molecular markers analysis, PmJZ was located on the long arm of chromosome 2B and flanked by markers BJK695-1 and CIT02g-20 with the genetic distances of 1.2 and 0.5 centimorgan, respectively, corresponding to the bread wheat genome of Chinese Spring (International Wheat Genome Sequencing Consortium RefSeq v2.1) 703.8 to 707.6 Mb. PmJZ is most likely different from the documented Pm genes on chromosome 2BL based on their physical positions, molecular markers analysis, and resistance spectrum. Based on the gene annotation information, five genes related to disease resistance could be considered as the candidate genes of PmJZ. To accelerate the application of PmJZ, the flanking markers BJK695-1 and CIT02g-20 can serve for marker-assisted selection of PmJZ in wheat disease-resistance breeding.

Modern Plant Breeding Techniques in Crop Improvement and Genetic Diversity: From Molecular Markers and Gene Editing to Artificial Intelligence-A Critical Review

With the development of new technologies in recent years, researchers have made significant progress in crop breeding. Modern breeding differs from traditional breeding because of great changes in technical means and breeding concepts. Whereas traditional breeding initially focused on high yields, modern breeding focuses on breeding orientations based on different crops' audiences or by-products. The process of modern breeding starts from the creation of material populations, which can be constructed by natural mutagenesis, chemical mutagenesis, physical mutagenesis transfer DNA (T-DNA), Tos17 (endogenous retrotransposon), etc. Then, gene function can be mined through QTL mapping, Bulked-segregant analysis (BSA), Genome-wide association studies (GWASs), RNA interference (RNAi), and gene editing. Then, at the transcriptional, post-transcriptional, and translational levels, the functions of genes are described in terms of post-translational aspects. This article mainly discusses the application of the above modern scientific and technological methods of breeding and the advantages and limitations of crop breeding and diversity. In particular, the development of gene editing technology has contributed to modern breeding research.

Wheat improvement through advances in single nucleotide polymorphism (SNP) detection and genotyping with a special emphasis on rust resistance

KEY MESSAGE

Single nucleotide polymorphism (SNP) markers in wheat and their prospects in breeding with special reference to rust resistance. Single nucleotide polymorphism (SNP)-based markers are increasingly gaining momentum for screening and utilizing vital agronomic traits in wheat. To date, more than 260 million SNPs have been detected in modern cultivars and landraces of wheat. This rapid SNP discovery was made possible through the release of near-complete reference and pan-genome assemblies of wheat and its wild relatives, coupled with whole genome sequencing (WGS) of thousands of wheat accessions. Further, genotyping customized SNP sites were facilitated by a series of arrays (9 to 820Ks), a cost effective substitute WGS. Lately, germplasm-specific SNP arrays have been introduced to characterize novel traits and detect closely linked SNPs for marker-assisted breeding. Subsequently, the kompetitive allele-specific PCR (KASP) assay was introduced for rapid and large-scale screening of specific SNP markers. Moreover, with the advances and reduction in sequencing costs, ample opportunities arise for generating SNPs artificially through mutations and in combination with next-generation sequencing and comparative genomic analyses. In this review, we provide historical developments and prospects of SNP markers in wheat breeding with special reference to rust resistance where over 50 genetic loci have been characterized through SNP markers. Rust resistance is one of the most essential traits for wheat breeding as new strains of the Puccinia fungus, responsible for rust diseases, evolve frequently and globally.

Isocitrate lyase promotes Puccinia striiformis f. sp. tritici susceptibility in wheat (Triticum aestivum) by suppressing accumulation of glyoxylate cycle intermediates

Plant fungal parasites manipulate host metabolism to support their own survival. Among the many central metabolic pathways altered during infection, the glyoxylate cycle is frequently upregulated in both fungi and their host plants. Here, we examined the response of the glyoxylate cycle in bread wheat (Triticum aestivum) to infection by the obligate biotrophic fungal pathogen Puccinia striiformis f. sp. tritici (Pst). Gene expression analysis revealed that wheat genes encoding the two unique enzymes of the glyoxylate cycle, isocitrate lyase (TaICL) and malate synthase, diverged in their expression between susceptible and resistant Pst interactions. Focusing on TaICL, we determined that the TaICL B homoeolog is specifically upregulated during early stages of a successful Pst infection. Furthermore, disruption of the B homoeolog alone was sufficient to significantly perturb Pst disease progression. Indeed, Pst infection of the TaICL-B disruption mutant (TaICL-BY400*) was inhibited early during initial penetration, with the TaICL-BY400* line also accumulating high levels of malic acid, citric acid, and aconitic acid. Exogenous application of malic acid or aconitic acid also suppressed Pst infection, with trans-aconitic acid treatment having the most pronounced effect by decreasing fungal biomass 15-fold. Thus, enhanced TaICL-B expression during Pst infection may lower accumulation of malic acid and aconitic acid to promote Pst proliferation. As exogenous application of aconitic acid and malic acid has previously been shown to inhibit other critical pests and pathogens, we propose TaICL as a potential target for disruption in resistance breeding that could have wide-reaching protective benefits for wheat and beyond.

Wheat NAC transcription factor NAC5-1 is a positive regulator of senescence

Wheat (Triticum aestivum L.) is an important source of both calories and protein in global diets, but there is a trade-off between grain yield and protein content. The timing of leaf senescence could mediate this trade-off as it is associated with both declines in photosynthesis and nitrogen remobilization from leaves to grain. NAC transcription factors play key roles in regulating senescence timing. In rice, OsNAC5 expression is correlated with increased protein content and upregulated in senescing leaves, but the role of the wheat ortholog in senescence had not been characterized. We verified that NAC5-1 is the ortholog of OsNAC5 and that it is expressed in senescing flag leaves in wheat. To characterize NAC5-1, we combined missense mutations in NAC5-A1 and NAC5-B1 from a TILLING mutant population and overexpressed NAC5-A1 in wheat. Mutation in NAC5-1 was associated with delayed onset of flag leaf senescence, while overexpression of NAC5-A1 was associated with slightly earlier onset of leaf senescence. DAP-seq was performed to locate transcription factor binding sites of NAC5-1. Analysis of DAP-seq and comparison with other studies identified putative downstream target genes of NAC5-1 which could be associated with senescence. This work showed that NAC5-1 is a positive transcriptional regulator of leaf senescence in wheat. Further research is needed to test the effect of NAC5-1 on yield and protein content in field trials, to assess the potential to exploit this senescence regulator to develop high-yielding wheat while maintaining grain protein content.

A complex receptor locus confers responsiveness to necrosis and ethylene-inducing like peptides in Brassica napus

Brassica crops are susceptible to diseases which can be mitigated by breeding for resistance. MAMPs (microbe-associated molecular patterns) are conserved molecules of pathogens that elicit host defences known as pattern-triggered immunity (PTI). Necrosis and Ethylene-inducing peptide 1-like proteins (NLPs) are MAMPs found in a wide range of phytopathogens. We studied the response to BcNEP2, a representative NLP from Botrytis cinerea, and showed that it contributes to disease resistance in Brassica napus. To map regions conferring NLP response, we used the production of reactive oxygen species (ROS) induced during PTI across a population of diverse B. napus accessions for associative transcriptomics (AT), and bulk segregant analysis (BSA) on DNA pools created from a cross of NLP-responsive and non-responsive lines. In silico mapping with AT identified two peaks for NLP responsiveness on chromosomes A04 and C05 whereas the BSA identified one peak on A04. BSA delimited the region for NLP-responsiveness to 3 Mbp, containing ~245 genes on the Darmor-bzh reference genome and four co-segregating KASP markers were identified. The same pipeline with the ZS11 genome confirmed the highest-associated region on chromosome A04. Comparative BLAST analysis revealed unannotated clusters of receptor-like protein (RLP) homologues on ZS11 chromosome A04. However, no specific RLP homologue conferring NLP response could be identified. Our results also suggest that BR-SIGNALLING KINASE1 may be involved with modulating the NLP response. Overall, we demonstrate that responsiveness to NLP contributes to disease resistance in B. napus and define the associated genomic location. These results can have practical application in crop improvement.

Identification of a Locus for High-Temperature Adult-Plant Resistance to Stripe Rust in the Wheat Yr8 Near-Isogenic Line Through Mutagenesis and Molecular Mapping

Aegilops species are wheat relatives that harbor valuable disease resistance genes for wheat breeding. The wheat Yr8 near-isogenic line AvSYr8NIL has long been believed to carry only Yr8 for race-specific all-stage resistance to stripe rust, caused by Puccinia striiformis f. sp. tritici, derived from Aegilops comosa. However, AvSYr8NIL has been found to have high-temperature adult-plant (HTAP) resistance in our field and greenhouse tests. To confirm both HTAP and Yr8 resistance, seeds from AvSYr8NIL were treated with ethyl methanesulfonate to generate mutant lines. The mutant lines with only Yr8 (M641) and only HTAP resistance (M488) were crossed with the susceptible recurrent parent Avocet S (AvS). The F1 and F4 lines of AvS/M641 were phenotyped with Yr8-avirulent races in the seedling stage at the low-temperature (4 to 20°C) profile, while the F1, F2, F4, and F5 lines of AvS/M488 were phenotyped with Yr8-virulent races in the adult-plant stage at the high-temperature (10 to 30°C) profile. Both Yr8 and the HTAP resistance gene (YrM488) were recessive. The F4 populations of AvS/M641 and AvS/M488 were genotyped using polymorphic Kompetitive allele-specific PCR markers converted from single-nucleotide polymorphisms. Yr8 was mapped to a 0.66-cM fragment, and YrM488 was mapped to a 1.22-cM interval on chromosome 2D. The physical distance between the two resistance genes was estimated to be more than 500 Mb, indicating their distinct loci. The mutant lines with separated resistance genes would be useful in enhancing our understanding of different types of resistance and in further studying the interactions between wheat and the stripe rust pathogen.

Sr65: a widely effective gene for stem rust resistance in wheat

KEY MESSAGE

Sr65 in chromosome 1A of Indian wheat landrace Hango-2 is a potentially useful all-stage resistance gene that currently protects wheat from stem rust in Australia, India, Africa and Europe. Stem rust, caused by Puccinia graminis f. sp. tritici (Pgt), threatened global wheat production with the appearance of widely virulent races that included TTKSK and TTRTF. Indian landrace Hango-2 showed resistance to Pgt races in India and Australia. Screening of a Hango-2/Avocet 'S' (AvS) recombinant inbred line population identified two stem rust resistance genes, a novel gene (temporarily named as SrH2) from Hango-2 and Sr26 from AvS. A mapping population segregating for SrH2 alone was developed from two recombinant lines. SrH2 was mapped on the short arm of chromosome 1A, where it was flanked by KASP markers KASP_7944 (proximal) and KASP_12147 (distal). SrH2 was delimited to an interval of 1.8-2.3 Mb on chromosome arm 1AS. The failure to detect candidate genes through MutRenSeq and comparative genomic analysis with the pan-genome dataset indicated the necessity to generate a Hango-2 specific assembly for detecting the gene sequence linked with SrH2 resistance. MutRenSeq however enabled identification of SrH2-linked KASP marker sunCS_265. Markers KASP_12147 and sunCS_265 showed 92% and 85% polymorphism among an Australian cereal cultivar diversity panel and can be used for marker-assisted selection of SrH2 in breeding programs. The effectiveness of SrH2 against Pgt races from Europe, Africa, India, and Australia makes it a valuable resource for breeding stem rust-resistant wheat cultivars. Since no wheat-derived gene was previously located in chromosome arm 1AS, SrH2 represents a new locus and named as SR65.

QTL Mapping Reveals Both All-Stage and Adult-Plant Resistance to Powdery Mildew in Chinese Elite Wheat Cultivars

Powdery mildew caused by Blumeria graminis f. sp. tritici is a threat to wheat production in China. Mapping quantitative trait loci (QTL) for resistance to powdery mildew and developing breeder-friendly markers are important initial steps in breeding resistant cultivars. An all-stage resistance gene and several QTL were identified using a population of 254 recombinant inbred lines developed from a Jingdong 8/Aikang 58 cross. The population was evaluated for powdery mildew resistance across six field environments over three consecutive growing seasons utilizing two different mixtures of B. graminis f. sp. tritici isolates, named #Bgt-HB and #Bgt-BJ. Using genotypic data obtained from the Wheat TraitBreed 50K single-nucleotide polymorphism array, seven stable QTL were identified on chromosome arms 1DL, 2AL, 2DS, 4DL, 5AL, 6BL.1, and 6BL.2. The QTL on 2AL conferred all-stage resistance to B. graminis f. sp. tritici race E20 in greenhouse tests and explained up to 52% of the phenotypic variance in field trials but was resistant only against #Bgt-HB. The gene involved in this QTL was predicted to be Pm4a based on genome location and gene sequence. QPmja.caas-1DL, QPmja.caas-4DL, and QPmja.caas-6BL.1 were identified as potentially new QTL for powdery mildew resistance. QPmja.caas-2DS and QPmja.caas-6BL.1 were effective against both B. graminis f. sp. tritici mixtures, indicating their probable broad-spectrum resistance. A Kompetitive allele-specific PCR marker closely linked to QPmja.caas-2DS was developed and validated in a panel of 286 wheat cultivars. Because both Jingdong 8 and Aikang 58 have been leading cultivars and breeding parents, the QTL and marker reported represent valuable resources for wheat researchers and breeders.

Intelligent reprogramming of wheat for enhancement of fungal and nematode disease resistance using advanced molecular techniques

Wheat (Triticum aestivum L.) diseases are major factors responsible for substantial yield losses worldwide, which affect global food security. For a long time, plant breeders have been struggling to improve wheat resistance against major diseases by selection and conventional breeding techniques. Therefore, this review was conducted to shed light on various gaps in the available literature and to reveal the most promising criteria for disease resistance in wheat. However, novel techniques for molecular breeding in the past few decades have been very fruitful for developing broad-spectrum disease resistance and other important traits in wheat. Many types of molecular markers such as SCAR, RAPD, SSR, SSLP, RFLP, SNP, and DArT, etc., have been reported for resistance against wheat pathogens. This article summarizes various insightful molecular markers involved in wheat improvement for resistance to major diseases through diverse breeding programs. Moreover, this review highlights the applications of marker assisted selection (MAS), quantitative trait loci (QTL), genome wide association studies (GWAS) and the CRISPR/Cas-9 system for developing disease resistance against most important wheat diseases. We also reviewed all reported mapped QTLs for bunts, rusts, smuts, and nematode diseases of wheat. Furthermore, we have also proposed how the CRISPR/Cas-9 system and GWAS can assist breeders in the future for the genetic improvement of wheat. If these molecular approaches are used successfully in the future, they can be a significant step toward expanding food production in wheat crops.

Identification of a Wheat Powdery Mildew Dominant Resistance Gene in the Pm5 Locus for High-Throughput Marker-Assisted Selection

Powdery mildew, caused by Blumeria graminis f. sp. tritici (Bgt), poses a severe threat to wheat yield and quality worldwide. Rapid identification and the accurate transference of effective resistance genes are important to the development of resistant cultivars and the sustainable control of this disease. In the present study, the wheat line AL11 exhibited high levels of resistance to powdery mildew at both the seedling and adult plant stages. Genetic analysis of the AL11 × 'Shixin 733' mapping population revealed that its resistance was controlled by a single dominant gene, tentatively designated PmAL11. Using bulked segregant RNA-Seq and molecular marker analysis, PmAL11 was mapped to the Pm5 locus on chromosome 7B where it cosegregated with the functional marker Pm5e-KASP. Sequence alignment analysis revealed that the Pm5e-homologous sequence in AL11 was identical to the reported recessive gene Pm5e in wheat landrace 'Fuzhuang 30'. It appears that PmAL11 was most probably Pm5e, but it was mediated by a dominant inheritance pattern, so it should provide a valuable resistance resource for both genetic study and wheat breeding. To efficiently use and trace PmAL11 in breeding, a new kompetitive allele-specific PCR marker AL11-K2488 that cosegregated with this gene was developed and confirmed to be applicable in the different wheat backgrounds, thus promoting its use in the marker-assisted selection of PmAL11.

Transcriptome profiling at the transition to the reproductive stage uncovers stage and tissue-specific genes in wheat

BACKGROUND

The transition from vegetative to floral phase is the result of complex crosstalk of exogenous and endogenous floral integrators. This critical physiological event is the response to environmental interaction, which causes biochemical cascades of reactions at different internal tissues, organs, and releases signals that make the plant moves from vegetative status to a reproductive phase. This network controlling flowering time is not deciphered largely in bread wheat. In this study, a comparative transcriptome analysis at a transition time in combination with genetic mapping was used to identify responsible genes in a stage and tissue-specific manner. For this reason, two winter cultivars that have been bred in Germany showing contrasting and stable heading time in different environments were selected for the analysis.

RESULTS

In total, 670 and 1075 differentially expressed genes in the shoot apical meristem and leaf tissue, respectively, could be identified in 23 QTL intervals for the heading date. In the transition apex, Histone methylation H3-K36 and regulation of circadian rhythm are both controlled by the same homoeolog genes mapped in QTL TaHd112, TaHd124, and TaHd137. TaAGL14 gene that identifies the floral meristem was mapped in TaHd054 in the double ridge. In the same stage, the homoeolog located on chromosome 7D of FLOWERING TIME LOCUS T mapped on chr 7B, which evolved an antagonist function and acts as a flowering repressor was uncovered. The wheat orthologue of transcription factor ASYMMETRIC LEAVES 1 (AS1) was identified in the late reproductive stage and was mapped in TaHd102, which is strongly associated with heading date. Deletion of eight nucleotides in the AS1 promoter could be identified in the binding site of the SUPPRESSOR OF CONSTANS OVEREXPRESSION 1 (SOC1) gene in the late flowering cultivar. Both proteins AS1 and SOC1 are inducing flowering time in response to gibberellin biosynthesis.

CONCLUSION

The global transcriptomic at the transition phase uncovered stage and tissue-specific genes mapped in QTL of heading date in winter wheat. In response to Gibberellin signaling, wheat orthologous transcription factor AS1 is expressed in the late reproductive phase of the floral transition. The locus harboring this gene is the strongest QTL associated with the heading date trait in the German cultivars. Consequently, we conclude that this is another indication of the Gibberellin biosynthesis as the mechanism behind the heading variation in wheat.

Screening of Candidate Genes Associated with Brown Stripe Resistance in Sugarcane via BSR-seq Analysis

Sugarcane brown stripe (SBS), caused by the fungal pathogen Helminthosporium stenospilum, is one of the most serious threats to sugarcane production. However, its outbreaks and epidemics require suitable climatic conditions, resulting in the inefficient improvement of the SBS resistance by phenotype selection. The sugarcane F1 population of SBS-resistant YT93-159 × SBS-susceptible ROC22 was used for constructing the bulks. Bulked segregant RNA-seq (BSR-seq) was then performed on the parents YT93-159 (T01) and ROC22 (T02), and the opposite bulks of 30 SBS-susceptible individuals mixed bulk (T03) and 30 SBS-resistant individuals mixed bulk (T04) collected from 287 F1 individuals. A total of 170.00 Gb of clean data containing 297,921 SNPs and 70,426 genes were obtained. Differentially expressed genes (DEGs) analysis suggested that 7787 and 5911 DEGs were identified in the parents (T01 vs. T02) and two mixed bulks (T03 vs. T04), respectively. In addition, 25,363 high-quality and credible SNPs were obtained using the genome analysis toolkit GATK for SNP calling. Subsequently, six candidate regions with a total length of 8.72 Mb, which were located in the chromosomes 4B and 7C of sugarcane wild species Saccharum spontaneum, were identified, and 279 genes associated with SBS-resistance were annotated by ED algorithm and ΔSNP-index. Furthermore, the expression profiles of candidate genes were verified by quantitative real-time PCR (qRT-PCR) analysis, and the results showed that eight genes (LRR-RLK, DHAR1, WRKY7, RLK1, BLH4, AK3, CRK34, and NDA2) and seven genes (WRKY31, CIPK2, CKA1, CDPK6, PFK4, CBL2, and PR2) of the 20 tested genes were significantly up-regulated in YT93-159 and ROC22, respectively. Finally, a potential molecular mechanism of sugarcane response to H. stenospilum infection is illustrate that the activations of ROS signaling, MAPK cascade signaling, Ca2+ signaling, ABA signaling, and the ASA-GSH cycle jointly promote the SBS resistance in sugarcane. This study provides abundant gene resources for the SBS resistance breeding in sugarcane.

Analysis of QTLs and Candidate Genes for Tassel Symptoms in Maize Infected with Sporisorium reilianum

Heat smut is a fungal soil-borne disease caused by Sporisorium reilianum, and affects the development of male and female tassels. Our previous research found that the tassel symptoms in maize infected with Sporisorium reilianum significantly differed in inbred lines with Sipingtou blood, and exhibited stable heredity over time at multiple locations. In this study, cytological analysis demonstrated that the cellular organization structures of three typical inbred lines (Huangzao4, Jing7, and Chang7-2) showed significant discrepancies at the VT stage. QTLs that control the different symptoms of maize tassels infected with Sporisorium reilianum were located in two F₂ populations, which were constructed using three typical inbred lines. The BSA (bulked segregation analysis) method was used to construct mixed gene pools based on typical tassel symptoms. The QTLs of different symptoms of maize tassels infected with Sporisorium reilianum were detected with 869 SSR markers covering the whole maize genome. The mixed gene pools were screened with polymorphic markers between the parents. Additional SSR markers were added near the above marker to detect genotypes in partially single plants in F₂ populations. The QTL controlling tassel symptoms in the Huangzao4 and Jing7 lines was located on the bin 1.06 region, between the markers of umc1590 and bnlg1598, and explained 21.12% of the phenotypic variation with an additive effect of 0.6524. The QTL controlling the tassel symptoms of the Jing7 and Chang7-2 lines was located on the bin 2.07 region, between the markers of umc1042 and bnlg1335, and explained 11.26% phenotypic variation with an additive effect of 0.4355. Two candidate genes (ZmABP2 and Zm00001D006403) were identified by a conjoint analysis of label-free quantification proteome sequencings.

Genetic identification of SNP markers and candidate genes associated with sugarcane smut resistance using BSR-Seq

Sugarcane smut caused by Sporisorium scitamineum is one of the most severe fungal diseases worldwide. In this study, a cross was made between a smut-resistant variety YT93-159 and a smut-susceptible variety ROC22, and 312 progenies were obtained. Two bulks of progenies were then constructed, one consisted of 27 highly smut resistant progenies and the other 24 smut susceptible progenies. Total RNAs of the progenies of each bulk, were pooled and subject to bulked segregant RNA-sequence analysis (BSR-Seq). A total of 164.44 Gb clean data containing 2,341,449 SNPs and 64,999 genes were obtained, 7,295 of which were differentially expressed genes (DEGs). These DEGs were mainly enriched in stress-related metabolic pathways, including carbon metabolism, phenylalanine metabolism, plant hormone signal transduction, glutathione metabolism, and plant-pathogen interactions. Besides, 45,946 high-quality, credible SNPs, a 1.27 Mb region at Saccharum spontaneum chromosome Chr5B (68,904,827 to 70,172,982), and 129 candidate genes were identified to be associated with smut resistance. Among them, twenty-four genes, either encoding key enzymes involved in signaling pathways or being transcription factors, were found to be very closely associated with stress resistance. RT-qPCR analysis demonstrated that they played a positive role in smut resistance. Finally, a potential molecular mechanism of sugarcane and S. scitamineum interaction is depicted that activations of MAPK cascade signaling, ROS signaling, Ca²⁺ signaling, and PAL metabolic pathway and initiation of the glyoxalase system jointly promote the resistance to S. scitamineum in sugarcane. This study provides potential SNP markers and candidate gene resources for smut resistance breeding in sugarcane.

Temporally coordinated expression of nuclear genes encoding chloroplast proteins in wheat promotes Puccinia striiformis f. sp. tritici infection

Targeting host processes that allow pathogens to thrive can be invaluable in resistance breeding. Here, we generated a deep-sequencing transcriptome time course for Puccinia striiformis f. sp. tritici (Pst) infection on wheat and compared datasets from three wheat varieties with different levels of susceptibility to two tested pathogen isolates. We sought genes specifically altered in a susceptible host as candidates that might support colonisation. Host responses differed between Pst-varietal pairs most prominently early during infection. Notably, however, nuclear genes encoding chloroplast-localised proteins (NGCPs) exhibited temporal coordination of expression profiles that differed at later time points in relation to Pst susceptibility. Disrupting one such NGCP, encoding the chloroplast-localised RNA binding protein TaCSP41a, led to lower Pst susceptibility. These analyses thus highlight NGCPs as prime targets for Pst manipulation during infection and point to TaCSP41a disruption as a potential source of Pst resistance for breeding programmes.

A transcriptome time course of Puccinia striiformis f. sp. tritici (Pst) infection reveals nuclear genes encoding chloroplast-localized proteins are manipulated during infection and highlights TaCSP41a disruption as a target for resistance breeding.

Genetically Dissecting the Novel Powdery Mildew Resistance Gene in Wheat Breeding Line PBDH1607

Powdery mildew is one of the most destructive diseases in wheat production. Identifying novel resistance genes and deploying them in new cultivars is the most effective approach to minimize wheat losses caused by powdery mildew. In this study, wheat breeding line PBDH1607 showed high resistance to powdery mildew at both the seedling and adult plant stages. Genetic analysis of the seedling data demonstrated that the resistance was controlled by a single dominant gene, tentatively designated PmPBDH. The ΔSNP index based on bulked segregant RNA sequencing indicated that PmPBDH was associated with an interval of about 30.8 Mb (713.5 to 744.3 Mb) on chromosome arm 4AL. Using newly developed markers, we mapped PmPBDH to a 3.2-cM interval covering 7.1 Mb (719,055,516 to 726,215,121 bp). This interval differed from those of Pm61 (717,963,176 to 719,260,469 bp), MlIW30 (732,769,506 to 732,790,522 bp), and MlNSF10 (729,275,816 to 731,365,462 bp) reported on the same chromosome arm. PmPBDH also differed from Pm61, MlIW30, and MlNSF10 by its response spectrum, origin, or inheritance mode, suggesting that PmPBDH should be a new Pm gene. In the candidate interval, five genes were found to be associated with PmPBDH via time course gene expression analysis, and thus they are candidate genes of PmPBDH. Six closely linked markers, including two kompetitive allele-specific PCR markers, were confirmed to be applicable for tracking PmPBDH in marker-assisted breeding.

Evaluation of nine statistics to identify QTLs in bulk segregant analysis using next generation sequencing approaches

Background

Bulk segregant analysis (BSA) combined with next generation sequencing is a powerful tool to identify quantitative trait loci (QTL). The impact of the size of the study population and the percentage of extreme genotypes analysed have already been assessed. But a good comparison of statistical approaches designed to identify QTL regions using next generation sequencing (NGS) technologies for BSA is still lacking.

Results

We developed an R code to simulate QTLs in bulks of F2 contrasted lines. We simulated a range of recombination rates based on estimations using different crop species. The simulations were used to benchmark the ability of statistical methods identify the exact location of true QTLs. A single QTL led to a shift in allele frequency across a large fraction of the chromosome for plant species with low recombination rate. The smoothed version of all statistics performed best notably the smoothed Euclidean distance-based statistics was always found to be more accurate in identifying the location of QTLs. We propose a simulation approach to build confidence interval statistics for the detection of QTLs.

Conclusion

We highlight the statistical methods best suited for BSA studies using NGS technologies in crops even when recombination rate is low. We also provide simulation codes to build confidence intervals and to assess the impact of recombination for application to other studies. This computational study will help select NGS-based BSA statistics that are useful to the broad scientific community.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08718-y.

Capturing Wheat Phenotypes at the Genome Level

Recent technological advances in next-generation sequencing (NGS) technologies have dramatically reduced the cost of DNA sequencing, allowing species with large and complex genomes to be sequenced. Although bread wheat (Triticum aestivum L.) is one of the world's most important food crops, efficient exploitation of molecular marker-assisted breeding approaches has lagged behind that achieved in other crop species, due to its large polyploid genome. However, an international public-private effort spanning 9 years reported over 65% draft genome of bread wheat in 2014, and finally, after more than a decade culminated in the release of a gold-standard, fully annotated reference wheat-genome assembly in 2018. Shortly thereafter, in 2020, the genome of assemblies of additional 15 global wheat accessions was released. As a result, wheat has now entered into the pan-genomic era, where basic resources can be efficiently exploited. Wheat genotyping with a few hundred markers has been replaced by genotyping arrays, capable of characterizing hundreds of wheat lines, using thousands of markers, providing fast, relatively inexpensive, and reliable data for exploitation in wheat breeding. These advances have opened up new opportunities for marker-assisted selection (MAS) and genomic selection (GS) in wheat. Herein, we review the advances and perspectives in wheat genetics and genomics, with a focus on key traits, including grain yield, yield-related traits, end-use quality, and resistance to biotic and abiotic stresses. We also focus on reported candidate genes cloned and linked to traits of interest. Furthermore, we report on the improvement in the aforementioned quantitative traits, through the use of (i) clustered regularly interspaced short-palindromic repeats/CRISPR-associated protein 9 (CRISPR/Cas9)-mediated gene-editing and (ii) positional cloning methods, and of genomic selection. Finally, we examine the utilization of genomics for the next-generation wheat breeding, providing a practical example of using in silico bioinformatics tools that are based on the wheat reference-genome sequence.

BSA-Seq and Fine Linkage Mapping for the Identification of a Novel Locus (qPH9) for Mature Plant Height in Rice (Oryza sativa)

BACKGROUND

Plant height is a key factor in the determination of rice yield since excessive height can easily cause lodging and reduce yield. Therefore, the identification and analysis of plant height-related genes to elucidate their physiological, biochemical, and molecular mechanisms have significant implications for rice breeding and production.

RESULTS

High-throughput quantitative trait locus (QTL) sequencing analysis of a 638-individual F2:3 mapping population resulted in the identification of a novel height-related QTL (qPH9), which was mapped to a 2.02-Mb region of Chromosome 9. Local QTL mapping, which was conducted using 13 single nucleotide polymorphism (SNP)-based Kompetitive allele-specific PCR (KASP) markers for the qPH9 region, and traditional linkage analysis, facilitated the localization of qPH9 to a 126-kb region that contained 15 genes. Subsequent haplotype and sequence analyses indicated that OsPH9 was the most probable candidate gene for plant height at this locus, and functional analysis of osph9 CRISPR/Cas9-generated OsPH9 knockout mutants supported this conclusion.

CONCLUSION

OsPH9 was identified as a novel regulatory gene associated with plant height in rice, along with a height-reducing allele in 'Dongfu-114' rice, thereby representing an important molecular target for rice improvement. The findings of the present study are expected to spur the investigation of genetic mechanisms underlying rice plant height and further the improvement of rice plant height through marker-assisted selection.

Study on Supergenus Rubus L.: Edible, Medicinal, and Phylogenetic Characterization

Rubus L. is one of the most diverse genera belonging to Rosaceae; it consists of more than 700 species with a worldwide distribution. It thus provides an ideal natural "supergenus" for studying the importance of its edible, medicinal, and phylogenetic characteristics for application in our daily lives and fundamental scientific studies. The Rubus genus includes many economically important species, such as blackberry (R. fruticosus L.), red raspberry (R. ideaus L.), black raspberry (R. occidentalis L.), and raspberry (R. chingii Hu), which are widely utilized in the fresh fruit market and the medicinal industry. Although Rubus species have existed in human civilization for hundreds of years, their utilization as fruit and in medicine is still largely inadequate, and many questions on their complex phylogenetic relationships need to be answered. In this review, we briefly summarize the history and progress of studies on Rubus, including its domestication as a source of fresh fruit, its medicinal uses in pharmacology, and its systematic position in the phylogenetic tree. Recent available evidence indicates that (1) thousands of Rubus cultivars were bred via time- and labor-consuming methods from only a few wild species, and new breeding strategies and germplasms were thus limited; (2) many kinds of species in Rubus have been used as medicinal herbs, though only a few species (R. ideaus L., R. chingii Hu, and R. occidentalis L.) have been well studied; (3) the phylogeny of Rubus is very complex, with the main reason for this possibly being the existence of multiple reproductive strategies (apomixis, hybridization, and polyploidization). Our review addresses the utilization of Rubus, summarizing major relevant achievements and proposing core prospects for future application, and thus could serve as a useful roadmap for future elite cultivar breeding and scientific studies.

SNP Discovery Using BSR-Seq Approach for Spot Blotch Resistance in Wheat (Triticum aestivum L.), an Essential Crop for Food Security

The pathogenic fungus, Bipolaris sorokiniana, that causes spot blotch (SB) disease of wheat, is a major production constraint in the Eastern Gangetic Plains of South Asia and other warm, humid regions of the world. A recombinant inbred line population was developed and phenotyped at three SB-prone locations in India. The single nucleotide polymorphism (SNP) for SB resistance was identified using a bulked segregant RNA-Seq-based approach, referred to as “BSR-Seq.” Transcriptome sequencing of the resistant parent (YS#24), the susceptible parent (YS#58), and their resistant and susceptible bulks yielded a total of 429.67 million raw reads. The bulk frequency ratio (BFR) of SNPs between the resistant and susceptible bulks was estimated, and selection of SNPs linked to resistance was done using sixfold enrichments in the corresponding bulks (BFR >6). With additional filtering criteria, the number of transcripts was further reduced to 506 with 1055 putative polymorphic SNPs distributed on 21 chromosomes of wheat. Based on SNP enrichment on chromosomal loci, five transcripts were found to be associated with SB resistance. Among the five SB resistance-associated transcripts, four were distributed on the 5B chromosome with putative 52 SNPs, whereas one transcript with eight SNPs was present on chromosome 3B. The SNPs linked to the trait were exposed to a tetra-primer ARMS-PCR assay, and an SNP-based allele-specific marker was identified for SB resistance. The in silico study of these five transcripts showed homology with pathogenesis-related genes; the metabolic pathway also exhibits similar results, suggesting their role in the plant defense mechanism.

Transcriptome profiling revealed candidate genes, pathways and transcription factors related to nitrogen utilization and excessive nitrogen stress in perennial ryegrass

Ryegrass (Lolium perenne L.), a high-quality forage grass, is a good nutrient source for herbivorous livestock. However, improving nitrogen use efficiency and avoiding nitrate toxicity caused by excessive nitrogen are continual challenges in ryegrass production. The molecular mechanism underlying the response of ryegrass to nitrogen, especially excessive nitrogen, remains unclear. In this study, the transcriptomic changes under different nitrogen levels were investigated in perennial ryegrass by high-throughput next-generation RNA sequencing. Phenotypic characterization showed that treatment with half of the standard N concentration (N0.5) led to a better growth state than the other three treatments. The treatments with the standard N concentration (N1) and treatments with ten times higher than the standard N concentration (N10) contained excessive nitrogen, which placed stress on plant growth. Analysis of differentially expressed genes indicated that 345 and 104 genes are involved in the regulation of nitrogen utilization and excessive nitrogen stress, respectively. KEGG enrichment analysis suggested that "photosynthesis-antenna proteins" may respond positively to appropriate nitrogen conditions, whereas "steroid biosynthesis", "carotenoid biosynthesis" and "C5-branched dibasic acid metabolism" were identified as the top significantly enriched pathways in response to excessive nitrogen. Additionally, 21 transcription factors (TFs) related to nitrogen utilization were classified into 10 families, especially the AP2-EREBP and MYB TF families. Four TFs related to excessive nitrogen stress were identified, including LOBs, NACs, AP2-EREBPs and HBs. The expression patterns of these selected genes were also analyzed. These results provide new insight into the regulatory mechanism of ryegrass in response to nitrogen utilization and excessive nitrogen stress.

Allelically and Differentially Expressed Genes After Infection of Edwardsiella ictaluri in Channel Catfish as Determined by Bulk Segregant RNA-Seq

Identification of genetic markers associated with resistance against enteric septicemia of catfish (ESC) is of great interest for genetic enhancement programs of catfish. In the present study, bulk segregant RNA-Seq analysis was applied to determine differentially expressed genes and alleles after ESC infection. Here we report three genomic regions on LG1, LG12, and LG26, containing significant single-nucleotide polymorphisms (SNPs). These genomic regions aligned well with quantitative trait loci (QTL) previously identified. Within the QTL regions, eleven genes were found to be differentially regulated between phenotypic bulks. Importantly, the QTL on linkage group 1 (LG1) were found to be expressed in the liver, whereas the QTL on LG12 and LG26 were expressed in the intestine, suggesting multiple mechanisms of ESC resistance. It is apparent that apolipoproteins may be important for ESC resistance as the QTL on LG1 included the 14-kDa apolipoprotein genes that are both allelically expressed and differentially expressed between the resistant and susceptible bulks. Traf2 and NCK-interacting protein kinase (TNIK) were found in the QTL on LG12, and it was downregulated in resistant fish, suggesting the importance of NCK downregulation in ESC resistance, as previously reported. In addition, we observed divergent gene expression patterns between the liver and intestine after infection. Immune/inflammatory-related processes were overrepresented from liver DEGs, while those DEGs identified from intestine were enriched for proteolysis and wounding processes. Taken together, the BSR-Seq analysis presented here advanced the knowledge of ESC resistance, providing information of not only positions of QTL but also genes and their differential expression between resistant and susceptible fish, making it one step closer to the identification of the causal genes for ESC resistance.

Long-read genome sequencing of bread wheat facilitates disease resistance gene cloning

The cloning of agronomically important genes from large, complex crop genomes remains challenging. Here we generate a 14.7 gigabase chromosome-scale assembly of the South African bread wheat (Triticum aestivum) cultivar Kariega by combining high-fidelity long reads, optical mapping and chromosome conformation capture. The resulting assembly is an order of magnitude more contiguous than previous wheat assemblies. Kariega shows durable resistance to the devastating fungal stripe rust disease1. We identified the race-specific disease resistance gene Yr27, which encodes an intracellular immune receptor, to be a major contributor to this resistance. Yr27 is allelic to the leaf rust resistance gene Lr13; the Yr27 and Lr13 proteins show 97% sequence identity2,3. Our results demonstrate the feasibility of generating chromosome-scale wheat assemblies to clone genes, and exemplify that highly similar alleles of a single-copy gene can confer resistance to different pathogens, which might provide a basis for engineering Yr27 alleles with multiple recognition specificities in the future.

Fine mapping and identification of the candidate gene BFS for fruit shape in wax gourd (Benincasa hispida)

Non-synonymous mutations in the BFS gene, which encodes the IQD protein, are responsible for the shape of wax gourd fruits. Fruit shape is an important agronomic trait in wax gourds. Therefore, in this study, we employed bulked segregant analysis (BSA) to identify a candidate gene for fruit shape in wax gourds within F₂ populations derived by crossing GX-71 (long cylindrical fruit, fruit shape index = 4.56) and MY-1 (round fruit, fruit shape index = 1.06) genotypes. According to BSA, the candidate gene is located in the 17.18 Mb region on chromosome 2. Meanwhile, kompetitive allele-specific PCR (KASP) markers were used to reduce it to a 19.6 Kb region. Only one gene was present within the corresponding region of the reference genome, namely Bch02G016830 (designated BFS). Subsequently, BFS was sequenced in six wax gourd varieties with different fruit shapes. Sequence analysis revealed two non-synonymous mutations in the round wax gourd and one non-synonymous mutation in the cylindrical wax gourd. Quantitative real‑time PCR (qRT-PCR) analysis further showed that the expression of BFS in round fruits was significantly higher than in long cylindrical fruits at the ovary formation stage. Therefore, BFS is a candidate gene for determination wax gourd shape. The predicted protein encoded by the BFS gene belongs to the IQ67-domain protein family, which have the structural characteristics of scaffold proteins and coordinate Ca²⁺ CaM signaling from the membrane to the nucleus. Ultimately, two derived cleaved amplified polymorphic sequence (dCAPS) markers were developed to facilitate marker-assisted selection for wax gourds breeding.

Identification of Two New Loci for Adult Plant Resistance to Leaf Rust and Stripe Rust in the Chinese Wheat Variety 'Neimai 836'

The characterization of leaf rust (caused by Puccinia triticina) and stripe rust (caused by Puccinia striiformis f. sp. tritici) resistance genes is the basis for breeding resistant wheat varieties and managing epidemics of these diseases in wheat. A cross between the susceptible wheat variety 'Apav#1' and resistant variety 'Neimai 836' was used to develop a mapping population containing 148 F₅ recombinant inbred lines (RILs). Leaf rust phenotyping was done in field trials at Ciudad Obregón, Mexico, in 2017 and 2018, and stripe rust data were generated at Toluca, Mexico, in 2017 and in Mianyang, Ezhou, and Gansu, China, in 2019. Inclusive complete interval mapping (ICIM) was used to create a genetic map and identify significant resistance quantitative trait loci (QTL) with 2,350 polymorphic markers from a 15K wheat single-nucleotide polymorphism (SNP) array and simple-sequence repeats (SSRs). The pleiotropic multipathogen resistance gene Lr46/Yr29 and four QTL were identified, including two new loci, QLr.hzau-3BL and QYr.hzau-5AL, which explained 3 to 16% of the phenotypic variation in resistance to leaf rust and 7 to 14% of that to stripe rust. The flanking SNP markers for the two loci were converted to Kompetitive Allele-Specific PCR (KASP) markers and used to genotype a collection of 153 wheat lines, indicating the Chinese origin of the loci. Our results suggest that Neimai 836, which has been used as a parent for many wheat varieties in China, could be a useful source of high-level resistance to both leaf rust and stripe rust.

RNA-seq bulked segregant analysis combined with KASP genotyping rapidly identified PmCH7087 as responsible for powdery mildew resistance in wheat

Powdery mildew causes considerable yield losses in common wheat (Triticum aestivum L.) production. Mapping and cloning powdery mildew-resistant quantitative trait loci can benefit stable yield production by facilitating the breeding of resistant varieties. In this study, we used the powdery mildew resistance introgression line 'CH7087' (harboring the resistance gene PmCH7087) and developed a large F₂ population and a corresponding F_2:3 segregation population containing 2,000 family lines for molecular mapping of PmCH7087. Genetic analysis demonstrated that the resistance phenotype was controlled by a single dominant gene. According to the performance exhibited by the F_2:3 lines, 50 resistant lines and 50 susceptible lines without phenotype segregation were chosen for pooling and bulked segregant RNA sequencing (BSR-Seq) analysis. A region spanning 42.77 Mb was identified, and genotyping of an additional 183 F_2:3 lines with extreme phenotypes using 20 kompetitive allele-specific polymerase chain reaction (KASP) markers in the BSR-Seq mapping regions confirmed this region and narrowed it to 9.68 Mb, in which 45 genes were identified and annotated. Five of these transcripts harbored nonsynonymous single nucleotide polymorphisms between the two parents, with the transcripts of TraesCS2B01G302800 being involved in signal transduction. Furthermore, TraesCS2B01G302800.2 was annotated as the closest homologue of serine/threonine-protein kinase PBS1, a typical participant in the plant disease immune response, indicating that TraesCS2B01G302800 was the candidate gene of PmCH7087. Our results may facilitate future research attempting to improve powdery mildew resistance in wheat and to identify candidate genes for further verification and gene cloning.

Identification of a major and stable QTL on chromosome 5A confers spike length in wheat (Triticum aestivum L.)

Spike length (SL) is the key determinant of plant architecture and yield potential. In this study, 193 recombinant inbred lines (RILs) derived from a cross between 13F10 and Chuanmai 42 (CM42) were evaluated for spike length in six environments. Sixty RILs consisting of 30 high and 30 low SLs were genotyped using the bulked segregant analysis exome sequencing (BSE-Seq) analysis for preliminary quantitative trait locus (QTL) mapping. A 6.69 Mb (518.43-525.12 Mb) region on chromosome 5AL was found to have a significant effect on the SL trait. Fifteen competitive allele-specific PCR (KASP) markers were successfully converted from the single nucleotide polymorphisms (SNPs) in the SL target region. Combined with four novel simple sequence repeat (SSR) markers, a genetic linkage map spanning 21.159 cM was constructed. The mapping result confirmed the identity of a major and stable QTL named QSl.cib-5A in the targeted region that explained 7.88-26.60% of the phenotypic variation in SL. QSl.cib-5A was narrowed to a region of 4.84 cM interval corresponding to a 4.67 Mb (516.60-521.27 Mb) physical region in the Chinese Spring RefSeq v2.0 containing 17 high-confidence genes with 25 transcripts. In addition, this QTL exhibited pleiotropic effects on spikelet density (SD), with the phenotypic variances proportion ranging from 11.34 to 19.92%. This study provides a foundational step for cloning the QSl.cib-5A, which is involved in the regulation of spike morphology in common wheat.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01249-6.

RNA-Seq-based DNA marker analysis of the genetics and molecular evolution of Triticeae species

The release of high-quality chromosome-level genome sequences of members of the Triticeae tribe has greatly facilitated genetic and genomic analyses of important crops such as wheat (Triticum aestivum) and barley (Hordeum vulgare). Due to the large diploid genome size of Triticeae plants (ca. 5 Gbp), transcript analysis is an important method for identifying genetic and genomic differences among Triticeae species. In this review, we summarize our results of RNA-Seq analyses of diploid wheat accessions belonging to the genera Aegilops and Triticum. We also describe studies of the molecular relationships among these accessions and provide insight into the evolution of common hexaploid wheat. DNA markers based on polymorphisms within species can be used to map loci of interest. Even though the genome sequence of diploid Aegilops tauschii, the D-genome donor of common wheat, has been released, the diploid barley genome continues to provide key information about the physical structures of diploid wheat genomes. We describe how a series of RNA-Seq analyses of wheat relatives has helped uncover the structural and evolutionary features of genomic and genetic systems in wild and cultivated Triticeae species.

The branched-chain amino acid aminotransferase TaBCAT1 modulates amino acid metabolism and positively regulates wheat rust susceptibility

Plant pathogens suppress defense responses to evade recognition and promote successful colonization. Although identifying the genes essential for pathogen ingress has traditionally relied on screening mutant populations, the post-genomic era provides an opportunity to develop novel approaches that accelerate identification. Here, RNA-seq analysis of 68 pathogen-infected bread wheat (Triticum aestivum) varieties, including three (Oakley, Solstice and Santiago) with variable levels of susceptibility, uncovered a branched-chain amino acid aminotransferase (termed TaBCAT1) as a positive regulator of wheat rust susceptibility. We show that TaBCAT1 is required for yellow and stem rust infection and likely functions in branched-chain amino acid (BCAA) metabolism, as TaBCAT1 disruption mutants had elevated BCAA levels. TaBCAT1 mutants also exhibited increased levels of salicylic acid (SA) and enhanced expression of associated defense genes, indicating that BCAA regulation, via TaBCAT1, has a key role in SA-dependent defense activation. We also identified an association between the levels of BCAAs and resistance to yellow rust infection in wheat. These findings provide insight into SA-mediated defense responses in wheat and highlight the role of BCAA metabolism in the defense response. Furthermore, TaBCAT1 could be manipulated to potentially provide resistance to two of the most economically damaging diseases of wheat worldwide.

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.

Comparative Transcriptomics and RNA-Seq-Based Bulked Segregant Analysis Reveals Genomic Basis Underlying Cronartium ribicola vcr2 Virulence

Breeding programs of five-needle pines have documented both major gene resistance (MGR) and quantitative disease resistance (QDR) to Cronartium ribicola (Cri), a non-native, invasive fungal pathogen causing white pine blister rust (WPBR). WPBR is one of the most deadly forest diseases in North America. However, Cri virulent pathotypes have evolved and can successfully infect and kill trees carrying resistance (R) genes, including vcr2 that overcomes MGR conferred by the western white pine (WWP, Pinus monticola) R gene (Cr2). In the absence of a reference genome, the present study generated a vcr2 reference transcriptome, consisting of about 20,000 transcripts with 1,014 being predicted to encode secreted proteins (SPs). Comparative profiling of transcriptomes and secretomes revealed vcr2 was significantly enriched for several gene ontology (GO) terms relating to oxidation-reduction processes and detoxification, suggesting that multiple molecular mechanisms contribute to pathogenicity of the vcr2 pathotype for its overcoming Cr2. RNA-seq-based bulked segregant analysis (BSR-Seq) revealed genome-wide DNA variations, including about 65,617 single nucleotide polymorphism (SNP) loci in 7,749 polymorphic genes shared by vcr2 and avirulent (Avcr2) pathotypes. An examination of the distribution of minor allele frequency (MAF) uncovered a high level of genomic divergence between vcr2 and Avcr2 pathotypes. By integration of extreme-phenotypic genome-wide association (XP-GWAS) analysis and allele frequency directional difference (AFDD) mapping, we identified a set of vcr2-associated SNPs within functional genes, involved in fungal virulence and other molecular functions. These included six SPs that were top candidate effectors with putative activities of reticuline oxidase, proteins with common in several fungal extracellular membrane (CFEM) domain or ferritin-like domain, polysaccharide lyase, rds1p-like stress responsive protein, and two Cri-specific proteins without annotation. Candidate effectors and vcr2-associated genes provide valuable resources for further deciphering molecular mechanisms of virulence and pathogenicity by functional analysis and the subsequent development of diagnostic tools for monitoring the virulence landscape in the WPBR pathosystems.

From Genetic Maps to QTL Cloning: An Overview for Durum Wheat

Durum wheat is one of the most important cultivated cereal crops, providing nutrients to humans and domestic animals. Durum breeding programs prioritize the improvement of its main agronomic traits; however, the majority of these traits involve complex characteristics with a quantitative inheritance (quantitative trait loci, QTL). This can be solved with the use of genetic maps, new molecular markers, phenotyping data of segregating populations, and increased accessibility to sequences from next-generation sequencing (NGS) technologies. This allows for high-density genetic maps to be developed for localizing candidate loci within a few Kb in a complex genome, such as durum wheat. Here, we review the identified QTL, fine mapping, and cloning of QTL or candidate genes involved in the main traits regarding the quality and biotic and abiotic stresses of durum wheat. The current knowledge on the used molecular markers, sequence data, and how they changed the development of genetic maps and the characterization of QTL is summarized. A deeper understanding of the trait architecture useful in accelerating durum wheat breeding programs is envisioned.

Identification of Solanum Immune Receptors by Bulked Segregant RNA-Seq and High-Throughput Recombinant Screening

The identification, understanding, and deployment of immune receptors are crucial to achieve high-level and durable resistance for crops against pathogens. In potato, many R genes have been identified using map-based cloning strategies. However, this is a challenging and laborious task that involves the development of a high number of molecular markers for the initial mapping, and the screening of thousands of plants for fine mapping. Bulked segregant RNA-Seq (BSR-Seq) has proven to be an efficient technique for the mapping of resistance genes. The RNA from two bulks of plants with contrasting phenotypes is sequenced and analyzed to identify single-nucleotide polymorphism (SNPs) markers linked to the target gene. Subsequently, the SNP markers that are identified can be used to delimit the mapping interval. Additionally, we designed an in vitro recombinant screening strategy that is advantageous for analyzing a large number of plants, in terms of time, space, and cost. Tips and detailed protocols, including BSR-Seq, bioinformatic analysis, and recombinant screening, are provided in this chapter.

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.

A haplotype-led approach to increase the precision of wheat breeding

Crop productivity must increase at unprecedented rates to meet the needs of the growing worldwide population. Exploiting natural variation for the genetic improvement of crops plays a central role in increasing productivity. Although current genomic technologies can be used for high-throughput identification of genetic variation, methods for efficiently exploiting this genetic potential in a targeted, systematic manner are lacking. Here, we developed a haplotype-based approach to identify genetic diversity for crop improvement using genome assemblies from 15 bread wheat (Triticum aestivum) cultivars. We used stringent criteria to identify identical-by-state haplotypes and distinguish these from near-identical sequences (~99.95% identity). We showed that each cultivar shares ~59 % of its genome with other sequenced cultivars and we detected the presence of extended haplotype blocks containing hundreds to thousands of genes across all wheat chromosomes. We found that genic sequence alone was insufficient to fully differentiate between haplotypes, as were commonly used array-based genotyping chips due to their gene centric design. We successfully used this approach for focused discovery of novel haplotypes from a landrace collection and documented their potential for trait improvement in modern bread wheat. This study provides a framework for defining and exploiting haplotypes to increase the efficiency and precision of wheat breeding towards optimising the agronomic performance of this crucial crop.

Brinton, Uauy and colleagues utilize genomic data from the 10+ Wheat Genome Project to develop a useful tool for studying and generating new wheat cultivars. This framework uses advanced exploitation of wheat haplotypes to bring newfound precision and efficiency to wheat breeding.

Evidence for the Accumulation of Nonsynonymous Mutations and Favorable Pleiotropic Alleles During Wheat Breeding

Plant breeding leads to the genetic improvement of target traits by selecting a small number of genotypes from among typically large numbers of candidate genotypes after careful evaluation. In this study, we first investigated how mutations at conserved nucleotide sites normally viewed as deleterious, such as nonsynonymous sites, accumulated in a wheat, Triticum aestivum, breeding lineage. By comparing a 150 year old ancestral and modern cultivar, we found recent nucleotide polymorphisms altered amino acids and occurred within conserved genes at frequencies expected in the absence of purifying selection. Mutations that are deleterious in other contexts likely had very small or no effects on target traits within the breeding lineage. Second, we investigated if breeders selected alleles with favorable effects on some traits and unfavorable effects on others and used different alleles to compensate for the latter. An analysis of a segregating population derived from the ancestral and modern parents provided one example of this phenomenon. The recent cultivar contains the Rht-B1b green revolution semi-dwarfing allele and compensatory alleles that reduce its negative effects. However, improvements in traits other than plant height were due to pleiotropic loci with favorable effects on traits and to favorable loci with no detectable pleiotropic effects. Wheat breeding appears to tolerate mutations at conserved nucleotide sites and to only select for alleles with both favorable and unfavorable effects on traits in exceptional situations.

Identification and validation of genetic loci for tiller angle in bread wheat

KEY MESSAGE

Two major QTL for tiller angle were identified on chromosomes 1AL and 5DL, and TaTAC-D1 might be the candidate gene for QTA.caas-5DL. An ideal plant architecture is important for achieving high grain yield in crops. Tiller angle (TA) is an important factor influencing yield. In the present study, 266 recombinant inbred lines (RILs) derived from a cross between Zhongmai 871 (ZM871) and its sister line Zhongmai 895 (ZM895) was used to map TA by extreme pool-genotyping and inclusive composite interval mapping (ICIM). Two quantitative trait loci (QTL) on chromosomes 1AL and 5DL were identified with reduced tiller angle alleles contributed by ZM895. QTA.caas-1AL was detected in six environments, explaining 5.4-11.2% of the phenotypic variances. The major stable QTL, QTA.caas-5DL, was identified in all eight environments, accounting for 13.8-24.8% of the phenotypic variances. The two QTL were further validated using BC1F4 populations derived from backcrosses ZM871/ZM895//ZM871 (121 lines) and ZM871/ZM895//ZM895 (175 lines). Gene TraesCS5D02G322600, located in the 5DL QTL and designated TaTAC-D1, had a SNP in the third exon with 'A' and 'G' in ZM871 and ZM895, respectively, resulting in a Thr169Ala amino acid change. A KASP marker based on this SNP was validated in two sets of germplasm, providing further evidence for the significant effects of TaTAC-D1 on TA. Thus extreme pool-genotyping can be employed to detect QTL for plant architecture traits and KASP markers tightly linked with the QTL can be used in wheat breeding programs targeting improved plant architecture.

The loss of function of HEL, which encodes a cellulose synthase interactive protein, causes helical and vine-like growth of tomato

Helical growth is an economical way for plant to obtain resources. The classic microtubule-microfibril alignment model of Arabidopsis helical growth involves restriction of the appropriate orientation of cellulose microfibrils appropriately in the cell walls. However, the molecular mechanism underlying tomato helical growth remains unknown. Here, we identified a spontaneous tomato helical (hel) mutant with right-handed helical cotyledons and petals but left-handed helical stems and true leaves. Genetic analysis revealed that the hel phenotype was controlled by a single recessive gene. Using map-based cloning, we cloned the HEL gene, which encodes a cellulose interacting protein homologous to CSI1 of Arabidopsis. We identified a 27 bp fragment replacement that generated a premature stop codon. Transgenic experiments showed that the helical growth phenotype could be restored by the allele of this gene from wild-type Pyriforme. In contrast, the knockout mutation of HEL in Pyriforme via CRISPR/Cas9 resulted in helical growth. These findings shed light on the molecular control of the helical growth of tomato.

A rare single nucleotide variant in Pm5e confers powdery mildew resistance in common wheat

Powdery mildew poses severe threats to wheat production. The most sustainable way to control this disease is through planting resistant cultivars. We report the map-based cloning of the powdery mildew resistance allele Pm5e from a Chinese wheat landrace. We applied a two-step bulked segregant RNA sequencing (BSR-Seq) approach in developing tightly linked or co-segregating markers to Pm5e. The first BSR-Seq used phenotypically contrasting bulks of recombinant inbred lines (RILs) to identify Pm5e-linked markers. The second BSR-Seq utilized bulks of genetic recombinants screened from a fine-mapping population to precisely quantify the associated genomic variation in the mapping interval, and identified the Pm5e candidate genes. The function of Pm5e was validated by transgenic assay, loss-of-function mutants and haplotype association analysis. Pm5e encodes a nucleotide-binding domain leucine-rich-repeat-containing (NLR) protein. A rare nonsynonymous single nucleotide variant (SNV) within the C-terminal leucine rich repeat (LRR) domain is responsible for the gain of powdery mildew resistance function of Pm5e, an allele endemic to wheat landraces of Shaanxi province of China. Results from this study demonstrate the value of landraces in discovering useful genes for modern wheat breeding. The key SNV associated with powdery mildew resistance will be useful for marker-assisted selection of Pm5e in wheat breeding programs.

Weeping candidate genes screened using comparative transcriptomic analysis of weeping and upright progeny in an F1 population of Prunus mume

Weeping is a specific plant architecture with high ornamental value. Despite the considerable importance of the weeping habit to landscaping applications and knowledge of plant architecture biology, little is known regarding the underlying molecular mechanisms. In this study, growth and phytohormone content were analyzed among the progeny of different branch types in an F1 mapping population of Prunus mume with varying plant architecture. Bulked segregant RNA sequencing was conducted to compare differences among progeny at a transcriptional level. The weeping habit appears to be a complex process regulated by a series of metabolic pathways, with photosynthesis and flavonoid biosynthesis highly enriched in differentially expressed genes between weeping and upright progeny. Based on functional annotation and homologous analyses, we identified 30 candidate genes related to weeping that merit further analysis, including 10 genes related to IAA and GA₃ biosynthesis, together with 6 genes related to secondary branch growth. The results of this study will facilitate further studies of the molecular mechanisms underlying the weeping habit in P. mume.

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.

RLP/K enrichment sequencing; a novel method to identify receptor-like protein (RLP) and receptor-like kinase (RLK) genes

The identification of immune receptors in crop plants is time-consuming but important for disease control. Previously, resistance gene enrichment sequencing (RenSeq) was developed to accelerate mapping of nucleotide-binding domain and leucine-rich repeat containing (NLR) genes. However, resistances mediated by pattern recognition receptors (PRRs) remain less utilized. Here, our pipeline shows accelerated mapping of PRRs. Effectoromics leads to precise identification of plants with target PRRs, and subsequent RLP/K enrichment sequencing (RLP/KSeq) leads to detection of informative single nucleotide polymorphisms that are linked to the trait. Using Phytophthora infestans as a model, we identified Solanum microdontum plants that recognize the apoplastic effectors INF1 or SCR74. RLP/KSeq in a segregating Solanum population confirmed the localization of the INF1 receptor on chromosome 12, and led to the rapid mapping of the response to SCR74 to chromosome 9. By using markers obtained from RLP/KSeq in conjunction with additional markers, we fine-mapped the SCR74 receptor to a 43-kbp G-LecRK locus. Our findings show that RLP/KSeq enables rapid mapping of PRRs and is especially beneficial for crop plants with large and complex genomes. This work will enable the elucidation and characterization of the nonNLR plant immune receptors and ultimately facilitate informed resistance breeding.

Identification of Stripe Rust Resistance Loci in U.S. Spring Wheat Cultivars and Breeding Lines Using Genome-Wide Association Mapping and Yr Gene Markers

Stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst), poses a major threat to wheat production worldwide, especially in the United States. To identify loci for effective stripe rust resistance in U.S. wheat, a genome-wide association study (GWAS) was conducted using a panel of 616 spring wheat cultivars and breeding lines. The accessions in this panel were phenotyped for stripe rust response in the greenhouse at seedling stage with five predominant and highly virulent races of Pst and in different field environments at adult-plant stage in 2017 and 2018. In total, 2,029 single-nucleotide polymorphism markers that cover the whole genome were generated with genotyping by multiplexed sequencing and used in GWAS. In addition, 23 markers of previously reported resistance genes or quantitative trait loci (QTLs) were used to genotype the population. This spring panel was grouped into three subpopulations based on principal component analysis. A total of 37 genes or QTLs including 10 potentially new QTLs for resistance to stripe rust were detected by GWAS and linked marker tests. The frequencies of the resistance genes or QTLs in various nurseries were determined, indicating different intensities of these genes or QTLs used in breeding programs of different regions. These resistance loci and the information on their markers, effectiveness, and distributions should be useful for improving stripe rust resistance in wheat cultivars.

Bulked Segregant RNA-Seq Reveals Distinct Expression Profiling in Chinese Wheat Cultivar Jimai 23 Responding to Powdery Mildew

Wheat powdery mildew, caused by Blumeria graminis f. sp. tritici (Bgt), is one of the most destructive fungal diseases threatening global wheat production. Host resistance is well known to be the most efficient method to control this disease. However, the molecular mechanism of wheat powdery mildew resistance (Pm) is still unclear. To analyze the molecular mechanism of Pm, we used the resistant wheat cultivar Jimai 23 to investigate its potential resistance components and profiled its expression in response to powdery mildew infection using bulked segregant RNA-Seq (BSR-Seq). We showed that the Pm of Jimai 23 was provided by a single dominant gene, tentatively designated PmJM23, and assigned it to the documented Pm2 region of chromosome 5DS. 3,816 consistently different SNPs were called between resistant and susceptible parents and the bulked pools derived from the combinations between the resistant parent Jimai23 and the susceptible parent Tainong18. 58 of the SNPs were assigned to the candidate region of PmJM23. Subsequently, 3,803 differentially expressed genes (DEGs) between parents and bulks were analyzed by GO, COG and KEGG pathway enrichment. The temporal expression patterns of associated genes following Bgt inoculation were further determined by RT-qPCR. Expression of six disease-related genes was investigated during Bgt infection and might serve as valuable genetic resources for the improvement of durable resistance to Bgt.

Characterization of Genetic and Allelic Diversity Amongst Cultivated and Wild Lentil Accessions for Germplasm Enhancement

Intensive breeding of cultivated lentil has resulted in a relatively narrow genetic base, which limits the options to increase crop productivity through selection. Assessment of genetic diversity in the wild gene pool of lentil, as well as characterization of useful and novel alleles/genes that can be introgressed into elite germplasm, presents new opportunities and pathways for germplasm enhancement, followed by successful crop improvement. In the current study, a lentil collection consisting of 467 wild and cultivated accessions that originated from 10 diverse geographical regions was assessed, to understand genetic relationships among different lentil species/subspecies. A total of 422,101 high-confidence SNP markers were identified against the reference lentil genome (cv. CDC Redberry). Phylogenetic analysis clustered the germplasm collection into four groups, namely, Lens culinaris/Lens orientalis, Lens lamottei/Lens odemensis, Lens ervoides, and Lens nigricans. A weak correlation was observed between geographical origin and genetic relationship, except for some accessions of L. culinaris and L. ervoides. Genetic distance matrices revealed a comparable level of variation within the gene pools of L. culinaris (Nei’s coefficient 0.01468–0.71163), L. ervoides (Nei’s coefficient 0.01807–0.71877), and L. nigricans (Nei’s coefficient 0.02188–1.2219). In order to understand any genic differences at species/subspecies level, allele frequencies were calculated from a subset of 263 lentil accessions. Among all cultivated and wild lentil species, L. nigricans exhibited the greatest allelic differentiation across the genome compared to all other species/subspecies. Major differences were observed on six genomic regions with the largest being on Chromosome 1 (c. 1 Mbp). These results indicate that L. nigricans is the most distantly related to L. culinaris and additional structural variations are likely to be identified from genome sequencing studies. This would provide further insights into evolutionary relationships between cultivated and wild lentil germplasm, for germplasm improvement and introgression.