Next Article in Journal
A Short ERAP2 That Binds IRAP Is Expressed in Macrophages Independently of Gene Variation
Next Article in Special Issue
Genome-Wide Identification of the Eucalyptus urophylla GATA Gene Family and Its Diverse Roles in Chlorophyll Biosynthesis
Previous Article in Journal
An In Silico Methodology That Facilitates Decision Making in the Engineering of Nanoscale Protein Materials
Previous Article in Special Issue
Functional Characterization of MdTAC1a Gene Related to Branch Angle in Apple (Malus x domestica Borkh.)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Insights into the Genetic Architecture and Genomic Prediction of Powdery Mildew Resistance in Flax (Linum usitatissimum L.)

1
Ottawa Research and Development Centre, Agriculture and Agri-Food Canada, Ottawa, ON K1A 0C6, Canada
2
Morden Research and Development Centre, Agriculture and Agri-Food Canada, Morden, MB R6M 1Y5, Canada
3
Department of Biology, University of Ottawa, 30 Marie Curie, Ottawa, ON K1N 6N5, Canada
4
Key Laboratory of Crop Genetics and Germplasm Enhancement, College of Agriculture, Nanjing Agricultural University/JCIC-MCP, Nanjing 210095, China
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(9), 4960; https://doi.org/10.3390/ijms23094960
Submission received: 8 April 2022 / Revised: 28 April 2022 / Accepted: 28 April 2022 / Published: 29 April 2022
(This article belongs to the Collection Genetics and Molecular Breeding in Plants)

Abstract

:
Powdery mildew (PM), caused by the fungus Oidium lini in flax, can cause defoliation and reduce seed yield and quality. To date, one major dominant gene (Pm1) and three quantitative trait loci (QTL) on chromosomes 1, 7 and 9 have been reported for PM resistance. To fully dissect the genetic architecture of PM resistance and identify QTL, a diverse flax core collection of 372 accessions augmented with an additional 75 breeding lines were sequenced, and PM resistance was evaluated in the field for eight years (2010–2017) in Morden, Manitoba, Canada. Genome-wide association studies (GWAS) were performed using two single-locus and seven multi-locus statistical models with 247,160 single nucleotide polymorphisms (SNPs) and the phenotypes of the 447 individuals for each year separately as well as the means over years. A total of 349 quantitative trait nucleotides (QTNs) were identified, of which 44 large-effect QTNs (R2 = 10–30%) were highly stable over years. The total number of favourable alleles per accession was significantly correlated with PM resistance (r = 0.74), and genomic selection (GS) models using all identified QTNs generated significantly higher predictive ability (r = 0.93) than those constructed using the 247,160 genome-wide random SNP (r = 0.69), validating the overall reliability of the QTNs and showing the additivity of PM resistance in flax. The QTNs were clustered on the distal ends of all 15 chromosomes, especially on chromosome 5 (0.4–5.6 Mb and 9.4–16.9 Mb) and 13 (4.7–5.2 Mb). To identify candidate genes, a dataset of 3230 SNPs located in resistance gene analogues (RGAs) was used as input for GWAS, from which an additional 39 RGA-specific QTNs were identified. Overall, 269 QTN loci harboured 445 RGAs within the 200 Kb regions spanning the QTNs, including 45 QTNs located within the RGAs. These RGAs supported by significant QTN/SNP allele effects were mostly nucleotide binding site and leucine-rich repeat receptors (NLRs) belonging to either coiled-coil (CC) NLR (CNL) or toll interleukin-1 (TIR) NLR (TNL), receptor-like kinase (RLK), receptor-like protein kinase (RLP), transmembrane-coiled-coil (TM-CC), WRKY, and mildew locus O (MLO) genes. These results constitute an important genomic tool for resistance breeding and gene cloning for PM in flax.

1. Introduction

Diseases are major constraints to crop production. Powdery mildew (PM) is a common fungal disease of many crops worldwide. Most PM pathogen species are host-specific or infect few hosts, implying that their genomes encode distinct “toolboxes” of pathogenesis-related genes [1]. Importantly, all PM fungi are obligate biotrophic plant pathogens; hence, their growth and reproduction are entirely reliant on the availability of water and nutrients from living host cells [2]. PM interferes with plant growth and reduces the crop’s quality. For example, PM can reduce grain yields by up to 40% in wheat [3]. PM occurs later in the season when temperatures are between 20 and 25 °C, and the extent of the damages is only determined after the crop has been harvested [4]. Once the host is infected, the infection spreads quickly under favourable conditions, resulting in PM outbreaks. Furthermore, the evolution of pathogen races with increased virulence may result in a “breakdown” of resistance. For instance, Pm17, Pm3a, and Pm4a were defeated in several Eastern and mid-Atlantic regions of the United States [5,6], as was Pm8 in China [7]. As a result, most resistance genes become ineffective after a period of time [8]. Thus, it is necessary to keep looking for new sources of resistance or create new combinations to stay abreast of new races. This can be achieved by identifying new genes/alleles that will allow us to build molecular tools to quickly and efficiently introduce them into breeding lines and to diversify the resistance sources against the rapidly evolving pathogen races.
Plant disease resistance is usually categorized as either “qualitative,” defined by the presence of major resistant (R) genes, or “quantitative” defined by the presence of resistance-related quantitative trait loci (QTL) [9,10]. R genes, for example, have been used successfully in wheat where they form the basis of resistance breeding programs that have produced many resistant commercial varieties [11]. The majority of the known R proteins have nucleotide-binding site (NBS) and leucine-rich repeat (LRR) domains and act as intracellular immune receptors that recognize their cognate effectors directly or indirectly [12]. Two other genetically specified R proteins are cell surface-localized receptor-like transmembrane proteins (RLPs) and receptor-like kinases (RLKs) [13]. Mildew locus O (MLO) is a well-studied PM susceptibility gene that was first discovered in barley in 1942 [14] and later identified in rice and wheat [15,16].
Qualitative resistance means “monogenic”, “vertical” or “race-specific”, while quantitative resistance refers to adult plant resistance (APR), “slow-mildewing”, or “partial resistance” [17]. Qualitative resistance genes are usually short-lived due to the frequent changes in the pathogen population [9]. Quantitative genetics approaches, such as estimating genetic elements, heritability, and efficient gene numbers, are typically used to investigate quantitative resistance. Several studies have identified more than 100 QTL on all chromosomes in wheat [11]. Quantitative resistance is more durable [18,19], and offers long-term defence to host plants. However, detecting such resistance is difficult, particularly when R genes are present. Traditionally, genetic studies of quantitative traits use segregating biparental populations that have been tested for the traits of interest and genotyped with DNA-based molecular markers. This method necessitates the costly and time-consuming creation of specific mapping populations. Additionally, the resolution is limited by the number of crossovers and the high linkage disequilibrium (LD), necessitating additional research to fine-map the QTL that often span several cM [20].
Genome-wide association studies (GWAS) are a powerful alternative strategy to the traditional linkage-based QTL mapping, focusing on LD obtained from unrelated genotypes of a collection that reflects historical recombination, thereby resulting in more accurate positioning of QTL and higher mapping resolution. However, GWAS are less effective at detecting alleles with small effects and rare alleles, and can produce a large number of false-positive associations. To improve the statistical mapping power of GWAS, a large population size and high marker density are needed. Several GWAS analyses for PM resistance have been performed in plant species such as wheat [11,21,22], barley [23] and oat [24]. These studies revealed that resistance to PM was mainly mediated by a few minor QTL with small effects and heavily influenced by the genetic background of the populations studied, the phenotyping conditions, and the genotype-by-environment interactions. These findings imply the limited efficiency of marker-assisted selection (MAS) to improve PM resistance by pyramiding small-effect favourable alleles. Genomic selection (GS), also known as genomic prediction (GP), is a promising alternative in crops, particularly for improving complex traits. In GS, genome-wide markers are used to predict the genomic-estimated breeding values (GEBVs) of individuals by capturing the benefits of both major- and minor-effect QTL [25,26]. As a result, GS captures a more significant proportion of the genetic variance of the selected traits than MAS, which is often restricted to a small number of markers linked to major QTL.
In flax, PM, caused by the obligate biotrophic ascomycete Oidium lini Skoric, is one of the most destructive and common flax foliar diseases [27,28]. It was first reported in Western Canada in 1997 [29] where it remains a sporadic disease. Infection with PM at an early stage of development will result in significant yield and seed quality reductions [30]. Most Canadian flax varieties are moderately resistant to PM under field conditions with natural inoculum [31]. One major dominant gene for resistance to PM, designated Pm1, has been identified from the Canadian varieties AC Watson, AC McDuff and AC Emerson as well as from the introduced varieties Atalante and Linda. Two additional dominant genes have also been postulated in Linda [32]. The use of resistant varieties in combination with a systematic disease management program is the most successful way to reduce the incidence of PM and save producers money.
In the present study, both GWAS and GS analyses were performed on 447 flax accessions. Field resistance assessment was performed for five to eight years (2010–2017) in Morden, Manitoba, Canada and genotyping was achieved through short-read sequencing of genomic DNA. The major objectives of this study were: (1) to discover major and minor QTNs associated with PM resistance in flax through GWAS using seven widely used multi-locus and two single-locus statistical methods; (2) to identify putative candidate genes for these QTNs; and (3) to evaluate the efficiency of the QTNs as markers in GS models.

2. Results

2.1. Evaluation of Powdery Mildew Resistance

The average PM rating of the 372 accessions of the core collection over five years (2012–2016) was 4.84 ± 1.35, while that of the 75 selected breeding lines was 2.24 ± 0.43 (Figure 1A, Table S1). The overall average was 4.5 ± 1.4 over the five years with large phenotypic variation (coefficient of variation (CV): 28.1–53.2%) (Table 1).
Based on the average PM ratings over five years, 85 genotypes were either highly resistant (HR) or resistant (R); this group included 70 of the 75 selected breeding lines and 15 accessions from the core collection. Of these 85 genotypes, 83 belong to the linseed morphotype while F_LTU_B_CN10111 and F_NLD_C_CN18983 are fibre-type flax. A total of 92 genotypes were moderately resistant (MR), including the five selected breeding lines; 114 genotypes were moderately susceptible (MS); 118 genotypes were susceptible (S); and 38 genotypes were highly susceptible (HS) (Figure 1B). The PM ratings averaged 2.2 ± 0.4, 3.4 ± 0.3, 4.4 ± 0.3, 5.8 ± 0.6 and 7.3 ± 0.3 for the R, MR, MS, S and HS groups, respectively (Figure 1B).
Similar performance for PM resistance between years was observed, although the yea-to-year PM ratings of the genotypes differed slightly. The PM ratings in 2014 and 2015 were generally lower than those in other years (Table 1), while the PM ratings of the 75 breeding lines were lower in 2010 (Table S1). Pearson correlation analysis corroborated the significant correlations across years (r = 0.35–0.81, p < 0.01) (Figure 2).

2.2. Genetic Structure of the Population

Principal component analysis (PCA) was performed for the flax core collection of 372 accessions and 75 selected breeding lines using all 247,160 SNPs identified from the 447 genotypes (Figure S1). The first three principal components (PCs) accounted for 11.8%, 8.1% and 6.1% of the total variation, respectively. These three PCs grouped all genotypes into either linseed or fibre morphotypes (Figure S1). Most of the 75 selected breeding lines were highly resistant to PM (Figure 1A), creating a spatial cluster that differed from the core collection (Figure S1). To capture more of the structural contributions of the PCs, the first ten PCs (38% of the total variation) were selected as a genetic structure matrix for the downstream GWAS analyses.

2.3. Identification of QTL

Two single-locus methods (GLM and MLM) and seven multi-locus methods (pLARmEB, pKWmEB, FASTmrMLM, ISIS EM-BLASSO, mrMLM, FASTmrEMMA, and FarmCPU) were used to identify PM resistance QTNs from all 247,160 SNPs and all 447 accessions. A total of 349 unique QTNs were identified for the six PM datasets (PM-2012, PM-2013, PM-2014, PM-2015, PM-2016, PM-Mean) based on the 47,564 haplotype blocks identified from the 247,160 SNPs using Plink (Table S2). QTNs located in the same haplotype blocks were treated as QTL, and the QTNs with the largest effect (R2) were chosen as the tag to represent these QTL. Singleton QTNs were considered independent QTL. Hereafter, we use the tag QTNs or singleton QTNs to represent the QTL. The 349 tag QTNs and their related information are listed in Table S3 and depicted on chromosomes in Figure 3.
The 349 unique QTNs were identified using a number of GWAS statistical models (Table S4). The single-locus model GLM detected only nine QTNs, but these had relatively large effects, ranging from 0.65 to 27.62% and averaging 11.34% of R2 (Table S4). No QTNs were detected with the single-locus model MLM. The multi-locus FarmCPU method detected 15 QTNs with relatively large effects (9.07% of R2). The remaining six mrMLM-based multi-locus models identified most of the QTNs, including some with large and minor effects accounting for 3.19–5.79% of R2. The QTN numbers varied from 68 with FASTmrEMMA to 108 with pLARmEB (Table S4). Only 2–4 QTNs were common between GLM and the seven multi-locus models; in contrast, 13–39 QTNs were shared by two of seven multi-locus models (Table S5).
The identified QTNs also differed across individual PM phenotypic datasets (Table S6). A total of 45–84 QTNs were identified from the six PM datasets, of which few were shared by more than one PM dataset (Table S7). However, when assessing the effects of the QTN alleles across PM datasets, most QTNs were found to be stable over years. Of the 349 unique QTNs, 265 had significant allele effects in at least three of the PM datasets (Table S3). A total of 122, 56, 54, 33, 32 and 52 QTNs had significant allele effects in six, five, four, three, two and one PM datasets, respectively. QTN effect R2 estimates were positively correlated with the number of PM datasets with significant allele effects and negatively correlated to the coefficients of variation (CVs) of the QTN effects, indicative of the stability of the QTNs (Figure 4). Large-effect QTNs were stable across years, while the small-effect QTNs only significant in one PM dataset were least stable (Figure 4, Table S3). Of the 122 QTNs stable across all six PM datasets, 44 were of large QTN effects ≥ 10% of R2, of which the following 11 had R2 values ≥ 20%: Lu2-1672205 (30.2%), Lu3-581507 (29.71%), Lu5-11130392 (27.8%), Lu5-12090990 (27.6%), Lu5-15697144 (26.9%), Lu5-16602027 (25.2%), Lu5-16840013 (22.6%), Lu7-17007593 (21.3%), Lu9-3920670 (20.9%), Lu9-20701159 (20.9%), Lu11-17188390 (20.7%).
Most large-effect QTNs were clustered on the distal ends of chromosomes, especially on chromosome 5 (0.4–5.6 Mb and 9.4–16.9 Mb) and 13 (2.6–4.9 Mb) (Figure 3, Table S3). Of the 44 large-effect QTNs with R2 ≥ 10%, 15 were located on chromosome 5 and 5 on chromosome 13.

2.4. Favourable Alleles

According to the effect direction (positive or negative) of two alleles for a QTN, the favourable allele (FA) composition of all identified QTNs existing within each genotype was determined. For all 349 unique QTNs, the total number of favourable alleles (NFAs) within a genotype ranged from 141 to 305. The NFAs were significantly negatively correlated with the PM ratings of the accessions (Figure 5A) (R2 = 0.62). The 75 selected breeding lines that were resistant (R) or moderately resistant (MR) were clearly distinguished through NFAs of accessions. The NFAs for five groups of genotypes were 262 ± 35, 204 ± 22, 193 ± 15, 180 ± 13 and 169 ± 12 for R, MR, MS, S and HS, respectively (Figure 5B), confirming the same linear correlation between NFAs and PM ratings displayed in Figure 5A.
Based on the number of favourable alleles in the genotypes of the population, three types of QTNs were observed (Figure 6, Table S3): type 1 QTNs with low favourable allele frequencies (FAFs) (Figure 6A), type 2 QTNs with high FAFs (Figure 6B), and type 3 QTNs with FAFs close to 0.5 (Figure 6C). Of the 447 genotypes in this study, 177 or 37% were R or MR; thus, the FAFs of good QTNs were expected to be approximately 0.37. We found that the 33 QTNs with R2 > 20% had an average FAF of 0.27 (0.15–0.33), and the 33 QTNs with R2 between 10% and 20% had an average FAF of 0.48 (0.19–0.86), compared to the remaining QTNs with R2 ≤ 10% that had an average FAF of 0.60 (0.06–0.94). The FAFs of most large-effect QTNs (>10%) ranged from 0.1 to 0.6 in all genotypes and 0.5 to 0.9 in resistant genotypes (Figure 7).

2.5. Relationship between PM Resistance and Flax Morphotypes

Of the 447 accessions used in this study, 367 and 80 belonged to the linseed and fibre morphotypes, respectively. Linseed and fibre genotypes had average PM ratings of 4.2 ± 1.5 and 5.3 ± 1.5, and average NFAs of 208.8 ± 38.1 and 175.7 ± 12.2, respectively (Figure 8). Significant differences in PM ratings and NFAs between linseed and fibre genotypes were observed (Wilcox test, p = 2.487 × 10−8 for PM rating and p < 2.2 × 10−16 for NFAs). The 10% most resistant accessions were all linseed types, while the 10% most susceptible included 18 fibre accessions (~38%) even though fibre accessions made up only 17.9% of the overall collection. Overall, a higher proportion of fibre accessions than expected based on their representation in the collection were susceptible, and this was congruent with the proportionally lower NFAs.

2.6. Candidate Genes

Of the 349 unique QTNs, 93 were located within genes (Table S3), of which 6 were RGAs: Lus10002249 on Lu8, Lus10015649 on Lu14, Lus10003971 on Lu7, Lus10006772 (RLK) on Lu12, and Lus10030587 and Lus10033608 (RLK) on Lu12. Lus10002249 and Lus10015649 are Toll/interleukin-1 receptor (TIR)-nucleotide-binding site (NBS)-leucine-rich repeats (LRR) (TNL) genes, and Lus10003971 is a transmembrane (TM) coiled-coil (CC) (TM-CC) gene, while the remaining three encode receptor-like protein kinase (RLK).
To further identify candidate genes for PM resistance, a subset of 3230 SNPs located on 838 flax RGAs were extracted from the overall SNP dataset and GWAS analysis was performed with the same models as mentioned above. Significant QTNs were identified from 42 RGAs, including Lus10030587, Lus10003971, and Lus10002249 that had already been detected with the 247,160 genome-wide SNP dataset. Combining the results from both datasets, a total of 388 QTNs were obtained, including 132 QTNs identified on 45 RGAs (Table 2) and 87 non-RGA genes (Table S8).
Of these 388 QTNs, 269 harboured 445 RGAs within 200 Kb regions of the QTNs (Table S8, Figure S2), encompassing 13 gene families: dirigent protein (DIR), disease resistance-zinc finger-chromosome condensation (DZC), extreme-drug-resistant (EDR), mildew resistance locus o (MLO), RLK, receptor-like protein (RLP), resistance to powdery mildew 8 (RPW8), TM-CC, TIR, nucleotide-binding site–leucine-rich repeats (NL), coiled-coil–nucleotide-binding site–leucine-rich repeats (CNL), TNL and WRKY (Table S8). Of these 445 candidate RGAs, 45 had QTNs identified within them, and 270 RGAs were supported by at least one SNP on each gene, which had significant SNP allele effects (Table S8, Figure S3). Candidate genes for large-effect QTNs were mostly TNL, CNL, WRKY, RLP and RLK genes (Table 3). For the 269 QTN regions harbouring RGAs, the average minimum distance between a QTN and a predicted candidate RGA was 32,346 bp with a range of 0–99,046 bp.
A total of 16 candidate RGA clusters that contained at least 3 candidate genes were observed on 11 of the 15 flax chromosomes (2, 3, 5, 7–10, 12–15) (Figure S2). The largest candidate gene cluster was associated with QTN Lu8-18351964 (R2 = 11.49%) located on Chr 8. It contained 15 tandemly duplicated TNL genes within a 126.5 Kb region (18,254,394–18,380,935 bp). The QTN Lu8-18351964 was identified within gene Lus10007812, while the remaining 14 genes had at least one SNP per gene with a significant allele effect on PM resistance. Another important gene cluster was associated with QTNs Lu5-1534998 (R2 = 27.83%) and Lu5-1535619 (R2 = 32.06%) and spanned a genomic region of 123.5 Kb (1,449,598–1,573,096 bp) containing four TNL genes: Lus10004726, Lus10004727, Lus10004719 and Lus10004747. The QTN Lu5-1535619 was identified within the coding region of Lus10004726, but the other three genes also had SNPs with significant allele effects on PM resistance (Table S8).
With the exception of Lus10040576, Lus10021001, and Lus10031043, the remaining 42/45 candidate RGA genes that co-located with QTNs were considered highly stable over years with CV < 100% (Table 2). These 45 RGAs belonged to the following gene families: TNL (16), RLK (16), WRKY (4), TM-CC (4), RLP (2), CNL (1), DIR (1), and MLO (1). Among those with the highest R2 (>20%) (Figure S2), Lus10004727, Lus10004726 and Lus10004719 represent tandemly duplicated TNL genes (R2 = 24–28%); Lus10029860 is another TNL gene (R2 = 35%); and Lus10032303 encodes a WRKY gene (R2 = 24%). All five genes were located on chromosome 5 (1.5–13.3 Mb). The large-effect candidate RGA Lus10027903 on chromosome 12 encodes an RLP gene (R2 = 20%).
A notable exception is the 108.6 Kb genomic region of chromosome 13 (Figure 9), which harboured three tandem duplicate RPW8 orthologous genes (Lus10000835, Lus10000836 and Lus10009328). Indeed, we identified four major QTNs at this locus, Lu13-4791823 (R2 = 13.79%), Lu13-4830850 (R2 = 7.17%), Lu13-4866704 (R2 = 12.24%), and Lu13-4900476 (R2 = 10.9%), but no SNPs were found within the three RPW8 orthologous genes, probably due to genome sequencing resulting in the low read coverage in the region. An alternative reason could be that the reads did not precisely map to the three orthologues and that the SNPs were not called because they did not pass the filtering criteria.
From the candidate RGA genes, 77 gene pairs were duplicated between or within chromosomes due to whole-genome duplication (Table S9). Duplicated genes involved 131 QTNs on all 15 chromosomes (Figure 3, Table S9). Out of 44 large-effect QTNs (R2 > 10%), 10 (Lu3-1644588, Lu4-11526385, Lu5-1552921, Lu5-1763832, Lu5-15697144, Lu7-17007593, Lu9-4948236, Lu10-11682031, Lu10-11695343, Lu15-46304) involved 13 pairs of duplicated genes, including 1 pair each of TNL genes, CNL genes, MLO genes, TM-CC genes and WRKY genes as well as 8 pairs of RLK genes (Figure 3, Table S9). Some QTNs harboured two or more candidate genes, and thus involved more than one pair of duplicated genes.

2.7. Genomic Prediction for PM Resistance

To evaluate the overall reliability of the identified QTNs, ten GS models were compared using the 349 unique QTNs and a five-fold cross-validation scheme was used to identify the best models for genomic prediction of PM resistance. All models performed similarly, except for RKHS and RFR (Figure S4). GBLUP, BRR and SVR generated the highest predictive ability of 0.93. Thus, GBLUP was used for the remaining comparisons.
Five marker sets were used to construct the GBLUP GS models: (1) 349 QTNs identified from the genome-wide dataset of 247,160 SNPs, (2) 388 QTNs obtained using the genome-wide and the RGA-derived SNP datasets, (3) 132 QTNs located within genes, including 45 RGAs and 87 non-RGA genes, (4) 294 QTNs or SNPs located on 294 candidate RGAs with significant allele effects, and (5) the genome-wide dataset of 247,160 SNPs. The 349 QTN and the 388 QTN datasets explained the highest genetic variation for PM resistance (h2 = 0.69–0.70) and generated the highest predictive ability (r = 0.925–0.926) (Table 4). The predictive ability obtained from the other datasets was significantly lower and the lowest was obtained from the genome-wide dataset of 247,160 SNPs (r = 0.690).
The prediction model constructed using all 447 accessions as a training population, the PM ratings over five years, and all 349 QTNs with GBLUP explained 96% of PM variation (R2 = 0.96), showing high predictive ability and the potential of this model in applied genomic prediction (Figure S5).

3. Discussion

Genetic studies of PM resistance in flax are few in number, limiting our understanding of the genetic architecture of this trait to a few major genes. Indeed, using traditional genetic analyses, the single dominant gene Pm1 that confers resistance to PM was identified in several Canadian (‘AC Watson’, ‘AC McDuff’, and ‘AC Emerson’) and introduced (‘Atalante’ and ‘Linda’) cultivars, and two putative dominant genes were additionally postulated in ‘Linda’ [32]. In agreement with the latter, a QTL mapping study using biparental F3 and F4 populations derived from a cross between the susceptible cultivar NorMan and the resistant cultivar Linda identified three PM resistance QTL [27]. Here, we exploited a large genetic panel that included 372 accessions from the diverse flax core collection [33,34] from which we defined a high density genome-wide dataset of 247,160 SNPs and an RGA-specific subset of 3230 SNPs, with the view of gaining greater insights into the genetic architecture of PM resistance in flax in order to design breeding methods for its improvement. Because the core collection contained few highly resistant lines, we augmented the germplasm with 75 selected breeding lines previously phenotyped and deemed highly resistant to PM. This significantly enhanced the detection power and the reliability of the identified QTNs through an increase in population size and the genetic variation of the PM ratings. We also used two single- and seven multi-locus statistical models to identify both large- and small-effect QTNs, resulting in a total of 349 unique QTNs from the genome-wide SNP dataset (Table S3). This has proven to be a good strategy to take advantage of the strength of each model as well as to palliate their shortcomings [35,36]. Post-identification QTN analysis (Figure S4) also contributed to the reliability of the QTNs by removing potentially redundant and false-positive QTNs. Overall, the methodology used herein constitutes a powerful strategy, combining different tools and methods to identify the most reliable QTN–trait association.
Three PM-resistant QTL located on LG1 (QPM-crc-LG1), 7 (QPM-crc-LG7), and 9 (QPM-crc-LG9) were identified using biparental F3 and F4 families derived from a cross between NorMan (PM-susceptible) and Linda (PM-resistant) [27]. The SSR markers defining these three QTL correspond to three genomic regions of 16,920,407–18,739,647 bp on Chr 1, 3,817,603–3,817,863 bp on Chr 7, and 357,191–357,510 bp on Chr 9 [37]. The genomic region of QPM-crc-LG1 spanned two QTNs identified in this study, Lu1-18139539 and Lu1-18452203, of which Lu1-18139539 was detected within a gene Lus10009706 (coordinated from 18,138,555 to 18,140,489 bp on Chr 1), encoding a tetratricopeptide repeat (TPR)-like superfamily protein. However, no QTNs identified here overlapped with QPM-crc-LG7 and QPM-crc-LG9. Because the Pm1 gene [32] was identified using phenotypic data only, we were not able to map this gene to a chromosome.
Current GWAS methods are limited to predicting genes or genetic features controlling traits. This may mostly depend on the density of genome-wide markers, the size of a genetic panel, the association of QTN with the trait or the extent of the QTN effect, and so on. To date, the main method for predicting candidate genes from QTNs identified by GWAS remains to scan the annotated genes in their vicinity. An optimal window size may be determined based on linkage disequilibrium decay [38,39]. In this study, we adopted a window of 200 Kb flanking a QTN [35,40]. To provide further support to the candidate RGAs as it relates to their function in disease resistance, we first narrowed the candidate genes to the flax RGAs identified from its annotated reference genome sequence [37,41,42]. Second, we reduced the set to the candidate RGAs located within the specified window size of 200 Kb and that had at least one SNP with significant allele effects on PM ratings within the RGAs. That is, even though no QTN was identified from a candidate gene, the SNP(s) on a candidate gene must have been significantly related to PM rating to be considered. Third, we performed GWAS using the 3230 SNPs exclusively located within flax RGAs in order to identify QTNs associated with PM resistance that were specifically located within RGAs. In this way, we identified 45 RGAs that possessed intragenic QTNs and an additional 270 RGAs supported by SNPs within these genes, all with significant allele effects.
QTL mapping provides a useful statistical genetics tool to identify QTNs and candidate genes associated with PM resistance. Molecular markers can be designed from some large-effect QTNs for MAS. Alternatively, or in addition to, all QTN markers can be used to establish GP models for predicting GEBVs of germplasm or breeding lines. In this study, we counted favourable alleles of all QTNs (ignoring QTN effects) in each accession of the genetic panel and we observed a significant correlation between the number of favourable alleles and PM ratings (R2 = 0.62, Figure 5A). PM-resistant accessions had significantly more favourable alleles than PM-susceptible accessions (Figure 5B), clearly indicating the significant additive feature of the identified QTNs. Consequently, the pyramiding of favourable alleles through cross prediction [43], hybridization and recombination, as well as GS is expected to translate into cultivars with improved resistance to PM. The 75 selected breeding lines with high PM resistance and a high number of favourable alleles demonstrate the potential of this approach. Genomic cross prediction is an advanced genomic tool that integrates computer simulation and GS to predict the genetic performance of different types of crosses by evaluating the expected breeding values and genetic variances of their segregating populations for the purpose of selecting superior crosses and consequently enhancing the potential for success [43]. GS has been evaluated via a cross-validation approach for agronomic, abiotic and biotic stress-related traits, including pasmo resistance in flax [36,44,45]. For example, the predictive ability of pasmo resistance in flax was 0.92 when 500 QTL were used for prediction [44], similar to our results of 0.925 obtained using 349 QTNs for PM resistance. Therefore, these QTNs offer the potential for PM resistance breeding using genomics-assisted breeding methods.
The potential candidate genes can be further validated using functional approaches. Once functionally validated, they can be genetically edited to improve cultivar PM resistance. For example, MLO is a resistance gene that confers PM resistance in many crops, such as grapevine [46], wheat [47], and barley [48]. Several technologies such as genome editing and targeting induced lesions in genomes (TILLING) have been used to create MLO mutants with enhanced PM resistance in bread wheat [47,49]. Nekrasov et al. (2017) [50] reported a non-transgenic tomato variety resistant to PM (Oidium neolycopersici) produced through editing the MLO gene (SlMlo1) using the CRISPR/Cas9 technology, which is based on the Cas9 DNA nuclease guided to a specific DNA target by a single guide-RNA (sgRNA). PMR4 encodes a callose synthase, and its loss-of-function mutants are resistant to PM in Arabidopsis and tomato. The CRISPR/Cas9-mediated knockout mutants of the PMR4 ortholog (SlPMR4) in tomato showed partial resistance against the PM pathogen O. neolycopersici [51]. RNA silencing of SlPMR4 also enhanced the resistance to PM in tomatoes [52]. In this study, we detected six MLO gene orthologues (Lus10036120, Lus10036121, Lus10023506, Lus10015461, Lus10012698, and Lus10001336) that had significant SNP effects with R2 of 1–10%. QTN Lu13-1749576 (R2 = 1.23%) was identified from the gene region of the MLO ortholog Lus10001336. These MLO and other identified candidate genes could be candidates for genome editing to improve PM resistance.
The RPW8 domain was found in several broad-spectrum powdery mildew resistance proteins from Arabidopsis thaliana and other dicots [53,54]. The A. thaliana locus RPW8 contains two naturally polymorphic, dominant R genes, RPW8.1 and RPW8.2, which individually control resistance to a broad range of powdery mildew pathogens [53]. Therefore, the three RPW8 genes identified on chromosome 13 are potentially important candidate genes. An additional PM resistance study based on a fibre flax recombinant inbred line (RIL) population derived from the Aurore/Adelie cross also showed that several RPW8 genes may explain the PM resistance difference between the two parents [55].
It is worth noting that some paired candidate duplicate genes co-located with some paired QTNs (Table S9). For instance, in Figure 9, three RPW8 genes (Lus10000835, Lus10000836 and Lus10009328) are paired duplicate genes with sequence similarity of 0.67–0.69 (Table S9). Correspondingly, four QTNs Lu13-4791823, Lu13-4830850, Lu13-4866704, and Lu13-4900476 are paired QTNs. This finding is worthy of further investigations to test whether these duplicated genes and QTNs contribute additively to PM resistance in flax.
We observed that in the combined population of 372 accessions from the core collection and 75 breeding lines, linseed accessions tended to be more resistant to PM than fibre accessions. The same result was also observed in the flax core collection alone [56]. This result is somewhat surprising because fibre flax is grown in coastal, humid environments that are conducive to the development of the disease and resistance to PM has long been a major objective of fibre flax breeding programs considering that PM has the potential to considerably reduce fibre quality. Therefore, crosses with resistant linseed varieties can be used to introduce new allelic diversity for PM resistance in fibre flax breeding programs and genomic selection can be used to assist in selecting the most resistant lines while preserving the fibre morphotype.

4. Materials and Methods

4.1. Genetic Population

A diverse genetic panel of 372 accessions from the previously described flax core collection [33,34] was used. The core collection was assembled from the world collection of 3378 flax accessions, collected from 39 countries and corresponding to 11 regions of the world: North America, South America, Eastern Asia, Western Asia, Southern Asia, Central and Eastern Europe, Western Europe, Southern Europe, Northern Europe, Oceania, and Africa. This panel contained 17 landraces, 84 breeding lines, 234 cultivars, and 37 accessions of unknown improvement status that were grouped into two morphotypes: 80 fibre and 292 linseed accessions [56].
The core collection [56] contained few highly resistant accessions. To empower the GWAS and GP analyses, an additional 75 resistant breeding lines were added to the core collection.

4.2. Phenotyping of Powdery Mildew Resistance and Statistical Analysis

The 372 accessions were evaluated for reaction to PM at Agriculture and Agri-Food Canada, Morden Research and Development Centre’s farm, Morden, Manitoba, Canada from 2012 to 2016. The experimental design was a randomized, complete block design with two replications. Each entry was seeded in 3 m rows spaced 30 cm apart during the 2nd or 3rd week of May every year. Inoculated susceptible plants were transplanted from the growth room into the field at the early flowering stage to serve as inoculum and ensure early disease infection and development in the field. One pot containing ten heavily infected plants was transplanted every ten rows. Disease ratings (PM severity) on leaves and stems were evaluated as the percentage of the leaf and stem areas covered by mycelium using the following 0 to 9 scale: 0 (HR) = no sign of PM—most vigorous plants, 1 (HR) = <10%, 2 (R) = 11–20%, 3 (MR) = 21–30%, 4 (MS) = 31–40%, 5 (S) = 41–50%, 6 (S) = 51–60%, 7 (S) = 61–70%, 8 (HS) = 71–80%, and 9 (HS) ≥ 80% [32]. This rating was converted into five resistance groups: (1) resistant (including highly resistant HR and resistant R, rating 1–2.9); (2) moderately resistant (MR, rating 3.0–3.9); (3) moderately susceptible (MS, rating 4.0–4.9); (4) susceptible (S, rating 5.0–6.9); and (5) highly susceptible (HS, rating 7.0–9.0). Field assessments were conducted at the early (PM1, 26 July) and late flowering stages (PM2, 7–10 days after PM1, around 2 August), the green boll stage (PM3, 7–10 days after PM2, around 10 August), and the early brown boll stage (PM4, 7–10 days after PM3, around 20 August). In 2012 and 2013, the observations from all four stages (PM1–PM4) were collected. Due to weather and other factors, only the data from PM1 and PM2 were available in subsequent years.
The 75 selected breeding lines were field-evaluated using the same procedure as the flax core collection, but this was accomplished during eight consecutive years (2010–2017). As both populations were assessed in the same PM nursery using the same procedure, the two datasets of the five common years (2012–2016) were combined for downstream analyses.

4.3. Genotyping and SNP Identification

Whole-genome resequencing methodology was employed to genotype all individuals of the core collection. The lines were grown in growth chambers with 20 h light at 22° and 4 h dark at 18° until they were approximately 7–8 cm tall. The compact tip of the plants (75–100 mg) was collected, flash-frozen in liquid nitrogen, and immediately lyophilized. Genomic DNA was extracted using the DNeasy 96 Plant kit (Qiagen, Mississauga, ON, Canada) and quantified using the Quant-iT PicoGreen dsDNA assay kit (Thermo Fisher Scientific, Waltham, MA, USA), both according to the manufacturer’s instructions. The Illumina HiSeq 2000 platform (Illumina Inc., San Diego, CA, USA) was used to generate 100 bp paired-end reads with ~15.5X genome coverage equivalents of the reference genome. All reads from each individual of the population were aligned to the scaffold sequences of the flax reference genome [57] using BWA v0.6.1 [58] with base-quality Q score in Phred scale >20 and other default parameters. The alignment file for each individual was used as input for SNP discovery using the software package SAMtools v1.12 [59]. All variants were further filtered to obtain a set of high-quality SNPs as previously described [60]. The SNP coordinates were then converted to the chromosome scale of the flax pseudomolecules v2.0 upon its release [42]. All procedures were implemented in the AGSNP pipeline [61,62] and its updated GBS version [60]. The detected SNPs were further filtered with minor allele frequency (MAF) > 0.05 and SNP call rate ≥ 60%. To minimize the contribution from regions of extensive strong linkage disequilibrium (LD), a single SNP was retained per 200 Kb window when pairwise correlation coefficients (r2) among neighbouring SNPs were greater than 0.8 [63,64], resulting in a total of 258,873 SNPs. Missing SNPs (on average, 14.13% of a missing data rate) were imputed using Beagle v.4.2 with default parameters [65].
A similar approach was used for the genotyping of the 75 selected breeding lines. Shotgun PCR-free library preparation (Lucigen, Middleton, WI, USA), library quality control (Illumina, San Diego, CA, USA) and sequencing were performed by the Centre d’expertise et de services Génome Québec (Montréal, QC, Canada). Twenty-one samples were indexed per lane, and sequencing on the HiSeqX platform (Illumina) generated 150 bp paired-end reads to an estimated average of 15X genome coverage per genotype. The same procedure was used to identify SNPs from the 75 selected breeding lines. The identified SNPs were combined with the SNP set from the core collection, resulting in a common set of 247,160 SNPs that was used in all downstream analyses.

4.4. Genome-Wide Association Studies (GWAS)

GWAS analyses were conducted separately for the five individual year datasets and the 5-year average dataset with the single-locus models GLM [66] and MLM [67] and the following seven multi-locus models: FarmCPU [68], pLARmEB [69], pKWmEB [70], FASTmrMLM [71], FASTmrEMMA [72], ISIS EM-BLASSO [73] and mrMLM [71]. The R package rMVP v1.0.4 [68] was used to run the GLM, MLM and FarmCPU models, while the R package mrMLM v5.0.1 [74] was used to run the remaining six multi-locus models. The kinship matrix used for each model was calculated using the module implemented in the corresponding software. The population structure of the 447 accessions was estimated using principal component analysis (PCA) and the first ten principal components (PCs) as cofactors in the models.
The threshold of significant associations for GLM, MLM and FarmCPU was determined by a critical p value (α = 0.05) subjected to Bonferroni correction, i.e., the corrected p-value = 2.02 × 10−7 (0.05/247,160 SNPs). A log of odds (LOD) score of 3.0 was used to detect significant associations for the six models implemented in the mrMLM package.
The putative QTNs identified were first analysed by testing the statistical significance of QTN alleles for PM resistance. Statistically significant differences between alleles provided validity to the QTNs. Wilcox non-parametric tests were performed using the R function wilcox.test to remove the non-significant QTNs at a 5% probability level. The direction (positive or negative) of QTN effects was subsequently determined. Only QTNs with consistent effect directions in all datasets were considered valid and were retained. Such QTNs were grouped into QTL by calculating haplotype blocks using plink v1.9 [75]. QTNs located in the same haplotype block were grouped into QTL. For each such defined QTL, the QTN with the largest average R2 over all datasets was chosen as the tag QTN representative of the QTL. R2 values were calculated based on simple regressions of QTNs on PM ratings, representing the proportion of the total variation for PM resistance explained by the QTNs/QTL.
To analyse the stability of the QTNs, we used the coefficient of variation (CV) of the allele effect values across the six PM datasets (five individual years and mean over years) for each QTN. Favourable alleles were determined based on the difference in PM rating of individuals with either of the two alleles at a given QTN. Stable QTNs have the same favourable allele in all PM datasets. QTNs with inconsistent effect values were considered unreliable and were removed. QTNs were declared highly stable if the effect differences across datasets were significant. Thus, for a significant QTN, its effect difference must be significant in at least one PM dataset. CV for effect values over datasets can be used to measure the stability of QTNs. The second criterion was R2 (%), i.e., the variance proportion of the phenotypic variation explained by the QTN. Thus, for each QTN, CV and R2 were used to describe its stability and the extent of its effect. A stable QTN is defined here as having R2 > 10% and CV < 100%.
To test QTN effect additivity, the number of QTNs with negative effect or favourable alleles (NFA) in all accessions was tallied. A simple regression of NFA on PM in the population was calculated. Correlations between the NFA and PM rating in the six datasets were calculated using the R function “cor”. A complete description of the post-identification QTN analysis pipeline is depicted in Figure S6.

4.5. Candidate Gene Prediction

A total of 1327 RGAs were identified using the RGAugury pipeline in the flax pseudomolecule [41,42]. An additional 117 disease-resistance-related genes were detected in the flax genome based on a homology search against the annotated Arabidopsis genes. Thus, a total of 1444 RGAs of known annotated functions to disease resistance were used for candidate gene prediction. To predict candidate resistance genes co-localized with QTNs, the RGAs located within 200 Kb of a QTN, i.e., 100 Kb on either side of the QTN, were investigated. In addition, GWAS were also performed using the models previously mentioned and the 3230 SNP subset present within 838 of the 1444 RGAs and that segregated in the germplasm under study.

4.6. Genomic Prediction (GP)

To gain an understanding of the best method for GP for PM resistance in flax, we evaluated the predictive ability of several marker sets: (1) all 247,160 SNPs, (2) QTNs identified by GWAS using all 247,160 SNPs, (3) QTNs within RGAs, (4) QTNs identified by GWAS using the 3230 SNPs located within RGAs, and (5) QTNs identified from all 247,160 SNPs and the 3230 SNPs located within RGAs. The following ten GP models were compared to find the optimal prediction models for PM resistance prediction: rrBLUP [76], GBLUP [77], RFR, BRR, BL, SVR [78], RKHS [76], BayesA, BayesB, and BayesC [25].
Five-fold cross-validation with 50 iterations was used to estimate the predictive ability of the models for the 372 accessions in the core collection and 75 selected breeding lines. The predictive ability (r) was calculated as the Pearson’s correlation coefficient between the mean genomic estimated breeding values (GEBVs) and the observed phenotypes. A custom genomic selection pipeline (GSPipeline v1.0) integrating the ten genomic prediction models implemented in the R packages rrBLUP v4.6.1 [76], BGLR v1.0.9 [79], BLR v1.6 [80], randomForest v4.7-1 [81] and sommer v4.1.6 [82] was used for GP model construction and cross-validation. Tukey’s multiple pairwise comparisons (HSD.test function) were performed to test the statistical significance of the predictive ability.

5. Conclusions

This study reaffirmed the quantitative nature of PM resistance in flax. The flax core collection is a valuable genetic panel that contains a broad range of genetic variation from diverse geographical origins and different improvement status, and has demonstrated, herein, that it also possesses considerable genetic variation for PM resistance. However, few accessions were highly resistant and its complementation with 75 highly resistant breeding lines was beneficial to our study. This large genetic panel with its high-density genome-wide SNPs combined with multiple single- and multi-locus GWAS models proved powerful to identify large- and minor-effect QTL, resulting in the identification of 349 QTNs and 445 candidate RGAs associated with PM resistance in flax. GWAS using the small RGA-derived SNP set further refined the candidate gene set down to 45 RGAs, which harboured QTNs within their coding sequences. Other candidate RGAs near QTNs were also supported by at least one significant SNP within their coding sequence. Significant additive features of the identified QTNs facilitate the application of these QTNs in marker-assisted and genomic selection, and a high predictive ability is expected for PM resistance. This large-scale QTL identification study provides great potential to use the identified QTNs and their potential candidate genes for PM resistance breeding and gene cloning.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms23094960/s1.

Author Contributions

F.M.Y., S.C. and K.Y.R. conceived and designed the study. K.Y.R. implemented the phenotyping of powdery mildew resistance. S.C. performed next-generation sequencing. F.M.Y. conducted SNP identification and analysis. F.M.Y., C.Z., N.K., Z.Y., P.L., J.X. and L.H. performed data analysis. F.M.Y. and C.Z. prepared tables and figures. F.M.Y. drafted the manuscript. All authors reviewed and edited the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Genome Canada and other industrial stakeholders for the Total Utilization Flax GENomics (TUFGEN) project, and by Agriculture Development Fund for Diverse Field Crop Cluster project in genomics and molecular markers to identify resistance genes in flax.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data are available in the Supplementary Materials.

Acknowledgments

We thank Ismael Moumen for his technical assistance to prepare part of the figures, manuscript review and English editing. We also thank Tara Edwards for English editing.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Schulze-Lefert, P.; Panstruga, R. A molecular evolutionary concept connecting nonhost resistance, pathogen host range, and pathogen speciation. Trends Plant Sci. 2011, 16, 117–125. [Google Scholar] [CrossRef] [PubMed]
  2. Panstruga, R. Establishing compatibility between plants and obligate biotrophic pathogens. Curr. Opin. Plant Biol. 2003, 6, 320–326. [Google Scholar] [CrossRef]
  3. Alam, M.A.; Mandal, M.S.N.; Wang, C.; Lim, W. Chromosomal location and SSR markers of a powdery mildew resistance gene in chromosome wheat line N0308. Afr. J. Micro. Res. 2013, 7, 477–482. [Google Scholar]
  4. Saharan, G.S.; Mehta, N.K.; Meena, P.D. Powdery Mildew Disease of Crucifers: Biology, Ecology and Disease Management; Springer Nature Singapore Pvt. Ltd.: Singapore, 2019. [Google Scholar]
  5. Niewoehner, A.S.; Leath, S. Virulence of Blumeria graminis f. sp. tritici on winter wheat in the Eastern United States. Plant Dis. 1998, 82, 64–68. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Cowger, C.; Parks, R.; Marshall, D. Appearance of powdery mildew of wheat caused by Blumeria graminis f. sp. tritici on Pm17-bearing cultivars in North Carolina. Plant Dis. 2009, 93, 1219. [Google Scholar] [CrossRef] [PubMed]
  7. Wang, Z.L.; Li, L.H.; He, Z.H.; Duan, X.Y.; Zhou, Y.L.; Chen, X.M.; Lillemo, M.; Singh, R.P.; Wang, H.; Xia, X.C. Seedling and adult plant resistance to powdery mildew in Chinese bread wheat cultivars and lines. Plant Dis. 2005, 89, 457–463. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. REX Consortium. Combining selective pressures to enhance the durability of disease resistance genes. Front. Plant Sci. 2016, 7, 1916. [Google Scholar] [CrossRef] [Green Version]
  9. Marone, D.; Russo, M.A.; Laido, G.; De Vita, P.; Papa, R.; Blanco, A.; Gadaleta, A.; Rubiales, D.; Mastrangelo, A.M. Genetic basis of qualitative and quantitative resistance to powdery mildew in wheat: From consensus regions to candidate genes. BMC Genom. 2013, 14, 562. [Google Scholar] [CrossRef] [Green Version]
  10. Yang, L.; Zhang, X.; Zhang, X.; Wang, J.; Luo, M.; Yang, M. Identification and evaluation of resistance to powdery mildew and yellow rust in a wheat mapping population. PLoS ONE 2017, 12, e0177905. [Google Scholar] [CrossRef] [Green Version]
  11. Kang, Y.; Zhou, M.; Merry, A.; Barry, K. Mechanisms of powdery mildew resistance of wheat—A review of molecular breeding. Plant Pathol. 2020, 69, 601–617. [Google Scholar] [CrossRef]
  12. Bonardi, V.; Cherkis, K.; Nishimura, M.T.; Dangl, J.L. A new eye on NLR proteins: Focused on clarity or diffused by complexity? Curr. Opin. Immunol. 2012, 24, 41–50. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Dangl, J.L.; Jones, J.D. Plant pathogens and integrated defence responses to infection. Nature 2001, 411, 826–833. [Google Scholar] [CrossRef] [PubMed]
  14. Jørgensen, J.H. Discovery, characterization and exploitation of Mlo powdery mildew resistance in barley. In Breeding for Disease Resistance. Developments in Plant Pathology; Johnson, R., Jellis, G.J., Eds.; Springer: Dordrecht, The Netherland, 1992; Volume 1. [Google Scholar] [CrossRef]
  15. Devoto, A.; Hartmann, H.A.; Piffanelli, P.; Elliott, C.; Simmons, C.; Taramino, G.; Goh, C.S.; Cohen, F.E.; Emerson, B.C.; Schulze-Lefert, P.; et al. Molecular phylogeny and evolution of the plant-specific seven-transmembrane MLO family. J. Mol. Evol. 2003, 56, 77–88. [Google Scholar] [CrossRef] [PubMed]
  16. Varallyay, E.; Giczey, G.; Burgyan, J. Virus-induced gene silencing of Mlo genes induces powdery mildew resistance in Triticum aestivum. Arch. Virol. 2012, 157, 1345–1350. [Google Scholar] [CrossRef] [PubMed]
  17. Marone, D.; Russo, M.A.; Laido, G.; De Leonardis, A.M.; Mastrangelo, A.M. Plant nucleotide binding site-leucine-rich repeat (NBS-LRR) genes: Active guardians in host defense responses. Int. J. Mol. Sci. 2013, 14, 7302–7326. [Google Scholar] [CrossRef] [Green Version]
  18. Brun, H.; Chevre, A.M.; Fitt, B.D.; Powers, S.; Besnard, A.L.; Ermel, M.; Huteau, V.; Marquer, B.; Eber, F.; Renard, M.; et al. Quantitative resistance increases the durability of qualitative resistance to Leptosphaeria maculans in Brassica napus. New Phytol. 2010, 185, 285–299. [Google Scholar] [CrossRef] [Green Version]
  19. Delourme, R.; Bousset, L.; Ermel, M.; Duffe, P.; Besnard, A.L.; Marquer, B.; Fudal, I.; Linglin, J.; Chadoeuf, J.; Brun, H. Quantitative resistance affects the speed of frequency increase but not the diversity of the virulence alleles overcoming a major resistance gene to Leptosphaeria maculans in oilseed rape. Infect. Genet. Evol. 2014, 27, 490–499. [Google Scholar] [CrossRef]
  20. Tanksley, S.D. Mapping polygenes. Annu. Rev. Genet. 1993, 27, 205–233. [Google Scholar] [CrossRef]
  21. Liu, W.; Maccaferri, M.; Chen, X.; Laghetti, G.; Pignone, D.; Pumphrey, M.; Tuberosa, R. Genome-wide association mapping reveals a rich genetic architecture of stripe rust resistance loci in emmer wheat (Triticum turgidum ssp. dicoccum). Theor. Appl. Genet. 2017, 130, 2249–2270. [Google Scholar] [CrossRef] [Green Version]
  22. Nordestgaard, N.V.; Thach, T.; Sarup, P.; Rodriguez-Algaba, J.; Andersen, J.R.; Hovmoller, M.S.; Jahoor, A.; Jorgensen, L.N.; Orabi, J. Multi-parental populations suitable for identifying sources of resistance to powdery mildew in winter wheat. Front. Plant Sci. 2020, 11, 570863. [Google Scholar] [CrossRef]
  23. Pogoda, M.; Liu, F.; Douchkov, D.; Djamei, A.; Reif, J.C.; Schweizer, P.; Schulthess, A.W. Identification of novel genetic factors underlying the host-pathogen interaction between barley (Hordeum vulgare L.) and powdery mildew (Blumeria graminis f. sp. hordei). PLoS ONE 2020, 15, e0235565. [Google Scholar] [CrossRef] [PubMed]
  24. Montilla-Bascon, G.; Rispail, N.; Sanchez-Martin, J.; Rubiales, D.; Mur, L.A.; Langdon, T.; Howarth, C.J.; Prats, E. Genome-wide association study for crown rust (Puccinia coronata f. sp. avenae) and powdery mildew (Blumeria graminis f. sp. avenae) resistance in an oat (Avena sativa) collection of commercial varieties and landraces. Front. Plant Sci. 2015, 6, 103. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Meuwissen, T.H.; Hayes, B.J.; Goddard, M.E. Prediction of total genetic value using genome-wide dense marker maps. Genetics 2001, 157, 1819–1829. [Google Scholar] [CrossRef] [PubMed]
  26. Desta, Z.A.; Ortiz, R. Genomic selection: Genome-wide prediction in plant improvement. Trends Plant Sci. 2014, 19, 592–601. [Google Scholar] [CrossRef]
  27. Asgarinia, P.; Cloutier, S.; Duguid, S.; Rashid, K.; Mirlohi, A.; Banik, M.; Saeidi, G. Mapping quantitative trait loci for powdery mildew resistance in flax (Linum usitatissimum L.). Crop Sci. 2013, 53, 2462–2472. [Google Scholar] [CrossRef] [Green Version]
  28. Aly, A. Examination of correlations between several biochemical components and powdery mildew resistance of flax cultivars. Korean Soc. Plant Pathol. 2012, 28, 149–155. [Google Scholar] [CrossRef] [Green Version]
  29. Rashid, K.Y. Powdery mildew on flax: A new disease in western Canada. Can. J. Plant Pathol. 1998, 20, 216. [Google Scholar]
  30. Beale, R.E. Studies of resistance in linseed cultivars to Oidium lini and Botrytis cinerea. In Production and Protection of Linseed; Aspects of Applied Biology Series; Association of Applied Biologists, Horticulture Research International, Wellesbourne: Warwick, UK, 1991; Volume 28, pp. 85–90. [Google Scholar]
  31. You, F.M.; Duguid, S.D.; Lam, I.; Cloutier, S.; Rashid, K.Y.; Booker, H. Pedigrees and genetic base of the flax varieties registered in Canada. Can. J. Plant Sci. 2016, 96, 837–852. [Google Scholar] [CrossRef] [Green Version]
  32. Rashid, K.; Duguid, S. Inheritance of resistance to powdery mildew in flax. Can. J. Plant Pathol. 2005, 27, 404–409. [Google Scholar] [CrossRef]
  33. Diederichsen, A.; Kusters, P.M.; Kessler, D.; Bainas, Z.; Gugel, R.K. Assembling a core collection from the flax world collection maintained by Plant Gene Resources of Canada. Genet. Resour. Crop Evol. 2013, 60, 1479–1485. [Google Scholar] [CrossRef]
  34. Soto-Cerda, B.J.; Diederichsen, A.; Ragupathy, R.; Cloutier, S. Genetic characterization of a core collection of flax (Linum usitatissimum L.) suitable for association mapping studies and evidence of divergent selection between fiber and linseed types. BMC Plant Biol. 2013, 13, 78. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. He, L.; Xiao, J.; Rashid, K.Y.; Yao, Z.; Li, P.; Jia, G.; Wang, X.; Cloutier, S.; You, F.M. Genome-wide association studies for pasmo resistance in flax (Linum usitatissimum L.). Front. Plant Sci. 2018, 9, 1982. [Google Scholar] [CrossRef] [PubMed]
  36. Lan, S.; Zheng, C.; Hauck, K.; McCausland, M.; Duguid, S.D.; Booker, H.M.; Cloutier, S.; You, F.M. Genomic prediction accuracy of seven breeding selection traits improved by QTL identification in flax. Int. J. Mol. Sci. 2020, 21, 1577. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. You, F.M.; Cloutier, S. Mapping quantitative trait loci onto chromosome-scale pseudomolecules in flax. Methods Protoc. 2020, 3, 28. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Sun, C.; Wang, B.; Wang, X.; Hu, K.; Li, K.; Li, Z.; Li, S.; Yan, L.; Guan, C.; Zhang, J.; et al. Genome-wide association study dissecting the genetic architecture underlying the branch angle trait in rapeseed (Brassica napus L.). Sci. Rep. 2016, 6, 33673. [Google Scholar] [CrossRef] [Green Version]
  39. Remington, D.L.; Thornsberry, J.M.; Matsuoka, Y.; Wilson, L.M.; Whitt, S.R.; Doebley, J.; Kresovich, S.; Goodman, M.M.; Buckler, E.S. Structure of linkage disequilibrium and phenotypic associations in the maize genome. Proc. Natl. Acad. Sci. USA 2001, 98, 11479–11484. [Google Scholar] [CrossRef] [Green Version]
  40. Kumar, S.; You, F.M.; Duguid, S.; Booker, H.; Rowland, G.; Cloutier, S. QTL for fatty acid composition and yield in linseed (Linum usitatissimum L.). Theor. Appl. Genet. 2015, 128, 965–984. [Google Scholar] [CrossRef]
  41. Li, P.; Quan, X.; Jia, G.; Xiao, J.; Cloutier, S.; You, F.M. RGAugury: A pipeline for genome-wide prediction of resistance gene analogs (RGAs) in plants. BMC Genom. 2016, 17, 852. [Google Scholar] [CrossRef] [Green Version]
  42. You, F.M.; Xiao, J.; Li, P.; Yao, Z.; Gao, J.; He, L.; Zhu, T.; Luo, M.-C.; Wang, X.; Deyholos, M.K.; et al. Chromosome-scale pseudomolecules refined by optical, physical, and genetic maps in flax. Plant J. 2018, 95, 371–384. [Google Scholar] [CrossRef]
  43. You, F.M.; Zheng, C.; Bartaula, S.; Khan, N.; Wang, J.; Cloutier, S. Genomic cross prediction for linseed improvement. In Accelerated Plant Breeding; Gosal, S.S., Wani, S.H., Eds.; Springer Nature: Berlin/Heidelberg, Germany, 2022; Volume 4, pp. 451–480. [Google Scholar]
  44. He, L.; Xiao, J.; Rashid, K.Y.; Jia, G.; Li, P.; Yao, Z.; Wang, X.; Cloutier, S.; You, F.M. Evaluation of genomic prediction for pasmo resistance in flax. Int. J. Mol. Sci. 2019, 20, 359. [Google Scholar] [CrossRef] [Green Version]
  45. You, F.M.; Booker, H.M.; Duguid, S.D.; Jia, G.; Cloutier, S. Accuracy of genomic selection in biparental populations of flax (Linum usitatissimum L.). Crop J. 2016, 4, 290–303. [Google Scholar] [CrossRef] [Green Version]
  46. Pessina, S.; Lenzi, L.; Perazzolli, M.; Campa, M.; Dalla Costa, L.; Urso, S.; Vale, G.; Salamini, F.; Velasco, R.; Malnoy, M. Knockdown of MLO genes reduces susceptibility to powdery mildew in grapevine. Hortic. Res. 2016, 3, 16016. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Acevedo-Garcia, J.; Spencer, D.; Thieron, H.; Reinstadler, A.; Hammond-Kosack, K.; Phillips, A.L.; Panstruga, R. mlo-based powdery mildew resistance in hexaploid bread wheat generated by a non-transgenic TILLING approach. Plant Biotechnol. J. 2017, 15, 367–378. [Google Scholar] [CrossRef]
  48. Buschges, R.; Hollricher, K.; Panstruga, R.; Simons, G.; Wolter, M.; Frijters, A.; van Daelen, R.; van der Lee, T.; Diergaarde, P.; Groenendijk, J.; et al. The barley Mlo gene: A novel control element of plant pathogen resistance. Cell 1997, 88, 695–705. [Google Scholar] [CrossRef] [Green Version]
  49. Wang, Y.; Cheng, X.; Shan, Q.; Zhang, Y.; Liu, J.; Gao, C.; Qiu, J.L. Simultaneous editing of three homoeoalleles in hexaploid bread wheat confers heritable resistance to powdery mildew. Nat. Biotechnol. 2014, 32, 947–951. [Google Scholar] [CrossRef] [PubMed]
  50. Nekrasov, V.; Wang, C.; Win, J.; Lanz, C.; Weigel, D.; Kamoun, S. Rapid generation of a transgene-free powdery mildew resistant tomato by genome deletion. Sci. Rep. 2017, 7, 482. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Koseoglou, E. The Study of SlPMR4 CRISPR/Cas9-Mediated Tomato Allelic Series for Resistance against Powdery Mildew; Wageningen University: Wageningen, The Netherlands, 2017. [Google Scholar]
  52. Huibers, R.P.; Loonen, A.E.; Gao, D.; Van den Ackerveken, G.; Visser, R.G.; Bai, Y. Powdery mildew resistance in tomato by impairment of SlPMR4 and SlDMR1. PLoS ONE 2013, 8, e67467. [Google Scholar] [CrossRef] [Green Version]
  53. Xiao, S.; Ellwood, S.; Calis, O.; Patrick, E.; Li, T.; Coleman, M.; Turner, J.G. Broad-spectrum mildew resistance in Arabidopsis thaliana mediated by RPW8. Science 2001, 291, 118–120. [Google Scholar] [CrossRef]
  54. Xiao, S.; Charoenwattana, P.; Holcombe, L.; Turner, J.G. The Arabidopsis genes RPW8.1 and RPW8.2 confer induced resistance to powdery mildew diseases in tobacco. Mol. Plant Microbe Interact. 2003, 16, 289–294. [Google Scholar] [CrossRef] [Green Version]
  55. Young, L.W.; Trouve, J.P.; Speck, A.; Kutcher, R.; Rashid, K.Y.; You, F.M.; Cloutier, S.; Booker, H.M. Mapping rust and powdery mildew resistant genes in flax. In Proceedings of the XXVII Plant and Animal Genome Conference, San Diego, CA, USA, 12–16 January 2019; p. PO0871. [Google Scholar]
  56. You, F.M.; Jia, G.; Xiao, J.; Duguid, S.D.; Rashid, K.Y.; Booker, H.M.; Cloutier, S. Genetic variability of 27 traits in a core collection of flax (Linum usitatissimum L.). Front. Plant Sci. 2017, 8, 1636. [Google Scholar] [CrossRef] [Green Version]
  57. Wang, Z.; Hobson, N.; Galindo, L.; Zhu, S.; Shi, D.; McDill, J.; Yang, L.; Hawkins, S.; Neutelings, G.; Datla, R.; et al. The genome of flax (Linum usitatissimum) assembled de novo from short shotgun sequence reads. Plant J. 2012, 72, 461–473. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Jo, H.; Koh, G. Faster single-end alignment generation utilizing multi-thread for BWA. Biomed. Mater. Eng. 2015, 26 (Suppl. 1), S1791–S1796. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R.; Genome Project Data Processing, Subgroup. The sequence alignment/map format and SAMtools. Bioinformatics 2009, 25, 2078–2079. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Kumar, S.; You, F.M.; Cloutier, S. Genome wide SNP discovery in flax through next generation sequencing of reduced representation libraries. BMC Genom. 2012, 13, 684. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. You, F.M.; Deal, K.R.; Wang, J.; Britton, M.T.; Fass, J.N.; Lin, D.; Dandekar, A.M.; Leslie, C.A.; Aradhya, M.; Luo, M.C.; et al. Genome-wide SNP discovery in walnut with an AGSNP pipeline updated for SNP discovery in allogamous organisms. BMC Genom. 2012, 13, 354. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. You, F.M.; Huo, N.; Deal, K.R.; Gu, Y.Q.; Luo, M.C.; McGuire, P.E.; Dvorak, J.; Anderson, O.D. Annotation-based genome-wide SNP discovery in the large and complex Aegilops tauschii genome using next-generation sequencing without a reference genome sequence. BMC Genom. 2011, 12, 59. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Huang, X.; Wei, X.; Sang, T.; Zhao, Q.; Feng, Q.; Zhao, Y.; Li, C.; Zhu, C.; Lu, T.; Zhang, Z.; et al. Genome-wide association studies of 14 agronomic traits in rice landraces. Nat. Genet. 2010, 42, 961–967. [Google Scholar] [CrossRef]
  64. Altshuler, D.; Donnelly, P.A.; International HapMap Consortium. A haplotype map of the human genome. Nature 2005, 437, 1299–1320. [Google Scholar] [CrossRef] [Green Version]
  65. Browning, S.R.; Browning, B.L. Rapid and accurate haplotype phasing and missing-data inference for whole-genome association studies by use of localized haplotype clustering. Am. J. Hum. Genet. 2007, 81, 1084–1097. [Google Scholar] [CrossRef] [Green Version]
  66. Price, A.L.; Patterson, N.J.; Plenge, R.M.; Weinblatt, M.E.; Shadick, N.A.; Reich, D. Principal components analysis corrects for stratification in genome-wide association studies. Nat. Genet. 2006, 38, 904–909. [Google Scholar] [CrossRef]
  67. Yu, J.; Pressoir, G.; Briggs, W.H.; Vroh Bi, I.; Yamasaki, M.; Doebley, J.F.; McMullen, M.D.; Gaut, B.S.; Nielsen, D.M.; Holland, J.B.; et al. A unified mixed-model method for association mapping that accounts for multiple levels of relatedness. Nat. Genet. 2006, 38, 203–208. [Google Scholar] [CrossRef] [PubMed]
  68. Liu, X.; Huang, M.; Fan, B.; Buckler, E.S.; Zhang, Z. Iterative usage of fixed and random effect models for powerful and efficient genome-wide association studies. PLoS Genet. 2016, 12, e1005767. [Google Scholar] [CrossRef] [PubMed]
  69. Zhang, J.; Feng, J.Y.; Ni, Y.L.; Wen, Y.J.; Niu, Y.; Tamba, C.L.; Yue, C.; Song, Q.; Zhang, Y.M. pLARmEB: Integration of least angle regression with empirical Bayes for multilocus genome-wide association studies. Heredity 2017, 118, 517–524. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Ren, W.L.; Wen, Y.J.; Dunwell, J.M.; Zhang, Y.M. pKWmEB: Integration of Kruskal-Wallis test with empirical Bayes under polygenic background control for multi-locus genome-wide association study. Heredity 2017, 120, 208–218. [Google Scholar] [CrossRef] [PubMed]
  71. Zhang, Y.M.; Tamba, C.L. A fast mrMLM algorithm for multi-locus genome-wide association studies. bioRxiv 2018, 341784. [Google Scholar] [CrossRef]
  72. Wen, Y.J.; Zhang, H.; Ni, Y.L.; Huang, B.; Zhang, J.; Feng, J.Y.; Wang, S.B.; Dunwell, J.M.; Zhang, Y.M.; Wu, R. Methodological implementation of mixed linear models in multi-locus genome-wide association studies. Brief. Bioinform. 2018, 19, 700–712. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Tamba, C.L.; Ni, Y.L.; Zhang, Y.M. Iterative sure independence screening EM-Bayesian LASSO algorithm for multi-locus genome-wide association studies. PLoS Comput. Biol. 2017, 13, e1005357. [Google Scholar] [CrossRef]
  74. Zhang, Y.-W.; Tamba, C.L.; Wen, Y.-J.; Li, P.; Ren, W.-L.; Ni, Y.-L.; Gao, J.; Zhang, Y.-M. mrMLM v4.0: An R platform for multi-locus genome-wide association studies. bioRxiv 2020, 18, 481–487. [Google Scholar] [CrossRef]
  75. Purcell, S.; Neale, B.; Todd-Brown, K.; Thomas, L.; Ferreira, M.A.; Bender, D.; Maller, J.; Sklar, P.; de Bakker, P.I.; Daly, M.J.; et al. PLINK: A tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 2007, 81, 559–575. [Google Scholar] [CrossRef] [Green Version]
  76. Endelman, J.B. Ridge regression and other kernels for genomic selection with R package rrBLUP. Plant Genome 2011, 4, 250–255. [Google Scholar] [CrossRef] [Green Version]
  77. VanRaden, P.M. Efficient methods to compute genomic predictions. J. Dairy Sci. 2008, 91, 4414–4423. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Tibshirani, R. Regression shrinkage and selection via the Lasso. J. R. Stat. Soc. Ser. B Methodol. 1996, 58, 267–288. [Google Scholar] [CrossRef]
  79. Pérez, P.; de los Campos, G. Genome-wide regression and prediction with the BGLR statistical package. Genetics 2014, 198, 483–495. [Google Scholar] [CrossRef] [PubMed]
  80. Pérez, P.; de Los Campos, G.; Crossa, J.; Gianola, D. Genomic-enabled prediction based on molecular markers and pedigree using the bayesian linear regression package in R. Plant Genome 2010, 3, 106–116. [Google Scholar] [CrossRef]
  81. Wu, Q.; Ye, Y.; Liu, Y.; Ng, M.K. SNP selection and classification of genome-wide SNP data using stratified sampling random forests. IEEE Trans. Nanobiosci. 2012, 11, 216–227. [Google Scholar] [CrossRef]
  82. Covarrubias-Pazaran, G. Genome-assisted prediction of quantitative traits using the R package sommer. PLoS ONE 2016, 11, e0156744. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Powdery mildew (PM) field evaluation over five years (2012–2016). Box plots of PM ratings for the 372 accessions of the core collection and the 75 selected breeding lines (A) and the number of accessions and mean PM ratings by resistance groups (B). R: resistant (PM ratings of 1 to <3); MR: moderately resistant (3 to <4); MS: moderately susceptible (4 to <5); S: susceptible (5 to <7); HS: highly susceptible (7 to 9).
Figure 1. Powdery mildew (PM) field evaluation over five years (2012–2016). Box plots of PM ratings for the 372 accessions of the core collection and the 75 selected breeding lines (A) and the number of accessions and mean PM ratings by resistance groups (B). R: resistant (PM ratings of 1 to <3); MR: moderately resistant (3 to <4); MS: moderately susceptible (4 to <5); S: susceptible (5 to <7); HS: highly susceptible (7 to 9).
Ijms 23 04960 g001
Figure 2. Distribution and correlation matrix of powdery mildew (PM) ratings in five consecutive years (2012–2016) and mean PM ratings over years. *** represents statistical significance at 0.001 probability level.
Figure 2. Distribution and correlation matrix of powdery mildew (PM) ratings in five consecutive years (2012–2016) and mean PM ratings over years. *** represents statistical significance at 0.001 probability level.
Ijms 23 04960 g002
Figure 3. Circos map of powdery mildew (PM) resistance quantitative trait nucleotides (QTNs) and their co-located candidate genes. Track A: 15 chromosomes of the flax genome, B: 1444 resistance gene analogues (RGAs), C: 445 RGA candidates with significant SNP allele effects for PM resistance, D: 45 RGA candidates containing PM resistance QTNs, E: 9 QTNs identified by GLM, F: 15 QTNs identified by FarmCPU, G: 349 unique QTNs identified by the six multi-locus models of the mrMLM package and another multi-locus model FarmCPU, H: 68 QTNs identified by FASTmrEMMA, I: 105 QTNs identified by FASTmrMLM, J: 83 QTNs identified by ISIS_EM_BLASSO, K: 72 QTNs identified by mrMLM, L: 90 QTNs identified by pKWmEB, M: 108 QTNs identified by mrMLM_pLARmEB. The central region links the duplicated RGA candidates for the 349 QTNs identified. Allele effects of QTNs (R2) are represented by different colours as indicated in the figure legend. QTN clusters are shaded in yellow across the tracks.
Figure 3. Circos map of powdery mildew (PM) resistance quantitative trait nucleotides (QTNs) and their co-located candidate genes. Track A: 15 chromosomes of the flax genome, B: 1444 resistance gene analogues (RGAs), C: 445 RGA candidates with significant SNP allele effects for PM resistance, D: 45 RGA candidates containing PM resistance QTNs, E: 9 QTNs identified by GLM, F: 15 QTNs identified by FarmCPU, G: 349 unique QTNs identified by the six multi-locus models of the mrMLM package and another multi-locus model FarmCPU, H: 68 QTNs identified by FASTmrEMMA, I: 105 QTNs identified by FASTmrMLM, J: 83 QTNs identified by ISIS_EM_BLASSO, K: 72 QTNs identified by mrMLM, L: 90 QTNs identified by pKWmEB, M: 108 QTNs identified by mrMLM_pLARmEB. The central region links the duplicated RGA candidates for the 349 QTNs identified. Allele effects of QTNs (R2) are represented by different colours as indicated in the figure legend. QTN clusters are shaded in yellow across the tracks.
Ijms 23 04960 g003
Figure 4. Relationship between R2, the coefficients of variation (CV, %) of quantitative trait nucleotide (QTN) effects and the number of powdery mildew (PM) datasets that showed significant QTN effects.
Figure 4. Relationship between R2, the coefficients of variation (CV, %) of quantitative trait nucleotide (QTN) effects and the number of powdery mildew (PM) datasets that showed significant QTN effects.
Ijms 23 04960 g004
Figure 5. Relationships between the number of favourable alleles and powdery mildew (PM) ratings (A) and the PM resistance groups (B). The average number of favourable alleles ( x ¯ ) per resistance group is indicated above the graph. A quadratic polynomial regression line is fitted in (A) with a confidence interval in grey colour.
Figure 5. Relationships between the number of favourable alleles and powdery mildew (PM) ratings (A) and the PM resistance groups (B). The average number of favourable alleles ( x ¯ ) per resistance group is indicated above the graph. A quadratic polynomial regression line is fitted in (A) with a confidence interval in grey colour.
Ijms 23 04960 g005
Figure 6. Quantitative trait nucleotide (QTN) effects of some large-effect and stable QTNs (R2 > 10%) showing QTNs with favourable allele frequencies (FAF) < 0.45 (A), QTNs with FAF > 0.55 (B) and QTNs with FAF ranging from 0.45 to 0.55 (C). The blue triangles indicate the favourable alleles. The QTN effects (R2) were calculated based on the PM-Mean dataset and are indicated above each box plot along with the actual FAF of the QTNs.
Figure 6. Quantitative trait nucleotide (QTN) effects of some large-effect and stable QTNs (R2 > 10%) showing QTNs with favourable allele frequencies (FAF) < 0.45 (A), QTNs with FAF > 0.55 (B) and QTNs with FAF ranging from 0.45 to 0.55 (C). The blue triangles indicate the favourable alleles. The QTN effects (R2) were calculated based on the PM-Mean dataset and are indicated above each box plot along with the actual FAF of the QTNs.
Ijms 23 04960 g006
Figure 7. Relationship between the favourable allele frequency (FAF) in all accessions and that in the 85 resistant accession subset (R, powdery mildew ratings < 3) based on their QTN effects (R2) as colour-coded in the legend.
Figure 7. Relationship between the favourable allele frequency (FAF) in all accessions and that in the 85 resistant accession subset (R, powdery mildew ratings < 3) based on their QTN effects (R2) as colour-coded in the legend.
Ijms 23 04960 g007
Figure 8. Relationships between flax morphotypes and powdery mildew (PM) ratings (A), and between flax morphotypes and the number of favourable alleles (B).
Figure 8. Relationships between flax morphotypes and powdery mildew (PM) ratings (A), and between flax morphotypes and the number of favourable alleles (B).
Ijms 23 04960 g008
Figure 9. RPW8 genes are potential candidates for four quantitative trait nucleotides (QTNs): Lu13-4791823 (R2 = 13.79%), Lu13-4830850 (R2 = 7.17%), Lu13-4866704 (R2 = 12.24%), and Lu13-4900476 (R2 = 10.9%) on chromosome 13.
Figure 9. RPW8 genes are potential candidates for four quantitative trait nucleotides (QTNs): Lu13-4791823 (R2 = 13.79%), Lu13-4830850 (R2 = 7.17%), Lu13-4866704 (R2 = 12.24%), and Lu13-4900476 (R2 = 10.9%) on chromosome 13.
Ijms 23 04960 g009
Table 1. Powdery mildew (PM) disease ratings for 447 flax genotypes (372 accessions from a core collection and 75 selected breeding lines) from 2012 to 2016.
Table 1. Powdery mildew (PM) disease ratings for 447 flax genotypes (372 accessions from a core collection and 75 selected breeding lines) from 2012 to 2016.
YearDataset CodeSample Size x ¯   ±   s RangeCV (%)
2012PM-20124474.9 ± 2.01.0–9.041.1
2013PM-20132535.2 ± 2.11.0–9.040.9
2014PM-20144474.1 ± 2.21.0–9.053.2
2015PM-20153913.6 ± 2.01.0–9.057.3
2016PM-20164474.9 ± 1.42.0–9.028.1
MeanPM-Mean4474.5 ± 1.41.6–8.031.1
x ¯ : population mean; s: standard deviation; CV: coefficient of variation. Due to missing data, only 253 and 391 accessions were available in 2013 and 2015, respectively.
Table 2. Quantitative trait nucleotides (QTNs) identified on 45 resistance gene analogues (RGAs) for flax powdery mildew (PM) resistance.
Table 2. Quantitative trait nucleotides (QTNs) identified on 45 resistance gene analogues (RGAs) for flax powdery mildew (PM) resistance.
GeneChrGene Start PositionGene End PositionGene FamilyTag QTN-PositionSNPFAFAFCV of QTN EffectsEffectR2
Lus1000605612868172728684082RLKLu1-28683876A/GA0.3927.12−0.414.44
Lus1003058722395541123958647RLKLu2-23956609A/CA0.2238.27−0.355.68
Lus10040576357460005750240TNLLu3-5748445G/AA0.78208.730.181.19
Lus1003689141242903512433792WRKYLu4-12432479G/TG0.3728.39−0.7011.51
Lus1004186041621294216215322RLKLu4-16213043A/GA0.1827.34−0.623.4
Lus10004727515342181535440TNLLu5-1534998A/GA0.2831.31−1.3327.83
Lus10004726515355021538672TNLLu5-1535619C/TC0.2627.43−1.5432.6
Lus10004719515687341573096TNLLu5-1569098G/TG0.3133.14−1.0627.63
Lus10032303530051883006765WRKYLu5-3006723T/CT0.4522.08−0.6824.08
Lus10032310530525663053123DIRLu5-3052714T/GG0.8472.850.643.12
Lus10032351532231373225286RLKLu5-3224350A/GA0.3939.99−0.6410.62
Lus10034795546439964646254RLKLu5-4646212C/TC0.4951.12−0.499.31
Lus1002986051327043013273781TNLLu5-13271207A/GA0.2228.62−1.7634.99
Lus10017649618786701885283RLKLu6-1883039A/GA0.0536.53−0.5511.79
Lus1002100161537694515378903RLPLu6-15378264G/AA0.78104.210.230.8
Lus1002521661666452316666663WRKYLu6-16666521C/TT0.5852.260.121.13
Lus1002319971735735517360952TM-CCLu7-17359522G/AG0.4461.38−0.150.73
Lus1000397171765700617660870TM-CCLu7-17659649A/GG0.9126.371.124.85
Lus10021849873863687423024RLKLu8-7394085A/GG0.7839.930.383.65
Lus1002226581599907916001545RLKLu8-15999956A/GG0.8727.700.471.36
Lus1000781281835070018354799TNLLu8-18351964C/TT0.6544.470.7311.49
Lus1000224981903748719042009TNLLu8-19040276G/AG0.1928.07−1.315.67
Lus100102219911825917007TNLLu9-916748T/GG0.9662.090.981.71
Lus10031043962659806269593RLKLu9-6266682T/CC0.74131.110.231.47
Lus1000677212719787720446RLKLu12-720013A/GG0.8061.750.52.61
Lus1000673212890823894178TNLLu12-892762C/GC0.1148.07−0.471.46
Lus100233231218938621897997RLKLu12-1896717A/TA0.1525.63−0.481.66
Lus100182891251103535112941TM-CCLu12-5111993G/CG0.4433.29−0.252.03
Lus10027903121661443116619009RLPLu12-16614785A/GG0.5353.260.2720.28
Lus10033608121912211919129113RLKLu12-19127670G/AA0.7793.090.342.78
Lus100013361317471601750173MLOLu13-1749576G/AG0.1441.08−0.441.23
Lus100197081345203774539065TNLLu13-4531367T/CT0.3929.66−0.7314.08
Lus100093641351416815142476RLKLu13-5142458C/TC0.1861.62−0.628.14
Lus10030845131821179418215062RLKLu13-18212664T/GG0.9252.020.884.89
Lus100286391411713451174215CNLLu14-1171479C/TC0.5336.69−0.4910.89
Lus100205341434571753462957TNLLu14-3458382C/GC0.5423.50−0.407.47
Lus100214481440188364024903TM-CCLu14-4021471A/TT0.9478.620.753.75
Lus100141501453754035382647RLKLu14-5382091A/GG0.9291.220.588.75
Lus100156481459556895959658TNLLu14-5959395G/CC0.8858.540.816.82
Lus100156491459600045963280TNLLu14-5960489G/AA0.7456.870.475.75
Lus10008320141054778710552007TNLLu14-10551333T/CC0.8045.780.755.78
Lus10035674141535723415366290TNLLu14-15360622T/CC0.5338.400.253.03
Lus10039211141720324817203932TNLLu14-17203266G/AG0.2442.02−1.0514.04
Lus10007610154790750779RLKLu15-50397G/AG0.4535.62−0.5215.63
Lus100126781539905884001681WRKYLu15-3991048T/GG0.8763.360.7518.51
Chr: chromosome; FA: favourable allele; FAF: favourable allele frequency; CV: coefficient of variation; RLK: receptor-like protein kinase; RLP: receptor-like protein; TM-CC: transmembrane coiled-coil protein; NBS: nucleotide-binding site domain; LRR: leucine-rich repeat; TIR: Toll/interleukin-1 receptor-like domain; CNL: CC–NBS–LRR; TNL: TIR-NBS-LRRs; TN: TIR–NBS; TX: TIR–unknown; MLO: mildew resistance locus o.
Table 3. Counts of candidate resistance-related genes based on gene families and R2 of QTNs. Only candidate genes with SNPs of significant allele effects on powdery mildew (PM) ratings are included.
Table 3. Counts of candidate resistance-related genes based on gene families and R2 of QTNs. Only candidate genes with SNPs of significant allele effects on powdery mildew (PM) ratings are included.
Candidate Gene FamilyR2 of QTNs (%)Total
<11–55–1010–20>20
DIR14 (1) 1 6 (1)
DZC 1 1
EDR 1 1
MLO 5 (1)1 6 (1)
RLK2665 (9)32 (4)13 (3)2138 (16)
RLP48 (1)622 (1)22 (2)
TM-CC418 (4)32 27 (4)
WRKY28 (1)43 (2)1 (1)18 (4)
CNL23 4 (1)110 (1)
TNL10 (1)25 (5)15 (2)7 (4)7 (4)64 (16)
TIR 1 1
Total49138 (21)62 (7)32 (10)13 (6)294 (45)
The numbers in parentheses are the number of candidate genes within which QTNs were identified. RLK: receptor-like protein kinase; RLP: receptor-like protein; TM-CC: transmembrane coiled-coil protein; NBS: nucleotide-binding site domain; LRR: leucine-rich repeat; TIR: Toll/interleukin-1 receptor-like domain; CNL: CC–NBS–LRR; TNL: TIR-NBS-LRR; TN: TIR–NBS; TX: TIR–unknown; MLO: mildew resistance locus o; EDR: extreme-drug-resistant; DZC: disease resistance, zinc finger, chromosome condensation; DIR: dirigent protein.
Table 4. Genomic heritability (h2) and genomic predictive ability (r) of five datasets for powdery mildew (PM) resistance. Predictive ability was estimated based on the PM-Mean dataset and the GBLUP model.
Table 4. Genomic heritability (h2) and genomic predictive ability (r) of five datasets for powdery mildew (PM) resistance. Predictive ability was estimated based on the PM-Mean dataset and the GBLUP model.
Marker Seth2r
349 QTNs (identified from 247,160 SNPs)0.691 ± 0.0400.925 ± 0.014 a
388 QTNs (identified from both 247,160 SNPs and 3230 RGA-derived SNPs)0.698 ± 0.0410.926 ± 0.013 a
132 QTNs located within genes (RGAs and non-RGAs)0.481 ± 0.0480.850 ± 0.024 b
279 QTNs/SNPs located within candidate RGAs0.587 ± 0.050 0.822 ± 0.031 c
247,160 SNPs0.600 ± 0.0700.690 ± 0.049 d
Letters on the right of the predictive ability (r) represent statistical significance at a 5% probability level using the Tukey HSD test.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

You, F.M.; Rashid, K.Y.; Zheng, C.; Khan, N.; Li, P.; Xiao, J.; He, L.; Yao, Z.; Cloutier, S. Insights into the Genetic Architecture and Genomic Prediction of Powdery Mildew Resistance in Flax (Linum usitatissimum L.). Int. J. Mol. Sci. 2022, 23, 4960. https://doi.org/10.3390/ijms23094960

AMA Style

You FM, Rashid KY, Zheng C, Khan N, Li P, Xiao J, He L, Yao Z, Cloutier S. Insights into the Genetic Architecture and Genomic Prediction of Powdery Mildew Resistance in Flax (Linum usitatissimum L.). International Journal of Molecular Sciences. 2022; 23(9):4960. https://doi.org/10.3390/ijms23094960

Chicago/Turabian Style

You, Frank M., Khalid Y. Rashid, Chunfang Zheng, Nadeem Khan, Pingchuan Li, Jin Xiao, Liqiang He, Zhen Yao, and Sylvie Cloutier. 2022. "Insights into the Genetic Architecture and Genomic Prediction of Powdery Mildew Resistance in Flax (Linum usitatissimum L.)" International Journal of Molecular Sciences 23, no. 9: 4960. https://doi.org/10.3390/ijms23094960

APA Style

You, F. M., Rashid, K. Y., Zheng, C., Khan, N., Li, P., Xiao, J., He, L., Yao, Z., & Cloutier, S. (2022). Insights into the Genetic Architecture and Genomic Prediction of Powdery Mildew Resistance in Flax (Linum usitatissimum L.). International Journal of Molecular Sciences, 23(9), 4960. https://doi.org/10.3390/ijms23094960

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop