Sex specific expression and distribution of small RNAs in papaya
© Aryal et al.; licensee BioMed Central Ltd. 2014
Received: 3 October 2013
Accepted: 30 December 2013
Published: 13 January 2014
Regulatory function of small non-coding RNAs (sRNA) in response to environmental and developmental cues has been established. Additionally, sRNA, also plays an important role in maintaining the heterochromatin and centromere structures of the chromosome. Papaya, a trioecious species with recently evolved sex chromosomes, has emerged as an excellent model system to study sex determination and sex chromosome evolution in plants. However, role of small RNA in papaya sex determination is yet to be explored.
We analyzed the high throughput sRNAs reads in the Illumina libraries prepared from male, female, and hermaphrodite flowers of papaya. Using the sRNA reads, we identified 29 miRNAs that were not previously reported from papaya. Including this and two previous studies, a total of 90 miRNAs has been identified in papaya. We analyzed the expression of these miRNAs in each sex types. A total of 65 miRNAs, including 31 conserved and 34 novel mirNA, were detected in at least one library. Fourteen of the 65 miRNAs were differentially expressed among different sex types. Most of the miRNA expressed higher in male flowers were related to the auxin signaling pathways, whereas the miRNAs expressed higher in female flowers were the potential regulators of the apical meristem identity genes. Aligning the sRNA reads identified the sRNA hotspots adjacent to the gaps of the X and Y chromosomes. The X and Y chromosomes sRNA hotspots has a 7.8 and 4.4 folds higher expression of sRNA, respectively, relative to the chromosome wide average. Approximately 75% of the reads aligned to the X chromosome hotspot was identical to that of the Y chromosome hotspot.
By analyzing the large-scale sRNA sequences from three sex types, we identified the sRNA hotspots flanking the gaps of papaya X, Y, and Yh chromosome. The sRNAs expression patterns in these regions were reminiscent of the pericentromeric region indicating that the only remaining gap in each of these chromosomes is likely the centromere. We also identified 14 differentially expressed miRNAs in male, female and hermaphrodite flowers of papaya. Our results provide valuable information toward understanding the papaya sex determination.
KeywordsCarica papaya Centromere miRNA siRNA Sex chromosome Sex determination
Micro RNA (miRNA) and small interfering RNA (siRNA) are two major classes of endogenous regulatory RNAs (sRNA) found in higher plants. These sRNAs are processed from RNA duplexes by a dicer family protein, which produces approximately 21-24nt final products. Many miRNAs regulate various developmental processes by sequence directed silencing of the mRNA at a posttranscriptional level [1–5]. The siRNAs regulate the genome function both at transcriptional and posttranscriptional levels through RNA guided DNA methylation and RNA guided RNA silencing, respectively [6–9].
Spatial and temporal pattern of organ formation in plants are regulated at various levels by sRNA mediated gene silencing. The development of the male and female gametophyte is regulated by the combinatorial action of various sRNAs [10, 11]. In Arabidopsis, early embryo patterning and transition from juvenile to adult plant is regulated by miR156 [12–14]. In rice, Arabidopsis, and maize, mir166 regulates the adaxial/abaxial patterning of the leaves . Micro-RNAs and other sRNAs are important regulators of flower development and floral organ identity in many plant species [13, 16–21]. Further analysis of sRNA transcriptome in various organs, tissues, and developmental phases will provide a better understanding on their function in plant development and organogenesis.
Previous reports suggest that miRNAs play a role in sex differentiation in plants [21–23]. Differential expression of miRNAs that regulates the development of carpels and stamens can ultimately lead to the development of unisexual flowers and to the sexual differentiation in plants. In maize, miR172 maintains the unisexual nature of the tassel by suppressing carpel development genes Ts6. In Petunia hybrida and Antirrhinum majus, the miR169 family genes, miRBL and miRFIS, restrict the expression of C class homeotic genes at the center of the flower. Plants defective in these miRNA genes produce the female flowers on genetically hermaphroditic individuals . Small RNAs are also important in maintaining the integrity of sex chromosomes by methylating the non-recombining regions during meiotic prophase .
Papaya (Carica papaya) is a model species to study the plant sex determination. It is a trioecious species with three sex types, hermaphrodite, male, and female. It belongs to the family Caricaceae, which comprises 35 species in six genera including one monoecious, 32 dioecious, and two trioecious species. Papaya diverged from its closest monoecious relative (Vasconcellea monoica) about 27.5 million years ago . Sex expression in papaya is controlled by a pair of recently evolved sex chromosomes. The male and hermaphrodite characteristics are determined by two slightly different Y-chromosomes, Y and Yh, respectively [26, 27]. The genotype XX determines female, XY determines male, and XYh determines hermaphrodite. All combinations of Y and Yh are embryonic lethal, indicating that the Y and Yh chromosomes have lost some genes necessary for embryo development . The Y and Yh chromosomes show 98.9% sequence similarity on average, differing mainly in intergenic and repetitive regions . The male specific region of Y chromosome (MSY) and hermaphrodite specific region of Yh (HSY) are highly methylated and heterochromatized compared to the corresponding region of X chromosome . Because sRNA is involved in establishing heterochromatic structures, analyzing the sRNA transcripts from the papaya sex chromosome will help to understand the epigenetic aspects of the recently evolved sex chromosomes.
Sequencing the HSY and its corresponding region in the X chromosome have produced 8.1 and 3.4 Mb pseudomolecules for the sex determining regions of the respective chromosomes [30, 31]. Sequence alignment between HSY and corresponding region of the X chromosome revealed two large-scale inversions and several intrachromosomal rearrangements were identified on the HSY, however precise location of these inversions are still unknown due to an unfilled gap on the physical map. In the X chromosome, the repeat density gradually increases towards the gap, suggesting that the gap is likely the centromere . Using the fluorescent in situ hybridization technique, the centromere of the Y chromosomes was identified near the knob 4 . The knob 4 region is now completely sequenced, however, analysis of repeat content in the region did not identify centromeric satellite elements in this region . These conflicting observations call for the search of centromere on the Y chromosomes. Identification of centromere will provide valuable information regarding the location of inversion and direction of sequence expansion at the early stage of papaya sex chromosome evolution.
We previously analyzed the genome-wide sRNA profile in papaya using the high throughput Illumina sRNA libraries prepared from flowers and leaves including Papaya Ringspot Virus (PRSV) infected leaves . Evidences from many plant species show that sex determination in plants are also regulated by epigenetic mechanisms including that of small noncoding RNAs . In this study, we analyzed the sRNAs libraries prepared from male, female, and hermaphrodite flowers of papaya. We used the high throughput sRNA sequences to identify the miRNAs in different sex types. Expression of miRNAs in the flowers of three papaya sex types is analyzed. Aligning the sRNA reads identified the sRNA hotspots on papaya sex chromosomes similar to that of pericentromeric region. We discuss the early events in papaya sex chromosome evolution in light of the identified centromeric location.
Sequencing sRNAs from papaya flowers
We previously reported a higher accumulation of purine-rich strands (sequences with more purine residue than pyrimidine residue) than the pyrimidine-rich strands (sequences with more pyrimidine residue than purine residue) in sRNA libraries from various plant species . To further confirm this phenomenon, we analyzed the purine-rich and pyrimidine-rich sequences in the three libraries. In congruence with the previous observation, all three libraries were overrepresented by purine-rich strands (Additional file 1: Figure S1). Purine-rich sequences constituted 62%, 65%, and 65% of the total reads in male, female and hermaphrodite libraries, respectively, whereas pyrimidine-rich sequences constituted 30%, 28% and 25% in the respective libraries. Sequences with equal number of purine and pyrimidine residues were found 8% in the male library, 7% in female the library, and 10% in the hermaphrodite library. To see whether the purine-pyrimidine bias is localized on a particular region or is spread across the sRNA sequences, we analyzed the frequency of nucleotides in each position of 21nt and 24nt sequences (Additional file 1: Figure S2). Adenine was the most frequent nucleotide throughout the 24nt sequence, whereas 21nt sequences were overrepresented by guanine nucleotides. Uracil and adenine was highly conserved at the 5’ end of 21nt and 24nt sequences, respectively.
Identification of miRNAs from the sRNA libraries
List of papaya miRNAs identified in this study
MiRNA abundance in papaya flowers
Normalized expression of miRNA in flowers of different sex types
Expression (counts per million reads)
In general, conserved miRNAs were expressed higher in all sex types than the novel miRNAs. Of the 34 novel miRNAs, only 7 were expressed more than 100 copies, whereas 24 of the 31 conserved miRNAs were expressed more than 100 copies. MiR166 showed highest overall expression followed by miR156 and miR159.
Small RNA landscape on papaya sex chromosomes
To analyze the sRNA transcriptome of the sex chromosomes, we mapped the sRNA reads from each library to the respective pseudomolecules . The adapter trimmed non-redundant (unique) reads were aligned to the male specific region of Y chromosome (MSY), HSY, and the corresponding region on X chromosome. Only the reads with 100% sequence alignment were taken for the downstream analysis. The number of reads mapped to the MSY was higher than that of HSY. Approximately 3.5% of the unique reads from each library were aligned to the HSY, 5.5% to the MSY, and 4% to the X chromosome (Additional file 1: Table S1).
The papaya MSY and HSY have been mapped near the centromere of both Y and Yh chromosomes [41–43]. A physical map of the MSY and HSY has been generated, but a large unfilled gap remains at the ‘border A’ region. Similarly, a gap remains in the middle of the X chromosome physical map [30, 39]. Analysis of repetitive sequences provided strong evidence that the centromere is located in the middle of the X chromosome where the gap remains on the physical map . We observed a sharp increase in sRNA expression at both edges of the gap in the X chromosome, which further strengthened the notion that the gap is the centromere. For the Y chromosomes, it was suggested that the centromere might be at either side of Knob 4 in the HSY based on fluorescent in situ hybridization of Knob-specific BACs on anaphase chromosomes . If this were the case, we would have mapped and sequenced the centromere as there is no gap in that region [30, 39], but there is no long tandem repeats nor the rise and fall of repetitive sequences, typical features of centromeres. Our evidence of sharp increase of sRNA at the edge of gap in the HSY and MSY strongly supports that the only remaining gap in the physical map is likely the centromere of the Y chromosomes (Figure 3). Identification of the Y chromosome centromere is a major advancement on papaya sex chromosome research, and set the stage for further characterization and eventual sequencing of the Y centromere. This finding also indicated that the first inversion that triggered sex chromosome evolution is pericentric.
The gap in HSY and MSY include a large genomic region corresponding to the Knob 1 that shared between X and Y chromosomes, and chromosome walking on the HSY kept landing on the X chromosome counterpart, which was mapped and sequenced. The next question would be how extensive of the pericentric region on the unmapped side of the gap in HSY and MSY. The question can be addressed by examine the X specific genes between the two evolutionary strata in the X chromosome. Analysis of X specific genes between genes corresponding to the edges of inversions 1 and 2 revealed only one X specific gene, suggesting that the unmapped side of gap in the HSY and MSY involved few genes in the Inversion 1, and consider faster rate of gene loss in the HSY , it could be only one (paired with the one X specific gene between the two strata) or none gene on the unmapped side of the centromere in inversion 1.
Higher expression sRNAs at the pericentromeric region of the chromosome has been observed in other plant species [44–46]. The sRNA sequences aligned to the pericentromeric regions were highly conserved between X and Y chromosomes with up to 76% of sequence mapped to these regions being common in both, providing further evidence for the location of the centromere in the respective gaps. High conservation among the sRNA sequences expressed from the pericentromeric regions of different chromosomes was also observed in the yeast centromeres . The centromeric sequences of the chromosomes remain elusive in the assemblies of many genomes due to their highly repetitive nature . To date, complete sequence for the centromeric region is available only for Chromosome 8 of rice that was recently evolved . The fact that only one gap each remaining in the papaya X, HSY and MSY physical maps validates the high quality of these physical maps and also indicated that these centromeres are likely the original centromere of the autosomes that these sex chromosomes evolved from. The alignment of candidate centromere position between papaya X chromosome and the orthologous autosome of V. monoica is in line with this conclusion. The inversion and degeneration of the Y chromosome did not destroy or caused turnover of the Y chromosome centromere.
Premature separation of the papaya sex chromosomes compared to the autosome has been detected in meiotic anaphase . This premature separation is likely the consequence of suppressed recombination in the sex specific region and improper alignment of the centromere due to their shifted position relative to each other caused by the pericentric inversion.
A total of 14 miRNAs were differentially expressed among male, female, and hermaphrodite flowers, indicating their potential function in papaya sex determination (Table 2, Figure 2). All conserved miRNAs that are expressed highest in male flowers (miR160, miR167a, miR167b, miR169, and miR393) regulates the genes in auxin signaling pathway [36–38, 50]. Experimental evidences indicate that auxin plays a central role in carpel development (see  for more review). Additionally, miR169 regulates floral development in Nocotiana benthamiana, Petunia hybrid and Antirrhinum majus by spatial restriction of ‘C’ class floral homeotic genes [18, 52]. In congruent with these studies, our result shows a gradual difference in miR169 expression pattern among different sexes – highest in male (>5 fold compared to female), intermediate in hermaphrodite (~2 fold compared to female) and lowest in female. The miRNAs expressed highest in the papaya female flowers (miR164, miR166, and miR394) are mainly involved in regulating the embryo patterning and floral meristem identity genes [53–55]. Although whether these differences are the cause or consequences of sexual dimorphism can not be concluded now, it is worth investigating these miRNAs for their roles in sex expression in papaya.
Using the large scale sRNA sequences from the male, female and hermaphrodite flowers of papaya, we have identified, 1) the sRNA hotspot on the papaya sex chromosome reminiscent of pericentromeric region, and 2) differentially expressed miRNAs in the flowers of different sex types. The identified pericentromeric regions of the sex chromosomes are located adjacent to the unfilled gap on the physical map, indicating that the centromere of these chromosomes lies in the gap. Relative position of the pericentromeric region on X and Y chromosomes revealed that the centromere of X and Y chromosomes are located 1.6 Mb apart from each other, indicating that the inversion on Y chromosome occurred at the pericentromeric region spanning the centromere. Our results provide valuable information for further characterization of papaya sex chromosomes.
By analyzing the miRNA expression in papaya flowers, we have identified 12 miRNAs differentially expressed among the three sex types. Majority of the miRNAs expressed higher in male flower targets the genes involved in auxin signaling pathway. We observed a higher expression of miR169 in male flowers, which has been implicated to repress ‘C’ class floral homeotic genes in other plant species. Further functional analysis of these miRNAs may reveal their role in papaya sex determination.
Small RNA library construction and sequencing
The small RNA libraries were constructed from male, female and hermaphrodite flowers of papaya. Mixture of flowers at different developmental stage, including young bud to fully mature flowers, was used for RNA extraction. Illumina sRNA library was prepared for each sex type as described previously . In brief, total RNA was isolated using the Trizol reagent (Invitrogen USA). Low molecular weight RNA was enriched from the total RNA by precipitating with 0.4 M NaCl and polyethylene glycol (PEG). The enriched low molecular weight RNA was then separated on a denaturing 15% polyacrylamide gel. The band corresponding to the RNA standard of 18–30 nucleotides was then excised from the gel and eluted overnight in 0.4 M NaCl at 4°C, and subsequently ligated with 5’ and 3’ Illumina small RNA adapters. The adapter-ligated library was converted into cDNA and amplified by using PCR and sequenced by Illumina Genome Analyzer II.
Identification of miRNAs from sRNA dataset
Computational analysis of sRNA reads obtained from the sRNA libraries was performed as reported previously [57, 58]. Briefly, sRNA reads were extracted after trimming the adaptor sequence and then created an unique small RNA data set along with their read counts. The unique sRNAs were aligned to REPBASE version 14, (http://www.girinst.org), and other classes of known noncoding RNAs (rRNAs, tRNAs, snRNAs, snoRNAs, etc.) obtained from RFAM database (http://www.sanger.ac.uk/Software/Rfam/ftp.shtml). Following the removal of sRNAs corresponding to repetitive elements and known ncRNAs (rRNAs, tRNAs and other messenger RNAs), conserved miRNAs were identified by aligning the unique sRNA reads to the miRBase version 15 (http://microrna.sanger.ac.uk/). For identification of novel miRNAs, unique sRNAs with more than 10 genomic hits were removed from further analysis. The flanking regions of the remaining genome-matched sequences were extracted and the fold-back structures were predicted using the RNAfold program . Next, we examined the resulting folding structures to choose those that had at least 18 base-pairs, one central loop, and folding energy not greater than 18 Kcal/mol. The fold-back structures that meet the accepted criteria for miRNA annotation, ≤ 6 mismatches, ≤2 bulged or asymmetrically unpaired nucleotides, and ≤2 continuous mismatches in the putative miRNA region, were selected using MIRCHECK program . Subsequently, a program was developed in house to check the existence of miRNA* sequence of the selected mature miRNA, based on the criterion that there were 2nt overhang(s) at the 3’ end(s) of either the miRNA or miRNA*.
Mapping of sRNAs reads on to the sex chromosomes
The sRNA reads were aligned to the papaya sex chromosome pseudomolecules by using the ‘bowtie’ short read aligner. The bowtie option was set to obtain only the reads that match 100% to the reference. The SAM output from bowtie was converted binary file (BAM) using Samtools. The BAM output was then used to generate the alignment map using the Bamview program. For the conservation analysis between the pericentromeric regions of the sex chromosomes, we manually extracted the sRNA hotspots from each chromosomes based on the alignment maps generated above. Subsequently, the reads mapped to the Y chromosome sRNA hotspot on Y chromosome was realigned to the X chromosome hotspot and vice versa. The papaya sex specific small RNA sequence data used in this study can be obtained from NCBI's Gene Expression Ombibus (GEO) using accession number GSE54097.
We would like to thank Dr. Stephen Downie, Marija Pushko, and Will Wadlington for editing and providing feedbacks on the manuscript. This work was supported by a grant from the National Science Foundation (NSF) Plant Genome Research Program to RM and QY (Award Nos. DBI 0922545). RS acknowledges the support from the Oklahoma Agricultural Experiment Station.
- Axtell MJ, Snyder JA, Bartel DP: Common functions for diverse small RNAs of land plants. Plant Cell. 2007, 19 (6): 1750-1769. 10.1105/tpc.107.051706.PubMed CentralPubMedView ArticleGoogle Scholar
- Bartel B, Bartel DP: MicroRNAs: at the root of plant development?. Plant Physiol. 2003, 132 (2): 709-717. 10.1104/pp.103.023630.PubMed CentralPubMedView ArticleGoogle Scholar
- Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116: 281-297. 10.1016/S0092-8674(04)00045-5.PubMedView ArticleGoogle Scholar
- Reinhart BJ, Weinstein EG, Rhoades MW, Bartel B, Bartel DP: MicroRNAs in plants. Genes Dev. 2002, 16 (13): 1616-1626. 10.1101/gad.1004402.PubMed CentralPubMedView ArticleGoogle Scholar
- Rhoades MW, Reinhart BJ, Lim LP, Burge CB, Bartel B, Bartel DP: Prediction of plant microRNA targets. Cell. 2002, 110: 513-520. 10.1016/S0092-8674(02)00863-2.PubMedView ArticleGoogle Scholar
- Bernstein E, Allis CD: RNA meets chromatin. Genes Dev. 2005, 19 (14): 1635-1655. 10.1101/gad.1324305.PubMedView ArticleGoogle Scholar
- Cantu D, Vanzetti LS, Sumner A, Dubcovsky M, Matvienko M, Distelfeld A, Michelmore RW, Dubcovsky J: Small RNAs, DNA methylation and transposable elements in wheat. BMC Genomics. 2010, 11: 408-10.1186/1471-2164-11-408.PubMed CentralPubMedView ArticleGoogle Scholar
- Axtell MJ: Classification and comparison of small RNAs from plants. Annu Rev Plant Biol. 2013, 64: 137-159. 10.1146/annurev-arplant-050312-120043.PubMedView ArticleGoogle Scholar
- Sunkar R, Chinnusamy V, Zhu J, Zhu JK: Small RNAs as big players in plant abiotic stress responses and nutrient deprivation. Trends Plant Sci. 2007, 12 (7): 301-309. 10.1016/j.tplants.2007.05.001.PubMedView ArticleGoogle Scholar
- Grant-Downton R, Hafidh S, Twell D, Dickinson HG: Small RNA pathways are present and functional in the angiosperm male gametophyte. Mol Plant. 2009, 2 (3): 500-512. 10.1093/mp/ssp003.PubMedView ArticleGoogle Scholar
- Olmedo-Monfil V, Duran-Figueroa N, Arteaga-Vazquez M, Demesa-Arevalo E, Autran D, Grimanelli D, Slotkin RK, Martienssen RA, Vielle-Calzada J-P: Control of female gamete formation by a small RNA pathway in Arabidopsis. Nature. 2010, 464 (11): 628-634.PubMed CentralPubMedView ArticleGoogle Scholar
- 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 (4): 750-759. 10.1016/j.cell.2009.06.031.PubMed CentralPubMedView ArticleGoogle Scholar
- Kim JJ, Lee JH, Kim W, Jung HS, Huijser P, Ahn JH: The microRNA156-SQUAMOSA PROMOTER BINDING PROTEIN-LIKE3 module regulates ambient temperature-responsive flowering via FLOWERING LOCUS T in Arabidopsis. Plant Physiol. 2012, 159 (1): 461-478. 10.1104/pp.111.192369.PubMed CentralPubMedView ArticleGoogle Scholar
- Nodine MD, Bartel DP: MicroRNAs prevent precocious gene expression and enable pattern formation during plant embryogenesis. Genes Dev. 2010, 24 (23): 2678-2692. 10.1101/gad.1986710.PubMed CentralPubMedView ArticleGoogle Scholar
- Juarez MT, Kui JS, Thomas J, Heller BA, Timmermans MCP: microRNA-mediated repression of rolled leaf1 specifies maize leaf polarity. Nature. 2004, 428: 84-88. 10.1038/nature02363.PubMedView ArticleGoogle Scholar
- Achard P, Herr A, Baulcombe DC, Harberd NP: Modulation of floral development by a gibberellin-regulated microRNA. Development. 2004, 131 (14): 3357-3365. 10.1242/dev.01206.PubMedView ArticleGoogle Scholar
- Aukerman MJ, Sakai H: Regulation of flowering time and floral organ identity by microRNA and its Apatala2-like target gene. Plant Cell. 2003, 15 (11): 2730-2741. 10.1105/tpc.016238.PubMed CentralPubMedView ArticleGoogle Scholar
- Cartolano M, Castillo R, Efremova N, Kuckenberg M, Zethof J, Gerats T, Schwarz-Sommer Z, Vandenbussche M: A conserved microRNA module exerts homeotic control over Petunia hybrida and Antirrhinum majus floral organ identity. Nat Genet. 2007, 39 (7): 901-905. 10.1038/ng2056.PubMedView ArticleGoogle Scholar
- Chen X: A microRNA as a translational repressor of APETALA2 in Arabidopsis flower development. Science. 2004, 303 (5666): 2022-2025. 10.1126/science.1088060.PubMedView ArticleGoogle Scholar
- Yamaguchi A, Wu MF, Yang L, Wu G, Poethig RS, Wagner D: The microRNA-regulated SBP-Box transcription factor SPL3 is a direct upstream activator of LEAFY, FRUITFULL, and APETALA1. Dev Cell. 2009, 17 (2): 268-278. 10.1016/j.devcel.2009.06.007.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhao L, Kim Y, Dinh TT, Chen X: miR172 regulates stem cell fate and defines the inner boundary of APETALA3 and PISTILLATA expression domain in Arabidopsis floral meristems. Plant J. 2007, 51 (5): 840-849. 10.1111/j.1365-313X.2007.03181.x.PubMed CentralPubMedView ArticleGoogle Scholar
- Banks JA: MicroRNA, sex determination and floral meristem determinacy in maize. Genome Biol. 2008, 9 (1): 204-10.1186/gb-2008-9-1-204.PubMed CentralPubMedView ArticleGoogle Scholar
- 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 (12): 1517-1521. 10.1038/ng.2007.20.PubMedView ArticleGoogle Scholar
- She X, Xu X, Fedotov A, Kelly WG, Maine EM: Regulation of heterochromatin assembly on unpaired chromosomes during Caenorhabditis elegans meiosis by components of a small RNA-mediated pathway. PLoS Genet. 2009, 5 (8): e1000624-10.1371/journal.pgen.1000624.PubMed CentralPubMedView ArticleGoogle Scholar
- Carvalho FA, Renner SS: A dated phylogeny of the papaya family (Caricaceae) reveals the crop's closest relatives and the family's biogeographic history. Mol Phylogen Evol. 2012, 65 (1): 46-53. 10.1016/j.ympev.2012.05.019.View ArticleGoogle Scholar
- Storey WB: Genetics of the papaya. J Hered. 1953, 44: 70-78.Google Scholar
- Ming R, Yu Q, Moore PH: Sex determination in papaya. Semin Cell Dev Biol. 2007, 18 (3): 401-408. 10.1016/j.semcdb.2006.11.013.PubMedView ArticleGoogle Scholar
- Yu Q, Navajas-Pérez R, Tong E, Robertson J, Moore PH, Paterson AH, Ming R: Recent origin of dioecious and gynodioecious Y chromosomes in papaya. Tropical Plant Biology. 2008, 1 (1): 49-57. 10.1007/s12042-007-9005-7.View ArticleGoogle Scholar
- Zhang W, Wang X, Yu Q, Ming R, Jiang J: DNA methylation and heterochromatinization in the male-specific region of the primitive Y chromosome of papaya. Genome Res. 2008, 18 (12): 1938-1943. 10.1101/gr.078808.108.PubMed CentralPubMedView ArticleGoogle Scholar
- Wang J, Na J-K, Yu Q, Gschwend AR, Han J, Zeng F, Aryal R, VanBuren R, Murray JE, Zhang W, et al: Sequencing papaya X and Yh chromosomes revealed molecular basis of incipient sex chromosome evolution. Proc Natl Acad Sci USA. 2012, 109 (34): 13710-13715. 10.1073/pnas.1207833109.PubMed CentralPubMedView ArticleGoogle Scholar
- Gschwend AR, Yu Q, Tong EJ, Zeng F, Han J, VanBuren R, Aryal R, Charlesworth D, Moore PH, Paterson AH, et al: Rapid divergence and expansion of the X chromosome in papaya. Proc Natl Acad Sci USA. 2012, 109 (34): 13716-13721. 10.1073/pnas.1121096109.PubMed CentralPubMedView ArticleGoogle Scholar
- Aryal R, Yang X, Yu Q, Sunkar R, Li L, Ming R: Asymmetric purine-pyrimidine distribution in cellular small RNA population of papaya. BMC Genomics. 2012, 13: 682-10.1186/1471-2164-13-682.PubMed CentralPubMedView ArticleGoogle Scholar
- Aryal R, Ming R: Sex determination in flowering plants: papaya as a model system. Plant Sci. 2013,http://dx.doi.org/10.1016/j.plantsci.2013.10.018,Google Scholar
- Meyers BC, Axtell MJ, Bartel B, Bartel DP, Baulcombe D, Bowman JL, Cao X, Carrington JC, Chen X, Green PJ, et al: Criteria for annotation of plant MicroRNAs. Plant Cell. 2008, 20 (12): 3186-3190. 10.1105/tpc.108.064311.PubMed CentralPubMedView ArticleGoogle Scholar
- Porter BW, Aizawa KS, Zhu YJ, Christopher DA: Differentially expressed and new non-protein-coding genes from a Carica papaya root transcriptome survey. Plant Sci. 2008, 174 (1): 38-50. 10.1016/j.plantsci.2007.09.013.View ArticleGoogle Scholar
- Chen ZH, Bao ML, Sun YZ, Yang YJ, Xu XH, Wang JH, Han N, Bian HW, Zhu MY: Regulation of auxin response by miR393-targeted transport inhibitor response protein 1 is involved in normal development in Arabidopsis. Plant Mol Biol. 2011, 77 (6): 619-629. 10.1007/s11103-011-9838-1.PubMedView ArticleGoogle Scholar
- Mallory AC, Bartel DP, Bartel B: MicroRNA-directed regulation of Arabidopsis AUXIN RESPONSE FACTOR17 is essential for proper development and modulates expression of early auxin response genes. Plant Cell. 2005, 17 (5): 1360-1375. 10.1105/tpc.105.031716.PubMed CentralPubMedView ArticleGoogle Scholar
- Liu Q, Zhang YC, Wang CY, Luo YC, Huang QJ, Chen SY, Zhou H, Qu LH, Chen YQ: Expression analysis of phytohormone-regulated microRNAs in rice, implying their regulation roles in plant hormone signaling. FEBS Lett. 2009, 583 (4): 723-728. 10.1016/j.febslet.2009.01.020.PubMedView ArticleGoogle Scholar
- Na JK, Wang J, Murray JE, Gschwend AR, Zhang W, Yu Q, Navajas-Perez R, Feltus FA, Chen C, Kubat Z, et al: Construction of physical maps for the sex-specific regions of papaya sex chromosomes. BMC Genomics. 2012, 13 (1): 176-10.1186/1471-2164-13-176.PubMedView ArticleGoogle Scholar
- Reinhart BJ, Bartel DP: Small RNAs correspond to centromere heterochromatic repeats. Science. 2002, 297 (5588): 1831-10.1126/science.1077183.PubMedView ArticleGoogle Scholar
- Yu Q, Hou S, Hobza R, Feltus FA, Wang X, Jin W, Skelton RL, Blas A, Lemke C, Saw JH, et al: Chromosomal location and gene paucity of the male specific region on papaya Y chromosome. Mol Genet Genomics. 2007, 278 (2): 177-185. 10.1007/s00438-007-0243-z.PubMedView ArticleGoogle Scholar
- Ma H, Moore PH, Liu Z, Kim MS, Yu Q, Fitch MMM, Sekioka T, Paterson AH, Ming R: High-density linkage mapping revealed suppression of recombination at the sex determination locus in papaya. Genetics. 2004, 166: 419-436. 10.1534/genetics.166.1.419.PubMed CentralPubMedView ArticleGoogle Scholar
- Liu Z, Moore PH, Ma H, Ackerman CM, Ragiba M, Yu Q, Pearl HM, Kim MS, Charlton JW, Stiles JI, et al: A primitive Y chromosome in papaya marks incipient sex chromosome evolution. Nature. 2004, 427: 348-352. 10.1038/nature02228.PubMedView ArticleGoogle Scholar
- The International Brachypodium Initiative: Genome sequencing and analysis of the model grass Brachypodium distachyon. Nature. 2010, 463 (7282): 763-768. 10.1038/nature08747.View ArticleGoogle Scholar
- Lu C, Tej SS, Luo S, Haudenschild CD, Meyers BC, Green PJ: Elucidation of the small RNA component of the transcriptome. Science. 2005, 309 (5740): 1567-1569. 10.1126/science.1114112.PubMedView ArticleGoogle Scholar
- Lu C, Kulkarni K, Souret FF, MuthuValliappan R, Tej SS, Poethig RS, Henderson IR, Jacobsen SE, Wang W, Green PJ, et al: MicroRNAs and other small RNAs enriched in the Arabidopsis RNA-dependent RNA polymerase-2 mutant. Genome Res. 2006, 16 (10): 1276-1288. 10.1101/gr.5530106.PubMed CentralPubMedView ArticleGoogle Scholar
- Hayden KE, Willard HF: Composition and organization of active centromere sequences in complex genomes. BMC Genomics. 2012, 13: 324-10.1186/1471-2164-13-324.PubMed CentralPubMedView ArticleGoogle Scholar
- Nagaki K, Cheng Z, Ouyang S, Talbert PB, Kim M, Jones KM, Henikoff S, Buell CR, Jiang J: Sequencing of a rice centromere uncovers active genes. Nat Genet. 2004, 36 (2): 138-145. 10.1038/ng1289.PubMedView ArticleGoogle Scholar
- Kumar LSS, Abraham A, Srinivasan VK: The cytology of Carica papaya Linn. Indian J Agr SCi. 1945, 15: 242-253.Google Scholar
- Liu Q, Chen YQ: Insights into the mechanism of plant development: interactions of miRNAs pathway with phytohormone response. Biochem Biophys Res Commun. 2009, 384 (1): 1-5. 10.1016/j.bbrc.2009.04.028.PubMedView ArticleGoogle Scholar
- Larsson E, Franks RG, Sundberg E: Auxin and the Arabidopsis thaliana gynoecium. J Exp Bot. 2013, 64 (9): 2619-2627. 10.1093/jxb/ert099.PubMedView ArticleGoogle Scholar
- Geuten K, Irish V: Hidden variability of floral homeotic B genes in Solanaceae provides a molecular basis for the evolution of novel functions. Plant Cell. 2010, 22 (8): 2562-2578. 10.1105/tpc.110.076026.PubMed CentralPubMedView ArticleGoogle Scholar
- 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 (12): 1035-1046. 10.1016/j.cub.2004.06.022.PubMedView ArticleGoogle Scholar
- Jung JH, Park CM: MIR166/165 genes exhibit dynamic expression patterns in regulating shoot apical meristem and floral development in Arabidopsis. Planta. 2007, 225 (6): 1327-1338. 10.1007/s00425-006-0439-1.PubMedView ArticleGoogle Scholar
- Knauer S, Holt AL, Rubio-Somoza I, Tucker EJ, Hinze A, Pisch M, Javelle M, Timmermans MC, Tucker MR, Laux T: A protodermal miR394 signal defines a region of stem cell competence in the Arabidopsis shoot meristem. Dev Cell. 2013, 24 (2): 125-132. 10.1016/j.devcel.2012.12.009.PubMedView ArticleGoogle Scholar
- Sunkar R, Zhou X, Zheng Y, Zhang W, Zhu JK: Identification of novel and candidate miRNAs in rice by high throughput sequencing. BMC Plant Biol. 2008, 8: 25-10.1186/1471-2229-8-25.PubMed CentralPubMedView ArticleGoogle Scholar
- Jagadeeswaran G, Nimmakayala P, Zheng Y, Gowdu K, Reddy UK, Sunkar R: Characterization of the small RNA component of leaves and fruits from four different cucurbit species. BMC Genomics. 2012, 13: 329-10.1186/1471-2164-13-329.PubMed CentralPubMedView ArticleGoogle Scholar
- Jagadeeswaran G, Zheng Y, Sumathipala N, Jiang H, Arrese EL, Soulages JL, Zhang W, Sunkar R: Deep sequencing of small RNA libraries reveals dynamic regulation of conserved and novel microRNAs and microRNA-stars during silkworm development. BMC Genomics. 2010, 11: 52-10.1186/1471-2164-11-52.PubMed CentralPubMedView ArticleGoogle Scholar
- Hofacker IL: Vienna RNA secondary structure server. Nucleic Acids Res. 2003, 31 (13): 3429-3431. 10.1093/nar/gkg599.PubMed CentralPubMedView ArticleGoogle Scholar
- Jones-Rhoades MW, Bartel DP: Computational identification of plant microRNAs and their targets, including a stress-induced miRNA. Mol Cell. 2004, 14 (6): 787-799. 10.1016/j.molcel.2004.05.027.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. 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.