- Research article
- Open Access
Characterization and expression profiles of miRNAs in the triploid hybrids of Brassica napus and Brassica rapa
BMC Genomicsvolume 20, Article number: 649 (2019)
Polyploidy provides a means of interspecific genome transfer to incorporate preferable traits from progenitor to progeny. However, few studies on miRNA expression profiles of interspecific hybrids of B. napus (AnAnCnCn) and B. rapa (ArAr) have been reported.
Here, we apply small RNA sequencing to explore miRNA expression patterns between B. napus, B. rapa and their F1 hybrid. Bioinformatics analysis identified 376, 378, 383 conserved miRNAs and 82, 76, 82 novel miRNAs in B. napus, B. rapa and the F1 hybrid, respectively. Moreover, 213 miRNAs were found to be differentially expressed between B. napus, B. rapa and the F1 hybrid. The present study also shows 211 miRNAs, including 77 upregulated and 134 downregulated miRNAs, to be nonadditively expressed in the F1 hybrid. Furthermore, miRNA synteny analysis revealed high genomic conservation between the genomes of B. napus, B. rapa and their F1 hybrid, with some miRNA loss and gain events in the F1 hybrid.
This study not only provides useful resources for exploring global miRNA expression patterns and genome structure but also facilitates genetic research on the roles of miRNAs in genomic interactions of Brassica allopolyploids.
Polyploidy has been recognized as an important force in plant evolution as well as in novel species formation . Genome duplication increases available genetic materials and results in stronger drought, disease and insect pest resistance, which enables polyploid species to have greater advantages in natural selection. [2, 3]. By comparing the cytology, molecular biology and molecular genetics of polyploids with their parents, very complex genetic variations, including methylation, gene silencing, gene activation and small RNA changes, have usually been found in the former [4, 5].
MicroRNAs (miRNAs) are a class of non-coding single-stranded RNA molecules encoded by endogenous genes with a length of approximately 22 nucleotides [6, 7]. It has been reported that miRNAs are involved in various stages of plant growth and development by silencing target gene(s), and they can be applied for molecular breeding of crop species through direct molecular modulation of specific traits . For example, miRNA controls multiple signal transduction pathways and participates in plant cell differentiation, flower development and seed formation [9,10,11,12,13]. It has been reported that small RNAs act as a genetic buffer in interspecific hybrids and allopolyploids of Arabidopsis . Additionally, multiple small RNAs, including miRNAs, are associated with the molecular mechanism of heterosis in maize [14,15,16,17]. Moreover, previous studies showed that various miRNAs were nonadditively expressed in Brassica hexaploids compared with their parents [18,19,20,21], suggesting a correlation of miRNA expression patterns with many phenotypes in plant allopolyploids. Brassica napus is an important polyploid oil crop and the largest oil-bearing crop in China, accounting for 45% of the total output of edible oil [22,23,24,25,26]. B. rapa belongs to the Brassicaceae family and is widely planted in China [27,28,29,30,31]. Additionally, Brassica species are widely used as model systems to study genomic changes in allopolyploidization. For example, a study of molecular markers and biomass heterosis in the interspecific hybrid between B. napus and B. rapa has been reported [32, 33]. Overall, natural triploid plants have great potential in polyploid breeding, and introgression of genetic resources of B. rapa (AA genome) into B. napus will broaden the genetic basis of the latter and have important application significance [34,35,36]. For instance, interspecific hybridization between B. napus and B. rapa was found to be an efficient and important method of introgressing B. rapa germplasm into B. napus . Triploids from crosses of B. napus and B. rapa are ideal models to assess homologous pairing and recombination occurring between their genomes. Moreover, these triploid hybrids have been used for the production of monosomic addition lines and for gene flow assessment after backcrossing to either parent . The availability of polyploids of B. napus and B. rapa allows us to explore the effects of polyploidization on global gene expression through more accurate comparisons of triploids and their progenitors. Furthermore, the genomes of B. napus and B. rapa have been recently sequenced [39, 40], providing available genome information to globally identify miRNAs and analyze miRNA expression patterns among B. napus, B. rapa and their hybrid.
miRNAs play important roles in the biological and metabolic processes of seeds, including embryogenesis, dormancy and germination . The embryo is a crucial tissue affecting seed development and germination. Although interspecific hybridization between B. napus and B. rapa has been widely investigated , the expression pattern and diversity of embryo miRNAs in triploids remain unclear. The purpose of this study was to investigate the expression profile of miRNAs and their potential targets in triploid hybrids and their progenitors, which may provide vital clues for further and detailed functional studies of embryo miRNAs. In this study, the embryo tissues of triploid hybrid (AnArCn, 21 days after pollination) and control embryo tissues of B. napus and B. rapa were isolated for small RNA sequencing. Conserved and novel miRNAs of B. rapa, B. napus and their hybrid were identified and characterized, and the miRNA expression profiles in the developing embryos of B. napus, B. rapa and their hybrid were investigated. The results showed 376, 378 and 383 conserved miRNAs and 82, 76 and 82 novel miRNAs in the B. napus, B. rapa and F1 hybrid, respectively. Moreover, 213 miRNAs were found to be differentially expressed between B. napus, B. rapa and the F1 hybrid, and 211 miRNAs, including 77 upregulated and 134 downregulated miRNAs, were nonadditively expressed in the F1 hybrid. miRNA synteny analysis revealed some miRNA loss and gain events in the genome of the F1 hybrid. Collectively, this study identified conserved and novel miRNAs and investigated genome-wide miRNA expression profiles between B. napus, B. rapa and their F1 triploid hybrid, and the findings will be very helpful in promoting research on miRNA roles in the genomic interactions of triploids.
Materials and methods
B. napus var. Huashuang 3 (Female parent, AnAnCnCn, 2n = 38) was obtained from Professor Jiangsheng Wu in Huazhong Agricultural University (Professor Jiangsheng Wu bred Huashuang 3 variety). B. rapa var. Tianmen Youcaibai (ArAr, 2n = 18) is a local rapeseed variety in Hubei province. The triploid hybrid (AnArCn, n = 29) was obtained by the hybridization between B. napus var. Huashuang 3 (Female parent, AnAnCnCn, 2n = 38) and B. rapa var. Tianmen Youcaibai (ArAr, 2n = 18). The developing embryos of these plants (21 days after pollination) were harvested and immediately frozen in liquid nitrogen. Total RNAs were isolated using Trizol reagent (Invitrogen, USA). The isolated RNAs were quantified and put in − 80 °C freezer.
Small RNAs sequencing
The total RNA of developing embryos of B. napus, B. rapa and their triploid hybrid (21 days after pollination) was isolated for small RNA sequencing as previously described . The small RNA libraries were sequenced on Illumina HiSeq 2000 platform. The adaptor sequences and low-quality reads were first removed before further bioinformatics analysis. The obtained clean reads were mapped to the non-coding RNA sequences in Rfam database. The mapped reads of snoRNAs, snRNAs, tRNAs, rRNAs and materials containing the poly(A) tail were removed. The filtered sRNAs-seq clean reads were submitted to the NCBI/SRA database with accession number.
Prediction of miRNAs and targets
The conserved and novel miRNAs of B. napus, B. rapa and the triploid hybrid were predicted according to published criteria . The sRNAs-seq clean reads of B. napus, B. rapa and their triploid hybrid were mapped to the genomes of B. napus and B. rapa, respectively. Mireap (http://sourceforge.net/projects/mireap/) was then used to predict miRNA hairpin structures with published parameters. By alignment with miRbase (http://www.mirbase.org/index.shtml), the miRNAs containing ≤2 nucleotide mismatches to known miRNAs were identified to be conserved miRNAs, while the miRNAs containing > 2 nucleotide mismatches to known miRNAs were considered to be novel miRNAs. The miRNA targets of B. napus, B. rapa and their triploid hybrid were predicted respectively using psRobot (http://omicslab.genetics.ac.cn/psRobot/.) according to the published rules [43, 44].
MiRNA differential expression analysis
In order to identify differentially expressed miRNAs between B. napus, B. rapa and their triploid hybrid, the expression level of miRNAs in B. napus or B. rapa was set as a control and the up- or down-regulated miRNAs in triploid hybrid were analyzed by edgeR . MiRNA expression quantification was normalized according to the expression of transcript per million (TPM). First, The TPM of the conserved and novel miRNA in B. napus, B. rapa and their hybrid were first calculated. The fold change values of selected miRNAs were then normalized based on their TPMs. The P-value significance threshold was set by the false discovery rate (FDR). For the comparison of miRNA expression level between B. napus, B. rapa and their hybrid, a fold change of miRNA expression level greater than 2 (with FDR < 0.001 and P value < 0.01) was considered as an indication of significant change and the miRNA was considered to be differentially expressed between B. napus, B. rapa and triploid hybrid. Hierarchical clustering analysis of differentially expressed miRNAs was performed using the Euclidean distance measurement with the R package pheatmap (https://cran.r-project.org/web/packages/pheatmap/index.html).
Quantitative RT-PCR analysis of differentially expressed miRNAs
To validate the differentially expressed miRNAs between B. napus, B. rapa and their triploid hybrid, qRT-PCR was used for expression level analysis. For cDNA synthesis, 1 μg of total RNA from B. napus, B. rapa and their triploid hybrid were reverse-transcribed a RevertAid First Strand cDNA Synthesis Kit (Thermo Scientific). Quantitative RT-PCR reactions were performed using miRNA qRT-PCR kit (#AMPR-0200, GeneCopoeia) with U6 snRNA as the internal control. The comparative CT method was used for the calculation of the relative expression fold changes of miRNAs. All reactions were performed in triplicate with three independent experiments. The primer sequences were listed in Additional file 1: Table S1.
MiRNA synteny analysis
The synteny among B. napus, B. rapa and their triploid hybrid was analyzed as previously described . Brassica napus genome resource (http://www.genoscope.cns.fr/brassicanapus/) and Brassica Database (http:// brassicadb.org/) were used for miRNA synteny analysis. The 10 flanking protein-coding loci of every miRNA were extracted from B. napus and B. rapa genomes. BLASTn was used to perform the homology test of miRNAs and their flanking genes and top-5 hits of every miRNA were selected. The miRNAs containing 1 identical upstream or downstream flanking coding genes were identified as syntenic miRNAs among B. napus, B. rapa and their triploid hybrid. The identified syntenic miRNAs were classified into 4 sets , in which the first 3 sets were used for the construction of Circos map .
Small RNA sequencing of B. napus, B. rapa and their triploid hybrid
Total RNA from developing embryos of B. napus, B. rapa and their F1 hybrid was isolated for the preparation of small RNA libraries. To investigate small RNA expression patterns in B. napus, B. rapa and their F1 hybrid, three biological replicates were sequenced using the Illumina Hi-Seq 2000 platform. After removing adaptor sequences and low-quality reads, ~ 15.2, 16.7 and 15.8 million clean reads were generated for B. napus, B. rapa and the F1 hybrid, respectively. The lengths of most clean reads were 21–24 nt according to statistical analysis (Fig. 1a). The richest small RNA reads were 24 nt and accounted for ~ 40% of the clean reads. Furthermore, we investigated the distribution of reads mapping on the chromosomes of B. napus, B. rapa and F1 hybrid genomes. As shown in Fig. 1b, Circos analysis showed that ~ 7.5 and 10.8 million clean reads mapped to the A and C genomes of B. napus, respectively. Additionally, we observed that ~ 15.7 million clean reads mapped to the A genome of B. rapa. As most of the small RNA reads were 24 nt, we examined the distribution of 24-nt clean reads and found that ~ 3.5 and 4.7 and 5.2 million clean reads of 24 nt mapped to the A and C genomes of B. napus and the A genome of B. rapa, respectively (Fig. 1b). Figure 1c shows that repeat RNA, NAT (Natural Antisense Transcript) RNA and intron sequences were the top 3 enriched types among all mapped small RNAs.
Identification of conserved and novel miRNAs
To identify conserved miRNAs in B. napus, B. rapa and the F1 hybrid, we first investigated correlation of the clean reads with TPM distribution in the B. napus, B. rapa and F1 hybrid samples. The correlation coefficients were 0.929, 0.822 and 0.825 between B. napus and the F1 hybrid, B. rapa and the F1 hybrid and B. napus and B. rapa, respectively (Fig. 2a). Additionally, the TPMs of B. napus, B. rapa and the F1 hybrid were mainly distributed around 1 and 2 (Log10 (TPM + 1)) (Fig. 2b). We next mapped the unique clean reads to known miRBase miRNAs, and 451 conserved miRNAs were identified, including 376 miRNAs in B. napus, 378 miRNAs in B. rapa and 383 miRNAs in the F1 hybrid (Fig. 2c, Additional file 2: Table S2). Further analysis revealed the 451 conserved miRNAs to be distributed in 69 families, including miRNA156, miR159, miRNA156, and miRNA171_1 (Additional file 3: Table S3). We observed that 21, 26 and 29 miRNAs were specifically expressed in B. napus, B. rapa and the F1 hybrid, respectively. In addition, the expression abundance of conserved miRNAs among B. napus, B. rapa and F1 hybrid was quite variable. The expression level of ~ 40% of the conserved miRNAs was less than 10 reads; in contrast, we observed from 10 to 1000 reads for ~ 43% of the conserved miRNAs and more than 1000 reads for ~ 17% of the conserved miRNAs (Additional file 4: Table S4). Novel miRNAs among B. napus, B. rapa and the F1 hybrid were identified using published criteria , with 88 identified (Fig. 2d, Additional file 5: Table S5) including 82, 76 and 82 novel miRNAs in B. napus, B. rapa and the F1 hybrid, respectively. Statistical analysis showed that ~ 77% of 88 novel miRNAs (68 miRNAs) were shared by B. napus, B. rapa and the F1 hybrid. Moreover, 70 novel miRNAs were coexpressed between B. napus and B. rapa, 80 novel miRNAs were coexpressed between B. napus and the F1 hybrid, and 70 novel miRNAs were coexpressed between B. rapa and the F1 hybrid. We also observed that the expression abundance of some novel miRNAs in the F1 hybrid was much higher than the average expression level in its progenitors. For instance, novel_82 had 3-fold higher expression in the F1 hybrid than the average level in its progenitors. In addition, expression abundance analysis of novel miRNAs had an expression level of less than 10 reads for ~ 33% of the novel miRNAs; 10–1000 reads was found for ~ 62% of the novel miRNAs and more than 1000 reads for ~ 5%. To investigate the global miRNA expression profile, we performed hierarchical clustering of expressed miRNAs between the triploid hybrid and its inbred parents, and the results showed 379 miRNAs to be expressed in all three samples (Fig. 2c, d). Euclidean distances based on hierarchical clustering analysis between B. napus and the F1 hybrid, B. rapa and the F1 hybrid, and B. napus and B. rapa were 85.46, 59.78 and 48.66, respectively, which indicates that the expression profile of miRNAs of the F1 hybrid was more similar to that of B. napus than B. rapa. (Fig. 2e).
Differential miRNA expression and miRNA target analysis
To investigate differentially expressed miRNAs between B. napus, B. rapa and the F1 hybrid, we analyzed the expression abundance of each miRNA by normalizing the read count to TPM (Fig. 2e). The results showed that a total of 213 miRNAs were differentially expressed between the B. napus, B. rapa and the F1 hybrid. Further analysis revealed 36 differentially expressed miRNAs (DE-miRNAs) between B. napus and the F1 hybrid, including 13 upregulated and 23 downregulated miRNAs (FDR < 0.001), respectively (Additional file 6: Table S6). Additionally, 155 DE-miRNAs were identified between B. rapa and the F1 hybrid, including 70 upregulated and 85 downregulated miRNAs, respectively (Additional file 6: Table S6). We also observed that 172 miRNAs were differentially expressed between B. napus and B. rapa, including 85 upregulated and 87 downregulated miRNAs, respectively (Additional file 6: Table S6). Collectively, more differentially expressed miRNAs were identified between B. rapa and the F1 hybrid or B. napus and B. rapa. To validate the differentially expressed miRNAs between B. napus, B. rapa and F1 hybrid, some miRNAs were selected for qRT-PCR analysis. Figure 3 showed that qRT-PCR experiment results were consistent with the small RNA sequencing data (Additional file 6: Table S6) and confirmed the differential expression of miRNAs between B. napus, B. rapa and F1 hybrid. The mRNA targets of the identified differentially expressed miRNAs among B. napus, B. rapa and the F1 hybrid were predicted using psRobot, which generated 9196 target genes for 276 miRNAs. All miRNA target genes were assigned Gene Ontology (GO) terms to examine the possible functions of these predicted miRNA target genes. These GO terms are distributed under three main categories: cellular component, biological process and molecular function. “Biological regulation”, “Regulation of cellular process” and “Regulation of biological process” were significantly enriched in the biological process category. “Cell”, “Cell part”, “Intracellular” and “Intracellular part” were significantly enriched in the cellular component category, and “Anion binding”, “Ribonucleotide binding”, “Nucleotide binding” and others were significantly enriched in the molecular function category. Furthermore, the predicted target genes of the miRNAs differentially expressed between the F1 hybrid and its progenitors were annotated using the KEGG database (ranked by P-value). The top 20 KEGG pathways are shown in Fig. 4. The top five pathways included “Cyanoamino acid metabolism” (33 genes, P = 5.68 × 10− 15), “Starch and sucrose metabolism” (42 genes, P = 1.68 × 10− 8), “Phenylpropanoid biosynthesis” (37 unigenes, P = 2.28 × 10− 8), “Selenocompound metabolism” (12 unigenes, P = 2.38 × 10− 6) and “Sulfur metabolism” (12 unigenes, P = 1.5 × 10− 3). As an essential process in plant growth and development, carbohydrate metabolism is closely related to biomass heterosis [48, 49]. We observed that various carbohydrate metabolism pathways, including “Cyanoamino acid metabolism”, “Starch and sucrose metabolism” and “Phenylpropanoid biosynthesis”, were significantly enriched among the target genes of differentially expressed miRNAs between B. napus, B. rapa and the F1 hybrid. Because these differentially expressed miRNAs mainly clustered in miR156, miR159 and miR166, miR167, miR171_1, miR395, miRNA396, and miR399 families, we conclude that these differentially expressed miRNAs might play important roles in biomass heterosis by controlling expression of target genes in the hybrid of B. napus and B. rapa.
Furthermore, we investigated the interaction network between miRNAs and their mRNA targets (Fig. 5). For instance, miR159c and miR156a and miR159a target a large number of genes in miRNA-mRNA interaction networks, suggesting that these three miRNAs might be the key regulators of multiple mRNA targets. As shown in Fig. 5, MYB65 and MYB101, common mRNA targets of miR159c and miR159a, are reported to be essential transcription factors in various plant development processes [50, 51]. Moreover, Fig. 5 indicates that ERS2 is one of the mRNA targets of miR408. A previous study showed that ERS2 is involved in the ethylene-activated signaling pathway and plays many vital roles in plant growth and development . Taken together, these results indicate that the target genes are linked to development and metabolism processes in the F1 hybrid and progenitors, which suggests essential roles of miRNA in allopolyploidization and heterosis.
Non-additive miRNA expression analysis
To investigate how polyploidization affects miRNA expression patterns, we performed non-additive miRNA expression analysis, as previously described . In brief, the miRNA expression values for the F1 hybrid were compared with the average miRNA expression values for B. napus and B. rapa, and miRNAs were identified as being nonadditively expressed if the fold change was greater than 2 (P ≤ 0.05) in comparison with the MPV (mid-parent value), the average miRNA expression value for B. napus and B. rapa. The results showed that 211 miRNAs, including 77 up- and 134 downregulated miRNAs, were nonadditively expressed in the F1 hybrid (Fig. 6a). Importantly, more downregulated miRNAs than upregulated miRNAs (134 vs 77) were nonadditively expressed in the F1 hybrid, which is consistent with previous nonadditive analysis in Brassica hexaploidy . We further analyzed the expression patterns of the 134 downregulated miRNAs in the F1 hybrid and found that 28 miRNAs with higher expression levels in B. napus than in B. rapa were downregulated in the F1 hybrid (Fig. 6b). Moreover, 66 miRNAs with higher expression in B. rapa than in B. napus were downregulated in the F1 hybrid, and 40 miRNAs with higher equal expression between B. napus and B. rapa were downregulated in the F1 hybrid. Taken together, our data indicate that more miRNAs with higher expression in B. rapa than in B. napus are repressed in the F1 hybrid.
MiRNA synteny analysis in allopolyploidization
Brassica species are good models for studying genome expansion and gene fractionation in polyploids [54,55,56,57,58]. The available genome sequences of B. napus and B. rapa are very helpful for investigating miRNA expansion and loss in the hybrid of B. napus and B. rapa. We first performed miRNA synteny analysis between B. napus, B. rapa and the F1 hybrid using the predicted miRNAs and their flanking sequences. As shown in Figs. 7, 95.5% of the B. napus and B. rapa miRNAs (715/749) are present in the genome of the F1 hybrid (Additional file 7: Table S7). Among them, 350 B. rapa miRNAs and 365 B. napus miRNAs mapped to the Ar and AnCn genomes of the F1 hybrid, respectively. We also observed some miRNAs to be located on the opposite genome of the F1 hybrid. For example, 238 B. rapa miRNAs are located on the AnCn genome and 135 B. napus miRNAs on the Ar genome in the F1 hybrid (Additional file 7: Table S7). Overall, the miRNA synteny analysis results indicated high conservation between the AnAnCnCn, ArAr genome and the AnArCn genome of the F1 hybrid. However, we did find that 74 miRNAs, including 68 conserved miRNAs and 6 novel miRNAs, have no homology with the AnArCn genome of the F1 hybrid. Comparative genome analysis between B. napus, B. rapa and their triploid hybrid revealed that some conserved miRNAs in miR166 and miR399 families have been lost from the AnArCn genome of the F1 hybrid, possibly due to gene fractionation and mutation. Furthermore, our results showed that 34 miRNAs, including 30 conserved miRNAs and 4 novel miRNAs (Additional file 7: Table S7), are not present in the syntenic regions of the AnAnCnCn and ArAr genomes, which indicates that these miRNAs may have been generated after the allopolyploidization event in the F1 hybrid. We hypothesize that these newly generated miRNAs were induced by random DNA segmental duplication or insertions in the F1 genome and that they are related to environmental adaptation under specific selection pressure. Taken together, the results indicate a certain amount of miRNA loss and insertions in the genome of the F1 hybrid in comparison with its progenitors B. napus and B. rapa.
Polyploidy is an important process of plant evolution. Interspecific hybridization usually results in hybrid vigor, disease resistance and some other traits in plant allopolyploids [59,60,61]. In the process of polyploid production, genetic and epigenetic effects are stimulated by distant hybridization and polyploidy. Various studies have shown that miRNAs play essential roles in the genomic interaction of many plant hybrids [14, 21, 62]. For instance, miRNA expression will increase, and transposons will be activated in interspecific hybrids, which will alter gene expression in the new polyploid and allow for adaptation to environmental changes [14, 19].
In this study, we used high-throughput sequencing technology to globally identify miRNAs, including conserved and novel miRNAs, in B. napus, B. rapa and their F1 hybrid. A total of 451 conserved miRNAs were identified, including 376, 378 and 383 miRNAs in the B. napus, B. rapa and F1 hybrid, respectively. In addition, 88 novel miRNAs were identified in B. napus, B. rapa and the F1 hybrid. Bioinformatics analysis showed that the expression level of these novel miRNAs was lower than that of conserved miRNAs, which suggested some important roles for these novel miRNAs in B. napus, B. rapa and the F1 hybrid. Moreover, bioinformatics analysis of differentially expressed miRNAs indicated that the genome-wide miRNA expression pattern of the F1 hybrid was more similar to that of B. napus than B. rapa, indicating that miRNAs play important roles in the genomic interaction of plant hybrids.
A total of 9196 target genes of the differentially expressed miRNAs were predicted by psRobot. Gene Ontology analysis showed that the target genes are mainly associated with “biological regulation”, “regulation of cellular process” and “regulation of biological process” in the biological process category. More importantly, we observed that some predicted miRNA target genes are important transcription factors (TFs), including TCP2, TCP4, MYB65 and MYB101 (Fig. 5). TCP transcription factors regulate the growth and development of leaves and flowers . Moreover, miR319 affects the formation of leaf shape and flowering time by specifically inhibiting expression of TCP mRNAs . In addition, a previous study showed that miR159 overexpression inhibits expression of MYB transcription factors and leads to male sterility . Taken together, compared with the roles of miRNAs in other plants, the miRNAs differentially expressed between the F1 hybrid and its progenitors might be potential determinants in the heterosis of the B. napus and B. rapa hybrid.
In polyploidy, nonadditive miRNA expression is usually induced by interspecific hybridization, resulting in gene expression and phenotype complexity. For example, comparative analysis between Arabidopsis thaliana, A. arenosa and their hybrid offspring showed that miRNAs are highly conserved among the species and that a large number of are are nonadditively expressed between the interspecific hybrid and its progenitors . In particular, comparative analysis between a hexaploid hybrid and its parents showed 68.8% of miRNAs to be nonadditively expressed. Moreover, most of the nonadditively expressed miRNAs were inhibited and displayed similar expression patterns in B. rapa and hexaploid hybrids . Our study showed that more downregulated miRNAs (134 miRNAs) than upregulated miRNAs (77 miRNAs) were nonadditively expressed in the F1 hybrid, which is consistent with a previous study investigating hexaploidy in Brassica . Our synteny analysis allowed exact miRNA localization in orthologous regions of the genome in the F1 hybrid, and the results showed that the B. napus and B. rapa genomes and F1 genome share very high homology. Nevertheless, no homology of 74 B. napus and B. rapa miRNAs was found for the F1 hybrid, suggesting that some miRNAs in B. napus and B. rapa have been lost from the F1 hybrid genome, possibly due to DNA fragment deletions and mutations. We also observed that 34 miRNAs in the F1 hybrid were not present in the syntenic regions of the B. napus and B. rapa genomes, indicating that these miRNAs might have been generated after the allopolyploidization event by random DNA segmental duplication or insertions in the F1 genome. Taken together, this study not only provides a useful resource to explore global miRNA expression patterns and evolution and genome structure in allopolyploid plants but also offers a promising perspective on the mechanism of miRNA regulation in allopolyploids, and the findings promote the genetic studies of genomic interaction in allopolyploids.
This study provides fundamental expression patterns of conserved and novel miRNAs between B. napus, B. rapa and F1 hybrid. Bioinformatics analysis showed that 376, 378 and 383 conserved miRNAs, as well as 82, 76 and 82 novel miRNAs, were identified in B. napus, B. rapa and F1 hybrid, respectively. Moreover, 213 miRNAs were found to be differentially expressed between B. napus, B. rapa and F1 hybrid. Meanwhile, the present study showed that a total of 211 miRNAs, including 77 up-regulated and 134 downregulated miRNAs, were non-additively expressed in F1 hybrid. Furthermore, miRNA synteny analysis revealed a high miRNA genomic conservation between the genomes of B. napus, B. rapa and their F1 hybrid genome, and some miRNA loss and gain events in the genome of F1 hybrid. This study provided useful resources to explore global miRNA expression pattern and genome structure, which will be very helpful to reveal the miRNA-mediated regulatory mechanisms in the genomic interactions of Brassica allopolyploid.
Availability of data and materials
The sequencing data of B. napus, B. rapa and F1 hybrid in this study were deposited in NCBI Sequence Read Archive (SRA) Sequence Database with accession number SRR9332479, SRR9332480, SRR9332481, SRR9332482, SRR9332475, SRR9332476, SRR9332477, SRR9332478, SRR9332483.
- A. thaliana :
- B. napus :
- B. rapa :
Differentially expressed miRNA
Transcripts per million
Otto SP, Whitton J. Polyploid incidence and evolution. Annu Rev Genet. 2000;34:401–37.
Ng DWK, Lu J, Chen ZJ. Big roles for small RNAs in polyploidy, hybrid vigor, and hybrid incompatibility. Curr Opin Plant Biol. 2012;15:154–61.
Ni Z, Kim ED, Ha M, Lackey E, Liu J, Zhang Y, Sun Q, Chen ZJ. Altered circadian rhythms regulate growth vigour in hybrids and allopolyploids. Nature. 2009;457:327.
Chen ZJ. Genetic and epigenetic mechanisms for gene expression and phenotypic variation in plant polyploids. Annu Rev Plant Biol. 2007;58:377–406.
Leitch AR, Leitch IJ. Genomic plasticity and the diversity of polyploidy plants. Science. 2008;320:481–3.
Chen XM. Small RNAs and their roles in plant development. Annu Rev Cell Dev Bi. 2009;25:21–44.
Rogers K, Chen X. Biogenesis, turnover, and mode of action of plant microRNAs. Plant Cell. 2013;25:2383–99.
Liu Q, Chen YQ. The potential roles of microRNAs in molecular breeding. Methods Mol Biol. 2012;877:303–11.
Khraiwesh B, Zhu JK, Zhu J. Role of miRNAs and siRNAs in biotic and abiotic stress responses of plants. Biochim Biophys Acta. 2012;1819:137–48.
Sunkar R, Li YF, Jagadeeswaran G. Functions of microRNAs in plant stress responses. Trends Plant Sci. 2012;17:196–203.
Shen D, Suhrkamp I, Wang Y, Liu S, Menkhaus J, Verreet JA, Fan L, Cai D. Identification and characterization of microRNAs in oilseed rape (Brassica napus) responsive to infection with the pathogenic fungus Verticillium longisporum using Brassica AA (Brassica rapa) and CC (Brassica oleracea) as reference genomes. New Phytol. 2014;204:577–94.
Sun G. MicroRNAs and their diverse functions in plants. Plant Mol Biol. 2012;80:17–36.
Li BS, Duan H, Li JG, Deng XW, Yin WL, Xia XL. Global identification of miRNAs and targets in Populus euphratica under salt stress. Plant Mol Biol. 2013;81:525–39.
Ha M, Lu J, Tian L, Ramachandran V, Kasschau KD, Chapman EJ, Carrington JC, Chen X, Wang XJ, Chen ZJ. Small RNAs serve as a genetic buffer against genomic shock in Arabidopsis interspecific hybrids and allopolyploids. Proc Natl Acad Sci U S A. 2009;106:17835–40.
Ding D, Wang Y, Han M, Fu Z, Li W, Liu Z, Hu Y, Tang J. MicroRNA transcriptomic analysis of heterosis during maize seed germination. PLoS One. 2012;7:e39578.
Barber WT, Zhang W, Win H, Varala KK, Dorweiler JE, Hudson ME, Moose SP. Repeat associated small RNAs vary among parents and following hybridization in maize. Proc Natl Acad Sci U S A. 2012;109:10444–9.
Thiebaut F, Rojas CA, Grativol C, Motta MR, Vieira T, Regulski M, Martienssen RA, Farinelli L, Hemerly AS, Ferreira PC. Genome-wide identification of microRNA and siRNA responsive to endophytic beneficial diazotrophic bacteria in maize. BMC Genomics. 2014;15:766.
He G, Zhu X, Elling AA, Chen L, Wang X, Guo L, Liang M, He H, Zhang H, Chen F, Qi Y, Chen R, Deng XW. Global epigenetic and transcriptional trends among two rice subspecies and their reciprocal hybrids. Plant Cell. 2010;22:17–33.
Kenan-Eichler M, Leshkowitz D, Tal L, Noor E, Melamed-Bessudo C, Feldman M, Levy AA. Wheat hybridization and polyploidization results in deregulation of small RNAs. Genetics. 2011;188:263–72.
Zanca AS, Vicentini R, Ortiz-Morea FA, Del Bem LE, da Silva MJ, Vincentz M, Nogueira FT. Identification and expression analysis of microRNAs and targets in the biofuel crop sugarcane. BMC Plant Biol. 2010;10:260.
Shivaprasad PV, Dunn RM, Santos BA, Bassett A, Baulcombe DC. Extraordinary transgressive phenotypes of hybrid tomato are influenced by epigenetics and small silencing RNAs. EMBO J. 2012;31:257–66.
Lu S, Sturtevant D, Aziz M, Jin C, Li Q, Chapman KD, Guo L. Spatial analysis of lipid metabolites and expressed genes reveals tissue-specific heterogeneity of lipid metabolism in high- and low-oil Brassica napus L. seeds. Plant J. 2018;94:915–32.
Bancroft I, Morgan C, Fraser F, Higgins J, Wells R, Clissold L, Baker D, Long Y, Meng J, Wang X, Liu S, Trick M. Dissecting the genome of the polyploid crop oilseed rape by transcriptome sequencing. Nat Biotechnol. 2011;29:762–6.
Xu YH, Zhong L, Wu XM, Fang XP, Wang JB. Rapid alterations of gene expression and cytosine methylation in newly synthesized Brassica napus allopolyploids. Planta. 2009;229:471–83.
Yan X, Dong C, Yu J, Liu W, Jiang C, Liu J, Hu Q, Fang X, Wei W. Transcriptome profile analysis of young floral buds of fertile and sterile plants from the self-pollinated offspring of the hybrid between novel restorer line NR1 and Nsa CMS line in Brassica napus. BMC Genomics. 2013;14:26.
Mason AS, Nelson MN, Castello MC, Yan G, Cowling WA. Genotypic effects on the frequency of homoeologous and homologous recombination in Brassica napus × B. carinata hybrids. Theor Appl Genet. 2011;122:543–53.
Huang X, Lei Y, Guan H, Hao Y, Liu H, Sun G, Chen R, Song S. Transcriptomic analysis of the regulation of stalk development in flowering Chinese cabbage (Brassicacampestris) by RNA sequencing. Sci Rep. 2017;7:15517.
Kakar KU, Nawaz Z, Kakar K, Ali E, Almoneafy AA, Ullah R, Ren XL, Shu QY. Comprehensive genomic analysis of the CNGC gene family in Brassica oleracea: novel insights into synteny, structures, and transcript profiles. BMC Genomics. 2017;18:869.
Chen J, Li R, Xia Y, Bai G, Guo P, Wang Z, Zhang H, Siddique KHM. Development of EST-SSR markers in flowering Chinese cabbage (Brassica campestris L. ssp. chinensis var. utilis Tsen et Lee) based on de novo transcriptomic assemblies. PLoS One. 2017;12:e0184736.
Khalid M, Hassani D, Bilal M, Liao J, Huang D. Elevation of secondary metabolites synthesis in Brassica campestris ssp. chinensis L. via exogenous inoculation of Piriformospora indica with appropriate fertilizer. PLoS One. 2017;12:e0177185.
Huang S, Liu Z, Li C, Yao R, Li D, Hou L, Li X, Liu W, Feng H. Transcriptome Analysis of a Female-sterile Mutant (fsm) in Chinese Cabbage (Brassica campestris ssp. pekinensis). Front Plant Sci. 2017;8:546.
Liu R, Qian W, Meng J. Association of RFLP markers and biomass heterosis in trigenomic hybrids of oilseed rape (Brassicanapus ╳ B. campestris). Theor Appl Genet. 2002;105:1050–7.
Li M, Qian W, Meng J, Li Z. Construction of novel Brassica napus genotypes through chromosomal substitution and elimination using interploid species hybridization. Chromosom Res. 2004;12:417–26.
Qian W, Chen X, Fu D, Zou J, Meng J. Intersubgenomic heterosis in seed yield potential observed in a new type of Brassica napus introgressed with partial Brassica rapa genome. Theor Appl Genet. 2005;110:1187–94.
Qian W, Meng J, Li M, Frauen M, Sass O, Noack J, Jung C. Introgression of genomic components from Chinese Brassica rapa contributes to widening the genetic diversity in rapeseed (B. napus L.), with emphasis on the evolution of Chinese rapeseed. Theor Appl Genet. 2006;113:49–54.
Li M, Li Z, Zhang C, Qian W, Meng J. Reproduction and cytogenetic characterization of interspecific hybrids derived from crosses between Brassica carinata and B. rapa. Theor Appl Genet. 2005;110:1284–9.
Liu RH, Meng JL. RFLP and AFLP analysis of inter- and intraspecific variation of Brassica rapa and B. napus shows that B. rapa is an important genetic resource for B. napus improvement. Yi Chuan Xue Bao. 2006;33:814–23.
Leflon M, Eber F, Letanneur JC, Chelysheva L, Coriton O, Huteau V, Ryder CD, Barker G, Jenczewski E, Chèvre AM. Pairing and recombination at meiosis of Brassica rapa (AA) x Brassica napus (AACC) hybrids. Theor Appl Genet. 2006;113:1467–80.
Chalhoub B, Denoeud F, Liu S, Parkin IA, Tang H, Wang X, et al. Plant genetics. Early allopolyploid evolution in the post-Neolithic Brassica napus oilseed genome. Science. 2014;345:950–3.
Wang X, Wang H, Wang J, Sun R, Wu J, Liu S, et al. The genome of the mesopolyploid crop species Brassica rapa. Nat Genet. 2011;43:1035–9.
Hafner M, Landgraf P, Ludwig J, Rice A, Ojo T, Lin C, Holoch D, Lim C, Tuschl T. Identification of microRNAs and other small regulatory RNAs using cDNA library sequencing. Methods. 2008;44:3–12.
Shen E, Zou J, Hubertus Behrens F, Chen L, Ye C, Dai S, Li R, Ni M, Jiang X, Qiu J, Liu Y, Wang W, Zhu QH, Chalhoub B, Bancroft I, Meng J, Cai D, Fan L. Identification, evolution, and expression partitioning of miRNAs in allopolyploid Brassica napus. J Exp Bot. 2015;66:7241–53.
Meyers BC, Axtell MJ, Bartel B, Bartel DP, Baulcombe D, Bowman JL, Cao X, Carrington JC, Chen X, Green PJ, Griffiths-Jones S, Jacobsen SE, Mallory AC, Martienssen RA, Poethig RS, Qi Y, Vaucheret H, Voinnet O, Watanabe Y, Weigel D, Zhu JK. Criteria for annotation of plant MicroRNAs. Plant Cell. 2008;20:3186–90.
HJ W, Ma YK, Chen T, Wang M, Wang XJ. PsRobot: a web-based plant small RNA meta-analysis toolbox. Nucleic Acids Res. 2012;40:W22–8.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;J26:139–40.
Shen D, Suhrkamp I, Wang Y, Liu S, Menkhaus J, Verreet JA, Fan L, Cai D. 2014. Identification and characterization of microRNAs in oilseed rape (Brassica napus) responsive to infection with the pathogenic fungus Verticillium longisporum using Brassica AA (Brassica rapa) and CC (Brassica oleracea) as reference genomes. New Phytol. 2014;204:577–94.
Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA. Circos: an information aesthetic for comparative genomics. Genome Res. 2009;19:1639–45.
Ko DK, Rohozinski D, Song Q, Taylor SH, Juenger TE, Harmon FG, Chen ZJ. Temporal Shift of Circadian-Mediated Gene Expression and Carbon Fixation Contributes to Biomass Heterosis in Maize Hybrids. PLoS Genet. 2016;12:e1006197.
Chen ZJ. Genomic and epigenetic insights into the molecular bases of heterosis. Nat Rev Genet. 2013;14:471–82.
Xue T, Liu Z, Dai X, Xiang F. Primary root growth in Arabidopsis thaliana is inhibited by the miR159 mediated repression of MYB33, MYB65 and MYB101. Plant Sci. 2017;262:182–9.
Schwab R, Palatnik JF, Riester M, Schommer C, Schmid M, Weigel D. Specific effects of microRNAs on the plant transcriptome. Dev Cell. 2005;8:517–27.
Tanase K, Ichimura K. Expression of ethylene receptors Dl-ERS1-3 and Dl-ERS2, and ethylene response during flower senescence in Delphinium. J Plant Physiol. 2006;163:1159–66.
Shen Y, Zhao Q, Zou J, Wang W, Gao Y, Meng J, Wang J. Characterization and expression patterns of small RNAs in synthesized Brassica hexaploids. Plant Mol Biol. 2014;85:287–99.
Cheng F, Wu J, Fang L, Sun S, Liu B, Lin K, Bonnema G, Wang X. Biased gene fractionation and dominant gene expression among the subgenomes of Brassica rapa. PLoS One. 2012;7:e36442.
Cheng F, Wu J, Cai X, Liang J, Freeling M, Wang X. Gene retention, fractionation and subgenome differences in polyploid plants. Nat Plants. 2018;4:258–68.
Parkin IA, Koh C, Tang H, Robinson SJ, Kagale S, Clarke WE, Town CD, Nixon J, Krishnakumar V, Bidwell SL, Denoeud F, Belcram H, Links MG, Just J, Clarke C, Bender T, Huebert T, Mason AS, Pires JC, Barker G, Moore J, Walley PG, Manoli S, Batley J, Edwards D, Nelson MN, Wang X, Paterson AH, King G, Bancroft I, Chalhoub B, Sharpe AG. Transcriptome and methylome profiling reveals relics of genome dominance in the mesopolyploid Brassica oleracea. Genome Biol. 2014;15:R77.
Sankoff D, Zheng C. Fractionation, rearrangement and subgenome dominance. Bioinformatics. 2012;28:i402–8.
Sankoff D, Zheng C. Whole genome duplication in plants: implications for evolutionary analysis. Methods Mol Biol. 2018;1704:291–315.
Tonosaki K, Michiba K, Bang SW, Kitashiba H, Kaneko Y, Nishio T. Genetic analysis of hybrid seed formation ability of Brassica rapa in intergeneric crossings with Raphanus sativus. Theor Appl Genet. 2013;126:837–46.
Li X, Guo W, Wang B, Li X, Chen H, Wei L, Wang Y, Wu J, Long H. Instability of chromosome number and DNA methylation variation induced by hybridization and amphidiploid formation between Raphanus sativus L. and Brassica alboglabra bailey. BMC Plant Biol. 2010;10:207.
Long H, Chen C, Wang B, Feng Y. rDNA genetic imbalance and nucleolar chromatin restructuring is induced by distant hybridization between Raphanus sativus and Brassica alboglabra. PLoS One. 2015;10:e0117198.
Shen Y, Sun S, Hua S, Shen E, Ye CY, Cai D, Timko MP, Zhu QH, Fan L. Analysis of transcriptional and epigenetic changes in hybrid vigor of allopolyploid Brassica napus uncovers key roles for small RNAs. Plant J. 2017;91:874–93.
Danisman S, van der Wal F, Dhondt S, Waites R, de Folter S, Bimbo A, van Dijk AD, Muino JM, Cutri L, Dornelas MC, Angenent GC, Immink RG. Arabidopsis class I and class II TCP transcription factors regulate jasmonic acid metabolism and leaf development antagonistically. Plant Physiol. 2012;159:1511–23.
Palatnik JF, Wollmann H, Schommer C, Schwab R, Boisbouvier J, Rodriguez R, Warthmann N, Allen E, Dezulian T, Huson D, Carrington JC, Weigel D. Sequence and expression differences underlie functional specialization of Arabidopsis microRNAs miR159 and miR319. Dev Cell. 2007;13:115–25.
Alonso-Peral MM, Li J, Li Y, Allen RS, Schnippenkoetter W, Ohms S, White RG, Millar AA. The microRNA159-regulated GAMYB-like genes inhibit growth and promote programmed cell death in Arabidopsis. Plant Physiol. 2010;154:757–71.
We thank Professor Jiangsheng Wu (Huazhong Agricultural University) for providing B. napus var. Huashuang 3. We also thank lab members of Dr. Maoteng Li (Huazhong University of Science and Technology) for their valuable comments and discussion on of the manuscript. We are also grateful to anonymous reviewers for their valuable suggestions and comments on the manuscript.
This study was supported by the National Key Research and Development Program (2017YFD0101701) and the New Century Talents Support Program by the Ministry of Education of China (NCET110172). The funders had no role in the study design, data analysis and collection, data interpretation and manuscript preparation.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no conflict of interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Primers for qRT-PCR validation of differentially expressed miRNAs between B. napus, B. rapa and the triploid hybrid. (DOCX 15 kb)
Table S2. The identified conserved miRNAs in B. napus, B. rapa and F1. (XLSX 23 kb)
Table S3. The identified miRNA families in B. napus, B. rapa and F1. (XLSX 35 kb)
Table S4. The expression level of the conserved miRNAs in B. napus, B. rapa and F1. (XLSX 23 kb)
Table S5. The identified novel miRNAs in B. napus, B. rapa and F1. (XLSX 13 kb)
Table S6. The differentially expressed miRNAs between B. napus, B. rapa and F1. (XLSX 12 kb)
Table S7. The miRNA synteny analysis of B. napus, B. rapa and F1. (XLSX 30 kb)