Next Article in Journal
Resveratrol Inhibition of the WNT/β-Catenin Pathway following Discogenic Low Back Pain
Previous Article in Journal
QSAR, Docking, and Molecular Dynamics Simulation Studies of Sigmacidins as Antimicrobials against Streptococci
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 8
10.3390/ijms23084090
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of Sex Differentiation-Related microRNAs in Spinach Female and Male Flower

by
Ning Li
,
Yueyuan Wang
,
Jiwen Wang
,
Wanqing Zhang
,
Ziwei Meng
,
Yuanshen Wang
,
Yulan Zhang
,
Shufen Li
,
Wujun Gao
and
Chuanliang Deng
*
College of Life Sciences, Henan Normal University, Xinxiang 453007, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(8), 4090; https://doi.org/10.3390/ijms23084090
Submission received: 3 March 2022 / Revised: 29 March 2022 / Accepted: 4 April 2022 / Published: 7 April 2022
(This article belongs to the Section Molecular Plant Sciences)
Browse Figures
Versions Notes

Abstract

:
Sex determination and differentiation is an important biological process for unisexual flower development. Spinach is a model plant to study the mechanism of sex determination and differentiation of dioecious plant. Till now, little is known about spinach sex determination and differentiation mechanism. MicroRNAs are key factors in flower development. Herein, small RNA sequencing was performed to explore the roles of microRNAs in spinach sex determination and differentiation. As a result, 92 known and 3402 novel microRNAs were identified in 18 spinach female and male flower samples. 74 differentially expressed microRNAs were identified between female and male flowers, including 20 female-biased and 48 male-biased expression microRNAs. Target prediction identified 22 sex-biased microRNA-target pairs, which may be involved in spinach sex determination or differentiation. Among the differentially expressed microRNAs between FNS and M03, 55 microRNAs were found to reside in sex chromosome; one of them, sol-miR2550n, was functionally studied via genetic transformation. Silencing of sol-miR2550n resulted in abnormal anther while overexpression of sol-miR2550n induced early flowering, indicating sol-miR2550n was a male-promoting factor and validating the reliability of our small RNA sequencing data. Conclusively, this work can supply valuable information for exploring spinach sex determination and differentiation and provide a new insight in studying unisexual flower development.

1. Introduction

Sex is a special character of biology. Most animals are unisexual and the mechanism of sex determination mammalian is already clearer. In plant kingdom, only 6% of angiosperms are dioecious plants and their sex determining mechanism are not well studied [1]. The gender origins of dioecious plant are not uniform, which enhances the difficulty to uncover plant sex determining mechanism. At present, there are two main ideas about the mechanism of sex determination in dioecious plants: one is the two-gene mutation model, and the other is the single-gene mutation model. For two-gene mutation model, the appearance of dioecious plants may be due to the mutation of two genes, M gene, which controls stamen fertility, and F gene, which regulates pistil development [2,3]. For single-gene mutation model, the emergence of dioecious plants may be caused by the two independent mutations of the same gene Q: one is loss-of-function mutation resulting in male sterility, and the other is gain-of-function mutation inducing female sterility [4,5,6]. From these two models, it can be seen that the mutation of flower development-related genes is the most important factor leading to the generation of plant sex. To date, the sex-determining genes or candidates have been identified in a few angiosperms such as Vitis vinifera [7,8], Diospyros lotus [9], Populus trichocarpa [10,11,12], Ficus carica [13], Date palm [14], Fragaria octoploids [15], Actinidia chinensis [16] and Asparagus officinalls [17]. The function of the sex-determining genes was studied only in few plants and the mechanism of sex determination remains largely elusive. Sex differentiation is an important process along with sex determination. Researches on sex differentiation may supply clues to uncover sex determination. Many genes related to sex differentiation have been identified in different plants such as in Silene latifolia [18,19], Carica papaya [20].
MicroRNAs (miRNAs) are a class of endogenous single-stranded non-coding small RNA molecules with a length of about 21–24 nt. The initial transcription product pri-miRNA is spliced into a single strand RNA precursor (pre-miRNA) with hairpin structure, and then formed the mature miRNA with regulatory function. MiRNAs negatively regulate their target messenger RNAs (mRNA) through mRNA degradation, translation inhibition or chromatin modification [21,22,23]. MiRNAs play important regulatory roles in plant morphogenesis, nutrient balance, biological and abiotic stress [24,25,26,27]. More and more studies have shown that miRNAs played important roles in plant flower development. Flowering is an important process in plant growth cycle and miRNAs are involved in flower formation processes. In Arabidopsis, miR164c regulated two NAC transcription factors CUP-SHAPED COTYLEDON1 (CUC1) and CUC2 to control the boundary foundation of floral organs [28]. MiR172 can regulate the formation of sepal and petal primordium via inhibiting the translation of floral homeotic gene APETALA2 [29,30]. MiR159, miR167 and miR319 can form a regulation module to influence flower development by targeting transcription factors MYB33, TCP4 and Auxin Response Factor6/8 [31]. In tomato, sly-miR160a was involved in flower development via regulating auxin response factor SlARF10A [32]. The flowers of corn, a hermaphrodite plant, are bisexual in the early development stage but become unisexual in the late stage as the development of pistil or stamen is blocked. It has been reported that zma-miR172E was a member of miR172 family and an insertional mutation in the promoter region of its coding gene resulted in a corn recessive mutant tasseleseed4, in which carpel developed in tassel with no stamen [33]. The sex determining gene of a dioecious plant, persimmon (Diospyros lotus) was identified as a miRNA, OGI, targeting a homeodomain transcription factor MeGI to determine plant sex [9,34]. These studies fully demonstrate that miRNAs participate in flower development, even in the formation of unisexual flowers and moreover miRNA can serve as sex determining gene.
Spinacia oleracea (2n = 12), a dioecious plant, has a pair of homomorphic sex chromosomes (XY) and is at the early stage of sex chromosome evolution [35]. Many sex-linked, sex chromosome-specific and sex-determining-gene-linked molecular markers were developed in spinach, such as T11A, V20A, 45 s rDNA, SP_0018, SpoX and so on [36,37,38,39]. These molecular markers are helpful for locating the accurate sex-determining region (SDR) and further identifying sex-determining genes. Qian et al. (2017) [40], Yu et al. (2021) [41] and Ma et al. (2022) [42] found the SDR locating on the sex chromosome using different method. Moreover, 166 sex-differentiation-related genes and 12 Y-specific genes were successively reported in spinach [40,43]. However, the function of these genes was still unclear and it is still controversial that spinach sex determination is regulated by single-gene mutation model or two-gene mutation model [43]. Spinach genome was firstly published in 2017 and constantly updated with the development of high throughput technology, which supply valuable information for functional genome study [41,42,44,45,46].
To verify whether miRNA can serve as sex determining gene in spinach, miRNAs resource should be firstly harvested. In 2017, spinach miRNAs were analyzed via in silico prediction in vegetative tissues of spinach [47]. Flower development-related genes are crucial for sex determination and differentiation, so we performed small RNA sequencing using 18 spinach flower samples at three early female and male flower development stages to identify some candidate miRNAs for spinach sex determination or differentiation. Virus-induced gene silencing (VIGS) technology is an efficient way to study gene function in plant [48]. In addition, this technology has been successfully applied in spinach [49,50]. Moreover, RNA interference technology and overexpression technology have been used to study miRNA function [51]. Hence, we performed VIGS and overexpression technology to analyze the function of sol-miR2550n, one of the candidate miRNAs, and to further validate the accuracy of our sequencing data.

2. Results

2.1. Small RNA Sequencing

The roles of miRNAs in flower development are not clear in spinach till now. To characterize miRNAs in spinach, we performed small RNA sequencing using 18 spinach flower samples and harvested 13,911,586 high quality reads (mean value of 18 samples data) (Table 1). After filtering, at least 91.37% of high quality reads were selected as small RNA clean tags, the primary length distribution of which is 24 nt (Figure 1). All small RNA clean tags were aligned in GeneBank database (Release 209.0), Rfam database (11.0), spinach genome (http://www.spinachbase.org/, accessed on 23 July 2021) and miRBase database (Release 21). Finally, 92 known miRNAs and 3402 novel miRNAs were identified in 18 spinach flower samples (Table 2).

2.2. Sex-Biased miRNAs of Spinach

According to the expression of miRNA, differentially expressed miRNAs (DE miRNAs) were analyzed among female and male flower samples. There were 431 DE miRNAs between FNS and M03, 451 DE miRNAs between FNB and M05, and 781 DE miRNAs between FYS and M10 (Figure 2a). Sex-biased genes, which exhibit significantly higher expression in one sex than in the other sex, always act downstream of sex-determining gene. Hence, identification of the sex-biased genes is helpful to uncover the sex-determination mechanism. Herein, 74 DE miRNAs (7 were known and 67 were novel) were identified between female and male flowers at three early developmental stages (Figure 2a; Table S1); and of them, 20 miRNAs displayed female-biased expression and 48 miRNAs displayed male-biased expression (Figure 2b). Moreover, there were 9 female-specific expression miRNAs and 17 male-specific expression miRNAs (Figure 2b). These sex-biased and sex-specific expression miRNAs may be involved in spinach sex determination or differentiation.

2.3. Target Genes of miRNAs

Target genes of miRNAs were analyzed according to our previous transcriptome data [52]. Totally, 20,460 genes were predicted as targets of 2855 miRNAs and 49,591 target sites were found in these target genes for miRNAs binding (Table 3). It can be seen that one miRNA can regulate more than one gene and one gene can be controlled by more than one miRNA. MiRNA and its target always display antagonistic expression profile. As previously reported [52], 1946 male-biased genes and 961 female-biased genes were found between female and male flowers at three early developmental stages. Herein, we found ten male-biased genes were targeted by five female-biased miRNAs and six female-biased genes were targeted by ten male-biased miRNAs, composing 22 miRNA-target pairs (Figure 3, Table S2). The 22 miRNA-target pairs may be involved in spinach sex determination or differentiation.

2.4. MiRNAs Residing in Sex Chromosome

In the published spinach genome data, chr4 is the sex chromosome where the sex-determining genes reside [43]. Blasting with spinach genome [44], we found 307 miRNAs (16 known and 291 novel) residing in Chr1, 275 miRNAs (6 known and 269 novel) residing in Chr2, 415 miRNAs (12 known and 403 novel) residing in Chr3, 496 miRNAs (8 known and 488 novel) residing in Chr4 (sex chromosome), 282 miRNAs (17 known and 265 novel) residing in Chr5, 240 miRNAs (2 known and 238 novel) residing in Chr6, and 1479 miRNAs (31 known and 1448 novel) residing in scaffolds (Figure 4). Among 496 miRNAs residing in sex chromosome, there were 55 DE miRNAs between FNS and M03, 56 DE miRNAs between FNB and M05, and 114 DE miRNAs between FYS and M10. However, no intersection was existed among these three data sets, i.e., no DE miRNA was identified between female and male flowers at three developmental stages. Considering that sex determination and differentiation occurs at the earlier flower development stage in spinach, so we further analyzed DE miRNAs between FNS and M03, the earliest stage of the three flower developmental stages. As a result, 14 differential expression genes were found as targets of ten DE miRNAs between FNS and M03 (Table 4).

2.5. Functional Analysis of Sol-miR172 and Sol-miR2550n

To uncover the potential function of these miRNAs residing in sex chromosome, we performed plant transformation experiments for the first time. sol-miR172 (miR172-y) is a DE miRNA between FNS and M03 and showed higher expression level in M03 stage; its homolog in Arabidopsis has been reported to regulate flower development [29], so we selected sol-miR172 as a positive control. sol-miR2550n (novel-m2550-5p) resided in sex chromosome and displayed higher expression level in M03 stage than in FNS stage (Table S3), so we selected it as a sex-determining candidate. Virus-induced-gene-silencing and heterologous-overexpression were performed for functional analysis of sol-miR172 and sol-miR2550n in spinach. Silencing of sol-miR172 or sol-miR2550n in spinach resulted in abnormal male flower (anther abortion) (Figure 5a,c,f). Overexpression of sol-miR172 or sol-miR2550n in Arabidopsis induced early flowering (Figure 5b,d,e,g,h). The phenotype of overexpression sol-miR172 in Arabidopsis was consistent with Aukerman and Sakai’s report [53]. These results verified the accuracy of the small RNA sequencing and analysis. Moreover, this work will supply abundant valuable information to explore the mechanism of spinach sex determination and differentiation.

3. Discussion

Sex determination and differentiation is an important question for monoecious and dioecious plant. In some species, the sex determining genes have been identified, such as asparagus, kiwifruit and persimmons. However, sex determining gene is still unclear in spinach; it may be transcription factor or noncoding RNA. Hence, we performed small RNA sequencing to search some clues.
It was the first time to sequencing small RNAs in spinach. Totally, 92 known miRNAs and 3402 novel miRNAs were identified in 18 spinach flower samples (Table 1), which enriched plant miRNA resource. To analyze the miRNAs related to spinach sex determination or differentiation, 74 differentially expressed miRNAs were screened out among female and male flower at three early developmental stages (Figure 2a), including 20 female-biased miRNAs, 9 female-specific miRNAs, 48 male-biased miRNAs and 17 male-specific expression miRNAs (Figure 2b). Moreover, 22 miRNA-target pairs were found through target prediction (Figure 3). These 74 sex-biased or sex-specific miRNAs can be served as candidates of spinach sex differentiation, but not sex determination as they were not reside in sex chromosome. Genes residing in sex chromosome, especially in sex-determining region, are important clues to explore spinach sex determination and differentiation. Hence, we analyzed the miRNAs residing in sex chromosome; and then we found 496 miRNAs in sex chromosome, but none was sex-biased miRNA (Figure 4). Considering that sex determinant plays a role in early spinach flower development stage, so we screened out ten DE miRNAs between FNS and M03, two earliest stages among samples (Table 4). 14 targets with opposite expression trend to these ten DE miRNAs were also identified (Table 4). The function of these ten DE miRNAs residing in sex chromosome will be studied in our future work.
Herein, one of the DE miRNAs residing in sex chromosome, sol-miR2550n, were firstly studied for its function. Through VIGS and overexpression method, we found that silencing of sol-miR2550n in spinach induced abnormal male flower (anther abortion) and overexpression of sol-miR2550n in Arabidopsis resulted in early flowering (Figure 5). sol-miR2550n showed higher expression level in male flower (M03) than in female flower (FNS) (Table S3), hence down-regulation of its expression influenced the male flower development. However, up-regulation of its expression didn’t affect flower structure but induced early flowering. In spinach, the flowering time of male plant is earlier than female plant [54]. Hence, an interesting hypothesis is that up-regulation of a male factor may promote the male trait, such as early flowering. Hence, sol-miR2550n may be a male-promoting factor. The exact molecular mechanism needs to be explored in future work. Meanwhile, sol-miR172, a homolog of the well-known flower-related factor miR172 [53,55], was also studied in this work as a positive control. In spinach, sol-miR172 also showed higher expression level in M03 (male flower) than in FNS (female flower) but not resided in sex chromosome. Overexpression of sol-miR172 in Arabidopsis resulted in early flowering, in accordance with ath-miR172 (Figure 5). However, down-regulation of sol-miR172 in spinach resulted in anther abortion of male flower, which is not similar with ath-miR172 (influencing flower whorls pattern) (Figure 5). Such difference may be due to the architecture difference between unisexual flower and bisexual flower and also indicated a new potential regulation pathway of miR172 in unisexual flower development. The molecular mechanism will be studied in our future work.

4. Materials and Methods

4.1. Plant Materials and Growth Conditions

All plants used in this study were Spinacia oleracea L. cv DA JIAN YE BO CAI, a dioecious plant. Seeds were obtained from U.S. National Plant Germplasm System (https://npgsweb.ars-grin.gov/gringlobal/search.aspx?, accessed on 23 July 2021, accession number is PI 527332) and grown in an experimental field at Henan Normal University, Xinxiang, China (113.90° E, 35.32° N). Female and male flowers collected at three stages were separately used for Small RNA Sequencing and qRT-PCR. Flower samples were same with our previous published paper [52].
Spinach plants used for VIGS treatment were planted in climate chamber (Jiangnanyiqi, Ningbo, China) with 16 h/8 h day/night, 18 °C, 60% humidity. Arabidopsis plants were planted in climate chamber (Jiangnanyiqi, Ningbo, China) with 16 h/8 h day/night, 22 °C, and 60% humidity.

4.2. Library Construction and Sequencing

After total RNA was extracted by TRIzol (Thermo Fisher Scientific, Waltham, MA, USA), the RNA molecules in a size range of 18–30 nt were enriched by polyacrylamide gel electrophoresis (CWBIO, Taizhou, China). Then the 3′ adapters were added and the 36–44 nt RNAs were enriched. The 5′ adapters were then ligated to the RNAs as well. The ligation products were reverse transcribed by PCR amplification and the 140–160 bp size PCR products were enriched to generate a cDNA library and sequenced using Illumina HiSeq TM 2500 (Illumina, San Diego, CA, USA) by Gene Denovo Biotechnology Co. (Guangzhou, China).
Reads obtained from the sequencing machine included dirty reads containing adapters or low quality bases which would affect the following assembly and analysis. Thus, to acquire clean tags, raw reads were further filtered using fastp software (https://github.com/OpenGene/fastp, version 0.12.4, 1 March 2022) according to the following rules:
(1)
Removing low quality reads containing more than one low quality (Q-value ≤ 20) base or containing unknown nucleotides (N);
(2)
Removing reads without 3′adapters;
(3)
Removing reads containing 5′adapters;
(4)
Removing reads containing 3′ and 5′ adapters but no small RNA fragment between them;
(5)
Removing reads containing ployA (the content of A base in a reads is higher than 70%) in small RNA fragment;
(6)
Removing reads shorter than 18 nt (not include adapters).

4.3. Alignment and Identification of Small RNA

All of the clean tags were aligned with small RNAs in GeneBank database (https://www.ncbi.nlm.nih.gov/genbank/, Release 209.0, 1 March 2022) using blast (https://blast.ncbi.nlm.nih.gov/Blast.cgi, version 2.2.25, 1 March 2022) (blastn, identity > 97%) to identify and remove rRNA, scRNA, snoRNA, snRNA and tRNA. Meanwhile all of the clean tags were aligned with small RNAs in Rfam database (http://rfam.xfam.org/, version 11.0, 1 March 2022) using blast (https://blast.ncbi.nlm.nih.gov/Blast.cgi, version 2.2.25, 1 March 2022) (blastn, identity > 97%) to identify and remove rRNA, scRNA, snoRNA, snRNA and tRNA. All of the clean tags were also aligned with reference genome [44] using bowtie (https://www.nature.com/articles/nmeth.1923, version 1.1.2, 1 March 2022) (parameters: -v 0 --best --strata –a). Those mapped to exons or introns might be fragments from mRNA degradation, so these tags were removed. The tags mapped to repeat sequences using RepeatMasker (http://www.repeatmasker.org/, 1 March 2022, version open-4.0.6, RepeatMasker Database: RepeatMaskerLib.embl Update 20150807, parameters: RepeatMasker -engine wublast -s -no_is -cutoff 255 -frag 20000) were also removed.
After removing rRNA, scRNA, snoRNA, snRNA, tRNA, mRNA degradation fragments and repeat sequences, the clean tags were then searched against miRBase database (http://www.mirbase.org/, 1 March 2022, Release 21) using bowtie (http://bowtie-bio.sourceforge.net/index.shtml, 1 March 2022, version 1.1.2, parameters: -v 0 --best --strata -a) to identify exist and known miRNAs. Then the unidentified clean tags were mapped to reference genome [44]. According to their genome positions and hairpin structures predicted by software Mireap_v0.2 (https://sourceforge.net/projects/mireap/, 1 March 2022), the novel miRNA candidates were identified. The default parameters of software Mireap_v0.2 were as follows:
(1)
Minimal miRNA sequence length is 18 nt;
(2)
Maximal miRNA sequence length is 25 nt;
(3)
Minimal miRNA reference sequence length is 20 nt;
(4)
Maximal miRNA reference sequence length is 23 nt;
(5)
Maximal copy number of miRNAs on reference is 20;
(6)
Maximal free energy allowed for a miRNA precursor is 18 kcal/mol;
(7)
Maximal space between miRNA and miRNA* is 300 nt;
(8)
Minimal space between miRNA and miRNA* is 16 nt;
(9)
Maximal bulge between miRNA and miRNA* is 4 nt;
(10)
Maximal asymmetry of miRNA/miRNA* duplex is 4 nt;
(11)
Flank sequence length of miRNA precursor is 20 nt.

4.4. MiRNA Expression Profiles

Total miRNAs consist of existing miRNAs, known miRNAs and novel miRNAs, based on their expression in each sample, the miRNA expression level was calculated and normalized to transcripts per million (TPM). The formula is as follows: TPM = Actual miRNA counts/Total counts of clean tags*106.
To identify DE miRNAs across samples or groups the formula was shown as follows:
p ( x | y ) = ( N 2 N 1 ) y ( x + y ) ! x ! y ! ( 1 + N 2 N 1 ) ( x + y + 1 ) C ( y ≤ y m i n | x ) = ∑ v = 0 y ≤ y m i n p ( y | x ) D ( y ≥ y m a x | x ) = ∑ y ≥ m a x ∞ p ( y | x )
We identified miRNAs with a fold change ≥ 2 and p value < 0.05 in a comparison as significant DE miRNAs.

4.5. Target Gene Prediction

Based on the sequences of the exist miRNAs, known miRNAs and novel miRNAs, the candidate target genes were predicted. The software patmatch (ftp://ftp.arabidopsis.org/home/tair/Software/Patmatch/, 1 March 2022, version 1.2) was used to predict target genes. The default parameters were as follows:
(1)
No more than four mismatches between sRNA & target (G-U bases count as 0.5 mismatches)
(2)
No more than two adjacent mismatches in the miRNA/target duplex
(3)
No adjacent mismatches in in positions 2–12 of the miRNA/target duplex (5’ of miRNA)
(4)
No mismatches in positions 10–11 of miRNA/target duplex
(5)
No more than 2.5 mismatches in positions 1–12 of the of the miRNA/target duplex (5’ of miRNA)
(6)
Minimum free energy (MFE) of the miRNA/target duplex should be > = 74% of the MFE of the miRNA bound to it’s perfect complement.

4.6. qRT-PCR

E.Z.N.A.® Plant miRNA Kit (OMEGA, Norcross, GA, USA) and E.Z.N.A.® Plant RNA Kit (OMEGA, Norcross, GA, USA) was used to isolate miRNA and total RNA, separately, according to manufacturer’s introduction. Mir-X miRNA First-Strand Synthesis Kit (TaKaRa, Shiga, Japan) was used to perform reverse transcription of miRNA. PrimeScript™ RT reagent Kit with gDNA Eraser (TaKaRa, Shiga, Japan) was used to synthesis the cDNA of total RNA. Mir-X miRNA qRT-PCR TB Green Kit (TaKaRa, Shiga, Japan) and TB Green ® Premix Ex Taq™ II (Tli RNaseH Plus) (TaKaRa, Shiga, Japan) were used to perform the qRT-PCR of miRNA and mRNA, respectively. qRT-PCR was carried out on LightCycler® 480 System (Roche, ) according to the manufacturer’s instructions. SpoEF, SpoUBQ, At5sRNA, AtACTIN2 were used as an internal control for normalization. The relative expression level of transcript was calculated with 2−ΔΔCT method. All primers sequences are listed in Table S4.

4.7. Vector Construction and Plant Transformation

The pre-miRNA was cloned using Trans-Start® FastPfu Fly DNA Polymerase (TRANSGEN, Beijing, China) and inserted into the pTRV2 vector (VIGS) or pLP100-35s vector (overexpression) with ClonExpressⅡOne Step Cloning Kit (Vazyme, Nanjing, China). These vectors were transformed into Agrobacterium competent cells (GV3101, WEIDI, Shanghai, China). Spinach seedlings with four leaves were injected with Agrobacterium infection buffer including pTRV2-pre-miRNA vector; flower phenotype was observed after about 30 days. Arabidopsis flower buds were dipped in Agrobacterium infection buffer including pLP100-35s-pre-miRNA vector; positive T1 plants were used to observe phenotype.

5. Conclusions

It was the first time to sequencing small RNAs in spinach. 92 known and 3402 novel miRNAs were identified, which enriched plant miRNA resource. 74 DE miRNAs including male-biased/specific and female-biased/specific expression miRNAs were identified between male and female flowers. Moreover, target prediction identified 22 sex-biased microRNA-target pairs, which may be involved in spinach sex determination or differentiation. Genes residing in sex chromosome, especially in sex-determining region, are very important for sex determination or differentiation. Hence, miRNAs residing in sex chromosome were analyzed and 55 DE microRNAs between FNS and M03 were identified; and one of them, sol-miR2550n, was analyzed via genetic transformation. Silencing of sol-miR2550n resulted in abnormal anther while overexpression of sol-miR2550n induced early flowering, indicating sol-miR2550n may be a male-promoting factor. Conclusively, our work can supply valuable information for exploring spinach sex determination and differentiation and provide a new insight in studying unisexual flower development.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms23084090/s1,

Author Contributions

Conceptualization, S.L., W.G. and C.D.; data curation, N.L. and Y.W. (Yueyuan Wang); formal analysis, J.W. and Z.M.; funding acquisition, N.L., S.L., W.G. and C.D.; investigation, Y.W. (Yueyuan Wang), W.Z. and Z.M.; methodology, Y.W. (Yueyuan Wang), J.W., W.Z. and Y.Z.; project administration, N.L. and C.D.; resources, S.L., W.G. and C.D.; software, J.W., W.Z., Y.W. (Yuanshen Wang) and Y.Z.; supervision, C.D.; validation, Y.W. (Yueyuan Wang), W.Z. and Z.M.; visualization, J.W. and Y.W. (Yuanshen Wang); writing–original draft, N.L.; writing–review and editing, N.L., S.L., W.G. and C.D. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by grants from the National Natural Science Foundation of China (31770346, 31970240 and 32000243), the key scientific research project plan for Henan provincial colleges and universities (20A180016), and the Natural Science foundation of Henan province (202300410242).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Sequencing data that support the findings of this study has been deposited in the NCBI SRA database [https://www.ncbi.nlm.nih.gov/sra/PRJNA796451, accessed on 17 March 2022]. All data generated or analyzed during this study are included in this published article.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Renner, S.S. The relative and absolute frequencies of angiosperm sexual systems: Dioecy, monoecy, gynodioecy, and an updated online database. Am. J. Bot. 2014, 101, 1588–1596. [Google Scholar] [CrossRef] [Green Version]
  2. Barrett, S.C. The evolution of plant sexual diversity. Nat. Rev. Genet. 2002, 3, 274–284. [Google Scholar] [CrossRef]
  3. Charlesworth, D. Plant contributions to our understanding of sex chromosome evolution. New Phytol. 2015, 208, 52–65. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Henry, I.M.; Akagi, T.; Tao, R.; Comai, L. One Hundred Ways to Invent the Sexes: Theoretical and Observed Paths to Dioecy in Plants. Annu. Rev. Plant Biol. 2018, 69, 553–575. [Google Scholar] [CrossRef] [Green Version]
  5. Akagi, T.; Charlesworth, D. Pleiotropic effects of sex-determining genes in the evolution of dioecy in two plant species. Proc. Biol. Sci. 2019, 286, 20191805. [Google Scholar] [CrossRef]
  6. Feng, G.; Sanderson, B.J.; Keefoverring, K.; Liu, J.; Ma, T.; Yin, T.; Smart, L.B.; Difazio, S.P.; Olson, M.S. Pathways to sex determination in plants: How many roads lead to Rome? Curr. Opin. Plant Biol. 2020, 54, 61–68. [Google Scholar] [CrossRef]
  7. Fechter, I.; Hausmann, L.; Daum, M.; Sorensen, T.R.; Viehover, P.; Weisshaar, B.; Topfer, R. Candidate genes within a 143 kb region of the flower sex locus in Vitis. Mol. Genet. Genom. 2012, 287, 247–259. [Google Scholar] [CrossRef] [PubMed]
  8. Picq, S.; Santoni, S.; Lacombe, T.; Latreille, M.; Weber, A.; Ardisson, M.; Ivorra, S.; Maghradze, D.; Arroyo-Garcia, R.; Chatelet, P.; et al. A small XY chromosomal region explains sex determination in wild dioecious V. vinifera and the reversal to hermaphroditism in domesticated grapevines. BMC Plant Biol. 2014, 14, 229. [Google Scholar] [CrossRef] [Green Version]
  9. Akagi, T.; Henry, I.M.; Tao, R.; Comai, L. A Y-chromosome-encoded small RNA acts as a sex determinant in persimmons. Science 2014, 346, 646–650. [Google Scholar] [CrossRef] [PubMed]
  10. Geraldes, A.; Hefer, C.A.; Capron, A.; Kolosova, N.; Martinez-Nunez, F.; Soolanayakanahally, R.Y.; Stanton, B.; Guy, R.D.; Mansfield, S.D.; Douglas, C.J.; et al. Recent Y chromosome divergence despite ancient origin of dioecy in poplars (Populus). Mol. Ecol. 2015, 24, 3243–3256. [Google Scholar] [CrossRef] [PubMed]
  11. Muller, N.A.; Kersten, B.; Leite Montalvao, A.P.; Mahler, N.; Bernhardsson, C.; Brautigam, K.; Carracedo Lorenzo, Z.; Hoenicka, H.; Kumar, V.; Mader, M.; et al. A single gene underlies the dynamic evolution of poplar sex determination. Nat. Plants 2020, 6, 630–637. [Google Scholar] [CrossRef] [PubMed]
  12. Xue, L.; Wu, H.; Chen, Y.; Li, X.; Hou, J.; Lu, J.; Wei, S.; Dai, X.; Olson, M.S.; Liu, J.; et al. Evidences for a role of two Y-specific genes in sex determination in Populus deltoides. Nat. Commun. 2020, 11, 5893. [Google Scholar] [CrossRef]
  13. Mori, K.; Shirasawa, K.; Nogata, H.; Hirata, C.; Tashiro, K.; Habu, T.; Kim, S.; Himeno, S.; Kuhara, S.; Ikegami, H. Identification of RAN1 orthologue associated with sex determination through whole genome sequencing analysis in fig (Ficus carica L.). Sci. Rep. 2017, 7, 46784. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Torres, M.F.; Mathew, L.S.; Ahmed, I.; Al-Azwani, I.K.; Krueger, R.; Rivera-Nunez, D.; Mohamoud, Y.A.; Clark, A.G.; Suhre, K.; Malek, J.A. Genus-wide sequencing supports a two-locus model for sex-determination in Phoenix. Nat. Commun. 2018, 9, 3969. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Tennessen, J.A.; Wei, N.; Straub, S.; Govindarajulu, R.; Liston, A.; Ashman, T.L. Repeated translocation of a gene cassette drives sex-chromosome turnover in strawberries. PLoS Biol. 2018, 16, e2006062. [Google Scholar] [CrossRef] [Green Version]
  16. Akagi, T.; Pilkington, S.M.; Varkonyi-Gasic, E.; Henry, I.M.; Sugano, S.S.; Sonoda, M.; Firl, A.; McNeilage, M.A.; Douglas, M.J.; Wang, T.C.; et al. Two Y-chromosome-encoded genes determine sex in kiwifruit. Nat. Plants 2019, 5, 801–809. [Google Scholar] [CrossRef]
  17. Harkess, A.; Huang, K.; van der Hulst, R.; Tissen, B.; Caplan, J.L.; Koppula, A.; Batish, M.; Meyers, B.C.; Leebens-Mack, J. Sex Determination by Two Y-Linked Genes in Garden Asparagus. Plant Cell 2020, 32, 1790–1796. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Bergero, R.; Charlesworth, D. Preservation of the Y Transcriptome in a 10-Million-Year-Old Plant Sex Chromosome System. Curr. Biol. 2011, 21, 1470–1474. [Google Scholar] [CrossRef] [Green Version]
  19. Kazama, Y.; Ishii, K.; Aonuma, W.; Ikeda, T.; Kawamoto, H.; Koizumi, A.; Filatov, D.A.; Chibalina, M.; Bergero, R.; Charlesworth, D.; et al. A new physical mapping approach refines the sex-determining gene positions on the Silene latifolia Y-chromosome. Sci. Rep. 2016, 6, 18917. [Google Scholar] [CrossRef] [Green Version]
  20. Zerpa-Catanho, D.; Wai, J.; Wang, M.L.; Yu, L.; Nguyen, J.; Ming, R. Differential gene expression among three sex types reveals a MALE STERILITY 1 (CpMS1) for sex differentiation in papaya. BMC Plant Biol. 2019, 19, 545. [Google Scholar] [CrossRef]
  21. Voinnet, O. Origin, Biogenesis, and Activity of Plant MicroRNAs. Cell 2009, 136, 669–687. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Chellappan, P.; Xia, J.; Zhou, X.F.; Gao, S.; Zhang, X.M.; Coutino, G.; Vazquez, F.; Zhang, W.X.; Jin, H.L. siRNAs from miRNA sites mediate DNA methylation of target genes. Nucleic Acids Res. 2010, 38, 6883–6894. [Google Scholar] [CrossRef] [PubMed]
  23. Wu, L.; Zhou, H.Y.; Zhang, Q.Q.; Zhang, J.G.; Ni, F.R.; Liu, C.; Qi, Y.J. DNA Methylation Mediated by a MicroRNA Pathway. Mol. Cell 2010, 38, 465–475. [Google Scholar] [CrossRef]
  24. Feng, H.; Zhang, Q.; Wang, Q.; Wang, X.; Liu, J.; Li, M.; Huang, L.; Kang, Z. Target of tae-miR408, a chemocyanin-like protein gene (TaCLP1), plays positive roles in wheat response to high-salinity, heavy cupric stress and stripe rust. Plant Mol. Biol. 2013, 83, 433–443. [Google Scholar] [CrossRef]
  25. Hackenberg, M.; Shi, B.-J.; Gustafson, P.; Langridge, P. Characterization of phosphorus-regulated miR399 and miR827 and their isomirs in barley under phosphorus-sufficient and phosphorus-deficient conditions. BMC Plant Biol. 2013, 13, 214. [Google Scholar] [CrossRef] [Green Version]
  26. Han, R.; Jian, C.; Lv, J.Y.; Yan, Y.; Chi, Q.; Li, Z.J.; Wang, Q.; Zhang, J.; Liu, X.L.; Zhao, H.X. Identification and characterization of microRNAs in the flag leaf and developing seed of wheat (Triticum aestivum L.). BMC Genom. 2014, 15, 289. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Song, X.; Li, Y.; Cao, X.; Qi, Y. MicroRNAs and Their Regulatory Roles in Plant–Environment Interactions. Annu. Rev. Plant Biol. 2019, 70, 489–525. [Google Scholar] [CrossRef]
  28. Baker, C.C.; Sieber, P.; Wellmer, F.; Meyerowitz, E.M. The early extra petals1 mutant uncovers a role for MicroRNA miR164c in regulating petal number in Arabidopsis. Curr. Biol. 2005, 15, 303–315. [Google Scholar] [CrossRef] [Green Version]
  29. Chen, X. A MicroRNA as a Translational Repressor of APETALA2 in Arabidopsis Flower Development. Science 2004, 303, 2022. [Google Scholar] [CrossRef] [Green Version]
  30. Zhao, L.; Kim, Y.J.; Dinh, T.T.; Chen, X.M. miR172 regulates stem cell fate and defines the inner boundary of APETALA3 and PISTILLATA expression domain in Arabidopsis floral meristems. Plant J. 2007, 51, 840–849. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Rubio-Somoza, I.; Weigel, D. Coordination of Flower Maturation by a Regulatory Circuit of Three MicroRNAs. PLoS Genet. 2013, 9, e1003374. [Google Scholar] [CrossRef] [Green Version]
  32. Damodharan, S.; Corem, S.; Gupta, S.K.; Arazi, T. Tuning of SlARF10A dosage by sly-miR160a is critical for auxin-mediated compound leaf and flower development. Plant J. 2018, 96, 855–868. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Chuck, G.; Meeley, R.; Irish, E.; Sakai, H.; Hake, S. The maize tasselseed4 microRNA controls sex determination and meristem cell fate by targeting Tasselseed6/indeterminate spikelet1. Nat. Genet. 2007, 39, 1517–1521. [Google Scholar] [CrossRef] [PubMed]
  34. Akagi, T.; Shirasawa, K.; Nagasaki, H.; Hirakawa, H.; Tao, R.; Comai, L.; Henry, I.M. The persimmon genome reveals clues to the evolution of a lineage-specific sex determination system in plants. PLoS Genet. 2020, 16, e1008845. [Google Scholar] [CrossRef]
  35. Ming, R.; Bendahmane, A.; Renner, S.S. Sex Chromosomes in Land Plants. Annu. Rev. Plant Biol. 2011, 62, 485–514. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Onodera, Y.; Yonaha, I.; Masumo, H.; Tanaka, A.; Niikura, S.; Yamazaki, S.; Mikami, T. Mapping of the genes for dioecism and monoecism in Spinacia oleracea L.: Evidence that both genes are closely linked. Plant Cell Rep. 2011, 30, 965–971. [Google Scholar] [CrossRef] [PubMed]
  37. Deng, C.; Qin, R.; Gao, J.; Cao, Y.; Li, S.; Gao, W.; Lu, L. Identification of sex chromosome of spinach by physical mapping of 45s rDNAs by FISH. Caryologia 2012, 65, 322–327. [Google Scholar] [CrossRef] [Green Version]
  38. Kudoh, T.; Takahashi, M.; Osabe, T.; Toyoda, A.; Hirakawa, H.; Suzuki, Y.; Ohmido, N.; Onodera, Y. Molecular insights into the non-recombining nature of the spinach male-determining region. Mol. Genet. Genom. 2018, 293, 557–568. [Google Scholar] [CrossRef]
  39. Wadlington, W.H.; Ming, R. Development of an X-specific marker and identification of YY individuals in spinach. Theor. Appl. Genet. 2018, 131, 1987–1994. [Google Scholar] [CrossRef]
  40. Qian, W.; Fan, G.Y.; Liu, D.D.; Zhang, H.L.; Wang, X.W.; Wu, J.; Xu, Z.S. Construction of a high-density genetic map and the X/Y sex-determining gene mapping in spinach based on large-scale markers developed by specific-locus amplified fragment sequencing (SLAF-seq). BMC Genom. 2017, 18, 276. [Google Scholar] [CrossRef] [Green Version]
  41. Yu, L.; Ma, X.; Deng, B.; Yue, J.; Ming, R. Construction of high-density genetic maps defined sex determination region of the Y chromosome in spinach. Mol. Genet. Genom. 2021, 296, 41–53. [Google Scholar] [CrossRef] [PubMed]
  42. Ma, X.; Yu, L.; Fatima, M.; Wadlington, W.H.; Hulse-Kemp, A.M.; Zhang, X.; Zhang, S.; Xu, X.; Wang, J.; Huang, H.; et al. The spinach YY genome reveals sex chromosome evolution, domestication, and introgression history of the species. Genome Biol. 2022, 23, 75. [Google Scholar] [CrossRef] [PubMed]
  43. Okazaki, Y.; Takahata, S.; Hirakawa, H.; Suzuki, Y.; Onodera, Y. Molecular evidence for recent divergence of X- and Y-linked gene pairs in Spinacia oleracea L. PLoS ONE 2019, 14, e0214949. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Xu, C.; Jiao, C.; Sun, H.; Cai, X.; Wang, X.; Ge, C.; Zheng, Y.; Liu, W.; Sun, X.; Xu, Y.; et al. Draft genome of spinach and transcriptome diversity of 120 Spinacia accessions. Nat. Commun. 2017, 8, 15275. [Google Scholar] [CrossRef] [PubMed]
  45. Hirakawa, H.; Toyoda, A.; Itoh, T.; Suzuki, Y.; Nagano, A.J.; Sugiyama, S.; Onodera, Y. A spinach genome assembly with remarkable completeness, and its use for rapid identification of candidate genes for agronomic traits. DNA Res. 2021, 28, dsab004. [Google Scholar] [CrossRef] [PubMed]
  46. Cai, X.; Sun, X.; Xu, C.; Sun, H.; Wang, X.; Ge, C.; Zhang, Z.; Wang, Q.; Fei, Z.; Jiao, C.; et al. Genomic analyses provide insights into spinach domestication and the genetic basis of agronomic traits. Nat. Commun. 2021, 12, 7246. [Google Scholar] [CrossRef]
  47. Avsar, B.; Aliabadi, D.E. In silico Analysis of MicroRNAs in Spinacia oleracea Genome and Transcriptome. Int. J. Biosci. Biochem. Bioinform. 2017, 7, 84–92. [Google Scholar] [CrossRef] [Green Version]
  48. Wang, P.; Liao, H.; Zhang, W.; Yu, X.; Zhang, R.; Shan, H.; Duan, X.; Yao, X.; Kong, H. Flexibility in the structure of spiral flowers and its underlying mechanisms. Nat. Plants 2015, 2, 15188. [Google Scholar] [CrossRef]
  49. West, N.W.; Golenberg, E.M. Gender-specific expression of GIBBERELLIC ACID INSENSITIVE is critical for unisexual organ initiation in dioecious Spinacia oleracea. New Phytol. 2018, 217, 1322–1334. [Google Scholar] [CrossRef] [Green Version]
  50. Kim, H.; Onodera, Y.; Masuta, C. Application of cucumber mosaic virus to efficient induction and long-term maintenance of virus-induced gene silencing in spinach. Plant Biotechnol. 2020, 37, 83–88. [Google Scholar] [CrossRef] [Green Version]
  51. Hu, Y.; Cheng, H.; Zhang, Y.; Zhang, J.; Niu, S.; Wang, X.; Li, W.; Zhang, J.; Yao, Y. The MdMYB16/MdMYB1-miR7125-MdCCR module regulates the homeostasis between anthocyanin and lignin biosynthesis during light induction in apple. New Phytol. 2021, 231, 1105–1122. [Google Scholar] [CrossRef] [PubMed]
  52. Li, N.; Meng, Z.; Tao, M.; Wang, Y.; Zhang, Y.; Li, S.; Gao, W.; Deng, C. Comparative transcriptome analysis of male and female flowers in Spinacia oleracea L. BMC Genom. 2020, 21, 850. [Google Scholar] [CrossRef]
  53. Aukerman, M.J.; Sakai, H. Regulation of flowering time and floral organ identity by a MicroRNA and its APETALA2-like target genes. Plant Cell 2003, 15, 2730–2741. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Rosa, J.T. Sex expression in spinach. Hilgardia 1925, 1, 259–274. [Google Scholar] [CrossRef] [Green Version]
  55. Zhang, B.; Chen, X. Secrets of the MIR172 family in plant development and flowering unveiled. PLoS Biol. 2021, 19, e3001099. [Google Scholar] [CrossRef] [PubMed]
Ijms 23 04090 g001 550
Figure 1. Frequency distribution of small RNA length. FNS, female flower without stigma (0.3 mm in diameter); FNB, female flower without stigma (0.5 mm in diameter); FYS, mature female flower with stigma; M03, male flower (0.3 mm in diameter); M05, male flower (0.5 mm in diameter); M10, male flower (1.0 mm in diameter).
Figure 1. Frequency distribution of small RNA length. FNS, female flower without stigma (0.3 mm in diameter); FNB, female flower without stigma (0.5 mm in diameter); FYS, mature female flower with stigma; M03, male flower (0.3 mm in diameter); M05, male flower (0.5 mm in diameter); M10, male flower (1.0 mm in diameter).
Ijms 23 04090 g001
Ijms 23 04090 g002 550
Figure 2. Differentially expressed miRNAs. (a) the Venn diagram of DE miRNAs among three flower stages; (b) the number of sex-biased expression miRNAs.
Figure 2. Differentially expressed miRNAs. (a) the Venn diagram of DE miRNAs among three flower stages; (b) the number of sex-biased expression miRNAs.
Ijms 23 04090 g002
Ijms 23 04090 g003 550
Figure 3. 22 pairs of sex-biased miRNA-target. The heatmaps represent the expression profile of 15 sex-biased miRNAs (left) and 16 sex-biased target genes (right); the dotted line links miRNA and its target.
Figure 3. 22 pairs of sex-biased miRNA-target. The heatmaps represent the expression profile of 15 sex-biased miRNAs (left) and 16 sex-biased target genes (right); the dotted line links miRNA and its target.
Ijms 23 04090 g003
Ijms 23 04090 g004 550
Figure 4. Distribution of known miRNAs and novel miRNAs in spinach chromosomes.
Figure 4. Distribution of known miRNAs and novel miRNAs in spinach chromosomes.
Ijms 23 04090 g004
Ijms 23 04090 g005 550
Figure 5. Phenotypes of silencing and overexpression of sol-miR2550n and sol-miR172. (a) anther phenotype of wildtype (WT), sol-miR2550n- or sol-miR172-silencing spinach male plant, scale bar = 0.5 mm; (b) phenotype of wildtype, sol-miR2550n- or sol-miR172-overexpressed Arabidopsis plant; (c) relative expression level of sol-miR2550n in wildtype and sol-miR2550n-silencing spinach plant; (d) relative expression level of sol-miR2550n in wildtype and sol-miR2550n-overexpressed Arabidopsis plant; (e) statistics of flowering time of sol-miR2550n-overexpressed Arabidopsis plants, n ≥ 20; (f) relative expression level of sol-miR172 in wildtype and sol-miR172-silencingspianch plant; (g) relative expression level of sol-miR172 in wildtype and sol-miR172-overexpressed Arabidopsis plant; (h) statistics of flowering time of sol-miR172-overexpressed Arabidopsis plants, n ≥ 20; “*” represents p-value < 0.05; “**” represents p-value < 0.01; “***” represents p-value < 0.001; error bar represents standard error.
Figure 5. Phenotypes of silencing and overexpression of sol-miR2550n and sol-miR172. (a) anther phenotype of wildtype (WT), sol-miR2550n- or sol-miR172-silencing spinach male plant, scale bar = 0.5 mm; (b) phenotype of wildtype, sol-miR2550n- or sol-miR172-overexpressed Arabidopsis plant; (c) relative expression level of sol-miR2550n in wildtype and sol-miR2550n-silencing spinach plant; (d) relative expression level of sol-miR2550n in wildtype and sol-miR2550n-overexpressed Arabidopsis plant; (e) statistics of flowering time of sol-miR2550n-overexpressed Arabidopsis plants, n ≥ 20; (f) relative expression level of sol-miR172 in wildtype and sol-miR172-silencingspianch plant; (g) relative expression level of sol-miR172 in wildtype and sol-miR172-overexpressed Arabidopsis plant; (h) statistics of flowering time of sol-miR172-overexpressed Arabidopsis plants, n ≥ 20; “*” represents p-value < 0.05; “**” represents p-value < 0.01; “***” represents p-value < 0.001; error bar represents standard error.
Ijms 23 04090 g005
Table
Table 1. Summary of small RNA sequencing.
Table 1. Summary of small RNA sequencing.
Sample IDClean ReadsHigh Quality3’Adapter NullInsert Null5’Adapter ContaminantsSmaller than 18 ntPolyAClean Tags
SpM10-114,183,245 (100%)13,907,336 (98.0547%)66,715 (0.4797%)39,362 (0.2830%)68,350 (0.4915%)193,229 (1.3894%)2113 (0.0152%)13,537,567 (97.3412%)
SpM10-214,464,387 (100%)14,164,634 (97.9276%)80,069 (0.5653%)67,309 (0.4752%)19,248 (0.1359%)553,683 (3.9089%)1861 (0.0131%)13,442,464 (94.9016%)
SpM10-314,506,554 (100%)14,227,721 (98.0779%)77,086 (0.5418%)117,860 (0.8284%)37,183 (0.2613%)1,254,040 (8.8141%)1497 (0.0105%)12,740,055 (89.5439%)
SpM05-115,642,311 (100%)15,300,380 (97.8141%)34,255 (0.2239%)129,458 (0.8461%)29,207 (0.1909%)778,380 (5.0873%)2248 (0.0147%)14,326,832 (93.6371%)
SpM05-215,687,586 (100%)15,366,874 (97.9556%)39,542 (0.2573%)104,255 (0.6784%)42,694 (0.2778%)514,753 (3.3498%)1978 (0.0129%)14,663,652 (95.4238%)
SpM05-313,426,285 (100%)13,140,154 (97.8689%)60,121 (0.4575%)96,948 (0.7378%)23,471 (0.1786%)566,578 (4.3118%)1664 (0.0127%)12,391,372 (94.3016%)
SpM03-114,134,390 (100%)13,839,374 (97.9128%)65,717 (0.4749%)87,932 (0.6354%)25,498 (0.1842%)619,950 (4.4796%)2138 (0.0154%)13,038,139 (94.2105%)
SpM03-214,974,564 (100%)14,661,075 (97.9065%)51,713 (0.3527%)49,690 (0.3389%)27,909 (0.1904%)716,623 (4.8879%)2337 (0.0159%)13,812,803 (94.2141%)
SpM03-314,763,485 (100%)14,459,587 (97.9416%)59,870 (0.4141%)102,487 (0.7088%)25,196 (0.1743%)1,058,011 (7.3170%)1825 (0.0126%)13,212,198 (91.3733%)
SpFYS-114,686,982 (100%)14,367,641 (97.8257%)42,045 (0.2926%)57,856 (0.4027%)18,191 (0.1266%)388,472 (2.7038%)1739 (0.0121%)13,859,338 (96.4622%)
SpFYS-216,785,272 (100%)16,407,458 (97.7491%)84,234 (0.5134%)85,457 (0.5208%)28,469 (0.1735%)790,536 (4.8182%)1885 (0.0115%)15,416,877 (93.9626%)
SpFYS-313,411,587 (100%)13,139,123 (97.9684%)59,803 (0.4552%)75,025 (0.5710%)19,904 (0.1515%)554,722 (4.2219%)1697 (0.0129%)12,427,972 (94.5875%)
SpFNB-112,270,405 (100%)12,010,945 (97.8855%)38,498 (0.3205%)53,131 (0.4424%)14,985 (0.1248%)314,233 (2.6162%)1141 (0.0095%)11,588,957 (96.4866%)
SpFNB-213,192,504 (100%)12,938,214 (98.0725%)26,417 (0.2042%)54,149 (0.4185%)21,619 (0.1671%)550,952 (4.2583%)1196 (0.0092%)12,283,881 (94.9426%)
SpFNB-312,539,936 (100%)12,279,735 (97.9250%)53,903 (0.4390%)38,428 (0.3129%)17,344 (0.1412%)364,801 (2.9708%)1132 (0.0092%)11,804,127 (96.1269%)
SpFNS-114,777,105 (100%)14,461,450 (97.8639%)44,934 (0.3107%)58,975 (0.4078%)19,961 (0.1380%)822,130 (5.6850%)1283 (0.0089%)13,514,167 (93.4496%)
SpFNS-213,104,189 (100%)12,831,211 (97.9169%)26,074 (0.2032%)66,202 (0.5159%)19,029 (0.1483%)426,303 (3.3224%)1346 (0.0105%)12,292,257 (95.7997%)
SpFNS-313,170,448 (100%)12,905,639 (97.9894%)97,757 (0.7575%)67,113 (0.5200%)18,541 (0.1437%)619,765 (4.8023%)1122 (0.0087%)12,101,341 (93.7679%)
Table
Table 2. Summary of miRNAs identified in each sample.
Table 2. Summary of miRNAs identified in each sample.
SampleTags
Total		Known miRNANovel miRNA
miRNA
Number		Tags
Number		Ratio amiRNA
Number		Tags
Number		Ratio b
SpM10-1135375671646552854.84%2434661030.49%
SpM10-2134424641797049225.24%2503814060.61%
SpM10-3127400551789132517.17%2216789220.62%
SpM05-1143268321836948894.85%2416614410.43%
SpM05-21466365216610056956.86%2338619390.42%
SpM05-3123913721736774575.47%2384584270.47%
SpM03-1130381391574498643.45%2581958570.74%
SpM03-2138128032196626344.80%25631050930.76%
SpM03-3132121981947346635.56%23861215330.92%
SpFYS-1138593381733475532.51%2772743580.54%
SpFYS-2154168772003588722.33%27901044830.68%
SpFYS-3124279721633604862.90%2499697760.56%
SpFNB-1115889571407339486.33%1980482370.42%
SpFNB-2122838811556569345.35%1768427690.35%
SpFNB-3118041271317680996.51%1796458680.39%
SpFNS-11351416716310107727.48%1857743670.55%
SpFNS-2122922571449777607.95%2107591930.48%
SpFNS-3121013411566876845.68%2032408740.34%
Total c 92 3402
Note:a, the ratio of known miRNA Tags number to Tags total; b, the ratio of novel miRNA Tags number to Tags total; c, the total number of known or novel miRNAs identified in all samples.
Table
Table 3. Statistics of target gene prediction of all miRNAs.
Table 3. Statistics of target gene prediction of all miRNAs.
Sample NamemiRNA NumberTarget Gene NumberTarget Site Number
SpM10-11726716114391
SpM10-21796761915036
SpM10-31598783215006
SpM05-11746835915876
SpM05-21679651813201
SpM05-31720774415111
SpM03-11822687413919
SpM03-21893945118542
SpM03-31712904517370
SpFYS-11952664614188
SpFYS-21976781716164
SpFYS-31744651813423
SpFNB-11430540110784
SpFNB-21323535510423
SpFNB-3129450919762
SpFNS-11361657911874
SpFNS-21523530010930
SpFNS-31481648112669
Total28552046049591
Table
Table 4. DE miRNAs with targets between FNS and M03.
Table 4. DE miRNAs with targets between FNS and M03.
miRNA IDTargetAnnotation
novel-m2307-5pIsoform014322ENTH/ANTH/VHS superfamily protein isoform 1
novel-m2430-5pIsoform009669protein IQ-DOMAIN 14
novel-m2524-5pIsoform010547F-box and leucine-rich repeat protein 2/20
novel-m2550-5pIsoform007857defensin Ec-AMP-D2-like
novel-m2554-3pIsoform007343probable prefoldin subunit 2
novel-m2572-3pSpo15344plant mobile domain family protein
novel-m2617-3pIsoform008852uncharacterized LOC104898421
novel-m2619-3pSpo05571Polygalacturonase (PG) (3.2.1.15) (Pectinase) (Precursor)
novel-m2641-5pIsoform0020287-hydroxymethyl chlorophyll a reductase, chloroplastic
Isoform002329ABC transporter A family member 7-like
Isoform007233cysteine-rich and transmembrane domain-containing protein A
Isoform014680transcription factor bHLH90
novel-m2763-3pIsoform000331probable xyloglucan endotransglucosylase/hydrolase protein 5
Isoform011756uncharacterized LOC104897309
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, N.; Wang, Y.; Wang, J.; Zhang, W.; Meng, Z.; Wang, Y.; Zhang, Y.; Li, S.; Gao, W.; Deng, C. Identification of Sex Differentiation-Related microRNAs in Spinach Female and Male Flower. Int. J. Mol. Sci. 2022, 23, 4090. https://doi.org/10.3390/ijms23084090

AMA Style

Li N, Wang Y, Wang J, Zhang W, Meng Z, Wang Y, Zhang Y, Li S, Gao W, Deng C. Identification of Sex Differentiation-Related microRNAs in Spinach Female and Male Flower. International Journal of Molecular Sciences. 2022; 23(8):4090. https://doi.org/10.3390/ijms23084090

Chicago/Turabian Style

Li, Ning, Yueyuan Wang, Jiwen Wang, Wanqing Zhang, Ziwei Meng, Yuanshen Wang, Yulan Zhang, Shufen Li, Wujun Gao, and Chuanliang Deng. 2022. "Identification of Sex Differentiation-Related microRNAs in Spinach Female and Male Flower" International Journal of Molecular Sciences 23, no. 8: 4090. https://doi.org/10.3390/ijms23084090

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop