Centromere Locations in Brassica A and C Genomes Revealed Through Half-Tetrad Analysis

The evolution and formation of centromeric repeats analysis in Vitis vinifera

MAIN CONCLUSION

Six grape centromere-specific markers for cytogenetics were mined by combining genetic and immunological assays, and the possible evolution mechanism of centromeric repeats was analyzed. Centromeric histone proteins are functionally conserved; however, centromeric repetitive DNA sequences may represent considerable diversity in related species. Therefore, studying the characteristics and structure of grape centromere repeat sequences contributes to a deeper understanding of the evolutionary process of grape plants, including their origin and mechanisms of polyploidization. Plant centromeric regions are mainly composed of repetitive sequences, including SatDNA and transposable elements (TE). In this research, the characterization of centromere sequences in the whole genome of grapevine (Vitis vinifera L.) has been conducted. Five centromeric tandem repeat sequences (Vv1, Vv2, Vv5, Vv6, and Vv8) and one long terminal repeat (LTR) sequence Vv24 were isolated. These sequences had different centromeric distributions, which indicates that grape centromeric sequences may undergo rapid evolution. The existence of extrachromosomal circular DNA (eccDNA) and gene expression in CenH3 subdomain region may provide various potential mechanisms for the generation of new centromeric regions.

Heterozygosity analysis of spontaneous 2n female gametes and centromere mapping of the diploid Hevea brasiliensis based on full-sib triploid populations

KEY MESSAGE

Unreduced megagametophytes via second-division restitution were confirmed through heterozygosity analysis, and four candidate physical centromeres of rubber were located for the first time. The evaluation of maternal heterozygosity restitution (MHR) is vital in identifying the mechanism of 2n gametogenesis and assessing the utilization value of 2n gametes. In this study, three full-sib triploid populations were employed to evaluate the MHR of 2n female gametes of rubber tree clone GT1 and to confirm their genetic derivation. The 2n female gametes of GT1 were derived from second-division restitution (SDR) and transmitted more than half of the parental heterozygosity. In addition, low recombination frequency markers were developed, and four candidate physical centromeres of rubber tree were located for the first time. The confirmation that 2n female gametes of rubber tree clone GT1 are derived from SDR provides insights into the molecular mechanisms of 2n gametogenesis. In addition, the identified centromere location will aid in the development of centromeric markers for the rapid identification of the 2n gametogenesis mechanism.

Genome composition in Brassica interspecific hybrids affects chromosome inheritance and viability of progeny

Interspecific hybridization is widespread in nature and can result in the formation of new hybrid species as well as the transfer of traits between species. However, the fate of newly formed hybrid lineages is relatively understudied. We undertook pairwise crossing between multiple genotypes of three Brassica allotetraploid species Brassica juncea (2n = AABB), Brassica carinata (2n = BBCC), and Brassica napus (2n = AACC) to generate AABC, BBAC, and CCAB interspecific hybrids and investigated chromosome inheritance and fertility in these hybrids and their self-pollinated progeny. Surprisingly, despite the presence of a complete diploid genome in all hybrids, hybrid fertility was very low. AABC and BBAC first generation (F1) hybrids both averaged ~16% pollen viability compared to 3.5% in CCAB hybrids: most CCAB hybrid flowers were male-sterile. AABC and CCAB F1 hybrid plants averaged 5.5 and 0.5 seeds per plant, respectively, and BBAC F1 hybrids ~56 seeds/plant. In the second generation (S1), all confirmed self-pollinated progeny resulting from CCAB hybrids were sterile, producing no self-pollinated seeds. Three AABC S1 hybrids putatively resulting from unreduced gametes produced 3, 14, and 182 seeds each, while other AABC S1 hybrids averaged 1.5 seeds/plant (0-8). BBAC S1 hybrids averaged 44 seeds/plant (range 0-403). We also observed strong bias towards retention rather than loss of the haploid genomes, suggesting that the subgenomes in the Brassica allotetraploids are already highly interdependent, such that loss of one subgenome is detrimental to fertility and viability. Our results suggest that relationships between subgenomes determine hybridization outcomes in these species.

Meta-analysis of seed weight QTLome using a consensus and highly dense genetic map in Brassica napus L

KEY MESSAGE

We report here the discovery of high-confidence MQTL regions and of putative candidate genes associated with seed weight in B. napus using a highly dense consensus genetic map and by comparing various large-scale multiomics datasets. Seed weight (SW) is a direct determinant of seed yield in Brassica napus and is controlled by many loci. To unravel the main genomic regions associated with this complex trait, we used 13 available genetic maps to construct a consensus and highly dense map, comprising 40,401 polymorphic markers and 9191 genetic bins, harboring a cumulative length of 3047.8 cM. Then, we performed a meta-analysis using 639 projected SW quantitative trait loci (QTLs) obtained from studies conducted since 1999, enabling the identification of 57 meta-QTLS (MQTLs). The confidence intervals of our MQTLs were 9.8 and 4.3 times lower than the average CIs of the original QTLs for the A and C subgenomes, respectively, resulting in the detection of some key genes and several putative novel candidate genes associated with SW. By comparing the genes identified in MQTL intervals with multiomics datasets and coexpression analyses of common genes, we defined a more reliable and shorter list of putative candidate genes potentially involved in the regulation of seed maturation and SW. As an example, we provide a list of promising genes with high expression levels in seeds and embryos (e.g., BnaA03g04230D, BnaC03g08840D, BnaA10g29580D and BnaA03g27410D) that can be more finely studied through functional genetics experiments or that may be useful for MQTL-assisted breeding for SW. The high-density genetic consensus map and the single nucleotide polymorphism (SNP) physical map generated from the latest B. napus cv. Darmor-bzh v10 assembly will be a valuable resource for further mapping and map-based cloning of other important traits.

Genome sequence analysis provides insights into the mode of 2n egg formation in Solanum malmeanum

KEY MESSAGE

Our genomic investigation confirms the mechanism of 2n eggs formation in S. malmeanum and aid in optimizing the use of wild germplasm. Wild potatoes are a valuable source of agronomic traits. However, substantial reproductive barriers limit gene flow into cultivated species. 2n gametes are instrumental in preventing endosperm abortion caused by genetic imbalances in the endosperm. However, little is known about the molecular mechanisms underlying the formation of 2n gametes. Here, the wild species Solanum malmeanum Bitter (2x, 1EBN, endosperm balance number) was used in inter- and intrapoloid crosses with other Solanum species, with viable seeds being produced only when S. malmeanum was used as the female parent to cross the 2EBN Solanum genus and with the likely involvement of 2n gametes. Subsequently, we substantiated the formation of 2n eggs in S. malmeanum using fluorescence in situ hybridization (FISH) and genomic sequencing technology. Additionally, the transmission rate of maternal heterozygous polymorphism sites was assessed from a genomic perspective to analyze the mode of 2n egg formation in S. malmeanum × S. tuberosum and S. malmeanum × S. chacoense crosses; each cross acquired an average of 31.12% and 22.79% maternal sites, respectively. This confirmed that 2n egg formation in S. malmeanum attributed to second-division restitution (SDR) coupled with the occurrence of exchange events. The high-throughput sequencing technology used in this study has strong advantages over traditional cytological analyses. Furthermore, S. malmeanum, which has a variety of excellent traits not available from present cultivated potato genepool, has received little research attention and has successfully achieved gene flow in cultivated species in the current study. These findings will facilitate the understanding and optimization of wild germplasm utilization in potatoes.

Novel gene loci associated with susceptibility or cryptic quantitative resistance to Pyrenopeziza brassicae in Brassica napus

Quantitative disease resistance (QDR) controls the association of the light leaf spot pathogen with Brassica napus; four QDR loci that were in linkage disequilibrium and eight gene expression markers were identified. Quantitative disease resistance (QDR) can provide durable control of pathogens in crops in contrast to resistance (R) gene-mediated resistance which can break down due to pathogen evolution. QDR is therefore a desirable trait in crop improvement, but little is known about the causative genes, and so it is difficult to incorporate into breeding programmes. Light leaf spot, caused by Pyrenopeziza brassicae, is an important disease of oilseed rape (canola, Brassica napus). To identify new QDR gene loci, we used a high-throughput screening pathosystem with P. brassicae on 195 lines of B. napus combined with an association transcriptomics platform. We show that all resistance against P. brassicae was associated with QDR and not R gene-mediated. We used genome-wide association analysis with an improved B. napus population structure to reveal four gene loci significantly (P = 0.0001) associated with QDR in regions showing linkage disequilibrium. On chromosome A09, enhanced resistance was associated with heterozygosity for a cytochrome P450 gene co-localising with a previously described locus for seed glucosinolate content. In addition, eight significant gene expression markers with a false discovery rate of 0.001 were associated with QDR against P. brassicae. For seven of these, expression was positively correlated with resistance, whereas for one, a HXXXD-type acyl-transferase, negative correlation indicated a potential susceptibility gene. The study identifies novel QDR loci for susceptibility and resistance, including novel cryptic QDR genes associated with heterozygosity, that will inform future crop improvement.

Frequent spontaneous structural rearrangements promote rapid genome diversification in a Brassica napus F1 generation

In a cross between two homozygous Brassica napus plants of synthetic and natural origin, we demonstrate that novel structural genome variants from the synthetic parent cause immediate genome diversification among F1 offspring. Long read sequencing in twelve F1 sister plants revealed five large-scale structural rearrangements where both parents carried different homozygous alleles but the heterozygous F1 genomes were not identical heterozygotes as expected. Such spontaneous rearrangements were part of homoeologous exchanges or segmental deletions and were identified in different, individual F1 plants. The variants caused deletions, gene copy-number variations, diverging methylation patterns and other structural changes in large numbers of genes and may have been causal for unexpected phenotypic variation between individual F1 sister plants, for example strong divergence of plant height and leaf area. This example supports the hypothesis that spontaneous de novo structural rearrangements after de novo polyploidization can rapidly overcome intense allopolyploidization bottlenecks to re-expand crops genetic diversity for ecogeographical expansion and human selection. The findings imply that natural genome restructuring in allopolyploid plants from interspecific hybridization, a common approach in plant breeding, can have a considerably more drastic impact on genetic diversity in agricultural ecosystems than extremely precise, biotechnological genome modifications.

R gene triplication confers European fodder turnip with improved clubroot resistance

Clubroot is one of the most important diseases for many important cruciferous vegetables and oilseed crops worldwide. Different clubroot resistance (CR) loci have been identified from only limited species in Brassica, making it difficult to compare and utilize these loci. European fodder turnip ECD04 is considered one of the most valuable resources for CR breeding. To explore the genetic and evolutionary basis of CR in ECD04, we sequenced the genome of ECD04 using de novo assembly and identified 978 candidate R genes. Subsequently, the 28 published CR loci were physically mapped to 15 loci in the ECD04 genome, including 62 candidate CR genes. Among them, two CR genes, CRA3.7.1 and CRA8.2.4, were functionally validated. Phylogenetic analysis revealed that CRA3.7.1 and CRA8.2.4 originated from a common ancestor before the whole-genome triplication (WGT) event. In clubroot susceptible Brassica species, CR-gene homologues were affected by transposable element (TE) insertion, resulting in the loss of CR function. It can be concluded that the current functional CR genes in Brassica rapa and non-functional CR genes in other Brassica species were derived from a common ancestral gene before WGT. Finally, a hypothesis for CR gene evolution is proposed for further discussion.

Modelling of gene loss propensity in the pangenomes of three Brassica species suggests different mechanisms between polyploids and diploids

Plant genomes demonstrate significant presence/absence variation (PAV) within a species; however, the factors that lead to this variation have not been studied systematically in Brassica across diploids and polyploids. Here, we developed pangenomes of polyploid Brassica napus and its two diploid progenitor genomes B. rapa and B. oleracea to infer how PAV may differ between diploids and polyploids. Modelling of gene loss suggests that loss propensity is primarily associated with transposable elements in the diploids while in B. napus, gene loss propensity is associated with homoeologous recombination. We use these results to gain insights into the different causes of gene loss, both in diploids and following polyploidization, and pave the way for the application of machine learning methods to understanding the underlying biological and physical causes of gene presence/absence.

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

Simple Summary

The selection of varieties more resilient to disease and climate change requires generating new genetic diversity for breeding. The main mechanism for reshuffling genetic information is through the recombination of chromosomes during meiosis. We showed in oilseed rape (Brassica napus, AACC, 2n = 4x = 38), which is a natural hybrid formed from a cross between turnip (B. rapa, AA, 2n = 2x = 20) and cabbage (B. oleracea, CC, 2n = 2x = 18), that there is significantly more crossovers occurring along the entire A chromosomes in allotriploid AAC (crossbetween B. napus and B. rapa) than in diploid AA or allotetraploid AACC hybrids. We demonstrated that these allotriploid AAC hybrids are highly efficient to introduce new variability within oilseed rape varieties, notably by enabling the introduction of small genomic regions carrying genes controlling agronomically interesting traits.

Abstract

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.

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 S₁ to S₉ 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.

Stable, fertile lines produced by hybridization between allotetraploids Brassica juncea (AABB) and Brassica carinata (BBCC) have merged the A and C genomes

Many flowering plant taxa contain allopolyploids that share one or more genomes in common. In the Brassica genus, crop species Brassica juncea and Brassica carinata share the B genome, with 2n = AABB and 2n = BBCC genome complements, respectively. Hybridization results in 2n = BBAC hybrids, but the fate of these hybrids over generations of self-pollination has never been reported. We produced and characterized B. juncea × B. carinata (2n = BBAC) interspecific hybrids over six generations of self-pollination under selection for high fertility using a combination of genotyping, fertility phenotyping, and cytogenetics techniques. Meiotic pairing behaviour improved from 68% bivalents in the F₁ to 98% in the S₅ /S₆ generations, and initially low hybrid fertility also increased to parent species levels. The S₅ /S₆ hybrids contained an intact B genome (16 chromosomes) plus a new, stable A/C genome (18-20 chromosomes) resulting from recombination and restructuring of A and C-genome chromosomes. Our results provide the first experimental evidence that two genomes can come together to form a new, restructured genome in hybridization events between two allotetraploid species that share a common genome. This mechanism should be considered in interpreting phylogenies in taxa with multiple allopolyploid species.

Enhancing backcross programs through increased recombination

Background

Introgression of a quantitative trait locus (QTL) by successive backcrosses is used to improve elite lines (recurrent parent) by introducing alleles from exotic material (donor parent). In the absence of selection, the proportion of the donor genome decreases by half at each generation. However, since selection is for the donor allele at the QTL, elimination of the donor genome around that QTL will be much slower than in the rest of the genome (i.e. linkage drag). Using markers to monitor the genome around the QTL and in the genetic background can accelerate the return to the recurrent parent genome. Successful introgression of a locus depends partly on the occurrence of crossovers at favorable positions. However, the number of crossovers per generation is limited and their distribution along the genome is heterogeneous. Recently, techniques have been developed to modify these two recombination parameters.

Results

In this paper, we assess, by simulations in the context of Brassicaceae, the effect of increased recombination on the efficiency of introgression programs by studying the decrease in linkage drag and the recovery of the recurrent genome. The simulated selection schemes begin by two generations of foreground selection and continue with one or more generations of background selection. Our results show that, when the QTL is in a region that initially lacked crossovers, an increase in recombination rate can decrease linkage drag by nearly ten-fold after the foreground selection and improves the return to the recurrent parent. However, if the QTL is in a region that is already rich in crossovers, an increase in recombination rate is detrimental.

Conclusions

Depending on the recombination rate in the region targeted for introgression, increasing it can be beneficial or detrimental. Thus, the simulations analysed in this paper help us understand how an increase in recombination rate can be beneficial. They also highlight the best methods that can be used to increase recombination rate, depending on the situation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00619-0.

Long-read assembly of the Brassica napus reference genome Darmor-bzh

BACKGROUND

The combination of long reads and long-range information to produce genome assemblies is now accepted as a common standard. This strategy not only allows access to the gene catalogue of a given species but also reveals the architecture and organization of chromosomes, including complex regions such as telomeres and centromeres. The Brassica genus is not exempt, and many assemblies based on long reads are now available. The reference genome for Brassica napus, Darmor-bzh, which was published in 2014, was produced using short reads and its contiguity was extremely low compared with current assemblies of the Brassica genus.

FINDINGS

Herein, we report the new long-read assembly of Darmor-bzh genome (Brassica napus) generated by combining long-read sequencing data and optical and genetic maps. Using the PromethION device and 6 flowcells, we generated ∼16 million long reads representing 93× coverage and, more importantly, 6× with reads longer than 100 kb. This ultralong-read dataset allows us to generate one of the most contiguous and complete assemblies of a Brassica genome to date (contig N50 > 10 Mb). In addition, we exploited all the advantages of the nanopore technology to detect modified bases and sequence transcriptomic data using direct RNA to annotate the genome and focus on resistance genes.

CONCLUSION

Using these cutting-edge technologies, and in particular by relying on all the advantages of the nanopore technology, we provide the most contiguous Brassica napus assembly, a resource that will be valuable to the Brassica community for crop improvement and will facilitate the rapid selection of agronomically important traits.

Localization and characterization of Citrus centromeres by combining half-tetrad analysis and CenH3-associated sequence profiling

The physical locations of citrus centromere are revealed by combining genetic and immunological assays for the first time and nine citrus centromere-specific markers for cytogenetics are mined. Centromere localization is challenging, because highly redundant repetitive sequences in centromeric regions make sequence assembly difficult. Although several citrus genomes have been released, the centromeric regions and their characteristics remain to be elucidated. Here, we mapped citrus centromeres through half-tetrad analysis (HTA) that included the genotyping of 54 tetraploid hybrids derived from 2n megagametophytes of Nadorcott tangor with 212 single nucleotide polymorphism (SNP) markers. The sizes of centromeric regions, which estimated based on the heterozygosity restitution rate pattern along the chromosomes, ranged from 1.12 to 18.19 Mb. We also profiled the binding sequences with the centromere-specific histone variant CenH3 by chromatin immunoprecipitation sequencing (ChIP-seq). Based on the positions of the top ten CenH3-enriched contigs, the sizes of centromeric regions were estimated to range from 0.01 to 7.60 Mb and were either adjacent to or included in the centromeric regions identified by HTA. We used DNA probes from two repeats selected from the centromeric regions and seven CenH3-binding centromeric repeats to verify centromeric locations by fluorescence in situ hybridization (FISH). Centromere localization in citrus will contribute to the mining of centromeric/pericentromeric markers, thus to facilitate the rapid identification of mechanisms underlying 2n gamete formation and serve the polyploidy breeding.

Assessing by Modeling the Consequences of Increased Recombination in Recurrent Selection of Oryza sativa and Brassica rapa

Meiotic recombination generates genetic diversity but in most species the number of crossovers per meiosis is limited. Previous modeling studies showed that increasing recombination can enhance response to selection. However, such studies did not assume a specific method of modifying recombination. Our objective was to test whether two methods used to increase recombination in plants could increase genetic gain in a population undergoing recurrent selection such as in genomic selection programs. The first method, in Oryza sativa, used a mutant of anti-crossover genes, increasing global recombination without affecting the recombination landscape shape. The second one used the ploidy level of a cross between Brassica rapa and Brassica napus, increasing recombination especially in pericentromeric regions. Our modeling framework used these recombination landscapes and sampled quantitative trait loci positions from the actual gene distributions. We simulated selection programs with initially a cross between two inbred lines, for two species. Increased recombination enhanced the response to selection. The amount of enhancement in the cumulative gain largely depended on the species and the number of quantitative trait loci (2, 10, 20, 50, 200 or 1000 per chromosome). Genetic gains were increased up to 30% after 20 generations. Furthermore, increasing recombination in cold regions was the most effective: the gain was larger by 25% with the first method and 34% with the second one in B. rapa, and 12% compared to 16% in O. sativa In summary, increased recombination enhances the genetic gain in long-term selection programs, with visible effects after four to five generations.

A quartet pollen phenotype identified in a population of Brassica interspecific hybrids shows incomplete penetrance and variable response to temperature

Quartet pollen, where pollen grains remain attached to each other post-meiosis, is useful for tetrad analysis, crossover assessment and centromere mapping. We observed the quartet pollen phenotype for the first time in the agriculturally significant Brassica genus, in an experimental population of allohexaploid Brassica hybrids derived from the cross (Brassica napus × B. carinata) × B. juncea followed by two self-pollination generations. Quartet pollen production was assessed in 144 genotypes under glasshouse conditions, following which a set of 16 genotypes were selected to further investigate the effect of environment (warm: 25 °C and cold: 10 °C temperatures) on quartet pollen production in growth cabinets. Under glasshouse phenotyping conditions, only 92 out of 144 genotypes produced enough pollen to score: of these, 30 did not produce any observable quartet pollen, while 62 genotypes produced quartet pollen at varying frequencies. Quartet pollen production appeared quantitative and did not clearly fall into phenotypic or qualitative categories indicative of major gene expression. No consistent effect of temperature on quartet pollen production was identified, with some genotypes producing more and some producing less quartet pollen under different temperature treatments. The genetic heterogeneity and frequent pollen infertility of this population prevents strong conclusions being made. However, it is clear that the quartet phenotype in this Brassica population does not show complete penetrance and shows variable (likely genotype-specific) response to temperature stress. In future, identification of quartet phenotypes in Brassica would perhaps best be carried out via screening of diploid (e.g. B. rapa) TILLING populations.

Genetic changes in a novel breeding population of Brassica napus synthesized from hundreds of crosses between B. rapa and B. carinata

Introgression of genomic variation between and within related crop species is a significant evolutionary approach for population differentiation, genome reorganization and trait improvement. Using the Illumina Infinium Brassica 60K SNP array, we investigated genomic changes in a panel of advanced generation new-type Brassica napus breeding lines developed from hundreds of interspecific crosses between 122 Brassica rapa and 74 Brassica carinata accessions, and compared them with representative accessions of their three parental species. The new-type B. napus population presented rich genetic diversity and abundant novel genomic alterations, consisting of introgressions from B. rapa and B. carinata, novel allelic combinations, reconstructed linkage disequilibrium patterns and haplotype blocks, and frequent deletions and duplications (nonrandomly distributed), particularly in the C subgenome. After a much shorter, but very intensive, selection history compared to traditional B. napus, a total of 15 genomic regions with strong selective sweeps and 112 genomic regions with putative signals of selective sweeps were identified. Some of these regions were associated with important agronomic traits that were selected for during the breeding process, while others were potentially associated with restoration of genome stability and fertility after interspecific hybridization. Our results demonstrate how a novel method for population-based crop genetic improvement can lead to rapid adaptation, restoration of genome stability and positive responses to artificial selection.

Genetic Properties of a Nested Association Mapping Population Constructed With Semi-Winter and Spring Oilseed Rapes

Nested association mapping (NAM) populations have been widely applied to dissect the genetic basis of complex quantitative traits in a variety of crops. In this study, we developed a Brassica napus NAM (BN-NAM) population consisting of 15 recombination inbred line (RIL) families with 2,425 immortal genotypes. Fifteen high-density genetic linkage maps were constructed by genotyping by sequencing (GBS) based on all RIL families, with further integration into a joint linkage map (JLM) having 30,209 unique markers in common with multiple linkage maps. Furthermore, an ultra-density whole-genome variation map was constructed by projecting 4,444,309 high-quality variants onto the JLM. The NAM population captured a total of 88,542 recombination events (REs). The uneven distribution of recombination rate along chromosomes is positively correlated with the densities of genes and markers, but negatively correlated with the density of transposable elements and linkage disequilibrium (LD). Analyses of population structure and principal components revealed that the BN-NAM population could be divided into three groups with weak stratification. The LD decay distance across genome varied between 170 and 2,400 Kb, with LD decay more rapid in the A than in the C sub-genome. The pericentromeric regions contained large LD blocks, especially in the C sub-genome. This NAM population provides a valuable resource for dissecting the genetic basis of important traits in rapeseed, especially in semi-winter oilseed rape.

Mapping of homoeologous chromosome exchanges influencing quantitative trait variation in Brassica napus

Genomic rearrangements arising during polyploidization are an important source of genetic and phenotypic variation in the recent allopolyploid crop Brassica napus. Exchanges among homoeologous chromosomes, due to interhomoeologue pairing, and deletions without compensating homoeologous duplications are observed in both natural B. napus and synthetic B. napus. Rearrangements of large or small chromosome segments induce gene copy number variation (CNV) and can potentially cause phenotypic changes. Unfortunately, complex genome restructuring is difficult to deal with in linkage mapping studies. Here, we demonstrate how high-density genetic mapping with codominant, physically anchored SNP markers can detect segmental homoeologous exchanges (HE) as well as deletions and accurately link these to QTL. We validated rearrangements detected in genetic mapping data by whole-genome resequencing of parental lines along with cytogenetic analysis using fluorescence in situ hybridization with bacterial artificial chromosome probes (BAC-FISH) coupled with PCR using primers specific to the rearranged region. Using a well-known QTL region influencing seed quality traits as an example, we confirmed that HE underlies the trait variation in a DH population involving a synthetic B. napus trait donor, and succeeded in narrowing the QTL to a small defined interval that enables delineation of key candidate genes.

Gene Introgression in Weeds Depends on Initial Gene Location in the Crop: Brassica napus-Raphanus raphanistrum Model

The effect of gene location within a crop genome on its transfer to a weed genome remains an open question for gene flow assessment. To elucidate this question, we analyzed advanced generations of intergeneric hybrids, derived from an initial pollination of known oilseed rape varieties (Brassica napus, AACC, 2n = 38) by a local population of wild radish (Raphanus raphanistrum, RrRr, 2n = 18). After five generations of recurrent pollination, 307 G5 plants with a chromosome number similar to wild radish were genotyped using 105 B. napus specific markers well distributed along the chromosomes. They revealed that 49.8% of G5 plants carried at least one B. napus genomic region. According to the frequency of B. napus markers (0-28%), four classes were defined: Class 1 (near zero frequency), with 75 markers covering ∼70% of oilseed rape genome; Class 2 (low frequency), with 20 markers located on 11 genomic regions; Class 3 (high frequency), with eight markers on three genomic regions; and Class 4 (higher frequency), with two adjacent markers detected on A10. Therefore, some regions of the oilseed rape genome are more prone than others to be introgressed into wild radish. Inheritance and growth of plant progeny revealed that genomic regions of oilseed rape could be stably introduced into wild radish and variably impact the plant fitness (plant height and seed number). Our results pinpoint that novel technologies enabling the targeted insertion of transgenes should select genomic regions that are less likely to be introgressed into the weed genome, thereby reducing gene flow.

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.

Segregation for fertility and meiotic stability in novel Brassica allohexaploids

Allohexaploid Brassica populations reveal ongoing segregation for fertility, while genotype influences fertility and meiotic stability. Creation of a new Brassica allohexaploid species is of interest for the development of a crop type with increased heterosis and adaptability. At present, no naturally occurring, meiotically stable Brassica allohexaploid exists, with little data available on chromosome behaviour and meiotic control in allohexaploid germplasm. In this study, 100 plants from the cross B. carinata × B. rapa (A2 allohexaploid population) and 69 plants from the cross (B. napus × B. carinata) × B. juncea (H2 allohexaploid population) were assessed for fertility and meiotic behaviour. Estimated pollen viability, self-pollinated seed set, number of seeds on the main shoot, number of pods on the main shoot, seeds per ten pods and plant height were measured for both the A2 and H2 populations and for a set of reference control cultivars. The H2 population had high segregation for pollen viability and meiotic stability, while the A2 population was characterised by low pollen fertility and a high level of chromosome loss. Both populations were taller, but had lower average fertility trait values than the control cultivar samples. The study also characterises fertility and meiotic chromosome behaviour in genotypes and progeny sets in heterozygous allotetraploid Brassica derived lines, and indicates that genotypes of the parents and H1 hybrids are affecting chromosome pairing and fertility phenotypes in the H2 population. The identification and characterisation of factors influencing stability in novel allohexaploid Brassica populations will assist in the development of this as a new crop species for food and agricultural benefit.

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.

The Impact of Open Pollination on the Structural Evolutionary Dynamics, Meiotic Behavior, and Fertility of Resynthesized Allotetraploid Brassica napus L

Allopolyploidy, which results from the merger and duplication of two divergent genomes, has played a major role in the evolution and diversification of flowering plants. The genomic changes that occur in resynthesized or natural neopolyploids have been extensively studied, but little is known about the effects of the reproductive mode in the initial generations that may precede its successful establishment. To truly reflect the early generations of a nascent polyploid, two resynthesized allotetraploid Brassica napus populations were obtained for the first time by open pollination. In these populations, we detected a much lower level of aneuploidy (third generation) compared with those previously published populations obtained by controlled successive selfing. We specifically studied 33 resynthesized B. napus individuals from our two open pollinated populations, and showed that meiosis was affected in both populations. Their genomes were deeply shuffled after allopolyploidization: up to 8.5 and 3.5% of the C and A subgenomes were deleted in only two generations. The identified deletions occurred mainly at the distal part of the chromosome, and to a significantly greater extent on the C rather than the A subgenome. Using Fluorescent In Situ Hybridization (BAC-FISH), we demonstrated that four of these deletions corresponded to fixed translocations (via homeologous exchanges). We were able to evaluate the size of the structural variations and their impact on the whole genome size, gene content, and allelic diversity. In addition, the evolution of fertility was assessed, to better understand the difficulty encountered by novel polyploid individuals before the putative formation of a novel stable species.

Genome-Wide Association Mapping Reveals the Genetic Control Underlying Branch Angle in Rapeseed (Brassica napus L.)

Plant architecture is vital not only for crop yield, but also for field management, such as mechanical harvesting. The branch angle is one of the key factors determining plant architecture. With the aim of revealing the genetic control underlying branch angle in rapeseed (Brassica napus L.), the positional variation of branch angles on individual plants was evaluated, and the branch angle increased with the elevation of branch position. Furthermore, three middle branches of individual plants were selected to measure the branch angle because they exhibited the most representative phenotypic values. An association panel with 472 diverse accessions was estimated for branch angle trait in six environments and genotyped with a 60K Brassica Infinium® SNP array. As a result of association mapping, 46 and 38 significantly-associated loci were detected using a mixed linear model (MLM) and a multi-locus random-SNP-effect mixed linear model (MRMLM), which explained up to 62.2 and 66.2% of the cumulative phenotypic variation, respectively. Numerous highly-promising candidate genes were identified by annotating against Arabidopsis thaliana homologous, including some first found in rapeseed, such as TAC1, SGR1, SGR3, and SGR5. These findings reveal the genetic control underlying branch angle and provide insight into genetic improvements that are possible in the plant architecture of rapeseed.

Population Genomic Analysis Reveals Differential Evolutionary Histories and Patterns of Diversity across Subgenomes and Subpopulations of Brassica napus L

The allotetraploid species Brassica napus L. is a global crop of major economic importance, providing canola oil (seed) and vegetables for human consumption and fodder and meal for livestock feed. Characterizing the genetic diversity present in the extant germplasm pool of B. napus is fundamental to better conserve, manage and utilize the genetic resources of this species. We used sequence-based genotyping to identify and genotype 30,881 SNPs in a diversity panel of 782 B. napus accessions, representing samples of winter and spring growth habits originating from 33 countries across Europe, Asia, and America. We detected strong population structure broadly concordant with growth habit and geography, and identified three major genetic groups: spring (SP), winter Europe (WE), and winter Asia (WA). Subpopulation-specific polymorphism patterns suggest enriched genetic diversity within the WA group and a smaller effective breeding population for the SP group compared to WE. Interestingly, the two subgenomes of B. napus appear to have different geographic origins, with phylogenetic analysis placing WE and WA as basal clades for the other subpopulations in the C and A subgenomes, respectively. Finally, we identified 16 genomic regions where the patterns of diversity differed markedly from the genome-wide average, several of which are suggestive of genomic inversions. The results obtained in this study constitute a valuable resource for worldwide breeding efforts and the genetic dissection and prediction of complex B. napus traits.

Genome-Wide Association Study Provides Insight into the Genetic Control of Plant Height in Rapeseed (Brassica napus L.)

Plant height is a key morphological trait of rapeseed. In this study, we measured plant height of a rapeseed population across six environments. This population contains 476 inbred lines representing the major Chinese rapeseed genepool and 44 lines from other countries. The 60K Brassica Infinium® SNP array was utilized to genotype the association panel. A genome-wide association study (GWAS) was performed via three methods, including a robust, novel, nonparametric Anderson-Darling (A-D) test. Consequently, 68 loci were identified as significantly associated with plant height (P < 5.22 × 10(-5)), and more than 70% of the loci (48) overlapped the confidence intervals of reported QTLs from nine mapping populations. Moreover, 24 GWAS loci were detected with selective sweep signals, which reflected the signatures of historical semi-dwarf breeding. In the linkage disequilibrium (LD) decay range up-and downstream of 65 loci (r (2) > 0.1), we found plausible candidates orthologous to the documented Arabidopsis genes involved in height regulation. One significant association found by GWAS colocalized with the established height locus BnRGA in rapeseed. Our results provide insights into the genetic basis of plant height in rapeseed and may facilitate marker-based breeding.