International Journal of Genomics

International Journal of Genomics / 2017 / Article

Research Article | Open Access

Volume 2017 |Article ID 6790478 | 13 pages |

MicroRNA and Putative Target Discoveries in Chrysanthemum Polyploidy Breeding

Academic Editor: Marco Gerdol
Received25 Jun 2017
Accepted30 Oct 2017
Published13 Dec 2017


MicroRNAs (miRNAs), around 22 nucleotides (nt) in length, are a class of endogenous and noncoding RNA molecule that play an essential role in plant development, either by suppressing the transcription of target genes at a transcriptional level or inhibiting translation at a posttranscriptional level. To understand the roles of miRNAs and their target genes in chrysanthemum polyploidy breeding, three sRNA libraries of normal and abnormal embryos after hybridization were performed by RNA-Seq. As a result, a total of 170 miRNAs were identified and there are 41 special miRNAs in cross of paternal chromosome doubling, such as miR169b, miR440, and miR528-5p. miR164c and miR159a were highly expressed in a normal embryo at 18 days after pollination, suggesting the regulatory role at the late stage of embryonic development. miR172c was only detected in the normal embryo at 18 days after pollination, which means that miR172c mainly mediates gene expression in postembryonic development and these genes may promote embryo maturation. Other miRNAs, including miR414, miR2661, and miR5021, may regulate the genes participated in pathways of auxin response and energy metabolism; then they regulate the complex embryonic development together.

1. Introduction

In plants, microRNAs (miRNAs) are a major class of small noncoding RNAs (sRNAs) with 20–22 nt. sRNAs have been identified to control the developmental processes in plants by regulating gene expression [1]. They have the potential to regulate a gene by two ways: (1) posttranscriptional gene silencing (PTGS) by binding to 3 untranslated region (UTR) of messenger RNAs (mRNAs) and repressing the translation of target mRNAs; (2) transcriptional gene silencing (TGS) by epigenetic modifications [2, 3]. PTGS is the main strategy used by miRNAs to regulate gene expression in plant development. Plant mature miRNAs are generated from miRNA precursors that are processed by a ribonuclease DICER-LIKE1, which negatively regulates specific target mRNAs [2].

A number of studies have investigated the key regulatory role of miRNAs in a wide range of growth and development processes in plants, including the regulation of embryonic development [4, 5]. Plant embryo development includes two stages, embryo morphogenesis and seed maturation. In the model plant Arabidopsis thaliana, embryonic pattern formation has been mainly concerned and studied [5]. Based on the embryo shape, the cell division goes through several stages from preglobular to mature embryo, tightly regulated by multiple genes [6]. As a class of small regulatory RNA, the function of miRNAs during plant embryo development has been considerably reported recently [7]. 28,645 mature miRNAs have been discovered and deposited in the public miRNA database miRBase (Release 21, and hundreds of miRNA target in plant embryo. A large number of Arabidopsis mutants of miRNA biogenesis genes have revealed the crucial roles of miRNA during seed morphogenesis and maturation. Willmann et al. reported the earlier timing of embryo maturation in Arabidopsis mutant for strong alleles of DCL1 (DICER-LIKE1) that are required for miRNA biogenesis and demonstrated the negative regulatory role of specific miRNAs during early embryogenesis and later in embryonic development and that miRNAs are key regulators during seed maturation program [8].

Embryonic miRNAs mediate plant embryo development by regulating transcription factors located in special spatial organization and other key developmental regulators [5]. In plants, miR165/166 is one of the best characterized miRNA families, which regulates five class III homeodomain leucine zipper genes (HD-ZIP IIIs) (PHB, PHV, REV, CNA, and ATHB8) [9, 10]. The HD-Zip III gene family regulates apical embryo patterning and organ polarity as well as controls shoot and root apical meristem (SAM and RAM) formation in embryogenesis [1113]. The family of miR160/miR167 regulates the target mRNAs of auxin response factors (ARFs) associated with auxin homeostasis [14, 15]. Auxin response is an important signaling pathway during embryonic pattern formation, embryo development, and seed maturation [16].

Chrysanthemum (Chrysanthemum morifolium) is an economically important flower around the world, with the increase of chrysanthemum consumption, breeders are driven to improve the cultivars’ traits, such as color, size, shape, and tolerance [17]. Artificial distant hybridization is one of the most effective methods to improve and create new cultivars. However, embryo abortion commonly happens during hybridization [18, 19]. Although previous studies analyzed many different reasons leading to plant embryo abortion, including maternal genotypes [20], parent ploidy [21, 22], and gene regulation [6, 23], the molecular mechanisms regulating embryo development are poorly understood. Zhang et al. revealed the gene, protein, and miRNA change in the stage of chrysanthemum embryo development [19, 24] and provided numerous information of embryo abortion in chrysanthemum hybridization breeding. Evidence supporting a major role for chromosome doubling in overcoming chrysanthemum embryo abortion has been obtained from some studies [22, 25], but the functions of miRNA remain largely uncharacterized.

Next generation sequencing (NGS) approaches have been employed to identify individual miRNAs in various samples, and bioinformatics analyses have offered the technical support to predict the miRNA targets [26]. In the present study, we sequenced three small RNA libraries from chrysanthemum embryo in cross C. morifolium “Yuhualuoying” × tetraploid C. nankingense and identified the key miRNAs and targets that may facilitate embryo development in cross of paternal chromosome doubling.

2. Materials and Methods

2.1. Plant Materials and Artificial Hybridization

Artificial hybridization was performed in cross C. morifolium “Yuhualuoying” (♀, 2n = 6X = 54) × tetraploid C. nankingense (♂, 2n = 4X = 36). Here, the male is an autopolyploid generated by colchicine doubling of the diploid C. nankingense (2n = 2X = 18) [19]. After pollination, three samples were collected, corresponding to the normal embryo at 12 days after pollination (DAP) (NE12), normal embryo at 18 DAP (NE18), and abnormal embryo at 18 DAP (AE18). For each sample, we collected 0.2 g independent biological replicates for verification and mixed the rest of triplicate samples (~0.5 g) for RNA-Seq. All of the samples were immediately frozen in liquid nitrogen and stored at −80°C.

2.2. RNA Extraction, Small RNA Library Preparation, and Illumina Sequencing

Total RNA was extracted with TRIzol reagent (Takara Bio Inc., Otsu, Japan) according to the manufacturer’s protocol. Small RNAs with 18–30 nt fragments were enriched by 15% denaturing polyacrylamide gel electrophoresis. After purification, they were ligated to 5 and 3 adaptors and reversed transcribed into cDNA by reverse transcription-PCR (RT-PCR). Three small RNA libraries were constructed and sequenced using the Illumina HiSeq™ 2000 by the Beijing Genomics Institute (BGI) (Shenzhen, Guangdong Province, China).

2.3. Conserved and Novel miRNA Prediction

Firstly, the data cleaning analysis was performed by getting rid of low-quality reads, reads with 5 primer contaminants and poly A, reads without 3 primer and the insert tag, and reads shorter than 18 nt. The small RNA tags with miRNA, rRNA, snRNA, snoRNA, and tRNA were annotated by aligning to GenBank ( and Rfam database ( using all clean reads of 18~30 nt. The number and proportion of each type of sRNAs were calculated in three libraries. Then, to identify the miRNAs in chrysanthemum embryo, miRBase 19.0 ( was used to search the conserved miRNAs by BlantN. Only 90% matched sequences were considered to be conserved miRNAs. To allow the unambiguous mapping of small RNAs to annotations, the priority rule was followed: rRNA (in which GenBank > Rfam) > known miRNA > repeat > exon > intron. Finally, the novel miRNAs were predicted using the Mireap software (; here, the chrysanthemum embryo transcriptome library obtained from the same sample with the present study (the NCBI accession number PRJNA315793) was used as the reference database. In the chrysanthemum embryo, a total of 99,119 unigenes were assembled with a mean length of 550–580 nt, which were used for prediction in the present study.

2.4. Target Prediction and Functional Annotation for miRNAs

The potential target genes of the known miRNAs were predicted by the web tool psRNATarget ( with parameters suggested by Allen et al. [27]. The target genes were identified in chrysanthemum embryo transcriptome dataset [22], and the function of these potential target genes was annotated using the two protein databases, Gene Ontology (GO) ( and Kyoto Encyclopedia of Genes and Genomes (KEGG) (

2.5. Differential Expression of Known miRNA and Their Targets

To find out the differentially expressed miRNAs in three samples, we normalized the expression of miRNA in three samples and obtained the expression of transcript per million (TPM). Normalized expression = actual miRNA count/total count of clean reads1000000. If the miRNA count was zero, it was revised to 0.01 for analysis of differential expression. The fold-change of log2 (sample 1/sample 2) and value from the normalized counts were calculated to determine significant expression changes. Finally, those miRNAs with fold-change > 1 and value <0.05 were considered to be differentially expressed in the two samples. Heatmap is a visual tool reflecting the expression of differential miRNAs. In the present study, the OmicShare tool, a free online platform for data analysis (, was plotted. Based on the transcriptome library, we searched the expression pattern of these target genes regulated by the differentially expressed miRNAs.

2.6. Real-Time Quantitative PCR (qRT-PCR) Validation of miRNAs and Target Genes

We randomly selected 12 miRNAs and 9 target genes to validate the reliability of sRNA sequencing by qRT-PCR. RNA samples used for sequencing were reverse transcribed using PrimeScript miRNA qPCR starter kit ver 2.0 (Takara Bio, Dalian, China). qRT-PCR was performed using the SYBR Premix EX Taq Kit (Takara, Dalian, China), and the PCR amplification was done as described by Zhang et al. [24]. Three biological replicates were performed for each sample, and relative expression levels were calculated by the 2−△△CT method. The chrysanthemum gene EF1α (elongation factor 1a) (GenBank accession number KF305681) was used as a reference, which is stably expressed in chrysanthemum [19]. Special primers were designed using PRIMER3 RELEASE 2.3.4. All of the primers were shown in Table S1.

3. Results

3.1. Sequencing Analysis of sRNA

To explore the regulation of miRNAs in chrysanthemum embryo development when paternal chromosomes were doubled, three cDNA libraries of small RNAs were constructed, which were named NE12, NE18, and AE18. All clean reads were obtained by filtering the low-quality sequences, adapter sequences, and poly-A sequences shorter than 18 nt, which altogether resulted in more than 99.79% of raw reads. When these small RNA tags were mapped to genome, about 20% of reads in each sample were matched. The unique sRNAs matched to genome in the three libraries were 3,864,037, 3,780,667, and 3,734,164, accounting for 9.62%, 10.04%, and 10.29% of all unique sRNAs (Table 1).

CategoryUnique sRNAsPercent (%)Total sRNAsPercent (%)Unique sRNAsPercent (%)Total sRNAsPercent (%)Unique sRNAsPercent (%)Total sRNAsPercent (%)

Mapping to genome371,5609.622,251,58219.52379,57910.042,194,21220.31384,29010.292,383,15021.24

When these unique sRNAs were aligned to the GenBank and miRBase, many types of sRNAs were identified, including miRNA, rRNA, snRNA, snoRNA, and tRNA, but the vast majority of sequences were unannotated. The length of different types of sRNAs was discrepant; since in chrysanthemum embryo, the most abundant classes of sRNA showed the length of 24 nt (dominant siRNA), then 21 nt (mainly miRNAs), and 22 nt (Figure 1). In the present study, the unique miRNAs in each library were taken into consideration for subsequent analysis, and the proportion of miRNAs was 0.26%, 0.24%, and 0.26% in libraries of NE12, NE18, and AE18 (Table 1).

3.2. miRNAs and Target Genes Identified in Three Libraries

In all samples, a total of 170 conversed miRNAs were identified in miRBase, and 130, 131, and 132 miRNAs were expressed in NE12, NE18, and AE18, respectively. 100 miRNAs (accounting for 58.8%) were detected in three samples; however, some were in two different samples, such as 10 miRNAs in NE12 and NE18, 7 miRNAs in NE18 and AE18, and 6 miRNAs in NE12 and AE18. Venn diagram (Figure 2) presented the quantity distribution of conserved miRNAs in chrysanthemum embryo.

Expression level is an important feature to explain the regulation function of miRNAs, and they commonly varied greatly in different samples. Here, the expression of 170 miRNAs with sequences was normalized and analyzed. The most abundant miRNAs in three samples were miR156a, miR157a, and miR166a, with the expression more than a thousand (Table S2). However, some miRNAs highly expressed in a particular sample, such as the expression of miR398b-5p, are 3524 in AE18 but neither in NE12 nor NE18, miR5721 only expressed in NE18, and miR5662 only in NE12.

To study the biological function of miRNA in chrysanthemum embryo development, the sequences of target mRNA were paid attention. A total of 770 target genes were identified and regulated by 88 miRNAs (51.76% of all miRNAs), and miRNA414 had the most targets (347 unigenes), followed were miR5293 (61 targets) and miR5021 (44 targets) (Table S3). 34 target genes were regulated by two miRNAs, such as CL1999. Contig2 was targeted by miR165a and miR166a; unigene10412 was the target of miR156a and miR157a (Table S3).

3.3. Functional Annotation of Target Genes

Gene Ontology (GO) is a standardized classification system of gene used to describe the characteristics of genes and gene products in organisms. The result showed that 517 target genes were classified into three gene ontology categories: biological processes, cellular components, and molecular functions (Figure 3). Of the 40 functional categories, the number of genes in each sample is different, but in most of categories, NE12 had the largest number but least in AE18. In the third ontology of molecular function, some target genes were expressed in a specific sample, such as the genes related to “electron carrier activity” only expressed in NE12, two genes in the category of “enzyme regulator activity” expressed in NE18, but there was no gene regulated “molecular transducer activity” in NE18 (Figure 3).

To understand the biological function of target genes, KEGG analysis, as with GO, was also used to analyze candidate targets. In NE12, a total of 337 target genes had the biological function on 207 KEGG pathways, but there were less target genes and pathways in NE18 and AE18. 211 target genes with 118 pathway annotation in NE18 and 179 annotated target genes with 113 pathways in AE18 (Table S4). Some target genes were annotated in dozens of pathways that are only in library of NE12, including some energy metabolism pathways, such as ko00280 (valine, leucine, and isoleucine degradation), ko00310 (lysine degradation), and ko00410 (beta-alanine metabolism).

3.4. Differentially Expressed miRNAs in Chrysanthemum Embryo

The aim of this study is to analyze the different miRNAs possibly involved in chrysanthemum embryo development, so we identified 112 differentially expressed miRNAs with at least 1.5-fold after standardized expression. The heatmap showed the expression pattern of these miRNAs (Figure 4). They were assembled in three groups depending on the expression trend. In group NE18/AE18, the number of upregulated miRNAs was almost the same as downregulated, and similar quantity distribution occurred in NE12/NE18. However, the difference is in NE12/AE18, in which only 1/3 of the miRNAs were upregulated in NE12 and means more negative genes in AE18.

3.5. Characteristics of Target Genes in Chrysanthemum Embryo Development

miRNAs regulate the plant development by mediating the expression of target genes. Apart from those redundancy regulated by several miRNAs, a total of 770 target genes were identified and annotated by chrysanthemum transcriptome dataset. The expression level and annotation were presented in Table S5. Most of them were regulated by one miRNA, and some were negatively regulated by miRNAs, such as the transcription factor MYB11 (Unigene2183) highly expressed in NE12, but the targeted miR858b was the lowest expression level in it. Another transcription factor WRKY48 (Unigene25838) was expressed to be the highest in AE18 and lowest in NE12, showing the negative regulation by miR414 (Figure 5, Tables S2 and S5). Several targets were regulated by two miRNAs, and there was no difference in the level of expression among samples. For example, unigene21756 was the target of miR5215 and miR5373, and its FPKM is near the three samples, without negative regulation by these two miRNAs.

In order to study the target genes that regulate the embryo development, some differentially expressed target genes (the fold of FPKM > 1.5) were selected according to their Nr/Nt annotation, which may be involve in chrysanthemum embryo development (Table 2). They contained some transcription factor, genes related to energy metabolism and protein synthesis, and some uncharacterized protein. The expression of transcription factor was various; WRKY and NAC were highest in AE18, but MYB was lowest. Some of the genes associated with auxin and ATP synthesis were downregulated in AE18 (Table 2 and Figure 5).

miRNAID of target genesNE12_FPKMNE18_FPKMAE18_FPKMNr/Nt annotation

miR159aUnigene13056_All22.877614.050114.212850S ribosomal protein
miR160aCL8743.Contig2_All5.26587.40085.3906Auxin response factor
miR160aUnigene9514_All17.71318.564412.0923Auxin response factor 10
miR164cUnigene24220_All1.03976.661211.3245NAC transcription factor
miR166aUnigene33904_All2.13793.85971.4223Homeobox-leucine zipper protein
miR167c-3pCL15007.Contig2_All16.057129.348435.2898AMP-activated protein kinase
miR169bCL3042.Contig2_All2.95118.09417.2676Nuclear transcription factor Y
miR169bCL3042.Contig3_All2.51077.06315.6611Nuclear transcription factor Y
miR169bUnigene26767_All11.851930.080426.5044Nuclear transcription factor Y
miR172cUnigene23005_All1.8273.09814.1985Floral homeotic protein APETALA2
miR2653aUnigene7445_All0.43874.618126.8109Solanum lycopersicum RAP2.6-like protein gene
miR2661Unigene35846_All1.88224.57980.9556ATP binding protein
miR2868Unigene46_All10.110817.979328.7056Protein NLP6-like
miR396a-3pCL936.Contig4_All10.61349.88297.7245MAK16 protein-like protein
miR397aUnigene30076_All29.732313.912513.4372Glucan endo-1,3-beta-glucosidase 7-like
miR397aUnigene4252_All1.97325.40822.1149Glycine max laccase-4-like
miR403CL516.Contig5_All26.117643.392545.8309Argonaute protein group
miR414CL4550.Contig3_All12.719519.8898.8036ADP-glucose pyrophosphorylase small subunit 2
miR414CL2586.Contig3_All2.68525.66277.2585ATP binding protein
miR414Unigene19498_All3.00412.92031.0519F-box protein At5g07610
miR414Unigene12916_All31.227455.686366.4644F-box/LRR-repeat protein
miR414Unigene10707_All0.1582.70394.4112Glutamate receptor 2.3
miR414CL7369.Contig4_All15.133111.867723.1584Glycine max surfeit locus protein 2-like
miR414CL4274.Contig1_All0.11415.5741.6227Histone H2B.2
miR414Unigene32715_All2.80179.38132.9311Minichromosome maintenance factor
miR414CL7491.Contig1_All40.176918.414917.7948Proline-rich protein
miR414CL11069.Contig1_All1.5942.20674.0463Protein NLP6-like
miR414Unigene30228_All1.79634.97349.0057Uncharacterized protein
miR414Unigene31458_All33.3706106.7705151.1735Uncharacterized protein
miR414Unigene28799_All0.66221.22.6895Uncharacterized protein
miR414CL15100.Contig1_All2.759911.790810.9824Uncharacterized protein
miR414CL5368.Contig4_All23.391871.71984.0485Uncharacterized protein
miR414CL5528.Contig1_All3.13584.97871.9595Uncharacterized protein
miR414Unigene35998_All2.94440.26950.8154Uncharacterized protein
miR414Unigene19911_All3.76452.68470.7964Uncharacterized protein
miR414Unigene25838_All2.17868.707524.971Transcription factor WRKY
miR414Unigene27592_All109.281189.670672.8937ATP synthase subunit O
miR5021CL13854.Contig1_All20.240911.18049.5598ATP binding protein
miR5021Unigene8421_All8.16326.98483.84PAPS-reductase-like protein
miR5021 miR5293CL11776.Contig1_All39.39934.238634.2098Uncharacterized protein
miR5049cCL559.Contig1_All33.847933.526845.5856Chrysanthemum x morifolium DREBa
miR5049c miR5052aUnigene1703_All1.45032.19024.4178NA
miR5077CL14985.Contig3_All2.67744.74555.6638Global transcription factor group
miR5227CL4864.Contig1_All58.291317.220911.1415ATP-citrate synthase beta chain protein
miR5270aUnigene3852_All2.52765.026113.2176Uncharacterized protein
miR5293CL7596.Contig2_All999.12822027.78282112.3254Dehydrin family protein
miR5293Unigene16818_All9.1697.10123.1716Glycine max reticulon-like protein B21-like
miR5293Unigene31007_All34.27532.3294106.0944Uncharacterized protein
miR5298bUnigene4147_All12.773420.358742.7602Uncharacterized protein
miR5298b miR5380aCL5088.Contig2_All22.250420.752920.0829Nucleosome assembly protein 1-like
miR5380aUnigene7407_All6.097788.9232149.3867Cytochrome P450
miR5380aUnigene19991_All0.60693.0551.0955Ferric reduction oxidase
miR5215CL965.Contig5_All00.56611.5225Disease-resistance protein NRSA1
miR5215 miR5373Unigene21756_All5.96455.66896.2018ATP binding protein
miR5373CL561.Contig3_All39.076227.12225.8912Glycyl-tRNA synthetase 1
miR5640Unigene11631_All1.352938.1984.7781Glutathione S-transferase 1
miR5641Unigene3971_All21.850214.301711.807Glycosyl transferase
miR5641Unigene7799_All1.28211.63845.4077RPL5A-related protein
miR6111-5pUnigene29644_All5.341311.375217.9804Cytokinin-regulated kinase 1
miR6220-3p miR818aCL1844.Contig1_All9.91433.94893.9549Benzoyl CoA benzoic acid benzoyltransferase
miR858bUnigene2183_All19.35489.29955.5141Transcription factor MYB11

3.6. Validation of qRT-PCR

qRT-PCR is an efficient and accurate way to examine the result of RNA-Seq. In the present study, a total of 12 miRNAs were double tested by qRT-PCR, suggesting the high compatibility of the expression pattern between the two methods (Figure 6). Seven of them had the highest expression in AE18, such as miR169b, miR159a, and miR858b. miR167c-3p was verified with the highest expression in NE18 by two methods. miR414 that regulated the most target genes with various expression levels was expressed dominantly in NE12 (Figure 6).

4. Discussion

4.1. Identified miRNAs in Chrysanthemum Embryo

miRNAs can regulate the plant embryonic development, which important and diverse roles have been studied in various species; however, the specific function of individual miRNA is still uncharacterized for stage-specific embryo [5]. Next-generation sequencing makes it easier to identify individual miRNA families from different organs or treated plants. In Pinus taeda [15] and Brassica napus [28], miRNAs were identified from zygotic embryos at late developmental stages. More and more researchers have provided abundant evidences that miRNAs are required for the majority of embryonic cell differentiation and development in Arabidopsis [29, 30]. In chrysanthemum cross breeding, embryonic development is a crucial stage for seed formation, but embryo abortion is prevalent in chrysanthemum distant hybridization [31, 32]. Previous studies explored this barrier at the level of cell structure, gene expression, and miRNA regulation [19, 24], in which 227 miRNAs were identified in hybrid embryos from cross C. morifolium and diploid C. nankingense. Because of the importance of chromosome doubling in cross breeding, in the present study, we performed the hybridization using C. morifolium and tetraploid C. nankingense; 179 miRNAs were identified in three embryo samples (Table S2). Compared with two crosses, less miRNAs were identified after paternal chromosome doubling; as a result, 135 miRNAs were identified simultaneously in two crosses, 44 new miRNAs in the present cross, such as miR169b, miR440, and miR528-5p, but 92 miRNAs expressed in the previous cross were not detected here, such as miR172a, miR172b, and miR391 [24]. These similarities and differences suggest that most of the identified miRNAs in two crosses are the main factor regulating embryonic development in chrysanthemum and those expressed in specific cross may regulate the development depending on whether the male chromosome doubled. To some extent, it implicated that chromosome doubling has an effect on embryo development in distant hybridization by regulation of miRNAs.

4.2. The Role of miRNA Chrysanthemum Embryo

In plant embryo, miRNAs can mediate their downstream targets and regulate the expression of some transcription factors or other key developmental regulators [5]. The overexpression of miRNAs and their target genes have allowed assignment of developmental roles in embryonic, vegetative, and floral organ boundary formation [33, 34], such as miR159, miR164, miR165/166, miR172, and miR319 families. It has demonstrated that miR164 directly regulates NAC domain genes for making the function of normal plant morphogenesis and normal embryonic [34]. In the present study, there was no miR164 detected, but miR164c was expressed in three samples and higher expression level in abnormal embryo (Table S2). Previous study showed the highest expression of miR164c in normal embryo at 18 DAP when the male was diploid C. nankingense [24], suggesting that miR164c may have the important regulatory role at the late stage of embryonic development, and the lower expression level of miR164c in normal embryo at 18 DAP may be beneficial to chrysanthemum embryo development normally. Similar situation happened in miR159a; there was no miR159 detected in chrysanthemum embryo, but the miR159a was differentially expressed in normal and abnormal embryo at 18 DAP, suggesting the connected function with miR164a during chrysanthemum embryonic development.

miR172 showed the regulatory function for early flowering and floral organ identity defects by regulating AP2 and TOE [35]. miR172c was detected in this study, and the expression was zero in samples of NE12 and AE18, which means that miR172c mainly mediates gene expression in postembryonic development and these genes may promote embryo maturation. DCL1 was a key gene for embryo development by embryo lethality, which was regulated by miR163 [36, 37]. However, there was no miR163 identified in chrysanthemum embryo whether the cross has paternal chromosome doubling or not. This result exhibits that DCL1 was not regulated by miR163 and participated in embryonic lethality in chrysanthemum.

4.3. miRNA-Mediated Target Genes in Chrysanthemum Embryo

Research in the last decade has demonstrated that miRNAs make the crucial roles during plant embryogenesis by regulating various genes and pathways [38]. It contains the process of spatial control of differentiation, regulation of auxin responses, and temporal control of differentiation [5]. A study on miR160-resistant ARF17 transgenes showed the defected cotyledons, suggesting that miR160 negatively regulated genes involved in auxin signaling that is critical for proper development of the embryo and cotyledons [39]. Auxin response is a critical biological pathway for embryonic development, and it has also been reported in chrysanthemum [22]. In the present study, miR160 was identified in NE12 and NE18, highly expressed in NE18 (Table S2). The result of target gene prediction showed that the targets of miR160 were auxin response factor with differential expression among samples (Table 2). In the late embryonic development (18 DAP), the target genes (unigene9514) of auxin response factor downregulated compared with NE12. Besides, there are a portion of target genes involved in energy metabolism according to Nr/Nt annotation, and their expression was downregulated in an abnormal embryo (Table 2). Transcriptome provided the evidence of importance of energy synthesis for normal embryo development [22]; here, these identified target genes were regulated by miRNAs, such as miR414, miR2661, and miR5021, also support the significance of energy metabolism for chrysanthemum embryo development.

5. Conclusion

Polyploid breeding will pay more attention from genomic research in the future as rapid advances in the next generation sequencing technology, which makes unprecedented opportunities to explore and understand the regulatory of genomic or transcriptomic changes. As a critical regulatory factor, the function of miRNAs has attracted a lot of attention during plant growth and development. In chrysanthemum distant hybridization, breeders always faced the barriers existed in embryo development. The present study provided some explanation about the embryo abortion and the miRNAs related to embryo development even their target genes. We propose that late embryonic miRNAs, especially miR164a, regulate NAC transcription factor and thereby affect the embryonic development. miR160 mediated the auxin response, and miR414, miR2661, and miR5021 regulate the genes involved in energy metabolism; together, they regulate the embryo development in chrysanthemum hybridization.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this article.

Authors’ Contributions

Fengjiao Zhang and Jingya Zhao contributed equally to this work.


This study was supported by the National Natural Science Foundation of China (31672182, 31171983), the Natural Science Foundation of Jiangsu Province (BK20161449), and the Fundamental Research Funds for the Central Universities (KYTZ201602).

Supplementary Materials

Supplementary 1. Table S1: primers used to verify the expression of miRNAs and target genes.

Supplementary 2. Table S2: expression level of miRNAs identified in chrysanthemum embryo.

Supplementary 3. Table S3: all identified target genes in chrysanthemum embryo.

Supplementary 4. Table S4: KEGG annotation of target genes.

Supplementary 5. Table S5: annotation of all target genes.


  1. B. J. Reinhart, E. G. Weinstein, M. W. Rhoades, B. Bartel, and D. P. Bartel, “MicroRNAs in plants,” Genes & Development, vol. 16, no. 13, pp. 1616–1626, 2002. View at: Publisher Site | Google Scholar
  2. J. C. Carrington and V. Ambros, “Role of microRNAs in plant and animal development,” Science, vol. 301, no. 5631, pp. 336–338, 2003. View at: Publisher Site | Google Scholar
  3. X. Chen, “Small RNAs and their roles in plant development,” Annual Review of Cell and Developmental Biology, vol. 25, no. 1, pp. 21–44, 2009. View at: Publisher Site | Google Scholar
  4. M. W. Jones-Rhoades, D. P. Bartel, and B. Bartel, “MicroRNAs and their regulatory roles in plants,” Annual Review of Plant Biology, vol. 57, no. 1, pp. 19–53, 2006. View at: Publisher Site | Google Scholar
  5. D. Vashisht and M. Nodine, “MicroRNA functions in plant embryos,” Biochemical Society Transactions, vol. 42, no. 2, pp. 352–357, 2014. View at: Publisher Site | Google Scholar
  6. I. Tzafrir, R. Pena-Muralla, A. Dickerman et al., “Identification of genes required for embryo development in Arabidopsis,” Plant Physiology, vol. 135, no. 3, pp. 1206–1220, 2004. View at: Publisher Site | Google Scholar
  7. W. F. Seefried, M. R. Willmann, R. L. Clausen, and P. D. Jenik, “Global regulation of embryonic patterning in Arabidopsis by microRNAs,” Plant Physiology, vol. 165, no. 2, pp. 670–687, 2014. View at: Publisher Site | Google Scholar
  8. M. R. Willmann, A. J. Mehalick, R. L. Packer, and P. D. Jenik, “MicroRNAs regulate the timing of embryo maturation in Arabidopsis,” Plant Physiology, vol. 155, no. 4, pp. 1871–1884, 2011. View at: Publisher Site | Google Scholar
  9. C. A. Kidner and R. A. Martienssen, “Spatially restricted microRNA directs leaf polarity through ARGONAUTE1,” Nature, vol. 428, no. 6978, pp. 81–84, 2004. View at: Publisher Site | Google Scholar
  10. G. K. Zhou, M. Kubo, R. Zhong, T. Demura, and Z. H. Ye, “Overexpression of miR165 affects apical meristem formation, organ polarity establishment and vascular development in Arabidopsis,” Plant and Cell Physiology, vol. 48, no. 3, pp. 391–404, 2007. View at: Publisher Site | Google Scholar
  11. M. J. Prigge, D. Otsuga, J. M. Alonso, J. R. Ecker, G. N. Drews, and S. E. Clark, “Class III homeodomain-leucine zipper gene family members have overlapping, antagonistic, and distinct roles in Arabidopsis development,” The Plant Cell, vol. 17, no. 1, pp. 61–76, 2005. View at: Publisher Site | Google Scholar
  12. L. Turchi, M. Carabelli, V. Ruzza et al., “Arabidopsis HD-Zip II transcription factors control apical embryo development and meristem function,” Development, vol. 140, no. 10, pp. 2118–2129, 2013. View at: Publisher Site | Google Scholar
  13. S. Miyashima, M. Honda, K. Hashimoto et al., “A comprehensive expression analysis of the Arabidopsis MICRORNA165/6 gene family during embryogenesis reveals a conserved role in meristem specification and a non-cell-autonomous function,” Plant and Cell Physiology, vol. 54, no. 3, pp. 375–384, 2013. View at: Publisher Site | Google Scholar
  14. J. W. Wang, L. J. Wang, Y. B. Mao, W. J. Cai, H. W. Xue, and X. Y. Chen, “Control of root cap formation by microRNA-targeted auxin response factors in Arabidopsis,” The Plant Cell, vol. 17, no. 8, pp. 2204–2216, 2005. View at: Publisher Site | Google Scholar
  15. T. J. Oh, R. M. Wartell, J. Cairney, and G. S. Pullman, “Evidence for stage-specific modulation of specific microRNAs (miRNAs) and miRNA processing components in zygotic embryo and female gametophyte of loblolly pine (Pinus taeda),” New Phytologist, vol. 179, no. 1, pp. 67–80, 2008. View at: Publisher Site | Google Scholar
  16. B. Möller and D. Weijers, “Auxin control of embryo patterning,” Cold Spring Harbor Perspectives in Biology, vol. 1, no. 5, article a001545, 2009. View at: Publisher Site | Google Scholar
  17. N. Anderson, “Chrysanthemum,” in Flower breeding and genetics: Issues, challenges and opportunities for the 21st century, pp. 389–437, Springer, Dordrecht, 2007. View at: Publisher Site | Google Scholar
  18. Y. Deng, N. Teng, S. Chen et al., “Reproductive barriers in the intergeneric hybridization between Chrysanthemum grandiflorum (Ramat.) Kitam. and Ajania przewalskii Poljak. (Asteraceae),” Euphytica, vol. 174, no. 1, pp. 41–50, 2010. View at: Publisher Site | Google Scholar
  19. F. Zhang, Z. Wang, W. Dong et al., “Transcriptomic and proteomic analysis reveals mechanisms of embryo abortion during chrysanthemum cross breeding,” Scientific Reports, vol. 4, p. 6536, 2014. View at: Publisher Site | Google Scholar
  20. K. Kärkkäinen, O. Savolainen, and V. Koski, “Why do plants abort so many developing seeds: bad offspring or bad maternal genotypes?” Evolutionary Ecology, vol. 13, no. 3, pp. 305–317, 1999. View at: Publisher Site | Google Scholar
  21. Z. J. Chen, “Molecular mechanisms of polyploidy and hybrid vigor,” Trends in Plant Science, vol. 15, no. 2, pp. 57–71, 2010. View at: Publisher Site | Google Scholar
  22. F. Zhang, L. Hua, J. Fei et al., “Chromosome doubling to overcome the chrysanthemum cross barrier based on insight from transcriptomic and proteomic analyses,” BMC Genomics, vol. 17, no. 1, p. 585, 2016. View at: Publisher Site | Google Scholar
  23. S. L. Stone, L. W. Kwong, K. M. Yee et al., “LEAFY COTYLEDON2 encodes a B3 domain transcription factor that induces embryo development,” Proceedings of the National Academy of Sciences of the United States of America, vol. 98, no. 20, pp. 11806–11811, 2001. View at: Publisher Site | Google Scholar
  24. F. Zhang, W. Dong, L. Huang et al., “Identification of microRNAs and their targets associated with embryo abortion during chrysanthemum cross breeding via high-throughput sequencing,” PloS One, vol. 10, no. 4, article e0124371, 2015. View at: Publisher Site | Google Scholar
  25. S. Liu, F. Zhang, S. Chen, and F. Chen, “Interspecific hybridization between the tetraploid Chrysanthemum nakingense and Ch. grandiflorum ‘Zhongshanzixing’ and the genetic performance of their F1 hybrids,” Scientia Agricultura Sinica, vol. 43, pp. 2500–2507, 2010. View at: Google Scholar
  26. M. Delseny, B. Han, and Y. Hsing, “High throughput DNA sequencing: the new sequencing revolution,” Plant Science, vol. 179, no. 5, pp. 407–422, 2010. View at: Publisher Site | Google Scholar
  27. E. Allen, Z. Xie, A. M. Gustafson, and J. C. Carrington, “microRNA-directed phasing during trans-acting siRNA biogenesis in plants,” Cell, vol. 121, no. 2, pp. 207–221, 2005. View at: Publisher Site | Google Scholar
  28. D. Huang, C. Koh, J. A. Feurtado, E. W. Tsang, and A. J. Cutler, “MicroRNAs and their putative targets in Brassica napus seed maturation,” BMC Genomics, vol. 14, no. 1, p. 140, 2013. View at: Publisher Site | Google Scholar
  29. M. D. Nodine and D. P. Bartel, “MicroRNAs prevent precocious gene expression and enable pattern formation during plant embryogenesis,” Genes & Development, vol. 24, no. 23, pp. 2678–2692, 2010. View at: Publisher Site | Google Scholar
  30. A. S. Rodrigues and C. M. Miguel, “The pivotal role of small non-coding RNAs in the regulation of seed development,” Plant Cell Reports, vol. 36, no. 5, pp. 653–667, 2017. View at: Publisher Site | Google Scholar
  31. C. Q. Sun, F. D. Chen, N. J. Teng, Z. L. Liu, W. M. Fang, and X. L. Hou, “Interspecific hybrids between Chrysanthemum grandiflorum (Ramat.) Kitamura and C. indicum (L.) Des Moul. and their drought tolerance evaluation,” Euphytica, vol. 174, no. 1, pp. 51–60, 2010. View at: Publisher Site | Google Scholar
  32. C. Q. Sun, Z. Z. Huang, Y. L. Wang et al., “Overcoming pre-fertilization barriers in the wide cross between Chrysanthemum grandiflorum (Ramat.) Kitamura and C. nankingense (Nakai) Tzvel. by using special pollination techniques,” Euphytica, vol. 178, no. 2, pp. 195–202, 2011. View at: Publisher Site | Google Scholar
  33. P. Laufs, A. Peaucelle, H. Morin, and J. Traas, “MicroRNA regulation of the CUC genes is required for boundary size control in Arabidopsis meristems,” Development, vol. 131, no. 17, pp. 4311–4322, 2004. View at: Publisher Site | Google Scholar
  34. A. C. Mallory, D. V. Dugas, D. P. Bartel, and B. Bartel, “MicroRNA regulation of NAC-domain targets is required for proper formation and separation of adjacent embryonic, vegetative, and floral organs,” Current Biology, vol. 14, no. 12, pp. 1035–1046, 2004. View at: Publisher Site | Google Scholar
  35. M. J. Axtell and D. P. Bartel, “Antiquity of microRNAs and their targets in land plants,” The Plant Cell, vol. 17, no. 6, pp. 1658–1673, 2005. View at: Publisher Site | Google Scholar
  36. Y. Kurihara and Y. Watanabe, “Arabidopsis micro-RNA biogenesis through Dicer-like 1 protein functions,” Proceedings of the National Academy of Sciences of the United States of America, vol. 101, no. 34, pp. 12753–12758, 2004. View at: Publisher Site | Google Scholar
  37. X. Chen, “MicroRNA biogenesis and function in plants,” FEBS Letters, vol. 579, no. 26, pp. 5923–5931, 2005. View at: Publisher Site | Google Scholar
  38. T. Yang, L. Xue, and L. An, “Functional diversity of miRNA in plants,” Plant Science, vol. 172, no. 3, pp. 423–432, 2007. View at: Publisher Site | Google Scholar
  39. A. C. Mallory, D. P. Bartel, and B. Bartel, “MicroRNA-directed regulation of Arabidopsis AUXIN RESPONSE FACTOR17 is essential for proper development and modulates expression of early auxin response genes,” The Plant Cell, vol. 17, no. 5, pp. 1360–1375, 2005. View at: Publisher Site | Google Scholar

Copyright © 2017 Fengjiao Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

696 Views | 365 Downloads | 6 Citations
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19. Sign up here as a reviewer to help fast-track new submissions.