Many previous studies have provided evidence for genome changes in polyploids, but there are little data on the overall population dynamics of genome change and whether it causes phenotypic variability. We analyzed genetic, epigenetic, gene expression, and phenotypic changes in ;50 resynthesized Brassica napus lines independently derived by hybridizing double haploids of Brassica oleracea and Brassica rapa. A previous analysis of the first generation (S0) found that genetic changes were rare, and cytosine methylation changes were frequent. Our analysis of a later generation found that most S0 methylation changes remained fixed in their S5 progeny, although there were some reversions and new methylation changes. Genetic changes were much more frequent in the S5 generation, occurring in every line with lines normally distributed for number of changes. Genetic changes were detected on 36 of the 38 chromosomes of the S5 allopolyploids and were not random across the genome. DNA fragment losses within lines often occurred at linked marker loci, and most fragment losses co-occurred with intensification of signal from homoeologous markers, indicating that the changes were due to homoeologous nonreciprocal transpositions (HNRTs). HNRTs between chromosomes A1 and C1 initiated in early generations, occurred in successive generations, and segregated, consistent with a recombination mechanism. HNRTs and deletions were correlated with qualitative changes in the expression of specific homoeologous genes and anonymous cDNA amplified fragment length polymorphisms and with phenotypic variation among S5 polyploids. Our data indicate that exchanges among homoeologous chromosomes are a major mechanism creating novel allele combinations and phenotypic variation in newly formed B. napus polyploids.

Key messageTheBrassica napusIllumina array provides genome-wide markers linked to the available genome sequence, a significant tool for genetic analyses of the allotetraploidB. napusand its progenitor diploid genomes.AbstractA high-density single nucleotide polymorphism (SNP) Illumina Infinium array, containing 52,157 markers, was developed for the allotetraploid Brassica napus. A stringent selection process employing the short probe sequence for each SNP assay was used to limit the majority of the selected markers to those represented a minimum number of times across the highly replicated genome. As a result approximately 60 % of the SNP assays display genome-specificity, resolving as three clearly separated clusters (AA, AB, and BB) when tested with a diverse range of B. napus material. This genome specificity was supported by the analysis of the diploid ancestors of B. napus, whereby 26,504 and 29,720 markers were scorable in B. oleracea and B. rapa, respectively. Forty-four percent of the assayed loci on the array were genetically mapped in a single doubled-haploid B. napus population allowing alignment of their physical and genetic coordinates. Although strong conservation of the two positions was shown, at least 3 % of the loci were genetically mapped to a homoeologous position compared to their presumed physical position in the respective genome, underlying the importance of genetic corroboration of locus identity. In addition, the alignments identified multiple rearrangements between the diploid and tetraploid Brassica genomes. Although mostly attributed to genome assembly errors, some are likely evidence of rearrangements that occurred since the hybridisation of the progenitor genomes in the B. napus nucleus. Based on estimates for linkage disequilibrium decay, the array is a valuable tool for genetic fine mapping and genome-wide association studies in B. napus and its progenitor genomes.Electronic supplementary materialThe online version of this article (doi:10.1007/s00122-016-2746-7) contains supplementary material, which is available to authorized users.

Publicly available genomic tools help researchers integrate information and make new discoveries. In this paper, we describe the development of immortal mapping populations of rapid cycling, self-compatible lines, molecular markers, and linkage maps for Brassica rapa and B. oleracea and make the data and germplasm available to the Brassica research community. The B. rapa population consists of 160 recombinant inbred (RI) lines derived from the cross of highly inbred lines of rapid cycling and yellow sarson B. rapa. The B. oleracea population consists of 155 double haploid (DH) lines derived from an F1 cross between two DH lines, rapid cycling and broccoli. A total of 120 RFLP probes, 146 SSR markers, and one phenotypic trait (flower color) were used to construct genetic linkage maps for both species. The B. rapa map consists of 224 molecular markers distributed along 10 linkage groups (A1-A10) with a total distance of 1125.3 cM and a marker density of 5.7 cM/marker. The B. oleracea genetic map consists of 279 molecular markers and one phenotypic marker distributed along nine linkage groups (C1-C9) with a total distance of 891.4 cM and a marker density of 3.2 cM/marker. A syntenic analysis with Arabidopsis thaliana identified collinear genomic blocks that are in agreement with previous studies, reinforcing the idea of conserved chromosomal regions across the Brassicaceae.

BackgroundFlowering time, plant height and seed yield are strongly influenced by climatic and day-length adaptation in crop plants. To investigate these traits under highly diverse field conditions in the important oilseed crop Brassica napus, we performed a genome-wide association study using data from diverse agroecological environments spanning three continents.MethodsA total of 158 European winter-type B.napus inbred lines were genotyped with 21,623 unique, single-locus single-nucleotide polymorphism (SNP) markers using the Brassica 60 K-SNP Illumina® Infinium consortium array. Phenotypic associations were calculated in the panel over the years 2010–2012 for flowering time, plant height and seed yield in 5 highly diverse locations in Germany, China and Chile, adding up to 11 diverse environments in total.ResultsWe identified 101 genome regions associating with the onset of flowering, 69 with plant height, 36 with seed yield and 68 cross-trait regions with potential adaptive value. Within these regions, B.napus orthologs for a number of candidate adaptation genes were detected, including central circadian clock components like CIRCADIAN CLOCK- ASSOCIATED 1 (Bna.CCA1) and the important flowering-time regulators FLOWERING LOCUS T (Bna.FT) and FRUITFUL (Bna.FUL).DiscussionGene ontology (GO) enrichment analysis of candidate regions suggested that selection of genes involved in post-transcriptional and epigenetic regulation of flowering time may play a potential role in adaptation of B. napus to highly divergent environments. The classical flowering time regulators Bna.FLC and Bna.CO were not found among the candidate regions, although both show functional variation. Allelic effects were additive for plant height and yield, but not for flowering time. The scarcity of positive minor alleles for yield in this breeding pool points to a lack of diversity for adaptation that could restrict yield gain in the face of environmental change.ConclusionsOur study provides a valuable framework to further improve the adaptability and yield stability of this recent allopolyploid crop under changing environments. The results suggest that flowering time regulation within an adapted B. napus breeding pool is driven by a high number of small modulating processes rather than major transcription factors like Bna.CO. In contrast, yield regulation appears highly parallel, therefore yield could be increased by pyramiding positively associated haplotypes.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-1950-1) contains supplementary material, which is available to authorized users.

The traditional development of simple sequence repeat (SSR) or microsatellite markers by probe hybridization can be time-consuming and requires the use of specialized laboratory equipment. In this study, probe hybridization was circumvented by using sequence information on 3,500 genomic clones mainly from Brassica oleracea to identify di, tri, tetra and penta-nucleotide repeats. A total of 587 primer pairs flanking SSR were developed using this approach. From these, 420 SSR markers amplified DNA in two parental lines of B. rapa (26% were polymorphic) and 523 in two parental lines of B. oleracea (32% were polymorphic). A diverse array of motif types was identified, characterized and compared with traditional SSR detection methods. The most abundant motifs found were di- (38%) and trinucleotides (33%) followed by penta- (16%) and tetranucleotide (13%) motifs. The type of motif class, motif length and repeat were not indicative of polymorphisms. The frequency of B. oleracea SSRs in genomic shotgun sequence was estimated to be 1 every 4 Kb. In general, the average motif length and repeat numbers were shorter than those obtained previously by probe hybridization, and they contained a more balanced representation of SSR motif types in the genome by identifying those that do not hybridize well to DNA probes. Brassica genomic DNA sequence information is a promising resource for developing a large number of SSR molecular markers in Brassica species.

Subfunctionalization is the process by which a pair of duplicated genes, or paralogs, experiences a reduction of individual expression patterns or function while still reproducing the complete expression pattern and function of the ancestral gene. Two germin-like protein (GLP)-encoding genes, GerB and GerF, are paralogs that belong to a small gene family in barley (Hordeum vulgare). Both genes share high nucleotide sequence similarity in coding and noncoding regions and encode identical apoplastic proteins. The use of RNA gel blots, coupled with single-stranded conformation polymorphism (SSCP) analysis of RT–PCR products, elucidated the developmental and tissue-specific expression patterns of each gene. Individual expression patterns provided evidence of both overlapping redundancy and early subfunctionalization. GerB is predominantly expressed in developing shoots, while GerF is predominantly expressed in seedling roots, developing spikes, and pericarp/testa. GerF promoter deletion studies located a region (−356/−97) responsible for high promoter activity and showed the ability of GerB and GerF upstream regions to drive gfp expression in coleoptiles, epicarps, and lemma/palea of developing spikes. The observed expression patterns are consistent with proposed roles in plant development and defense mechanisms for this gene family. These roles may explain why redundancy has been selectively maintained in this duplicate gene pair.

Targeted genomic selection methodologies, or sequence capture, allow for DNA enrichment and large-scale resequencing and characterization of natural genetic variation in species with complex genomes, such as rapeseed canola (Brassica napus L., AACC, 2n=38). The main goal of this project was to combine sequence capture with next generation sequencing (NGS) to discover single nucleotide polymorphisms (SNPs) in specific areas of the B. napus genome historically associated (via quantitative trait loci –QTL– analysis) to traits of agronomical and nutritional importance. A 2.1 million feature sequence capture platform was designed to interrogate DNA sequence variation across 47 specific genomic regions, representing 51.2 Mb of the Brassica A and C genomes, in ten diverse rapeseed genotypes. All ten genotypes were sequenced using the 454 Life Sciences chemistry and to assess the effect of increased sequence depth, two genotypes were also sequenced using Illumina HiSeq chemistry. As a result, 589,367 potentially useful SNPs were identified. Analysis of sequence coverage indicated a four-fold increased representation of target regions, with 57% of the filtered SNPs falling within these regions. Sixty percent of discovered SNPs corresponded to transitions while 40% were transversions. Interestingly, fifty eight percent of the SNPs were found in genic regions while 42% were found in intergenic regions. Further, a high percentage of genic SNPs was found in exons (65% and 64% for the A and C genomes, respectively). Two different genotyping assays were used to validate the discovered SNPs. Validation rates ranged from 61.5% to 84% of tested SNPs, underpinning the effectiveness of this SNP discovery approach. Most importantly, the discovered SNPs were associated with agronomically important regions of the B. napus genome generating a novel data resource for research and breeding this crop species.

Single nucleotide polymorphisms (SNPs) have rapidly become the molecular marker of choice in plant and animal association mapping (AM) studies. In this work, a genome-wide association study (GWAS) and candidate quantitative trait loci (cQTL) approaches were used to identify SNP markers associated with seed quality traits, in a Brassica napus L. association panel composed of 89 adapted winter oilseed rape accessions. Six seed quality traits (oil and protein content, linolenic acid, total glucosinolates, hemicellulose and cellulose content) were evaluated in two different locations for two seasons. For GWAS, 4025 SNP markers evenly distributed along the B. napus genome were genotyped using a 6K Illumina array platform. For cQTL, 100 SNP markers previously discovered in genomic regions underlying seed quality QTL were genotyped using a competitive allele-specific PCR (KASPar). Analysis of the population structure revealed the presence of two weakly differentiated subpopulations (F ST = 0.037), with 82 % of the pairwise kinship comparisons ranging from 0 to 0.1. The GWAS approach resulted in the identification of 17 and 5 significant associations for seed glucosinolate content and seed hemicellulose content, respectively. The cQTL approach identified 4 significant associations for seed glucosinolate content and 6 significant associations for seed hemicellulose content. The associated SNPs were consistently identified across environments and were mapped to previously reported QTL. These results illustrate the suitability of AM to identify SNP markers associated with seed quality traits in B. napus.