Polyploid formation pathways have an impact on genetic rearrangements in resynthesized Brassica napus

Unambiguous chromosome identification reveals the factors impacting irregular chromosome behaviors in allotriploid AAC Brassica

KEY MESSAGE

The major irregular chromosome pairing and mis-segregation were detected during meiosis through unambiguous chromosome identification and found that allotriploid Brassica can undergo meiosis successfully and produce mostly viable aneuploid gametes. Triploids have played a crucial role in the evolution of species by forming polyploids and facilitating interploidy gene transfer. It is widely accepted that triploids cannot undergo meiosis normally and predominantly produce nonfunctional aneuploid gametes, which restricts their role in species evolution. In this study, we demonstrated that natural and synthetic allotriploid Brassica (AAC), produced by crossing natural and synthetic Brassica napus (AACC) with Brassica rapa (AA), exhibits basically normal chromosome pairing and segregation during meiosis. Homologous A chromosomes paired faithfully and generally segregated equally. Monosomic C chromosomes were largely retained as univalents and randomly entered daughter cells. The primary irregular meiotic behaviors included associations of homoeologs and 45S rDNA loci at diakinesis, as well as homoeologous chromosome replacement and premature sister chromatid separation at anaphase I. Preexisting homoeologous arrangements altered meiotic behaviors in both chromosome irregular pairing and mis-segregation by increasing the formation of A-genomic univalents and A-C bivalents, as well as premature sister chromatid separation and homologous chromosome nondisjunction. Meiotic behaviors depended significantly on the genetic background and heterozygous homoeologous rearrangement. AAC triploids mainly generated aneuploid gametes, most of which were viable. These results demonstrate that allotriploid Brassica containing an intact karyotype can proceed through meiosis successfully, broadening our current understanding of the inheritance and role in species evolution of allotriploid.

Admixture of divergent genomes facilitates hybridization across species in the family Brassicaceae

Hybridization and polyploidization are pivotal to plant evolution. Genetic crosses between distantly related species are rare in nature due to reproductive barriers but how such hurdles can be overcome is largely unknown. Here we report the hybrid genome structure of xBrassicoraphanus, a synthetic allotetraploid of Brassica rapa and Raphanus sativus. We performed cytogenetic analysis and de novo genome assembly to examine chromosome behaviors and genome integrity in the hybrid. Transcriptome analysis was conducted to investigate expression of duplicated genes in conjunction with epigenome analysis to address whether genome admixture entails epigenetic reconfiguration. Allotetraploid xBrassicoraphanus retains both parental chromosomes without genome rearrangement. Meiotic synapsis formation and chromosome exchange are avoided between nonhomologous progenitor chromosomes. Reconfiguration of transcription network occurs, and less divergent cis-elements of duplicated genes are associated with convergent expression. Genome-wide DNA methylation asymmetry between progenitors is largely maintained but, notably, B. rapa-originated transposable elements are transcriptionally silenced in xBrassicoraphanus through gain of DNA methylation. Our results demonstrate that hybrid genome stabilization and transcription compatibility necessitate epigenome landscape adjustment and rewiring of cis-trans interactions. Overall, this study suggests that a certain extent of genome divergence facilitates hybridization across species, which may explain the great diversification and expansion of angiosperms during evolution.

Genetic Variability in Balkan Paleoendemic Resurrection Plants Ramonda serbica and R. nathaliae Across Their Range and in the Zone of Sympatry

The genus Ramonda includes three Paleoendemic and Tertiary relict species that survived in refugial habitats of the Balkan Peninsula (R. nathaliae and R. serbica) and the Iberian Peninsula (R. myconi). They are all "resurrection plants," a rare phenomenon among flowering plants in Europe. Ramonda myconi and R. nathaliae are diploids (2n = 2x = 48), while R. serbica is a hexaploid (2n = 6x = 144). The two Balkan species occur in sympatry in only two localities in eastern Serbia, where tetraploid potential hybrids (2n = 4x = 96) were found. This observation raised questions about the existence of gene flow between the two species and, more generally, about the evolutionary processes shaping their genetic diversity. To address this question, genetic markers (AFLP) and an estimate of genome size variation were used in a much larger sample and at a larger geographic scale than previously. The combination of AFLP markers and genome size results suggested ongoing processes of interspecific and interploidy hybridization in the two sites of sympatry. The data also showed that interspecific gene flow was strictly confined to sympatry. Elsewhere, both Ramonda species were characterized by low genetic diversity within populations and high population differentiation. This is consistent with the fact that the two species are highly fragmented into small and isolated populations, likely a consequence of their postglacial history. Within sympatry, enormous variability in cytotypes was observed, exceeding most reported cases of mixed ploidy in complex plant species (from 2x to >8x). The AFLP profiles of non-canonical ploidy levels indicated a diversity of origin pathways and that backcrosses probably occur between tetraploid interspecific hybrids and parental species. The question arises whether this diversity of cytotypes corresponds to a transient situation. If not, the question arises as to the genetic and ecological mechanisms that allow this diversity to be maintained over time.

Rapid and Synchronous Breeding of Cytoplasmic Male Sterile and Maintainer Line Through Mitochondrial DNA Rearrangement Using Doubled Haploid Inducer in Brassica napus

When homozygously fertile plants were induced using doubled haploid (DH) induction lines Y3380 and Y3560, the morphology of the induced F1 generation was basically consistent with the female parent, but the fertility was separated, showing characteristics similar to cytoplasmic male sterile (CMS) and maintainer lines. In this study, the morphology, fertility, ploidy, and cytoplasm genotype of the induced progeny were identified, and the results showed that the sterile progeny was polima cytoplasm sterile (pol CMS) and the fertile progeny was nap cytoplasm. The molecular marker and test-cross experimental results showed that the fertile progeny did not carry the restorer gene of pol CMS and the genetic distance between the female parent and the offspring was 0.002. This suggested that those inductions which produced sterile and fertile progeny were coordinated to CMS and maintainer lines. Through the co-linearity analysis of the mitochondrial DNA (mtDNA), it was found that the rearrangement of mtDNA by DH induction was the key factor that caused the transformation of fertility (nap) into sterility (pol). Also, when heterozygous females were induced with DH induction lines, the induction F2 generation also showed the segregation of fertile and sterile lines, and the genetic distance between sterile and fertile lines was approximately 0.075. Therefore, the induction line can induce different types of female parents, and the breeding of the sterile line and the maintainer line can be achieved through the rapid synchronization of sister crosses and self-crosses. The induction of DH inducer in B. napus can provide a new model for the innovation of germplasm resources and open up a new way for its application.

Novel Brassica hybrids with different resistance to Leptosphaeria maculans reveal unbalanced rDNA signal patterns

Hybridization of Brassica napus with various Brassicaceae species can result in obtaining new forms with increased resistance to blackleg, a dangerous disease caused mainly by Leptosphaeria maculans. In this study, we aimed to correlate the field resistance of selected Brassica hybrids to blackleg with chromosomal structure revealed by Fluorescence in situ hybridization. Tested genotypes varied in the number of chromosomes and rDNA signals. The greatest variation was observed for A1-type chromosomes. Field evaluation also revealed significant differences in L. maculans resistance. Performed analyses allowed to distinguish three B. napus × Brassica fruticulosa genotypes in which variable patterns of chromosomal structure might be connected to field resistance. However, a more thorough study, including the detection of all A-genome chromosomes, is required.

Androgenesis-Based Doubled Haploidy: Past, Present, and Future Perspectives

Androgenesis, which entails cell fate redirection within the microgametophyte, is employed widely for genetic gain in plant breeding programs. Moreover, androgenesis-responsive species provide tractable systems for studying cell cycle regulation, meiotic recombination, and apozygotic embryogenesis within plant cells. Past research on androgenesis has focused on protocol development with emphasis on temperature pretreatments of donor plants or floral buds, and tissue culture optimization because androgenesis has different nutritional requirements than somatic embryogenesis. Protocol development for new species and genotypes within responsive species continues to the present day, but slowly. There is more focus presently on understanding how protocols work in order to extend them to additional genotypes and species. Transcriptomic and epigenetic analyses of induced microspores have revealed some of the cellular and molecular responses required for or associated with androgenesis. For example, microRNAs appear to regulate early microspore responses to external stimuli; trichostatin-A, a histone deacetylase inhibitor, acts as an epigenetic additive; ά-phytosulfokine, a five amino acid sulfated peptide, promotes androgenesis in some species. Additionally, present work on gene transfer and genome editing in microspores suggest that future endeavors will likely incorporate greater precision with the genetic composition of microspores used in doubled haploid breeding, thus likely to realize a greater impact on crop improvement. In this review, we evaluate basic breeding applications of androgenesis, explore the utility of genomics and gene editing technologies for protocol development, and provide considerations to overcome genotype specificity and morphogenic recalcitrance in non-model plant systems.

Genome-wide unbalanced expression bias and expression level dominance toward Brassica oleracea in artificially synthesized intergeneric hybrids of Raphanobrassica

Raphanobrassica (RrRrCrCr, 2n = 4x = 36), which is generated by distant hybridization between the maternal parent Raphanus sativus (RsRs, 2n = 2x = 18) and the paternal parent Brassica oleracea (C°C°, 2n = 2x = 18), displays intermediate silique phenotypes compared to diploid progenitors. However, the hybrid shares much more similarities in silique phenotypes with those of B. oleracea than those of R. sativus. Strikingly, the silique of Raphanobrassica is obviously split into two parts. To investigate the gene expression patterns behind these phenomena, transcriptome analysis was performed on the upper, middle, and lower sections of pods (RCsiu, RCsim, and RCsil), seeds in the upper and lower sections of siliques (RCseu and RCsel) from Raphanobrassica, whole pods (Rsi and Csi) and all seeds in the siliques (Rse and Cse) from R. sativus and B. oleracea. Transcriptome shock was observed in all five aforementioned tissues of Raphanobrassica. Genome-wide unbalanced biased expression and expression level dominance were also discovered, and both of them were toward B. oleracea in Raphanobrassica, which is consistent with the observed phenotypes. The present results reveal the global gene expression patterns of different sections of siliques of Raphanobrassica, pods, and seeds of B. oleracea and R. sativus, unraveling the tight correlation between global gene expression patterns and phenotypes of the hybrid and its parents.

Genetic characteristics and ploidy trigger the high inducibility of double haploid (DH) inducer in Brassica napus

BACKGROUND

Our recently reported doubled haploid (DH) induction lines e.g., Y3380 and Y3560 are allo-octoploid (AAAACCCC, 2n = 8× ≈ 76), which can induce the maternal parent to produce DH individuals. Whether this induction process is related to the production of aneuploid gametes form male parent and genetic characteristics of the male parent has not been reported yet.

RESULTS

Somatic chromosome counts of DH inducer parents, female wax-less parent (W1A) and their F₁ hybrid individuals revealed the reliability of flow cytometry analysis. Y3560 has normal chromosome behavior in metaphase I and anaphase I, but chromosome division was not synchronized in the tetrad period. Individual phenotypic identification and flow cytometric fluorescence measurement of F₁ individual and parents revealed that DH individuals can be distinguished on the basis of waxiness trait. The results of phenotypic identification and flow cytometry can identify the homozygotes or heterozygotes of F₁ generation individuals. The data of SNP genotyping coupled with phenotypic waxiness trait revealed that the genetic distance between W1A and F₁ homozygotes were smaller as compared to their heterozygotes. It was found that compared with allo-octoploids, aneuploidy from allo-octoploid segregation did not significantly increase the DH induction rate, but reduced male infiltration rate and heterozygous site rate of induced F₁ generation. The ploidy, SNP genotyping and flow cytometry results cumulatively shows that DH induction is attributed to the key genes regulation from the parents of Y3560 and Y3380, which significantly increase the induction efficiency as compared to ploidy.

CONCLUSION

Based on our findings, we hypothesize that genetic characteristics and aneuploidy play an important role in the induction of DH individuals in Brassca napus, and the induction process has been explored. It provides an important insight for us to locate and clone the genes that regulate the inducibility in the later stage.

Challenges and prospects for a potential allohexaploid Brassica crop

The production of a new allohexaploid Brassica crop (2n = AABBCC) is increasingly attracting international interest: a new allohexaploid crop could benefit from several major advantages over the existing Brassica diploid and allotetraploid species, combining genetic diversity and traits from all six crop species with additional allelic heterosis from the extra genome. Although early attempts to produce allohexaploids showed mixed results, recent technological and conceptual advances have provided promising leads to follow. However, there are still major challenges which exist before this new crop type can be realized: (1) incorporation of sufficient genetic diversity to form a basis for breeding and improvement of this potential crop species; (2) restoration of regular meiosis, as most allohexaploids are genetically unstable after formation; and (3) improvement of agronomic traits to the level of "elite" breeding material in the diploid and allotetraploid crop species. In this review, we outline these major prospects and challenges and propose possible plans to produce a stable, diverse and agronomically viable allohexaploid Brassica crop.

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.

The tight genome size of ants: diversity and evolution under ancestral state reconstruction and base composition

The mechanisms and processes driving change and variation in the genome size (GS) are not well known, and only a small set of ant species has been studied. Ants are an ecologically successful insect group present in most distinct ecosystems worldwide. Considering their wide distribution and ecological plasticity in different environmental contexts, we aimed to expand GS estimation within Formicidae to examine distribution patterns and variation in GS and base composition and to reconstruct the ancestral state of this character in an attempt to elucidate the generalized pattern of genomic expansions. Genome size estimates were generated for 99 ant species, including new GS estimates for 91 species of ants, and the mean GS of Formicidae was found to be 0.38 pg. The AT/GC ratio was 62.40/37.60. The phylogenetic reconstruction suggested an ancestral GS of 0.38 pg according to the Bayesian inference/Markov chain Monte Carlo method and 0.37 pg according to maximum likelihood and parsimony methods; significant differences in GS were observed between the subfamilies sampled. Our results suggest that the evolution of GS in Formicidae occurred through loss and accumulation of non-coding regions, mainly transposable elements, and occasionally by whole genome duplication. However, further studies are needed to verify whether these changes in DNA content are related to colonization processes, as suggested at the intraspecific level.

QTL Mapping for Resistance to Cankers Induced by Pseudomonas syringae pv. actinidiae (Psa) in a Tetraploid Actinidia chinensis Kiwifruit Population

Polyploidy is a key driver of significant evolutionary changes in plant species. The genus Actinidia (kiwifruit) exhibits multiple ploidy levels, which contribute to novel fruit traits, high yields and resistance to the canker-causing dieback disease incited by Pseudomonas syringae pv. actinidiae (Psa) biovar 3. However, the genetic mechanism for resistance to Psa observed in polyploid kiwifruit is not yet known. In this study we performed detailed genetic analysis of a tetraploid Actinidia chinensis var. chinensis population derived from a cross between a female parent that exhibits weak tolerance to Psa and a highly Psa-resistant male parent. We used the capture-sequencing approach across the whole kiwifruit genome and generated the first ultra-dense maps in a tetraploid kiwifruit population. We located quantitative trait loci (QTLs) for Psa resistance on these maps. Our approach to QTL mapping is based on the use of identity-by-descent trait mapping, which allowed us to relate the contribution of specific alleles from their respective homologues in the male and female parent, to the control of Psa resistance in the progeny. We identified genes in the diploid reference genome whose function is suggested to be involved in plant defense, which underly the QTLs, including receptor-like kinases. Our study is the first to cast light on the genetics of a polyploid kiwifruit and suggest a plausible mechanism for Psa resistance in this species.

Genome-Wide Duplication of Allotetraploid Brassica napus Produces Novel Characteristics and Extensive Ploidy Variation in Self-Pollinated Progeny

Whole genome duplications (WGDs) have played a major role in angiosperm species evolution. Polyploid plants have undergone multiple cycles of ancient WGD events during their evolutionary history. However, little attention has been paid to the additional WGD of the existing allopolyploids. In this study, we explored the influences of additional WGD on the allopolyploid Brassica napus Compared to tetraploid B. napus, octoploid B. napus (AAAACCCC, 2n = 8x =76) showed significant differences in phenotype, reproductive ability and the ploidy of self-pollinated progeny. Genome duplication also altered a key reproductive organ feature in B. napus, that is, increased the number of pollen apertures. Unlike autopolyploids produced from the diploid Brassica species, the octoploid B. napus produced from allotetraploid B. napus had a relatively stable meiotic process, high pollen viability and moderate fertility under self-pollination conditions, indicating that sub-genomic interactions may be important for the successful establishment of higher-order polyploids. Doubling the genome of B. napus provided us with an opportunity to gain insight into the flexibility of the Brassica genomes. The genome size of self-pollinated progeny of octoploid B. napus varied greatly, and was accompanied by extensive genomic instability, such as aneuploidy, mixed-ploidy and mitotic abnormality. The octoploid B. napus could go through any of genome reduction, equilibrium or expansion in the short-term, thus providing a novel karyotype library for the Brassica genus. Our results reveal the short-term evolutionary consequences of recurrent polyploidization events, and help to deepen our understanding of polyploid plant evolution.

Development and Characteristics of Interspecific Hybrids between Brassica oleracea L. and B. napus L

Interspecific hybridization between B. oleracea inbred lines of head cabbage, Brussels sprouts, kale and B. taurica and inbred lines of rapeseed (B. napus L.) were performed aiming at the development of the new sources of genetic variability of vegetable Brassicas. Using conventional crossings and the embryo-rescue techniques the following interspecific hybrids were developed: 11 genotypes of F1 generation, 18 genotypes of F2 and F1 × F2 generations (produced after self- and cross-pollination of interspecific F1 hybrids), 10 plants of the BC1 generation (resulted from crossing head cabbage cytoplasmic male-sterile lines with interspecific hybrids of the F2 and F1 generations) and 8 plants of BC1 × (F1 × F2). No viable seeds of the BC2 generation (B. oleracea) were obtained due to the strong incompatibility and high mortality of embryos. The morphological characteristics during the vegetative and generative stages, pollen characteristics, seed development and propagation, nuclear DNA contents and genome compositions of interspecific hybrids were analyzed. All the interspecific F1 hybrids were male-fertile with a majority of undeveloped and malformed pollen grains. They showed intermediate values for morphological traits and nuclear DNA contents and had nearly triploid chromosomal numbers (27 to 29) compared with parental lines. The F2 generation had a doubled nuclear DNA content, with 52 and 56 chromosomes, indicating their allohexaploid nature. F2 hybrids were characterized by a high heterosis of morphological characteristics, viable pollen and good seed development. F1 × F2 hybrids were male-fertile with a diversified DNA content and intermediate pollen viability. BC1 plants were male-sterile with an intermediate nuclear DNA content between the F2 and head cabbage, having 28 to 38 chromosomes. Plants of the BC1 × (F1 × F2) generation were in majority male-fertile with 38–46 chromosomes, high seed set, high heterosis and intermediate values for morphological traits. The obtained interspecific hybrids are valuable as new germplasm for improving Brassica-breeding programs.

Response of Brassica napus to Plasmodiophora brassicae Involves Salicylic Acid-Mediated Immunity: An RNA-Seq-Based Study

Clubroot, caused by the obligate parasite Plasmodiophora brassicae, is an important disease of the Brassicaceae and poses a significant threat to the $26.7 billion canola/oilseed rape (Brassica napus) industry in western Canada. While clubroot is managed most effectively by planting resistant host varieties, new pathotypes of P. brassicae have emerged recently that can overcome this resistance. Whole genome analyses provide both a toolbox and a systemic view of molecular mechanisms in host-pathogen interactions, which can be used to design new breeding strategies to increase P. brassicae resistance. We used RNA-seq to evaluate differential gene expression at 7, 14 and 21 days after inoculation (dai) of two B. napus genotypes with differential responses to P. brassicae pathotype 5X. Gall development was evident at 14 dai in the susceptible genotype (the oilseed rape 'Brutor'), while gall development in the resistant genotype (the rutabaga (B. napus) 'Laurentian') was limited and not visible until 21 dai. Immune responses were better sustained through the time-course in 'Laurentian', and numerous genes from immune-related functional categories were associated with salicylic acid (SA)-mediated responses. Jasmonic acid (JA)-mediated responses seemed to be mostly inhibited, especially in the resistant genotype. The upregulation of standard defense-related proteins, like chitinases and thaumatins, was evident in 'Laurentian'. The enrichment, in both host genotypes, of functional categories for syncytium formation and response to nematodes indicated that cell enlargement during P. brassicae infection, and the metabolic processes therein, share similarities with the response to infection by nematodes that produce similar anatomical symptoms. An analysis of shared genes between the two genotypes at different time-points, confirmed that the nematode-like responses occurred earlier for 'Brutor', along with cell metabolism and growth changes. Additionally, the susceptible cultivar turned off defense mechanisms earlier than 'Laurentian'. Collectively, this study showed the importance of SA in triggering immune responses and suggested some key resistance and susceptibility factors that can be used in future studies for resistance breeding through gene-editing approaches.

Meiotic Chromosome Stability and Suppression of Crossover Between Non-homologous Chromosomes in xBrassicoraphanus, an Intergeneric Allotetraploid Derived From a Cross Between Brassica rapa and Raphanus sativus

Hybridization and polyploidization are major driving forces in plant evolution. Allopolyploids can be occasionally formed from a cross between distantly related species but often suffer from chromosome instability and infertility. xBrassicoraphanus is an intergeneric allotetraploid (AARR; 2n = 38) derived from a cross between Brassica rapa (AA; 2n = 20) and Raphanus sativus (RR; 2n = 18). xBrassicoraphanus is fertile and genetically stable, while retaining complete sets of both B. rapa and R. sativus chromosomes. Precise control of meiotic recombination is essential for the production of balanced gametes, and crossovers (COs) must occur exclusively between homologous chromosomes. Many interspecific hybrids have problems with meiotic division at early generations, in which interactions between non-homologous chromosomes often bring about aneuploidy and unbalanced gamete formation. We analyzed meiotic chromosome behaviors in pollen mother cells (PMCs) of allotetraploid and allodiploid F1 individuals of newly synthesized xBrassicoraphanus. Allotetraploid xBrassicoraphanus PMCs showed a normal diploid-like meiotic behavior. By contrast, allodiploid xBrassicoraphanus PMCs displayed abnormal segregation of chromosomes mainly due to the absence of homologous pairs. Notably, during early stages of meiosis I many of allodiploid xBrassicoraphanus chromosomes behave independently with few interactions between B. rapa and R. sativus chromosomes, forming many univalent chromosomes before segregation. Chromosomes were randomly assorted at later stages of meiosis, and tetrads with unequal numbers of chromosomes were formed at completion of meiosis. Immunolocalization of HEI10 protein mediating meiotic recombination revealed that COs were more frequent in synthetic allotetraploid xBrassicoraphanus than in allodiploid, but less than in the stabilized line. These findings suggest that structural dissimilarity between B. rapa and R. sativus chromosomes prevents non-homologous interactions between the parental chromosomes in allotetraploid xBrassicoraphanus, allowing normal diploid-like meiosis when homologous pairing partners are present. This study also suggests that CO suppression between non-homologous chromosomes is required for correct meiotic progression in newly synthesized allopolyploids, which is important for the formation of viable gametes and reproductive success in the hybrid progeny.

Whole genome duplication does not promote common modes of reproductive isolation in Trifolium pratense

PREMISE

Although polyploidy has been studied since the early 1900s, fundamental aspects of polyploid ecology and evolution remain unexplored. In particular, surprisingly little is known about how newly formed polyploids (neopolyploids) become demographically established. Models predict that most polyploids should go extinct within the first few generations as a result of reproductive disadvantages associated with being the minority in a primarily diploid population (i.e., the minority cytotype principle), yet polyploidy is extremely common. Therefore, a key goal in the study of polyploidy is to determine the mechanisms that promote polyploid establishment in nature. Because premating isolation is critical in order for neopolylpoids to avoid minority cytotype exclusion and thus facilitate establishment, we examined floral morphology and three common premating barriers to determine their importance in generating reproductive isolation of neopolyploids from diploids.

METHODS

We induced neopolyploidy in Trifolium pratense and compared their floral traits to the diploid progenitors. In addition to shifts in floral morphology, we examined three premating barriers: isolation by self-fertilization, flowering-time asynchrony, and pollinator-mediated isolation.

RESULTS

We found significant differences in the morphology of diploid and neopolyploid flowers, but these changes did not facilitate premating barriers that would generate reproductive isolation of neopolyploids from diploids. There was no difference in flowering phenology, pollinator visitation, or selfing between the cytotypes.

CONCLUSIONS

Our results indicate that barriers other than the ones tested in this study-such as geographic isolation, vegetative reproduction, and pistil-stigma incompatibilities-may be more important in facilitating isolation and establishment of neopolyploid T. pratense.

Homoeolog expression bias and expression level dominance (ELD) in four tissues of natural allotetraploid Brassica napus

Background

Allopolyploidy is widespread in angiosperms, and they can coordinate two or more different genomes through genetic and epigenetic modifications to exhibit stronger vigor and adaptability. To explore the changes in homologous gene expression patterns in the natural allotetraploid Brassica napus (A_nA_nC_nC_n) relative to its two diploid progenitors, B. rapa (A_rA_r) and B. oleracea (C_oC_o), after approximately 7500 years of domestication, the global gene pair expression patterns in four major tissues (stems, leaves, flowers and siliques) of these three species were analyzed using an RNA sequencing approach.

Results

The results showed that the ‘transcriptomic shock’ phenomenon was alleviated in natural B. napus after approximately 7500 years of natural domestication, and most differentially expressed genes (DEGs) in B. napus were downregulated relative to those in its two diploid progenitors. The KEGG analysis indicated that three pathways related to photosynthesis were enriched in both comparison groups (A_nA_nC_nC_n vs A_rA_r and A_nA_nC_nC_n vs C_oC_o), and these pathways were all downregulated in four tissues of B. napus. In addition, homoeolog expression bias and expression level dominance (ELD) in B. napus were thoroughly studied through analysis of expression levels of 27,609 B. rapa-B. oleracea orthologous gene pairs. The overwhelming majority of gene pairs (an average of 86.7%) in B. napus maintained their expression pattern in two diploid progenitors, and approximately 78.1% of the gene pairs showed expression bias with a preference toward the A subgenome. Overall, an average of 48, 29.7 and 22.3% homologous gene pairs exhibited additive expression, ELD and transgressive expression in B. napus, respectively. The ELD bias varies from tissue to tissue; specifically, more gene pairs in stems and siliques showed ELD-A, whereas the opposite was observed in leaves and flowers. More transgressive upregulation, rather than downregulation, was observed in gene pairs of B. napus.

Conclusions

In general, these results may provide a comprehensive understanding of the changes in homologous gene expression patterns in natural B. napus after approximately 7500 years of evolution and domestication and may enhance our understanding of allopolyploidy.

Extensive changes in gene expression and alternative splicing due to homoeologous exchange in rice segmental allopolyploids

We report rampant homoeologous exchanges in progenies of a newly synthesized rice segmental allotetraploid and demonstrate their consequences to changes of gene expression and alternative splicing. Allopolyploidization is recurrent across the tree of angiosperms and known as a driving evolutionary force in both plants and animals. A salient feature of allopolyploidization is the induction of homoeologous exchange (HE) events between the constituent subgenomes, which may in turn cause changes in gene expression, transcript alternative splicing, and phenotypic novelty. However, this issue has been poorly studied, largely because lack of a system in which the exact parentage donating the subgenomes is known and the HE events are occurring in real time. Here, we employed whole-genome re-sequencing and RNA-seq-based transcriptome profiling in four randomly chosen progeny individuals (at the 10th-selfed generation) of segmental allotetraploids that were constructed by colchicine-mediated whole-genome doubling of F1 hybrids between the two subspecies (japonica and indica) of Asian cultivated Oryza sativa. We show that rampant HE events occurred in these tetraploid individuals, which converted most of the otherwise heterozygous genomic regions into a homogenized state of one parental subgenome. We demonstrate that genes within these homogenized genomic regions in the tetraploids showed high frequencies of altered expression and enhanced alternative splicing relative to their counterparts in the corresponding diploid parents in the embryo tissue. Intriguingly, limited overlaps between the differentially expressed genes and the differential alternative spliced genes were identified, which were partitioned to distinctly enriched gene ontology terms. Together, our results indicate that HE is a major mechanism to rapidly generate novelty in gene expression and transcriptome diversity, which may facilitate phenotypic innovation in nascent allopolyploids and relevant to allopolyploid crop breeding.

Transcriptional regulation of anthocyanin biosynthesis in a high-anthocyanin resynthesized Brassica napus cultivar

BACKGROUND

Anthocyanins are plant secondary metabolites with key roles in attracting insect pollinators and protecting against biotic and abiotic stresses. They have potential health-promoting effects as part of the human diet. Anthocyanin biosynthesis has been elucidated in many species, enabling the development of anthocyanin-enriched fruits, vegetables, and grains; however, few studies have investigated Brassica napus anthocyanin biosynthesis.

RESULTS

We developed a high-anthocyanin resynthesized B. napus line, Rs035, by crossing anthocyanin-rich B. rapa (A genome) and B. oleracea (C genome) lines, followed by chromosome doubling. We identified and characterized 73 and 58 anthocyanin biosynthesis genes in silico in the A and C genomes, respectively; these genes showed syntenic relationships with 41 genes in Arabidopsis thaliana and B. napus. Among the syntenic genes, twelve biosynthetic and six regulatory genes showed transgressively higher expression in Rs035, and eight structural genes and one regulatory gene showed additive expression. We identified three early-, four late-biosynthesis pathways, three transcriptional regulator genes, and one transporter as putative candidates enhancing anthocyanin accumulation in Rs035. Principal component analysis and Pearson's correlation coefficients corroborated the contribution of these genes to anthocyanin accumulation.

CONCLUSIONS

Our study lays the foundation for producing high-anthocyanin B. napus cultivars. The resynthesized lines and the differentially expressed genes we have identified could be used to transfer the anthocyanin traits to other commercial rapeseed lines using molecular and conventional breeding.

Homoeolog expression bias and expression level dominance in resynthesized allopolyploid Brassica napus

Background

Allopolyploids require rapid genetic and epigenetic modifications to reconcile two or more sets of divergent genomes. To better understand the fate of duplicate genes following genomic mergers and doubling during allopolyploid formation, in this study, we explored the global gene expression patterns in resynthesized allotetraploid Brassica napus (AACC) and its diploid parents B. rapa (AA) and B. oleracea (CC) using RNA sequencing of leaf transcriptomes.

Results

We found that allopolyploid B. napus formation was accompanied by extensive changes (approximately one-third of the expressed genes) in the parental gene expression patterns (‘transcriptome shock’). Interestingly, the majority (85%) of differentially expressed genes (DEGs) were downregulated in the allotetraploid. Moreover, the homoeolog expression bias (relative contribution of homoeologs to the transcriptome) and expression level dominance (total expression level of both homoeologs) were thoroughly investigated by monitoring the expression of 23,766 B. oleracea-B. rapa orthologous gene pairs. Approximately 36.5% of the expressed gene pairs displayed expression bias with a slight preference toward the A-genome. In addition, 39.6, 4.9 and 9.0% of the expressed gene pairs exhibited expression level dominance (ELD), additivity expression and transgressive expression, respectively. The genome-wide ELD was also biased toward the A-genome in the resynthesized B. napus. To explain the ELD phenomenon, we compared the individual homoeolog expression levels relative to those of the diploid parents and found that ELD in the direction of the higher-expression parent can be explained by the downregulation of homoeologs from the dominant parent or upregulation of homoeologs from the nondominant parent; however, ELD in the direction of the lower-expression parent can be explained only by the downregulation of the nondominant parent or both homoeologs. Furthermore, Gene Ontology (GO) enrichment analysis suggested that the alteration in the gene expression patterns could be a prominent cause of the phenotypic variation between the newly formed B. napus and its parental species.

Conclusions

Collectively, our data provide insight into the rapid repatterning of gene expression at the beginning of Brassica allopolyploidization and enhance our knowledge of allopolyploidization processes.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4966-5) contains supplementary material, which is available to authorized users.

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.

Speciation Success of Polyploid Plants Closely Relates to the Regulation of Meiotic Recombination

Polyploidization is a widespread phenomenon, especially in flowering plants that have all undergone at least one event of whole genome duplication during their evolutionary history. Consequently, a large range of plants, including many of the world's crops, combines more than two sets of chromosomes originating from the same (autopolyploids) or related species (allopolyploids). Depending on the polyploid formation pathway, different patterns of recombination will be promoted, conditioning the level of heterozygosity. A polyploid population harboring a high level of heterozygosity will produce more genetically diverse progenies. Some of these individuals may show a better adaptability to different ecological niches, increasing their chance for successful establishment through natural selection. Another condition for young polyploids to survive corresponds to the formation of well-balanced gametes, assuring a sufficient level of fertility. In this review, we discuss the consequences of polyploid formation pathways, meiotic behavior and recombination regulation on the speciation success and maintenance of polyploid species.

Mixed-Ploidy Species: Progress and Opportunities in Polyploid Research

Mixed-ploidy species harbor a unique form of genomic and phenotypic variation that influences ecological interactions, facilitates genetic divergence, and offers insights into the mechanisms of polyploid evolution. However, there have been few attempts to synthesize this literature. We review here research on the cytotype distribution, diversity, and dynamics of intensively studied mixed-ploidy species and consider the implications for understanding mechanisms of polyploidization such as cytotype formation, establishment, coexistence, and post-polyploid divergence. In general, mixed-ploidy species are unevenly represented among families: they exhibit high cytotype diversity, often within populations, and frequently comprise rare and odd-numbered ploidies. Odd-ploidies often occur in association with asexuality. We highlight research hypotheses and opportunities that take advantage of the unique properties of ploidy variation.

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.

Evolutionary Dynamics of Unreduced Gametes

Unreduced gametes, which have the somatic (2n) chromosome number, are an important precursor to polyploid formation and apomixis. The product of irregularities in meiosis, 2n gametes are expected to be rare and deleterious in most natural populations, contrary to their wide taxonomic distribution and the prevalence of polyploidy. To better understand this discrepancy, we review contemporary evidence related to four aspects of 2n gamete dynamics in natural populations: (i) estimates of their frequency; (ii) their environmental and genetic determinants; (iii) adaptive and nonadaptive processes regulating their evolution; and (iv) factors regulating their union and production of polyploids in diploid populations. Aided by high-throughput methods of detection, these foci will advance our understanding of variation in 2n gametes within and among species, and their role in polyploid evolution.

Homoeologous chromosome pairing across the eukaryote phylogeny

During the past quarter century, molecular phylogenetic inferences have significantly resolved evolutionary relationships spanning the eukaryotic tree of life. With improved phylogenies in hand, the focus of systematics will continue to expand from estimating species relationships toward examining the evolution of specific, fundamental traits across the eukaryotic tree. Undoubtedly, this will expose knowledge gaps in the evolution of key traits, particularly with respect to non-model lineages. Here, we examine one such trait across eukaryotes-the regulation of homologous chromosome pairing during meiosis-as an illustrative example. Specifically, we present an overview of the breakdown of homologous chromosome pairing in model eukaryotes and provide a discussion of various meiotic aberrations that result in the failure of homolog recognition, with a particular focus on lineages with a history of hybridization and polyploidization, across major eukaryotic clades. We then explore what is known about these processes in natural and non-model eukaryotic taxa, thereby exposing disparities in our understanding of this key trait among non-model groups.

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.

Gene conversion events and variable degree of homogenization of rDNA loci in cultivars of Brassica napus

BACKGROUND AND AIMS

Brassica napus (AACC, 2n = 38, oilseed rape) is a relatively recent allotetraploid species derived from the putative progenitor diploid species Brassica rapa (AA, 2n = 20) and Brassica oleracea (CC, 2n = 18). To determine the influence of intensive breeding conditions on the evolution of its genome, we analysed structure and copy number of rDNA in 21 cultivars of B. napus, representative of genetic diversity.

METHODS

We used next-generation sequencing genomic approaches, Southern blot hybridization, expression analysis and fluorescence in situ hybridization (FISH). Subgenome-specific sequences derived from rDNA intergenic spacers (IGS) were used as probes for identification of loci composition on chromosomes.

KEY RESULTS

Most B. napus cultivars (18/21, 86 %) had more A-genome than C-genome rDNA copies. Three cultivars analysed by FISH ('Darmor', 'Yudal' and 'Asparagus kale') harboured the same number (12 per diploid set) of loci. In B. napus 'Darmor', the A-genome-specific rDNA probe hybridized to all 12 rDNA loci (eight on the A-genome and four on the C-genome) while the C-genome-specific probe showed weak signals on the C-genome loci only. Deep sequencing revealed high homogeneity of arrays suggesting that the C-genome genes were largely overwritten by the A-genome variants in B. napus 'Darmor'. In contrast, B. napus 'Yudal' showed a lack of gene conversion evidenced by additive inheritance of progenitor rDNA variants and highly localized hybridization signals of subgenome-specific probes on chromosomes. Brassica napus 'Asparagus kale' showed an intermediate pattern to 'Darmor' and 'Yudal'. At the expression level, most cultivars (95 %) exhibited stable A-genome nucleolar dominance while one cultivar ('Norin 9') showed co-dominance.

CONCLUSIONS

The B. napus cultivars differ in the degree and direction of rDNA homogenization. The prevalent direction of gene conversion (towards the A-genome) correlates with the direction of expression dominance indicating that gene activity may be needed for interlocus gene conversion.

Flowering Time Gene Variation in Brassica Species Shows Evolutionary Principles

Flowering time genes have a strong influence on successful reproduction and life cycle adaptation. However, their regulation is highly complex and only well understood in diploid model systems. For crops with a polyploid background from the genus Brassica, data on flowering time gene variation are scarce, although indispensable for modern breeding techniques like marker-assisted breeding. We have deep-sequenced all paralogs of 35 Arabidopsis thaliana flowering regulators using Sequence Capture followed by Illumina sequencing in two selected accessions of the vegetable species Brassica rapa and Brassica oleracea, respectively. Using these data, we were able to call SNPs, InDels and copy number variations (CNVs) for genes from the total flowering time network including central flowering regulators, but also genes from the vernalisation pathway, the photoperiod pathway, temperature regulation, the circadian clock and the downstream effectors. Comparing the results to a complementary data set from the allotetraploid species Brassica napus, we detected rearrangements in B. napus which probably occurred early after the allopolyploidisation event. Those data are both a valuable resource for flowering time research in those vegetable species, as well as a contribution to speciation genetics.

Recreating Stable Brachypodium hybridum Allotetraploids by Uniting the Divergent Genomes of B. distachyon and B. stacei

Brachypodium hybridum (2n = 30) is a natural allopolyploid with highly divergent sub-genomes derived from two extant diploid species, B. distachyon (2n = 10) and B. stacei (2n = 20) that differ in chromosome evolution and number. We created synthetic B. hybridum allotetraploids by hybridizing various lines of B. distachyon and B. stacei. The initial amphihaploid F1 interspecific hybrids were obtained at low frequencies when B. distachyon was used as the maternal parent (0.15% or 0.245% depending on the line used) and were sterile. No hybrids were obtained from reciprocal crosses or when autotetraploids of the parental species were crossed. Colchicine treatment was used to double the genome of the F1 amphihaploid lines leading to allotetraploids. The genome-doubled F1 plants produced a few S1 (first selfed generation) seeds after self-pollination. S1 plants from one parental combination (Bd3-1×Bsta5) were fertile and gave rise to further generations whereas those of another parental combination (Bd21×ABR114) were sterile, illustrating the importance of the parental lineages crossed. The synthetic allotetraploids were stable and resembled the natural B. hybridum at the phenotypic, cytogenetic and genomic levels. The successful creation of synthetic B. hybridum offers the possibility to study changes in genome structure and regulation at the earliest stages of allopolyploid formation in comparison with the parental species and natural B. hybridum.

Extraction of the Constituent Subgenomes of the Natural Allopolyploid Rapeseed (Brassica napus L.)

As the dynamic nature of progenitor genomes accompanies the speciation by interspecific hybridization, the extraction of the constituent subgenome(s) from a natural allopolyploid species of long history and then restitution of the progenitor(s) provides the unique opportunity to study the genome evolution and interplay. Herein, the A subgenome from the allotetraploid oilseed rape (Brassica napus L., AACC) was extracted through inducing the preferential elimination of C-subgenome chromosomes in intertribal crosses and the progenitor B. rapa was restituted (RBR). Then by crossing and backcrossing RBR with B. napus donor, the C subgenome was in situ dissected by adding each of its nine chromosomes to the extracted A subgenome and establishing the whole set of monosonic alien addition lines (MAALs). RBR from spring-type B. napus genotype "Oro" expressed a phenotype resembling some type of B. rapa never observed before, but showed a winter-type flowering habit. This RBR had weaker growth vigor and suffered more seriously from biotic and abiotic stresses compared with Oro. The phenotypes specific for these MAALs showed the location of the related genes on the particular C-subgenome chromosomes. These MAALs exhibited obviously different frequencies in homeologous pairing and transmission of additional C-subgenome chromosomes, which were associated with the distinct degrees of their relatedness, and even with the possible genetic regulation for meiotic pairing evolved in B. napus Finally, large scaffolds undetermined for sequence assembly of B. napus were anchored to specific C-subgenome chromosomes using MAALs.

Polyploidy: Pitfalls and paths to a paradigm

Investigators have long searched for a polyploidy paradigm-rules or principles that might be common following polyploidization (whole-genome duplication, WGD). Here we attempt to integrate what is known across the more thoroughly investigated polyploid systems on topics ranging from genetics to ecology. We found that while certain rules may govern gene retention and loss, systems vary in the prevalence of gene silencing vs. homeolog loss, chromosomal change, the presence of a dominant genome (in allopolyploids), and the relative importance of hybridization vs. genome doubling per se. In some lineages, aspects of polyploidization are repeated across multiple origins, but in other species multiple origins behave more stochastically in terms of genetic and phenotypic change. Our investigation also reveals that the path to synthesis is hindered by numerous gaps in our knowledge of even the best-known systems. Particularly concerning is the absence of linkage between genotype and phenotype. Moreover, most recent studies have focused on the genetic and genomic attributes of polyploidy, but rarely is there an ecological or physiological context. To promote a path to a polyploidy paradigm (or paradigms), we propose a major community goal over the next 10-20 yr to fill the gaps in our knowledge of well-studied polyploids. Before a meaningful synthesis is possible, more complete data sets are needed for comparison-systems that include comparable genetic, genomic, chromosomal, proteomic, as well as morphological, physiological, and ecological data. Also needed are more natural evolutionary model systems, as most of what we know about polyploidy continues to come from a few crop and genetic models, systems that often lack the ecological context inherent in natural systems and necessary for understanding the drivers of biodiversity.

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.

Direct vs. indirect effects of whole-genome duplication on prezygotic isolation in Chamerion angustifolium: Implications for rapid speciation

PREMISE OF THE STUDY

The depiction of polyploid speciation as instantaneous implies that strong prezygotic and postzygotic isolation form as a direct result of whole-genome duplication. However, the direct vs. indirect contributions of genome duplication to phenotypic divergence and prezygotic isolation are rarely quantified across multiple reproductive barriers.

METHODS

We compared the phenotypic differences between diploid and both naturally occurring and synthesized tetraploids (neotetraploids) of the plant Chamerion angustifolium. Using this information and additional published values for this species, we compared the magnitude of isolation (ecological, flowering, pollinator, and gametic) between diploids and natural-occurring tetraploids to that between diploids and neotetraploids.

KEY RESULTS

Differences among ploidy cytotypes were observed for eight of 12 vegetative and reproductive traits measured. Neotetraploids resembled diploids but differed from natural tetraploids with respect to four traits, including flowering time and plant height. Diploid-neotetraploid (2x-4xneo) experimental arrays exhibited lower pollinator fidelity to cytotype and seed set compared with 2x-4xnat arrays. Based on these results and published evidence, reproductive isolation between diploids and neotetraploids across all four life stages averaged 0.48 and deviated significantly from that between diploids and natural tetraploids (RI = 0.96).

CONCLUSIONS

Genome duplication causes phenotypic shifts and contributes directly to prezygotic isolation for some barriers (gametic isolation) but cannot account for the cumulative isolation from diploids observed in natural tetraploids. Therefore, conditions for species formation through genome duplication are not necessarily instantaneous and selection to strengthen prezygotic barriers in young polyploids is critical for the establishment of polyploid species in sympatry.

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.

Linkage and regional association analysis reveal two new tightly-linked major-QTLs for pod number and seed number per pod in rapeseed (Brassica napus L.)

To facilitate the pseudochromosomes assembly and gene cloning in rapeseed, we developed a reference genetic population/map (named BnaZNF2) from two sequenced cultivars, Zhongshuang11 and No.73290, those exhibit significant differences in many traits, particularly yield components. The BnaZNF2 genetic map exhibited perfect collinearity with the physical map of B. napus, indicating its high quality. Comparative mapping revealed several genomic rearrangements between B. napus and B. rapa or B. oleracea. A total of eight and 16 QTLs were identified for pod number and seed number per pod, respectively, and of which three and five QTLs are identical to previously identified ones, whereas the other five and 11 are novel. Two new major QTL respectively for pod number and seed number per pod, qPN.A06-1 and qSN.A06-1 (R(2 )= 22.8% and 32.1%), were colocalised with opposite effects, and only qPN.A06-1 was confirmed and narrowed by regional association analysis to 180 kb including only 33 annotated genes. Conditional QTL analysis and subsequent NILs test indicated that tight linkage, rather than pleiotropy, was the genetic causation of their colocalisation. Our study demonstrates potential of this reference genetic population/map for precise QTL mapping and as a base for positional gene cloning in rapeseed.

Digital gene expression analysis of gene expression differences within Brassica diploids and allopolyploids

BACKGROUND

Brassica includes many successfully cultivated crop species of polyploid origin, either by ancestral genome triplication or by hybridization between two diploid progenitors, displaying complex repetitive sequences and transposons. The U's triangle, which consists of three diploids and three amphidiploids, is optimal for the analysis of complicated genomes after polyploidization. Next-generation sequencing enables the transcriptome profiling of polyploids on a global scale.

RESULTS

We examined the gene expression patterns of three diploids (Brassica rapa, B. nigra, and B. oleracea) and three amphidiploids (B. napus, B. juncea, and B. carinata) via digital gene expression analysis. In total, the libraries generated between 5.7 and 6.1 million raw reads, and the clean tags of each library were mapped to 18547-21995 genes of B. rapa genome. The unambiguous tag-mapped genes in the libraries were compared. Moreover, the majority of differentially expressed genes (DEGs) were explored among diploids as well as between diploids and amphidiploids. Gene ontological analysis was performed to functionally categorize these DEGs into different classes. The Kyoto Encyclopedia of Genes and Genomes analysis was performed to assign these DEGs into approximately 120 pathways, among which the metabolic pathway, biosynthesis of secondary metabolites, and peroxisomal pathway were enriched. The non-additive genes in Brassica amphidiploids were analyzed, and the results indicated that orthologous genes in polyploids are frequently expressed in a non-additive pattern. Methyltransferase genes showed differential expression pattern in Brassica species.

CONCLUSION

Our results provided an understanding of the transcriptome complexity of natural Brassica species. The gene expression changes in diploids and allopolyploids may help elucidate the morphological and physiological differences among Brassica species.

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.

Transposon variation by order during allopolyploidisation between Brassica oleracea and Brassica rapa

Although many studies have shown that transposable element (TE) activation is induced by hybridisation and polyploidisation in plants, much less is known on how different types of TE respond to hybridisation, and the impact of TE-associated sequences on gene function. We investigated the frequency and regularity of putative transposon activation for different types of TE, and determined the impact of TE-associated sequence variation on the genome during allopolyploidisation. We designed different types of TE primers and adopted the Inter-Retrotransposon Amplified Polymorphism (IRAP) method to detect variation in TE-associated sequences during the process of allopolyploidisation between Brassica rapa (AA) and Brassica oleracea (CC), and in successive generations of self-pollinated progeny. In addition, fragments with TE insertions were used to perform Blast2GO analysis to characterise the putative functions of the fragments with TE insertions. Ninety-two primers amplifying 548 loci were used to detect variation in sequences associated with four different orders of TE sequences. TEs could be classed in ascending frequency into LTR-REs, TIRs, LINEs, SINEs and unknown TEs. The frequency of novel variation (putative activation) detected for the four orders of TEs was highest from the F1 to F2 generations, and lowest from the F2 to F3 generations. Functional annotation of sequences with TE insertions showed that genes with TE insertions were mainly involved in metabolic processes and binding, and preferentially functioned in organelles. TE variation in our study severely disturbed the genetic compositions of the different generations, resulting in inconsistencies in genetic clustering. Different types of TE showed different patterns of variation during the process of allopolyploidisation.

Genome-wide disruption of gene expression in allopolyploids but not hybrids of rice subspecies

Hybridization and polyploidization are prominent processes in plant evolution. Hybrids and allopolyploids typically exhibit radically altered gene expression patterns relative to their parents, a phenomenon termed “transcriptomic shock.” To distinguish the effects of hybridization from polyploidization on coregulation of divergent alleles, we analyzed expression of parental copies (homoeologs) of 11,608 genes using RNA-seq-based transcriptome profiling in reciprocal hybrids and tetraploids constructed from subspecies japonica and indica of Asian rice (Oryza sativa L.). The diploid hybrids and their derived allopolyploids differ dramatically in morphology, despite having the same suite of genes and genic proportions. Allelic and homoeolog-specific transcripts were unequivocally diagnosed in the hybrids and tetraploids based on parent-specific SNPs. Compared with the in silico hybrid (parental mix), the range of progenitor expression divergence was significantly reduced in both reciprocally generated F1 hybrids, presumably due to the ameliorating effects of a common trans environment on divergent cis-factors. In contrast, parental expression differences were greatly elaborated at the polyploid level, which we propose is a consequence of stoichiometric disruptions associated with the numerous chromosomal packaging and volumetric changes accompanying nascent polyploidy. We speculate that the emergent property of “whole genome doubling” has repercussions that reverberate throughout the transcriptome and downstream, ultimately generating altered phenotypes. This perspective may yield insight into the nature of adaptation and the origin of evolutionary novelty accompanying polyploidy.

Extensive tRNA gene changes in synthetic Brassica napus

Allopolyploidization, where two species come together to form a new species, plays a major role in speciation and genome evolution. Transfer RNAs (abbreviated tRNA) are typically 73-94 nucleotides in length, and are indispensable in protein synthesis, transferring amino acids to the cell protein synthesis machinery (ribosome). To date, the regularity and function of tRNA gene sequence variation during the process of allopolyploidization have not been well understood. In this study, the inter-tRNA gene corresponding to tRNA amplification polymorphism method was used to detect changes in tRNA gene sequences in the progeny of interspecific hybrids between Brassica rapa and B. oleracea, mimicking the original B. napus (canola) species formation event. Cluster analysis showed that tRNA gene variation during allopolyploidization did not appear to have a genotypic basis. Significant variation occurred in the early generations of synthetic B. napus (F1 and F2 generations), but fewer alterations were observed in the later generation (F3). The variation-prone tRNA genes tended to be located in AT-rich regions. BlastN analysis of novel tRNA gene variants against a Brassica genome sequence database showed that the variation of these tRNA-gene-associated sequences in allopolyploidization might result in variation of gene structure and function, e.g., metabolic process and transport.

Nucleolar dominance and different genome behaviors in hybrids and allopolyploids

Many plants are allopolyploids with different nuclear genomes from two or more progenitors, but cytoplasmic genomes typically inherited from the female parent. The importance of this speciation mechanism has stimulated the extensive investigations of genetic consequences of genome mergers in several experimental systems during last 20 years. The dynamic nature of polyploid genomes is recognized, and widespread changes to gene expression are revealed by transcriptomic analysis. These progresses show different stabilities of parental genomes and their unequal contributions to the transcriptome, proteome, and phenotype. We review the results in systems where extensive genetic analyses have been conducted and propose possible mechanisms for biased behavior of parental genomes in allopolyploids, including the role of nucleolar dominance. It is hypothesized that the novel ribosomes with rRNAs from uniparental genome and the ribosomal proteins of biparental origins have some impacts on the biased cellular and genetic behaviors of parental genomes in hybrids and allopolyploids.

A large-scale introgression of genomic components of Brassica rapa into B. napus by the bridge of hexaploid derived from hybridization between B. napus and B. oleracea

Brassica rapa (AA) has been used to widen the genetic basis of B. napus (AACC), which is a new but important oilseed crop worldwide. In the present study, we have proposed a strategy to develop new type B. napus carrying genomic components of B. rapa by crossing B. rapa with hexaploid (AACCCC) derived from B. napus and B. oleracea (CC). The hexaploid exhibited large flowers and high frequency of normal chromosome segregation, resulting in good seed set (average of 4.48 and 12.53 seeds per pod by self and open pollination, respectively) and high pollen fertility (average of 87.05 %). It was easy to develop new type B. napus by crossing the hexaploid with 142 lines of B. rapa from three ecotype groups, with the average crossability of 9.24 seeds per pod. The genetic variation of new type B. napus was diverse from that of current B. napus, especially in the A subgenome, revealed by genome-specific simple sequence repeat markers. Our data suggest that the strategy proposed here is a large-scale and highly efficient method to introgress genomic components of B. rapa into B. napus.

Revisiting the dioecy-polyploidy association: alternate pathways and research opportunities

The evolutionary transition from hermaphroditism (combined sexes) to dioecy (separate sexes) is associated with whole genome duplication (polyploidy) in several flowering plant genera. Moreover, there is evidence for transitions in the opposite direction, i.e. a loss of dioecy with an increase in ploidy. Here, we review evidence for these associations, synthesize previous ideas on the mechanism underlying the patterns and explore alternative pathways. Specifically, we examine potential ecological and genetic mechanisms, differentiated by whether ploidy or gender (functional sex expression of the plant) changes are the primary cause and whether the effect is direct or indirect. An analysis of 22 genera variable for both ploidy and gender indicates that gender monomorphism (hermaphroditism, monoecy) is more common among diploid than polyploid species, whereas gender dimorphism (dioecy, gynodioecy, subdioecy) is more frequent among polyploid species. The transition from diploid hermaphroditic to polyploid gender-dimorphic taxa may arise directly through changes in gender as a result of genome duplication through genomic rearrangements or homeologous recombination, or changes in gender may result in increased unreduced gamete production leading to polyploid formation. Alternatively, the transition may occur through the indirect effects of genome duplication on mating system and inbreeding depression, which favor selection for unisexuality, or habitat shifts associated with unisexuality may simultaneously cause increased unreduced gamete production. Novel mechanisms for transitions in the opposite direction (from dioecy to hermaphroditism with increase in ploidy) include disruption of genetic sex determination and the benefits of reproductive assurance. We highlight key questions requiring further attention and promising approaches for answering them and better clarifying the genesis of sexual system polyploidy associations. See also the sister article focusing on animals by Wertheim et al. in this themed issue.

Lessons from natural and artificial polyploids in higher plants

Polyploidy in higher plants is a major source of genetic novelty upon which selection may act to drive evolution, as evidenced by the widespread success of polyploid species in the wild. However, research into the effects of polyploidy can be confounded by the entanglement of several processes: genome duplication, hybridisation (allopolyploidy is frequent in plants) and subsequent evolution. The discovery of the chemical agent colchicine, which can be used to produce artificial polyploids on demand, has enabled scientists to unravel these threads and understand the complex genomic changes involved in each. We present here an overview of lessons learnt from studies of natural and artificial polyploids, and from comparisons between the 2, covering basic cellular and metabolic consequences through to alterations in epigenetic gene regulation, together with 2 in-depth case studies in Senecio and Glycine. See also the sister article focusing on animals by Arai and Fujimoto in this themed issue.