Epigenomic and structural events preclude recombination in Brassica napus.

Meiotic recombination is a major evolutionary process generating genetic diversity at each generation in sexual organisms. However, this process is highly regulated, with the majority of crossovers lying in the distal chromosomal regions that harbor low DNA methylation levels. Even in these regions, some islands without recombination remain, for which we investigated the underlying causes. Genetic maps were established in two Brassica napus hybrids to detect the presence of such large nonrecombinant islands. The role played by DNA methylation and structural variations in this local absence of recombination was determined by performing bisulfite sequencing and whole genome comparisons. Inferred structural variations were validated using either optical mapping or oligo fluorescence in situ hybridization. Hypermethylated or inverted regions between Brassica genomes were associated with the absence of recombination. Pairwise comparisons of nine B. napus genome assemblies revealed that such inversions occur frequently and may contain key agronomic genes such as resistance to biotic stresses. We conclude that such islands without recombination can have different origins, such as DNA methylation or structural variations in B. napus. It is thus essential to take into account these features in breeding programs as they may hamper the efficient combination of favorable alleles in elite varieties.

Untangling structural factors driving genome stabilization in nascent Brassica napus allopolyploids.

Allopolyploids have globally higher fitness than their diploid progenitors; however, by comparison, most resynthesized allopolyploids have poor fertility and highly unstable genome. Elucidating the evolutionary processes promoting genome stabilization and fertility is thus essential to comprehend allopolyploid success. Using the Brassica model, we mimicked the speciation process of a nascent allopolyploid species by resynthesizing allotetraploid Brassica napus and systematically selecting for euploid individuals over eight generations in four independent allopolyploidization events with contrasted genetic backgrounds, cytoplasmic donors, and polyploid formation type. We evaluated the evolution of meiotic behavior and fertility and identified rearrangements in S1 to S9 lineages to explore the positive consequences of euploid selection on B. napus genome stability. Recurrent selection of euploid plants for eight generations drastically reduced the percentage of aneuploid progenies as early as the fourth generation, concomitantly with a decrease in number of newly fixed homoeologous rearrangements. The consequences of homoeologous rearrangements on meiotic behavior and seed number depended strongly on the genetic background and cytoplasm donor. The combined use of both self-fertilization and recurrent euploid selection allowed identification of genomic regions associated with fertility and meiotic behavior, providing complementary evidence to explain B. napus speciation success.

Amplifying recombination genome-wide and reshaping crossover landscapes in Brassicas.

Meiotic recombination by crossovers (COs) is tightly regulated, limiting its key role in producing genetic diversity. However, while COs are usually restricted in number and not homogenously distributed along chromosomes, we show here how to disrupt these rules in Brassica species by using allotriploid hybrids (AAC, 2n = 3x = 29), resulting from the cross between the allotetraploid rapeseed (B. napus, AACC, 2n = 4x = 38) and one of its diploid progenitors (B. rapa, AA, 2n = 2x = 20). We produced mapping populations from different genotypes of both diploid AA and triploid AAC hybrids, used as female and/or as male. Each population revealed nearly 3,000 COs that we studied with SNP markers well distributed along the A genome (on average 1 SNP per 1.25 Mbp). Compared to the case of diploids, allotriploid hybrids showed 1.7 to 3.4 times more overall COs depending on the sex of meiosis and the genetic background. Most surprisingly, we found that such a rise was always associated with (i) dramatic changes in the shape of recombination landscapes and (ii) a strong decrease of CO interference. Hybrids carrying an additional C genome exhibited COs all along the A chromosomes, even in the vicinity of centromeres that are deprived of COs in diploids as well as in most studied species. Moreover, in male allotriploid hybrids we found that Class I COs are mostly responsible for the changes of CO rates, landscapes and interference. These results offer the opportunity for geneticists and plant breeders to dramatically enhance the generation of diversity in Brassica species by disrupting the linkage drag coming from limits on number and distribution of COs.

The poor lonesome A subgenome of Brassica napus var. Darmor (AACC) may not survive without its mate.

Constitutive genomes of allopolyploid species evolve throughout their life span. However, the consequences of long-term alterations on the interdependency between each original genome have not been established. Here, we attempted an approach corresponding to subgenome extraction from a previously sequenced natural allotetraploid, offering a unique opportunity to evaluate plant viability and structural evolution of one of its diploid components. We employed two different strategies to extract the diploid AA component of the Brassica napus variety 'Darmor' (AACC, 2n = 4x = 38) and we assessed the genomic structure of the latest AA plants obtained (after four to five rounds of selection), using a 60K single nucleotide polymorphism Illumina array. Only one strategy was successful and the diploid AA plants that were structurally characterized presented a lower proportion of the B. napus A subgenome extracted than expected. In addition, our analyses revealed that some genes lost in a polyploid context appeared to be compensated for plant survival, either by conservation of genomic regions from B. rapa, used in the initial cross, or by some introgressions from the B. napus C subgenome. We conclude that as little as c. 7500 yr of coevolution could lead to subgenome interdependency in the allotetraploid B. napus as a result of structural modifications.

Homoeologous Chromosome Sorting and Progression of Meiotic Recombination in Brassica napus: Ploidy Does Matter!

Meiotic recombination is the fundamental process that produces balanced gametes and generates diversity within species. For successful meiosis, crossovers must form between homologous chromosomes. This condition is more difficult to fulfill in allopolyploid species, which have more than two sets of related chromosomes (homoeologs). Here, we investigated the formation, progression, and completion of several key hallmarks of meiosis in Brassica napus (AACC), a young polyphyletic allotetraploid crop species with closely related homoeologous chromosomes. Altogether, our results demonstrate a precocious and efficient sorting of homologous versus homoeologous chromosomes during early prophase I in two representative B. napus accessions that otherwise show a genotypic difference in the progression of homologous recombination. More strikingly, our detailed comparison of meiosis in near isogenic allohaploid and euploid plants showed that the mechanism(s) promoting efficient chromosome sorting in euploids is adjusted to promote crossover formation between homoeologs in allohaploids. This suggests that, in contrast to other polyploid species, chromosome sorting is context dependent in B. napus.

Crossover rate between homologous chromosomes and interference are regulated by the addition of specific unpaired chromosomes in Brassica.

Recombination is a major mechanism generating genetic diversity, but the control of the crossover rate remains a key question. In Brassica napus (AACC, 2n = 38), we can increase the homologous recombination between A genomes in AAC hybrids. Hypotheses for this effect include the number of C univalent chromosomes, the ratio between univalents and bivalents and, finally, which of the chromosomes are univalents. To test these hypotheses, we produced AA hybrids with zero, one, three, six or nine additional C chromosomes and four different hybrids carrying 2n = 32 and 2n = 35 chromosomes. The genetic map lengths for each hybrid were established to compare their recombination rates. The rates were 1.4 and 2.7 times higher in the hybrids having C6 or C9 alone than in the control (0C). This enhancement reached 3.1 and 4.1 times in hybrids carrying six and nine C chromosomes, and it was also higher for each pair of hybrids carrying 2n = 32 or 2n = 35 chromosomes, with a dependence on which chromosomes remained as univalents. We have shown, for the first time, that the presence of one chromosome, C9 , affects significantly the recombination rate and reduces crossover interference. This result will have fundamental implications on the regulation of crossover frequency.

Non-random distribution of extensive chromosome rearrangements in Brassica napus depends on genome organization.

Chromosome rearrangements are common, but their dynamics over time, mechanisms of occurrence and the genomic features that shape their distribution and rate are still poorly understood. We used allohaploid Brassica napus (AC, n = 19) as a model to analyze the effect of genomic features on the formation and diversity of meiotically driven chromosome rearrangements. We showed that allohaploid B. napus meiosis leads to extensive new structural diversity. Almost every allohaploid offspring carried a unique combination of multiple rearrangements throughout the genome, and was thus structurally differentiated from both its haploid parent and its sister plants. This large amount of genome reshuffling was remarkably well-tolerated in the heterozygous state, as neither male nor female fertility were strongly reduced, and meiosis behavior was normal in most cases. We also used a quantitative statistical model, which accounted for 75% of the observed variation in rearrangement rates, to show that the distribution of meiotically driven chromosome rearrangements was not random but was shaped by three principal genomic features. In descending order of importance, the rate of marker loss increased strongly with genetic distance from the centromere, the degree of collinearity between chromosomes, and the genome of origin (A < C). Overall, our results demonstrate that B. napus accumulates a large number of genetic changes, but these rearrangements are not randomly distributed in the genome. The structural genetic diversity produced by the allohaploid pathway and its role in the evolution of polyploid species compared to diploid meiosis are discussed.

Allopolyploidy has a moderate impact on restructuring at three contrasting transposable element insertion sites in resynthesized Brassica napus allotetraploids.

The role played by whole-genome duplication (WGD) in evolution and adaptation is particularly well illustrated in allopolyploids, where WGD is concomitant with interspecific hybridization. This 'Genome Shock', usually accompanied by structural and functional modifications, has been associated with the activation of transposable elements (TEs). However, the impact of allopolyploidy on TEs has been studied in only a few polyploid species, and not in Brassica, which has been marked by recurrent polyploidy events. Here, we developed sequence-specific amplification polymorphism (SSAP) markers for three contrasting TEs, and compared profiles between resynthesized Brassica napus allotetraploids and their diploid Brassica progenitors. To evaluate restructuring at TE insertion sites, we scored changes in SSAP profiles and analysed a large set of differentially amplified SSAP bands. No massive structural changes associated with the three TEs surveyed were detected. However, several transposition events, specific to the youngest TE originating from the B. oleracea genome, were identified. Our study supports the hypothesis that TE responses to allopolyploidy are highly specific. The changes observed in SSAP profiles lead us to hypothesize that they may partly result from changes in DNA methylation, questioning the role of epigenetics during the formation of a new allopolyploid genome.

The dispensable chromosome of Leptosphaeria maculans shelters an effector gene conferring avirulence towards Brassica rapa.

Phytopathogenic fungi frequently contain dispensable chromosomes, some of which contribute to host range or pathogenicity. In Leptosphaeria maculans, the stem canker agent of oilseed rape (Brassica napus), the minichromosome was previously suggested to be dispensable, without evidence for any role in pathogenicity. Using genetic and genomic approaches, we investigated the inheritance and molecular determinant of an L. maculans-Brassica rapa incompatible interaction. Single gene control of the resistance was found, while all markers located on the L. maculans minichromosome, absent in the virulent parental isolate, co-segregated with the avirulent phenotype. Only one candidate avirulence gene was identified on the minichromosome, validated by complementation experiments and termed AvrLm11. The minichromosome was frequently lost following meiosis, but the frequency of isolates lacking it remained stable in field populations sampled at a 10-yr time interval, despite a yearly sexual stage in the L. maculans life cycle. This work led to the cloning of a new 'lost in the middle of nowhere' avirulence gene of L. maculans, interacting with a B. rapa resistance gene termed Rlm11 and introgressed into B. napus. It demonstrated the dispensability of the L. maculans minichromosome and suggested that its loss generates a fitness deficit.

Repeated polyploidy drove different levels of crossover suppression between homoeologous chromosomes in Brassica napus allohaploids.

Allopolyploid species contain more than two sets of related chromosomes (homoeologs) that must be sorted during meiosis to ensure fertility. As polyploid species usually have multiple origins, one intriguing, yet largely underexplored, question is whether different mechanisms suppressing crossovers between homoeologs may coexist within the same polyphyletic species. We addressed this question using Brassica napus, a young polyphyletic allopolyploid species. We first analyzed the meiotic behavior of 363 allohaploids produced from 29 accessions, which represent a large part of B. napus genetic diversity. Two main clear-cut meiotic phenotypes were observed, encompassing a twofold difference in the number of univalents at metaphase I. We then sequenced two chloroplast intergenic regions to gain insight into the maternal origins of the same 29 accessions; only two plastid haplotypes were found, and these correlated with the dichotomy of meiotic phenotypes. Finally, we analyzed genetic diversity at the PrBn locus, which was shown to determine meiotic behavior in a segregating population of B. napus allohaploids. We observed that segregation of two alleles at PrBn could adequately explain a large part of the variation in meiotic behavior found among B. napus allohaploids. Overall, our results suggest that repeated polyploidy resulted in different levels of crossover suppression between homoeologs in B. napus allohaploids.

Crossovers get a boost in Brassica allotriploid and allotetraploid hybrids.

Meiotic crossovers are necessary to generate balanced gametes and to increase genetic diversity. Even if crossover number is usually constrained, recent results suggest that manipulating karyotype composition could be a new way to increase crossover frequency in plants. In this study, we explored this hypothesis by analyzing the extent of crossover variation in a set of related diploid AA, allotriploid AAC, and allotetraploid AACC Brassica hybrids. We first used cytogenetic methods to describe the meiotic behavior of the different hybrids. We then combined a cytogenetic estimation of class I crossovers in the entire genome by immunolocalization of a key protein, MutL Homolog1, which forms distinct foci on meiotic chromosomes, with genetic analyses to specifically compare crossover rates between one pair of chromosomes in the different hybrids. Our results showed that the number of crossovers in the allotriploid AAC hybrid was higher than in the diploid AA hybrid. Accordingly, the allotetraploid AACC hybrid showed an intermediate behavior. We demonstrated that this increase was related to hybrid karyotype composition (diploid versus allotriploid versus allotetraploid) and that interference was maintained in the AAC hybrids. These results could provide another efficient way to manipulate recombination in traditional breeding and genetic studies.

Immediate unidirectional epigenetic reprogramming of NORs occurs independently of rDNA rearrangements in synthetic and natural forms of a polyploid species Brassica napus.

The dynamics of genome modification that occurred from the initial hybridization event to the stabilization of allopolyploid species remains largely unexplored. Here, we studied inheritance and expression of rDNA loci in the initial generations of Brassica napus allotetraploids (2n = 38, AACC) resynthesized from Brassica oleracea (2n = 18, CC) and B. rapa (2n = 20, AA) and compared the patterns to natural forms. Starting already from F1 generation, there was a strong uniparental silencing of B. oleracea genes. The epigenetic reprogramming was accompanied with immediate condensation of C-genome nucleolar organizer region (NOR) and progressive transgeneration hypermethylation of polymerase I promoters, mainly at CG sites. No such changes were observed in the A-genome NORs. Locus loss and gains affecting mainly non-NOR loci after the first allotetraploid meiosis did not influence established functional status of NORs. Collectively, epigenetic and genetic modifications in synthetic lines resemble events that accompanied formation of natural allopolyploid species.

Genome structure affects the rate of autosyndesis and allosyndesis in AABC, BBAC and CCAB Brassica interspecific hybrids.

Gene introgression into allopolyploid crop species from diploid or polyploid ancestors can be accomplished through homologous or homoeologous chromosome pairing during meiosis. We produced trigenomic Brassica interspecific hybrids (genome complements AABC, BBAC and CCAB) from the amphidiploid species Brassica napus (AACC), Brassica juncea (AABB) and Brassica carinata (BBCC) in order to test whether the structure of each genome affects frequencies of homologous and homoeologous (both allosyndetic and autosyndetic) pairing during meiosis. AABC hybrids produced from three genotypes of B. napus were included to assess the genetic control of homoeologous pairing. Multi-colour fluorescent in situ hybridisation was used to quantify homologous pairing (e.g. A-genome bivalents in AABC), allosyndetic associations (e.g. B-C in AABC) and autosyndetic associations (e.g. B-B in AABC) at meiosis. A high percentage of homologous chromosomes formed pairs (97.5-99.3%), although many pairs were also involved in autosyndetic and allosyndetic associations. Allosyndesis was observed most frequently as A-C genome associations (mean 4.0 per cell) and less frequently as A-B genome associations (0.8 per cell) and B-C genome associations (0.3 per cell). Autosyndesis occurred most frequently in the haploid A genome (0.75 A-A per cell) and least frequently in the haploid B genome (0.13 B-B per cell). The frequency of C-C autosyndesis was greater in BBAC hybrids (0.75 per cell) than in any other hybrid. The frequency of A-B, A-C and B-C allosyndesis was affected by the genomic structure of the trigenomic hybrids. Frequency of allosyndesis was also influenced by the genotype of the B. napus paternal parent for the three AABC (B. juncea × B. napus) hybrid types. Homoeologous pairing between the Brassica A, B and C genomes in interspecific hybrids may be influenced by complex interactions between genome structure and allelic composition.

A Modified Meiotic Recombination in Brassica napus Largely Improves Its Breeding Efficiency.

Meiotic recombination is the main tool used by breeders to generate biodiversity, allowing genetic reshuffling at each generation. It enables the accumulation of favorable alleles while purging deleterious mutations. However, this mechanism is highly regulated with the formation of one to rarely more than three crossovers, which are not randomly distributed. In this study, we showed that it is possible to modify these controls in oilseed rape (Brassica napus, AACC, 2n = 4x = 38) and that it is linked to AAC allotriploidy and not to polyploidy per se. To that purpose, we compared the frequency and the distribution of crossovers along A chromosomes from hybrids carrying exactly the same A nucleotide sequence, but presenting three different ploidy levels: AA, AAC and AACC. Genetic maps established with 202 SNPs anchored on reference genomes revealed that the crossover rate is 3.6-fold higher in the AAC allotriploid hybrids compared to AA and AACC hybrids. Using a higher SNP density, we demonstrated that smaller and numerous introgressions of B. rapa were present in AAC hybrids compared to AACC allotetraploid hybrids, with 7.6 Mb vs. 16.9 Mb on average and 21 B. rapa regions per plant vs. nine regions, respectively. Therefore, this boost of recombination is highly efficient to reduce the size of QTL carried in cold regions of the oilseed rape genome, as exemplified here for a QTL conferring blackleg resistance.

Quantitative resistance increases the durability of qualitative resistance to Leptosphaeria maculans in Brassica napus.

It has frequently been hypothesized that quantitative resistance increases the durability of qualitative (R-gene mediated) resistance but supporting experimental evidence is rare. To test this hypothesis, near-isogenic lines with/without the R-gene Rlm6 introduced into two Brassica napus cultivars differing in quantitative resistance to Leptosphaeria maculans were used in a 5-yr field experiment. Recurrent selection of natural fungal populations was done annually on each of the four plant genotypes, using crop residues from each genotype to inoculate separately the four series of field trials for five consecutive cropping seasons. Severity of phoma stem canker was measured on each genotype and frequencies of avirulence alleles in L. maculans populations were estimated. Recurrent selection of virulent isolates by Rlm6 in a susceptible background rendered the resistance ineffective by the third cropping season. By contrast, the resistance was still effective after 5 yr of selection by the genotype combining this gene with quantitative resistance. No significant variation in the performance of quantitative resistance alone was noted over the course of the experiment. We conclude that quantitative resistance can increase the durability of Rlm6. We recommend combining quantitative resistance with R-gene mediated resistance to enhance disease control and crop production.

Corrigendum: FANCM Limits Meiotic Crossovers in Brassica Crops.

[This corrects the article DOI: 10.3389/fpls.2018.00368.].

Homeologous recombination plays a major role in chromosome rearrangements that occur during meiosis of Brassica napus haploids.

Chromosomal rearrangements can be triggered by recombination between distinct but related regions. Brassica napus (AACC; 2n = 38) is a recent allopolyploid species whose progenitor genomes are widely replicated. In this article, we analyze the extent to which chromosomal rearrangements originate from homeologous recombination during meiosis of haploid B. napus (n = 19) by genotyping progenies of haploid x euploid B. napus with molecular markers. Our study focuses on three pairs of homeologous regions selected for their differing levels of divergence (N1/N11, N3/N13, and N9/N18). We show that a high number of chromosomal rearrangements occur during meiosis of B. napus haploid and are transmitted by first division restitution (FDR)-like unreduced gametes to their progeny; half of the progeny of Darmor-bzh haploids display duplications and/or losses in the chromosomal regions being studied. We demonstrate that half of these rearrangements are due to recombination between regions of primary homeology, which represents a 10- to 100-fold increase compared to the frequency of homeologous recombination measured in euploid lines. Some of the other rearrangements certainly result from recombination between paralogous regions because we observed an average of one to two autosyndetic A-A and/or C-C bivalents at metaphase I of the B. napus haploid. These results are discussed in the context of genome evolution of B. napus.

FANCM Limits Meiotic Crossovers in Brassica Crops.

Meiotic crossovers (COs) are essential for proper chromosome segregation and the reshuffling of alleles during meiosis. In WT plants, the number of COs is usually small, which limits the genetic variation that can be captured by plant breeding programs. Part of this limitation is imposed by proteins like FANCM, the inactivation of which results in a 3-fold increase in COs in Arabidopsis thaliana. Whether the same holds true in crops needed to be established. In this study, we identified EMS induced mutations in FANCM in two species of economic relevance within the genus Brassica. We showed that CO frequencies were increased in fancm mutants in both diploid and tetraploid Brassicas, Brassica rapa and Brassica napus respectively. In B. rapa, we observed a 3-fold increase in the number of COs, equal to the increase observed previously in Arabidopsis. In B. napus we observed a lesser but consistent increase (1.3-fold) in both euploid (AACC) and allohaploid (AC) plants. Complementation tests in A. thaliana suggest that the smaller increase in crossover frequency observed in B. napus reflects residual activity of the mutant C copy of FANCM. Altogether our results indicate that the anti-CO activity of FANCM is conserved across the Brassica, opening new avenues to make a wider range of genetic diversity accessible to crop improvement.

Numerous and rapid nonstochastic modifications of gene products in newly synthesized Brassica napus allotetraploids.

Polyploidization is a widespread process that results in the merger of two or more genomes in a common nucleus. To investigate modifications of gene expression occurring during allopolyploid formation, the Brassica napus allotetraploid model was chosen. Large-scale analyses of the proteome were conducted on two organs, the stem and root, so that >1600 polypeptides were screened. Comparative proteomics of synthetic B. napus and its homozygous diploid progenitors B. rapa and B. oleracea showed that very few proteins disappeared or appeared in the amphiploids (<1%), but a strikingly high number (25-38%) of polypeptides displayed quantitative nonadditive pattern. Nonstochastic gene expression repatterning was found since 99% of the detected variations were reproducible in four independently created amphiploids. More than 60% of proteins displayed a nonadditive pattern closer to the paternal parent B. rapa. Interspecific hybridization triggered the majority of the deviations (89%), whereas very few variations (approximately 3%) were associated with genome doubling and more significant alterations arose from selfing (approximately 9%). Some nonadditive proteins behaved similarly in both organs, while others exhibited contrasted behavior, showing rapid organ-specific regulation. B. napus formation was therefore correlated with immediate and directed nonadditive changes in gene expression, suggesting that the early steps of allopolyploidization repatterning are controlled by nonstochastic mechanisms.

Mapping PrBn and other quantitative trait loci responsible for the control of homeologous chromosome pairing in oilseed rape (Brassica napus L.) haploids.

In allopolyploid species, fair meiosis could be challenged by homeologous chromosome pairing and is usually achieved by the action of homeologous pairing suppressor genes. Oilseed rape (Brassica napus) haploids (AC, n=19) represent an attractive model for studying the mechanisms used by allopolyploids to ensure the diploid-like meiotic pairing pattern. In oilseed rape haploids, homeologous chromosome pairing at metaphase I was found to be genetically based and controlled by a major gene, PrBn, segregating in a background of polygenic variation. In this study, we have mapped PrBn within a 10-cM interval on the C genome linkage group DY15 and shown that PrBn displays incomplete penetrance or variable expressivity. We have identified three to six minor QTL/BTL that have slight additive effects on the amount of pairing at metaphase I but do not interact with PrBn. We have also detected a number of other loci that interact epistatically, notably with PrBn. Our results support the idea that, as in other polyploid species, metaphase I homeologous pairing in oilseed rape haploids is controlled by an integrated system of several genes, which function in a complex manner.