BACKGROUND
wide hybridisation and polyploidisation is a common phenomenon in plant evolution that results in the formation of new species  <cit> .
wide hybrids often exhibit more vigorous growth than their parents, and this effect is mainly demonstrated by increases in drought tolerance, pest resistance, flowering time, organ size and biomass, among other factors  <cit> .
furthermore, polyploids show novel traits that are not present in their diploid progenitors  <cit> .
for example, allotetraploid cotton  produces more abundant and higher-quality fibres, and this effect is derived from their aa and dd extant diploid species  <cit> .
however, wide hybrids also exhibit the disadvantage of infertility, as was observed with the wide hybridisation combination of b. rapa and b. oleracea  <cit> .
however, despite this finding, the growth and adaptability advantages of polyploidy have always been a puzzling phenomenon, and the underlying molecular mechanisms are among the most interesting subjects in plant breeding.
allopolyploid formation consists of two processes, namely, hybridisation and chromosome doubling.
hybridisation involves a different genomic combination in the same cells and genomes experience “genomic shock”, whereas chromosome doubling doubles and restructures the genome dosage  <cit> .
recent studies have shown that small rnas, particularly the 24-nt sirnas, play an important role in genome reconstruction and stabilisation  <cit> .
the role of 24-nt sirnas is primarily reflected in two aspects.
the first aspect is the modification of transposons and repetitive sequences for the maintenance of genome stability, which is mediated by rna-dependent dna methylation, and the other aspect is the cis-regulation of gene expression via the transposon gene fragment in the gene near the region   <cit> .
the reduction of sirnas in f <dig> may induce phenotypic and genomic instability in arabidopsis allotetraploids  <cit> .
in interspecific hybrids of arabidopsis, the sirna populations underwent rapid changes in the f <dig> allotetraploid, but stability was maintained in the f <dig> allotetraploid, in which the dna and chromatin were significantly modified  <cit> .
however, in wheat, the number of 24-nt sirna transcripts significantly decreased in the hexaploid compared with that obtained in the parents, and this decrease was accompanied by a decrease in the dna chain cpg island methylation levels  <cit> .
thus, these studies concluded that a decrease in the dna methylation levels may be one of the causes of genomic instability in allopolyploids during the early stage  <cit> .
this phenomenon of changes to the sirnas with an increase in the polyploidy level suggests that sirnas may play a key role in genomic reconstruction and stability after chromosome doubling.
however, little is known regarding how sirnas play this role during wide hybridisation  and chromosome doubling .
mirnas are the primary mediators of the trans-regulation of gene expression  <cit> .
gene silencing mediated by mirnas is an important strategy used at the post-transcriptional level of gene regulation  <cit> .
furthermore, mirnas are conserved in evolution but become active in polyploidisation  <cit> .
changes in the mirna expression levels can affect the expression of the target genes, and this effect is considered to be one of the main causes that results in the phenotypic variation of the polyploidy  <cit> .
in arabidopsis, the number of mirnas in allotetraploids is higher than those observed in their autotetraploid parents  <cit> .
a similar phenomenon was also found in the synthesis of hexaploid wheat  <cit> : the number of mirnas increased with an increase in the ploidy level, and this effect was not associated with the gene dosage balance hypothesis.
this phenomenon showed that other mechanisms, such as cis-epigenetic regulation, may exist.
thus, small rnas, which are a product of non-coding rnas, are involved in regulating gene expression and have become an important factor of gene expression during allopolyploidisation.
however, little is known regarding the changes in mirnas and their regulation of gene expression and phenotypic variation during wide hybridisation and chromosome doubling, and these data are essential for elucidating the mechanism of heterosis.
in a previous study, we performed wide hybridisation between b. rapa  and b. nigra  and obtained an allodiploid  and allotetraploid .
in addition, we showed that chromosome doubling resulted in higher levels of genetic and phenotypic variation compared with wide hybridisation  <cit> .
in this study, we first analysed the allodiploids and allotetraploids using sequence-related amplified polymorphism  and methylation-sensitive amplification polymorphism  to determine the differences in the genetic changes and epigenetic alterations between wide hybridisation and chromosome doubling.
second, the allodiploids and allotetraploids were analysed through the high-throughput sequencing of small rnas to determine how the changes in small rnas occur during these two processes.
different genomes were combined, the genome dosage was doubled, and the correlation between the sirna and dna methylation at different polyploid levels was assessed.
third, the different expression levels of known mirnas and their target genes were analysed to explore how mirnas and their target genes affect the different phenotypes of the allodiploids and allotetraploids.
RESULTS
phenotypic analysis of the parents and their wide hybrids
in our previous study, the wide hybridisation of b. rapa  and b. nigra  was performed, and the allodiploid  and allotetraploid  were obtained.
we aimed to determine the phenotypic differences between the wide hybrid and their parents.
the characteristics of the allodiploid and allotetraploid and their parents were compared .
the results showed that allotetraploids had a high leaf length and flower size compared with their parents and the allodiploids .
in our previous study, we found that allotetraploids had greater vigour than their parents and allodiploids  <cit> .
thus, allotetraploids had high heterosis compared with their parents and the allodiploids.figure  <dig>
the layout of the experiment plants; the parents and their allodiploid and allotetraploid.
a. leaves and b. flower of the parents and their allodiploid and allotetraploid.
allotetraploids had a high leaf length and flower size compared with their parents and the allodiploids.
dna methylation patterns of the parents and their wide hybrids
to elucidate the epigenetic mechanisms related to the processes of hybridisation and polyploidisation, methylation-sensitive amplification polymorphism  analysis was used to analyse the parents and their allodiploid  and allotetraploid .
after treatment with double-restriction ecori/mspi or ecori /hpaii, the amplified fragments were classified as one of four types:  non-methylated in all of the samples,  methylated in all of the samples,  demethylated in the hybrids compared with the parents, and  hyper-methylated in the hybrids compared with the parents .
in this study,  <dig> reproducible and clear loci were obtained using  <dig> primer pairs .
these  <dig> loci were classified into four major groups , including  <dig> categories according to the variation model between the parents and their allodiploids and allotetraploids .
group a consisted of  <dig> % and  <dig> % of the monomorphic loci in ab and aabb, respectively.
in group b,  <dig> % of the polymorphic loci were specifically found in ab and aabb.
of the loci in group c,  <dig> % and  <dig> % of the polymorphic loci were specifically found in ab and aabb, respectively.
compared with the parents, the group d loci displayed alterations in dna methylation only in  <dig> % and  <dig> % of those found in ab and aabb, respectively.
there was a significant difference in the methylation patterns between ab and aabb.
with respect to the dna methylation status, the ratios between the allodiploids, the allotetraploids, and their parents demonstrated significant differences .
in addition, the cg methylation was high  in ab compared with aabb .
thus, the dna methylation alteration in ab was significantly higher compared with that in aabb.
moreover, the genetic study revealed that two types of the fragments can be used to estimate the genomic changes in the allodiploids and allotetraploids, which have a loss of fragments compared with the parents in addition to novel fragments.
the percentages of genetic changes were  <dig> % in ab and  <dig> % in aabb .
thus, the percentage of genetic changes was significantly higher in aabb compared with ab.figure  <dig>
dna methylation in the parents and their allodiploid and allotetraploid.
the percentage of mirnas, sirnas, genetic changes, and cg methylation in the parents and their allodiploid and allotetraploid.
high-throughput sequencing of small rnas
a small rna library was prepared for the analysis of the parents and their allodiploids and allotetraploids.
sixty million reads were obtained from the small rna sequencing of the four libraries .
a total of  <dig> , <dig> reads were obtained, and  <dig> , <dig> of these reads were of high quality and corresponded to  <dig> , <dig> unique sequence tags.
the small rna sizes ranged from  <dig> to 30 nt, which included the two prominent classes of 21-nt and 24-nt long small rnas .
the 21-nt class corresponded predominantly to mirnas, whereas the 24-nt class corresponded to sirnas  <cit> .
the 24-nt long small rnas were most abundant within the two parental and ab libraries, whereas the 21-nt long small rnas were most prevalent in the aabb library.
interestingly, the amount of mirnas relative to the total small rnas increased with increasing levels of polyploidy: the lowest percentages were  <dig> % in aa,  <dig> % in bb,  <dig> % in ab, and  <dig> % in aabb.
however, the amount of sirnas relative to the total small rna decreased with the increasing levels of polyploidy: the highest percentages were  <dig> % in aa,  <dig> % in bb,  <dig> % in ab, and  <dig> % in aabb .
moreover, the total small rnas demonstrated a high interaction of  <dig> % in aa/ab, and the unique small rnas showed a high level of interaction of  <dig> % in aa/ab .
the genome-matched small rna tags were then clustered into several rna categories  in the four libraries .
in addition, a high percentage of small rnas were sorted as unann rnas .
the repeat-associated srnas were matched based on ltr/copia:  <dig>  ltr/copia:  <dig>  ltr/gypsy:  <dig>  and ltr/gypsy:  <dig> in both the unique tags and total tags.
unexpectedly, all four types of repeat-associated srnas accumulated in lower levels in aabb compared with ab .
among the four types, the accumulation of 21-nt srnas was higher compared with that of 24-nt srnas.figure  <dig>
the length of the small rnas in the parents and their allodiploid and allotetraploid.
known mirnas
to identify conserved mirnas in the parents and their allodiploids and allotetraploids, small rnas that were 18– <dig> nucleotides in length were searched using blastn against mirbase version v <dig> .
the  <dig>  and  <dig>  unique sequences  that were found in the aa and bb libraries, respectively, were annotated as mirna candidates.
totals of  <dig>  and  <dig>  unique sequences  were found in the ab and aabb libraries, respectively .
the expression of known mirnas in the four samples was demonstrated by plotting the log2-ratio .
these results showed  <dig> differentially expressed mirnas:  <dig>   <dig>   <dig>   <dig>   <dig>  and  <dig> up-regulated mirnas and  <dig>   <dig>   <dig>   <dig>   <dig>  and  <dig> down-regulated mirnas in bb/aa, aa/ab, aa/aabb, bb/ab, bb/aabb, and ab/aabb, respectively .
these findings were significantly different between the four libraries.
in our study, further analysis identified totals of  <dig>   <dig>   <dig>  and  <dig> conserved mirnas, which belong to  <dig> plants families, in the aa, bb, ab, and aabb libraries, respectively .figure  <dig>
the expression levels of the known mirnas in the parents and their allodiploid and allotetraploid.
novel mirnas
the expression of novel mirnas in the four samples was demonstrated by plotting the log2-ratio .
these results showed  <dig> novel mirnas:  <dig>   <dig>   <dig>   <dig>   <dig>  and  <dig> up-regulated mirnas and  <dig>   <dig>   <dig>   <dig>   <dig>  and  <dig> down-regulated mirnas in bb/aa, aa/ab, aa/aabb, bb/ab, bb/aabb, and ab/aabb, respectively .
these findings were significantly different between the four libraries.figure  <dig>
the expression levels of the novel mirnas in the parents and their allodiploid and allotetraploid.
qrt-pcr analysis of mirnas and their target genes
many mirna targets play important roles in the regulation of their expression.
the targets of the differentially expressed mirnas were predicted to elucidate the relationship between functions and phenotypes.
the prediction revealed a total of  <dig> known mirnas and  <dig> target genes:  <dig>   <dig>   <dig>  and  <dig> known mirnas and  <dig>   <dig>   <dig>  and  <dig> target genes in aa, bb, ab, and aabb, respectively .
moreover, a total of  <dig> novel mirnas and  <dig> target genes were detected:  <dig>   <dig>   <dig>  and  <dig> novel mirnas and  <dig>   <dig>   <dig>  and  <dig> target genes in aa, bb, ab, and aabb, respectively .
furthermore, we performed a quantitative analysis of nine mirnas and  <dig> targets genes involved in different vegetative and reproductive functions.
moreover, the expression levels of the target genes  were higher in ab compared with aabb.
however, the expression levels of sbp and ago <dig> were higher in aabb compared with ab .
although the corresponding accumulated abundances of various mirnas  were lower in ab compared with aabb, the levels of mir156a and mir157a were higher in ab compared with aabb .
moreover, only mir165a, mir166a, and mir395a were not inversely correlated with the targeted genes  in aabb and ab, respectively.
thus, the expression levels of the targets were negatively corrected with the abundance of significantly expressed mirnas in this study.figure  <dig>
the expression of mirna targets in the parents and their allodiploid and allotetraploid.
the expression levels of the mirnas in the parents and their allodiploid and allotetraploid.
discussion
wide hybridisation and polyploidisation showed different patterns of dna methylation
in this study, the sirnas and dna methylation patterns were significantly different between the allodiploid  and the allotetraploid .
dna methylation was mediated by 24-nt sirnas, which are derived from repetitive dna, transposons, and intergenic and genic regions  <cit> .
based on the length distribution of small rnas, 24-nt sirnas ranged from  <dig> %  to  <dig> %  .
the repeat-associated sirnas were matched on the ltr  and showed higher levels in ab compared with aabb .
the ltr can be reactivated by hybridisation, as has been demonstrated by several previous studies  <cit> .
thus, a hypothesis may be reached that the extent of retrotransponson activation varies depending on the wide hybridisation.
thus, the regulation of genomic dosage may display different patterns in ab and aabb.
dna methylation combined with the activation of a transposable element has been proposed as the stabilising mechanism underlying the epigenetic changes mediated by sirnas during hybridisation and polyploidisation  <cit> .
in response to genomic shock, the sirnas maintain genomic stability in allopolyploids  <cit> .
thus, it can be speculated that sirnas and dna methylation are conducive to the maintenance of genome stability in aabb when faced with genomic shock in ab.
moreover, reduced sirnas levels are largely associated with genes, and genes are associated with altered sirnas levels and have correlations with changes in the dna methylation and expression levels  <cit> .
in the present study, the detected sirnas and dna methylation were low in the allotetraploids compared with the allodiploids.
thus, it has been proposed that sirnas play a key role in maintaining the genomic stability of allotetraploids.
changes in sirna and mirna during wide hybridisation and polyploidisation
sirnas and mirnas induce rapid and dynamic changes during the early stage of allopolyploid formation  <cit> .
the most abundant small rnas found in this study through high-throughput sequencing were identified as mirnas and sirnas, as assessed by the composition of mirnas and sirnas, which were often 21-nt and 24-nt in length, respectively  <cit> .
the relative amount of small rnas corresponding to mirnas increased with an increase in the polyploidy level .
conversely, the relative amount of sirnas corresponding to transposons decreased with an increase in the polyploidy level .
similarly, cantu et al.
<cit>  demonstrated that sirnas corresponding to transposons are observed at lower levels in hexaploids compared with tetraploids.
moreover, the relative amount of epigenetic alterations decreased with an increase in the polyploidy level .
kenan-eichler et al.
<cit>  also reported that the percentage of sirnas and epigenetic changes decreased and the mirna levels increased in hexaploids.
thus, the mirna expression increased with an increase in the polyploidy level, whereas the sirna and epigenetic alteration levels decreased with an increase in the polyploidy level.
furthermore, this ploidy dependence was insensitive to genomic composition but was sensitive to dosage, such as that of ab and aabb, which featured the same two genomes  at varying doses , resulting in their expression of divergent small rnas profiles.
the expression of mirnas and their target genes in the allodiploids and allotetraploids
mirnas function as negative regulators of gene expression and are known to play markedly expanded roles in a variety of developmental processes affecting meristems, leaves, roots, and inflorescences  <cit> .
the evolutionary conserved mir <dig>  mir <dig>  and mir <dig> regulate the development of leaves and contribute to the construction of leaf morphology.
the overexpression of mir <dig>  mir <dig> and mir <dig> reduce the levels of all cuc <dig>  phb, phv, and rev genes and increase the development of sam, which has effects on leaf development  <cit> .
similarly, the mir164a, mir165a, and mir166a levels were high and the cuc <dig>  phb, phv, and rev levels were low in aabb compared with ab .
in our previous study, ab and aabb presented significantly different phenotypes with respect to leaf length and width   <cit> .
in the vegetative to the reproductive phase, the targets of mir <dig>  mir <dig>  mir <dig> and mir <dig> participate in the activation of floral meristem identity genes.
in arabidopsis, as development proceeds, the decrease in the mir156/mir <dig> levels and the increase in spls in the sam result in the activation of floral meristem identity genes  <cit> .
a high expression level of mir <dig> is known to have a positive effect on the timing and development of flowers downstream of nfya <dig>  <cit> .
similarly, in our phenotypic study, flower development was high in aabb compared with ab .
moreover, the overexpression of mir <dig> results in a delay in flowering, which is associated with a reduction in the levels of myb  <cit> .
in our findings, flower timing was delayed in aabb compared with ab .
a large number of mirnas from diverse plants have been identified in the response to metal stress  <cit> .
recent studies have shown that mir <dig> is induced by sulphur starvation and regulates a low-affinity sulphate transporter  and three atp sulphurylases   <cit> .
furthermore, transgenic plants over-expressing mir <dig> accumulate more sulphate in the plant shoots, which suggests that mir <dig> is involved in the regulation of sulphate allocation by targeting aps genes and sultr2; <dig>  <cit> .
in our study, aabb plants were found to overexpress mir <dig> and presented low expression of its target genes  compared with ab plants .
these results showed that aabb plants exhibit tolerance to sulphate deficiency and heavy metal stress.
previous studies have shown that ago <dig> mrna is targeted by mir <dig>  <cit> .
moreover, the regulation of mirnas and their targets may result in novel phenotypes in allopolyploids  <cit> .
in addition, aabb plants presented a higher level of ago <dig> compared with the ab plants .
taken together, the data suggest that mirnas regulate gene expression and induce phenotype variation, such as heterosis, in allotetraploids.
this study explored the role of small rnas in wide hybridisation and polyploidisation between b. rapa and b nigra.
the different mirnas showed different expression levels in allodiploids and allotetraploids, which performed the phenotypic variation.
however, some questions remain elusive, including how nonadditively expressed mirnas and sirnas affect growth and developmental traits, such as leaf shape, plant stature, biomass, flowering time, and fitness in allodiploids and allotetraploids, and whether and how sirnas and dna methylation play roles in this process.
therefore, further studies of the allodiploids and allotetraploids will be necessary.
CONCLUSIONS
this study explored the role of small rnas in wide hybridisation and allopolyploidisation between brassica rapa  and brassica nigra .
when the a genome was crossed with the b genome, sirna levels increased and decreased relative to their parents; b. rapa and b. nigra, respectively, while the dna methylation levels increased relative to their parents.
when the genome ab was doubled, the sirna and dna methylation levels of the allotetraploid decreased compared with its allodiploid.
when the a genome was crossed with the b genome, the mirna levels increased relative to their parents.
when the genome ab was doubled, the mirna levels of the allotetraploid increased compared with its allodiploid.
this result showed that sirnas, dna methylation and mirna play key roles in maintaining the genomic stability through the regulation of small rna levels.
moreover, most mirnas were highly overexpressed in the allotetraploid, which might be induced by the heterosis, such as mir <dig>  mir <dig>  and mir <dig>  mir <dig>  and mir <dig>  which have a major role in flower and leaf development in the allotetraploid.
taken together, the findings of this study demonstrated that sirnas and mirnas maintain the genomic and phenotypic stability in the allotetraploid.
therefore, the present findings may provide new information for elucidating the effects of small rnas on the formation of allopolyploidisation.

 wide hybridisation and polyploidisation is a common phenomenon in plant evolution that results in the formation of new species. wide hybrids often exhibit more vigorous growth than their parents, and this effect is mainly demonstrated by increases in drought tolerance, pest resistance, flowering time, organ size and biomass, among other factors. furthermore, polyploids show novel traits that are not present in their diploid progenitors. for example, allotetraploid cotton produces more abundant and higher-quality fibres, and this effect is derived from their aa and dd extant diploid species. however, wide hybrids also exhibit the disadvantage of infertility, as was observed with the wide hybridisation combination of b. rapa and b. oleracea. however, despite this finding, the growth and adaptability advantages of polyploidy have always been a puzzling phenomenon, and the underlying molecular mechanisms are among the most interesting subjects in plant breeding. allopolyploid formation consists of two processes, namely, hybridisation and chromosome doubling. hybridisation involves a different genomic combination in the same cells and genomes experience “genomic shock”, whereas chromosome doubling doubles and restructures the genome dosage. recent studies have shown that small rnas, particularly the 24-nt sirnas, play an important role in genome reconstruction and stabilisation. the role of 24-nt sirnas is primarily reflected in two aspects. the first aspect is the modification of transposons and repetitive sequences for the maintenance of genome stability, which is mediated by rna-dependent dna methylation, and the other aspect is the cis-regulation of gene expression via the transposon gene fragment in the gene near the region. the reduction of sirnas in f may induce phenotypic and genomic instability in arabidopsis allotetraploids. in interspecific hybrids of arabidopsis, the sirna populations underwent rapid changes in the f allotetraploid, but stability was maintained in the f allotetraploid, in which the dna and chromatin were significantly modified. however, in wheat, the number of 24-nt sirna transcripts significantly decreased in the hexaploid compared with that obtained in the parents, and this decrease was accompanied by a decrease in the dna chain cpg island methylation levels. thus, these studies concluded that a decrease in the dna methylation levels may be one of the causes of genomic instability in allopolyploids during the early stage. this phenomenon of changes to the sirnas with an increase in the polyploidy level suggests that sirnas may play a key role in genomic reconstruction and stability after chromosome doubling. however, little is known regarding how sirnas play this role during wide hybridisation and chromosome doubling. mirnas are the primary mediators of the trans-regulation of gene expression. gene silencing mediated by mirnas is an important strategy used at the post-transcriptional level of gene regulation. furthermore, mirnas are conserved in evolution but become active in polyploidisation. changes in the mirna expression levels can affect the expression of the target genes, and this effect is considered to be one of the main causes that results in the phenotypic variation of the polyploidy. in arabidopsis, the number of mirnas in allotetraploids is higher than those observed in their autotetraploid parents. a similar phenomenon was also found in the synthesis of hexaploid wheat: the number of mirnas increased with an increase in the ploidy level, and this effect was not associated with the gene dosage balance hypothesis. this phenomenon showed that other mechanisms, such as cis-epigenetic regulation, may exist. thus, small rnas, which are a product of non-coding rnas, are involved in regulating gene expression and have become an important factor of gene expression during allopolyploidisation. however, little is known regarding the changes in mirnas and their regulation of gene expression and phenotypic variation during wide hybridisation and chromosome doubling, and these data are essential for elucidating the mechanism of heterosis. in a previous study, we performed wide hybridisation between b. rapa and b. nigra and obtained an allodiploid and allotetraploid. in addition, we showed that chromosome doubling resulted in higher levels of genetic and phenotypic variation compared with wide hybridisation. in this study, we first analysed the allodiploids and allotetraploids using sequence-related amplified polymorphism and methylation-sensitive amplification polymorphism to determine the differences in the genetic changes and epigenetic alterations between wide hybridisation and chromosome doubling. second, the allodiploids and allotetraploids were analysed through the high-throughput sequencing of small rnas to determine how the changes in small rnas occur during these two processes. different genomes were combined, the genome dosage was doubled, and the correlation between the sirna and dna methylation at different polyploid levels was assessed. third, the different expression levels of known mirnas and their target genes were analysed to explore how mirnas and their target genes affect the different phenotypes of the allodiploids and allotetraploids.  phenotypic analysis of the parents and their wide hybrids in our previous study, the wide hybridisation of b. rapa and b. nigra was performed, and the allodiploid and allotetraploid were obtained. we aimed to determine the phenotypic differences between the wide hybrid and their parents. the characteristics of the allodiploid and allotetraploid and their parents were compared. the results showed that allotetraploids had a high leaf length and flower size compared with their parents and the allodiploids. in our previous study, we found that allotetraploids had greater vigour than their parents and allodiploids. thus, allotetraploids had high heterosis compared with their parents and the allodiploids.figure the layout of the experiment plants; the parents and their allodiploid and allotetraploid. a. leaves and b. flower of the parents and their allodiploid and allotetraploid. allotetraploids had a high leaf length and flower size compared with their parents and the allodiploids. dna methylation patterns of the parents and their wide hybrids to elucidate the epigenetic mechanisms related to the processes of hybridisation and polyploidisation, methylation-sensitive amplification polymorphism analysis was used to analyse the parents and their allodiploid and allotetraploid. after treatment with double-restriction ecori/mspi or ecori /hpaii, the amplified fragments were classified as one of four types: non-methylated in all of the samples, methylated in all of the samples, demethylated in the hybrids compared with the parents, and hyper-methylated in the hybrids compared with the parents. in this study, reproducible and clear loci were obtained using primer pairs. these loci were classified into four major groups, including categories according to the variation model between the parents and their allodiploids and allotetraploids. group a consisted of % and % of the monomorphic loci in ab and aabb, respectively. in group b, % of the polymorphic loci were specifically found in ab and aabb. of the loci in group c, % and % of the polymorphic loci were specifically found in ab and aabb, respectively. compared with the parents, the group d loci displayed alterations in dna methylation only in % and % of those found in ab and aabb, respectively. there was a significant difference in the methylation patterns between ab and aabb. with respect to the dna methylation status, the ratios between the allodiploids, the allotetraploids, and their parents demonstrated significant differences. in addition, the cg methylation was high in ab compared with aabb. thus, the dna methylation alteration in ab was significantly higher compared with that in aabb. moreover, the genetic study revealed that two types of the fragments can be used to estimate the genomic changes in the allodiploids and allotetraploids, which have a loss of fragments compared with the parents in addition to novel fragments. the percentages of genetic changes were % in ab and % in aabb. thus, the percentage of genetic changes was significantly higher in aabb compared with ab.figure dna methylation in the parents and their allodiploid and allotetraploid. the percentage of mirnas, sirnas, genetic changes, and cg methylation in the parents and their allodiploid and allotetraploid. high-throughput sequencing of small rnas a small rna library was prepared for the analysis of the parents and their allodiploids and allotetraploids. sixty million reads were obtained from the small rna sequencing of the four libraries. a total of, reads were obtained, and, of these reads were of high quality and corresponded to, unique sequence tags. the small rna sizes ranged from to 30 nt, which included the two prominent classes of 21-nt and 24-nt long small rnas. the 21-nt class corresponded predominantly to mirnas, whereas the 24-nt class corresponded to sirnas. the 24-nt long small rnas were most abundant within the two parental and ab libraries, whereas the 21-nt long small rnas were most prevalent in the aabb library. interestingly, the amount of mirnas relative to the total small rnas increased with increasing levels of polyploidy: the lowest percentages were % in aa, % in bb, % in ab, and % in aabb. however, the amount of sirnas relative to the total small rna decreased with the increasing levels of polyploidy: the highest percentages were % in aa, % in bb, % in ab, and % in aabb. moreover, the total small rnas demonstrated a high interaction of % in aa/ab, and the unique small rnas showed a high level of interaction of % in aa/ab. the genome-matched small rna tags were then clustered into several rna categories in the four libraries. in addition, a high percentage of small rnas were sorted as unann rnas. the repeat-associated srnas were matched based on ltr/copia: ltr/copia: ltr/gypsy: and ltr/gypsy: in both the unique tags and total tags. unexpectedly, all four types of repeat-associated srnas accumulated in lower levels in aabb compared with ab. among the four types, the accumulation of 21-nt srnas was higher compared with that of 24-nt srnas.figure the length of the small rnas in the parents and their allodiploid and allotetraploid. known mirnas to identify conserved mirnas in the parents and their allodiploids and allotetraploids, small rnas that were 18– nucleotides in length were searched using blastn against mirbase version v. the and unique sequences that were found in the aa and bb libraries, respectively, were annotated as mirna candidates. totals of and unique sequences were found in the ab and aabb libraries, respectively. the expression of known mirnas in the four samples was demonstrated by plotting the log2-ratio. these results showed differentially expressed mirnas: and up-regulated mirnas and and down-regulated mirnas in bb/aa, aa/ab, aa/aabb, bb/ab, bb/aabb, and ab/aabb, respectively. these findings were significantly different between the four libraries. in our study, further analysis identified totals of and conserved mirnas, which belong to plants families, in the aa, bb, ab, and aabb libraries, respectively.figure the expression levels of the known mirnas in the parents and their allodiploid and allotetraploid. novel mirnas the expression of novel mirnas in the four samples was demonstrated by plotting the log2-ratio. these results showed novel mirnas: and up-regulated mirnas and and down-regulated mirnas in bb/aa, aa/ab, aa/aabb, bb/ab, bb/aabb, and ab/aabb, respectively. these findings were significantly different between the four libraries.figure the expression levels of the novel mirnas in the parents and their allodiploid and allotetraploid. qrt-pcr analysis of mirnas and their target genes many mirna targets play important roles in the regulation of their expression. the targets of the differentially expressed mirnas were predicted to elucidate the relationship between functions and phenotypes. the prediction revealed a total of known mirnas and target genes: and known mirnas and and target genes in aa, bb, ab, and aabb, respectively. moreover, a total of novel mirnas and target genes were detected: and novel mirnas and and target genes in aa, bb, ab, and aabb, respectively. furthermore, we performed a quantitative analysis of nine mirnas and targets genes involved in different vegetative and reproductive functions. moreover, the expression levels of the target genes were higher in ab compared with aabb. however, the expression levels of sbp and ago were higher in aabb compared with ab. although the corresponding accumulated abundances of various mirnas were lower in ab compared with aabb, the levels of mir156a and mir157a were higher in ab compared with aabb. moreover, only mir165a, mir166a, and mir395a were not inversely correlated with the targeted genes in aabb and ab, respectively. thus, the expression levels of the targets were negatively corrected with the abundance of significantly expressed mirnas in this study.figure the expression of mirna targets in the parents and their allodiploid and allotetraploid. the expression levels of the mirnas in the parents and their allodiploid and allotetraploid. discussion wide hybridisation and polyploidisation showed different patterns of dna methylation in this study, the sirnas and dna methylation patterns were significantly different between the allodiploid and the allotetraploid. dna methylation was mediated by 24-nt sirnas, which are derived from repetitive dna, transposons, and intergenic and genic regions. based on the length distribution of small rnas, 24-nt sirnas ranged from % to %. the repeat-associated sirnas were matched on the ltr and showed higher levels in ab compared with aabb. the ltr can be reactivated by hybridisation, as has been demonstrated by several previous studies. thus, a hypothesis may be reached that the extent of retrotransponson activation varies depending on the wide hybridisation. thus, the regulation of genomic dosage may display different patterns in ab and aabb. dna methylation combined with the activation of a transposable element has been proposed as the stabilising mechanism underlying the epigenetic changes mediated by sirnas during hybridisation and polyploidisation. in response to genomic shock, the sirnas maintain genomic stability in allopolyploids. thus, it can be speculated that sirnas and dna methylation are conducive to the maintenance of genome stability in aabb when faced with genomic shock in ab. moreover, reduced sirnas levels are largely associated with genes, and genes are associated with altered sirnas levels and have correlations with changes in the dna methylation and expression levels. in the present study, the detected sirnas and dna methylation were low in the allotetraploids compared with the allodiploids. thus, it has been proposed that sirnas play a key role in maintaining the genomic stability of allotetraploids. changes in sirna and mirna during wide hybridisation and polyploidisation sirnas and mirnas induce rapid and dynamic changes during the early stage of allopolyploid formation. the most abundant small rnas found in this study through high-throughput sequencing were identified as mirnas and sirnas, as assessed by the composition of mirnas and sirnas, which were often 21-nt and 24-nt in length, respectively. the relative amount of small rnas corresponding to mirnas increased with an increase in the polyploidy level. conversely, the relative amount of sirnas corresponding to transposons decreased with an increase in the polyploidy level. similarly, cantu et al. demonstrated that sirnas corresponding to transposons are observed at lower levels in hexaploids compared with tetraploids. moreover, the relative amount of epigenetic alterations decreased with an increase in the polyploidy level. kenan-eichler et al. also reported that the percentage of sirnas and epigenetic changes decreased and the mirna levels increased in hexaploids. thus, the mirna expression increased with an increase in the polyploidy level, whereas the sirna and epigenetic alteration levels decreased with an increase in the polyploidy level. furthermore, this ploidy dependence was insensitive to genomic composition but was sensitive to dosage, such as that of ab and aabb, which featured the same two genomes at varying doses, resulting in their expression of divergent small rnas profiles. the expression of mirnas and their target genes in the allodiploids and allotetraploids mirnas function as negative regulators of gene expression and are known to play markedly expanded roles in a variety of developmental processes affecting meristems, leaves, roots, and inflorescences. the evolutionary conserved mir mir and mir regulate the development of leaves and contribute to the construction of leaf morphology. the overexpression of mir mir and mir reduce the levels of all cuc phb, phv, and rev genes and increase the development of sam, which has effects on leaf development. similarly, the mir164a, mir165a, and mir166a levels were high and the cuc phb, phv, and rev levels were low in aabb compared with ab. in our previous study, ab and aabb presented significantly different phenotypes with respect to leaf length and width. in the vegetative to the reproductive phase, the targets of mir mir mir and mir participate in the activation of floral meristem identity genes. in arabidopsis, as development proceeds, the decrease in the mir156/mir levels and the increase in spls in the sam result in the activation of floral meristem identity genes. a high expression level of mir is known to have a positive effect on the timing and development of flowers downstream of nfya. similarly, in our phenotypic study, flower development was high in aabb compared with ab. moreover, the overexpression of mir results in a delay in flowering, which is associated with a reduction in the levels of myb. in our findings, flower timing was delayed in aabb compared with ab. a large number of mirnas from diverse plants have been identified in the response to metal stress. recent studies have shown that mir is induced by sulphur starvation and regulates a low-affinity sulphate transporter and three atp sulphurylases. furthermore, transgenic plants over-expressing mir accumulate more sulphate in the plant shoots, which suggests that mir is involved in the regulation of sulphate allocation by targeting aps genes and sultr2;. in our study, aabb plants were found to overexpress mir and presented low expression of its target genes compared with ab plants. these results showed that aabb plants exhibit tolerance to sulphate deficiency and heavy metal stress. previous studies have shown that ago mrna is targeted by mir. moreover, the regulation of mirnas and their targets may result in novel phenotypes in allopolyploids. in addition, aabb plants presented a higher level of ago compared with the ab plants. taken together, the data suggest that mirnas regulate gene expression and induce phenotype variation, such as heterosis, in allotetraploids. this study explored the role of small rnas in wide hybridisation and polyploidisation between b. rapa and b nigra. the different mirnas showed different expression levels in allodiploids and allotetraploids, which performed the phenotypic variation. however, some questions remain elusive, including how nonadditively expressed mirnas and sirnas affect growth and developmental traits, such as leaf shape, plant stature, biomass, flowering time, and fitness in allodiploids and allotetraploids, and whether and how sirnas and dna methylation play roles in this process. therefore, further studies of the allodiploids and allotetraploids will be necessary.  this study explored the role of small rnas in wide hybridisation and allopolyploidisation between brassica rapa and brassica nigra. when the a genome was crossed with the b genome, sirna levels increased and decreased relative to their parents; b. rapa and b. nigra, respectively, while the dna methylation levels increased relative to their parents. when the genome ab was doubled, the sirna and dna methylation levels of the allotetraploid decreased compared with its allodiploid. when the a genome was crossed with the b genome, the mirna levels increased relative to their parents. when the genome ab was doubled, the mirna levels of the allotetraploid increased compared with its allodiploid. this result showed that sirnas, dna methylation and mirna play key roles in maintaining the genomic stability through the regulation of small rna levels. moreover, most mirnas were highly overexpressed in the allotetraploid, which might be induced by the heterosis, such as mir mir and mir mir and mir which have a major role in flower and leaf development in the allotetraploid. taken together, the findings of this study demonstrated that sirnas and mirnas maintain the genomic and phenotypic stability in the allotetraploid. therefore, the present findings may provide new information for elucidating the effects of small rnas on the formation of allopolyploidisation. 
recent studies have demonstrated that small rnas, play an important role in maintaining the genome reconstruction and stability. the study showed that sirnas and dna methylation play an important role in maintaining the genome stability in the formation of an allotetraploid. the mirnas regulate gene expression and induce the phenotype variation, which may play an important role in the occurrence of heterosis in the allotetraploid. hybridisation makes a different genome combined in the same cell, and genome “shock” and instability occur during this process, whereas chromosome doubling results in doubling and reconstructing the genome dosage. the small rnas results showed that the expression levels of mirnas increased in the allodiploid and allotetraploid compared with the parents, and the expression levels of sirnas increased and decreased compared with the parents b. rapa and b. nigra, respectively. the phenotypic analysis exhibited that the allotetraploid had high heterosis compared with their parents and the allodiploid. furthermore, qrt-pcr analysis showed that the expression levels of the target genes were negatively corrected with the expressed mirnas. however, to date, little is known regarding the role of small rnas during the process of wide hybridisation and chromosome doubling, which is essential to elucidate the mechanism of polyploidisation. wide hybridisationallopolyploidisationdna methylationsmall rnasissue-copyright-statement© the author 2014 therefore, the genetic and dna methylation alterations and changes in the sirna and mirna were assessed during the formation of an allodiploid and its allotetraploid between brassica rapa and brassica nigra in the present study. an allopolyploid formation consists of the two processes of hybridisation and chromosome doubling. the findings of this study may provide new information for elucidating that the allotetraploids have a growth advantage over the parents and the allodiploids.
 an allopolyploid formation consists of the two processes of hybridisation and chromosome doubling. hybridisation makes a different genome combined in the same cell, and genome “shock” and instability occur during this process, whereas chromosome doubling results in doubling and reconstructing the genome dosage. recent studies have demonstrated that small rnas, play an important role in maintaining the genome reconstruction and stability. however, to date, little is known regarding the role of small rnas during the process of wide hybridisation and chromosome doubling, which is essential to elucidate the mechanism of polyploidisation. therefore, the genetic and dna methylation alterations and changes in the sirna and mirna were assessed during the formation of an allodiploid and its allotetraploid between brassica rapa and brassica nigra in the present study.  the phenotypic analysis exhibited that the allotetraploid had high heterosis compared with their parents and the allodiploid. the methylation-sensitive amplification polymorphism analysis indicated that the proportion of changes in the methylation pattern of the allodiploid was significantly higher than that found in the allotetraploid, while the dna methylation ratio was higher in the parents than the allodiploid and allotetraploid. the small rnas results showed that the expression levels of mirnas increased in the allodiploid and allotetraploid compared with the parents, and the expression levels of sirnas increased and decreased compared with the parents b. rapa and b. nigra, respectively. moreover, the percentages of mirnas increased with an increase in the polyploidy levels, but the percentages of sirnas and dna methylation alterations decreased with an increase in the polyploidy levels. furthermore, qrt-pcr analysis showed that the expression levels of the target genes were negatively corrected with the expressed mirnas.  the study showed that sirnas and dna methylation play an important role in maintaining the genome stability in the formation of an allotetraploid. the mirnas regulate gene expression and induce the phenotype variation, which may play an important role in the occurrence of heterosis in the allotetraploid. the findings of this study may provide new information for elucidating that the allotetraploids have a growth advantage over the parents and the allodiploids. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords wide hybridisationallopolyploidisationdna methylationsmall rnasissue-copyright-statement© the author 2014 