Skip to main content

Identification and characterization of miRNAs and PHAS loci related to the early development of the embryo and endosperm in Fragaria × ananassa

Abstract

Background

The strawberry fleshy fruit is actually enlarged receptacle tissue, and the successful development of the embryo and endosperm is essential for receptacle fruit set. MicroRNAs (miRNAs) and phased small interfering RNAs (phasiRNAs) play indispensable regulatory roles in plant growth and development. However, miRNAs and phasiRNAs participating in the regulation of strawberry embryo and endosperm development have yet to be explored.

Results

Here, we performed genome-wide identification of miRNA and phasiRNA-producing loci (PHAS) in strawberry seeds with a focus on those involved in the development of the early embryo and endosperm. We found that embryos and endosperm have different levels of small RNAs. After bioinformatics analysis, the results showed that a total of 404 miRNAs (352 known and 52 novel) and 156 PHAS genes (81 21-nt and 75 24-nt genes) could be found in strawberry seed-related tissues, of which four and nine conserved miRNA families displayed conserved expression in the endosperm and embryo, respectively. Based on refined putative annotation of PHAS loci, some auxin signal-related genes, such as CM3, TAR2, AFB2, ASA1, NAC and TAS3, were found, which demonstrates that IAA biosynthesis is important for endosperm and embryo development during early fruit growth. Additionally, some auxin signal-related conserved (miR390-TAS3) and novel (miR156-ASA1) trigger-PHAS pairs were identified.

Conclusions

Taken together, these results expand our understanding of sRNAs in strawberry embryo and endosperm development and provide a genomic resource for early-stage fruit development.

Peer Review reports

Background

Cultivated strawberry (Fragaria × ananassa Duch.) is one of most popular berry crops in the Rosaceae family, and China has been the world’s largest strawberry producer since 1994 [1]. In the botanical sense, the strawberry fruit is neither a berry nor a true fruit; it is actually a group of achenes dotting the receptacle’s surface. Total or partial removal of the achenes from the receptacle resulted in failure of the fleshy fruit to grow and ripen or an abnormal shape unless exogenous auxin was supplied [2, 3]. Achenes are rich sources of auxin, whose concentration peaks during early development and then declines quickly as fruit growth and ripening progress. This is consistent with the levels of the auxin-related amino acid tryptophan (Trp) and the lack of free auxin in the receptacles during all developmental stages [2, 4,5,6]. This auxin level pattern correlates with auxin promoting receptacle enlargement and inhibiting strawberry ripening [7, 8]. In contrast to the auxin level, the abundance of abscisic acid (ABA) in achenes remained constant during early development and then increased at the onset of fruit ripening [9]. Moreover, ABA is self-induced in achenes and coordinates with IAA in the regulation of early fruit development and later ripening in strawberry [9]. As a consequence, achenes are essential to strawberry receptacle enlargement and ripening. Each achene consists of four tissues: a hard and relatively thick ovary wall, a thin seed coat, a single-layer endosperm, and a small embryo [6, 10]. Among them, the endosperm is a nourishing tissue that supports embryo growth and germination, and failure of endosperm development will cause embryo arrest [11]. In addition, endosperm may play a more prominent role than the embryo in the biosynthesis of auxin and gibberellin for strawberry fruit set [10]. In this regard, the role of endosperm is widely accepted as being altruistic. However, the embryo also regulates the development of the endosperm; for example, scanning electron microscopy (SEM) photos that show the decomposition of starch granules in the endosperm cells facing the scutellum of the embryo [12], and VP1 genes expressed in or GA/ABA signalling released from the embryo modulate endosperm development [13,14,15]. Therefore, successful communication between the endosperm and embryo is necessary for seed development and has a profound impact on fruit set and quality.

At present, many studies have revealed that small RNAs (sRNAs), such as miRNAs and short interfering RNAs (siRNAs), are key genetic and epigenetic regulators in diverse biological processes related to plant growth, development, and stress responses [16]. miRNAs are the most abundant sRNAs and show wide conservation in plants. Mature miRNAs are 20–22 nucleotides in length, are derived from single-stranded RNAs, and regulate their target genes through either transcript cleavage or translation inhibition [17]. Distinct from miRNA biogenesis, siRNAs are generated from long double-stranded RNAs, and distinct siRNA species from both strands depend on the Dicer proteins involved. Accordingly, siRNAs can be further grouped into several subclasses, including heterochromatic small interfering RNAs (hc-siRNAs), phasiRNA and other endogenous siRNAs [18]. Among them, hc-siRNAs are predominantly 24 nt long and are mainly derived from repetitive sequences or transposons. They direct de novo cytosine DNA methylation at homologous loci through an RNA-directed DNA methylation (RdDM) pathway [19]. PhasiRNAs are a special type of siRNA that rely on the function of miRNAs in their biogenesis. After cleavage of the target miRNA, the 5’ fragment of the target mRNA is rapidly degraded and the 3’ fragment is converted into dsRNA via the activity of RNA-DEPENDENT RNA POLYMERASE6 (RDR6). Then, the dsRNA is iteratively cleaved by a Dicer protein (DCL4 or DCL5) from the strand containing the cleavage site, producing duplexes of phasiRNAs (DCL4 forms 21-nt phasiRNAs, DCL5 forms 24-nt phasiRNAs). The biogenesis of phasiRNA is guided by miRNAs through either one or two miRNA binding sites [17]. PhasiRNA-producing loci are named PHAS loci, and the first PHAS loci, TAS3, was identified in Arabidopsis as a conserved noncoding genes among land plants. The TAS3 gene targeted by miR390 generates several tasiRNAs that target AUXIN RESPONSIVE FACTOR (ARF) genes [20]. The three components, miR390, TAS3 and ARF, are indispensable regulators in auxin signalling, and the miR390-TAS3-ARF pathway plays an important role in leaf morphogenesis and flower development [21, 22]. PhasiRNAs are also generated from protein-coding genes; for instance, the PENTATRICOPEPTIDE REPEAT (PPR), nucleotide-binding leucine-rich repeat receptor (NB-LRR) and v-myb avian myeloblastosis viral oncogene homologue (MYB) families were shown to generate phasiRNAs in Arabidopsis [23], Medicago [24], apple (Malus × domestica ‘Golden delicious’) [25], litchi (Litchi chinensis ‘Huaizhi’) [26], woodland strawberry (Fragaria vesca) [27] and Rosa rugosa [28].

To date, research on sRNAs in strawberry has mainly focused on fruit ripening [29, 30] and stress [31], but there is little research on seed (embryo and endosperm) development. In this study, we focused on the developing embryo and endosperm at the early development stage (12 days after pollination (DAP)) by genome-wide characterization of miRNAs, their target genes and the PHAS pathway. Our results demonstrate the large miRNA population and some miRNA-PHAS pairs in the embryo and endosperm, which enhances our understanding of plant sRNAs in strawberry seed development.

Results

Library construction and sequencing

To identify known and novel miRNAs in the embryo and endosperm of strawberry, twelve sRNA libraries derived from reciprocally crossed embryo and endosperm tissues of F. × ananassa cultivar ‘Benihoppe’ and ‘Sweet Charlie’ at 12 DAP were constructed and sequenced. Initially, a total of 9,935,499 to 54,785,143 raw reads were obtained for each library (Table S1). After filtering the low-quality reads, adaptor and contaminant sequences, reads with lengths < 16 nt and reads containing polyA in the small RNA fragment, 8,100,514 to 44,609,527 clean reads were obtained for each library (Table S1). We analysed the length distribution of the sequences in the 12 libraries and found that the majority of sRNAs ranged from 20 to 25 nt in length. The concentrated length distribution peaked at 24-nt sRNAs, while 23-nt and 22-nt sRNAs were the second and third most abundant populations in each library, respectively. Furthermore, there were more 24-nt sRNAs present in the embryo than in the endosperm, whereas the 21–22-nt sRNA showed a reversed trend in the embryo and endosperm of BS, which may be due to tissue-specific expression of sRNA in strawberry development (Fig. 1). Moreover, all clean reads were aligned against the GenBank database [32], Rfam 11.0 database [33], and cultivated strawberry genome [34]. An average of 37.16%, 0.76%, 0.14%, and 2.58% reads matched other noncoding RNAs, including rRNAs, tRNAs, snRNAs, and snoRNAs, and these sRNA populations contained more reads in the endosperm than in the embryo at the 12 DAP stage. In addition, a small portion of reads was mapped to protein coding sequences and intron sequences, which are likely to be RNA degradation products (Figure S1).

Fig. 1
figure 1

Length distribution and abundance of sRNAs in the reciprocally crossed embryo and endosperm of strawberry. BS and SB represent ‘Benihoppe’ (♀) × ‘Sweet Charlie’ (♂) and ‘Sweet Charlie’ (♀) × ‘Benihoppe’ (♂), respectively. em: embryo; en: endosperm. For each tissue, three biological replicates (em1, em2, and em3; en1, en2, and en3) were performed

Identification of known and novel miRNAs in the developing embryo and endosperm

Clean reads (excluding reads mapped to snRNAs, snoRNAs, tRNAs and coding sequences or intron sequences) from 12 libraries were aligned against the miRbase 22 database [35]. A total of 352 known miRNAs were identified (Table S2); these miRNAs belong to 46 miRNA families. Twenty-two out of the 46 families are highly conserved in plants, and a few less conserved miRNAs, including miR477, miR482/2118, and miR535, were identified. As shown Fig. 2A, the most abundant miRNA families are miR156 (38 members), miR482 (25 members), and miR171 (24 members). In addition, 52 novel mature potential candidate miRNAs from 12 libraries were identified, of which most were 21 nt in length (Fig. 2B & Table S2). Pri-miRNAs with stem‒loop structures are cleaved into miRNA::miRNA* duplexes, and in this study, we also analysed the miRNA* structure (Fig. 2C & Table S2). The abundance of miRNA* was significantly lower than that of mature miRNAs; in contrast, we found that fan-miR408* and fan-miRN1* were more abundant than their references, fan-miR408 and fan-miRN1. Interestingly, miR408* was also found to have a higher abundance than rco-miR408 in caster been [36].

Fig. 2
figure 2

The number of members in known and novel miRNA families. A A total of 46 known miRNA families were identified in developing strawberry embryos and endosperm. B Fifty-two novel miRNAs belonging to 30 novel miRNA families were identified in embryos and endosperm. Different colours are used to present the length of each miRNA family member. denotes 20 nt miRNA,  21 nt miRNA and  22 nt miRNA. C The secondary structures of fan-miR166a, fan-miR408b, fan-miRN1a, and fan-miRN25 identified from strawberry. Mature miRNA and predicted miRNA* sequences are marked in  and  , respectively

Expression profiling of miRNAs in the developing embryo and endosperm

To explore known miRNA families that are highly expressed in the development of embryos and endosperm in strawberry, the TPM reads of miRNA families were counted by summing. Among them, miR1511 was the most abundant family, followed by miR156 and miR159 (Figure S2). With these expression data, we can examine how highly expressed miRNA families are involved in the embryo and endosperm. These high expression levels of known highly expressed miRNA families were divided into five groups according to their expression profile (Fig. 3A). Among them, group I and group V were highly expressed in the endosperm of SB and BS, respectively, and we considered that this was related to variety specificity. The expression levels of group III miRNA families, including miR164, miR167, miR168, miR396, miR477, miR1511, miR6275, miR7125, miR11289, miR11292 and miR11312, decreased in the embryo but increased in the endosperm, indicating that their functions are potentially involved in the development of the endosperm. In contrast, group IV had a high expression level in the embryo from the reciprocal cross, which suggested that these miRNA families play an important role in the embryo. Novel miRNA families with more than 50 TPM reads were found to be highly expressed in the development of strawberry seeds. Of them, 26 novel miRNA families were classified, and miRN12 was the most abundant family, followed by miRN10 and miRN8. According to the different expression levels between the embryo and endosperm, these novel families were divided into three groups (Fig. 3B). Ten miRNA families showed upregulated expression in the endosperm from the reciprocal cross relative to the embryo, so we placed these in group I. In group II, miRN23 showed tissue-specific expression in the embryo of BS, and miRN24 showed strain-specific expression in SB. Group III miRNA families have upregulated expression in the embryo from reciprocal crosses, showing their regulatory function in heart embryo development. To verify the results of sRNA sequencing, tailing-reaction quantitative reverse transcription PCR (qRT‒PCR) was performed to determine the expression level of 12 differentially expressed miRNAs in the development of strawberry embryos and endosperm. As shown in Fig. 3C, the expression patterns of miRNAs determined by qRT‒PCR were broadly consistent with those from sRNA sequencing data.

Fig. 3
figure 3

Expression profiling of known and novel miRNA families in the developing strawberry embryo and endosperm. A Forty-six known miRNA families were divided into five groups according to their expression patterns in the embryo and endosperm. B Twenty-six novel miRNA families were classified into three groups. C Tailing-reaction qRT‒PCR verification of 12 differentially expressed miRNAs in the development of strawberry embryos and endosperm. The error bar represents the error values of three biological replicates

Prediction and analysis of miRNA target genes in the embryo and endosperm

To analyse the function of miRNAs in strawberry embryo and endosperm development, we predicted targets of the 404 miRNA candidates, including 352 known miRNAs and 52 novel miRNAs; 2984 target genes were detected (Table S3). The majority of conserved miRNA targets are various transcription factors, including SPL, MYB, ARF, NFYA, AP2, TCP2, and HD-ZIP, which play an important role in regulating plant development [16, 37]. Other conserved miRNA targets include NAC (miR164), scarecrow-like protein (miR171), TAS3 (miR390), sulfate transporter (miR395), laccase genes (miR397), umecyanin-like (miR398), ubiquitin-conjugating enzyme or inorganic phosphate transporter (miR399), uclacyanin (miR408) and SPX domain-containing membrane protein/protein TORNADO 2 (miR827). miR162 and miR168 were found to target genes encoding DCL1 and AGO1, respectively, which are key components of the general process of miRNA biogenesis and activity. miR391, miR482, miR2109, miR2118 and three novel miRNA families containing miRN19, miRN25 and miRN29 were found to target disease resistance proteins. F-box genes are a large family in plants; approximately 800 of these genes in F. vesca are involved in different biological events [16, 38] and have been identified to be targeted by the conserved miRNA of miR394, the less conserved miRNAs including miR11288 and miR11315, and two novel miRNA families, miRN1 and miRN20. Two less conserved (miR11285 and miR11301) and two novel (miRN24 and miRN26) miRNA families targeted the PPR protein. Receptor-like serine/threonine protein kinases are a large family of biological enzymes that participate in the phosphorylation of proteins [39]; these kinases were targeted by several miRNA families, including miR390, miR482, miR535, miR7125, miR11284 and miR11308. Moreover, novel miRNAs targeted a number of functional genes, such as NEDD8-specific protease 1 (miRN3), RNA polymerase II degradation factor 1-like (miRN4), sulfate transporter (miRN5), elongation of fatty acids protein (miRN11), ACT domain-containing protein (miRN14), piriformospora indica-insensitive protein (miRN17), protection of telomeres protein (miRN18), and phosphoenolpyruvate carboxylase (miRN31). To further demonstrate that these identified genes were targeted by miRNA, qRT‒PCR was performed to analyse the fan-miR159a-FaMYB101, fan-miR160a-FaARF18, fan-miR164a-FaNAC, fan-miR167a-FaARF8, fan-miR167h-FaARF6, and fan-miR172b-FaAP2 pairs. The qRT‒PCR results showed that these 6 target gene expression profiles were opposite those of their corresponding miRNAs (Fig. 4), which was consistent with previous studies; for instance, miR167 directly targets ARF6 and ARF8 and decreases their transcript levels, which is essential for the fertility of both ovules and anthers [40]. miR172 inhibits the expression of AP2, which limits cell division and expansion in Arabidopsis [41], and overaccumulation of miR172 leads to the silencing of AP2 in apple [42].

Fig. 4
figure 4

Quantitative real-time PCR validation of six identified miRNA‒mRNA pairs in strawberry. Blue and green bars represent the miRNAs and their corresponding target mRNAs, respectively. The error bar represents the error values of three biological replicates

Identification of PHAS loci in the developing embryo and endosperm

To investigate phasiRNA biogenesis and regulation in embryos and endosperm, PHAS loci were identified based on small RNA reads by using sRNAanno [43]. Hundreds of PHAS loci (251 twenty-one-nucleotide PHAS loci and 798 twenty-four-nucleotide PHAS loci) were identified. When using a stringent phase score (phase score >  = 10), 81 21-nt PHAS loci could still be identified. We identified 75 24-nt PHAS loci by using a phase score of 5 (phase score >  = 5) (Table S4).

To refine the putative annotation of the predicted PHAS loci, we aligned these PHAS loci to the NCBI Nucleotide Collection (nr/nt) database and the PlantRep database (Table S4). Among the 21-nt PHAS loci, some PHAS loci were protein coding genes, such as PPR, TAR2, AFB2, NAC and kinase, which are conserved in plants [44, 45] (Fig. 5A). We found that six CM3 genes and one ASA1 gene in strawberry also generated 21-nt phasiRNAs, which play an important role in the shikimic acid pathway [46]. In addition, 36% of the PHAS loci originated from ncRNA, and two identified 21-nt PHAS loci were TAS3 loci; we did not find TAS1 and TAS2 loci in strawberry. The 81 PHAS loci were mapped to the genome using TBtools [47] software, and among them, the majority of PHAS loci were clustered on chromosomes Fvb4-3 and Fvb6-1 (Fig. 5C). Up to 55 PHAS loci overlapped with annotated protein-coding genes of cultured strawberry, including CM1/3, PPR, peptidyl-prolyl cis–trans isomerase FKBP62, TAR2, DNA repair protein, E3 ubiquitin-protein ligase, AFB2, and other proteins (Table S4). A large number (33/55) of PHAS loci overlapped with genes of unknown function, many of which were single copies in the genome, suggesting cis- rather than trans-activity. The remaining 26 PHAS loci were found in intergenic regions of the genome (Table S4).

Fig. 5
figure 5

Classification of identified PHAS loci. A The 21-nt PHAS loci were classified into three categories according to their functional annotation. CM3, chorismate mutase 3; PPR, pentatricopeptide repeat; TAR2, tryptophan aminotransferase-related 2; AFB2, auxin signalling F-box 2; ASA1, anthranilate synthase alpha subunit 1; NAC, NAM, ATAF and CUC proteins. B Classification of 24 identified PHAS loci. Coding (blue) genes were classified according to their annotation. LEC2, B3 domain-containing transcription factor LEC2; RKS1, G-type lectin S-receptor-like serine/threonine-protein kinase RKS1; FT1, flowering locus T1. Noncoding (orange) loci are grouped based on the length of the PHAS loci. C 21-nt PHAS loci located in the strawberry genome. D 24-nt PHAS loci located in the strawberry genome

Many 24-nt PHAS loci originated from protein coding genes, such as LEC2, RKS1, CSC1, FT1 and MYB (Table S4). In addition to protein coding genes, 9 predicted noncoding genes produced abundant phasiRNAs, which were of short length (100 bp ~ 400 bp) based on sRNA data (Fig. 5B). Interestingly, we found that most 24-nt PHAS loci (61.33%) were located in the repeating sequence region, which is different from the 21-nt PHAS loci (Table S4). The majority of PHAS loci were clustered on chromosomes Fvb1-1, Fvb4-2, Fvb6-3 and Fvb7-1 (Fig. 5D). Seventy-two percent (54/75) of the PHAS loci were located in intergenic regions of the genome, and 28% (21/75) of the PHAS loci overlapped with annotated protein-coding genes in cultured strawberry. Interestingly, these annotated protein-coding genes are enzymes, except for uncharacterized proteins, including DNA ligase, isoleucine–tRNA ligase, RNA helicase, protein phosphatase 2C, methylesterase 17, xyloglucan galactosyltransferase KATAMARI1 and key players in epigenetic modification, including the ENHANCED DOWNY MILDEW 2 (EDM2) protein, the helicase protein MOM1, and the methyltransferase PMT9 (Table S4).

The expression patterns of PHAS loci in the developing embryo and endosperm

We also examined the expression of each PHAS locus (Fig. 6 and Figure S3). Among 81 of the 21-nt PHAS loci, 56 showed differential expression in the embryo and endosperm, where 25 and 31 loci accumulated at higher levels in the embryo and endosperm, respectively (Fig. 6A). In this study, we found four PPR loci in strawberry, and only one PPR locus was specifically expressed in embryos (Fig. 6B). Conversely, other conserved PHAS loci, including two TAR2 and one AFB2, showed higher expression in the endosperm (Fig. 6C). In addition, among six CM3 loci, strawberry also had a higher level in the endosperm (Fig. 6D). However, we observed that most of the 24-nt PHAS loci showed embryo-specific expression, and only nine had higher expression in the endosperm, indicating that 24-nt phasiRNAs play important roles in embryo growth and development (Figure S3A). To understand the specific expression patterns of 24-nt PHAS loci, sRNA abundance and phasing score at an endosperm-specific zinc transporter protein locus and embryo-specific MYB locus were examined, as shown in Figure S3B and C, respectively.

Fig. 6
figure 6

Expression patterns of 21-nt PHAS loci in strawberry embryos and endosperm. A Cluster analysis and heatmap of the expression patterns of all 21-nt PHAS loci in embryos and endosperm. B The distribution of sRNA reads in the sRNA profiles from the embryo (left) and endosperm (right) tissues at the PPR locus. Two tracks showing the abundance (RP10M) of sRNAs in both strands of the PPR locus (upper); sRNA sizes are denoted by different colours. The lower track shows the phasing score of each read in this locus. C sRNA abundance and phasing score at the TAR2 locus. The patterns of another TAR2 locus are similar to those of this TAR2 locus, as indicated. D sRNA abundance and phasing score at the CM3 locus. The patterns of other CM3 loci are similar to this CM3 locus, as indicated

Identification of miRNA triggers for PHAS loci

Seventeen families of miRNAs, including 14 known miRNAs and 3 novel miRNAs, were found to trigger phasiRNA production from 21-nt PHAS loci. Among the 40 PHAS loci triggered by miRNAs, 22 overlapped with annotated genes in the strawberry genome, such as E3 ubiquitin-protein ligase, PPR, anthranilate synthase alpha subunit 1 (ASA1), and CM3. Among the miRNA triggers, most miRNA triggers were 21 nt long, including conserved miRNAs, such as miR390, miR319, miR395 and miR482, and novel miRNAs, such as miRN3 and miRN14, and the remaining miRNA triggers were 20 and 22 nt long. All miRNA trigger-PHAS pairs are shown in Table 1. Some conserved miRNA trigger-PHAS pairs were found in strawberry, such as miR390-TAS3 (Fig. 7A) and miR482/2118-ncRNA. In addition, we also found some novel miRNA trigger-PHAS pairs, such as miR11301/miR11285-PPR/CSC1-like (Fig. 7B, C), miR156-ASA1 (Fig. 7D), miR390-HLCS (Fig. 7E), miR319/miR159-MIR159 (Fig. 7F), miRN11-ncRNA/40S ribosomal protein S29 and miRN14-ncRNA (Table 1).

Table 1 Triggers of PHAS loci in strawberry
Fig. 7
figure 7

miRNA-triggered 21-nt phasiRNAs in strawberry. A sRNA abundance and phasing score are viewed at putative TAS3 genes and are targeted by miR390 via a “two-hit”-like model in strawberry. B Novel PHAS locus CSC1-like targeted by miR11285 and miR11301. C Conserved PHAS locus PPR was identified in Fvb3-4 and targeted by miR11301. D miR156n-triggered phasiRNAs produced from the ASA1 locus. E miR390 also triggered phasiRNA production from the novel PHAS locus HLCS. F Novel PHAS locus MIR159c targeted by miR159 and miR390

Twelve of the 24-nt PHAS loci were triggered by miRNA, including 7 known miRNAs and 3 novel miRNAs. Among the miRNA triggers, most miRNA triggers are 21 nt long, such as miRNA395, miRNA482, miRNA6275, and miRNA11288 and the novel miRNAs miRN14 and miRN29. The remaining miRNA triggers were 20 and 22 nt long, including miR156, miR172, miR11285 and miRN16, respectively. Table 1 also shows miRNA trigger-24-nt PHAS pairs, such as miRN29-ncRNA (Figure S4A), miR11285-CSC1-like (Figure S4B), miRN14-ncRNA (Figure S4C) and miR11288/miR395-uncharacterized mRNA (Figure S4D).

Discussion

Bidirectional dialogue between the endosperm and embryo is important for fruit set and development [10, 12,13,14,15, 48]. However, studies that focused on the sRNAs involved in endosperm and embryo development have not been found for Fragaria.

Abundance of sRNA in the developing embryo and endosperm

In this study, we constructed 12 sRNA libraries from reciprocally crossed embryo and endosperm tissues from strawberry and sequenced the libraries using high-throughput sequencing technology to gain insight into the identified miRNAs and PHAS loci that might be involved in strawberry embryo and endosperm development. Here, we found that the majority of the sRNAs ranged from 20 ~ 25 nt in length, with the 24-nt class being the most abundant in the embryo and endosperm (Fig. 1). The result was similar to previous reports from castor bean [36] and maize [49, 50], while in the green fruit stage and red fruit stage of cultivated strawberry, 21-nt sRNAs represented the most abundant sRNAs [51], suggesting that 24-nt sRNA biogenesis is conserved in embryo and endosperm development.

Features of the miRNAs in the developing embryo and endosperm

miRNAs are the most functionally important and most widely studied class of sRNA in plants. Our current study confirmed that 352 known miRNAs belonging to 46 miRNA families were correlated with the embryo and endosperm development processes in strawberry. We also found 52 novel miRNAs that have not been reported elsewhere. Among the miRNA families, the miR156 family was the largest family, with 28 members, followed by miR482 and miR171, with 25 and 24 members, respectively (Fig. 2A). These results were similar to those of previous studies, albeit with some differences. For example, Huang et al. [49] found that the miR166 and miR167 families were the largest families, with 15 members in maize endosperm, followed by the miR156 (13 members) and miR171 (11 members) families. The differences in these miRNA family members in different plants revealed that they could function in common or unique regulatory pathways. Compared to the miRNA target genes identified in other plants, many conserved miRNAs, such as miR156, miR160, miR166, miR167, miR395, miR397, and miR408, etc., showed similar targets as their homologues in strawberry [16, 37]. These conserved miRNAs and their targets usually play fundamental regulatory role in plant growth and development. For instance, the conserved miR156 family is involved in regulation of flowering and the vegetative-to-reproductive transition, by cleaving SPL transcripts and promoting their translational repression [52]. Moreover, miR156 and miR172 are positively regulated by the transcription factors they target, participates in the stability of juvenile-to-adult phase transition in plants [53, 54]. In plants, miR165 and miR166 share the same function of targeting and cleaving HD-ZIP III genes, which act against meristem maintenance [55]. miR397 regulated fruit cell lignification in pear by inhibiting expression of laccase genes which is useful for improving pear fruit quality [56]. However, some conserved miRNAs and most novel miRNAs exhibited species-specific targets in cultivated strawberry. For example, ARF transcription factors, which are essential in the early auxin response, were targeted by the miR160 and miR167 families, which is consistent with our results [16, 57]. Our study showed that miR160 family members were all more highly expressed in the embryo than in the endosperm and targeted ARF17 and ARF18. In the miR167 family, many members were highly expressed in the embryo, while miR167j, miR167k and miR167l were highly expressed in the endosperm. Their target genes were ARF6 and ARF8. Moreover, we found that ARFs were also targeted by miR159. The results of our study showed that miR159 family members, including miR159m, miR159n, miR159o, miR159p, miR159q and miR159r, showed endosperm-specific expression (Figure S5) and targeted ARF23 (Table S3). Previous research showed that ARF23 is a pseudogene in A. thaliana, and its biological function and ability to bind DNA is not known [58,59,60,61]. Among the various functional mechanisms of pseudogenes, antisense transcription-related regulation and competition with miRNA and RNA-binding proteins appear to be common in humans [62, 63]. The genetic functions and evolutionary trends for ARF23 in the strawberry genome are unknown. Additionally, it is unclear how the miR159 family that was predicted to target ARF23 regulates the development of strawberry endosperm and embryos. These are interesting issues that require further investigation.

In addition, we classified the miRNA families according to their expression level profiles to explore diverse molecular and physiological pathways of miRNAs in strawberry embryos and endosperm. We found that some of the conserved miRNA families, such as miR164, miR167, miR168 and miR395, showed conserved expression in the developing reciprocal cross endosperm, and other conserved miRNA families, including miR159, miR160, miR162, miR166, miR171, miR319, miR395, miR399, and miR535, showed conserved expression in the developing reciprocal cross embryo (Fig. 3A). This result is similar to the research on embryos and endosperm in maize [64].

PhasiRNAs involved in the developing embryo and endosperm

In addition to miRNAs, phasiRNAs have been shown to be involved in developmental regulation and stress responses in flowering plants. Genome-wide identification of phasiRNAs is an important tool that can be used to investigate gene regulation involving small RNAs in legumes [24], soybean [65], tomato [66], litchi [26], Panax notoginseng [67], wheat [68], and pummelo [69]. In the woodland strawberry F. vesca, a total of 515 PHAS loci, including 192 21-nt PHAS loci and 323 24-nt PHAS loci, were identified [27]. In our study, 81 21-nt PHAS loci and 75 24-nt PHAS loci were identified in the cultivated strawberry F. × ananassa genome (Table S4). However, only 49.38% (40/81) of the 21-nt PHAS loci and 16% (12/75) of the 24-nt PHAS loci were identified as targets of miRNA triggers. Many studies have revealed that a few phasiRNA pathways are highly conserved alongside conserved miRNAs in land plants [50]. TAS3 is arguably the most well-conserved PHAS locus, as it has been identified across a broad range of species, and the miR390-TAS3-ARF pathway is the most conserved and archetypal [44, 70]. Xia et al. [38] identified two TAS3 genes, named TAS3L and TAS3S, in F. vesca. Cultivated strawberry also has two TAS3 genes targeted by miR390, but we did not find other TAS genes in the embryo and endosperm, which suggested that TAS3 is primarily conserved in Fragaria. A number of PPRs were shown to generate phasiRNAs and to be conserved, especially in eudicots. The pathway for the generation of PPR phasiRNAs includes a small number of diverse trigger miRNAs. For instance, miR161, miR400, and miR173 in Arabidopsis [20, 71, 72], miR1509 in legumes [24], miR7122 in Rosaceae [73], and miR11285 in F. vesca [38]. In this study, our research showed that PPR generated phasiRNAs by triggering miR11285 or miR11301 (Fig. 7C and Table 1); this result was similar to that in F. vesca. We also identified some conserved PHAS loci in strawberry, i.e., TAR2, NAC, and AFB2 (Fig. 5A). As an auxin biosynthetic gene, TAR2 plays a role in the indole-3-pyruvic acid (IPA)-dependent route for the conversion of Trp to IAA, which is the main mechanism for the de novo synthesis of auxin in plants and is required for proper embryo development, lateral root growth and several developmental processes in Arabidopsis [74,75,76,77,78]. As previously described in Arabidopsis, FaTAR2 works as a Trp aminotransferase in strawberry and is involved in auxin biosynthesis. Transient silencing of this gene was accompanied by a diminished responsiveness to auxin in ripening receptacles [48]. NACs are a family of genes specific to plants that play a role in a diverse set of developmental processes; their functions are partially similar to that of TAR2. For example, NAC1 is also involved in the auxin-signalling pathway to promote lateral root development, and in Arabidopsis, miR164-directed regulation of the NAC-domain gene CUC1 is necessary for normal embryonic development [79, 80]. In addition, NAC1, as a transcriptional activator, regulates the expression of two auxin-responsive genes, AIR3 and DBP [79]. AFB2 has positive regulatory roles as an auxin receptor in the regulation of fruit development compared with the AFB gene family [81, 82] and has the highest expression level in the early fruit development stage in cucumber and the diploid strawberry F. vesca [10, 83]. Moreover, AFB2 plays a major role in the control of jasmonic acid conjugation through the ARF6/ARF8 auxin signalling module during Arabidopsis adventitious rooting initiation [84]. In addition, the PHAS genes ASA1 and CM3 were newly identified in strawberry and found to play an important role in the shikimate pathway [46]. The shikimate pathway serves as an important precursor for auxin biosynthesis, and Trp accumulation occurs via this pathway to sustain auxin production [76]. Overall, we suggest that the biosynthesis of auxin-related phasiRNAs is an indispensable regulator of early endosperm and embryo development and of receptacle fruit set in strawberry.

Conclusions

In this study, a total of 404 miRNAs (352 known and 52 novel) could be found in strawberry embryo and endosperm tissues, of which four and nine conserved miRNA families displayed conserved expression in the endosperm and embryo, respectively. We also identified 81 PHAS loci producing 21-nt phasiRNAs and 75 PHAS loci producing 24-nt phasiRNAs. Among them, some auxin signal-related genes, CM3, TAR2, AFB2, ASA1, NAC and TAS3, were found, which demonstrates that PHAS locus-mediated IAA biosynthesis is important for endosperm and embryo during early fruit growth. Additionally, some auxin signal-related conserved (miR390-TAS3) and novel (miR156-ASA1) trigger-PHAS pairs were identified. The above results improve our knowledge of sRNAs in strawberry embryo and endosperm development and provide candidate PHAS loci for the early developmental stage of fruit in F. × ananassa.

Materials and methods

Plant materials and growth conditions

Two cultivated strawberry cultivars, F. × ananassa cultivars ‘Benihoppe’ and ‘Sweet Charlie’, were used as materials for sRNA identification and were grown in the Baima Teaching and Scientific Research Base of Nanjing Agricultural University, Nanjing, China. Developing fruits (receptacles with achenes), BS achenes (‘Benihoppe’ as the female, pollinated with ‘Sweet Charlie’) and SB achenes (‘Sweet Charlie’ as the female, pollinated with ‘Benihoppe’) from reciprocal crosses between ‘Benihoppe’ and ‘Sweet Charlie’ were collected at 12 DAP. The embryo and endosperm of a developing ovary (i.e., achene) were hand-dissected under a stereomicroscope, washed with tissue separation buffer [85], and then frozen immediately in liquid nitrogen. Approximately 2500 embryos and endosperm samples were dissected and pooled into one biological replicate, and three biological replicates were collected.

Small RNA library construction and sequencing

After total RNA of embryos and endosperm from each reciprocal cross achene was extracted by a TRIzol reagent kit (Invitrogen, Carlsbad, CA, USA), the small RNAs in a size range of 16–30 nt were enriched by 15% denaturing polyacrylamide gel electrophoresis (PAGE). Then, 3’ adapters were added, and the 36–44 nt RNAs were enriched. Next, 5’ adapters were ligated to the RNAs. The ligation products were reverse transcribed by PCR amplification, and the 140–160 bp PCR products were enriched to generate a cDNA library and sequenced using Illumina HiSeq™ 2500.

Small RNA sequencing analysis

Clean tags were obtained by removing the adaptor sequences, filtering out the low-quality tags and eliminating short lengths (< 14 nt). Then, all the clean reads were aligned with small RNA sequences in GenBank databases [32] and the Rfam 11.0 database [33] to identify and filter the ribosomal RNA (rRNA), transfer RNA (tRNA), small cytoplasmic RNA (scRNA), small nuclear RNA (snRNA) and small nucleolar RNA (snoRNA). After removing noncoding RNAs, the clean small RNA sequences were mapped to the F. × ananassa ‘Camarosa’ genome [34] to obtain unique reads with abundance and location information based on the genome using the SOAP 2.0 program (http://soap.genomics.org.cn/). The unique RNA sequences that perfectly matched the genome were subjected to subsequent analysis. Sequence reads that mapped to exons and introns of mRNA were removed to avoid mRNA degradation products, and the tags that mapped to repeat sequences were also removed.

Identification of known and novel miRNAs and their targets

To identify the known miRNAs, the unique RNA sequences that perfectly matched the F. × ananassa ‘Camarosa’ genome [34] were searched against the mature or precursor miRNAs of cultivated strawberry and mature plant miRNA sequences in the miRBase 22 database [35] by using the sRNAanno database [43]. The small RNA sequences that perfectly matched known miRNAs were considered to be known miRNAs. Potential miRNA precursor sequences were obtained by aligning the reads to the strawberry reference genome, and their secondary structure were predicted by using the RNAfold [86]. The results of secondary structure were analysed to verify characteristic miRNAs based on the previous criteria [28]. Sequences that did not find matches in the miRBase were identified as novel miRNAs by sRNAanno [43], which uses a set of well-established criteria, as previously reported [87]. Secondary structures of miRNAs are shown by RNAfold [86].

Based on the sequences of the known miRNAs and novel miRNAs, the software PatMatch (Version 1.2) [88] was used to predict target genes. The following criteria were applied: (1) no more than four mismatches between the sRNA and the target; (2) no more than two adjacent mismatches in the miRNA/target duplex; and (3) no mismatches in positions 10–11 of the miRNA/target duplex. (4) The minimum free energy (MFE) of the miRNA/target duplex should be ≥ 74% of the MFE of the miRNA bound to its perfect complement.

Identification of differentially expressed miRNAs in the embryo and endosperm of strawberry

To screen differentially expressed miRNAs during embryo and endosperm development in strawberry, the expression level of miRNAs in each sample was calculated and normalized to transcripts per million (TPM) using miRNA counts. The significant differentially expressed miRNAs were determined by a fold change ≥ 2 and P value < 0.05 for each comparison. The heatmap of miRNA expression was generated using TBtools [47].

qRT‒PCR analysis of differentially expressed miRNAs and their target genes

For quantification of mature miRNAs, we used the tailing-reaction qRT‒PCR method. Tailing-reaction RT‒PCR was performed using a miRNA 1st strand cDNA synthesis kit (Accurate Biotechnology, China) with general RT tailing-reaction primers. TB Green® Premix Ex Taq™ (TaKaRa, Japan) and tailing-reaction 3’ primer (Accurate Biotechnology, China) were used for qPCR on a QuantStudio 5 system (ABI, USA). For expression analysis of target genes, qRT‒PCR was performed using a QuantStudio 5 system (ABI, USA), with a PrimeScript™ RT reagent Kit (TaKaRa, Japan) for reverse transcription and TB Green® Premix Ex Taq™ (TaKaRa, Japan) for qPCR. Three independent biological replicates were used, and 26S rRNA and EF-1α were used as the reference genes for the analysis of the expression of miRNAs and target genes, respectively. The relative expression levels of miRNAs and target genes were calculated by using the 2−ΔΔCt method. All primers used in this experiment are listed in Table S5.

Identification and analyses of putative PHAS loci and miRNA triggers

To investigate phasiRNA biogenesis and regulation in embryos and endosperm, a genome-wide search for PHAS loci was conducted based on small RNA reads by using the sRNAanno database [43]. Loci with a P value < 0.001 were considered to be reliable PHAS loci. These PHAS loci were refined for putative annotation by aligning their sequences to the NCBI Nucleotide Collection (nr/nt) database (https://www.ncbi.nlm.nih.gov/) and the PlantRep database [89]. A potential miRNA trigger of PHAS loci was also predicted using the sRNAanno database [43].

Availability of data and materials

The raw sRNA sequencing data have been submitted to the NCBI SRA database under the accession numbers SAMN28078473, SAMN28078474, SAMN28078475 and SAMN28078476.

Abbreviations

ABA:

Abscisic acid

AFB2:

Auxin signalling F-Box 2

ARF:

AUXIN RESPONSIVE FACTOR

ASA1:

Anthranilate synthase alpha subunit 1

BS:

‘Benihoppe’ as the female, pollinated with ‘Sweet Charlie’

CM3:

Chorismate mutase 3

DAP:

Days after pollination

EDM2:

ENHANCED DOWNY MILDEW 2

FT1:

Flowering locus T1

hc-siRNAs:

Heterochromatic small interfering RNAs

IPA:

Indole-3-pyruvic acid

LEC2:

B3 domain-containing transcription factor LEC2

MFE:

Minimum free energy

miRNAs:

MicroRNAs

MYB:

V-myb avian myeloblastosis viral oncogene homologue

NAC:

(NAM, ATAF and CUC) proteins

NB-LRR:

Nucleotide-binding leucine-rich repeat receptors

PAGE:

Polyacrylamide gel electrophoresis

PHAS loci:

PhasiRNA-producing loci

phasiRNAs:

Phased small interfering RNAs

PPR:

PENTATRICOPEPTIDE REPEAT

qRT‒PCR:

Quantitative reverse transcription PCR

RdDM:

RNA-directed DNA methylation

RDR6:

RNA-DEPENDENT RNA POLYMERASE6

rRNAs:

Ribosomal RNAs

SB:

‘Sweet Charlie’ as the female, pollinated with ‘Benihoppe’

scRNAs:

Small cytoplasmic RNAs

SEM:

Scanning electron microscopy

siRNAs:

Short interfering RNAs

snoRNAs:

Small nucleolar RNAs

snRNAs:

Small nuclear RNAs

sRNAs:

Small RNAs

TAR2:

Tryptophan aminotransferase-related 2

TPM:

Transcripts per million

tRNAs:

Transfer RNAs

Trp:

Tryptophan

References

  1. Lei JJ, Jiang S, Ma RY, Xue L, Zhao J, Dai HP. Current status of strawberry industry in China. Acta Hortic. 2021;1309:349–52.

    Article  Google Scholar 

  2. Nitsch JP. Growth and morphogenesis of the strawberry as related to auxin. Am J Bot. 1950;37:211–5.

    Article  CAS  Google Scholar 

  3. Veluthambi K, Rhee JK, Mizrahi Y, Poovaiah BW. Correlation between lack of receptacle growth in response to auxin and accumulation of a specific polypeptide in a strawberry (Fragaria ananassa Duch.) variant genotype. Plant Cell Physiol. 1985;26:317–24.

    CAS  Google Scholar 

  4. Nitsch JP. Free auxins and free tryptophane in the strawberry. Plant Physiol. 1955;30:33–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Dreher T, Poovaiah BW. Changes in auxin content during development in strawberry fruits. J Plant Growth Regul. 1982;1:267–76.

    Google Scholar 

  6. Fait A, Hanhineva K, Beleggia R, Dai N, Rogachev I, Nikiforova VJ, Fernie AR, Aharoni A. Reconfiguration of the achene and receptacle metabolic networks during strawberry fruit development. Plant Physiol. 2008;148:730–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Manning K. Soft fruit. In: Seymour GB, Taylor JE, Tucker GA, editors. Biochemistry of fruit ripening. London: Chapman & Hall; 1993. p. 347–78.

    Chapter  Google Scholar 

  8. Symons GM, Chua YJ, Ross JJ, Quittenden LJ, Davies NW, Reid JB. Hormonal changes during non-climacteric ripening in strawberry. J Exp Bot. 2012;63:4741–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Li T, Dai Z, Zeng B, Li J, Ouyang J, Kang L, Wang W, Jia W. Autocatalytic biosynthesis of abscisic acid and its synergistic action with auxin to regulate strawberry fruit ripening. Hortic Res. 2022.https://doi.org/10.1093/hr/uhab076.

  10. Kang C, Darwish O, Geretz A, Shahan R, Alkharouf N, Liu Z. Genome-scale transcriptomic insights into early-stage fruit development in woodland strawberry Fragaria vesca. Plant Cell. 2013;25:1960–78.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Chaudhury AM, Koltunow A, Payne T, Luo M, Tucker MR, Dennis ES, Peacock WJ. Control of early seed development. Annu Rev Cell Dev Biol. 2001;17:677–99.

    Article  CAS  PubMed  Google Scholar 

  12. An L, Tao Y, Chen H, He M, Xiao F, Li G, Ding Y, Liu Z. Embryo-endosperm interaction and its agronomic relevance to rice quality. Front Plant Sci. 2020;11:587641.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Zheng X, Li Q, Li C, An D, Xiao Q, Wang W, Wu Y. Intra-kernel reallocation of proteins in maize depends on VP1-mediated scutellum development and nutrient assimilation. Plant Cell. 2019;31:2613–35.

    CAS  PubMed  PubMed Central  Google Scholar 

  14. Zeeman SC. “Carbohydrate metabolism” in Biochemistry and molecular biology of plants. In: Buchanan BB, Gruissem W, Jones RL, editors. John Wiley & Sons. Chichester: Academic; 2015. p. 567–609.

    Google Scholar 

  15. Doll NM, Just J, Brunaud V, Caïus J, Grimault A, Depège-Fargeix N, Esteban E, Pasha A, Provart NJ, Ingram GC, Rogowsky PM, Widiez T. Transcriptomics at maize embryo/endosperm interfaces identifies a transcriptionally distinct endosperm subdomain adjacent to the embryo scutellum. Plant Cell. 2020;32:833–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Chen C, Zeng Z, Liu Z, Xia R. Small RNAs, emerging regulators critical for the development of horticultural traits. Hortic Res. 2018;5:63.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Chen X. Small RNAs and their roles in plant development. Annu Rev Cell Dev Biol. 2009;25:21–44.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  18. Axtell MJ. Classification and comparison of small RNAs from plants. Annu Rev Plant Biol. 2013;64:137–59.

    Article  CAS  PubMed  Google Scholar 

  19. Ye R, Wang W, Iki T, Liu C, Wu Y, Ishikawa M, Zhou X, Qi Y. Cytoplasmic assembly and selective nuclear import of Arabidopsis ARGONAUTE4/siRNA complexes. Mol Cell. 2012;46:859–70.

    Article  CAS  PubMed  Google Scholar 

  20. Axtell MJ, Jan C, Rajagopalan R, Bartel DP. A two-hit trigger for siRNA biogenesis in plants. Cell. 2006;127:565–77.

    Article  CAS  PubMed  Google Scholar 

  21. Adenot X, Elmayan T, Lauressergues D, Boutet S, Bouché N, Gasciolli V, Vaucheret H. DRB4-dependent TAS3 trans-acting siRNAs control leaf morphology through AGO7. Curr Biol. 2006;16:927–32.

    Article  CAS  PubMed  Google Scholar 

  22. Zhou C, Han L, Fu C, Wen J, Cheng X, Nakashima J, Ma J, Tang Y, Tan Y, Tadege M, Mysore KS, Xia G, Wang ZY. The trans-acting short interfering RNA3 pathway and NO APICAL MERISTEM antagonistically regulate leaf margin development and lateral organ separation, as revealed by analysis of an argonaute7/lobed leaflet1 mutant in Medicago truncatula. Plant Cell. 2013;25:4845–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. López-Márquez D, Del-Espino Á, López-Pagán N, Rodríguez-Negrete EA, Rubio-Somoza I, Ruiz-Albert J, Bejarano ER, Beuzón CR. MiRNA and phasiRNAs-mediated regulation of TIR-NBS-LRR defense genes in Arabidopsis thaliana. bioRxiv. 2020. https://doi.org/10.1101/2020.03.02.972620.

  24. Zhai J, Jeong DH, de Paoli E, Park S, Rosen BD, Li Y, González AJ, Yan Z, Kitto SL, Grusak MA, Jackson SA, Stacey G, Cook DR, Green PJ, Sherrier DJ, Meyers BC. MicroRNAs as master regulators of the plant NB-LRR defense gene family via the production of phased, trans-acting siRNAs. Gene Dev. 2011;25:2540–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Xia R, Zhu H, An YQ, Beers EP, Liu Z. Apple miRNAs and tasiRNAs with novel regulatory networks. Genome Biol. 2012;13:R47.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Ma W, Chen C, Liu Y, Zeng M, Meyers BC, Li J, Xia R. Coupling of microRNA-directed phased small interfering RNA generation from long noncoding genes with alternative splicing and alternative polyadenylation in small RNA-mediated gene silencing. New Phytol. 2018;217:1535–50.

    Article  CAS  PubMed  Google Scholar 

  27. Feng L, Xia R, Liu Y. Comprehensive characterization of miRNA and PHAS loci in the diploid strawberry (Fragaria vesca) genome. Hortic Plant J. 2019;5:255–67.

    Article  Google Scholar 

  28. Guo J, Wang Q, Liu L, Ren S, Li S, Liao P, Zhao Z, Lu C, Jiang B, Sunkar R, Zheng Y. Analysis of microRNAs, phased small interfering RNAs and their potential targets in Rosa rugosa Thunb. BMC Genomics. 2019;19:983.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  29. Li H, Mao W, Liu W, Dai H, Liu Y, Ma Y, Zhang Z. Deep sequencing discovery of novel and conserved microRNAs in wild type and a white-flesh mutant strawberry. Planta. 2013;238:695–713.

    Article  CAS  PubMed  Google Scholar 

  30. Li D, Mou W, Xia R, Li L, Zawora C, Ying T, Mao L, Liu Z, Luo Z. Integrated analysis of high-throughput sequencing data shows abscisic acid-responsive genes and miRNAs in strawberry receptacle fruit ripening. Hortic Res. 2019;6:26.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  31. Xu X, Ma X, Lei H, Yin L, Shi X, Song H. MicroRNAs play an important role in the regulation of strawberry fruit senescence in low temperature. Postharvest Biol Tec. 2015;108:39–47.

    Article  CAS  Google Scholar 

  32. Benson DA, Cavanaugh M, Clark K, Karsch-Mizrachi I, Lipman DJ, Ostell J, Sayers EW. GenBank. Nucleic Acids Res. 2013;41:D36-42.

    Article  CAS  PubMed  Google Scholar 

  33. Burge SW, Daub J, Eberhardt R, Tate J, Barquist L, Nawrocki EP, Eddy SR, Gardner PP, Bateman A. Rfam 11.0: 10 years of RNA families. Nucleic Acids Res. 2013;41:D226-32.

    Article  CAS  PubMed  Google Scholar 

  34. Edger PP, Poorten TJ, VanBuren R, Hardigan MA, Colle M, McKain MR, Smith RD, Teresi SJ, Nelson ADL, Wai CM, Alger EI, Bird KA, Yocca AE, Pumplin N, Ou S, Ben-Zvi G, Brodt A, Baruch K, Swale T, Shiue L, Acharya CB, Cole GS, Mower JP, Childs KL, Jiang N, Lyons E, Freeling M, Puzey JR, Knapp SJ. Origin and evolution of the octoploid strawberry genome. Nat Genet. 2019;51:541–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Kozomara A, Birgaoanu M, Griffiths-Jones S. miRBase: from microRNA sequences to function. Nucleic Acids Res. 2019;47:D155–62.

    Article  CAS  PubMed  Google Scholar 

  36. Xu W, Cui Q, Li F, Liu A. Transcriptome-wide identification and characterization of microRNAs from castor bean (Ricinus communis L.). PLoS One. 2013;8:e69995.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. D’Ario M, Griffiths-Jones S, Kim M. Small RNAs: big impact on plant development. Trends Plant Sci. 2017;22:1056–68.

    Article  PubMed  CAS  Google Scholar 

  38. Xia R, Ye S, Liu Z, Meyers BC, Liu Z. Novel and recently evolved microRNA clusters regulate expansive F-box gene networks through phased small interfering RNAs in wild diploid strawberry. Plant Physiol. 2015;169:594–610.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  39. Schenk PW, Snaar-Jagalska BE. Signal perception and transduction: the role of protein kinases. BBA-Mol Cell Res. 1999;1449:1–24.

    CAS  Google Scholar 

  40. Wu MF, Tian Q, Reed JW. Arabidopsis microRNA167 controls patterns of ARF6 and ARF8 expression, and regulates both female and male reproduction. Development. 2006;133:4211–8.

    Article  CAS  PubMed  Google Scholar 

  41. José Ripoll J, Bailey LJ, Mai QA, Wu SL, Hon CT, Chapman EJ, Ditta GS, Estelle M, Yanofsky MF. microRNA regulation of fruit growth. Nat Plants. 2015;1:15036.

    Article  PubMed  CAS  Google Scholar 

  42. Yao JL, Xu J, Cornille A, Tomes S, Karunairetnam S, Luo Z, Bassett H, Whitworth C, Rees-George J, Ranatunga C, Snirc A, Crowhurst R, de Silva N, Warren B, Deng C, Kumar S, Chagné D, Bus VG, Volz RK, Rikkerink EH, Gardiner SE, Giraud T, MacDiarmid R, Gleave AP. A microRNA allele that emerged prior to apple domestication may underlie fruit size evolution. Plant J. 2015;84:417–27.

    Article  CAS  PubMed  Google Scholar 

  43. Chen C, Li J, Feng J, Liu B, Feng L, Yu X, Li G, Zhai J, Meyers BC, Xia R. sRNAanno—a database repository of uniformly annotated small RNAs in plants. Hortic Res. 2021;8:45.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Fei Q, Rui X, Meyers BC. Phased, secondary, small interfering RNAs in posttranscriptional regulatory networks. Plant Cell. 2013;25:2400–15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Liu Y, Teng C, Xia R, Meyers BC. PhasiRNAs in plants: their biogenesis, genic sources, and roles in stress responses, development, and reproduction. Plant Cell. 2020;32:3059–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Herrmann KM. The shikimate pathway as an entry to aromatic secondary metabolism. Plant Physiol. 1995;107:7–12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Chen C, Chen H, Zhang Y, Thomas HR, Frank MH, He Y, Xia R. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol Plant. 2020;13:1194–202.

    Article  CAS  PubMed  Google Scholar 

  48. Estrada-Johnson E, Csukasi F, Pizarro CM, Vallarino JG, Kiryakova Y, Vioque A, Brumos J, Medina-Escobar N, Botella MA, Alonso JM, Fernie AR, Sánchez-Sevilla JF, Osorio S, Valpuesta V. Transcriptomic analysis in strawberry fruits reveals active auxin biosynthesis and signaling in the ripe receptacle. Front Plant Sci. 2017;8:889.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Huang H, Long J, Zheng L, Li Y, Hu Y, Yu G, Liu H, Liu Y, Huang Z, Zhang J, Huang Y. Identification and characterization of microRNA in maize endosperm response to exogenous sucrose using small RNA sequencing. Genomics. 2016;108:216–23.

    Article  CAS  PubMed  Google Scholar 

  50. Liu J, Guo X, Zhai T, Shu A, Zhao L, Liu Z, Zhang S. Genome-wide identification and characterization of microRNA responding to ABA and GA in maize embryos during seed germination. Plant Biol. 2020;22:949–57.

    Article  CAS  PubMed  Google Scholar 

  51. Cheng J, Niu Q, Zhang B, Chen K, Yang R, Zhu JK, Zhang Y, Lang Z. Downregulation of RdDM during strawberry fruit ripening. Genome Biol. 2018;19:212.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Xu M, Hu T, Zhao J, Park MY, Earley KW, Wu G, Yang L, Poethig RS. Developmental functions of miR156-regulated SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) genes in Arabidopsis thaliana. PLoS Genet. 2016;12:e1006263.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  53. Wu G, Park MY, Conway SR, Wang JW, Weigel D, Poethig RS. The sequential action of miR156 and miR172 regulates developmental timing in Arabidopsis. Cell. 2009;138:750–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Huang D, Koh C, Feurtado JA, Tsang EW, Cutler AJ. MicroRNAs and their putative targets in Brassica napus seed maturation. BMC Genomics. 2013;14:140.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Jung JH, Park CM. MIR166/165 genes exhibit dynamic expression patterns in regulating shoot apical meristem and floral development in Arabidopsis. Planta. 2007;225:1327–38.

    Article  CAS  PubMed  Google Scholar 

  56. Xue C, Yao JL, Qin MF, Zhang MY, Allan AC, Wang DF, Wu J. PbrmiR397a regulates lignification during stone cell development in pear fruit. Plant Biotechnol J. 2019;17:103–17.

    Article  CAS  PubMed  Google Scholar 

  57. Jodder J. miRNA-mediated regulation of auxin signaling pathway during plant development and stress responses. J Biosciences. 2020;45:91.

    Article  CAS  Google Scholar 

  58. Guilfoyle TJ, Hagen G. Auxin response factors. J Plant Growth Regul. 2001;20:281–91.

    Article  CAS  Google Scholar 

  59. Okushima Y, Overvoorde PJ, Arima K, Alonso JM, Chan A, Chang C, Ecker JR, Hughes B, Lui A, Nguyen D, Onodera C, Quach H, Smith A, Yu G, Theologis A. Functional genomic analysis of the AUXIN RESPONSE FACTOR gene family members in Arabidopsis thaliana: unique and overlapping functions of ARF7 and ARF19. Plant Cell. 2005;17:444–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Chandler JW. Auxin response factors. Plant Cell Environ. 2016;39:1014–28.

    Article  CAS  PubMed  Google Scholar 

  61. Roosjen M, Paque S, Weijers D. Auxin response factors: output control in auxin biology. J Exp Bot. 2018;69:179–88.

    Article  CAS  PubMed  Google Scholar 

  62. Kovalenko TF, Patrushev LI. Pseudogenes as functionally significant elements of the genome. Biochemistry Moscow. 2018;83:1332–49.

    Article  CAS  PubMed  Google Scholar 

  63. Ma Y, Liu S, Gao J, Chen C, Zhang X, Yuan H, Chen Z, Yin X, Sun C, Mao Y, Zhou F, Shao Y, Liu Q, Xu J, Cheng L, Yu D, Li P, Yu J. Genome-wide analysis of pseudogenes reveals HBBP1’s human-specific essentiality in erythropoiesis and implication in β-thalassemia. Dev Cell. 2021;56:478–93.

    Article  CAS  PubMed  Google Scholar 

  64. Li D, Liu Z, Gao L, Wang L, Gao M, Jiao Z, Qiao H, Yang J, Chen M, Yao L, Kan Y. Genome-wide identification and characterization of microRNAs in developing grains of Zea mays L. PLoS One. 2016;11:e0153168.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  65. Arikit S, Xia R, Kakrana A, Huang K, Zhai J, Yan Z, Valdés-López O, Prince S, Musket TA, Nguyen HT, Stacey C, Meyers BC. An atlas of soybean small RNAs identifies phased siRNAs from hundreds of coding genes. Plant Cell. 2014;26:4584–601.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Wu F, Chen Y, Tian X, Zhu X, Jin W. Genome-wide identification and characterization of phased small interfering RNA genes in response to Botrytis cinerea infection in Solanum lycopersicum. Sci Rep-UK. 2017;7:3019.

    Article  CAS  Google Scholar 

  67. Chen K, Liu L, Zhang X, Yuan Y, Ren S, Guo J, Wang Q, Liao P, Li S, Gui X, Li YF, Zheng Y. Phased secondary small interfering RNAs in Panax notoginseng. BMC Genomics. 2018;19:41.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  68. Zhang R, Huang S, Li S, Song G, Li Y, Li W, Li J, Gao J, Gu T, Li D, Zhang S, Li G. Evolution of PHAS loci in the young spike of allohexaploid wheat. BMC Genomics. 2020;21:200.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Fang YN, Yang XM, Jiang N, Wu XM, Guo WW. Genome-wide identification and expression profiles of phased siRNAs in a male-sterile somatic cybrid of pummelo (Citrus grandis). Tree Genet Genomes. 2020;16:46.

    Article  Google Scholar 

  70. Swetha C, Narjala A, Pandit A, Tirumalai V, Shivaprasad PV. Degradome comparison between wild and cultivated rice identifies differential targeting by miRNAs. BMC Genomics. 2022;23:53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  71. Chen HM, Li YH, Wu SH. Bioinformatic prediction and experimental validation of a microRNA-directed tandem trans-acting siRNA cascade in Arabidopsis. Proc Natl Acad Sci U S A. 2007;104:3318–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  72. Howell MD, Fahlgren N, Chapman EJ, Cumbie JS, Sullivan CM, Givan SA, Kasschau KD, Carrington JC. Genome-wide analysis of the RNA-DEPENDENT RNA POLYMERASE6/ DICER-LIKE4 pathway in Arabidopsis reveals dependency on miRNA- and tasiRNA-directed targeting. Plant Cell. 2007;19:926–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Xia R. MicroRNAs and trans-acting siRNA pathways in apple (Malus x domestica Borkh.) and peach (Prunus persica). Dissertation, Virginia Polytechnic Institute and State University, Blacksburg, Virginia, USA. 2013.

  74. Bartel B. Auxin biosynthesis. Annu Rev Plant Biol. 1997;48:51–66.

    Article  CAS  Google Scholar 

  75. Stepanova AN, Robertson-Hoyt J, Yun J, Benavente LM, Xie DY, Doležal K, Schlereth A, Jürgens G, Alonso JM. TAA1-mediated auxin biosynthesis is essential for hormone crosstalk and plant development. Cell. 2008;133:177–91.

    Article  CAS  PubMed  Google Scholar 

  76. Mashiguchi K, Tanaka K, Sakai T, Sugawara S, Kawaide H, Natsume M, Hanada A, Yaeno T, Shirasu K, Yao H, McSteend P, Zhao Y, Hayashif K, Kamiyaa Y, Kasaharaet H. The main auxin biosynthesis pathway in Arabidopsis. Proc Natl Acad Sci U S A. 2011;108:18512–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  77. Ma W, Li J, Qu B, Xue H, Zhao X, Li B, Fu X, Tong Y. Auxin biosynthetic gene TAR2 is involved in low nitrogen-mediated reprogramming of root architecture in Arabidopsis. Plant J. 2014;78:70–9.

    Article  CAS  PubMed  Google Scholar 

  78. Zhao Y. Essential roles of local auxin biosynthesis in plant development and in adaptation to environmental changes. Annu Rev Plant Biol. 2018;69:417–35.

    Article  CAS  PubMed  Google Scholar 

  79. Xie Q, Frugis G, Colgan D, Chua N. Arabidopsis NAC1 transduces auxin signal downstream of TIR1 to promote lateral root development. Gene Dev. 2000;14:3024–36.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  80. Mallory AC, Dugas DV, Bartel DP, Bartel B. MicroRNA regulation of NAC-domain targets is required for proper formation and separation of adjacent embryonic, vegetative, and floral organs. Curr Biol. 2004;14:1035–46.

    Article  CAS  PubMed  Google Scholar 

  81. Parry G, Calderon-Villalobos LI, Prigge M, Peret B, Dharmasiri S, Itoh H, Lechner E, Gray WM, Bennett M, Estelle M. Complex regulation of the TIR1/AFB family of auxin receptors. Proc Natl Acad Sci U S A. 2009;106:22540–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  82. Hu Z, Keçeli MA, Piisilä M, Li J, Survila M, Heino P, Brader G, Palva ET, Li J. F-box protein AFB4 plays a crucial role in plant growth, development and innate immunity. Cell Res. 2012;22:777–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  83. Cui L, Zhang T, Li J, Lou Q, Chen J. Cloning and expression analysis of Cs-TIR1/AFB2: the fruit development-related genes of cucumber (Cucumis sativus L.). Acta Physiol Plant. 2014;36:139–49.

    Article  CAS  Google Scholar 

  84. Lakehal A, Chaabouni S, Cavel E, Hir RL, Ranjan A, Raneshan Z, Novák O, Pacurar DI, Perrone I, Jobert F, Gutierrez L, Bakò L, Bellini C. A molecular framework for the control of adventitious rooting by TIR1/AFB2-Aux/IAA-dependent auxin signaling in Arabidopsis. Mol Plant. 2019;12:1499–514.

    Article  CAS  PubMed  Google Scholar 

  85. Liu Y, Jing X, Zhang H, Xiong J, Qiao Y. Identification of imprinted genes based on homology: an example of Fragaria vesca. Genes. 2021;12:380.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  86. Hofacker IL. Vienna RNA secondary structure server. Nucleic Acids Res. 2003;31:3429–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  87. Axtell MJ, Meyers BC. Revisiting criteria for plant microRNA annotation in the era of big data. Plant Cell. 2018;30:272–84.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  88. Yan T, Yoo D, Berardini TZ, Mueller LA, Weems DC, Weng S, Cherry JM, Rhee SY. PatMatch: a program for finding patterns in peptide and nucleotide sequences. Nucleic Acids Res. 2005;33:W262–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  89. Luo X, Chen S, ZhangY. PlantRep: a resource of plant repeats. Research Square. 2021. https://doi.org/10.21203/rs.3.rs-498874/v1.

Download references

Acknowledgements

We thank Dr. Yuehua Ma (Central Laboratory of College of Horticulture, Nanjing Agricultural University) for assistance in using quantitative real-time PCR.

Funding

This research was funded by the National Natural Science Foundation of China, grant number 31872056, and the Fundamental Research Funds for the Central Universities, grant number KYZZ2022004.

Author information

Authors and Affiliations

Authors

Contributions

YQ, XJ and HZ designed the experiment. HZ, XJ, QA and XH collected the samples and prepared sRNA samples for sequencing. XJ and XH carried out the qRT‒PCR. XJ analysed the data and wrote the manuscript. YQ and XJ reviewed and proofread the manuscript. All authors have read and agreed to the published version of the manuscript.

Corresponding author

Correspondence to Yushan Qiao.

Ethics declarations

Ethics approval and consent to participate

All relevant permissions were obtained for the collection of strawberry lines. The methods carried out in the manuscript were in accordance with the relevant local, national and international guidelines and regulations.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

Figure S1. Statistical analysis of sequencing reads in twelve libraries in strawberry reciprocally crossed embryo and endosperm. BS and SB represent ‘Benihoppe’ (♀) × ‘Sweet Charlie’ (♂) and ‘Sweet Charlie’ (♀) × ‘Benihoppe’ (♂), respectively. em: embryo; en: endosperm; rRNA: ribosomal RNA; tRNA: transfer RNA; snoRNA: small nucleolar RNA; snRNA: small nuclear RNA. For each tissue, three biological replicates (em1, em2, and em3; en1, en2, and en3) were performed.

Additional file 2:

Figure S2. The abundant of miRNA families in strawberry reciprocally crossed embryo and endosperm.

Additional file 3: Figure S3.

Expression patterns of 24nt PHAS loci in strawberry embryos and endosperm. A. Heatmap of all 24nt PHAS loci expression patterns in embryos and endosperm. B. sRNA abundance and phasing score at the zinc transporter protein locus. C. sRNA abundance and phasing score at the MYB locus.

Additional file 4: Figure S4.

The miRNA-triggered 24nt phasiRNAs in strawberry. A. sRNA abundance and phasing score are viewed at ncRNA is targeted by miRN29. B. Novel PHAS locus CSC1-like targeted by miR11285. C. ncRNA targeted by miRN14. D. The novel PHAS locus targeted by miR395 and miR11288.

Additional file 5: Figure S5.

Some miR159 family members was predicted targeted ARF23 that showed an endosperm-specific expression in strawberry.

Additional file 6: Table S1.

Distribution of the sRNA sequencing in the 12 libraries.

Additional file 7: Table S2.

Known and novel miRNAs in the developing embryo and endosperm.

Additional file 8: Table S3.

Prediction and analysis miRNA target genes.

Additional file 9: Table S4.

The basic information and annotation of the predicted PHAS loci in the strawberry.

Additional file 10: Table S5.

The primers for qRT-PCR.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Jing, X., Zhang, H., Huai, X. et al. Identification and characterization of miRNAs and PHAS loci related to the early development of the embryo and endosperm in Fragaria × ananassa. BMC Genomics 23, 638 (2022). https://doi.org/10.1186/s12864-022-08864-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-022-08864-3

Keywords