Sex specific expression and distribution of small RNAs in papaya

Complementing model species with model clades

Model species continue to underpin groundbreaking plant science research. At the same time, the phylogenetic resolution of the land plant tree of life continues to improve. The intersection of these 2 research paths creates a unique opportunity to further extend the usefulness of model species across larger taxonomic groups. Here we promote the utility of the Arabidopsis thaliana model species, especially the ability to connect its genetic and functional resources, to species across the entire Brassicales order. We focus on the utility of using genomics and phylogenomics to bridge the evolution and diversification of several traits across the Brassicales to the resources in Arabidopsis, thereby extending scope from a model species by establishing a "model clade." These Brassicales-wide traits are discussed in the context of both the model species Arabidopsis and the family Brassicaceae. We promote the utility of such a "model clade" and make suggestions for building global networks to support future studies in the model order Brassicales.

Combined analysis of mRNA and miRNA transcriptomes reveals the regulatory mechanism of Xanthomonas arboricola pv pruni resistance in Prunus persica

BACKGROUND

Peach bacterial shot hole, caused by Xanthomonas arboricola pv pruni (Xap), is a global bacterial disease that poses a threat to the yield and quality of cultivated peach trees (Prunus persica).

RESULTS

This study compared the mRNA and miRNA profiles of two peach varieties, 'Yanbao' (resistant) and 'Yingzui' (susceptible), after inoculation with Xap to identify miRNAs and target genes associated with peach tree resistance. mRNA sequencing results revealed that in the S0-vs-S3 comparison group, 1574 genes were upregulated and 3975 genes were downregulated. In the R0-vs-R3 comparison group, 1575 genes were upregulated and 3726 genes were downregulated. Through miRNA sequencing, a total of 112 known miRNAs belonging to 70 miRNA families and 111 new miRNAs were identified. Notably, some miRNAs were exclusively expressed in either resistant or susceptible varieties. Additionally, 59 miRNAs were downregulated and 69 miRNAs were upregulated in the R0-vs-R3 comparison group, while 46 miRNAs were downregulated and 52 miRNAs were upregulated in the S0-vs-S3 comparison group. Joint analysis of mRNA and miRNA identified 79 relationship pairs in the S0-vs-S3 comparison group, consisting of 48 miRNAs and 51 target genes. In the R0-vs-R3 comparison group, there were 58 relationship pairs, comprising 28 miRNAs and 20 target genes. Several target genes related to resistance, such as SPL6, TIFY6B, and Prupe.4G041800_v2.0.a1 (PPO), were identified through literature reports and GO/KEGG enrichment analysis.

CONCLUSION

In conclusion, this study discovered several candidate genes involved in peach tree resistance by analyzing differential expression of mRNA and miRNA. These findings provide valuable insights into the mechanisms underlying resistance to Xap in peach trees.

QTL analysis of femaleness in monoecious spinach and fine mapping of a major QTL using an updated version of chromosome-scale pseudomolecules

Although spinach is predominantly dioecious, monoecious plants with varying proportions of female and male flowers are also present. Recently, monoecious inbred lines with highly female and male conditions have been preferentially used as parents for F1-hybrids, rather than dioecious lines. Accordingly, identifying the loci for monoecism is an important issue for spinach breeding. We here used long-read sequencing and Hi-C technology to construct SOL_r2.0_pseudomolecule, a set of six pseudomolecules of spinach chromosomes (total length: 879.2 Mb; BUSCO complete 97.0%) that are longer and more genetically complete than our previous version of pseudomolecules (688.0 Mb; 81.5%). Three QTLs, qFem2.1, qFem3.1, and qFem6.1, responsible for monoecism were mapped to SOL_r2.0_pseudomolecule. qFem3.1 had the highest LOD score and corresponded to the M locus, which was previously identified as a determinant of monoecious expression, by genetic analysis of progeny from female and monoecious plants. The other QTLs were shown to modulate the ratio of female to male flowers in monoecious plants harboring a dominant allele of the M gene. Our findings will enable breeders to efficiently produce highly female- and male-monoecious parental lines for F1-hybrids by pyramiding the three QTLs. Through fine-mapping, we narrowed the candidate region for the M locus to a 19.5 kb interval containing three protein-coding genes and one long non-coding RNA gene. Among them, only RADIALIS-like-2a showed a higher expression in the reproductive organs, suggesting that it might play a role in reproductive organogenesis. However, there is no evidence that it is involved in the regulation of stamen and pistil initiation, which are directly related to the floral sex differentiation system in spinach. Given that auxin is involved in reproductive organ formation in many plant species, genes related to auxin transport/response, in addition to floral organ formation, were identified as candidates for regulators of floral sex-differentiation from qFem2.1 and qFem6.1.

Sex determination in papaya: Current status and perspectives

Papaya (Carica papaya L.) is an economically significant plant that produces fruit consumed worldwide due to its organoleptic characteristics. Since their commercial production, papaya fruits have faced several problems, such as pests, which have been partly resolved using transgenic varieties. Nevertheless, a principal challenge in this cultivation is the plant's sex determination. The sex issue in papaya is complex because papaya flowers can bear three sex forms: male, female, and hermaphrodite, which affects their fruit production, shape, and yield. Fruits from hermaphrodite plants are preferred more by consumers than female ones, and male plants rarely produce fruits without commercial value. Chromosomes are responsible for sex determination in papaya, denoted as XY for male, XX for female, and XYh for hermaphrodite. However, genes related to sex have been reported but are not conclusive. Factors such as the environment, hormones, and genetic and epigenetic background can also affect sex expression. Therefore, in this review, we will discuss recent research on the sex of papaya, from reported genes to date, their biology, and sexing approaches using molecular markers and their advantages.

Identification of miRNAs and Target Genes at Key Stages of Sexual Differentiation in Androdioecious Osmanthus fragrans

Androdioecy is the crucial transition state in the evolutionary direction of hermaphroditism to dioecy, however, the molecular mechanisms underlying the formation of this sex system remain unclear. While popular in China for its ornamental and cultural value, Osmanthus fragrans has an extremely rare androdioecy breeding system, meaning that there are both male and hermaphroditic plants in a population. To unravel the mechanisms underlying the formation of androdioecy, we performed small RNA sequencing studies on male and hermaphroditic O. fragrans. A total of 334 miRNAs were identified, of which 59 were differentially expressed. Functional categorization revealed that the target genes of differentially expressed miRNAs were mainly involved in the biological processes of reproductive development and the hormone signal transduction pathway. We speculated that the miRNA160, miRNA167, miRNA393 and miRNA396 families may influence the sex differentiation in O. fragrans. Overall, our study is the first exploration of miRNAs in the growth and development process of O. fragrans, and is also the first study of androdioecious plants from the miRNA sequencing perspective. The analysis of miRNAs and target genes that may be involved in the sex differentiation process lay a foundation for the ultimate discovery of the androdioecious molecular mechanism in O. fragrans.

Identification and verification of seed development related miRNAs in kernel almond by small RNA sequencing and qPCR

Many studies have investigated the role of miRNAs on the yield of various plants, but so far, no report is available on the identification and role of miRNAs in fruit and seed development of almonds. In this study, preliminary analysis by high-throughput sequencing of short RNAs of kernels from the crosses between almond cultivars 'Sefid' × 'Mamaee' (with small and large kernels, respectively) and 'Sefid' × 'P. orientalis' (with small kernels) showed that the expressions of several miRNAs such as Pdu-miR395a-3p, Pdu-miR8123-5p, Pdu-miR482f, Pdu-miR6285, and Pdu-miR396a were significantly different. These miRNAs targeted genes encoding different proteins such as NYFB-3, SPX1, PGSIP3 (GUX2), GH3.9, and BEN1. The result of RT-qPCR revealed that the expression of these genes showed significant differences between the crosses and developmental stages of the seeds, suggesting that these genes might be involved in controlling kernel size because the presence of these miRNAs had a negative effect on their target genes. Pollen source can influence kernel size by affecting hormonal signaling and metabolic pathways through related miRNAs, a phenomenon known as xenia.

Transcriptome-based analysis of the hormone regulation mechanism of gender differentiation in Juglans mandshurica Maxim

Juglans mandshurica Maxim is a hermaphroditic plant belonging to the genus Juglans in the family Juglandaceae. The pollination period of female flowers is different from the loose powder period of male flowers on the same tree. In several trees, female flowers bloom first, whereas in others, male flowers bloom first. In this study, male and female flower buds of J. mandshurica at the physiological differentiation stage were used. Illumina-based transcriptome sequencing was performed, and the quality of the sequencing results was evaluated and analyzed. A total of 138,138 unigenes with an average length of 788 bp were obtained. There were 8,116 differentially expressed genes (DEGs); 2,840 genes were upregulated, and 5,276 genes were downregulated. The DEGs were classified by Gene Ontology and analyzed by Kyoto Encyclopedia of Genes and Genomes. The signal transduction factors involved in phytohormone synthesis were selected. The results displayed that ARF and SAUR were expressed differently in the auxin signaling pathway. Additionally, DELLA protein (a negative regulator of gibberellin), the cytokinin synthesis pathway, and A-ARR were downregulated. On April 2nd, the contents of IAA, GA, CTK, ETH and SA in male and female flower buds of two types of J. mandshurica were opposite, and there were obvious genes regulating gender differentiation. Overall, we found that the sex differentiation of J. mandshurica was related to various hormone signal transduction pathways, and hormone signal transduction plays a leading role in regulation.

Sex biased expression of hormone related genes at early stage of sex differentiation in papaya flowers

Sex types of papaya are controlled by a pair of nascent sex chromosomes, but molecular genetic mechanisms of sex determination and sex differentiation in papaya are still unclear. We performed comparative analysis of transcriptomic profiles of male and female floral buds at the early development stage before the initiation of reproductive organ primordia at which there is no morphological difference between male and female flowers. A total of 1734 differentially expressed genes (DEGs) were identified, of which 923 showed female-biased expression and 811 showed male-biased expression. Functional annotation revealed that genes related to plant hormone biosynthesis and signaling pathways, especially in abscisic acid and auxin pathways, were overrepresented in the DEGs. Transcription factor binding motifs, such as MYB2, GAMYB, and AP2/EREBP, were enriched in the promoters of the hormone-related DEGs, and transcription factors with those motifs also exhibited differential expression between sex types. Among these DEGs, we also identified 11 genes in the non-recombining region of the papaya sex chromosomes and 9 genes involved in stamen and carpel development. Our results suggested that sex differentiation in papaya may be regulated by multiple layers of regulation and coordination and involved transcriptional, epigenetic, and phytohormone regulation. Hormones, especially ABA and auxin, transcription factors, and genes in the non-recombination region of the sex chromosome could be involved in this process. Our findings may facilitate the elucidation of signal transduction and gene interaction in sex differentiation of unisexual flowers in papaya.

Sex and the flower - developmental aspects of sex chromosome evolution

Background

The evolution of dioecious plants is occasionally accompanied by the establishment of sex chromosomes: both XY and ZW systems have been found in plants. Structural studies of sex chromosomes are now being followed up by functional studies that are gradually shedding light on the specific genetic and epigenetic processes that shape the development of separate sexes in plants.

Scope

This review describes sex determination diversity in plants and the genetic background of dioecy, summarizes recent progress in the investigation of both classical and emerging model dioecious plants and discusses novel findings. The advantages of interspecies hybrids in studies focused on sex determination and the role of epigenetic processes in sexual development are also overviewed.

Conclusions

We integrate the genic, genomic and epigenetic levels of sex determination and stress the impact of sex chromosome evolution on structural and functional aspects of plant sexual development. We also discuss the impact of dioecy and sex chromosomes on genome structure and expression.

Small RNA and degradome profiling involved in seed development and oil synthesis of Brassica napus

MicroRNAs (miRNAs) play a prominent role in post-transcriptional gene expression regulation and have been involved in various biological and metabolic processes to regulate gene expression. For Brassica napus, improving seed-weight and oil-content is the main breeding goal. In order to better understand the regulation mechanism of miRNAs during seed-weight formation and oil-content accumulation in B. napus, in this study, a high-throughput sequencing technology was used to profile miRNAs expression of Brassica napus immature seeds from one to six weeks after flowering. A total of 1,276 miRNAs, including 1,248 novel and 28 known miRNAs, were obtained from both the high-seed-weight with low-oil-content RNA pool (S03) and the low-seed-weight with high-oil-content RNA pool (S04). Analysis of their expression profiles disclosed that 300 novel and two known miRNAs were differentially expressed between S03 and S04. For degradome analysis, 57 genes with 64 degradation sites were predicted to be targeted for degradation by these miRNAs. Further bioinformatics analysis indicated that these differentially expressed miRNAs might participate in regulation of myriad cellular and molecular processes, during seed development and oil synthesis. Finally, 6 target genes with potential roles in regulation of seed development and 9 other targets in seed oil synthesis, were further confirmed as candidate genes from small RNA and degradome sequencing.

The development of functional mapping by three sex-related loci on the third whorl of different sex types of Carica papaya L

Carica papaya L. is an important economic crop worldwide and is used as a model plant for sex-determination research. To study the different flower sex types, we screened sex-related genes using alternative splicing sequences (AS-seqs) from a transcriptome database of the three flower sex types, i.e., males, females, and hermaphrodites, established at 28 days before flowering using 15 bacterial artificial chromosomes (BACs) of C. papaya L. After screening, the cDNA regions of the three sex-related loci, including short vegetative phase-like (CpSVPL), the chromatin assembly factor 1 subunit A-like (CpCAF1AL), and the somatic embryogenesis receptor kinase (CpSERK), which contained eight sex-related single-nucleotide polymorphisms (SNPs) from the different sex types of C. papaya L., were genotyped using high-resolution melting (HRM). The three loci were examined regarding the profiles of the third whorl, as described below. CpSVPL, which had one SNP associated with the three sex genotypes, was highly expressed in the male and female sterile flowers (abnormal hermaphrodite flowers) that lacked the fourth whorl structure. CpCAF1AL, which had three SNPs associated with the male genotype, was highly expressed in male and normal hermaphrodite flowers, and had no AS-seqs, whereas it exhibited low expression and an AS-seqs in intron 11 in abnormal hermaphrodite flowers. Conversely, carpellate flowers (abnormal hermaphrodite flowers) showed low expression of CpSVPL and AS-seqs in introns 5, 6, and 7 of CpSERK, which contained four SNPs associated with the female genotype. Specifically, the CpSERK and CpCAF1AL loci exhibited no AS-seq expression in the third whorl of the male and normal hermaphrodite flowers, respectively, and variance in the AS-seq expression of all other types of flowers. Functional mapping of the third whorl of normal hermaphrodites indicated no AS-seq expression in CpSERK, low CpSVPL expression, and, for CpCAF1AL, high expression and no AS-seq expression on XY^h-type chromosomes.

Genome-wide identification of Hami melon miRNAs with putative roles during fruit development

MicroRNAs represent a family of small endogenous, non-coding RNAs that play critical regulatory roles in plant growth, development, and environmental stress responses. Hami melon is famous for its attractive flavor and excellent nutritional value, however, the mechanisms underlying the fruit development and ripening remains largely unknown. Here, we performed small RNA sequencing to investigate the roles of miRNAs during Hami melon fruit development. Two batches of flesh samples were collected at four fruit development stages. Small RNA sequencing yielded a total of 54,553,424 raw reads from eight libraries. 113 conserved miRNAs belonging to 30 miRNA families and nine novel miRNAs comprising nine miRNA families were identified. The expression of 42 conserved miRNAs and three Hami melon-specific miRNAs significantly changed during fruit development. Furthermore, 484 and 124 melon genes were predicted as putative targets of 29 conserved and nine Hami melon-specific miRNA families, respectively. GO enrichment analysis were performed on target genes, "transcription, DNA-dependent", "rRNA processing", "oxidation reduction", "signal transduction", "regulation of transcription, DNA-dependent", and "metabolic process" were the over-represented biological process terms. Cleavage sites of six target genes were validated using 5' RACE. Our results present a comprehensive set of identification and characterization of Hami melon fruit miRNAs and their potential targets, which provide valuable basis towards understanding the regulatory mechanisms in programmed process of normal Hami fruit development and ripening. Specific miRNAs could be selected for further research and applications in breeding practices.

Overexpression of a flower-specific aerolysin-like protein from the dioecious plant Rumex acetosa alters flower development and induces male sterility in transgenic tobacco

Sex determination in Rumex acetosa, a dioecious plant with a complex XY1 Y2 sex chromosome system (females are XX and males are XY1 Y2 ), is not controlled by an active Y chromosome but depends on the ratio between the number of X chromosomes and autosomes. To gain insight into the molecular mechanisms of sex determination, we generated a subtracted cDNA library enriched in genes specifically or predominantly expressed in female floral buds in early stages of development, when sex determination mechanisms come into play. In the present paper, we report the molecular and functional characterization of FEM32, a gene encoding a protein that shares a common architecture with proteins in different plants, animals, bacteria and fungi of the aerolysin superfamily; many of these function as β pore-forming toxins. The expression analysis, assessed by northern blot, RT-PCR and in situ hybridization, demonstrates that this gene is specifically expressed in flowers in both early and late stages of development, although its transcripts accumulate much more in female flowers than in male flowers. The ectopic expression of FEM32 under both the constitutive promoter 35S and the flower-specific promoter AP3 in transgenic tobacco showed no obvious alteration in vegetative development but was able to alter floral organ growth and pollen fertility. The 35S::FEM32 and AP3::FEM32 transgenic lines showed a reduction in stamen development and pollen viability, as well as a diminution in fruit set, fruit development and seed production. Compared with other floral organs, pistil development was, however, enhanced in plants overexpressing FEM32. According to these effects, it is likely that FEM32 functions in Rumex by arresting stamen and pollen development during female flower development. The aerolysin-like pore-forming proteins of eukaryotes are mainly involved in defence mechanisms against bacteria, fungi and insects and are also involved in apoptosis and programmed cell death (PCD), a mechanism that could explain the role of FEM32 in Rumex sex determination.

A Century of Sex Determination in Flowering Plants

Plants have evolved a diverse array of strategies for sexual reproduction, particularly through the modification of male and female organs at distinct points in development. The immense variation in sexual systems across the land plants provides a unique opportunity to study the genetic, epigenetic, phylogenetic, and ecological underpinnings of sex determination. Here, we reflect on more than a century of research into flowering plant sex determination, placing a particular focus on the foundational genetic and cytogenetic observations, experiments, and hypotheses. Building on the seminal work on the genetics of plant sex, modern comparative genomic analyses now allow us to address longstanding questions about sex determination and the origins of sex chromosomes.

Structure-based discovery of potentially active semiochemicals for Cydia pomonella (L.)

The development of physiologically active semiochemicals is largely limited by the labor-consuming searching process. How to screen active semiochemicals efficiently is of significance to the extension of behavior regulation in pest control. Here pharmacophore modeling and shape-based virtual screening were combined to predict candidate ligands for Cydia pomonella pheromone binding protein 1 (CpomPBP1). Out of the predicted compounds, ETrME displayed the highest affinity to CpomPBP1. Further studies on the interaction between CpomPBP1 and ETrME, not only depicted the binding mode, but also revealed residues providing negative and positive contributions to the ETrME binding. Moreover, key residues involved in interacting with ETrME of CpomPBP1 were determined as well. These findings were significant to providing insights for the future searching and optimization of active semiochemicals.

Genomic identification, rapid evolution, and expression of Argonaute genes in the tilapia, Oreochromis niloticus

Argonaute proteins are key components of the small RNA-induced silencing complex and have multiple roles in RNA-directed regulatory pathways. Argonaute genes can be divided into two subfamilies: the Ago (interacting with microRNA/small interfering RNA) and Piwi subfamilies (interacting with piwi-interacting RNAs (piRNAs)). In the present study, genome-wide analyses firstly yielded the identification of different members of Agos and Piwis in the tilapia, coelacanth, spotted gar, and elephant shark. The additional teleost Ago3b was generated following the fish-specific genome duplication event. Selective pressure analysis on Agos and Piwis between cichlids and other teleosts showed an accelerated evolution of Piwil1 in the cichlid lineages, and the positive selected sites were located in the region of PIWI domain, suggesting that these amino acid substitutions are adapt to targeted cleavage of messenger RNA (mRNA) in cichlids. Ago1 and Ago4 were detected at higher levels at 5 days after hatching (dah) in both ovaries and testes compared with other stages, supporting the previously reported requirement of Ago-mediated pathways to clear the maternal mRNAs during the early embryogenesis. The Piwis were abundantly expressed in tilapia testes, indicating their essential roles in male germline, especially in spermatogenesis. Notable expression of Piwis was also detected in skeletal muscle, indicating that piRNA pathway may not only be confined to development and maintenance of the germline but may also play important roles in somatic tissues. The expression of Piwil1 and Piwil2 was examined by quantitative PCR (qPCR) and in situ hybridization (ISH) to validate the spatial and temporal expression profiles. Taken together, these results present a thorough overview of tilapia Argonaute family and provide a new perspective on the evolution and function of this family in teleosts.

Novel and conserved microRNAs in soybean floral whorls

MicroRNAs (miRNAs) correspond to a class of endogenous small non-coding RNAs (19-24 nt) that regulates the gene expression, through mRNA target cleavage or translation inhibition. In plants, miRNAs have been shown to play pivotal roles in a wide variety of metabolic and biological processes like plant growth, development, and response to biotic and abiotic stress. Soybean is one of the most important crops worldwide, due to the production of oil and its high protein content. The reproductive phase is considered the most important for soybean yield, which is mainly intended to produce the grains. The identification of miRNAs is not yet saturated in soybean, and there are no studies linking them to the different floral organs. In this study, three different mature soybean floral whorls were used in the construction of sRNA libraries. The sequencing of petal, carpel and stamen libraries generated a total of 10,165,661 sequences. Subsequent analyses identified 200 miRNAs sequences, among which, 41 were novel miRNAs, 80 were conserved soybean miRNAs, 31 were new antisense conserved soybean miRNAs and 46 were soybean miRNAs isoforms. We also found a new miRNA conserved in other plant species, and finally one miRNA-sibling of a soybean conserved miRNA. Conserved and novel miRNAs were evaluated by RT-qPCR. We observed a differential expression across the three whorls for six miRNAs. Computational predicted targets for miRNAs analyzed by RT-qPCR were identified and present functions related to reproductive process in plants. In summary, the increased accumulation of specific and novel miRNAs in different whorls indicates that miRNAs are an important part of the regulatory network in soybean flower.

Identification of miRNAs with potential roles in regulation of anther development and male-sterility in 7B-1 male-sterile tomato mutant

BACKGROUND

The 7B-1 tomato line (Solanum lycopersicum cv. Rutgers) is a photoperiod-sensitive male-sterile mutant, with potential application in hybrid seed production. Small RNAs (sRNAs) in tomato have been mainly characterized in fruit development and ripening, but none have been studied with respect to flower development and regulation of male-sterility. Using sRNA sequencing, we identified miRNAs that are potentially involved in anther development and regulation of male-sterility in 7B-1 mutant.

RESULTS

Two sRNA libraries from 7B-1 and wild type (WT) anthers were sequenced and thirty two families of known miRNAs and 23 new miRNAs were identified in both libraries. MiR390, miR166, miR159 were up-regulated and miR530, miR167, miR164, miR396, miR168, miR393, miR8006 and two new miRNAs, miR#W and miR#M were down-regulated in 7B-1 anthers. Ta-siRNAs were not differentially expressed and likely not associated with 7B-1 male-sterility. miRNA targets with potential roles in anther development were validated using 5'-RACE. QPCR analysis showed differential expression of miRNA/target pairs of interest in anthers and stem of 7B-1, suggesting that they may regulate different biological processes in these tissues. Expression level of most miRNA/target pairs showed negative correlation, except for few. In situ hybridization showed predominant expression of miR159, GAMYBL1, PMEI and cystatin in tapetum, tetrads and microspores.

CONCLUSION

Overall, we identified miRNAs with potential roles in anther development and regulation of male-sterility in 7B-1. A number of new miRNAs were also identified from tomato for the first time. Our data could be used as a benchmark for future studies of the molecular mechanisms of male-sterility in other crops.

microRNAs and the evolution of complex multicellularity: identification of a large, diverse complement of microRNAs in the brown alga Ectocarpus

There is currently convincing evidence that microRNAs have evolved independently in at least six different eukaryotic lineages: animals, land plants, chlorophyte green algae, demosponges, slime molds and brown algae. MicroRNAs from different lineages are not homologous but some structural features are strongly conserved across the eukaryotic tree allowing the application of stringent criteria to identify novel microRNA loci. A large set of 63 microRNA families was identified in the brown alga Ectocarpus based on mapping of RNA-seq data and nine microRNAs were confirmed by northern blotting. The Ectocarpus microRNAs are highly diverse at the sequence level with few multi-gene families, and do not tend to occur in clusters but exhibit some highly conserved structural features such as the presence of a uracil at the first residue. No homologues of Ectocarpus microRNAs were found in other stramenopile genomes indicating that they emerged late in stramenopile evolution and are perhaps specific to the brown algae. The large number of microRNA loci in Ectocarpus is consistent with the developmental complexity of many brown algal species and supports a proposed link between the emergence and expansion of microRNA regulatory systems and the evolution of complex multicellularity.

Genome-wide comparison of microRNAs and their targeted transcripts among leaf, flower and fruit of sweet orange

BACKGROUND

In plants, microRNAs (miRNAs) regulate gene expression mainly at the post-transcriptional level. Previous studies have demonstrated that miRNA-mediated gene silencing pathways play vital roles in plant development. Here, we used a high-throughput sequencing approach to characterize the miRNAs and their targeted transcripts in the leaf, flower and fruit of sweet orange.

RESULTS

A total of 183 known miRNAs and 38 novel miRNAs were identified. An in-house script was used to identify all potential secondary siRNAs derived from miRNA-targeted transcripts using sRNA and degradome sequencing data. Genome mapping revealed that these miRNAs were evenly distributed across the genome with several small clusters, and 69 pre-miRNAs were co-localized with simple sequence repeats (SSRs). Noticeably, the loop size of pre-miR396c was influenced by the repeat number of CUU unit. The expression pattern of miRNAs among different tissues and developmental stages were further investigated by both qRT-PCR and RNA gel blotting. Interestingly, Csi-miR164 was highly expressed in fruit ripening stage, and was validated to target a NAC transcription factor. This study depicts a global picture of miRNAs and their target genes in the genome of sweet orange, and focused on the comparison among leaf, flower and fruit tissues.

CONCLUSIONS

This study provides a global view of miRNAs and their target genes in different tissue of sweet orange, and focused on the identification of miRNA involved in the regulation of fruit ripening. The results of this study lay a foundation for unraveling key regulators of orange fruit development and ripening on post-transcriptional level.