Qi Pan
Bin Zhu
Dawei Zhang1†
Chaobo Tong
Xianhong Ge
Shengyi Liu
Zaiyun Li- 1National Key Laboratory of Crop Genetic Improvement, Key Laboratory of Rapeseed Genetics and Breeding of Agriculture Ministry of China, College of Plant Science and Technology, Huazhong Agricultural University, Wuhan, China
- 2Key Laboratory of Biology and Genetic Improvement of Oil Crops, Ministry of Agriculture, Wuhan, China
- 3Oil Crops Research Institute of the Chinese Academy of Agricultural Sciences, Wuhan, China
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, AeAe, 2n = 20) obtained from Brassica napus (AnAnCnCn, 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 (ArAr), 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 An-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 An and Cn subgenomes. Enrichment of distinct gene ontology (GO) categories of the above sets of genes further indicated potential functional cooperation of the An and Cn subgenome of B. napus. Whole genome methylation analysis revealed a small number of DEGs were identified in the differentially methylated regions.
Introduction
Polyploidy (whole-genome duplication) is an important speciation mechanism for all eukaryotes, especially in higher plant evolution. Allopolyploid speciation results from interspecific hybridization and genome doubling. During the formation and evolution of allopolyploids, genome-wide genetic, epigenetic, and gene expression changes might occur, when compared with the parents (Chen, 2007; Jackson and Chen, 2010; Madlung and Wendel, 2013). On one hand, these alternations undoubtedly promote new traits formation, thereby increasing adaptability. On the other hand, heterozygosity and intergenomic interactions in allopolyploids may also result in heterosis, giving rise to phenotypic variations and increased growth vigor (Comai, 2005; Chen, 2010).
It is now realized that the genome of the allopolyploid is not just the sum of the two parental genomes. Rapid and extensive genomic changes have been revealed in the past decades in various synthetic allopolyploids or newly-formed natural polyploids (Madlung and Wendel, 2013), including sequence loss (Shaked et al., 2001; Han et al., 2005; Ozkan and Feldman, 2009; Tate et al., 2009a), gene conversion (Kovarik et al., 2008; Salmon et al., 2010), rDNA loci changes (Joly et al., 2004; Pontes et al., 2004; Lim et al., 2008), activation of transposons (Kashkush et al., 2003; Madlung et al., 2005; Chen et al., 2008; Khasdan et al., 2010; Parisod et al., 2010; Hegarty et al., 2011), and recombination between homoeologous chromosomes (Gaeta and Chris Pires, 2010; Salmon et al., 2010; Szadkowski et al., 2010). Specifically, studies have provided strong evidence that homoeologous recombination contributes much to genetic changes and might be an important driver of the genome evolution of B. napus (Gaeta and Chris Pires, 2010; Xiong and Pires, 2011) and coffee (Lashermes et al., 2014). Associated with these genomic changes, variations in DNA methylation between parental species and allopolyploids has also been reported in Arabidopsis (Madlung et al., 2005; Beaulieu et al., 2009), wheat (Zhao et al., 2011), Brassica (Lukens et al., 2006; Xu et al., 2009; Ksiazczyk et al., 2011), dandelion (Verhoeven et al., 2010), and many other species.
Genetic and epigenetic changes will lead to new expression patterns of duplicated genes in polyploids that might be different from those in each parental genome. The total expression level of the duplicated genes can show deviations from parental additivity (the value of the gene in the polyploids is equal to the mid-parent value, MPV). Such nonadditive gene expression patterns have been widely studied in the past in various species and are key factors determining the formation and evolution of allopolyploids (Yoo et al., 2014). For example, using microarrays, it was found that the total expression of duplicated genes in allopolyploid cotton may be similar to that of one of its parents, or may be lower or higher than those of both parents (Rapp et al., 2009; Flagel and Wendel, 2010). The first scenario was then termed as expression level dominance, which was also found in many other species (Grover et al., 2012). By directly comparing the contribution of homeologs of two parents with the transcripts pool of allopolyploids, it was found that duplicate gene pairs may express unequally, and display homeolog expression bias that usually varies among tissues and was associated with progenitor diploid expression levels (Flagel et al., 2008; Chaudhary et al., 2009; Buggs et al., 2011; Dong and Adams, 2011; Combes et al., 2013). In addition to the genetic and epigenetic changes, nonadditive gene expression could also result from many other underlying causal factors, such as maternal–paternal influence (Gong et al., 2012), gene dosage balance (Birchler and Veitia, 2010), and cis- and trans regulation (Shi et al., 2012).
However, the knowledge of genetic remodeling and gene expression changes of polyploid formation is largely derived from the studies of resynthesized polyploids that mimic natural counterparts, obtained from hybridization between the limited genotypes of extant presumed diploids in several generations. These synthetics or neo-allopolyploids commonly display phenotypic and genetic instability, low fertility, and poor seed yields. These genetic changes and resulting genome instability were possibly common during the formation of existing polyploids, which were eliminated during thousands of years of natural and artificial selection favoring stability and fertility. As natural polyploidy species formed over widely different timeframes (thousands to millions of years), except those in Tragopogon (Ownbey, 1950), Spartina (Ainouche et al., 2004), and Senecio (Ashton and Abbott, 1992) within the last 150 years, it is generally impossible to find their precise diploid progenitors and to accurately determine the genetic changes and interactions of parental genomes during the evolutionary process. The extant presumed diploids most closely related to the original diploid progenitors also evolved in parallel with the polyploids under the same or different environments. Also, recurrent or multiple origins for a significant proportion of polyploids (Soltis and Soltis, 2009; Tate et al., 2009b) further complicates the precise detection of genetic changes for each parental genome in natural or domesticated polyploidy species. However, the extraction of the AABB component from bread wheat (AABBDD) and the investigation of the modifications in its phenotype, karyotype, and gene expression provide novel clues for exploring the evolutionary history events after wheat formation (Kerber, 1964; Zhang et al., 2014).
B. napus (genomes AACC, 2n = 4× = 38), canola or oilseed rape in the mustard family, is the most important edible oil crop. This young species (emerging ~7,500 years ago), without known wild populations, most likely originated from multiple spontaneous and independent interspecific hybridizations between turnip rape (B. rapa; AA, 2n = 2× = 20) and cabbage or kale (B. oleracea; CC, 2n = 2× = 18) growing side by side in southern Europe or the Mediterranean region (Prakash and Hinata, 1980; Kimber and McGregor, 1995; Pratap et al., 2007). In contrast to the resynthesized B. napus, which shows genetic and epigenetic changes, extensive chromosomal rearrangements, and changes in gene expression and phenotype (Song et al., 1995; Pires et al., 2004; Lukens et al., 2006; Gaeta et al., 2007; Cifuentes et al., 2010; Gaeta and Chris Pires, 2010; Szadkowski et al., 2010; Xiong et al., 2011), comparative mapping and the whole gene sequence revealed that naturally domesticated B. napus is essentially the sum of their diploid progenitors with few genomic rearrangements (Parkin et al., 1995; Axelsson et al., 2000; Rana et al., 2004; Cheung et al., 2009; Wang et al., 2011; Chalhoub et al., 2014; Liu et al., 2014; Parkin et al., 2014). However, without access to exact parental lines, it is still unclear what changes occurred in the expression of parental genes, what factors cause these changes, and what roles these changes play in the evolution and domestication of B. napus.
In previous studies, we first isolated a B. rapa type plant from the B. napus “Oro” variety (Tu et al., 2010). The process of isolation included wide hybridization between B. napus and Iastis indigotica followed by continuous backcrosses between the hybrid with AAC chromosome complements (2n = 29) and I. indigotica. This process might include subsequent genetic events during hybrid embryo development and the growth of young plants, such as (1) the rapid elimination of paternal chromosomes; (2) chromosome doubling of the hybrid genome; and, (3) genome-specific retention or loss of maternal chromosomes. Despite the exact process, C genome chromosomes were lost step-wise during the backcrosses. Using a similar method, another restituted B. rapa (designated as RBR, AeAe, 2n = 20) containing the An sub-genome of B. napus ZS11 (AnAnCnCn, 2n = 38) was also developed (Zhu et al., 2016). Step-wise cytological investigations have revealed the complete retained A genome chromosome and the lost C genome chromosome (Zhu et al., 2016). RBR manifested many significant phenotypes at multiple growth and developmental stages that differed from those of their parental B. napus lines, such as small leaf lobes, short siliques, and greater susceptibility to diseases during and after flowering in humid and high-temperature environments. Particularly, as with many other B. rapa lines, RBR showed severe self-incompatibility that resulted in poor seed set by self-pollination, although with a high pollen stainability (Zhu et al., 2016). Because its parent ZS11 is self-compatible, this result indicated that the self-incompatibility of RBR might result from the loss of C genome interaction when the A genome was extracted. The availability RBR and representative natural B. rapa and B. oleracea provide a unique opportunity to investigate the gene expression changes of the A genome during the extraction process or those that occurred over time at the allotetraploid level. Here, we detail these changes and provide a potential mechanism.
Materials and Methods
Plant Materials and Nomenclature of Different Genomes
A B. rapa line (designated as RBR, restituted B. rapa) with an “extracted” genome (AeAe) of the allotetraploid B. napus (cultivar Zhongshuang 11; designated as ZS11) was produced previously (Zhu et al., 2016). The A and C subgenome within nucleus of B. napus was named as An and Cn, respectively. RBR was propagated via self-pollination for one more generation after its extraction and the second generation was used in this study. B. napus ZS11, four varieties of B. rapa, one belonging to var. pekinensis (cultivar Chiifu-401) and three belonging to var. campestris (TRA1, RA2 and BY1) and three varieties of B. oleracea belonging to var. capitata (Ganlan cv02-12), var. alboglabra (Cjielan) and wild type (TRC) were also used in this study (Supplementary Table 1). The A genome within nucleus of four natural B. rapa was named as Ar. The relationships of all materials used in this study were descripted in Supplementary Figure 1. All plants were grown in a greenhouse under controlled conditions at a 22/20°C day/night cycle, with 16-h day length. For DNA and RNA extraction, young leaves were collected at four-leaf development stages.
RNA Extraction and Transcriptome Analysis
Total RNA was extracted from the young leaves of ZS11, RBR, four B. rapa varieties and three varieties of B. oleracea using Trizol reagent (Invitrogen, CA, USA) and purified using RNeasy Mini Spin Columns (Qiagen, Germany). For each variety used in this study, RNA was extracted independently at least from two plants (Supplementary Table 1). The integrity of RNA was checked with the Agilent Bioanalyzer 2100 Eukaryote Total RNA Nano Series II system. Individual mRNA-Seq libraries (Strand-Specific) with insert sizes of 300 bp for each sample were constructed and sequenced with paired-end reads of 150 bp on the Illumina Hiseq 4000 platform. FastQC (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/) and the NGSQC Toolkit (v2.3) (Patel and Jain, 2012) was used to check, visualize, and control the quality of RNA-seq reads.
To make comparisons of gene expression across species, all clean mRNA-Seq reads were then mapped to the B. napus ‘Darmor-Bzh' reference genome using BWA (Li and Durbin, 2009) with default parameters (Version: 0.7.8-r455), and the uniquely mapped reads were used for further analysis. Differentially expressed genes or homoeologous gene pairs were analyzed using the DESeq package (Anders and Huber, 2010). The raw read counts for each gene were calculated using HTSeq-count (v0.6.1) (Anders et al., 2015) under union mode and were normalized by respective size factors, which were estimated with the function estimateSizeFactors in the DESeq package. Two-fold changes and P values <0.05 were characterized as DEGs according to a negative binomial model of DESeq and multi-factor analyses relying on the generalized linear model (GLM) formulation.
DNA Methylation Analysis Based on Bisulfite-Treated DNA Sequencing
Qualified genomic DNA of ZS11 and RBR, with two biological duplicates (Supplementary Table 1), along with 0.5% w/w of unmethylated lambda DNA (Promega, Madison, WI, US) included to assess bisulfite conversion efficiency, was sheared using the Covaris S220 sonicator to produce an average fragment size of 200–300 bp. Library construction was performed using the Illumina Paired-End Sample Prep Kit and QIAGEN PCR purification kits according to the Illumina whole-genome bisulfite sequencing for methylation analysis protocol. Methylated adapters from the Illumina TruSeq DNA Sample Prep Kit v2 were substituted for the unmethylated adapters within the Illumina Paired-End Sample Prep Kit. Ligation products equivalent to 450 bp were purified using the BluePippin DNA size selection system (Sage Science). Bisulfite conversion of this DNA was performed using the EZ DNA Methylation Gold Kit (Zymo Research). The libraries were quantified using a BioAnalyzer (Agilent) and a KAPA library quantification kit for Illumina (KAPA Biosystems) and sequenced from both ends (paired-end) for 100 cycles on either a 2,000 or 2,500 platform, and 125 bp paired-end reads were generated.
Trimmomatic (version 0.32) (Bolger et al., 2014) was used for filtering adapters and low-quality reads for raw bisulfite-treated DNA sequencing data. Bismark software (Krueger and Andrews, 2011) was used to align the clean reads to the reference genome using the default parameters. The methylation level for each C cytosine was calculated using the Bismark methylation extractor based on uniquely mapped reads. Fisher's exact test was used for the identification of differentially methylated regions (DMRs) with a false discovery rate (FDR) < 0.01 in 200 bp windows of the B. napus genome.
Gene Ontology (Go) Enrichment Analysis
GO annotations were performed using agriGO (http://bioinfo.cau.edu.cn/agriGO/index.php). The Singular Enrichment Analysis tool was used to do the GO annotations and significant GO term enrichment analysis (Du et al., 2010), which computed GO term enrichment in one set of genes by comparing it with another set, named the target and reference list, respectively. Enrichment was calculated by Fisher's exact test with Hochberg's multitest adjustment (FDR, P <0.05).
Results
Transcriptome Changes in the An Subgenome of B. napus Reflect Effects of Allohexaploidization/Domestication
Because we do not know the exact original progenitor B. rapa of ZS11, we compared the transcriptome of An-subgenome and RBR with those of natural B. rapa (ArAr, 2n = 20), which represented by four cultivars of two subspecies, pekinensis (Chiifu-401) and campestris (TRA1,TRA2, BY1) using RNA-seq (Supplementary Figure 1). We found 15,614(47.3%), 11,642(36.9%), 13,574(42.7%) and 6,218(19.3%) genes to be differentially expressed in An-subgenome compared to those in Chiifu-401, TRA1, TRA2 and BY1, respectively (Supplementary Table 2). Similarly, there were 14,184 (43.8%), 10,366 (33.7%), 13,034 (43.9%), and 3,142 (10.0%) DEGs between RBR and each B. rapa, respectively (Supplementary Table 3). Obviously, transcriptome divergence between An-subgenome of ZS11 and BY1 is the least and the divergence between An-subgenome and Chiifu-401 is the greatest. The same is true with RBR. These results were consistent with the known breeding history of ZS11 and the evolution and/or domestication histories of these B. rapa lines: (1) BY1 has been used as one of the cross parents during the breeding of ZS11 (Supplementary Figure 1); (2) Chiifu-401 is domesticated for use as a vegetable, whereas the other B. rapa are used as oilseeds.
Next, we analyzed the proportions of upregulated versus downregulated genes in paired comparisons. When comparing An-subgenome with each B. rapa, BY1, TRA1, TRA2 and Chiifu-401, 3,281(52.8%), 6,139(52.7%), 7,259(53.5%) and 8,803(56.4%) DEGs, respectively, were found to be upregulated, whereas 2,937(47.2%), 5,503(47.3%), 6,315(46.5%) and 6,811(43.6%) DEGs, respectively, were found to be downregulated (Supplementary Table 2; Figure 1A). When comparing RBR with each B. rapa, 1,573 (50.1%), 5,542 (53.5%), 6,964 (53.4%), and 7,997 (56.4%) DEGs, respectively, were found to be upregulated, whereas 1,569 (49.9%), 4,824 (46.5%), 6,070 (46.6%), and 6,187(43.6%) DEGs, respectively, were found to be downregulated (Supplementary Table 3; Figure 1B). For each comparison, the numbers of upregulated genes were significantly higher than the numbers of downregulated genes in An-subgenome and RBR, except for that in the comparison between RBR and BY1, according to a binomial test (Supplementary Tables 2 and 3). These results indicated that the gene expression in An-subgenome of B. napus and Ae of RBR changed towards a generally higher level than those in Ar genome of natural B. rapa. Notably, this trend was most pronounced between An-subgenome/RBR and Chiifu-401, indicating the combined effects of allopolyploidization and selection under domestication.
Figure 1 Histograms of the up- and downregulated gene numbers in each of the four pairwise comparisons between An-subgenome (A) and RBR (B) and four varieties of B. rapa: Chiifu-401, BY1, TRA1, and TRA2. An represented An-subgenome of B. napus ZS11 and Ae represented RBR.
To test this speculation further, we analyzed shared DEGs between An-subgenome and four natural B. rapa as well as between RBR and four natural B. rapa. We found that 1,247 common genes were upregulated and 740 genes were downregulated in the An-subgenome of ZS11 (Supplementary Table 2; Figure 1A). We also found that 507 common genes were upregulated and 242 genes were downregulated in RBR (Supplementary Table 3; Figure 1B). In both cases, the numbers of upregulated genes were significantly higher than those of downregulated genes. These changes may be considered as the results of the effects of allopolyploidization/domestication occurring since the origin of B. napus. Together, these transcriptome comparisons suggest that allopolyploidization and domestication at the tetraploid level might have significantly increased the expression level of a substantial proportion of genes in the An-subgenome of B. napus relative to its progenitor B. rapa.
Reversible Gene Expression Changes During the Allo-/Deallopolyploidization of B. napus
Using four natural B. rapa representing the original diploid progenitor of B. napus, together with RBR, we were able to analyze the reversible gene expression changes occurring during the allo-/deallopolyploidization process in B. napus. We compared the transcriptome of RBR with that of B. napus An subgenome by RNA-seq analysis with three biological replicates. We found that 5,995 genes (18.8% of a total of 31,867 genes) were differentially expressed between Ae of RBR and An subgenome of ZS11. Among these genes, 3,273 (53.7%) were downregulated and 2,722 genes were upregulated in RBR. The numbers of down regulated genes were significantly higher than that of the upregulated genes (Supplementary Table 3). In order to investigate whether these changes are in fact indicative of the recovery of antecedent changes occurring during the formation and evolution of B. napus, we analyzed the expression of shared DEGs between An-subgenome and four natural B. rapa in RBR. We found that 837 (67.1%) of 1,247 upregulated genes and 475 (64.3%) of 740 downregulated genes in An subgenome of ZS11 were down- and up-regulated in RBR, respectively. These genes represented 25.6% of the total 3,273 downregulated genes and 17.4% of the total 2722 upregulated genes observed during B. rapa extraction (Figure 2). Obviously, by the unique restitution of B. rapa RBR, expression changes of a set of genes (21.5%) were revealed to be reversible during the allo-/deallopolyploidization of B. napus. In other words, these genes are highly or lowly expressed specifically in the A genome at the tetraploid level in B. napus. So, it is reasonable that the pattern of gene expression of the A genome in RBR was changed rapidly towards that in the natural B. rapa after extraction from the nucleus of B. napus.
Figure 2 Diagram of the total number of genes specially up- and downregulated in the A subgenome of ZS11 (An) in comparison with RBR (Ae) and four natural B. rapa lines (Ar). The arrows indicate the genomes in which genes are down-regulated, and the numbers inside represent the number of downregulated genes. Red and green wavy lines indicate a continuous trend of gene expression from Ar to An and Ae. The numbers represent the numbers of genes with this expression trend, and the percentages represent the proportions of these genes in the total up- and downregulated genes in each pairwise comparison. Ar vs An represented DEGs shared by all B. rapa varieties versus An-subgenome of B. napus. Ar vs Ae represented DEGs shared by all four B. rapa versus RBR. An vs Ae represented DEGs between RBR and those in An-subgenome of B. napus ZS11.
B. napus-Specific Upregulated Genes Are Enriched in Distinct Gene Ontology Categories
To explore whether those genes with reversible expression patterns in the allo-/deallopolyploidization are related to specific functional features, we performed a gene ontology (GO) analysis. We found that both down- and upregulated genes in RBR versus those in the An subgenome of B. napus showed over-representation of several categories of genes involved in metabolism and responses to stress stimulus (Supplementary Datasheet 2; Figure 3A, B). Particularly, B. napus-specific upregulated genes, either compared with natural B. rapa or with RBR, showed over-representation for two categories: heterocyclic metabolic processes and responses to abiotic stimuli (Supplementary Datasheet 1; Figure 3C). In order to elucidate the origins of these changes of An-subgenome of B. napus, we compared the homoeologous gene (Supplementary Datasheet 1) expression in two parental species of B. napus: B. oleracea and B. rapa, which were represented by four varieties of B. rapa and three varieties of B. oleracea. In total, 12 pairs of comparisons were made (Supplementary Table 4). Among the DEGs in each comparison, 455 and 518 upregulated genes were shared by all B. rapa and all B. oleracea lines, respectively. GO enrichment analysis found that 455 highly expressed genes in B. rapa showed several enriched GO categories involved in metabolic processes, whereas 518 highly expressed genes in B. oleracea lines showed several enriched GO categories involved in stress responses (Supplementary Datasheet 2; Figure 3D, E). These results indicate that there are differences in capabilities of metabolism and disease resistance between natural B. rapa and B. oleracea. An subgenome-specific up-/downregulated genes shown enrichment in metabolism and responses to stress might result from trans regulation, either positive or negative, by the Cn subgenome.
Figure 3 Gene ontology (GO) enrichment (> 50%) for genes that are differentially expressed before and after RBR extraction, specially up- and downregulated in B. napus as well as between natural B. rapa and B. oleracea. The x-axis is the GO category terms, and the y-axis is the percentages of genes mapped by the GO category terms, which are calculated by the number of genes mapped to a given GO category divided by the number of all mapped genes (Vermilion bars denote the percentages of each GO category in all expressed genes, 26,539 in total). (A) Gene ontology analysis of the 2722 upregulated genes in RBR comparing to those in An-subgenome of B. napus. (B) Gene ontology analysis of the 3,273 downregulated genes in RBR comparing to those in An-subgenome of B. napus. (C) Gene ontology analysis of the upregulated 837 genes in An-subgenome of B. napus comparing to those in natural B. rapa as well as in RBR. (D) Gene ontology analysis of the upregulated 455 homoeologous genes in B. rapa comparing to those in B. oleracea. (E) Gene ontology analysis of the 518 upregulated homoeologous genes in B. oleracea comparing to those in B. rapa. Note: only the Ontology of biological process was shown here for each enrichment.
Immediate and Dramatic Gene Expression Changes of An Subgenome During Its Extraction Are Associated With An and Cn Homoeologous Expression Bias
Studies in B. napus lines showed that no pronounced genome dominance was present, but there are significantly biased contributions to the gene expression of the An and Cn homeologs in different tissues (Chalhoub et al., 2014). We also analyzed the contribution of homeolog An and Cn (Chalhoub et al., 2014) to the transcriptome of B. napus in young leaves used here. ZS11 showed significantly higher expression for 3,323 (12.9%) An homeologs and 3,542 (13.7%) Cn homeologs in 25,699 total gene pairs. There are 18,834 (73.3%) gene pairs in ZS11, which showed equal contributions to the transcriptome. These results are consistent with those found in Darmor-bzh (Chalhoub et al., 2014).
We next analyzed the expression of these differentially expressed gene pairs in RBR and natural B. rapa. A total of 683 (20.5%) and 302 (9.1%) of the 3,323 highly expressed An homeologs in B. napus ZS11 were downregulated in RBR and in all natural B. rapa, respectively (Figure 4A). Moreover, 193 common genes were identified simultaneously in gene sets of An > Cn, An > Ar and An > Ae (Figure 4A). In contrast, only 121, and 25 of 3,323 highly expressed An homeologs were upregulated in RBR and all natural B. rapa, respectively. Only 12 common genes were found simultaneously in above three gene sets (Figure 4B). Similarly, 769 (21.7%) and 290 (8.2%) of 3,542 An homeologs with lower expression were upregulated in RBR and all natural B. rapa, respectively, and 203 common genes were identified simultaneously in gene sets of An < Cn, An < Ar and An < Ae (Figure 4C). However, only 167 and 34 genes were downregulated in RBR and four natural B. rapa, respectively, and 21 common genes were found in gene sets of An < Cn, An < Ar and An < Ae (Figure 4D). Taken together, these results indicated that the An homeologs that showed higher expression than Cn homeologs in B. napus tended to be downregulated in extracted B. rapa and natural B. rapa, while the An homeologs with lower expression tended to be upregulated.
Figure 4 Diagram shows the potential relationships between gene expression changes in pairwise comparisons (Ae vs An and An vs Ar) and homoeologous expression bias of An and Cn. A, B: Specially up- (A) and downregulated (B) genes under An > Cn; C, D: Specially up (C) and down (D) regulated genes under An < Cn.
Alternations in DNA Methylation in AA Components After Its Extraction From B. napus
To investigate how DNA methylation changes occurred in RBR immediately after its extraction from B. napus, which might contribute to the gene expression changes, we performed whole genome bisulfite sequencing analysis of RBR, as well as of ZS11. For RBR and ZS11, 59,473,375 and 106,938,717 raw reads were produced, respectively. After removal of those low-quality data, 57,134,049 and 102,602,342 clean reads were produced, respectively. Using the B. napus genome as a reference, 30,586,180 (53.53%) and 46,804,759 (45.62%) reads were uniquely mapped. In total, 43.21%, 18.34%, and 3.96% of the total CG, CHG, and CHH (where H = A, C or T) loci were methylated in RBR, and 55.50%, 27.90%, and 5.81% of the total loci were methylated in ZS11. The overall methylation at different loci in the later is consistent with those in sequenced B. napus leaves, in which 53% CG, 22% CHG, and 7% CHH were methylated (Chalhoub et al., 2014). To compare methylation levels between RBR and the An-subgenome within B. napus, the DMR (differentially methylated region) for each was identified. In total, there are 1,523 (500 bp-4.5 kb) and 4,654 (500 bp-4.7 kb) regions, which span about 0.98 Mb and 3.07 Mb of the genome, respectively, show hypo- and hyper-methylation in RBR versus ZS11 (Supplementary Datasheet 3). Obviously, RBR showed higher numbers of hyper-methylated regions than that of hypo-methylated regions immediately after its extraction. We found 274 (6.4% of total 3,726) genes were down-regulated in RBR exactly within the hyper-methylated regions, while 149 (4.7% of total 3,204) upregulated genes were located in hypo-methylated regions. In both cases, most of the DEGs within DMR have altered methylation involved in promoter and exon regions (Figure 5). Interestingly, there are many more genes in hyper-methylated regions (207) than in hypo-methylated regions (93), with changed methylation status involved in exon (Figure 5). Overall, more regions were changed toward hyper-methylation in RBR, which is consistent with the prevalence of more downregulated genes than upregulated genes in RBR compared to those in An subgenome of B. napus. GO enrichment analysis revealed that 274 down-regulated genes within hyper-methylated regions in RBR involved biological process mainly in meiotic chromosome organization and in molecular function in calcium ion binding, DNA binding, transferase activity, nucleic acid binding and methyltransferase activity. Meanwhile, 149 upregulated genes located in hypo-methylated regions showed enrichment in biological of disaccharide metabolic process, oligosaccharide metabolic process, hormone biosynthetic and metabolic process, etc. and molecular function of lyase activity, phosphoric ester hydrolase activity and phosphatase activity (Supplementary Datasheet 4).
Figure 5 Histograms of the numbers of genes with alternations in methylation status at different locations. P+E represents methylation changes involved in promoter and exon regions and P+I represent methylation changes involved in promoter and intron regions.
Discussion
Over the past few decades, researchers have revealed a wide range of changes in parental gene expression during the formation and evolution of allopolyploids, using different methods investigating natural and synthetic allopolyploids (Adams et al., 2003; Hegarty et al., 2006; Wang et al., 2006; Chaudhary et al., 2009; Pumphrey et al., 2009; Rapp et al., 2009; Akhunova et al., 2010; Chague et al., 2010; Buggs et al., 2011; Grover et al., 2012; Li et al., 2014). These changes occur either immediately after polyploid formation or over long periods of evolution. Although the molecular mechanisms are still not well understood, genetic and epigenetic changes of parental genomes, cis and trans regulation between different genomes, gene dosage, and stochastic effects are the likely main mediators of these changes (Wendel, 2000; Osborn et al., 2003; Riddle and Birchler, 2003; Adams and Wendel, 2005; Comai, 2005; Chen, 2007; Doyle et al., 2008; Soltis and Soltis, 2009; Jackson and Chen, 2010; Birchler, 2012; Madlung and Wendel, 2013; Yoo et al., 2014). However, gene expression changes of a given subgenome(s) within an allopolyploid species has been studied only in polyploid genome environment, thus obstructing any assessment of transcriptome changes to a specific subgenome(s) without confounding effects (e.g., trans-acting factors) from the other subgenome(s). RBR, which was obtained by a deallopolyploidization process, is thus a fascinating plant for analysis of genomic and/or gene expression changes in the An-subgenome of B. napus accumulated since its origin. Here, without the affects from Cn-subgenome, we found that RBR lines have significantly more upregulated genes than downregulated genes in comparison with natural B. rapa. This result indicated that some gene expression of A genome had been changed directionally and permanently since it merged with the C genome in the allotetraploid genome of B. napus. However, it should be pointed out that the exact founder genotype of B. rapa that donated the An sub-genome of B. napus ZS11 might be different from those we analyzed, an issue largely prevalent in all studies using model extant accessions to represent actual progenitors of naturally formed polyploids (Gottlieb, 2004; Zhang et al., 2014).
Gene expression changes originating from the genetic and heritable epigenetic will be “irreversible,” while changes arising from the other mechanism might be “reversible” regarding the subgenome(s) of a given allopolyploid. RBR which was extracted from B. napus, a neopolyploid, also provide the opportunity to identify genes with reversible expression pattern. This way, we found that 837 and 475 genes were specifically up- and downregulated in B. napus, respectively, suggesting that the changes of these genes are “reversible” and might be caused by trans factors or gene dosage effects, for example by small RNA (Palacios et al., 2019). Because these genes are common to both RBR and the An sub-genome of B. napus in comparison with the four natural B. rapa, they would be underestimated in number if the investigation was performed on the same line before polyploid formation and after the extraction. These results also indicated that genes of Ae genome tend to be expressed at levels as in Ar genome of natural B. rapa after its extraction from the AnAnCnCn genome, despite coevolved with the C genome for thousands years.
Gene expression changes might result from both genetic and epigenetic mechanisms (Madlung and Wendel, 2013). Due to the close relationship between the A and C sub-genomes and the lack of genetic mechanisms suppressing homoeologous pairing (Cui et al., 2012), homoeologous exchanges (HE) is considered a very important mechanism in B. napus evolution and domestication, and has been widely investigated in B. napus in particular (Udall et al., 2005; Szadkowski et al., 2010; Xiong et al., 2011; Chalhoub et al., 2014; Stein et al., 2017; Sun et al., 2017; Hurgobin et al., 2018). Upon analysis of the DEGs in each comparison in this study, the genes involved in HE (Homoeologous exchanges) had been removed; thus, those alternations might stem from small fragment changes such as deletions, insertions, and SNPs. In resynthesized B. napus, initial generations showed extensive changes at the genomic level, including deletion bands and novel bands (Song et al., 1995; Lukens et al., 2006; Gaeta et al., 2007; Cui et al., 2013). These changes cannot be entirely attributed to the homoeologous changes driven by meiosis because these changes are detected at the S0 generation (Lukens et al., 2006; Cui et al., 2013). Accordingly, by AFLP analysis with 57 pairs of primers, it was found that RBR not only lost 34.51% of the B. napus-specific fragments which due to the C genome chromosome lost but also have 13.41% of fragments novel for two parents. Also, changes in DNA methylation are usually associated with changes in gene expression, and heritable epigenetic alterations may contribute (Lukens et al., 2006; Xu et al., 2009; Cui et al., 2013). In bread wheat, permanent silencing of gene homeologs can be caused by altered DNA methylation (Shitsukawa et al., 2007; Hu et al., 2013). Similarly, we found that about 423 (6.1% of total 6930) DEGs between RBR and the An-subgenome exclusively occurred within DMR regions.
Our GO enrichment analysis indicated that highly expressed genes in B. rapa and B. oleracea showed several enriched GO categories involved in metabolic processes and defense responses, respectively (Figure 3). We showed that An subgenome-specific up-/downregulated genes shown enrichment in metabolism and responses to stress might result from trans regulation, either positive or negative, by the Cn subgenome. These results are consistent with the observation that B. oleracea usually shows greater disease resistance than B. rapa but grows more slowly. In China, B. napus has replaced the native B. rapa and B. juncea owing to its higher seed yield and stronger resistance to biotic and abiotic stresses. Empirically, it was considered that the improved resistance of B. napus was largely contributed by Cn subgenome which was origin from B. oleracea. The diminished vigor and strong susceptibility to diseases of RBR also suggested that the functional partitioning of disease resistance was enhanced for the C subgenome but attenuated for the A subgenome in B. napus. Based on the above results, we suggest that during B. napus formation and evolution, genes from B. rapa related to defense responses might have been transiently upregulated, whereas those associated with biosynthetic and metabolic processes were transiently downregulated by trans regulation from the C genome. These types of regulation also led to the differential expression between the An and Cn homeologs formed in B. napus. However, it must be pointed out that the current results are mainly based on transcriptome analysis, and experiments on specific gene expression, regulation, and function analysis is needed in future for further confirmation.
Conclusions
Through unique restitution of B. rapa, we were able to investigate the gene expression changes of the AA component of B. napus at allotetraploid level accumulated since its formation, and revealed reversible gene expression changes that occurred during the process of allo-/deallopolyploidization. Functional analysis of involved genes indicated the potential cooperation of parental genomes, which will shed lights on subgenome-specific contribution to polyploids.
Data Availability Statement
All raw data from this study have been submitted to the NCBI SRA database (http://www.ncbi.nlm.nih.gov/sra/) under the project number PRJNA480641.
Author Contributions
ZL, XG, and SL conceived and designed this study. BZ prepared the materials and plants. BZ, XG, CT, QP, and DZ analyzed the data. CT, XG, ZL, and SL wrote and revised the manuscript. All authors have read and approved the final version of the manuscript.
Funding
This research was supported by National Key Research and Development Program of China (grant 2016YFD0100202, 2018YFE0108000, 2016YFD0100305, 2016YFD0101007, 2016YFD0100602), National Natural Science Foundation of China (grant 31471530, 31801391, 31770250), Natural Science Foundation of Hubei Province (2019CFB628) and the Basic CAAS Fundation for International Cooperation (Y2019GH08).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2019.01279/full#supplementary-material
References
Adams, K. L., Wendel, J. F. (2005). Polyploidy and genome evolution in plants. Curr. Opin. Plant Biol. 8, 135–141. doi: 10.1016/j.pbi.2005.01.001
Adams, K. L., Cronn, R., Percifield, R., Wendel, J. F. (2003). Genes duplicated by polyploidy show unequal contributions to the transcriptome and organ-specific reciprocal silencing. Proc. Natl. Acad. Sci. U. S. A. 100, 4649–4654. doi: 10.1073/pnas.0630618100
Ainouche, M. L., Baumel, A., Salmon, A. (2004). Spartina anglica C. E. Hubbard: a natural model system for analysing early evolutionary changes that affect allopolyploid genomes. Biol. J. Linn. Soc. 82, 475–484. doi: 10.1111/j.1095-8312.2004.00334.x
Akhunova, A. R., Matniyazov, R. T., Liang, H., Akhunov, E. D. (2010). Homoeolog-specific transcriptional bias in allopolyploid wheat. BMC Genomics 11, 505. doi: 10.1186/1471-2164-11-505
Anders, S., Huber, W. (2010). Differential expression analysis for sequence count data. Genome Biol. 11, R106. doi: 10.1186/gb-2010-11-10-r106
Anders, S., Pyl, P. T., Huber, W. (2015). HTSeq–a Python framework to work with high-throughput sequencing data. Bioinformatics 31, 166–169. doi: 10.1093/bioinformatics/btu638
Ashton, P. A., Abbott, R. J. (1992). Multiple origins and genetic diversity in the newly arisen allopolyploid species, Senecio cambrensis Rosser (Compositae). Heredity 68, 25–32. doi: 10.1038/hdy.1992.3
Axelsson, T., Bowman, C. M., Sharpe, A. G., Lydiate, D. J., Lagercrantz, U. (2000). Amphidiploid Brassica juncea contains conserved progenitor genomes. Genome 43, 679–688. doi: 10.1139/g00-026
Beaulieu, J., Jean, M., Belzile, F. (2009). The allotetraploid Arabidopsis thaliana-Arabidopsis lyrata subsp. petraea as an alternative model system for the study of polyploidy in plants. Mol. Genet. Genomics 281, 421–435. doi: 10.1007/s00438-008-0421-7
Birchler, J. A., Veitia, R. A. (2010). The gene balance hypothesis: implications for gene regulation, quantitative traits and evolution. New Phytol. 186, 54–62. doi: 10.1111/j.1469-8137.2009.03087.x
Birchler, J. A. (2012). Insights from paleogenomic and population studies into the consequences of dosage sensitive gene expression in plants. Curr. Opin. In Plant Biol. 15, 544–548. doi: 10.1016/j.pbi.2012.08.005
Bolger, A. M., Lohse, M., Usadel, B. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120. doi: 10.1093/bioinformatics/btu170
Buggs, R. J., Zhang, L., Miles, N., Tate, J. A., Gao, L., Wei, W., et al. (2011). Transcriptomic shock generates evolutionary novelty in a newly formed, natural allopolyploid plant. Curr. Biology: CB 21, 551–556. doi: 10.1016/j.cub.2011.02.016
Chagué, V., Just, J., Mestiri, I., Balzergue, S., Tanguy, A. M., Huneau, C., et al. (2010). Genome-wide gene expression changes in genetically stable synthetic and natural wheat allohexaploids. New Phytol. 187, 1181–1194. doi: 10.1111/j.1469-8137.2010.03339.x
Chalhoub, B., Denoeud, F., Liu, S., Parkin, I. A., Tang, H., Wang, X., et al. (2014). Early allopolyploid evolution in the post-Neolithic Brassica napus oilseed genome. Science 345, 950–953. doi: 10.1126/science.1253435
Chaudhary, B., Flagel, L., Stupar, R. M., Udall, J. A., Verma, N., Springer, N. M., et al. (2009). Reciprocal silencing, transcriptional bias and functional divergence of homeologs in polyploid cotton (Gossypium). Genetics 182, 503–517. doi: 10.1534/genetics.109.102608
Chen, M., Ha, M., Lackey, E., Wang, J., Chen, Z. J. (2008). RNAi of met1 reduces DNA methylation and induces genome-specific changes in gene expression and centromeric small RNA accumulation in Arabidopsis allopolyploids. Genetics 178, 1845–1858. doi: 10.1534/genetics.107.086272
Chen, Z. J. (2007). Genetic and epigenetic mechanisms for gene expression and phenotypic variation in plant polyploids. Annu. Rev. Plant Biol. 58, 377–406. doi: 10.1146/annurev.arplant.58.032806.103835
Chen, Z. J. (2010). Molecular mechanisms of polyploidy and hybrid vigor. Trends Plant Sci. 15, 57–71. doi: 10.1016/j.tplants.2009.12.003
Cheung, F., Trick, M., Drou, N., Lim, Y. P., Park, J. Y., Kwon, S. J., et al. (2009). Comparative analysis between homoeologous genome segments of Brassica napus and its progenitor species reveals extensive sequence-level divergence. Plant Cell 21, 1912–1928. doi: 10.1105/tpc.108.060376
Cifuentes, M., Eber, F., Lucas, M. O., Lode, M., Chevre, A. M., Jenczewski, E. (2010). Repeated polyploidy drove different levels of crossover suppression between homoeologous chromosomes in Brassica napus allohaploids. Plant Cell 22, 2265–2276. doi: 10.1105/tpc.109.072991
Comai, L. (2005). The advantages and disadvantages of being polyploid. Nat. Rev. Genet. 6, 836–846. doi: 10.1038/nrg1711
Combes, M. C., Dereeper, A., Severac, D., Bertrand, B., Lashermes, P. (2013). Contribution of subgenomes to the transcriptome and their intertwined regulation in the allopolyploid Coffea arabica grown at contrasted temperatures. New Phytol. 200, 251–260. doi: 10.1111/nph.12371
Cui, C., Ge, X., Zhou, Y., Li, M., Li, Z. (2012). Cytoplasmic and genomic effects on meiotic pairing in Brassica hybrids and allotetraploids from pair crosses of three cultivated diploids. Genetics 191, 725–738. doi: 10.1534/genetics.112.140780
Cui, C., Ge, X., Zhou, Y., Li, M., Li, Z. (2013). Cytoplasmic and genomic effects on non-meiosis-driven genetic changes in Brassica hybrids and allotetraploids from pairwise crosses of three cultivated diploids. PloS One 8, e65078. doi: 10.1371/journal.pone.0065078
Dong, S., Adams, K. L. (2011). Differential contributions to the transcriptome of duplicated genes in response to abiotic stresses in natural and synthetic polyploids. New Phytol. 190, 1045–1057. doi: 10.1111/j.1469-8137.2011.03650.x
Doyle, J. J., Flagel, L. E., Paterson, A. H., Rapp, R. A., Soltis, D. E., Soltis, P. S., et al. (2008). Evolutionary genetics of genome merger and doubling in plants. Annu. Rev. Genet. 42, 443–461. doi: 10.1146/annurev.genet.42.110807.091524
Du, Z., Zhou, X., Ling, Y., Zhang, Z., Su, Z. (2010). agriGO, a GO analysis toolkit for the agricultural community. Nucleic Acids Res. 38, W64–W70. doi: 10.1093/nar/gkq310
Flagel, L. E., Wendel, J. F. (2010). Evolutionary rate variation, genomic dominance and duplicate gene expression evolution during allotetraploid cotton speciation. New Phytol. 186, 184–193. doi: 10.1111/j.1469-8137.2009.03107.x
Flagel, L., Udall, J., Nettleton, D., Wendel, J. (2008). Duplicate gene expression in allopolyploid Gossypium reveals two temporally distinct phases of expression evolution. BMC Biol. 6, 16. doi: 10.1186/1741-7007-6-16
Gaeta, R. T., Chris Pires, J. (2010). Homoeologous recombination in allopolyploids: the polyploid ratchet. New Phytol. 186, 18–28. doi: 10.1111/j.1469-8137.2009.03089.x
Gaeta, R. T., Pires, J. C., Iniguez-Luy, F., Leon, E., Osborn, T. C. (2007). Genomic changes in resynthesized Brassica napus and their effect on gene expression and phenotype. Plant Cell 19, 3403–3417. doi: 10.1105/tpc.107.054346
Gong, L., Salmon, A., Yoo, M. J., Grupp, K. K., Wang, Z., Paterson, A. H. (2012). The cytonuclear dimension of allopolyploid evolution: an example from cotton using rubisco. Mol. Biol. Evol. 29, 3023–3036. doi: 10.1093/molbev/mss110
Gottlieb, L. D. (2004). Rethinking classic examples of recent speciation in plants. New Phytol. 161, 71–82. doi: 10.1046/j.1469-8137.2003.00922.x
Grover, C. E., Gallagher, J. P., Szadkowski, E. P., Yoo, M. J., Flagel, L. E., Wendel, J. F. (2012). Homoeolog expression bias and expression level dominance in allopolyploids. New Phytol. 196, 966–971. doi: 10.1111/j.1469-8137.2012.04365.x
Han, F., Fedak, G., Guo, W., Liu, B. (2005). Rapid and repeatable elimination of a parental genome-specific DNA repeat (pGc1R-1a) in newly synthesized wheat allopolyploids. Genetics 170, 1239–1245. doi: 10.1534/genetics.104.039263
Hegarty, M. J., Barker, G. L., Wilson, I. D., Abbott, R. J., Edwards, K. J., Hiscock, S. J. (2006). Transcriptome shock after interspecific hybridization in senecio is ameliorated by genome duplication. Curr. Biol. 16, 1652–1659. doi: 10.1016/j.cub.2006.06.071
Hegarty, M. J., Batstone, T., Barker, G. L., Edwards, K. J., Abbott, R. J., Hiscock, S. J. (2011). Nonadditive changes to cytosine methylation as a consequence of hybridization and genome duplication in Senecio (Asteraceae). Mol. Ecol. 20, 105–113. doi: 10.1111/j.1365-294X.2010.04926.x
Hu, Z., Han, Z., Song, N., Chai, L., Yao, Y., Peng, H., et al. (2013). Epigenetic modification contributes to the expression divergence of three TaEXPA1 homoeologs in hexaploid wheat (Triticum aestivum). New Phytologist. 197, 1344–1352. doi: 10.1111/nph.12131
Hurgobin, B., Golicz, A. A., Bayer, P. E., Chan, C. K., Tirnaz, S., Dolatabadian, A., et al. (2018). Homoeologous exchange is a major cause of gene presence/absence variation in the amphidiploid Brassica napus. Plant Biotechnol. J. 16, 1265–1274. doi: 10.1111/pbi.12867
Jackson, S., Chen, Z. J. (2010). Genomic and expression plasticity of polyploidy. Curr. Opin. In Plant Biol. 13, 153–159. doi: 10.1016/j.pbi.2009.11.004
Joly, S., Rauscher, J. T., Sherman-Broyles, S. L., Brown, A. H., Doyle, J. J. (2004). Evolutionary dynamics and preferential expression of homeologous 18S-5.8S-26S nuclear ribosomal genes in natural and artificial glycine allopolyploids. Mol. Biol. Evol. 21, 1409–1421. doi: 10.1093/molbev/msh140
Kashkush, K., Feldman, M., Levy, A. A. (2003). Transcriptional activation of retrotransposons alters the expression of adjacent genes in wheat. Nat. Genet. 33, 102–106. doi: 10.1038/ng1063
Kerber, E. R. (1964). Wheat: reconstitution of the Tetraploid component (AABB) of hexaploids. Science 143, 253–255. doi: 10.1126/science.143.3603.253
Khasdan, V., Yaakov, B., Kraitshtein, Z., Kashkush, K. (2010). Developmental timing of DNA elimination following allopolyploidization in wheat. Genetics 185, 387–390. doi: 10.1534/genetics.110.116178
Kimber, D. S., McGregor, D. I. (1995). “The species and their origin, cultivation and world production,” in Brassica Oilseeds: Production and Utilization. Eds. Kimber, D., McGregor, D. I. (Wallingford, UK19: CABI Publishing).
Kovarik, A., Dadejova, M., Lim, Y. K., Chase, M. W., Clarkson, J. J., Knapp, S., et al. (2008). Evolution of rDNA in Nicotiana allopolyploids: a potential link between rDNA homogenization and epigenetics. Ann. Bot. 101, 815–823. doi: 10.1093/aob/mcn019
Krueger, F., Andrews, S. R. (2011). Bismark: a flexible aligner and methylation caller for Bisulfite-Seq applications. Bioinformatics 27, 1571–1572. doi: 10.1093/bioinformatics/btr167
Książczyk, T., Kovarik, A., Eber, F., Huteau, V., Khaitova, L., Tesarikova, Z., et al. (2011). Immediate unidirectional epigenetic reprogramming of NORs occurs independently of rDNA rearrangements in synthetic and natural forms of a polyploid species Brassica napus. Chromosoma 120, 557–571. doi: 10.1007/s00412-011-0331-z
Lashermes, P., Combes, M. C., Hueber, Y., Severac, D., Dereeper, A. (2014). Genome rearrangements derived from homoeologous recombination following allopolyploidy speciation in coffee. Plant J. 78, 674–685. doi: 10.1111/tpj.12505
Li, H., Durbin, R. (2009). Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754–1760. doi: 10.1093/bioinformatics/btp324
Li, A., Liu, D., Wu, J., Zhao, X., Hao, M., Geng, S., et al. (2014). mRNA and small RNA transcriptomes reveal insights into dynamic homoeolog regulation of allopolyploid heterosis in nascent hexaploid wheat. Plant Cell 26, 1878–1900. doi: 10.1105/tpc.114.124388
Lim, K. Y., Soltis, D. E., Soltis, P. S., Tate, J., Matyasek, R., Srubarova, H., et al. (2008). Rapid chromosome evolution in recently formed polyploids in Tragopogon (Asteraceae). PloS One 3, e3353. doi: 10.1371/journal.pone.0003353
Liu, S., Liu, Y., Yang, X., Tong, C., Edwards, D., Parkin, I. A., et al. (2014). The Brassica oleracea genome reveals the asymmetrical evolution of polyploid genomes. Nat. Commun. 5, 3930. doi: 10.1038/ncomms4930
Lukens, L. N., Pires, J. C., Leon, E., Vogelzang, R., Oslach, L., Osborn, T. (2006). Patterns of sequence loss and cytosine methylation within a population of newly resynthesized Brassica napus allopolyploids. Plant Physiol. 140, 336–348. doi: 10.1104/pp.105.066308
Madlung, A., Wendel, J. F. (2013). Genetic and epigenetic aspects of polyploid evolution in plants. Cytogenetic Genome Res. 140, 270–285. doi: 10.1159/000351430
Madlung, A., Tyagi, A. P., Watson, B., Jiang, H., Kagochi, T., Doerge, R. W., et al. (2005). Genomic changes in synthetic Arabidopsis polyploids. Plant J. 41, 221–230. doi: 10.1111/j.1365-313X.2004.02297.x
Martinez Palacios, P., Jacquemot, M. P., Tapie, M., Rousselet, A., Diop, M., Remoué, C., et al. (2019). Assessing the response of small RNA populations to allopolyploidy using resynthesized Brassica napus Allotetraploids. Mol. Biol. Evol. 36, 709–726. doi: 10.1093/molbev/msz007
Osborn, T. C., Pires, J. C., Birchler, J. A., Auger, D. L., Chen, Z. J., Lee, H. S., et al. (2003). Understanding mechanisms of novel gene expression in polyploids. Trends Genet. 19, 141–147. doi: 10.1016/S0168-9525(03)00015-5
Ownbey, M. (1950). Natural Hybridization and Amphiploidy in the Genus Tragopogon. Am. J. Bot. 37, 487–499. doi: 10.1002/j.1537-2197.1950.tb11033.x
Ozkan, H., Feldman, M. (2009). Rapid cytological diploidization in newly formed allopolyploids of the wheat (Aegilops-Triticum) group. Genome 52, 926–934. doi: 10.1139/G09-067
Parisod, C., Alix, K., Just, J., Petit, M., Sarilar, V., Mhiri, C., et al. (2010). Impact of transposable elements on the organization and function of allopolyploid genomes. New Phytol. 186, 37–45. doi: 10.1111/j.1469-8137.2009.03096.x
Parkin, I. A., Sharpe, A. G., Keith, D. J., Lydiate, D. J. (1995). Identification of the A and C genomes of amphidiploid Brassica napus (oilseed rape). Genome 38, 1122–1131. doi: 10.1139/g95-149
Parkin, I. A., Koh, C., Tang, H., Robinson, S. J., Kagale, S., Clarke, W. E., et al. (2014). Transcriptome and methylome profiling reveals relics of genome dominance in the mesopolyploid Brassica oleracea. Genome Biol. 15, R77. doi: 10.1186/gb-2014-15-6-r77
Patel, R. K., Jain, M. (2012). NGS QC Toolkit: a toolkit for quality control of next generation sequencing data. PloS One 7, e30619. doi: 10.1371/journal.pone.0030619
Pires, J. C., Zhao, J., Schranz, M. E., Leon, E. J., Osborn, T. C., et al. (2004). Flowering time divergence and genomic rearrangements in resynthesized Brassica polyploids (Brassicaceae). Biol. J. Linn. Soc. 82, 675–688. doi: 10.1111/j.1095-8312.2004.00350.x
Pontes, O., Neves, N., Silva, M., Lewis, M. S., Madlung, A., Comai, L., et al. (2004). Chromosomal locus rearrangements are a rapid response to formation of the allotetraploid Arabidopsis suecica genome. Proc. Natl. Acad. Sci. U. S. A. 101, 18240–18245. doi: 10.1073/pnas.0407258102
Prakash, S., Hinata, K. (1980). Taxonomy, cytogenetics and origin of crop Brassicas, a review. Opera Bot. 55, 1–57.
Pratap, A., Kumar, J., Pratap, A., Kumar, J. (2007). History, origin and evolution. Adv. In Botanical Res. 45, 1–20. doi: 10.1016/S0065-2296(07)45001-7
Pumphrey, M., Bai, J., Laudencia-Chingcuanco, D., Anderson, O., Gill, B. S. (2009). Nonadditive expression of homoeologous genes is established upon polyploidization in hexaploid wheat. Genetics 181, 1147–1157. doi: 10.1534/genetics.108.096941
Rana, D., van den Boogaart, T., O’Neill, C. M., Hynes, L., Bent, E., Macpherson, L., et al. (2004). Conservation of the microstructure of genome segments in Brassica napus and its diploid relatives. Plant J. 40, 725–733. doi: 10.1111/j.1365-313X.2004.02244.x
Rapp, R. A., Udall, J. A., Wendel, J. F. (2009). Genomic expression dominance in allopolyploids. BMC Biol. 7, 18. doi: 10.1186/1741-7007-7-18
Riddle, N. C., Birchler, J. A. (2003). Effects of reunited diverged regulatory hierarchies in allopolyploids and species hybrids. Trends In Genet. 19, 597–600. doi: 10.1016/j.tig.2003.09.005
Salmon, A., Flagel, L., Ying, B., Udall, J. A., Wendel, J. F. (2010). Homoeologous nonreciprocal recombination in polyploid cotton. New Phytol. 186, 123–134. doi: 10.1111/j.1469-8137.2009.03093.x
Shaked, H., Kashkush, K., Ozkan, H., Feldman, M., Levy, A. A. (2001). Sequence elimination and cytosine methylation are rapid and reproducible responses of the genome to wide hybridization and allopolyploidy in wheat. Plant Cell 13, 1749–1759. doi: 10.1105/TPC.010083
Shi, X., Ng, D. W., Zhang, C., Comai, L., Ye, W., Chen, Z. J. (2012). Cis- and trans-regulatory divergence between progenitor species determines gene-expression novelty in Arabidopsis allopolyploids. Nat. Commun. 3, 950. doi: 10.1038/ncomms1954
Shitsukawa, N., Tahira, C., Kassai, K., Hirabayashi, C., Shimizu, T., Takumi, S., et al. (2007). Genetic and epigenetic alteration among three homoeologous genes of a class E MADS box gene in hexaploid wheat. Plant Cell 19, 1723–1737. doi: 10.1105/tpc.107.051813
Soltis, P. S., Soltis, D. E. (2009). The role of hybridization in plant speciation. Annu. Rev. Plant Biol. 60, 561–588. doi: 10.1146/annurev.arplant.043008.092039
Song, K., Lu, P., Tang, K., Osborn, T. C. (1995). Rapid genome change in synthetic polyploids of Brassica and its implications for polyploid evolution. Proc. Natl. Acad. Sci. U. S. A. 92, 7719–7723. doi: 10.1073/pnas.92.177719
Stein, A., Coriton, O., Rousseau-Gueutin, M., Samans, B., Schiessl, S. V., Obermeier, C. 1., et al. (2017). Mapping of homoeologous chromosome exchanges influencing quantitative trait variation in Brassica napus. Plant Biotechnol. J. 15, 1478–1489. doi: 10.1111/pbi.12732
Sun, F., Fan, G., Hu, Q., Zhou, Y., Guan, M., Tong, C., et al. (2017). The high-quality genome of Brassica napus cultivar ‘ZS11' reveals the introgression history in semi-winter morphotype. Plant J. 92, 452–468. doi: 10.1111/tpj.13669
Szadkowski, E., Eber, F., Huteau, V., Lodé, M., Huneau, C., Belcram, H., et al. (2010). The first meiosis of resynthesized Brassica napus, a genome blender. New Phytol. 186, 102–112. doi: 10.1111/j.1469-8137.2010.03182.x
Tate, J. A., Joshi, P., Soltis, K. A., Soltis, P. S., Soltis, D. E. (2009a). On the road to diploidization? Homoeolog loss in independently formed populations of the allopolyploid Tragopogon miscellus (Asteraceae). BMC Plant Biol. 9, 80. doi: 10.1186/1471-2229-9-80
Tate, J. A., Symonds, V. V., Doust, A. N., Buggs, R. J., Mavrodiev, E., Majure, L. C., et al. (2009b). Synthetic polyploids of Tragopogon miscellus and T. mirus (Asteraceae): 60 Years after Ownbey's discovery. Am. J. Bot. 96, 979–988. doi: 10.3732/ajb.0800299
Tu, Y. Q., Sun, J., Ge, X. H., Li, Z. Y. (2010). Production and genetic analysis of partial hybrids from intertribal sexual crosses between Brassica napus and isatis indigotica and progenies. Genome 53, 1–11. doi: 10.1139/G09-093
Udall, J. A., Quijada, P. A., Osborn, T. C. (2005). Detection of chromosomal rearrangements derived from homologous recombination in four mapping populations of Brassica napus L. Genetics 169, 967–979. doi: 10.1534/genetics.104.033209
Verhoeven, K. J., Jansen, J. J., van Dijk, P. J., Biere, A. (2010). Stress-induced DNA methylation changes and their heritability in asexual dandelions. New Phytol. 185, 1108–1118. doi: 10.1111/j.1469-8137.2009.03121.x
Wang, J., Tian, L., Lee, H. S., Wei, N. E., Jiang, H., Watson, B., et al. (2006). Genomewide nonadditive gene regulation in Arabidopsis allotetraploids. Genetics 172, 507–517. doi: 10.1534/genetics.105.047894
Wang, X., Wang, H., Wang, J., Sun, R., Wu, J., Liu, S., et al. (2011). The genome of the mesopolyploid crop species Brassica rapa. Nat. Genet. 43, 1035–1039. doi: 10.1038/ng.919
Wendel, J. F. (2000). Genome evolution in polyploids. Plant Mol. Biol. 42, 225–249. doi: 10.1023/A:1006392424384
Xiong, Z., Pires, J. C. (2011). Karyotype and identification of all homoeologous chromosomes of allopolyploid Brassica napus and its diploid progenitors. Genetics 187, 37–49. doi: 10.1534/genetics.110.122473
Xiong, Z., Gaeta, R. T., Pires, J. C. (2011). Homoeologous shuffling and chromosome compensation maintain genome balance in resynthesized allopolyploid Brassica napus. Proc. Natl. Acad. Sci. U. S. A. 108, 7908–7913. doi: 10.1073/pnas.1014138108
Xu, Y., Zhong, L., Wu, X., Fang, X., Wang, J. (2009). Rapid alterations of gene expression and cytosine methylation in newly synthesized Brassica napus allopolyploids. Planta 229, 471–483. doi: 10.1007/s00425-008-0844-8
Yoo, M. J., Liu, X., Pires, J. C., Soltis, P. S., Soltis, D. E. (2014). Nonadditive gene expression in polyploids. Annu. Rev. Genet. 48, 485–517. doi: 10.1146/annurev-genet-120213-092159
Zhang, H., Zhu, B., Qi, B., Gou, X., Dong, Y., Xu, C., et al. (2014). Evolution of the BBAA component of bread wheat during its history at the allohexaploid level. Plant Cell 26, 2761–2776. doi: 10.1105/tpc.114.128439
Zhao, N., Zhu, B., Li, M., Wang, L., Xu, L., Zhang, H., et al. (2011). Extensive and heritable epigenetic remodeling and genetic stability accompany allohexaploidization of wheat. Genetics 188, 499–510. doi: 10.1534/genetics.111.127688
Keywords: Brassica napus, B. rapa, polyploidization, gene expression, DNA methylation
Citation: Pan Q, Zhu B, Zhang D, Tong C, Ge X, Liu S and Li Z (2019) Gene Expression Changes During the Allo-/Deallopolyploidization Process of Brassica napus. Front. Genet. 10:1279. doi: 10.3389/fgene.2019.01279
Received: 12 June 2019; Accepted: 21 November 2019;
Published: 19 December 2019.
Edited by:
Longjiang Fan, Zhejiang University, ChinaReviewed by:
Jinjin Jiang, Yangzhou University, ChinaChangwei Shao, Yellow Sea Fisheries Research Institute (CAFS), China
Copyright © 2019 Pan, Zhu, Zhang, Tong, Ge, Liu and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Bin Zhu emh1Z2cxMzBAMTI2LmNvbQ==; Chaobo Tong dG9uZ2NoYW9ib0AxMjYuY29t; Xianhong Ge Z2V4aWFuaG9uZ0BtYWlsLmh6YXUuZWR1LmNu
†Present address: Bin Zhu, School of Life Science, Guizhou Normal University, Guiyang, China
Dawei Zhang, College of Life Science, Hunan University of Science and Technology, Xiangtan, China