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

Negative regulators of grain yield and mineral contents in rice: potential targets for CRISPR-Cas9-mediated genome editing

Rice is a major global staple food crop, and improving its grain yield and nutritional quality has been a major thrust research area since last decades. Yield and nutritional quality are complex traits which are controlled by multiple signaling pathways. Sincere efforts during past decades of research have identified several key genetic and molecular regulators that governed these complex traits. The advent of clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated protein 9 (Cas9)-mediated gene knockout approaches has accelerated the development of improved varieties; however, finding out target gene with negative regulatory function in particular trait without giving any pleiotropic effect remains a challenge. Here, we have reviewed past and recent literature and identified important negative regulators of grain yield and mineral contents which could be potential targets for CRISPR-Cas9-mediated gene knockout. Additionally, we have also compiled a list of microRNAs (miRNAs), which target positive regulators of grain yield, plant stress tolerance, and grain mineral contents. Knocking out these miRNAs could help to increase expression of such positive regulators and thus improve the plant trait. The knowledge presented in this review would help to further accelerate the CRISPR-Cas9-mediated trait improvement in rice.

Functional role of microRNA in the regulation of biotic and abiotic stress in agronomic plants

The increasing demand for food is the result of an increasing population. It is crucial to enhance crop yield for sustainable production. Recently, microRNAs (miRNAs) have gained importance because of their involvement in crop productivity by regulating gene transcription in numerous biological processes, such as growth, development and abiotic and biotic stresses. miRNAs are small, non-coding RNA involved in numerous other biological functions in a plant that range from genomic integrity, metabolism, growth, and development to environmental stress response, which collectively influence the agronomic traits of the crop species. Additionally, miRNA families associated with various agronomic properties are conserved across diverse plant species. The miRNA adaptive responses enhance the plants to survive environmental stresses, such as drought, salinity, cold, and heat conditions, as well as biotic stresses, such as pathogens and insect pests. Thus, understanding the detailed mechanism of the potential response of miRNAs during stress response is necessary to promote the agronomic traits of crops. In this review, we updated the details of the functional aspects of miRNAs as potential regulators of various stress-related responses in agronomic plants.

Characterization and identification of grapevine heat stress-responsive microRNAs revealed the positive regulated function of vvi-miR167 in thermostability

High temperature stress is one of the primary abiotic stresses that restrict fruit tree production. Grapevine (Vitis vinifera) with high economic value throughout the world is a cultivated fruit crop, and its growth and development is often influenced by high temperature stress. Studying the heat stress-response mechanism of grapevine has great significance for understanding the acclimation to heat stress. In this study, we identified a series of heat stress responsive miRNAs and analyzed their function during the heat tolerance response. CK (control group, 25 °C) and heat treatment stress (TS, 45 °C) small RNA (sRNA) libraries were constructed and sequenced by high-throughput sequencing in 'Thompson seedless' grapevine. 873 known-miRNAs and 86 novel-miRNAs were identified, of which 88 known and three novel miRNAs were expressed differentially under heat stress. 322 genes were predicted to be targeted by the miRNAs. Eight selected miRNAs and its targets were confirmed by real time quantitative PCR (RT - qPCR), indicating that these "miRNA - target" were responsive to heat stress. In addition, most of the predicted target genes were negatively regulated by corresponding miRNAs. Gene function and pathway analyses indicated that these genes probably play crucial roles in heat stress tolerance. Vvi-miR167b transiently overexpression in grapevine leaves decreased target gene vvARF6, vvARF6-like and vvARF8 expression. The function of vvi-miR167 was verified by ectopic transformation in Arabidopsis thaliana, and the heat tolerance in transgenic lines was enhanced significantly, suggesting that the vvi-miR167 plays a positive regulatory role in grape thermostability. Comparison of miRNA expression patterns between heat treatment stress and CK can help elucidate the heat stress response and resistance mechanisms in grapes. In conclusion, these results gave us useful information to better understand the heat stress-response during domestication as well as for breeding new cultivars with heat stress resistance in fruit trees.

What Do We Know about Barley miRNAs?

Plant miRNAs are powerful regulators of gene expression at the post-transcriptional level, which was repeatedly proved in several model plant species. miRNAs are considered to be key regulators of many developmental, homeostatic, and immune processes in plants. However, our understanding of plant miRNAs is still limited, despite the fact that an increasing number of studies have appeared. This systematic review aims to summarize our current knowledge about miRNAs in spring barley (Hordeum vulgare), which is an important agronomical crop worldwide and serves as a common monocot model for studying abiotic stress responses as well. This can help us to understand the connection between plant miRNAs and (not only) abiotic stresses in general. In the end, some future perspectives and open questions are summarized.

miR171 modulates induction of somatic embryogenesis in citrus callus

Overexpression of miR171 restored SE competence in the recalcitrant citrus callus, and inhibition of miR171 function weakened SE competence in the strongly embryogenic citrus callus. Somatic embryogenesis (SE) is an important way of in vitro regeneration for plants. For perennial woody crops such as citrus, embryogenic callus is usually induced from unfertilized aborted ovules and widely used in biotechnology aided breeding. However, SE capacity always declines in callus during subculture, which makes regeneration difficult and hinders the application of biotechnology. We previously found that miR171 may be a regulator of SE in citrus, based on the abundant expression of csi-miR171c in the embryogenic callus and during SE of citrus. Here, we report that miR171 promotes SE and is required for SE in citrus. Overexpression of miR171 restored SE competence in the recalcitrant callus of 'Guoqing No.1' Satsuma mandarin (G1), whereas inhibition of miR171 function by Short Tandem Target Mimic (STTM) weakened SE competence in the strongly embryogenic callus of 'Valencia' sweet orange (V). The comparative transcriptomic analysis in miR171 overexpressed callus line (OE) and the wild type callus (WT) indicated that overexpression of miR171 decreased the expression level of its SCARECROW-LIKE (CsSCL) targets, and activated stress response related biological processes and metabolic processes that are required for cell differentiation. However, CsSCLs were up-regulated in the OE callus during SE induction process, which activated the cell division and developmental processes that are required for embryogenesis progress. Our results validate the function of miR171 in regulation of SE and reveal the biological responses provoked by miR171 in citrus that may promote SE.

Identification of miRNA-mRNA Regulatory Modules Involved in Lipid Metabolism and Seed Development in a Woody Oil Tree (Camellia oleifera)

Tea oil camellia (Camellia oleifera), an important woody oil tree, is a source of seed oil of high nutritional and medicinal value that is widely planted in southern China. However, there is no report on the identification of the miRNAs involved in lipid metabolism and seed development in the high- and low-oil cultivars of tea oil camellia. Thus, we explored the roles of miRNAs in the key periods of oil formation and accumulation in the seeds of tea oil camellia and identified miRNA–mRNA regulatory modules involved in lipid metabolism and seed development. Sixteen small RNA libraries for four development stages of seed oil biosynthesis in high- and low-oil cultivars were constructed. A total of 196 miRNAs, including 156 known miRNAs from 35 families, and 40 novel miRNAs were identified, and 55 significantly differentially expressed miRNAs were found, which included 34 upregulated miRNAs, and 21 downregulated miRNAs. An integrated analysis of the miRNA and mRNA transcriptome sequence data revealed that 10 miRNA–mRNA regulatory modules were related to lipid metabolism; for example, the regulatory modules of ath-miR858b–MYB82/MYB3/MYB44 repressed seed oil biosynthesis, and a regulation module of csi-miR166e-5p–S-ACP-DES6 was involved in the formation and accumulation of oleic acid. A total of 23 miRNA–mRNA regulatory modules were involved in the regulation of the seed size, such as the regulatory module of hpe-miR162a_L-2–ARF19, involved in early seed development. A total of 12 miRNA–mRNA regulatory modules regulating growth and development were identified, such as the regulatory modules of han-miR156a_L+1–SPL4/SBP2, promoting early seed development. The expression changes of six miRNAs and their target genes were validated using quantitative real-time PCR, and the targeting relationship of the cpa-miR393_R-1–AFB2 regulatory module was verified by luciferase assays. These data provide important theoretical values and a scientific basis for the genetic improvement of new cultivars of tea oil camellia in the future.

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.

CRISPR/Cas9 mutants of tomato MICRORNA164 genes uncover their functional specialization in development

Plant MICRORNA164 (miR164) plays diverse regulatory functions by post-transcriptional repression of certain NAM/ATAF/CUC-domain transcription factors. However, the involvement of miR164 in fleshy fruit development and ripening remains poorly understood. Here, de novo prediction of tomato (Solanum lycopersicum) MIR164 genes identified four genes (SlMIR164a-d), of which SlMIR164d has an atypically long pre-miRNA. The roles of the fruit expressed SlMIR164a, b, and d were studied by analysis of their Clustered Regularly Interspaced Short Palindromic Repeats mutants. The slmir164bCR mutant plants exhibited shoot and flower abnormalities characteristic of ectopic boundary specification, whereas the shoot and flower development of slmir164aCR and slmir164dCR mutants were indistinguishable from wild-type. Strikingly, the knockout of SlMIR164a practically eliminated sly-miR164 from the developing and ripening fruit pericarp. The sly-miR164-deficient slmir164aCR fruits were smaller than the wild-type, due to reduced pericarp cell division and expansion, and displayed intense red color and matte, instead of glossy appearance, upon ripening. We found that the fruit skin phenotypes were associated with morphologically abnormal outer epidermis and thicker cuticle. Quantitation of sly-miR164 target transcripts in slmir164aCR ripening fruits demonstrated the upregulation of SlNAM3 and SlNAM2. Specific expression of their miR164-resistant versions in the pericarp resulted in the formation of extremely small fruits with abnormal epidermis, highlighting the importance of their negative regulation by sly-miR164a. Taken together, our results demonstrate that SlMIR164a and SlMIR164b play specialized roles in development: SlMIR164b is required for shoot and flower boundary specification, and SlMIR164a is required for fruit growth including the expansion of its outer epidermis, which determines the properties of the fruit skin.

MicroRNA Profiling During Mulberry (Morus atropurpurea Roxb) Fruit Development and Regulatory Pathway of miR477 for Anthocyanin Accumulation

To understand the mechanism of small non-coding RNAs (miRNA)-mediated development and ripening of mulberry fruits, three small RNA libraries from mulberry fruits at different development stages were constructed, and 159 conserved miRNAs as well as 86 novel miRNAs were successfully identified. Among the miRNAs identified, there were 90 miRNAs which showed differential expression patterns at different stages of fruit development and ripening. The target genes of these differential expressed (DE) miRNAs were involved in growth and development, transcription and regulation of transcription, metabolic processes, and etc. Interestingly, it was found that the expression level of mul-miR477 was increased with fruit ripening, and it can target the antisense lncRNA (Mul-ABCB19AS) of the ATP binding cassette (ABC) transporter B 19 gene (Mul-ABCB19). Our results showed that mul-miR477 can repress the expression of Mul-ABCB19AS and increase the expression of Mul-ABCB19, and it acted as a positive regulator participating anthocyanin accumulation through the regulatory network of mul-miR477-Mul-ABCB19AS-Mul-ABCB19.

MicroRNAs in Woody Plants

MicroRNAs (miRNAs) are small (∼21-nucleotides) non-coding RNAs found in plant and animals. MiRNAs function as critical post-transcriptional regulators of gene expression by binding to complementary sequences in their target mRNAs, leading to mRNA destabilization and translational inhibition. Plant miRNAs have some distinct characteristics compared to their animal counterparts, including greater evolutionary conservation and unique miRNA processing methods. The lifecycle of a plant begins with embryogenesis and progresses through seed germination, vegetative growth, reproductive growth, flowering and fruiting, and finally senescence and death. MiRNAs participate in the transformation of plant growth and development and directly monitor progression of these processes and the expression of certain morphological characteristics by regulating transcription factor genes involved in cell growth and differentiation. In woody plants, a large and rapidly increasing number of miRNAs have been identified, but their biological functions are largely unknown. In this review, we summarize the progress of miRNA research in woody plants to date. In particular, we discuss the potential roles of these miRNAs in growth, development, and biotic and abiotic stresses responses in woody plants.

Small RNAs, Degradome, and Transcriptome Sequencing Provide Insights into Papaya Fruit Ripening Regulated by 1-MCP

As an inhibitor of ethylene receptors, 1-methylcyclopropene (1-MCP) can delay the ripening of papaya. However, improper 1-MCP treatment will cause a rubbery texture in papaya. Understanding of the underlying mechanism is still lacking. In the present work, a comparative sRNA analysis was conducted after different 1-MCP treatments and identified a total of 213 miRNAs, of which 44 were known miRNAs and 169 were novel miRNAs in papaya. Comprehensive functional enrichment analysis indicated that plant hormone signal pathways play an important role in fruit ripening. Through the comparative analysis of sRNAs and transcriptome sequencing, a total of 11 miRNAs and 12 target genes were associated with the ethylene and auxin signaling pathways. A total of 1741 target genes of miRNAs were identified by degradome sequencing, and nine miRNAs and eight miRNAs were differentially expressed under the ethylene and auxin signaling pathways, respectively. The network regulation diagram of miRNAs and target genes during fruit ripening was drawn. The expression of 11 miRNAs and 12 target genes was verified by RT-qPCR. The target gene verification showed that cpa-miR390a and cpa-miR396 target CpARF19-like and CpERF RAP2-12-like, respectively, affecting the ethylene and auxin signaling pathways and, therefore, papaya ripening.

Noncoding RNAs in Medicinal Plants and their Regulatory Roles in Bioactive Compound Production

BACKGROUND

Noncoding RNAs (ncRNAs), such as microRNAs (miRNAs), small interfering RNAs (siRNAs) and long noncoding RNAs (lncRNAs), play significant regulatory roles in plant development and secondary metabolism and are involved in plant response to biotic and abiotic stresses. They have been intensively studied in model systems and crops for approximately two decades and massive amount of information have been obtained. However, for medicinal plants, ncRNAs, particularly their regulatory roles in bioactive compound biosynthesis, are just emerging as a hot research field.

OBJECTIVE

This review aims to summarize current knowledge on herbal ncRNAs and their regulatory roles in bioactive compound production.

RESULTS

So far, scientists have identified thousands of miRNA candidates from over 50 medicinal plant species and 11794 lncRNAs from Salvia miltiorrhiza, Panax ginseng, and Digitalis purpurea. Among them, more than 30 miRNAs and five lncRNAs have been predicted to regulate bioactive compound production.

CONCLUSION

The regulation may achieve through various regulatory modules and pathways, such as the miR397-LAC module, the miR12112-PPO module, the miR156-SPL module, the miR828-MYB module, the miR858-MYB module, and other siRNA and lncRNA regulatory pathways. Further functional analysis of herbal ncRNAs will provide useful information for quality and quantity improvement of medicinal plants.

Integrated sRNAome and RNA-Seq analysis reveals miRNA effects on betalain biosynthesis in pitaya

BACKGROUND

MicroRNAs (miRNAs) and their regulatory functions in anthocyanin, carotenoid, and chlorophyll accumulation have been extensively characterized in many plant species. However, the miRNA regulatory mechanism in betalain biosynthesis remains mostly unknown.

RESULTS

In this study, 126 conserved miRNAs and 41 novel miRNAs were first isolated from Hylocereus monacanthus, among which 95 conserved miRNAs belonged to 53 miRNA families. Thirty-four candidate miRNAs related to betalain biosynthesis were differentially expressed. The expression patterns of those differential expressed miRNAs were analyzed in various pitaya tissues by RT-qPCR. A significantly negative correlation was detected between the expression levels of half those miRNAs and corresponding target genes. Target genes of miRNAs i.e. Hmo-miR157b-HmSPL6-like, Hmo-miR160a-Hpcyt P450-like3, Hmo-miR6020-HmCYP71A8-like, Hmo-novel-2-HmCYP83B1-like, Hmo-novel-15-HmTPST-like, Hmo-miR828a-HmTT2-like, Hmo-miR858-HmMYB12-like, Hmo-miR858-HmMYBC1-like and Hmo-miR858-HmMYB2-like were verified by 5'RACE and transient expression system in tobacco.

CONCLUSIONS

Hmo-miR157b, Hmo-miR160a, Hmo-miR6020 Hmo-novel-2, Hmo-novel-15, Hmo-miR828a and Hmo-miR858 play important roles in pitaya fruit coloration and betalain accumulation. Our findings provide new insights into the roles of miRNAs and their target genes of regulatory functions involved in betalain biosynthesis of pitaya.

Genome-wide identification and characterization of gibberellin metabolic and signal transduction (GA MST) pathway mediating seed and berry development (SBD) in grape (Vitis vinifera L.)

BACKGROUND

Grape is highly sensitive to gibberellin (GA), which is crucial during seed and berry development (SBD) either by itself or by interacting with other hormones, such as auxin, Abscisic acid (ABA), and Cytokinin (CK). However, no systematic analysis of GA metabolic and signal transduction (MST) pathway has been undertaken in grapevine.

RESULTS

In this study, total endogenous GA3 content significantly decreased during SBD, and a total of 48 known genes in GA metabolic (GAM; 31) and signal transduction (ST; 17) pathways were identified in this process. In the GAM pathway, out of 31 genes, VvGA20ox1-1, VvGA3ox4-1, and VvGA2ox1-1 may be the major factors interacting at the green-berry stage (GBS) accompanied with higher accumulation rate. GA biosynthesis was greater than GA inactivation at GBS, confirming the importance of seeds in GA synthesis. The visible correlation between endogenous GA3 content and gene expression profiles suggested that the transcriptional regulation of GA biosynthesis pathway genes was a key mechanism of GA accumulation at the stone-hardening stage (SHS). Interestingly, we observed a negative feedback regulation between VvGA3oxs-VvGAI1-4, VvGA2oxs-VvGAI1-4, and VvGID1B-VvGAI1-4 in maintaining the balance of GA3 content in berries. Moreover, 11 miRNAs may be involved in the modulation of GA MST pathway by mediating their target genes, such as VvGA3ox, VvGID1B, and VvGAMYB. Many genes in auxin, ABA, and CK MST pathways were further identified and found to have a special pattern in the berry, and the crosstalk between GA and these hormones may modulate the complex process during SBD through the interaction gene network of the multihormone pathway. Lastly, based on the expression characterization of multihormone MST pathway genes, a proposed model of the GA-mediated multihormone regulatory network during SBD was proposed.

CONCLUSIONS

Our results provided novel insights into GA-mediated regulatory networks during SBD in grape. The complexity of GA-mediated multihormone ST in SBD was also elucidated, thereby providing valuable information for future functional characterizations of specific genes in grape.

Regulation of DNA Damage Response and Homologous Recombination Repair by microRNA in Human Cells Exposed to Ionizing Radiation

Ionizing radiation may be of both artificial and natural origin and causes cellular damage in living organisms. Radioactive isotopes have been used significantly in cancer therapy for many years. The formation of DNA double-strand breaks (DSBs) is the most dangerous effect of ionizing radiation on the cellular level. After irradiation, cells activate a DNA damage response, the molecular path that determines the fate of the cell. As an important element of this, homologous recombination repair is a crucial pathway for the error-free repair of DNA lesions. All components of DNA damage response are regulated by specific microRNAs. MicroRNAs are single-stranded short noncoding RNAs of 20–25 nt in length. They are directly involved in the regulation of gene expression by repressing translation or by cleaving target mRNA. In the present review, we analyze the biological mechanisms by which miRNAs regulate cell response to ionizing radiation-induced double-stranded breaks with an emphasis on DNA repair by homologous recombination, and its main component, the RAD51 recombinase. On the other hand, we discuss the ability of DNA damage response proteins to launch particular miRNA expression and modulate the course of this process. A full understanding of cell response processes to radiation-induced DNA damage will allow us to develop new and more effective methods of ionizing radiation therapy for cancers, and may help to develop methods for preventing the harmful effects of ionizing radiation on healthy organisms.

Abundant expression of maternal siRNAs is a conserved feature of seed development

Small RNAs are abundant in plant reproductive tissues, especially 24-nucleotide (nt) small interfering RNAs (siRNAs). Most 24-nt siRNAs are dependent on RNA Pol IV and RNA-DEPENDENT RNA POLYMERASE 2 (RDR2) and establish DNA methylation at thousands of genomic loci in a process called RNA-directed DNA methylation (RdDM). In Brassica rapa, RdDM is required in the maternal sporophyte for successful seed development. Here, we demonstrate that a small number of siRNA loci account for over 90% of siRNA expression during B. rapa seed development. These loci exhibit unique characteristics with regard to their copy number and association with genomic features, but they resemble canonical 24-nt siRNA loci in their dependence on RNA Pol IV/RDR2 and role in RdDM. These loci are expressed in ovules before fertilization and in the seed coat, embryo, and endosperm following fertilization. We observed a similar pattern of 24-nt siRNA expression in diverse angiosperms despite rapid sequence evolution at siren loci. In the endosperm, siren siRNAs show a marked maternal bias, and siren expression in maternal sporophytic tissues is required for siren siRNA accumulation. Together, these results demonstrate that seed development occurs under the influence of abundant maternal siRNAs that might be transported to, and function in, filial tissues.

Genome-wide analysis of the citrus B3 superfamily and their association with somatic embryogenesis

BACKGROUND

In citrus, genetic improvement via biotechnology is hindered by the obstacle of in vitro regeneration via somatic embryogenesis (SE). Although a few B3 transcription factors are reported to regulate embryogenesis, little is known about the B3 superfamily in citrus, and which members might be involved in SE.

RESULTS

Genome-wide sequence analysis identified 72 (CsB3) and 69 (CgB3) putative B3 superfamily members in the genomes of sweet orange (Citrus sinensis, polyembryonic) and pummelo (C. grandis, monoembryonic), respectively. Genome duplication analysis indicated that segmental and tandem duplication events contributed to the expansion of the B3 superfamily in citrus, and that the B3 superfamily evolved under the effect of purifying selection. Phylogenetic relationships were well supported by conserved gene structure and motifs outside the B3 domain, which allowed possible functions to be inferred by comparison with homologous genes from Arabidopsis. Expression analysis identified 23 B3 superfamily members that were expressed during SE in citrus and 17 that may play functional roles at late SE stages. Eight B3 genes were identified that were specific to the genome of polyembryonic sweet orange compared to monoembryonic pummelo. Of these eight B3 genes, CsARF19 was found to be specifically expressed at higher levels in embryogenic callus (EC), implying its possible involvement in EC initiation.

CONCLUSIONS

This study provides a genome-wide analysis of the citrus B3 superfamily, including its genome organization, evolutionary features and expression profiles, and identifies specific family members that may be associated with SE.

The Critical Role of miRNAs in Regulation of Flowering Time and Flower Development

Flowering is an important biological process for plants that ensures reproductive success. The onset of flowering needs to be coordinated with an appropriate time of year, which requires tight control of gene expression acting in concert to form a regulatory network. MicroRNAs (miRNAs) are non-coding RNAs known as master modulators of gene expression at the post-transcriptional level. Many different miRNA families are involved in flowering-related processes such as the induction of floral competence, floral patterning, and the development of floral organs. This review highlights the diverse roles of miRNAs in controlling the flowering process and flower development, in combination with potential biotechnological applications for miRNAs implicated in flower regulation.

Small RNA and mRNA Sequencing Reveal the Roles of microRNAs Involved in Pomegranate Female Sterility

Female sterility is a key factor restricting plant reproduction. Our previous studies have revealed that pomegranate female sterility mainly arose from the abnormality of ovule development. MicroRNAs (miRNAs) play important roles in ovule development. However, little is known about the roles of miRNAs in female sterility. In this study, a combined high-throughput sequencing approach was used to investigate the miRNAs and their targeted transcripts involved in female development. A total of 103 conserved and 58 novel miRNAs were identified. Comparative profiling indicated that the expression of 43 known miRNAs and 14 novel miRNAs were differentially expressed between functional male flowers (FMFs) and bisexual flowers (BFs), 30 known miRNAs and nine novel miRNAs showed significant differences among different stages of BFs, and 20 known miRNAs and 18 novel miRNAs exhibited remarkable expression differences among different stages of FMFs. Gene ontology (GO) analyses of 144 predicted targets of differentially expressed miRNAs indicated that the “reproduction process” and “floral whorl development” processes were significantly enriched. The miRNA–mRNA interaction analyses revealed six pairs of candidate miRNAs and their targets associated with female sterility. Interestingly, pg-miR166a-3p was accumulated, whereas its predicted targets (Gglean012177.1 and Gglean013966.1) were repressed in functional male flowers (FMFs), and the interaction between pg-miR166a-3p and its targets (Gglean012177.1 and Gglean013966.1) were confirmed by transient assay. A. thaliana transformed with 35S-pre-pg-miR166a-3p verified the role of pg-miR166a-3p in ovule development, which indicated pg-miR166a-3p’s potential role in pomegranate female sterility. The results provide new insights into molecular mechanisms underlying the female sterility at the miRNA level.

Up-regulation of microRNA targets correlates with symptom severity in Citrus sinensis plants infected with two different isolates of citrus psorosis virus

MAIN CONCLUSION: miRNA targets from Citrus sinensis are predicted and validated using degradome data. They show an up-regulation upon infection with CPsV, with a positive correlation between target expression and symptom severity. Sweet orange (Citrus sinensis) may suffer from disease symptoms induced by virus infections, thus resulting in drastic economic losses. Infection of sweet orange plants with two isolates of citrus psorosis virus (CPsV), expressing different symptomatologies, alters the accumulation of a set of endogenous microRNAs (miRNAs). Here, we predicted ten putative targets from four down-regulated miRNAs: three belonging to the CCAAT-binding transcription factor family (CBFAs); an Ethylene-responsive transcription factor (RAP2-7); an Integrase-type DNA-binding superfamily protein (AP2B); Transport inhibitor response 1 (TIR1); GRR1-like protein 1-related (GRR1); Argonaute 2-related (AGO2), Argonaute 7 (AGO7), and a long non-coding RNA (ncRNA). We validated six of them through analysis of leaf degradome data. Expressions of the validated targets increase in infected samples compared to healthy tissue, showing a more striking up-regulation those samples with higher symptom severity. This study contributes to the understanding of the miRNA-mediated regulation of important transcripts in Citrus sinensis through target validation and shed light in the manner a virus can alter host regulatory mechanisms leading to symptom expression.

Integrated analysis of miRNAs and their targets reveals that miR319c/TCP2 regulates apical bud burst in tea plant (Camellia sinensis)

MAIN CONCLUSION

The roles of microRNA-mediated epigenetic regulation were highlighted in the bud dormancy-activity cycle, implying that certain differentially expressed miRNAs play crucial roles in apical bud burst, such as csn-miR319c/TCP2. microRNAs (miRNAs) are a class of small non-coding RNAs that regulate gene expression by targeting mRNA transcripts for cleavage or directing translational inhibition. To investigate whether miRNAs regulate bud dormancy-activation transition in tea plant, which largely affects the yield and price of tea products and adaptability of tea trees, we constructed small RNA libraries from three different periods of bud dormancy-burst transition. Through sequencing analysis, 262 conserved and 83 novel miRNAs were identified, including 118 differentially expressed miRNAs. Quantitative RT-PCR results for randomly selected miRNAs exhibited that our comprehensive analysis is highly reliable and accurate. The content of caffeine increased continuously from the endodormancy bud to flushing bud, and differentially expressed miRNAs coupling with their targets associated with bud burst were identified. Remarkably, csn-miR319c was downregulated significantly from the quiescent bud to burst bud, while its target gene CsnTCP2 (TEOSINTE BRANCHED/CYCLOIDEA/PROLIFERATING CELL FACTOR 2) displayed opposite expression patterns. Co-transformation experiment in tobacco demonstrated that csn-miR319c can significantly suppress the functions of CsnTCP2. This study on miRNAs and the recognition of target genes could provide new insights into the molecular mechanism of the bud dormancy-activation transition in tea plant.

Identification of miRNAs involved in fruit ripening by deep sequencing of Olea europaea L. transcriptome

BACKGROUND

The ripening process of olive fruits is associated with chemical and/or enzymatic specific transformations, making them particularly attractive to animals and humans. In olive drupes, including 'Cassanese' ones, ripening is usually accompanied by progressive chromatic change, resulting in a final red-brown colourization of both epidermis and mesocarp. This event has an exception in the 'Leucocarpa', in which we observed the destabilization in the equilibrium between the chlorophyll metabolism and that of the other pigments, particularly the anthocyanins, whose switch-off during maturation promotes the white colouration of the fruits. Recently, transcription profiling of 'Leucocarpa' and 'Cassanese' olives along ripening, performed through an Illumina RNA-seq approach, has provided useful insights on genes functions involved in fruit maturation such as those related to the biosynthesis of flavonoids and anthocyanins.

METHODOLOGY

To assess expression alterations of genes involved in flavonoids and anthocyanins biosynthetic pathways during ripening, possibly caused by small nuclear RNA (snRNA) in olive drupes, snRNA libraries from 'Leucocarpa' and 'Cassanese' were constructed with RNAs extracted at 100 and 130 Days After Flowering (DAF) and sequenced by an Illumina approach. 130 conserved microRNAs (miRNA) in the Viridiplantae belonging to 14 miRNA families were identified. Regarding the 130 conserved miRNAs, approximately the 48% were identified in all libraries, 5 and 18 miRNAs were shared between the "Cassanese" (C100, C130) and "Leucocarpa" (L100, L130) libraries, respectively.

CONCLUSION

For the remaining reads not-matching with known miRNAs in the Viridiplantae, we combined secondary structure and minimum free energy to discover novel olive miRNAs. Based on these analyses, 492 sequences were considered as putative novel miRNAs. The putative target genes of identified miRNA were computationally predicted by alignment with the olive drupe transcripts obtained from the same samples. A total of 218 transcripts were predicted as targets of 130 known and 492 putative novel miRNAs. Interestingly, some identified target genes are involved in negative regulation of anthocyanin metabolic process. Quantification of the expression pattern of three miRNA and their target transcripts by qRT-PCR assay confirmed the results of Illumina sequencing.

Comparative Analysis of miRNA Abundance Revealed the Function of Vvi-miR828 in Fruit Coloring in Root Restriction Cultivation Grapevine (Vitis vinifera L.)

Root restriction cultivation leads to early maturation and quality improvement, especially in the anthocyanin content in grapevine. However, the molecular mechanisms that underlie these changes have not been thoroughly elucidated. In this study, four small RNA libraries were constructed, which included the green soft stage (GS) and ripe stage (RS) of ‘Muscat’ (Vitis vinifera L.) grape berries that were grown under root restriction (RR) and in traditional cultivation (no root restriction, CK). A total of 162 known miRNAs and 14 putative novel miRNAs were detected from the four small RNA libraries by high-throughput sequencing. An analysis of differentially expressed miRNAs (DEMs) revealed that 13 miRNAs exhibited significant differences in expression between RR and CK at the GS and RS stages, respectively. For different developmental stages of fruit, 23 and 34 miRNAs showed expression differences between the GS and RS stages in RR and CK, respectively. The expression patterns of the eight DEMs and their targets were verified by qRT-PCR, and the expression profiles of target genes were confirmed to be complementary to the corresponding miRNAs in RR and CK. The function of Vvi-miR828, which showed the down regulated expression in the RS stage under root restriction, was identified by gene transformation in Arabidopsis. The anthocyanin content significantly decreased in transgenic lines, which indicates the regulatory capacity of Vvi-miR828 in fruit coloration. The miRNA expression pattern comparison between RR and CK might provide a means of unraveling the miRNA-mediated molecular process regulating grape berry development under root restricted cultivation.

Transcriptome analysis of strawberry (Fragaria × ananassa) fruits under osmotic stresses and identification of genes related to ascorbic acid pathway

Strawberry (Fragaria ananassa Duch.) is an economically important fruit with a high demand owing to its good taste and medicinal properties. However, its cultivation is affected by various biotic and abiotic stresses. Plants exhibit several intrinsic mechanisms to deal with stresses. In the case of strawberry, the mechanisms highlighting the response against these stresses remain to be elucidated, which has hampered the efforts to develop and cultivate strawberry plants with high yield and quality. Although a virtual reference genome of F. ananassa has recently been published, there is still a lack of information on the expression of genes in response to various stresses. Therefore, to provide molecular information for further studies with strawberry plants, we present the reference transcriptome dataset of F. ananassa, assembled and annotated from deep RNA-Seq data of fruits cultivated under salinity and drought stresses. We also systematically arranged a series of transcripts differentially expressed during these stresses, with an emphasis on genes related to the accumulation of ascorbic acid (AsA). Ascorbic acid is the most potent antioxidant present in these fruits and highly considered during biofortification. A comparison of the expression profile of these genes by RT-qPCR with the content of AsA in the fruits verified a tight regulation and balance between the expression of genes, from biosynthesis, degradation and recycling pathways, resulting in the reduced content of AsA in fruits under these stresses. These results provide a useful repertoire of genes for metabolic engineering, thereby improving the tolerance to stresses.

Evolutionary dynamics of lincRNA transcription in nine citrus species

Long intergenic non-coding RNAs (lincRNAs) play important roles in various biological processes in plants. However, little information is known about the evolutionary characteristics of lincRNAs among closely related plant species. Here, we present a large-scale comparative study of lincRNA transcription patterns in nine citrus species. By strand-specific RNA-sequencing, we identified 18 075 lincRNAs (14 575 lincRNA loci) from 34 tissue samples. The results indicated that the evolution of lincRNA transcription is more rapid than that of mRNAs. In total, 82.8-97.6% of sweet orange (Citrus sinensis) lincRNA genes were shown to have homologous sequences in other citrus genomes. However, only 15.5-28.8% of these genes had transcribed homologous lincRNAs in these citrus species, presenting a strong contrast to the high conservation of mRNA transcription (81.6-84.7%). Moreover, primitive and modern citrus lincRNAs were preferentially expressed in reproductive and vegetative organs, respectively. Evolutionarily conserved lincRNAs showed higher expression levels and lower tissue specificity than species-specific lincRNAs. Notably, we observed a similar tissue expression pattern of homologous lincRNAs in sweet orange and pummelo (Citrus grandis), suggesting that these lincRNAs may be functionally conserved and selectively maintained. We also identified and validated a lincRNA with the highest expression in fruit that acts as an endogenous target mimic (eTM) of csi-miR166c, and two lincRNAs that act as a precursor and target of csi-miR166c, respectively. These lincRNAs together with csi-miR166c could form an eTM166-miR166c-targeted lincRNA regulatory network that possibly affects citrus fruit development.

MicroRNA profiling analysis of developing berries for 'Kyoho' and its early-ripening mutant during berry ripening

BACKGROUND

'Fengzao' is an early-ripening bud mutant of 'Kyoho', which matures nearly 30 days earlier than 'Kyoho'. To gain a better understanding of the regulatory role of miRNAs in early-ripening of grape berry, high-throughput sequencing approach and quantitative RT-PCR validation were employed to identify miRNAs at the genome-wide level and profile the expression patterns of the miRNAs during berry development in 'Kyho' and 'Fengzao', respectively.

RESULTS

Nine independent small RNA libraries were constructed and sequenced in two varieties from key berry development stages. A total of 108 known miRNAs and 61 novel miRNAs were identified. Among that, 159 miRNAs identified in 'Fengzao' all completely expressed in 'Kyoho' and there were 10 miRNAs specifically expressed in 'Kyoho'. The expression profiles of known and novel miRNAs were quite similar between two varieties. As the major differentially expressed miRNAs, novel_144, vvi-miR3626-3p and vvi-miR3626-5p only expressed in 'Kyoho', vvi-miR399b and vvi-miR399e were down-regulated in 'Fengzao', while vvi-miR477b-3p up-regulated in 'Fengzao'. According to the expression analysis and previous reports, miR169-NF-Y subunit, miR398-CSD, miR3626-RNA helicase, miR399- phosphate transporter and miR477-GRAS transcription factor were selected as the candidates for further investigations of miRNA regulation role in the early-ripening of grape. The qRT-PCR analyses validated the contrasting expression patterns for these miRNAs and their target genes.

CONCLUSIONS

The miRNAome of the grape berry development of 'Kyoho', and its early-ripening bud mutant, 'Fengzao' were compared by high-throughput sequencing. The expression pattern of several key miRNAs and their target genes during grape berry development and ripening stages was examined. Our results provide valuable basis towards understanding the regulatory mechanisms of early-ripening of grape berry.

New insights into tomato microRNAs

Cultivated tomato, Solanum lycopersicum, is one of the most common fruits in the global food industry. Together with the wild tomato Solanum pennellii, it is widely used for developing better cultivars. MicroRNAs affect mRNA regulation, inhibiting its translation and/or promoting its degradation. Important proteins involved in these processes are ARGONAUTE and DICER. This study aimed to identify and characterize the genes involved in the miRNA processing pathway, miRNA molecules and target genes in both species. We validated the presence of pathway genes and miRNA in different NGS libraries and 6 miRNA families using quantitative RT-PCR. We identified 71 putative proteins in S. lycopersicum and 108 in S. pennellii likely involved in small RNAs processing. Of these, 29 and 32 participate in miRNA processing pathways, respectively. We identified 343 mature miRNAs, 226 pre-miRNAs in 87 families, including 192 miRNAs, which were not previously identified, belonging to 38 new families in S. lycopersicum. In S. pennellii, we found 388 mature miRNAs and 234 pre-miRNAs contained in 85 families. All miRNAs found in S. pennellii were unpublished, being identified for the first time in our study. Furthermore, we identified 2471 and 3462 different miRNA target in S. lycopersicum and S. pennellii, respectively.

Despacito: the slow evolutionary changes in plant microRNAs

MicroRNAs (miRNAs) are key regulators of gene expression. A handful of miRNAs are broadly conserved in land plants, while the majority are lineage specific; this review describes the processes by which new miRNAs are hypothesized to have emerged. Two major models describe miRNA origins, firstly, de novo emergence via inverted duplication of target gene fragments, and secondly, the expansion and neofunctionalization of existing miRNA families. The occasional acquisition of target sites by previously un-targeted genes adds further dynamism to the process by which miRNAs may shift roles during evolution. Additional factors guiding miRNA evolution include functional constraints on their length and the importance of precursor conservation that is observed in regions above or below the mature miRNA duplex; these regions represent recognition sites for components of biogenesis machinery and direct precursor processing. Insights into the mechanisms of miRNA emergence and divergence are important for understanding plant genome evolution and the impact of miRNA regulatory networks.

Coupling of microRNA-directed phased small interfering RNA generation from long noncoding genes with alternative splicing and alternative polyadenylation in small RNA-mediated gene silencing

MicroRNAs (miRNAs) and phased small interfering RNAs (phasiRNAs) play vital regulatory roles in plant growth and development. Little is known about these small RNAs in litchi (Litchi chinensis), an economically important fruit crop widely cultivated in Southeast Asia. We profiled the litchi small RNA population with various deep-sequencing techniques and in-depth bioinformatic analyses. The genome-wide identification of miRNAs, their target genes, and phasiRNA-generating (PHAS) genes/loci showed that the function of miR482/2118 has expanded, relative to its canonical function. We also discovered that, for 29 PHAS loci, miRNA-mediated phasiRNA production was coupled with alternative splicing (AS) and alternative polyadenylation (APA). Most of these loci encoded long noncoding RNAs. An miR482/2118 targeted locus gave rise to four main transcript isoforms through AS/APA, and diverse phasiRNAs generated from these isoforms appeared to target long terminal repeat (LTR) retrotransposons and other unrelated genes. This coupling enables phasiRNA production from different exons of noncoding PHAS genes and yields diverse phasiRNA populations, both broadening and altering the range of downstream phasiRNA-regulated genes. Our results reveal the diversity of miRNA and phasiRNA in litchi, and demonstrate AS/APA as a new layer of regulation in small RNA-mediated gene silencing.

Genome-wide identification of conserved and novel microRNAs in one bud and two tender leaves of tea plant (Camellia sinensis) by small RNA sequencing, microarray-based hybridization and genome survey scaffold sequences

BACKGROUND

MicroRNAs (miRNAs) are important for plant growth and responses to environmental stresses via post-transcriptional regulation of gene expression. Tea, which is primarily produced from one bud and two tender leaves of the tea plant (Camellia sinensis), is one of the most popular non-alcoholic beverages worldwide owing to its abundance of secondary metabolites. A large number of miRNAs have been identified in various plants, including non-model species. However, due to the lack of reference genome sequences and/or information of tea plant genome survey scaffold sequences, discovery of miRNAs has been limited in C. sinensis.

RESULTS

Using small RNA sequencing, combined with our recently obtained genome survey data, we have identified and analyzed 175 conserved and 83 novel miRNAs mainly in one bud and two tender leaves of the tea plant. Among these, 93 conserved and 18 novel miRNAs were validated using miRNA microarray hybridization. In addition, the expression pattern of 11 conserved and 8 novel miRNAs were validated by stem-loop-qRT-PCR. A total of 716 potential target genes of identified miRNAs were predicted. Further, Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that most of the target genes were primarily involved in stress response and enzymes related to phenylpropanoid biosynthesis. The predicted targets of 4 conserved miRNAs were further validated by 5'RLM-RACE. A negative correlation between expression profiles of 3 out of 4 conserved miRNAs (csn-miR160a-5p, csn-miR164a, csn-miR828 and csn-miR858a) and their targets (ARF17, NAC100, WER and MYB12 transcription factor) were observed.

CONCLUSION

In summary, the present study is one of few such studies on miRNA detection and identification in the tea plant. The predicted target genes of majority of miRNAs encoded enzymes, transcription factors, and functional proteins. The miRNA-target transcription factor gene interactions may provide important clues about the regulatory mechanism of these miRNAs in the tea plant. The data reported in this study will make a huge contribution to knowledge on the potential miRNA regulators of the secondary metabolism pathway and other important biological processes in C. sinensis.

Genome-wide identification of microRNAs responsive to Ectropis oblique feeding in tea plant (Camellia sinensis L.)

The tea plant (Camellia sinensis L.) is vulnerable to the geometrid Ectropis oblique; although microRNAs (miRNAs) are important for plant growth, development and stress response, the function of miRNAs in the response of C. sinensis to stress from E. oblique is unclear. To identify E. oblique stress-responsive miRNAs and their target genes in tea plant, three small RNA libraries were constructed from leaves subjected to mechanical wounding (MW), geometrid attack (GA) and from healthy control (CK) leaves. Using high-throughput sequencing, 130 known miRNAs and 512 novel miRNAs were identified; of these, differential expression under GA stress was observed for 36 known and 139 novel miRNAs. Furthermore, 169 GA-responsive and 173 MW-responsive miRNAs were detected by miRNA microarray. The expression patterns of six GA-responsive miRNAs were validated by qRT-PCR. Several target genes for these miRNAs encode various transcription factors, including ethylene-responsive transcription factors and squamosa promoter-binding-like proteins, which suggests that these miRNAs may regulate stress-responsive transcriptional processes in tea plant. The present findings provide novel insights into miRNA-mediated regulatory mechanisms underlying the response to GA stress, and also offer valuable information for development of pest resistance using RNA interference-based strategies in tea plants.

miR3954 is a trigger of phasiRNAs that affects flowering time in citrus

In plant, a few 22-nt miRNAs direct cleavages of their targets and trigger the biogenesis of phased small interfering RNAs (phasiRNAs) in plant. In this study, we characterized a miRNA triggering phasiRNAs generation, miR3954, and explored its downstream target genes and potential function. Our results demonstrated that miR3954 showed specific expression in the flowers of citrus species, and it targeted a NAC transcription factor (Cs7 g22460) and two non-coding RNA transcripts (lncRNAs, Cs1 g09600 and Cs1 g09635). The production of phasiRNAs was detected from transcripts targeted by miR3954, and was further verified in both sequencing data and transient expression experiments. PhasiRNAs derived from the two lncRNAs targeted not only miR3954-targeted NAC gene but also additional NAC homologous genes. No homologous genes of these two lncRNAs were found in plants other than citrus species, implying that this miR3954-lncRNAs-phasiRNAs-NAC pathway is likely citrus-specific. Transgenic analysis indicated that the miR3954-overexpressing lines showed decreased transcripts of lncRNA, elevated abundance of phasiRNAs and reduced expression of NAC genes. Interestingly, the overexpression of miR3954 leads to early flowering in citrus plants. In summary, our results illustrated a model of the regulatory network of miR3954-lncRNA-phasiRNAs-NAC, which may be functionally involved in flowering in citrus.

Identification of protoplast-isolation responsive microRNAs in Citrus reticulata Blanco by high-throughput sequencing

Protoplast isolation is a stress-inducing process, during which a variety of physiological and molecular alterations take place. Such stress response affects the expression of totipotency of cultured protoplasts. MicroRNAs (miRNAs) play important roles in plant growth, development and stress responses. However, the underlying mechanism of miRNAs involved in the protoplast totipotency remains unclear. In this study, high-throughput sequencing technology was used to sequence two populations of small RNA from calli and callus-derived protoplasts in Citrus reticulata Blanco. A total of 67 known miRNAs from 35 families and 277 novel miRNAs were identified. Among these miRNAs, 18 known miRNAs and 64 novel miRNAs were identified by differentially expressed miRNAs (DEMs) analysis. The expression patterns of the eight DEMs were verified by qRT-PCR. Target prediction showed most targets of the miRNAs were transcription factors. The expression levels of half targets showed a negative correlation to those of the miRNAs. Furthermore, the physiological analysis showed high levels of antioxidant activities in isolated protoplasts. In short, our results indicated that miRNAs may play important roles in protoplast-isolation response.

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.

Identification and characterization of miRNAs in two closely related C4 and C3 species of Cleome by high-throughput sequencing

Cleome gynandra and Cleome hassleriana, which are C4 and C3 plants, respectively, are two species of Cleome. The close genetic relationship between C. gynandra and C. hassleriana provides advantages for discovering the differences in leaf development and physiological processes between C3 and C4 plants. MicroRNAs (miRNAs) are a class of important regulators of various biological processes. In this study, we investigate the differences in the characteristics of miRNAs between C. gynandra and C. hassleriana using high-throughput sequencing technology. In total, 94 and 102 known miRNAs were identified in C. gynandra and C. hassleriana, respectively, of which 3 were specific for C. gynandra and 10 were specific for C. hassleriana. Ninety-one common miRNAs were identified in both species. In addition, 4 novel miRNAs were detected, including three in C. gynandra and three in C. hassleriana. Of these miRNAs, 67 were significantly differentially expressed between these two species and were involved in extensive biological processes, such as glycol-metabolism and photosynthesis. Our study not only provided resources for C. gynandra and C. hassleriana research but also provided useful clues for the understanding of the roles of miRNAs in the alterations of biological processes in leaf tissues during the evolution of the C4 pathway.

Transcriptome Analysis of Leaves, Flowers and Fruits Perisperm of Coffea arabica L. Reveals the Differential Expression of Genes Involved in Raffinose Biosynthesis

Coffea arabica L. is an important crop in several developing countries. Despite its economic importance, minimal transcriptome data are available for fruit tissues, especially during fruit development where several compounds related to coffee quality are produced. To understand the molecular aspects related to coffee fruit and grain development, we report a large-scale transcriptome analysis of leaf, flower and perisperm fruit tissue development. Illumina sequencing yielded 41,881,572 high-quality filtered reads. De novo assembly generated 65,364 unigenes with an average length of 1,264 bp. A total of 24,548 unigenes were annotated as protein coding genes, including 12,560 full-length sequences. In the annotation process, we identified nine candidate genes related to the biosynthesis of raffinose family oligossacarides (RFOs). These sugars confer osmoprotection and are accumulated during initial fruit development. Four genes from this pathway had their transcriptional pattern validated by quantitative reverse transcription polymerase chain reaction (qRT-PCR). Furthermore, we identified ~24,000 putative target sites for microRNAs (miRNAs) and 134 putative transcriptionally active transposable elements (TE) sequences in our dataset. This C. arabica transcriptomic atlas provides an important step for identifying candidate genes related to several coffee metabolic pathways, especially those related to fruit chemical composition and therefore beverage quality. Our results are the starting point for enhancing our knowledge about the coffee genes that are transcribed during the flowering and initial fruit development stages.

Identification of Conserved and Novel MicroRNAs in Blueberry

MicroRNAs (miRNAs) are a class of small endogenous RNAs that play important regulatory roles in cells by negatively affecting gene expression at both transcriptional and post-transcriptional levels. There have been extensive studies aiming to identify miRNAs and to elucidate their functions in various plant species. In the present study, we employed the high-throughput sequencing technology to profile miRNAs in blueberry fruits. A total of 9,992,446 small RNA tags with sizes ranged from 18 to 30 nt were obtained, indicating that blueberry fruits have a large and diverse small RNA population. Bioinformatic analysis identified 412 conserved miRNAs belonging to 29 families, and 35 predicted novel miRNAs that are likely to be unique to blueberries. Among them, expression profiles of five conserved miRNAs were validated by stem loop qRT-PCR. Furthermore, the potential target genes of conserved and novel miRNAs were predicted and subjected to Gene Ontology (GO) annotation. Enrichment analysis of the GO-represented biological processes and molecular functions revealed that these target genes were potentially involved in a wide range of metabolic pathways and developmental processes. Particularly, anthocyanin biosynthesis has been predicted to be directly or indirectly regulated by diverse miRNA families. This study is the first report on genome-wide miRNA profile analysis in blueberry and it provides a useful resource for further elucidation of the functional roles of miRNAs during fruit development and ripening.

Integrated mRNA and miRNA transcriptome reveal a cross-talk between developing response and hormone signaling for the seed kernels of Siberian apricot

Recently, our transcriptomic analysis has identified some functional genes responsible for oil biosynthesis in developing SASK, yet miRNA-mediated regulation for SASK development and oil accumulation is poorly understood. Here, 3 representative periods of 10, 30 and 60 DAF were selected for sRNA sequencing based on the dynamic patterns of growth tendency and oil content of developing SASK. By miRNA transcriptomic analysis, we characterized 296 known and 44 novel miRNAs in developing SASK, among which 36 known and 6 novel miRNAs respond specifically to developing SASK. Importantly, we performed an integrated analysis of mRNA and miRNA transcriptome as well as qRT-PCR detection to identify some key miRNAs and their targets (miR156-SPL, miR160-ARF18, miR164-NAC1, miR171h-SCL6, miR172-AP2, miR395-AUX22B, miR530-P2C37, miR393h-TIR1/AFB2 and psi-miRn5-SnRK2A) potentially involved in developing response and hormone signaling of SASK. Our results provide new insights into the important regulatory function of cross-talk between development response and hormone signaling for SASK oil accumulation.

Genome-scale mRNA and small RNA transcriptomic insights into initiation of citrus apomixis

Nucellar embryony (NE) is an adventitious form of apomixis common in citrus, wherein asexual embryos initiate directly from nucellar cells surrounding the embryo sac. NE enables the fixation of desirable agronomic traits and the production of clonal offspring of virus-free rootstock, but impedes progress in hybrid breeding. In spite of the great importance of NE in citrus breeding and commercial production, little is understood about the underlying molecular mechanisms. In this study, the stages of nucellar embryo initiation (NEI) were determined for two polyembryonic citrus cultivars via histological observation. To explore the genes and regulatory pathways involved in NEI, we performed mRNA-seq and sRNA-seq analyses of ovules immediately prior to and at stages during NEI in the two pairs of cultivars. A total of 305 differentially expressed genes (DEGs) were identified between the poly- and monoembryonic ovules. Gene ontology (GO) analysis revealed that several processes are significantly enriched based on DEGs. In particular, response to stress, and especially response to oxidative stress, was over-represented in polyembryonic ovules. Nearly 150 miRNAs, comprising ~90 conserved and ~60 novel miRNAs, were identified in the ovules of either cultivar pair. Only two differentially expressed miRNAs (DEMs) were identified, of which the novel miRN23-5p was repressed whereas the targets accumulated in the polyembryonic ovules. This integrated study on the transcriptional and post-transcriptional regulatory profiles between poly- and monoembryonic citrus ovules provides new insights into the mechanism of NE, which should contribute to revealing the regulatory mechanisms of plant apomixis.

High-throughput sequencing and degradome analysis reveal altered expression of miRNAs and their targets in a male-sterile cybrid pummelo (Citrus grandis)

Background

G1 + HBP is a male sterile cybrid line with nuclear genome from Hirado Buntan pummelo (C. grandis Osbeck) (HBP) and mitochondrial genome from “Guoqing No.1” (G1, Satsuma mandarin), which provides a good opportunity to study male sterility and nuclear-cytoplasmic cross talk in citrus. High-throughput sRNA and degradome sequencing were applied to identify miRNAs and their targets in G1 + HBP and its fertile type HBP during reproductive development.

Results

A total of 184 known miRNAs, 22 novel miRNAs and 86 target genes were identified. Some of the targets are transcription factors involved in floral development, such as auxin response factors (ARFs), SQUAMOSA promoter binding protein box (SBP-box), MYB, basic region-leucine zipper (bZIP), APETALA2 (AP2) and transport inhibitor response 1 (TIR1). Eight target genes were confirmed to be sliced by corresponding miRNAs using 5’ RACE technology. Based on the sequencing abundance, 42 differentially expressed miRNAs between sterile line G1 + HBP and fertile line HBP were identified. Differential expression of miRNAs and their target genes between two lines was validated by quantitative RT-PCR, and reciprocal expression patterns between some miRNAs and their targets were demonstrated. The regulatory mechanism of miR167a was investigated by yeast one-hybrid and dual-luciferase assays that one dehydrate responsive element binding (DREB) transcription factor binds to miR167a promoter and transcriptionally repress miR167 expression.

Conclusion

Our study reveals the altered expression of miRNAs and their target genes in a male sterile line of pummelo and highlights that miRNA regulatory network may be involved in floral bud development and cytoplasmic male sterility in citrus.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2882-0) contains supplementary material, which is available to authorized users.

Microsatellite markers: what they mean and why they are so useful

Microsatellites or Single Sequence Repeats (SSRs) are extensively employed in plant genetics studies, using both low and high throughput genotyping approaches. Motivated by the importance of these sequences over the last decades this review aims to address some theoretical aspects of SSRs, including definition, characterization and biological function. The methodologies for the development of SSR loci, genotyping and their applications as molecular markers are also reviewed. Finally, two data surveys are presented. The first was conducted using the main database of Web of Science, prospecting for articles published over the period from 2010 to 2015, resulting in approximately 930 records. The second survey was focused on papers that aimed at SSR marker development, published in the American Journal of Botany's Primer Notes and Protocols in Plant Sciences (over 2013 up to 2015), resulting in a total of 87 publications. This scenario confirms the current relevance of SSRs and indicates their continuous utilization in plant science.

Identification of Novel and Conserved miRNAs in Leaves of In vitro Grown Citrus reticulata "Lugan" Plantlets by Solexa Sequencing

MicroRNAs (miRNAs) play essential roles in plant development, but the roles in the in vitro plant development are unknown. Leaves of ponkan plantlets derived from mature embryos at in vitro culture conditions were used to sequence small RNA fraction via Solexa sequencing, and the miRNAs expression was analyzed. The results showed that there were 3,065,625 unique sequences in ponkan, of which 0.79% were miRNAs. The RNA sequences with lengths of 18-25 nt derived from the library were analyzed, leading to the identification of 224 known miRNAs, of which the most abundant were miR157, miR156, and miR166. Three hundred and fifty-eight novel miRNA candidates were also identified, and the number of reads of ponkan novel miRNAs varied from 5 to 168,273. The expression of the most known miRNAs obtained was at low levels, which varied from 5 to 4,946,356. To better understand the role of miRNAs during the preservation of ponkan in vitro plantlet, the expression patterns of cre-miR156a/159b/160a/166a/167a/168a/171/398b were validated by quantitative real-time PCR (qPCR). The results showed that not only the development-associated miRNAs, e.g., cre-miR156/159/166/396, expressed highly at the early preservation period in the in vitro ponkan plantlet leaves but also the stress-related miRNAs, e.g., cre-miR171 and cre-miR398b, expressed highly at the same time. The expression levels of most tested miRNAs were found to decrease after 6 months and the amounts of these miRNAs were kept at low levels at 18 months. After analyzing the expression level of their targets during the reservation of the ponkan in vitro plantlet, development-associated cre-ARF6 and stress-related cre-CSD modules exhibited negative correlation with miR167 and miR398, respectively, indicating an involvement of the miRNAs in the in vitro development of ponkan and function in the conservation of ponkan germplasm.

Comparative Analysis of miRNAs and Their Target Transcripts between a Spontaneous Late-Ripening Sweet Orange Mutant and Its Wild-Type Using Small RNA and Degradome Sequencing

Fruit ripening in citrus is not well-understood at the molecular level. Knowledge of the regulatory mechanism of citrus fruit ripening at the post-transcriptional level in particular is lacking. Here, we comparatively analyzed the miRNAs and their target genes in a spontaneous late-ripening mutant, "Fengwan" sweet orange (MT) (Citrus sinensis L. Osbeck), and its wild-type counterpart ("Fengjie 72-1," WT). Using high-throughput sequencing of small RNAs and RNA degradome tags, we identified 107 known and 21 novel miRNAs, as well as 225 target genes. A total of 24 miRNAs (16 known miRNAs and 8 novel miRNAs) were shown to be differentially expressed between MT and WT. The expression pattern of several key miRNAs and their target genes during citrus fruit development and ripening stages was examined. Csi-miR156k, csi-miR159, and csi-miR166d suppressed specific transcription factors (GAMYBs, SPLs, and ATHBs) that are supposed to be important regulators involved in citrus fruit development and ripening. In the present study, miRNA-mediated silencing of target genes was found under complicated and sensitive regulation in citrus fruit. The identification of miRNAs and their target genes provide new clues for future investigation of mechanisms that regulate citrus fruit ripening.

Identification, characterization and expression analysis of lineage-specific genes within sweet orange (Citrus sinensis)

Background

With the availability of rapidly increasing number of genome and transcriptome sequences, lineage-specific genes (LSGs) can be identified and characterized. Like other conserved functional genes, LSGs play important roles in biological evolution and functions.

Results

Two set of citrus LSGs, 296 citrus-specific genes (CSGs) and 1039 orphan genes specific to sweet orange, were identified by comparative analysis between the sweet orange genome sequences and 41 genomes and 273 transcriptomes. With the two sets of genes, gene structure and gene expression pattern were investigated. On average, both the CSGs and orphan genes have fewer exons, shorter gene length and higher GC content when compared with those evolutionarily conserved genes (ECs). Expression profiling indicated that most of the LSGs expressed in various tissues of sweet orange and some of them exhibited distinct temporal and spatial expression patterns. Particularly, the orphan genes were preferentially expressed in callus, which is an important pluripotent tissue of citrus. Besides, part of the CSGs and orphan genes expressed responsive to abiotic stress, indicating their potential functions during interaction with environment.

Conclusion

This study identified and characterized two sets of LSGs in citrus, dissected their sequence features and expression patterns, and provided valuable clues for future functional analysis of the LSGs in sweet orange.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2211-z) contains supplementary material, which is available to authorized users.

Identification of miRNAs involved in fruit ripening in Cavendish bananas by deep sequencing

BACKGROUND

MicroRNAs (miRNAs) are a family of non-coding small RNAs that play an important regulatory role in various biological processes. Previous studies have reported that miRNAs are closely related to the ripening process in model plants. However, the miRNAs that are closely involved in the banana fruit ripening process remain unknown.

METHODS

Here, we investigated the miRNA populations from banana fruits in response to ethylene or 1-MCP treatment using a deep sequencing approach and bioinformatics analysis combined with quantitative RT-PCR validation.

RESULTS

A total of 125 known miRNAs and 26 novel miRNAs were identified from three libraries. MiRNA profiling of bananas in response to ethylene treatment compared with 1-MCP treatment showed differential expression of 82 miRNAs. Furthermore, the differentially expressed miRNAs were predicted to target a total of 815 target genes. Interestingly, some targets were annotated as transcription factors and other functional proteins closely involved in the development and the ripening process in other plant species. Analysis by qRT-PCR validated the contrasting expression patterns between several miRNAs and their target genes.

CONCLUSIONS

The miRNAome of the banana fruit in response to ethylene or 1-MCP treatment were identified by high-throughput sequencing. A total of 82 differentially expressed miRNAs were found to be closely associated with the ripening process. The miRNA target genes encode transcription factors and other functional proteins, including SPL, APETALA2, EIN3, E3 ubiquitin ligase, β-galactosidase, and β-glucosidase. These findings provide valuable information for further functional research of the miRNAs involved in banana fruit ripening.

Transcriptome-Wide Identification of miRNA Targets under Nitrogen Deficiency in Populus tomentosa Using Degradome Sequencing

miRNAs are endogenous non-coding small RNAs with important regulatory roles in stress responses. Nitrogen (N) is an indispensable macronutrient required for plant growth and development. Previous studies have identified a variety of known and novel miRNAs responsive to low N stress in plants, including Populus. However, miRNAs involved in the cleavage of target genes and the corresponding regulatory networks in response to N stress in Populus remain largely unknown. Consequently, degradome sequencing was employed for global detection and validation of N-responsive miRNAs and their targets. A total of 60 unique miRNAs (39 conserved, 13 non-conserved, and eight novel) were experimentally identified to target 64 mRNA transcripts and 21 precursors. Among them, we further verified the cleavage of 11 N-responsive miRNAs identified previously and provided empirical evidence for the cleavage mode of these miRNAs on their target mRNAs. Furthermore, five miRNA stars (miRNA*s) were shown to have cleavage function. The specificity and diversity of cleavage sites on the targets and miRNA precursors in P. tomentosa were further detected. Identification and annotation of miRNA-mediated cleavage of target genes in Populus can increase our understanding of miRNA-mediated molecular mechanisms of woody plants adapted to low N environments.

Identification of fruit related microRNAs in cucumber (Cucumis sativus L.) using high-throughput sequencing technology

MicroRNAs (miRNAs) are approximately 21 nt noncoding RNAs that influence the phenotypes of different species through the post-transcriptional regulation of gene expression. Although many miRNAs have been identified in a few model plants, less is known about miRNAs specific to cucumber (Cucumis sativus L.). In this study, two libraries of cucumber RNA, one based on fruit samples and another based on mixed samples from leaves, stems, and roots, were prepared for deep-sequencing. A total of 110 sequences were matched to known miRNAs in 47 families, while 56 sequences in 46 families are newly identified in cucumber. Of these, 77 known and 44 new miRNAs were differentially expressed, with a fold-change of at least 2 and p-value < 0.05. In addition, we predicted the potential targets of known and new miRNAs. The identification and characterization of known and new miRNAs will enable us to better understand the role of these miRNAs in the formation of cucumber fruit.

Genomewide analysis of small RNAs in nonembryogenic and embryogenic tissues of citrus: microRNA- and siRNA-mediated transcript cleavage involved in somatic embryogenesis

Somatic embryogenesis (SE) is a process of somatic cells becoming dedifferentiated and generating embryos. SE has been widely used in biotechnology as a powerful way of regeneration and a model system for studying plant embryogenesis, but the controlling mechanisms of SE are far from clear. Here, we show the genomewide profiles of miRNAs/siRNAs and their target genes in nonembryogenic and embryogenic tissues of 'Valencia' sweet orange. By high-throughput sequencing (HTS) of small RNAs and RNA degradome tags, we identified 50 known and 45 novel miRNAs, 130 miniature inverted-repeat transposable elements (MITEs) derived, 94 other and 235 phased small interfering RNAs (siRNAs), as well as 203 target genes. The majority of the abundantly expressed miRNAs/siRNAs exhibit lower expression levels in embryogenic callus (EC) or during SE process than in nonembryogenic callus (NEC), which is supposed to derepress the target genes that are involved in development and stress response, thus to activate the biological processes required for cell differentiation. However, the conserved csi-miR156a/b, miR164b and 171c directed suppression of specific transcription factors (TFs) are supposed to inactivate the postembryonic growth thus to maintain normal SE. In this study, miRNA- and siRNA-mediated silencing of target genes was found under sophisticated regulation in citrus SE system; the enhancement effect of specific conserved miRNAs on SE was discussed, providing new clues for future investigation of mechanisms that control SE.

miRNAome analysis associated with anatomic and transcriptomic investigations reveal the polar exhibition of corky split vein in boron deficient Citrus sinensis

Corky split vein can develop under long-term boron deficient conditions in Citrus sinensis L. Osbeck cv. Newhall. This symptom only occurs in the upper rather than the lower epidermis of old leaves. Our previous study demonstrated that vascular hypertrophy was involved in the symptoms, and the 3rd developmental stage of corky split vein (BD3) was the critical stage for phenotype formation. Here, we performed an intensive study on the BD3 vein and its control sample (CK3 vein). A lignin test demonstrated that the lignin content in BD3 vein was approximately 1.7 times more than the CK3 vein. Anatomical investigation of the corky split vein indicated that the upper epidermis was destroyed by overgrowing vascular cells, and the increased lignin may contribute to vascular cell differentiation and wounding-induced lignification. In a subsequent small RNA sequencing of the BD3 and CK3 veins, 99 known miRNAs and 22 novel miRNAs were identified. Comparative profiling of these miRNAs demonstrated that the 57 known miRNAs and all novel miRNAs exhibited significant expression differences between the two small RNAs libraries of the BD3 and CK3 veins. Associated with our corresponding digital gene expression data, we propose that the decreased expression of two miRNAs, csi-miR156b and csi-miR164, which leads to the up-regulation of their target genes, SPLs (csi-miR156b-targeted) and CUC2 (csi-miR164-targeted), may promote vascular cell division and orderless stage transition in old leaves.