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

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

Molecular characterization of intergeneric hybrids between Malus and Pyrus

Apple (Malus) and pear (Pyrus) are economically important fruit crops well known for their unique textures, flavours, and nutritional qualities. Both genera are characterised by a distinct pattern of secondary metabolites, which directly affect not only resistance to certain diseases, but also have significant impacts on the flavour and nutritional value of the fruit. The identical chromosome numbers, similar genome size, and their recent divergence date, together with DNA markers have shown that apple and pear genomes are highly co-linear. This study utilized comparative genomic approaches, including simple sequence repeats, high resolution single nucleotide polymorphism melting analysis, and single nucleotide polymorphism chip analysis to identify genetic differences among hybrids of Malus and Pyrus, and F2 offspring. This research has demonstrated and validated that these three marker types, along with metabolomics analysis are very powerful tools to detect and confirm hybridity of progeny derived from crosses between apple and pear in both cross directions. Furthermore, this work analysed the genus-specific metabolite patterns and the resistance to fire blight (Erwinia amylovora) in progeny. The findings of this work will enhance and accelerate the breeding of novel tree fruit crops that benefit producers and consumers, by enabling marker assisted selection of desired traits introgressed between pear and apple.

Exploring the gene pool of Brassica napus by genomics-based approaches

De novo allopolyploidization in Brassica provides a very successful model for reconstructing polyploid genomes using progenitor species and relatives to broaden crop gene pools and understand genome evolution after polyploidy, interspecific hybridization and exotic introgression. B. napus (AACC), the major cultivated rapeseed species and the third largest oilseed crop in the world, is a young Brassica species with a limited genetic base resulting from its short history of domestication, cultivation, and intensive selection during breeding for target economic traits. However, the gene pool of B. napus has been significantly enriched in recent decades that has been benefit from worldwide effects by the successful introduction of abundant subgenomic variation and novel genomic variation via intraspecific, interspecific and intergeneric crosses. An important question in this respect is how to utilize such variation to breed crops adapted to the changing global climate. Here, we review the genetic diversity, genome structure, and population-level differentiation of the B. napus gene pool in relation to known exotic introgressions from various species of the Brassicaceae, especially those elucidated by recent genome-sequencing projects. We also summarize progress in gene cloning, trait-marker associations, gene editing, molecular marker-assisted selection and genome-wide prediction, and describe the challenges and opportunities of these techniques as molecular platforms to exploit novel genomic variation and their value in the rapeseed gene pool. Future progress will accelerate the creation and manipulation of genetic diversity with genomic-based improvement, as well as provide novel insights into the neo-domestication of polyploid crops with novel genetic diversity from reconstructed genomes.

Gene Expression Changes During the Allo-/Deallopolyploidization Process of Brassica napus

Gene expression changes due to allopolyploidization have been extensively studied in plants over the past few decades. Nearly all these studies focused on comparing the changes before and after genome merger. In this study, we used the uniquely restituted Brassica rapa (RBR, A_eA_e, 2n = 20) obtained from Brassica napus (A_nA_nC_nC_n, 2n = 38) to analyze the gene expression changes and its potential mechanism during the process of allo-/deallopolyploidization. RNA-seq-based transcriptome profiling identified a large number of differentially expressed genes (DEGs) between RBR and natural B. rapa (A_rA_r), suggesting potential effects of allopolyploidization/domestication of AA component of B. napus at the tetrapolyploid level. Meanwhile, it was revealed that up to 20% of gene expressions were immediately altered when compared with those in the A_n-subgenome. Interestingly, one fifth of these changes are in fact indicative of the recovery of antecedent gene expression alternations occurring since the origin of B. napus and showed association with homoeologous expression bias between A_n and C_n subgenomes. Enrichment of distinct gene ontology (GO) categories of the above sets of genes further indicated potential functional cooperation of the A_n and C_n subgenome of B. napus. Whole genome methylation analysis revealed a small number of DEGs were identified in the differentially methylated regions.

Transcriptional Aneuploidy Responses of Brassica rapa-oleracea Monosomic Alien Addition Lines (MAALs) Derived From Natural Allopolyploid B. napus

Establishing the whole set of aneuploids, for one naturally evolved allopolyploid species, provides a unique opportunity to elucidate the transcriptomic response of the constituent subgenomes to serial aneuploidy. Previously, the whole set of monosomic alien addition lines (MAALs, C1-C9) with each of the nine C subgenome chromosomes, added to the extracted A subgenome, was developed in the context of the allotetraploid Brassica napus donor “Oro,” after the restitution of the ancestral B. rapa (RBR Oro) was realized. Herein, transcriptomic analysis using high-throughput technology was conducted to detect gene expression alterations in these MAALs and RBR. Compared to diploid RBR, the genes of all of the MAALs showed various degrees of dysregulated expressions that resulted from cis effects and more prevailing trans effects. In addition, the trans-effect on gene expression in MAALs increased with higher levels of homology between the recipient A subgenome and additional C subgenome chromosomes, instead of gene numbers of extra chromosomes. A total of 10 trans-effect dysregulated genes, among all pairwise comparisons, were mainly involved in the function of transporter activity. Furthermore, highly expressed genes were more prone to downregulation and vice-versa, suggesting a common trend for transcriptional pattern responses to aneuploidy. These results provided a comprehensive insight of the impact of gene expression of individual chromosomes, in one subgenome, on another intact subgenome for one allopolyploid with a long evolutionary history.

Genome-Wide Gene Expression Disturbance by Single A1/C1 Chromosome Substitution in Brassica rapa Restituted From Natural B. napus

Alien chromosome substitution (CS) lines are treated as vital germplasms for breeding and genetic mapping. Previously, a whole set of nine Brassica rapa-oleracea monosonic alien addition lines (MAALs, C1-C9) was established in the background of natural B. napus genotype "Oro," after the restituted B. rapa (RBR) for Oro was realized. Herein, a monosomic substitution line with one alien C1 chromosome (Cs1) in the RBR complement was selected in the progenies of MAAL C1 and RBR, by the PCR amplification of specific gene markers and fluorescence in situ hybridization. Cs1 exhibited the whole plant morphology similar to RBR except for the defective stamens without fertile pollen grains, but it produced some seeds and progeny plants carrying the C1 chromosome at high rate besides those without the alien chromosome after pollinated by RBR. The viability of the substitution and its progeny for the RBR diploid further elucidated the functional compensation between the chromosome pairs with high homoeology. To reveal the impact of such aneuploidy on genome-wide gene expression, the transcriptomes of MAAL C1, Cs1 and euploid RBR were analyzed. Compared to RBR, Cs1 had sharply reduced gene expression level across chromosome A1, demonstrating the loss of one copy of A1 chromosome. Both additional chromosome C1 in MAAL and substitutional chromosome C1 in Cs1 caused not only cis-effect but also prevalent trans-effect differentially expressed genes. A dominant gene dosage effects prevailed among low expressed genes across chromosome A1 in Cs1, and moreover, dosage effects for some genes potentially contributed to the phenotype deviations. Our results provided novel insights into the transcriptomic perturbation and gene dosage effects on phenotype in CS related to one naturally evolved allopolyploid.

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.

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.

Genome-Wide Gene Expressions Respond Differently to A-subgenome Origins in Brassica napus Synthetic Hybrids and Natural Allotetraploid

The young allotetraploid Brassica napus (2n = 38, AACC) is one of models to study genomic responses to allopolyploidization. The extraction of AA component from natural B. napus and then restitution of progenitor B. rapa should provide a unique opportunity to reveal the genome interplay for gene expressions during the evolution. Herein, B. napus hybrids (2n = 19, AC) between the extracted and extant B. rapa (2n = 20, AA) and the same B. oleracea genotype (2n = 18, CC) were studied by RNA-seq and compared with natural B. napus donor, to reveal the gene expression changes from hybridization and domestication and the effects of A genome with different origins. Upon the initial merger of two diploid genomes, additive gene expression was prevalent in these two hybrids, for non-additively expressed genes only represented a small portion of total expressed genes. A high proportion of genes exhibited expression level dominance, with no preference to either of the parental genomes. Comparison of homoeolog expressions also showed no bias toward any genomes and the parental expression patterns were often maintained in the hybrids and natural allotetraploids. Although, the overall patterns of gene expression were highly conserved between two hybrids, the extracted B. rapa responded less and appeared more compatible for hybridization than the extant B. rapa. Our results suggested that expression level dominance and homoeolog expressions bias were balanced at the initial stage of genome merger, and such balance were largely maintained during the domestication of B. napus, despite the increased extent over time.

Transcriptomic analysis reveals differential gene expressions for cell growth and functional secondary metabolites in induced autotetraploid of Chinese woad (Isatis indigotica Fort.)

The giant organs and enhanced concentrations of secondary metabolites realized by autopolyploidy are attractive for breeding the respective medicinal and agricultural plants and studying the genetic mechanisms. The traditional medicinal plant Chinese woad (Isatis indigotica Fort., 2n = 2x = 14) is now still largely used for the diseases caused by bacteria and viruses in China. In this study, its autopolyploids (3x, 4x) were produced and characterized together with the 2x donor for their phenotype and transcriptomic alterations by using high-throughput RNA sequencing. With the increase of genome dosage, the giantism in cells and organs was obvious and the photosynthetic rate was higher. The 4x plants showed predominantly the normal meiotic chromosome pairing (bivalents and quadrivalents) and equal segregation and then produced the majority of 4x progeny. The total 70136 All-unigenes were de novo assembled, and 56,482 (80.53%) unigenes were annotated based on BLASTx searches of the public databases. From pair-wise comparisons between transcriptomic data of 2x, 3x, 4x plants, 1856 (2.65%)(2x vs 4x), 693(0.98%)(2x vs 3x), 1045(1.48%)(3x vs 4x) unigenes were detected to differentially expressed genes (DEGs), including both up- and down-regulated ones. These DEGs were mainly involved in cell growth (synthesis of expansin and pectin), cell wall organization, secondary metabolite biosynthesis, response to stress and photosynthetic pathways. The up-regulation of some DEGs for metabolic pathways of functional compounds in the induced autotetraploids substantiates the promising new type of this medicinal plant with the increased biomass and targeted metabolites.

Genome-wide gene expression perturbation induced by loss of C2 chromosome in allotetraploid Brassica napus L

Aneuploidy with loss of entire chromosomes from normal complement disrupts the balanced genome and is tolerable only by polyploidy plants. In this study, the monosomic and nullisomic plants losing one or two copies of C2 chromosome from allotetraploid Brassica napus L. (2n = 38, AACC) were produced and compared for their phenotype and transcriptome. The monosomics gave a plant phenotype very similar to the original donor, but the nullisomics had much smaller stature and also shorter growth period. By the comparative analyses on the global transcript profiles with the euploid donor, genome-wide alterations in gene expression were revealed in two aneuploids, and their majority of differentially expressed genes (DEGs) resulted from the trans-acting effects of the zero and one copy of C2 chromosome. The higher number of up-regulated genes than down-regulated genes on other chromosomes suggested that the genome responded to the C2 loss via enhancing the expression of certain genes. Particularly, more DEGs were detected in the monosomics than nullisomics, contrasting with their phenotypes. The gene expression of the other chromosomes was differently affected, and several dysregulated domains in which up- or downregulated genes obviously clustered were identifiable. But the mean gene expression (MGE) for homoeologous chromosome A2 reduced with the C2 loss. Some genes and their expressions on C2 were correlated with the phenotype deviations in the aneuploids. These results provided new insights into the transcriptomic perturbation of the allopolyploid genome elicited by the loss of individual chromosome.

Development of a complete set of monosomic alien addition lines between Brassica napus and Isatis indigotica (Chinese woad)

A complete set of monosomic alien addition lines of Brassica napus with one of the seven chromosomes of Isatis indigotica and the recombinant mitochondria was developed and characterized. Monosomic alien addition lines (MAALs) are valuable for elucidating the genome structure and transferring the useful genes and traits in plant breeding. Isatis indigotica (Chinese woad, 2n = 14, II) in Isatideae tribe of Brassicaceae family has been widely cultivated as a medicinal and dye plant in China. Herein, the intertribal somatic hybrid (2n = 52, AACCII) between B. napus cv. Huashuang 3 (2n = 38, AACC) and I. indigotica produced previously was backcrossed recurrently to parental B. napus, and 32 MAAL plants were isolated. Based on their phenotype, 5S and 45S rDNA loci and chromosome-specific SSR markers, these MAALs were classified into seven groups corresponding to potential seven types of MAALs carrying one of the seven I. indigotica chromosomes. One of the MAALs could be distinguishable by expressing the brown anthers of I. indigotica, other two hosted the chromosome with 5S or 45S rDNA locus, but the remaining four were identifiable by SSR markers. The simultaneous detection of the same SSR maker and gene locus in different MAALs revealed the paralogs on the chromosomes involved. The recombinant mitochondrial genome in MAALs was likely related with their male sterility with carpellody stamens, while the MAAL with normal brown anthers probably carried the restoring gene for the male sterility. The complete set of MAALs should be useful for exploiting the I. indigotica genome and for promoting the introgression of valuable genes to B. napus.

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.