Skip to main content

Homoeologous duplicated regions are involved in quantitative resistance of Brassica napus to stem canker

Abstract

Background

Several major crop species are current or ancient polyploids. To better describe the genetic factors controlling traits of agronomic interest (QTL), it is necessary to understand the structural and functional organisation of these QTL regions in relation to genome duplication. We investigated quantitative resistance to the fungal disease stem canker in Brassica napus, a highly duplicated amphidiploid species, to assess the proportion of resistance QTL located at duplicated positions.

Results

Genome-wide association analysis on a panel of 116 oilseed rape varieties genotyped with 3228 SNP indicated that 321 markers, corresponding to 64 genomic regions, are associated with resistance to stem canker. These genomic regions are relatively equally distributed on the A (53%) and C (47%) genomes of B. napus. Overall, 44% of these regions (28/64) are duplicated homoeologous regions. They are located in duplications of six (E, J, R, T, U and W) of the 24 ancestral blocks that constitute the B. napus genome. Overall, these six ancestral blocks have 34 duplicated copies in the B.napus genome. Almost all of the duplicated copies (82% of the 34 regions) harboured resistance associated markers for stem canker resistance, which suggests structural and functional conservation of genetic factors involved in this trait in B. napus.

Conclusions

Our study provides information on the involvement of duplicated loci in the control of stem canker resistance in B. napus. Further investigation of the similarity/divergence in sequence and gene content of these duplicated regions will provide insight into the conservation and allelic diversity of the underlying genes.

Background

Polyploidy or whole genome duplication (WGD) is an important phenomenon that has occurred during speciation and diversification of most plant species [1, 2]. It is estimated that 70% of angiosperms are polyploid [1, 2]. Polyploidy results from either the duplication of a same genome (autopolyploidy) or the interspecific hybridisation of genomes of two related species (allopolyploidy). The duplicated regions undergo important changes that can cause structural, functional and/or regulatory modifications to the duplicated genes. Various reports have suggested mechanisms to explain the evolution of duplicated genes including the gene balance hypothesis, diploidization, and neo- / sub- functionalization [14]. It is expected that these evolutionary processes might play a fundamental role in the diversification of the genes underlying complex traits. Various studies, mainly in allopolyploid plant species, have highlighted the involvement of duplicated loci (genes or QTL, for Quantitative Trait Loci), located at homoeologous positions, in the control of agronomic traits. Examples were reported in hexaploid wheat for a large number of traits (flowering time, glutenin synthesis and resistance to cyst nematodes) [57] and in soybean for flowering time [8]. Such duplicated homoeologous QTL may often represent a large proportion of the genetic factors controlling a complex trait. In cotton, a meta-analysis of the QTL involved in lint fibre quality showed that 21% of the QTL controlling the trait are located at homoeologous positions [9]. Another recent example of genetic analysis of fruit quality in strawberry revealed that 23% of the QTL detected are homoeologous [10].

Brassica napus is a suitable model for studying the effects of WGD on genetic factors involved in the control of complex traits. B. napus (2n = 4x = 38, genome: AACC) is an allotetraploid species formed from the hybridization between B. rapa (2n = 2x = 20, A genome) and B. oleracea (2n = 2x = 18, C genome) [11]. Brassica ancestors have undergone two duplication events (α and β) [12, 13] and two triplication events, one ancient event shared by a large majority of Angiosperms (γ) and a more recent event specific to the Brassica clade [14]. These WGD events, along with the merger of the two parental genomes, have resulted in a large number of duplicated regions in the B. napus genome. Information about the homoeologous relationships between the two genomes of B. napus is available. Indeed, Parkin et al. [15] identified regions of homoeology through the analysis of intra- and intergenomic duplications in the A and C genomes of B. napus. Moreover, the structural organisation of these duplicated genomic regions was studied by taking advantage of the close relatedness of Brassica sp. and the model plant Arabidopsis thaliana [16, 17]. Twenty four conserved blocks of colinearity were identified between the Arabidopsis and Brassica genomes [17, 18] in comparative mapping studies which showed that these blocks were highly duplicated in the B. napus genome [17, 19]. These duplications are regions that are either orthologous between (located on homoeologous linkage groups) or paralogous within the A and C genomes. As the B. rapa and B. oleracea genomes contain three subgenomes resulting from a meso-triplication [14], paralogous regions within the B. napus A and C genomes could also correspond to ancient homoeologous regions in the B. rapa and B. oleracea genomes. For this reason, in the present study all these duplicated regions will be referred to as homoeologous/duplicated regions.

Several genetic analyses for various complex agronomic traits such as oil content, seed yield, flowering time [20, 21] and disease resistance [22, 23] have been carried out in B. napus, mostly using linkage mapping approaches. Few studies reported duplicated QTL. For example, the involvement of duplicated QTL at homoeologous positions was described for seed glucosinolates [24], flowering time [25], yield-related traits [26] and resistance to sclerotinia stem rot [27]. Due to advances in genome sequencing, comparative mapping and computational technologies, it is now possible to further characterise the duplicated regions involved in the control of complex traits by more precise mapping of QTL and better estimate the proportion of QTL at homoeologous positions. In the present study, we focused on the quantitative resistance to stem canker, caused by the fungal pathogen Leptosphaeria maculans, which is one of the major diseases of B. napus crops worldwide [28, 29]. Linkage mapping approaches have been used to identify resistance QTL in B. napus segregating populations [3033] and to date, more than twenty QTL have been reported [33, 34]. In our laboratory, QTL detected with linkage mapping studies in bi-parental populations and in one connected population were analysed in relation to the homoeology information provided by Parkin et al. [15, 16]. This preliminary analysis suggested that several QTL for resistance to stem canker are localised in homoeologous regions, in particular on linkage groups A1/C1, A2/C2 and A3/C3 (unpublished data). The fact that some QTL are located at homoeologous positions can reflect functional redundancy or be involved in increased allelic diversity of the genes controlling the trait [1, 2, 35, 36]. Knowledge of the resulting diversity in the regions involved in stem canker resistance would help to construct resistant varieties with improved durability with the hypothesis that increasing the diversity of genetic factors controlling the resistance would result in an increase of the potential durability of the resistance [34]. A first step would be to increase the accuracy of detection of homoeologous duplicated QTL using new methods and precise data on the location of duplicated regions in the B. napus genome.

The aim of the present study was to assess the proportion of resistance QTL to stem canker located at homoeologous duplicated positions in the B. napus genome. To identify genetic regions involved in the resistance, we carried out a genome-wide association analysis. This approach allows the larger genetic diversity present in collections of varieties, compared to bi-parental or connected populations derived from a few parental lines, to be taken into account. Thus QTL detection is more exhaustive than in linkage analysis. The genomic position of the resistance-associated markers was then analysed in relation to the structural organisation of the duplicated regions in the B. napus genome in order to estimate the proportion of homoeologous duplicated QTL.

Results

Single Nucleotide Polymorphism marker analysis and genetic diversity

The panel of oilseed rape (OSR) varieties was first genotyped with 4329 single nucleotide polymorphism (SNP) markers. Rare alleles were eliminated and a set of 3228 SNP markers with a major allele frequency (MajAF) less than 0.95 was retained. These 3228 markers were well distributed on the B. napus linkage map with an average of one SNP every 0.62 centimorgan (cM). A total of 1986 SNP markers was located on the A genome (1 SNP each 0.50 cM) and 1242 SNPs on the C genome (1 SNP each 0.83 cM). The mean polymorphism information content (PIC) value was 0.27 and 0.28 on the A and C genomes, respectively. Out of the 3228 SNP markers, more than 87% showed less than 5% of heterozygous genotyping data and 0.4% had 10 to 13% of heterozygous genotyping data.

Linkage disequilibrium

The linkage disequilibrium (LD) was evaluated by calculating the r2 coefficient. Genotyping data of the 3228 SNP markers with a MajAF < 0.95 were used to calculate the r2 coefficient for the entire panel of 116 OSR varieties. Pairs of markers located on the same linkage group (LG) were defined as linked markers and those located on different LGs were defined as unlinked markers. The mean r2 was equal to 0.023 for the whole genome. For linked and unlinked markers, the mean r2 was equal to 0.078 and 0.020, respectively. A total of 8.62% and 0.3% of the linked and unlinked tested pairs of SNPs had a r2 coefficient superior to 0.2, respectively. The percentage of pairs of markers in LD (r2 > 0.2) was 7.77% on the A genome and of 9.57% on the C genome. However, the LD decayed a bit more rapidly on the C genome (up to 1.11 cM) compared with the A genome (up to 1.36 cM) (Additional file 1). On the whole genome, the LD extended up to 1.28 cM, with some variation depending on the LG and the region of the LG considered (Additional file 1, Additional file 2). Overall, markers in strong LD were mapped close to each others on the LGs. However, some exceptions were observed on LG A2, A8 and C8 (Additional file 3) on which LD extended further.

Population structure and kinship

A set of 727 SNP markers out of the 3228 available was used for Principal Component Analysis (PCA). These 727 SNPs were selected by eliminating close markers that were in strong LD to limit bias in the structure analysis. Results of the PCA (Additional file 4) did not show a strong structure within the population. The variance on the principal eigenvector was mainly due to the accession “Yudal” located at the extreme right of the vector, which is a spring type OSR variety. The two first principal components explained 13.9% of the variation in the panel. Sixteen significant axes that explained 54.7% of the variation within the panel were retained in the Tracy-Widom test. The matrix of coordinates of the accessions on these 16 axes was used as the P matrix for the association analysis.

The distribution of the kinship coefficients, calculated as the proportion of shared alleles between pairs of varieties, is shown in Additional file 5. The mean kinship coefficient value was 0.64. Up to 82% of the pairs of varieties tested had a kinship coefficient value gathered around the mean kinship value, i.e. between 0.55 and 0.70. Only a few coefficients showed extreme values. The matrix of kinship coefficients between pairs of varieties was used as the K matrix for the association analysis.

Marker-trait association analysis

Marker-trait association analysis was carried out with the 3228 SNPs (MajAF < 0.95) on the collection of 116 OSR varieties. Three linear models were used to test the marker/trait associations: (i) a General Linear Model (GLM), which does not take into account the relatedness between individuals, or the population structure; (ii) a K Compressed Mixed Linear Model (K CMLM), which only takes into account the relatedness between the varieties by using the K matrix; (iii) and a KP CMLM, which takes into account both the relatedness and the population structure, by using both the K and P matrices. A total of 1009 markers associated with resistance to stem canker was identified with the three models (type I error α = 5%). Of these 1009 markers, 192 (19%) were identified with at least two models, 689 (68%) were identified with the GLM only, 108 (11%) with the KP CMLM only and 20 (2%) with the K CMLM only (Figure 1). The markers identified with the GLM had overall lower p values than those identified with the K or KP CMLM. Indeed, the p value for approximately 50% of the markers identified with the GLM model was less than 0.01 whereas only 17% and 20% of the markers identified with the KP and K CMLM, respectively, had a p value less than 0.01. When a FDR test was applied at 0.25, no significant associations (α = 5%) were calculated with the K and KP models, whereas 300 (35% of the 855) associations remained significant with the GLM. Of these, 112 were common to the K and/or KP CML models and the remaining 188 markers were located close to markers associated with the K and or KP CML models.

Figure 1
figure 1

Number of resistance-associated markers identified with three different linear models. A panel of 116 winter oilseed rape varieties was used for marker-resistance association analysis. The General Linear Model is a simple linear model, the K Compressed Mixed Linear Model (CMLM) includes the relatedness between individuals of the panel and the KP CML model includes the relatedness and structure of the panel. The number in italics below the model name indicates the total number of markers identified with that model. Marker-trait associations were identified with a type I error of 0.05.

QQ plots of the observed p-values against the expected p-values under the null hypothesis obtained for the three models suggested that false-positive associations were more likely to be identified with the GLM than with the K and KP models (Figure 2). For the GLM (Figure 2A), deviation from the null hypothesis across the entire distribution indicated bias due to population stratification and/or relatedness [37, 38]. When the relatedness and structure in the panel were included in the model (Figure 2B and C), the deviation from the null hypothesis was corrected, indicating improved control of this type I error.

Figure 2
figure 2

Regression of the expected –log 10 (p value) on the observed –log 10 (p value). A Q-Q plot was used to compare the distribution of observed associations (blue dots) with statistics expected under the null hypothesis of no associations (represented by the red line) for the three linear models: a General Linear model (A), K Compressed Mixed Linear model (B) and a KP Compressed Mixed Linear model (C). The negative log10(p) was used instead of the p value so that the most significant markers are located at the top right corner of the graphs. A distribution close to the null hypothesis for most of the markers indicated a good control for false positives since it is expected that only a small number of tested markers are truly associated with the resistance.

Because of their capacity to control false-positive associations, in further analyses we focused on results from the K and KP models. The K and KP CML models identified a total of 321 markers (Additional file 6) significantly associated with resistance to stem canker (p-value < 0.05) (Figure 3). Eighty four markers (26%) were identified with both models, 120 (37%) with the K CMLM only and 117 (36%) with the KP CMLM only. The associated markers are located on all LGs, except for the LG C5 where no association was identified (Figure 3). Among these resistance associated markers, 217 markers are located on the A genome and 104 markers on the C genome. This represents 10.9% of the 1986 markers mapped to the A genome and 8.4% of the 1242 markers mapped to the C genome.

Figure 3
figure 3

Genome-wide association analysis of resistance to stem canker in a population of 116 oilseed rape cultivars. The K Compressed Mixed Linear model (A) takes into account the kinship between the varieties while the KP Compressed Mixed Linear model (B) takes into account the kinship and the structure of the varieties within the panel. Negative log10 of p-values were plotted against the genetic distance in cM on the 19 B. napus linkage groups (from 1 to 10: LGs A01 to A10; from 11 to 19: LGs C01 to C09). The horizontal line indicates the genome-wide association significance threshold.

Structural organisation of the markers associated with resistance to stem canker

The structural organisation of the resistance-associated markers was then analysed in relation to their location on the 24 conserved blocks of colinearity between B. napus and A. thaliana (henceforth named AK blocks) and to their anchorage on A. thaliana genes. The results are presented in Table 1 and in Additional file 7. Out of the 321 associated markers, 279 were assigned to one unique block, 25 to two or three different blocks and 17 could not be assigned (Additional file 6). The 279 markers with unique anchorage were located on 23 AK blocks corresponding to 64 distinct genomic regions (Additional file 7). For nine of these genomic regions, resistance associated markers were identified only on one copy of the corresponding nine blocks. For the other 55 genomic regions, resistance associated markers were identified on at least two copies of the 14 corresponding blocks. This suggests that more than 85% (55 out of 64) of the genomic regions associated with resistance to stem canker are duplicated regions. A more detailed analysis, based on the correspondence between resistance associated markers and A. thaliana genes, was carried out to investigate these 55 duplicated regions. In 28 genomic regions, corresponding to six of the 14 blocks, resistance-associated markers had significant hit against the same interval of the A.thaliana sequence suggesting that these are located in strictly duplicated regions (Table 2). In the other 27 regions corresponding to the eight other blocks, resistance-associated markers aligned with neighbouring intervals of the A. thaliana sequence suggesting that these resistance-associated are not strictly located in the same duplicated region (Additional file 7). We then examined more closely the resistance-associated markers in the six strictly duplicated regions on the E, J, R, T U and W blocks. Within the B. napus genome, the E block is duplicated four times, and the J, R, T, U and W blocks are duplicated six times (Figure 4). Resistance-associated markers were identified on the four duplicated E blocks, on the six duplicated J and U blocks, on two of the six duplicated T blocks and on five of the six duplicated R and W blocks. Thus, out of the 34 genomic regions corresponding to the duplication of the six investigated blocks, 28 (82%) carried resistance-associated markers. These 28 genomic regions grouped 69% of the 279 associated markers assigned to one unique AK block and 60% of the total 321 associated markers identified. Details of the organisation of the 28 regions with resistance associated markers on the B. napus genome are shown in Figure 4. Twenty-three of the 28 genomic regions were located on five collinear homoeologous LGs: A2/C2, A3/C3, A5/C4, A7/C6 and A10/C9. In four of the 28 genomic regions, QTL were identified on only one of the two homoeologous LGs. Finally, one of the 28 QTL (on block U, LG A8) was located in a region for which no homoeologous region has been identified in the B. napus genome to date. The number of resistance-associated markers detected and their location in the six duplicated blocks of interest are presented in Table 2.

Table 1 Distribution of genomic regions associated to stem canker resistance in relation to the organisation of duplicated blocks in the B. napus genome
Table 2 Resistance-associated markers detected in strictly duplicated regions
Figure 4
figure 4

Distribution of resistance-associated markers located in strictly duplicated regions of the B. napus genome. The coloured segments correspond to the six duplicated blocks, as defined by Schranz et al. [17], on which associated markers were identified in regions collinear to the same A. thaliana sequence interval. All the known duplications of these six blocks in the B. napus genome are shown. The number of resistance-associated markers identified on each block is indicated in parenthesis. The absence of a number indicates that no association was identified in the region. The resistance QTL identified in bi-parental or connected populations by linkage mapping on the copies of the 6 duplicated blocks are represented by stars on the right of the linkage groups. Orange stars indicate co-localisation of associated markers and previously identified QTL, black stars indicate loci where QTL were detected but where no resistance associated marker was identified. Linkage groups are organized according to the colinearity between the A and C genomes. The orientation of the arrows indicates the orientation of the colinearity between the genomic regions.

Discussion

This study highlighted that a large proportion of genomic regions involved in resistance to stem canker is located in duplicated homoeologous regions of the B. napus genome. These results were obtained using a genome-wide association analysis with a large number of SNPs on a panel of oilseed rape (OSR) varieties and by exploiting information available on the structural organisation of duplicated blocks in the B. napus genome.

Genome-wide association studies (GWAS) allow the identification of QTL in polyploid species

Advances in genome sequencing and computational technologies have led to high throughput SNP discovery in polyploid species [39] including B. napus [40]. Thus a large number of markers mapping at single loci are now available that can be used to conduct a precise GWAS as in diploid species. The power of association detection by GWAS partly depends on marker density and the extent of LD [41]. In the present study the average genome wide LD decayed within 1.28 cM, supporting previous studies in winter oilseed rape collections [40, 42, 43]. The average marker density (1 SNP every 0.62 cM) and the extent of LD in our panel allowed high-power detection of resistance-associated markers over the whole genome. Two LGs where the extent of LD was less than the marker density were the exception (Additional file 1). We limited the detection of false positive associations by taking into account the population structure and the varieties relatedness as previously reported [44]. Results from the QQ-plots showed that the K and KP models were better at limiting the false positive association rate than the GLM. This conclusion has been made in several other association mapping studies which also took into account the panel structure [44, 45], including that carried out by Jestin et al. [46] on the expanded OSR panel.

Some of the resistance-associated markers identified here are localised in genomic regions where we previously detected QTL by linkage mapping. A total of 28% of the associated markers identified with the K and KP CML models co-localised with 55% and 78% of the QTL detected in the biparental DYDH segregating population and in a connected multiparental population, respectively [30, 31, 47]. Some other regions were detected with only one or the other method. Differences between genomic regions identified by association mapping and linkage analysis have been observed in B. napus for other traits such as oil content [48] as well as in other species [8]. These differences are due to the type of population used, the allele frequency in the population, the type and number of markers used or the genetic map coverage. Our results show that a large proportion of QTL identified in the segregating populations can be detected in the association panel but also that new sources of resistance are present in the panel. This highlights the need to use both QTL detection methods in combination to exhaustively dissect the genetic architecture of complex traits for plant breeding.

In this study, we focused on the quantitative resistance to stem canker. Qualitative resistance controlled by specific major resistance genes have also been identified in B. napus [34]. These genes are Rlm2 on LG A10 and Rlm1 and the putative cluster Rlm3-Rlm4-Rlm7-Rlm9 on LG A07 [34]. To maximise chances to only study the quantitative resistance, the OSR varieties in our panel were chosen because they did not carry any effective specific resistance gene to L. maculans in our field conditions. Thus, we excluded OSR varieties known as carrying Rlm7, a highly effective gene that was recently introduced in OSR varieties in France. However, we cannot exclude that some associated markers might be linked to genes with a lower level of effectiveness. Indeed, associated markers located on the E bloc on linkage group A07 are in the region of the Rlm1 gene indicating that these markers could be linked to this specific resistance gene. In contrast, no associated markers were detected in the vicinity of the Rlm cluster region on A07.

Organisation of QTL for resistance to stem canker in the B. napusgenome

All of the resistance associated markers identified by GWAS corresponded to 64 genomic regions involved in the control of stem canker. These genomic regions are quite equally distributed on the A and C genomes of B. napus indicating that both genomes contribute to stem canker quantitative resistance. The number of genomic regions involved in resistance identified in this study is higher than that found by Jestin et al. [46] in the expanded panel of WOSR, mainly due to the use of a denser map. By exploiting the collinearity between the A. thaliana and B. napus genomes, we found that 28 (44%) of the 64 regions associated to resistance correspond to strictly duplicated genomic regions. Some other genomic associated regions correspond to neighbouring duplicated genomic regions of some of the AK blocks. If the accuracy of the GWAS analysis can be improved, we may find that some of these regions also correspond to strictly duplicated genomic regions. Thus it appears that more than 44% of the genomic regions involved in resistance to stem canker are strictly duplicated. This percentage of homoeologous genomic regions controlling a complex trait is higher than that reported so far in other polyploid species. Indeed, in allotetraploid cotton, Rong et al. [9] estimated that 21% of the 196 QTL involved in the control of cotton lint fibre quality are homoeologous QTL. In allooctoploid strawberry, Lerceteau-Kholer et al. [10], found that 23% of the 87 QTL involved in controlling strawberry fruit quality are homoeologous QTL. It is obvious that the proportion of homoeologous duplicated QTL that can be identified depends on the ability to exhaustively detect causative variants and the accuracy in identifying homoeologous regions (which requires a complete dissection of the genome structure). In addition, the proportion of homoeologous duplicated QTLs that can be observed also depends on the divergence time between the parental species and the evolutionary process (mainly artificial or natural selection) taking place in the history of these species for the studied trait. Thus it is difficult to conclude if there are really more duplicated genomic regions involved in the control of a quantitative trait in B. napus as compared with other polyploid species.

However, the high percentage of identified homoeologous genomic regions involved in resistance to stem canker in B. napus could be tightly linked to the significant structural conservation between B. rapa, B. oleracea and the A and C genomes of B. napus [49]. In addition, the two parental species have diverged more recently compared to other polyploid species such as cotton [14, 50]. Li et al. [51] recently identified QTL for morphological and yield component traits in duplicated regions within the A genome of B. rapa and in homoeologous regions between the A, B, and C genomes of B. rapa, B. juncea and B. napus, which supports the strong level of conservation between Brassica genomes.

Another factor affecting the relatively large proportion of duplicated homoeologous genomic regions involved in resistance to stem canker may be the nature of the trait. It is assumed that polyploidization events result in redundancy of duplicated genes but this gene redundancy can be reduced in a long term evolutionary process through various diploidization mechanisms among which genome fractionation is a major force [2, 3, 52]. Fractionation has been shown in polyploid species such as Arabidopsis [53], maize [54, 55] or B. rapa [56]. However, various studies highlighted that some functional gene categories are preferentially preserved from this reduction and are over-retained in duplicated copies. These include transcriptional factors, and protein kinases in the Arabidopsis lineage [57, 58], genes involved in signal transduction and some responses to external stimuli in rice [59], and genes involved in networks with a high level of connectivity in soybean [60]. In B. rapa, a recent study showed that genes involved in resistance to pathogens and especially in broad spectrum defence, were over-retained after the ancestral triplication event [61]. Given that the genomes of B. napus and B. rapa are highly conserved, it can be hypothesized that the large number of identified genomic regions involved in resistance to stem canker is tightly linked to this over-retention phenomenon. However, some genome fractionation or modification of gene expression/regulation could have occurred in some of the duplicated regions which would explain why, in some cases, resistance associated markers or QTL were not found on all the homoeologous copies in B. napus genome. Indeed, small-scale deletions have been observed at the genome microstructure level and in sequences in Brassica [49]. Examples of gene neofunctionalisation and subfunctionalisation have also been observed in B. rapa [6264]. These modifications were also extensively demonstrated in synthetic B. napus allotetraploids after the polyploidization process [65, 66] and in stabilized natural B. napus [49]. Nevertheless, the absence of resistance associated markers at some homoeologous loci could also be due to our experimental set-up. The causative variants at the “missing” QTL may be too rare, or even absent, in the panel to be detected through association mapping. Indeed, some of these “missing” QTL were identified in previous studies by linkage mapping, which is a more efficient method for detecting rare alleles [30, 46, 47].

It will now be interesting to further investigate the regions carrying the homoeologous QTL identified in this study. Characterisation of their gene content and expression, and sequence similarities will help elucidate the consequences of their conservation and evolution on the diversity of genetic factors involved in quantitative resistance to stem canker.

Conclusions

Our GWA study showed that many homoeologous duplicated regions are involved in the control of resistance of B. napus to stem canker. The resistance related genomic regions identified with GWAS overlapped with those identified by linkage mapping but this study also provided new information. Thus both methods can be valuable for a complete dissection of the architecture of genomic regions controlling complex traits. Our results showed that both the A and C genomes equally contribute to stem canker resistance and that 44% of the regions involved in the control of the resistance corresponded to homoeologous duplicated regions. In these regions, the associated markers were located in strictly duplicated intervals of six blocks which are syntenic with A. thaliana. Further studies to characterise the similarity/divergence in gene content and sequence of these duplicated regions are needed to gain insight into the conservation and allelic diversity of the underlying genes.

Methods

Plant material

A panel consisting of 115 European winter oilseed rape (OSR) varieties and one Asiatic spring OSR variety was used. This panel is a subset of that used in Jestin et al. [46]. Most of the varieties in the panel are double low, i.e. without erucic acid and with low glucosinolate seed content, and are registered in the French and the European catalogues (Additional file 8). These varieties were chosen because they did not carry any effective specific resistance gene to L. maculans in our field conditions. Thus, we excluded OSR varieties known as carrying Rlm7, a highly effective gene that was recently introduced in OSR varieties in France . Moreover, varieties of the panel showed a large range of responses to the infection, according to previous results obtained in different trials from different years [46].

Phenotypic evaluation

The material was previously evaluated in a field trial in 2006 at Le Rheu (France) [46]. For the phenotypic evaluation a disease index called the “G2 index” was calculated to classify the genotypes on a scale running from 0 to 9; from the most resistant (G2 index = 0) to the most susceptible varieties (G2 index = 9). The G2 index indicates the area of necrosis at the base stem section of the plant [44]. Statistical analysis of field trial results are presented in Jestin et al. [46].

Statistical analysis of genotyping data

The whole panel was genotyped with 5685 SNP markers. Within these markers, 4329 were mapped on an integrated map generated in our laboratory [40]. The map covers 2027.7 cM, which represents a marker density of 2.56 SNP per cM. A total of 2839 (65.6%) and 1490 (34.4%) of these SNP markers were mapped on the A and the C genome of B. napus, respectively. Only the genotyping data obtained with the mapped 4329 SNP markers were conserved for the following analysis.

The major allele frequency (MajAF), percentage of heterozygosity and polymorphic information content (PIC) were estimated for each SNP marker using PowerMarker v3.25 software [67]. Markers with a MajAF ≥ 0.95 and varieties with more than 15% of heterozygous genotyping data were removed from the data set to prevent bias during the association analysis.

Linkage disequilibrium

LD was evaluated on the whole genome and on each LG by calculating the correlation coefficient r2 between each pair of markers located on the same LG (linked markers) and for pairs of markers located on different LGs (unlinked markers). This analysis was performed using PLINK v1.07 software [68]. LD decay was evaluated using a nonlinear regression of the expected r2 as described by Sved [69] with E[r2] = 1/(1 + 4Nec) where c is the recombination rate in morgans and Ne the effective population size. The E[r2] was plotted against the genetic distance between molecular markers (in cM) in order to estimate the extend of LD with R software [70]. The r2 threshold value of 0.2 was chosen as the value below which the LD is no longer significant. Heat maps of LD between markers of the same LG were generated using the LDHeatMap package developed on R software [71].

Population structure and kinship

To investigate population structure, close markers that were in LD were first eliminated because these might cause a bias in the estimation of population structure. PLINK v1.07 [68] was used to select the markers. A sliding window of 20 SNP was defined with a 2-SNP step. In each window, the r2 coefficient was calculated for each possible marker pair. When the r2 value was greater than 0.2, one of the markers of the pair was eliminated so that only the markers that are not in LD were retained in the window. The set of markers obtained after this first analysis was then used to lead a principal component analysis (PCA) on the basis of genotyping data using the software EIGENSTRAT [72]. A Tracy-Widom test (alpha = 5%) was applied to identify significant axes of the PCA. The matrix of coordinates of the varieties on these significant axes, called the P matrix, was used to control the population structure in the association analysis.

An identity by state (IBS) kinship matrix was computed as a genotype similarity between the different pairs of markers using the algorithm “EMMA” in the GAPIT package [73]. The matrix of kinship values between pairs of varieties was used as a K matrix in the association analysis.

Genome-Wide Association analysis

Association mapping analysis was carried out using the GAPIT package [73]. All the markers with a MajAF < 0.95 were used in three different linear models to test the association between molecular markers and the G2 disease index trait. First, the GLM was applied. In this model, each marker is considered as an independent variable and there is no control over the population structure or kinship. Then, the K CMLM was used. In this model, a kinship matrix (K) is used to take into account the relationship between individual genotypes. Finally, the KP CMLM allows both the population kinship and structure to be controlled using the K and P matrices. Compared with a simple mixed linear model, compressed mixed linear models increase the statistical power and reduce computing time for large samples [74]. The K and KP CM linear models should allow better control of false associations caused by population structure.

Structural organisation of stem canker QTL

To study the organisation of the genetic factors controlling resistance to L. maculans in relation to the duplications of B. napus genome, our integrated genetic map was anchored on the Arabidopsis genome using the homology with the SNP context sequences [40]. A map with the duplicated blocks on our integrated map was constructed and blocks were presented in relation to the 24 blocks defined by Schranz et al. [17] and recently refined by Cheng et al. [18].

Availability of supporting data

The context sequences of the SNP markers mentioned as private in Additional file 6 can be requested from S. Faure (Sebastien.Faure@biogemma.com).

References

  1. Jackson S, Chen ZJ: Genomic and expression plasticity of polyploidy. Curr Opin Plant Biol. 2010, 13 (2): 153-159.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  2. Wang Y, Wang X, Paterson AH: Genome and gene duplications and gene expression divergence: a view from plants. Ann NY Acad Sci. 2012, 1256: 1-14.

    Article  PubMed  Google Scholar 

  3. Doyle JJ, Flagel LE, Paterson AH, Rapp RA, Soltis DE, Soltis PS, Wendel JF: Evolutionary Genetics of Genome Merger and Doubling in Plants. Annu Rev Genet. 2008, 42: 443-461.

    Article  CAS  PubMed  Google Scholar 

  4. Paterson AH: Polyploidy, evolutionary opportunity, and crop adaptation. Genetica. 2005, 123 (1–2): 191-196.

    Article  CAS  PubMed  Google Scholar 

  5. Borner A, Schumann E, Furste A, Coster H, Leithold B, Roder MS, Weber WE: Mapping of quantitative trait loci determining agronomic important characters in hexaploid wheat (Triticum aestivum L.). Theor Appl Genet. 2002, 105 (6–7): 921-936.

    PubMed  Google Scholar 

  6. Anderson OD, Abraham-Pierce FA, Tam A: Conservation in wheat high-molecular-weight glutenin gene promoter sequences: comparisons among loci and among alleles of the GLU-B1-1 locus. Theor Appl Genet. 1998, 96: 568-576.

    Article  CAS  Google Scholar 

  7. De Majnik J, Ogbonnaya FC, Moullet O, Lagudah ES: The Cre1 and Cre3 nematode resistance genes are located at homeologous loci in the wheat genome. Mol Plant-Microbe Interact. 2003, 16: 1129-1134.

    Article  CAS  PubMed  Google Scholar 

  8. Liu W, Kim MY, Kang YJ, Van K, Lee YH, Srinives P, Yuan DL, Lee SH: QTL identification of flowering time at three different latitudes reveals homeologous genomic regions that control flowering in soybean. Theor Appl Genet. 2011, 123 (4): 545-553.

    Article  PubMed  Google Scholar 

  9. Rong J, Feltus EA, Waghmare VN, Pierce GJ, Chee PW, Draye X, Saranga Y, Wright RJ, Wilkins TA, May OL, Smith CW, Gannaway JR, Wendel JR, Paterson AH: Meta-analysis of polyploid cotton QTL shows unequal contributions of subgenomes to a complex network of genes and gene clusters implicated in lint fiber development. Genetics. 2007, 176 (4): 2577-2588.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  10. Lerceteau-Kohler E, Moing A, Guerin G, Renaud C, Petit A, Rothan C, Denoyes B: Genetic dissection of fruit quality traits in the octoploid cultivated strawberry highlights the role of homoeo-QTL in their control. Theor Appl Genet. 2012, 124 (6): 1059-1077.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  11. U N: Genome analysis in Brassica with special reference to the experimental formation of B. napus and peculiar mode of fertilization. Jpn J Bot. 1935, 7: 389-452.

    Google Scholar 

  12. De Bodt S, Maere S, Van de Peer Y: Genome duplication and the origin of angiosperms. Trends Ecol Evol. 2005, 20 (11): 591-597.

    Article  PubMed  Google Scholar 

  13. Simillion C, Vandepoele K, Van Montagu MCE, Zabeau M, Van de Peer Y: The hidden duplication past of Arabidopsis thaliana. Proc Natl Acad Sci USA. 2002, 99 (21): 13627-13632.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  14. Rana D, Boogaart T, O’Neill CM, Hynes L, Bent E, Macpherson L, Park JY, Lim YP, Bancroft I: Conservation of the microstructure of genome segments in Brassica napus and its diploid relatives. Plant J. 2004, 40 (5): 725-733.

    Article  CAS  PubMed  Google Scholar 

  15. Parkin IAP, Sharpe AG, Lydiate DJ: Patterns of genome duplication within the Brassica napus genome. Genome. 2003, 46 (2): 291-303.

    Article  CAS  PubMed  Google Scholar 

  16. Parkin IAP, Gulden SM, Sharpe AG, Lukens L, Trick M, Osborn TC, Lydiate DJ: Segmental structure of the Brassica napus genome based on comparative analysis with Arabidopsis thaliana. Genetics. 2005, 171 (2): 765-781.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  17. Schranz ME, Lysak MA, Mitchell-Olds T: The ABC’s of comparative genomics in the Brassicaceae: building blocks of crucifer genomes. Trends Plant Sci. 2006, 11 (11): 535-542.

    Article  CAS  PubMed  Google Scholar 

  18. Cheng F, Mandakova T, Wu J, Xie Q, Lysak MA, Wang X: Deciphering the Diploid Ancestral Genome of the Mesohexaploid Brassica rapa. Plant Cell. 2013, 25: 1541-1554.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  19. Parkin I: Chasing Ghosts : Comparative Mapping in the Brassicaceae. Plant Genetics and Genomics : Crops and Models. Edited by: Schmidt R, Bancroft I. 2011, New York: Springer Science + Business Media, 153-170.

    Google Scholar 

  20. Chen G, Geng JF, Rahman M, Liu XP, Tu JX, Fu TD, Li GY, McVetty PBE, Tahir M: Identification of QTL for oil content, seed yield, and flowering time in oilseed rape (Brassica napus). Euphytica. 2010, 175 (2): 161-174.

    Article  CAS  Google Scholar 

  21. Long Y, Shi J, Qiu D, Li R, Zhang C, Wang J, Hou J, Zhao J, Shi L, Park BS, Choi SR, Lim YP, Meng J: Flowering time quantitative trait loci analysis of oilseed Brassica in multiple environments and genomewide alignment with Arabidopsis. Genetics. 2007, 177 (4): 2433-2444.

    CAS  PubMed Central  PubMed  Google Scholar 

  22. Delourme R, Brun H, Ermel M, Lucas MO, Vallee P, Domin C, Walton G, Li H, Sivasithamparam K, Barbetti MJ: Expression of resistance to Leptosphaeria maculans in Brassica napus double haploid lines in France and Australia is influenced by location. Ann Appl Biol. 2008, 153 (2): 259-269.

    Google Scholar 

  23. Yin XR, Yi B, Chen W, Zhang WJ, Tu JX, Fernando WGD, Fu TD: Mapping of QTLs detected in a Brassica napus DH population for resistance to Sclerotinia sclerotiorum in multiple environments. Euphytica. 2010, 173 (1): 25-35.

    Article  CAS  Google Scholar 

  24. Howell PM, Sharpe AG, Lydiate DJ: Homoeologous loci control the accumulation of seed glucosinolates in oilseed rape (Brassica napus). Genome. 2003, 46 (3): 454-460.

    Article  CAS  PubMed  Google Scholar 

  25. Wang J, Long Y, Wu BD, Liu J, Jiang CC, Shi L, Zhao JW, King GJ, Meng JL: The evolution of Brassica napus FLOWERING LOCUS T paralogues in the context of inverted chromosomal duplication blocks. BMC Evol Biol. 2009, 9: 271-

    Article  PubMed Central  PubMed  Google Scholar 

  26. Chen W, Zhang Y, Liu XP, Chen BY, Tu JX, Fu TD: Detection of QTL for six yield-related traits in oilseed rape (Brassica napus) using DH and immortalized F-2 populations. Theor Appl Genet. 2007, 115 (6): 849-858.

    Article  CAS  PubMed  Google Scholar 

  27. Zhao JW, Udall JA, Quijada PA, Grau CR, Meng JL, Osborn TC: Quantitative trait loci for resistance to Sclerotinia sclerotiorum and its association with a homeologous non-reciprocal transposition in Brassica napus L. Theor Appl Genet. 2006, 112 (3): 509-516.

    Article  CAS  PubMed  Google Scholar 

  28. Fitt BDL, Brun H, Barbetti MJ, Rimmer SR: World-wide importance of phoma stem canker (Leptosphaeria maculans and L. biglobosa) on oilseed rape (Brassica napus). Eur J Plant Pathol. 2006, 114 (1): 3-15.

    Article  Google Scholar 

  29. West JS, Kharbanda PD, Barbetti MJ, Fitt BDL: Epidemiology and management of Leptosphaeria maculans (phoma stem canker) on oilseed rape in Australia, Canada and Europe. Plant Pathol. 2001, 50 (1): 10-27.

    Article  Google Scholar 

  30. Jestin C, Vallée P, Domin C, Manzanares-Dauleux M, Delourme R: Assessment of a new strategy for selective phenotyping applied to complex traits in Brassica napus. Open J Genet Cytol. 2012, 2: 190-201.

    Article  CAS  Google Scholar 

  31. Pilet ML, Delourme R, Foisset N, Renard M: Identification of loci contributing to quantitative field resistance to blackleg disease, causal agent Leptosphaeria maculans (Desm.) Ces. et de Not., in Winter rapeseed (Brassica napus L.). Theor Appl Genet. 1998, 96 (1): 23-30.

    Article  Google Scholar 

  32. Pilet ML, Duplan G, Archipiano H, Barret P, Baron C, Horvais R, Tanguy X, Lucas MO, Renard M, Delourme R: Stability of QTL for field resistance to blackleg across two genetic backgrounds in oilseed rape. Crop Sci. 2001, 41 (1): 197-205.

    Article  CAS  Google Scholar 

  33. Kaur S, Cogan NOI, Ye G, Baillie RC, Hand ML, Ling AE, McGearey AK, Kaur J, Hopkins CJ, Todorovic M, Mountford H, Edwards D, Batley J, Burton W, Salisbury P, Gororo N, Marcroft S, Kearney G, Smith KF, Forster JW, Spangenberg GC: Genetic map construction and QTL mapping of resistance to blackleg (Leptosphaeria maculans) disease in Australian canola (Brassica napus L.) cultivars. Theor Appl Genet. 2009, 120: 71-83.

    Article  CAS  PubMed  Google Scholar 

  34. Delourme R, Chevre AM, Brun H, Rouxel T, Balesdent MH, Dias JS, Salisbury P, Renard M, Rimmer SR: Major gene and polygenic resistance to Leptosphaeria maculans in oilseed rape (Brassica napus). Eur J Plant Pathol. 2006, 114 (1): 41-52.

    Article  Google Scholar 

  35. Comai L: The advantages and disadvantages of being polyploid. Nat Rev Genet. 2005, 6 (11): 836-846.

    Article  CAS  PubMed  Google Scholar 

  36. Osborn TC, Pires JC, Birchler JA, Auger DL, Chen ZJ, Lee HS, Comai L, Madlung A, Doerge RW, Colot V, Martienssen RA: Understanding mechanisms of novel gene expression in polyploids. Trends Genet. 2003, 19 (3): 141-147.

    Article  CAS  PubMed  Google Scholar 

  37. Balding DJ: A tutorial on statistical methods for population association studies. Nat Rev Genet. 2006, 7 (10): 781-791.

    Article  CAS  PubMed  Google Scholar 

  38. McCarthy MI, Hirschhorn JN: Genome-wide association studies: potential next steps on a genetic journey. Hum Mol Genet. 2008, 17: R156-R165.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  39. Kaur S, Francki MG, Forster JW: Identification, characterization and interpretation of single-nucleotide sequence variation in allopolyploid crop species. Plant Biotechnol J. 2012, 10 (2): 125-138.

    Article  CAS  PubMed  Google Scholar 

  40. Delourme R, Falentin C, Fopa Fomeju B, Boillot M, Lassalle G, André I, Duarte J, Gauthier V, Lucante N, Marty A, Pauchon M, Pichon JP, Ribière N, Trotoux G, Blanchard P, Rivière N, Martinant JP, Pauquet J: High-density SNP-based genetic map development and linkage disequilibroum assessment in Brassica napus L. BMC Genomics. 2013, 14: 120-

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  41. Bradbury P, Parker T, Hamblin MT, Jannink J-L: Assessment of Power and False Discovery Rate in Genome-Wide Association Studies using the BarleyCAP Germplasm. Crop Sci. 2011, 51 (1): 52-59.

    Article  Google Scholar 

  42. Bus A, Koerber N, Snowdon RJ, Stich B: Patterns of molecular variation in a species-wide germplasm set of Brassica napus. Theor Appl Genet. 2011, 123 (8): 1413-1423.

    Article  PubMed  Google Scholar 

  43. Ecke W, Clemens R, Honsdorf N, Becker HC: Extent and structure of linkage disequilibrium in canola quality winter rapeseed (Brassica napus L.). Theor Appl Genet. 2010, 120 (5): 921-931.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  44. Yu JM, Pressoir G, Briggs WH, Bi IV, Yamasaki M, Doebley JF, McMullen MD, Gaut BS, Nielsen DM, Holland JB, Kresovich S, Buckler ES: A unified mixed-model method for association mapping that accounts for multiple levels of relatedness. Nat Genet. 2006, 38 (2): 203-208.

    Article  CAS  PubMed  Google Scholar 

  45. Atwell S, Huang YS, Vilhjalmsson BJ, Willems G, Horton M, Li Y, Meng D, Platt A, Tarone AM, Hu TT, Jiang R, Muliyati NW, Zhang X, Amer MA, Baxter I, Brachi B, Chory J, Dean C, Debieu M, de Meaux J, Ecker JR, Faure N, Kniskern JM, Jones JDG, Michael T, Nemri A, Roux F, Salt DE, Tang C, Todesco M, et al: Genome-wide association study of 107 phenotypes in Arabidopsis thaliana inbred lines. Nature. 2010, 465: 7298-

    Article  Google Scholar 

  46. Jestin C, Lode M, Vallee P, Domin C, Falentin C, Horvais R, Coedel S, Manzanares-Dauleux MJ, Delourme R: Association mapping of quantitative resistance for Leptosphaeria maculans in oilseed rape (Brassica napus L.). Mol Breed. 2011, 27 (3): 271-287.

    Article  Google Scholar 

  47. Jestin C: Approche conjointe de cartographie par analyse de liaison et génétique d’association pour appréhender l’architecture génétique de la résistance quantitative du colza à la nécrose du collet. PhD Thesis Agrocampus ouest. 2010, 224-

    Google Scholar 

  48. Zou J, Jiang C, Cao Z, Li R, Long Y, Chen S, Meng J: Association mapping of seed oil content in Brassica napus and comparison with quantitative trait loci identified from linkage mapping. Genome. 2010, 53: 908-916.

    Article  CAS  PubMed  Google Scholar 

  49. Cheung F, Trick M, Drou N, Lim YP, Park JY, Kwon SJ, Kim JA, Scott R, Pires JC, Paterson AH, Town C, Bancroft I: Comparative Analysis between Homoeologous Genome Segments of Brassica napus and Its Progenitor Species Reveals Extensive Sequence-Level Divergence. Plant Cell. 2009, 21 (7): 1912-1928.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  50. Wendel JF, Cronn RC: Polyploidy and the evolutionary history of cotton. Adv Agron. 2003, 78: 139-186.

    Article  Google Scholar 

  51. Li XN, Ramchiary N, Dhandapani V, Choi SR, Hur Y, Nou IS, Yoon MK, Lim YP: Quantitative Trait Loci Mapping in Brassica rapa Revealed the Structural and Functional Conservation of Genetic Loci Governing Morphological and Yield Component Traits in the A, B, and C Subgenomes of Brassica Species. DNA Res. 2013, 20 (1): 1-16.

    Article  PubMed Central  PubMed  Google Scholar 

  52. Adams KL: Evolution of duplicate gene expression in polyploid and hybrid plants. J Hered. 2007, 98 (2): 136-141.

    Article  CAS  PubMed  Google Scholar 

  53. Thomas BC, Pedersen B, Freeling M: Following tetraploidy in an Arabidopsis ancestor, genes were removed preferentially from one homeolog leaving clusters enriched in dose-sensitive genes. Genome Res. 2006, 16 (7): 934-946.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  54. Woodhouse MR, Schnable JC, Pedersen BS, Lyons E, Lisch D, Subramaniam S, Freeling M: Following Tetraploidy in Maize, a Short Deletion Mechanism Removed Genes Preferentially from One of the Two Homeologs. PLoS Biol. 2010, 8: 6-

    Article  Google Scholar 

  55. Schnable JC, Springer NM, Freeling M: Differentiation of the maize subgenomes by genome dominance and both ancient and ongoing gene loss. Proc Natl Acad Sci USA. 2011, 108 (10): 4069-4074.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  56. Tang H, Woodhouse MR, Cheng F, Schnable JC, Pedersen BS, Conant G, Wang X, Freeling M, Pires JC: Altered Patterns of Fractionation and Exon Deletions in Brassica rapa Support a Two-Step Model of Paleohexaploidy. Genetics. 2012, 190: 1563-1574.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  57. Blanc G, Wolfe KH: Functional divergence of duplicated genes formed by polyploidy during Arabidopsis evolution. Plant Cell. 2004, 16 (7): 1679-1691.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  58. Seoighe C, Gehring C: Genome duplication led to highly selective expansion of the Arabidopsis thaliana proteome. Trends Genet. 2004, 20 (10): 461-464.

    Article  CAS  PubMed  Google Scholar 

  59. Wu YF, Zhu ZG, Ma LG, Chen MS: The preferential retention of starch synthesis genes reveals the impact of whole-genome duplication on grass evolution. Mol Biol Evol. 2008, 25 (6): 1003-1006.

    Article  CAS  PubMed  Google Scholar 

  60. Severin AJ, Cannon SB, Graham MM, Grant D, Shoemaker RC: Changes in Twelve Homoeologous Genomic Regions in Soybean following Three Rounds of Polyploidy. Plant Cell. 2011, 23 (9): 3129-3136.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  61. Wang XW, Wang HZ, Wang J, Sun RF, Wu J, Liu SY, Bai YQ, Mun JH, Bancroft I, Cheng F, Huang SW, Li XX, Hua W, Wang JY, Wang XY, Freeling M, Pires JC, Paterson AH, Chalhoub B, Wang B, Hayward A, Sharpe AG, Park BS, Weisshaar B, Liu BH, Li B, Liu B, Tong CB, Song C, Duran C, et al: The genome of the mesopolyploid crop species Brassica rapa. Nat Genet. 2011, 43 (10): 1035-U1157.

    Article  CAS  PubMed  Google Scholar 

  62. Li J, Wu XM: Genome-wide identification, classification and expression analysis of genes encoding putative fasciclin-like arabinogalactan proteins in Chinese cabbage (Brassica rapa L.). Mol Biol Rep. 2012, 39 (12): 10541-10555.

    Article  CAS  Google Scholar 

  63. Mun JH, Yu HJ, Shin JY, Oh M, Hwang HJ, Chung H: Auxin response factor gene family in Brassica rapa: genomic organization, divergence, expression, and evolution. Mol Gen Genomics. 2012, 287 (10): 765-784.

    Article  CAS  Google Scholar 

  64. Zhao JW, Buchwaldt L, Rimmer SR, Brkic M, Bekkaoui D, Hegedus D: Differential expression of duplicated peroxidase genes in the allotetraploid Brassica napus. Plant Physiol Biochem. 2009, 47 (7): 653-656.

    Article  CAS  PubMed  Google Scholar 

  65. Albertin W, Balliau T, Brabant P, Chevre AM, Eber F, Malosse C, Thiellement H: Numerous and rapid nonstochastic modifications of gene products in newly synthesized Brassica napus allotetraploids. Genetics. 2006, 173 (2): 1101-1113.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  66. Gaeta RT, Pires JC, Iniguez-Luy F, Leon E, Osborn TC: Genomic changes in resynthesized Brassica napus and their effect on gene expression and phenotype. Plant Cell. 2007, 19 (11): 3403-3417.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  67. Liu KJ, Muse SV: PowerMarker: an integrated analysis environment for genetic marker analysis. Bioinformatics. 2005, 21 (9): 2128-2129.

    Article  CAS  PubMed  Google Scholar 

  68. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, Maller J, Sklar P, De Bakker PIW, Daly MJ, Sham PC: PLINK: A tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007, 81 (3): 559-575.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  69. Sved J: Linkage disequilibrium and homozygosity of chromosome segments in finite populations. Theor Popul Biol. 1971, 2: 125-141.

    Article  CAS  PubMed  Google Scholar 

  70. R Development Core team: R: A language and environment for statistical computing. 2.14 edn: R Foundation for Statistical Computing. 2012

    Google Scholar 

  71. Shin J-H, Blay S, McNeney B, Graham J: LDheatmap: An R Function for Graphical Display of Pairwise Linkage Disequilibria Between Single Nucleotide Polymorphism. J Stat Softw. 2006, 16: 1-10.

    Article  Google Scholar 

  72. Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D: Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 2006, 38 (8): 904-909.

    Article  CAS  PubMed  Google Scholar 

  73. Lipka AE, Tian F, Wang QS, Peiffer J, Li M, Bradbury PJ, Gore MA, Buckler ES, Zhang ZW: GAPIT: genome association and prediction integrated tool. Bioinformatics. 2012, 28 (18): 2397-2399.

    Article  CAS  PubMed  Google Scholar 

  74. Zhang Z, Ersoz E, Lai C-Q, Todhunter RJ, Tiwari HK, Gore MA, Bradbury PJ, Yu J, Arnett DK, Ordovas JM, Buckler ES: Mixed linear model approach adapted for genome-wide association studies. Nat Genet. 2010, 42: 355-U118.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

Download references

Acknowledgements

This work was supported by the French Institut National de la Recherche Agronomique (INRA) and the Région Bretagne. We thank the team of the INRA Experimental Unit at Le Rheu for performing the disease evaluation trials.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Régine Delourme.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

BF performed linkage disequilibrium, structure and GWA analyses. CF carried out genetic mapping and BLAST analyses for the construction of the integrated map used in this study and the positioning of duplicated blocks on this map. GL contributed to genetic mapping and BLAST analyses. RD and MMD helped analyse the results and coordinated the project. BF, RD and MMD wrote the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

12864_2013_6190_MOESM1_ESM.xlsx

Additional file 1: Table S1: Extent of linkage disequilibrium on each linkage group of B. napus. Extent of linkage disequilibrium was calculated according to the formula E[r2] = 1/(1 + 4Nec) as described in Sved et al. [69]. The effective population size (Ne) and the distance to which the LD extended (c) at a r2 threshold of 0.2 are given for each linkage group, for the A and C genomes and for the whole genome. (XLSX 10 KB)

12864_2013_6190_MOESM2_ESM.pdf

Additional file 2: Figure S1: Linkage disequilibrium related to the genetic distance between markers along B. napus genome. The linkage disequilibrium decay was measured in the whole panel of 116 winter oilseed rape cultivars. The red curve shows the non-linear regression trend line of r2 versus genetic distance. (PDF 125 KB)

12864_2013_6190_MOESM3_ESM.pdf

Additional file 3: Figure S2: The linkage disequilibrium pattern and results for association analysis of resistance to stem canker for each linkage group. Negative log10 p values obtained from the K and KP CML models were plotted against the genetic distance (in cM) for each linkage group. The corresponding linkage disequilibrium pattern calculated between pairs of tested markers for association is presented below the Manhattan plot. The more the colour is closer to red, the higher the linkage disequilibrium is. (PDF 1 MB)

12864_2013_6190_MOESM4_ESM.pdf

Additional file 4: Figure S3: Principal Component Analysis of 116 winter oilseed rape cultivars based on 727 SNP markers. The proportion of variance explained by the two principal components (Eigenvector 1 and 2) is indicated in parentheses. Varieties were tagged according to their euric acid and glucosinolate contents, “00” indicates a double low in erucic acid and glucosinolate content, “0+” indicated low erucic acid and high glucosinolate content, “+0” indicates high erucic acid and low glucosinolate content and “++” indicated high erucic acid and high glucosinolate content. The spring type variety ‘Yudal’ is circled in red on the right of the graph. (PDF 167 KB)

12864_2013_6190_MOESM5_ESM.pdf

Additional file 5: Figure S4: Distribution of the kinship coefficients calculated between 116 winter oilseed rape cultivars. (PDF 171 KB)

12864_2013_6190_MOESM6_ESM.xlsx

Additional file 6: Table S2: List of the 321 resistance associated markers. The positions (in cM) of the markers on our linkage map, including information about their assignation to AK block are provided. When possible, context sequences of associated SNP are given. For a set of markers mentioned “private”, context sequences are available under request (see ‘Availability of supporting data’ section). (XLSX 33 KB)

12864_2013_6190_MOESM7_ESM.xlsx

Additional file 7: Table S3: Position of resistance-associated markers in relation to the duplicated blocks in B. napus genome. P values given in the table for the K and KP linear models were those not corrected with the FDR test. Information about the co-localisation of resistance-associated markers with previously identify QTL in our laboratory (in a double haploid -DH- and/or a connected -CP- populations) is indicated. (XLSX 20 KB)

12864_2013_6190_MOESM8_ESM.xlsx

Additional file 8: Table S4: List of Brassica napus varieties used for the genome-wide association analysis. *The G2 index is a mean of three replicates. **“00” indicates a double low in erucic acid and glucosinolate content, “0+” indicated low erucic acid and high glucosinolate content, “+0” indicates high erucic acid and low glucosinolate content and “++” indicated high erucic acid and high glucosinolate content. The dash indicate missing data. IHAR: Instytut Hodowli I Aklimatyzacji Roslin, Poznan; INRA: Institut National de la Recherche Agronomique; JD: Jouffray-Drillaud; KWS: KWS Saat AG; Momont: Sarl Adrien Momont & Fils; NK: Syngenta Seeds SAS; NPZ: NPZ Lembke Semences Sarl; RAPS Gbr: Raps GbR Saatzucht Lundsgaard; SW: Svalöf Weibull AB. (XLSX 15 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Fopa Fomeju, B., Falentin, C., Lassalle, G. et al. Homoeologous duplicated regions are involved in quantitative resistance of Brassica napus to stem canker. BMC Genomics 15, 498 (2014). https://doi.org/10.1186/1471-2164-15-498

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-15-498

Keywords