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

Dosage-sensitivity shapes how genes transcriptionally respond to allopolyploidy and homoeologous exchange in resynthesized Brassica napus

The gene balance hypothesis proposes that selection acts on the dosage (i.e. copy number) of genes within dosage-sensitive portions of networks, pathways, and protein complexes to maintain balanced stoichiometry of interacting proteins, because perturbations to stoichiometric balance can result in reduced fitness. This selection has been called dosage balance selection. Dosage balance selection is also hypothesized to constrain expression responses to dosage changes, making dosage-sensitive genes (those encoding members of interacting proteins) experience more similar expression changes. In allopolyploids, where whole-genome duplication involves hybridization of diverged lineages, organisms often experience homoeologous exchanges that recombine, duplicate, and delete homoeologous regions of the genome and alter the expression of homoeologous gene pairs. Although the gene balance hypothesis makes predictions about the expression response to homoeologous exchanges, they have not been empirically tested. We used genomic and transcriptomic data from 6 resynthesized, isogenic Brassica napus lines over 10 generations to identify homoeologous exchanges, analyzed expression responses, and tested for patterns of genomic imbalance. Groups of dosage-sensitive genes had less variable expression responses to homoeologous exchanges than dosage-insensitive genes, a sign that their relative dosage is constrained. This difference was absent for homoeologous pairs whose expression was biased toward the B. napus A subgenome. Finally, the expression response to homoeologous exchanges was more variable than the response to whole-genome duplication, suggesting homoeologous exchanges create genomic imbalance. These findings expand our knowledge of the impact of dosage balance selection on genome evolution and potentially connect patterns in polyploid genomes over time, from homoeolog expression bias to duplicate gene retention.

Genetic factors inherited from both diploid parents interact to affect genome stability and fertility in resynthesized allotetraploid Brassica napus

Established allopolyploids are known to be genomically stable and fertile. However, in contrast, most newly resynthesized allopolyploids are infertile and meiotically unstable. Identifying the genetic factors responsible for genome stability in newly formed allopolyploid is key to understanding how 2 genomes come together to form a species. One hypothesis is that established allopolyploids may have inherited specific alleles from their diploid progenitors which conferred meiotic stability. Resynthesized Brassica napus lines are often unstable and infertile, unlike B. napus cultivars. We tested this hypothesis by characterizing 41 resynthesized B. napus lines produced by crosses between 8 Brassica rapa and 8 Brassica oleracea lines for copy number variation resulting from nonhomologous recombination events and fertility. We resequenced 8 B. rapa and 5 B. oleracea parent accessions and analyzed 19 resynthesized lines for allelic variation in a list of meiosis gene homologs. SNP genotyping was performed using the Illumina Infinium Brassica 60K array for 3 individuals per line. Self-pollinated seed set and genome stability (number of copy number variants) were significantly affected by the interaction between both B. rapa and B. oleracea parental genotypes. We identified 13 putative meiosis gene candidates which were significantly associated with frequency of copy number variants and which contained putatively harmful mutations in meiosis gene haplotypes for further investigation. Our results support the hypothesis that allelic variants inherited from parental genotypes affect genome stability and fertility in resynthesized rapeseed.

Patterns, mechanisms, and consequences of homoeologous exchange in allopolyploid angiosperms: a genomic and epigenomic perspective

Allopolyploids result from hybridization between different evolutionary lineages coupled with genome doubling. Homoeologous chromosomes (chromosomes with common shared ancestry) may undergo recombination immediately after allopolyploid formation and continue over successive generations. The outcome of this meiotic pairing behavior is dynamic and complex. Homoeologous exchanges (HEs) may lead to the formation of unbalanced gametes, reduced fertility, and selective disadvantage. By contrast, HEs could act as sources of novel evolutionary substrates, shifting the relative dosage of parental gene copies, generating novel phenotypic diversity, and helping the establishment of neo-allopolyploids. However, HE patterns vary among lineages, across generations, and even within individual genomes and chromosomes. The causes and consequences of this variation are not fully understood, though interest in this evolutionary phenomenon has increased in the last decade. Recent technological advances show promise in uncovering the mechanistic basis of HEs. Here, we describe recent observations of the common patterns among allopolyploid angiosperm lineages, underlying genomic and epigenomic features, and consequences of HEs. We identify critical research gaps and discuss future directions with far-reaching implications in understanding allopolyploid evolution and applying them to the development of important phenotypic traits of polyploid crops.

The final piece of the Triangle of U: Evolution of the tetraploid Brassica carinata genome

Ethiopian mustard (Brassica carinata) is an ancient crop with remarkable stress resilience and a desirable seed fatty acid profile for biofuel uses. Brassica carinata is one of six Brassica species that share three major genomes from three diploid species (AA, BB, and CC) that spontaneously hybridized in a pairwise manner to form three allotetraploid species (AABB, AACC, and BBCC). Of the genomes of these species, that of B. carinata is the least understood. Here, we report a chromosome scale 1.31-Gbp genome assembly with 156.9-fold sequencing coverage for B. carinata, completing the reference genomes comprising the classic Triangle of U, a classical theory of the evolutionary relationships among these six species. Our assembly provides insights into the hybridization event that led to the current B. carinata genome and the genomic features that gave rise to the superior agronomic traits of B. carinata. Notably, we identified an expansion of transcription factor networks and agronomically important gene families. Completion of the Triangle of U comparative genomics platform has allowed us to examine the dynamics of polyploid evolution and the role of subgenome dominance in the domestication and continuing agronomic improvement of B. carinata and other Brassica species.

Reconstruction of ancestral karyotype illuminates chromosome evolution in the genus Cucumis

Karyotype dynamics driven by complex chromosome rearrangements constitute a fundamental issue in evolutionary genetics. The evolutionary events underlying karyotype diversity within plant genera, however, have rarely been reconstructed from a computed ancestral progenitor. Here, we developed a method to rapidly and accurately represent extant karyotypes with the genus, Cucumis, using highly customizable comparative oligo-painting (COP) allowing visualization of fine-scale genome structures of eight Cucumis species from both African-origin and Asian-origin clades. Based on COP data, an evolutionary framework containing a genus-level ancestral karyotype was reconstructed, allowing elucidation of the evolutionary events that account for the origin of these diverse genomes within Cucumis. Our results characterize the cryptic rearrangement hotspots on ancestral chromosomes, and demonstrate that the ancestral Cucumis karyotype (n = 12) evolved to extant Cucumis genomes by hybridizations and frequent lineage- and species-specific genome reshuffling. Relative to the African species, the Asian species, including melon (Cucumis melo, n = 12), Cucumis hystrix (n = 12) and cucumber (Cucumis sativus, n = 7), had highly shuffled genomes caused by large-scale inversions, centromere repositioning and chromothripsis-like rearrangement. The deduced reconstructed ancestral karyotype for the genus allowed us to propose evolutionary trajectories and specific events underlying the origin of these Cucumis species. Our findings highlight that the partitioned evolutionary plasticity of Cucumis karyotype is primarily located in the centromere-proximal regions marked by rearrangement hotspots, which can potentially serve as a reservoir for chromosome evolution due to their fragility.

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.

Chromosome inheritance and meiotic stability in allopolyploid Brassica napus

Homoeologous recombination, aneuploidy, and other genetic changes are common in resynthesized allopolyploid Brassica napus. In contrast, the chromosomes of cultivars have long been considered to be meiotically stable. To gain a better understanding of the underlying mechanisms leading to stabilization in the allopolyploid, the behavior of chromosomes during meiosis can be compared by unambiguous chromosome identification between resynthesized and natural B. napus. Compared with natural B. napus, resynthesized lines show high rates of nonhomologous centromere association, homoeologous recombination leading to translocation, homoeologous chromosome replacement, and association and breakage of 45S rDNA loci. In both natural and resynthesized B. napus, we observed low rates of univalents, A–C bivalents, and early sister chromatid separations. Reciprocal homoeologous chromosome exchanges and double reductions were photographed for the first time in meiotic telophase I. Meiotic errors were non-uniformly distributed across the genome in resynthesized B. napus, and in particular homoeologs sharing synteny along their entire length exhibited multivalents at diakinesis and polysomic inheritance at telophase I. Natural B. napus appeared to resolve meiotic errors mainly by suppressing homoeologous pairing, resolving nonhomologous centromere associations and 45S rDNA associations before diakinesis, and reducing homoeologous cross-overs.

Homoeologous Exchanges, Segmental Allopolyploidy, and Polyploid Genome Evolution

Polyploidy is a major force in plant evolution and speciation. In newly formed allopolyploids, pairing between related chromosomes from different subgenomes (homoeologous chromosomes) during meiosis is common. The initial stages of allopolyploid formation are characterized by a spectrum of saltational genomic and regulatory alterations that are responsible for evolutionary novelty. Here we highlight the possible effects and roles of recombination between homoeologous chromosomes during the early stages of allopolyploid stabilization. Homoeologous exchanges (HEs) have been reported in young allopolyploids from across the angiosperms. Although all lineages undergo karyotype change via chromosome rearrangements over time, the early generations after allopolyploid formation are predicted to show an accelerated rate of genomic change. HEs can also cause changes in allele dosage, genome-wide methylation patterns, and downstream phenotypes, and can hence be responsible for speciation and genome stabilization events. Additionally, we propose that fixation of duplication - deletion events resulting from HEs could lead to the production of genomes which appear to be a mix of autopolyploid and allopolyploid segments, sometimes termed "segmental allopolyploids." We discuss the implications of these findings for our understanding of the relationship between genome instability in novel polyploids and genome evolution.

"Doubled-haploid" allohexaploid Brassica lines lose fertility and viability and accumulate genetic variation due to genomic instability

Microspore culture stimulates immature pollen grains to develop into plants via tissue culture and is used routinely in many crop species to produce "doubled haploids": homozygous, true-breeding lines. However, microspore culture is also often used on material that does not have stable meiosis, such as interspecific hybrids. In this case, the resulting progeny may lose their "doubled haploid" homozygous status as a result of chromosome missegregation and homoeologous exchanges. However, little is known about the frequency of these effects. We assessed fertility, meiosis and genetic variability in self-pollinated progeny sets (the MDL2 population) resulting from first-generation plants (the MDL1 population) derived from microspores of a near-allohexaploid interspecific hybrid from the cross (Brassica napus × B. carinata) × B. juncea. Allelic inheritance and copy number variation were predicted using single nucleotide polymorphism marker data from the Illumina Infinium 60K Brassica array. Seed fertility and viability decreased substantially from the MDL1 to the MDL2 generation. In the MDL2 population, 87% of individuals differed genetically from their MDL1 parent. These genetic differences resulted from novel homoeologous exchanges between chromosomes, chromosome loss and gain, and segregation and instability of pre-existing karyotype abnormalities. Novel karyotype change was extremely common, with 2.2 new variants observed per MDL2 individual. Significant differences between progeny sets in the number of novel genetic variants were also observed. Meiotic instability clearly has the potential to dramatically change karyotypes (often without detectable effects on the presence or absence of alleles) in putatively homozygous, microspore-derived lines, resulting in loss of fertility and viability.

Detecting de Novo Homoeologous Recombination Events in Cultivated Brassica napus Using a Genome-Wide SNP Array

The heavy selection pressure due to intensive breeding of Brassica napus has created a narrow gene pool, limiting the ability to produce improved varieties through crosses between B. napus cultivars. One mechanism that has contributed to the adaptation of important agronomic traits in the allotetraploid B. napus has been chromosomal rearrangements resulting from homoeologous recombination between the constituent A and C diploid genomes. Determining the rate and distribution of such events in natural B. napus will assist efforts to understand and potentially manipulate this phenomenon. The Brassica high-density 60K SNP array, which provides genome-wide coverage for assessment of recombination events, was used to assay 254 individuals derived from 11 diverse cultivated spring type B. napus. These analyses identified reciprocal allele gain and loss between the A and C genomes and allowed visualization of de novo homoeologous recombination events across the B. napus genome. The events ranged from loss/gain of 0.09 Mb to entire chromosomes, with almost 5% aneuploidy observed across all gametes. There was a bias toward sub-telomeric exchanges leading to genome homogenization at chromosome termini. The A genome replaced the C genome in 66% of events, and also featured more dominantly in gain of whole chromosomes. These analyses indicate de novo homoeologous recombination is a continuous source of variation in established Brassica napus and the rate of observed events appears to vary with genetic background. The Brassica 60K SNP array will be a useful tool in further study and manipulation of this phenomenon.

Homoeologous exchange is a major cause of gene presence/absence variation in the amphidiploid Brassica napus

Homoeologous exchanges (HEs) have been shown to generate novel gene combinations and phenotypes in a range of polyploid species. Gene presence/absence variation (PAV) is also a major contributor to genetic diversity. In this study, we show that there is an association between these two events, particularly in recent Brassica napus synthetic accessions, and that these represent a novel source of genetic diversity, which can be captured for the improvement of this important crop species. By assembling the pangenome of B. napus, we show that 38% of the genes display PAV behaviour, with some of these variable genes predicted to be involved in important agronomic traits including flowering time, disease resistance, acyl lipid metabolism and glucosinolate metabolism. This study is a first and provides a detailed characterization of the association between HEs and PAVs in B. napus at the pangenome level.

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.

Homoeologous exchanges cause extensive dosage-dependent gene expression changes in an allopolyploid crop

Structural variation is a major source of genetic diversity and an important substrate for selection. In allopolyploids, homoeologous exchanges (i.e. between the constituent subgenomes) are a very frequent type of structural variant. However, their direct impact on gene content and gene expression had not been determined. Here, we used a tissue-specific mRNA-Seq dataset to measure the consequences of homoeologous exchanges (HE) on gene expression in Brassica napus, a representative allotetraploid crop. We demonstrate that expression changes are proportional to the change in gene copy number triggered by the HEs. Thus, when homoeologous gene pairs have unbalanced transcriptional contributions before the HE, duplication of one copy does not accurately compensate for loss of the other and combined homoeologue expression also changes. These effects are, however, mitigated over time. This study sheds light on the origins, timing and functional consequences of homeologous exchanges in allopolyploids. It demonstrates that the interplay between new structural variation and the resulting impacts on gene expression, influences allopolyploid genome evolution.

The high-quality genome of Brassica napus cultivar 'ZS11' reveals the introgression history in semi-winter morphotype

Allotetraploid oilseed rape (Brassica napus L.) is an agriculturally important crop. Cultivation and breeding of B. napus by humans has resulted in numerous genetically diverse morphotypes with optimized agronomic traits and ecophysiological adaptation. To further understand the genetic basis of diversification and adaptation, we report a draft genome of an Asian semi-winter oilseed rape cultivar 'ZS11' and its comprehensive genomic comparison with the genomes of the winter-type cultivar 'Darmor-bzh' as well as two progenitors. The integrated BAC-to-BAC and whole-genome shotgun sequencing strategies were effective in the assembly of repetitive regions (especially young long terminal repeats) and resulted in a high-quality genome assembly of B. napus 'ZS11'. Within a short evolutionary period (~6700 years ago), semi-winter-type 'ZS11' and the winter-type 'Darmor-bzh' maintained highly genomic collinearity. Even so, certain genetic differences were also detected in two morphotypes. Relative to 'Darmor-bzh', both two subgenomes of 'ZS11' are closely related to its progenitors, and the 'ZS11' genome harbored several specific segmental homoeologous exchanges (HEs). Furthermore, the semi-winter-type 'ZS11' underwent potential genomic introgressions with B. rapa (A_r ). Some of these genetic differences were associated with key agronomic traits. A key gene of A03.FLC3 regulating vernalization-responsive flowering time in 'ZS11' was first experienced HE, and then underwent genomic introgression event with A_r , which potentially has led to genetic differences in controlling vernalization in the semi-winter types. Our observations improved our understanding of the genetic diversity of different B. napus morphotypes and the cultivation history of semi-winter oilseed rape in Asia.

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.

Targeted deep sequencing of flowering regulators in Brassica napus reveals extensive copy number variation

Gene copy number variation (CNV) is increasingly implicated in control of complex trait networks, particularly in polyploid plants like rapeseed (Brassica napus L.) with an evolutionary history of genome restructuring. Here we performed sequence capture to assay nucleotide variation and CNV in a panel of central flowering time regulatory genes across a species-wide diversity set of 280 B. napus accessions. The genes were chosen based on prior knowledge from Arabidopsis thaliana and related Brassica species. Target enrichment was performed using the Agilent SureSelect technology, followed by Illumina sequencing. A bait (probe) pool was developed based on results of a preliminary experiment with representatives from different B. napus morphotypes. A very high mean target coverage of ~670x allowed reliable calling of CNV, single nucleotide polymorphisms (SNPs) and insertion-deletion (InDel) polymorphisms. No accession exhibited no CNV, and at least one homolog of every gene we investigated showed CNV in some accessions. Some CNV appear more often in specific morphotypes, indicating a role in diversification.

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.

Post-polyploidisation morphotype diversification associates with gene copy number variation

Genetic models for polyploid crop adaptation provide important information relevant for future breeding prospects. A well-suited model is Brassica napus, a recent allopolyploid closely related to Arabidopsis thaliana. Flowering time is a major adaptation trait determining life cycle synchronization with the environment. Here we unravel natural genetic variation in B. napus flowering time regulators and investigate associations with evolutionary diversification into different life cycle morphotypes. Deep sequencing of 35 flowering regulators was performed in 280 diverse B. napus genotypes. High sequencing depth enabled high-quality calling of single-nucleotide polymorphisms (SNPs), insertion-deletions (InDels) and copy number variants (CNVs). By combining these data with genotyping data from the Brassica 60 K Illumina® Infinium SNP array, we performed a genome-wide marker distribution analysis across the 4 ecogeographical morphotypes. Twelve haplotypes, including Bna.FLC.A10, Bna.VIN3.A02 and the Bna.FT promoter on C02_random, were diagnostic for the diversification of winter and spring types. The subspecies split between oilseed/kale (B. napus ssp. napus) and swedes/rutabagas (B. napus ssp. napobrassica) was defined by 13 haplotypes, including genomic rearrangements encompassing copies of Bna.FLC, Bna.PHYA and Bna.GA3ox1. De novo variation in copies of important flowering-time genes in B. napus arose during allopolyploidisation, enabling sub-functionalisation that allowed different morphotypes to appropriately fine-tune their lifecycle.

Distinct subgenome stabilities in synthesized Brassica allohexaploids

Trigenomic Brassica allohexaploids synthesized from three crossing strategies showed diploidized and non-diploidized meiotic behaviors and produced both euploid and aneuploid progenies during successive generations, revealing the distinct subgenome stabilities (B > A> C). Three cultivated allotetraploid Brassica species (Brassica napus, B. juncea, B. carinata) represent the model system of speciation through interspecific hybridization and allopolyploidization, but no Brassica species at higher ploidy level exists in nature. In this study, Brassica allohexaploids (2n = 54, AABBCC) were artificially synthesized using three crossing strategies, and had combinations of the genomes from the extant allotetraploids and diploids (B. rapa, B. oleracea and B. nigra). The chromosome numbers and complements of these allohexaploids and the self-pollinated progenies of successive generations (S0-S7) were determined using multicolor fluorescent in situ hybridization that distinguished the chromosomes of three constituent genomes from each other. Both euploid and aneuploid progenies were identified. The most aneuploids maintained all B- and A-genome chromosomes and variable number of C-genome chromosomes, suggesting that genome stability was B > A > C. In the extreme case, loss of whole set of C-genome chromosomes led to the production of B. juncea-type progeny. Some aneuploid progenies had the same number of chromosomes (2n = 54) as the euploid, but the simultaneous loss and gain of A- and C-genome chromosomes. The diploidized and non-diploidized meiotic behaviors co-occurred in all allohexaploid individuals of consecutive generations. The aberrant chromosome pairing and segregation mainly involved the chromosomes of A and C genomes, which resulted in aneuploidy in self-pollinated progenies. The mechanisms for the differential stability of three genomes and the stabilization of the new allohexaploids are discussed.

Genome-wide association analyses reveal complex genetic architecture underlying natural variation for flowering time in canola

Optimum flowering time is the key to maximize canola production in order to meet global demand of vegetable oil, biodiesel and canola-meal. We reveal extensive variation in flowering time across diverse genotypes of canola under field, glasshouse and controlled environmental conditions. We conduct a genome-wide association study and identify 69 single nucleotide polymorphism (SNP) markers associated with flowering time, which are repeatedly detected across experiments. Several associated SNPs occur in clusters across the canola genome; seven of them were detected within 20 Kb regions of a priori candidate genes; FLOWERING LOCUS T, FRUITFUL, FLOWERING LOCUS C, CONSTANS, FRIGIDA, PHYTOCHROME B and an additional five SNPs were localized within 14 Kb of a previously identified quantitative trait loci for flowering time. Expression analyses showed that among FLC paralogs, BnFLC.A2 accounts for ~23% of natural variation in diverse accessions. Genome-wide association analysis for FLC expression levels mapped not only BnFLC.C2 but also other loci that contribute to variation in FLC expression. In addition to revealing the complex genetic architecture of flowering time variation, we demonstrate that the identified SNPs can be modelled to predict flowering time in diverse canola germplasm accurately and hence are suitable for genomic selection of adaptative traits in canola improvement programmes.

Interspecific hybridization, polyploidization, and backcross of Brassica oleracea var. alboglabra with B. rapa var. purpurea morphologically recapitulate the evolution of Brassica vegetables

Brassica oleracea and B. rapa are two important vegetable crops. Both are composed of dozens of subspecies encompassing hundreds of varieties and cultivars. Synthetic B. napus with these two plants has been used extensively as a research model for the investigation of allopolyploid evolution. However, the mechanism underlying the explosive evolution of hundreds of varieties of B. oleracea and B. rapa within a short period is poorly understood. In the present study, interspecific hybridization between B. oleracea var. alboglabra and B. rapa var. purpurea was performed. The backcross progeny displayed extensive morphological variation, including some individuals that phenocopied subspecies other than their progenitors. Numerous interesting novel phenotypes and mutants were identified among the backcross progeny. The chromosomal recombination between the A and C genomes and the chromosomal asymmetric segregation were revealed using Simple Sequence Repeats (SSR) markers. These findings provide direct evidence in support of the hypothesis that interspecific hybridization and backcrossing have played roles in the evolution of the vast variety of vegetables among these species and suggest that combination of interspecific hybridization and backcrossing may facilitate the development of new mutants and novel phenotypes for both basic research and the breeding of new vegetable crops.

Effects of genome structure variation, homeologous genes and repetitive DNA on polyploid crop research in the age of genomics

Compared to diploid species, allopolyploid crop species possess more complex genomes, higher productivity, and greater adaptability to changing environments. Next generation sequencing techniques have produced high-density genetic maps, whole genome sequences, transcriptomes and epigenomes for important polyploid crops. However, several problems interfere with the full application of next generation sequencing techniques to these crops. Firstly, different types of genomic variation affect sequence assembly and QTL mapping. Secondly, duplicated or homoeologous genes can diverge in function and then lead to emergence of many minor QTL, which increases difficulties in fine mapping, cloning and marker assisted selection. Thirdly, repetitive DNA sequences arising in polyploid crop genomes also impact sequence assembly, and are increasingly being shown to produce small RNAs to regulate gene expression and hence phenotypic traits. We propose that these three key features should be considered together when analyzing polyploid crop genomes. It is apparent that dissection of genomic structural variation, elucidation of the function and mechanism of interaction of homoeologous genes, and investigation of the de novo roles of repeat sequences in agronomic traits are necessary for genomics-based crop breeding in polyploids.

Microspore culture reveals complex meiotic behaviour in a trigenomic Brassica hybrid

BACKGROUND

Development of synthetic allohexaploid Brassica (2n = AABBCC) would be beneficial for agriculture, as allelic contributions from three genomes could increase hybrid vigour and broaden adaptation. Microspore culture of a near-allohexaploid hybrid derived from the cross (B. napus × B. carinata) × B. juncea was undertaken in order to assess the frequency and distribution of homologous and homoeologous crossovers in this trigenomic hybrid. SNP and SSR molecular markers were used to detect inheritance of A, B and C genome alleles in microspore-derived (MD) progeny. SNP allele copy number was also assessed. The MD progeny were also compared to progeny derived by self-pollination and open-pollination for fertility (estimated by self-pollinated seed set and pollen viability) and DNA ploidy (measured by flow cytometry).

RESULTS

In the trigenomic hybrid, homologous chromosome pairs A(j)-A(n), B(j)-B(c) and C(n)-C(c) had similar meiotic crossover frequencies and segregation to that previously observed in established Brassica species, as demonstrated by marker haplotype analysis of the MD population. Homoeologous pairing between chromosomes A1-C1, A2-C2 and A7-C6 was detected at frequencies of 12-18 %, with other homoeologous chromosome regions associating from 8 % (A3-C3) to 0-1 % (A8-C8, A8-C9) of the time. Copy number analysis revealed eight instances of additional chromosomes and 20 instances of chromosomes present in one copy in somatically doubled MD progeny. Presence of chromosome A6 was positively correlated with self-pollinated seed set and pollen viability in the MD population. Many MD progeny were unable to produce self-pollinated seed (76 %) or viable pollen (53 %), although one MD plant produced 198 self-pollinated seeds. Average fertility was significantly lower in progeny obtained by microspore culture than progeny obtained by self-pollination or open-pollination, after excluding MD progeny which had not undergone chromosome doubling.

CONCLUSIONS

Based on SNP data analysis of the microspore-derived progeny, crossover frequency per chromosome in the allohexaploid hybrid was found to be similar to that in established Brassica species, suggesting that the higher chromosome number did not significantly disrupt cellular regulation of meiosis. SNP allele copy number analysis revealed the occurrence not only of homoeologous duplication/deletion events but also other cryptic duplications and deletions that may have been the result of mitotic instability. Microspore culture simplified the assessment of chromosome behaviour in the allohexaploid hybrid but yielded progeny with lower fertility and a greater range of ploidy levels compared to progeny obtained by self- or open-pollination.

Genome-specific differential gene expressions in resynthesized Brassica allotetraploids from pair-wise crosses of three cultivated diploids revealed by RNA-seq

Polyploidy is popular for the speciation of angiosperms but the initial stage of allopolyploidization resulting from interspecific hybridization and genome duplication is associated with different extents of changes in genome structure and gene expressions. Herein, the transcriptomes detected by RNA-seq in resynthesized Brassica allotetraploids (Brassica juncea, AABB; B. napus, AACC; B. carinata, BBCC) from the pair-wise crosses of the same three diploids (B. rapa, AA; B. nigra, BB; B. oleracea, CC) were compared to reveal the patterns of gene expressions from progenitor genomes and the effects of different types of genome combinations and cytoplasm, upon the genome merger and duplication. From transcriptomic analyses for leaves and silique walls, extensive expression alterations were revealed in these resynthesized allotetraploids relative to their diploid progenitors, as well as during the transition from vegetative to reproductive development, for differential and transgressive gene expressions were variable in numbers and functions. Genes involved in glucosinolates and DNA methylation were transgressively up-regulated among most samples, suggesting that gene expression regulation was immediately established after allopolyploidization. The expression of ribosomal protein genes was also tissue-specific and showed a similar expression hierarchy of rRNA genes. The balance between the co-up and co-down regulation was observed between reciprocal B. napus with different types of the cytoplasm. Our results suggested that gene expression changes occurred after initial genome merger and such profound alterations might enhance the growth vigor and adaptability of Brassica allotetraploids.

SNP markers-based map construction and genome-wide linkage analysis in Brassica napus

An Illumina Infinium array comprising 5306 single nucleotide polymorphism (SNP) markers was used to genotype 175 individuals of a doubled haploid population derived from a cross between Skipton and Ag-Spectrum, two Australian cultivars of rapeseed (Brassica napus L.). A genetic linkage map based on 613 SNP and 228 non-SNP (DArT, SSR, SRAP and candidate gene markers) covering 2514.8 cM was constructed and further utilized to identify loci associated with flowering time and resistance to blackleg, a disease caused by the fungus Leptosphaeria maculans. Comparison between genetic map positions of SNP markers and the sequenced Brassica rapa (A) and Brassica oleracea (C) genome scaffolds showed several genomic rearrangements in the B. napus genome. A major locus controlling resistance to L. maculans was identified at both seedling and adult plant stages on chromosome A07. QTL analyses revealed that up to 40.2% of genetic variation for flowering time was accounted for by loci having quantitative effects. Comparative mapping showed Arabidopsis and Brassica flowering genes such as Phytochrome A/D, Flowering Locus C and agamous-Like MADS box gene AGL1 map within marker intervals associated with flowering time in a DH population from Skipton/Ag-Spectrum. Genomic regions associated with flowering time and resistance to L. maculans had several SNP markers mapped within 10 cM. Our results suggest that SNP markers will be suitable for various applications such as trait introgression, comparative mapping and high-resolution mapping of loci in B. napus.

High-resolution molecular karyotyping uncovers pairing between ancestrally related Brassica chromosomes

How do chromosomal regions with differing degrees of homology and homeology interact at meiosis? We provide a novel analytical method based on simple genetics principles which can help to answer this important question. This method interrogates high-throughput molecular marker data in order to infer chromosome behavior at meiosis in interspecific hybrids. We validated this method using high-resolution molecular marker karyotyping in two experimental Brassica populations derived from interspecific crosses among B. juncea, B. napus and B. carinata, using a single nucleotide polymorphism chip. This method of analysis successfully identified meiotic interactions between chromosomes sharing different degrees of similarity: full-length homologs; full-length homeologs; large sections of primary homeologs; and small sections of secondary homeologs. This analytical method can be applied to any allopolyploid species or fertile interspecific hybrid in order to detect meiotic associations. This genetic information can then be used to identify which genomic regions share functional homeology (i.e., retain enough similarity to allow pairing and segregation at meiosis). When applied to interspecific hybrids for which reference genome sequences are available, the question of how differing degrees of homology and homeology affect meiotic interactions may finally be resolved.

The fate of chromosomes and alleles in an allohexaploid Brassica population

Production of allohexaploid Brassica (2n = AABBCC) is a promising goal for plant breeders due to the potential for hybrid heterosis and useful allelic contributions from all three of the Brassica genomes present in the cultivated diploid (2n = AA, 2n = BB, 2n = CC) and allotetraploid (2n = AABB, 2n = AACC, and 2n = BBCC) crop species (canola, cabbages, mustards). We used high-throughput SNP molecular marker assays, flow cytometry, and fluorescent in situ hybridization (FISH) to characterize a population of putative allohexaploids derived from self-pollination of a hybrid from the novel cross (B. napus × B. carinata) × B. juncea to investigate whether fertile, stable allohexaploid Brassica can be produced. Allelic segregation in the A and C genomes generally followed Mendelian expectations for an F2 population, with minimal nonhomologous chromosome pairing. However, we detected no strong selection for complete 2n = AABBCC chromosome complements, with weak correlations between DNA content and fertility (r(2) = 0.11) and no correlation between missing chromosomes or chromosome segments and fertility. Investigation of next-generation progeny resulting from one highly fertile F2 plant using FISH revealed general maintenance of high chromosome numbers but severe distortions in karyotype, as evidenced by recombinant chromosomes and putative loss/duplication of A- and C-genome chromosome pairs. Our results show promise for the development of meiotically stable allohexaploid lines, but highlight the necessity of selection for 2n = AABBCC karyotypes.

A consensus map of rapeseed (Brassica napus L.) based on diversity array technology markers: applications in genetic dissection of qualitative and quantitative traits

Background

Dense consensus genetic maps based on high-throughput genotyping platforms are valuable for making genetic gains in Brassica napus through quantitative trait locus identification, efficient predictive molecular breeding, and map-based gene cloning. This report describes the construction of the first B. napus consensus map consisting of a 1,359 anchored array based genotyping platform; Diversity Arrays Technology (DArT), and non-DArT markers from six populations originating from Australia, Canada, China and Europe. We aligned the B. napus DArT sequences with genomic scaffolds from Brassica rapa and Brassica oleracea, and identified DArT loci that showed linkage with qualitative and quantitative loci associated with agronomic traits.

Results

The integrated consensus map covered a total of 1,987.2 cM and represented all 19 chromosomes of the A and C genomes, with an average map density of one marker per 1.46 cM, corresponding to approximately 0.88 Mbp of the haploid genome. Through in silico physical mapping 2,457 out of 3,072 (80%) DArT clones were assigned to the genomic scaffolds of B. rapa (A genome) and B. oleracea (C genome). These were used to orientate the genetic consensus map with the chromosomal sequences. The DArT markers showed linkage with previously identified non-DArT markers associated with qualitative and quantitative trait loci for plant architecture, phenological components, seed and oil quality attributes, boron efficiency, sucrose transport, male sterility, and race-specific resistance to blackleg disease.

Conclusions

The DArT markers provide increased marker density across the B. napus genome. Most of the DArT markers represented on the current array were sequenced and aligned with the B. rapa and B. oleracea genomes, providing insight into the Brassica A and C genomes. This information can be utilised for comparative genomics and genomic evolution studies. In summary, this consensus map can be used to (i) integrate new generation markers such as SNP arrays and next generation sequencing data; (ii) anchor physical maps to facilitate assembly of B. napus genome sequences; and (iii) identify candidate genes underlying natural genetic variation for traits of interest.

Persistent whole-chromosome aneuploidy is generally associated with nascent allohexaploid wheat

Allopolyploidization has been a driving force in plant evolution. Formation of common wheat (Triticum aestivum L.) represents a classic example of successful speciation via allopolyploidy. Nevertheless, the immediate chromosomal consequences of allopolyploidization in wheat remain largely unexplored. We report here an in-depth investigation on transgenerational chromosomal variation in resynthesized allohexaploid wheats that are identical in genome constitution to common wheat. We deployed sequential FISH, genomic in situ hybridization (GISH), and homeolog-specific pyrosequencing, which enabled unequivocal identification of each of the 21 homologous chromosome pairs in each of >1,000 individual plants from 16 independent lines. We report that whole-chromosome aneuploidy occurred ubiquitously in early generations (from selfed generation S(1) to >S(20)) of wheat allohexaploidy although at highly variable frequencies (20-100%). In contrast, other types of gross structural variations were scant. Aneuploidy included an unexpected hidden type, which had a euploid chromosome number of 2n = 42 but with simultaneous loss and gain of nonhomeologous chromosomes. Of the three constituent subgenomes, B showed the most lability for aneuploidy, followed by A, but the recently added D subgenome was largely stable in most of the studied lines. Chromosome loss and gain were also unequal across the 21 homologous chromosome pairs. Pedigree analysis showed no evidence for progressive karyotype stabilization even with multigenerational selection for euploidy. Profiling of two traits directly related to reproductive fitness showed that although pollen viability was generally reduced by aneuploidy, the adverse effect of aneuploidy on seed-set is dependent on both aneuploidy type and synthetic line.