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

Advances in the Study of the Transcriptional Regulation Mechanism of Plant miRNAs

MicroRNAs (miRNA) are a class of endogenous, non-coding, small RNAs with about 22 nucleotides (nt), that are widespread in plants and are involved in various biological processes, such as development, flowering phase transition, hormone signal transduction, and stress response. The transcriptional regulation of miRNAs is an important process of miRNA gene regulation, and it is essential for miRNA biosynthesis and function. Like mRNAs, miRNAs are transcribed by RNA polymerase II, and these transcription processes are regulated by various transcription factors and other proteins. Consequently, the upstream genes regulating miRNA transcription, their specific expression, and the regulating mechanism were reviewed to provide more information for further research on the miRNA regulatory mechanism and help to further understand the regulatory networks of plant miRNAs.

Evolutionary assessment of SQUAMOSA PROMOTER BINDING PROTEIN-LIKE genes in citrus relatives with a specific focus on flowering

Phase transition and floral induction in citrus requires several years of juvenility after germination. Such a long period of juvenility has been a major hindrance to its genetic improvement program. Studies have shown that miR156 along with its downstream genes SQUAMOSA PROMOTER BINDING PROTEINS (SBP) and SBP-LIKE (SPL) mediate the phase transition and floral induction process in plants. Our current study has systematically analyzed SPLs in 15 different citrus-related species, systematically annotated them based on their close homology to their respective Arabidopsis orthologs, and confirmed the functional attributes of the selected members in floral precocity. The majority of the species harbored 15 SPLs. Their cis-element assessment suggested the involvement of the SPLs in diverse developmental and physiological processes in response to different biotic and abiotic cues. Among all, SPL5, SPL9, and SPL11 stood out as consistently differentially expressed SPLs in the adult and young tissues of different citrus-related species. Independent overexpression of their F. hindsii orthologs (FhSPL5, FhSPL9, and FhSPL11) brought an enhanced expression of endogenous FLOWERING LOCUS T leading to the significantly precocious flowering in transgenic Arabidopsis lines. Future study of the genes in the citrus plant itself is expected to conclude the assessments made in the current study.

Spatiotemporal profiles of gene activity in stamen delineate nucleo-cytoplasmic interaction in a male-sterile somatic cybrid citrus

Cytoplasmic male sterility (CMS) has long been used to produce seedless fruits in perennial woody crops like citrus. A male-sterile somatic cybrid citrus (G1 + HBP) was generated by protoplast fusion between a CMS callus parent 'Guoqing No. 1' Satsuma mandarin (Citrus unshiu, G1) and a fertile mesophyll parent Hirado Buntan pummelo (Citrus grandis, HBP). To uncover the male-sterile mechanism of G1 + HBP, we compared the transcriptome profiles of stamen organ and cell types at five stages between G1 + HBP and HBP, including the initial stamen primordia, enlarged stamen primordia, pollen mother cells, tetrads, and microspores captured by laser microdissection. The stamen organ and cell types showed distinct gene expression profiles. A majority of genes involved in stamen development were differentially expressed, especially CgAP3.2, which was downregulated in enlarged stamen primordia and upregulated in tetrads of G1 + HBP compared with HBP. Jasmonic acid- and auxin-related biological processes were enriched among the differentially expressed genes of stamen primordia, and the content of jasmonic acid biosynthesis metabolites was higher in flower buds and anthers of G1 + HBP. In contrast, the content of auxin biosynthesis metabolites was lower in G1 + HBP. The mitochondrial tricarboxylic acid cycle and oxidative phosphorylation processes were enriched among the differentially expressed genes in stamen primordia, meiocytes, and microspores, indicating the dysfunction of mitochondria in stamen organ and cell types of G1 + HBP. Taken together, the results indicate that malfunction of mitochondria-nuclear interaction might cause disorder in stamen development, and thus lead to male sterility in the citrus cybrid.

Protoplast Technology and Somatic Hybridisation in the Family Apiaceae

Species of the family Apiaceae occupy a major market share but are hitherto dependent on open pollinated cultivars. This results in a lack of production uniformity and reduced quality that has fostered hybrid seed production. The difficulty in flower emasculation led breeders to use biotechnology approaches including somatic hybridization. We discuss the use of protoplast technology for the development of somatic hybrids, cybrids and in-vitro breeding of commercial traits such as CMS (cytoplasmic male sterility), GMS (genetic male sterility) and EGMS (environment-sensitive genic male sterility). The molecular mechanism(s) underlying CMS and its candidate genes are also discussed. Cybridization strategies based on enucleation (Gamma rays, X-rays and UV rays) and metabolically arresting protoplasts with chemicals such as iodoacetamide or iodoacetate are reviewed. Differential fluorescence staining of fused protoplast as routinely used can be replaced by new tagging approaches using non-toxic proteins. Here, we focused on the initial plant materials and tissue sources for protoplast isolation, the various digestion enzyme mixtures tested, and on the understanding of cell wall re-generation, all of which intervene in somatic hybrids regeneration. Although there are no alternatives to somatic hybridization, various approaches also discussed are emerging, viz., robotic platforms, artificial intelligence, in recent breeding programs for trait identification and selection.

MicroRNAs Involved in Regulatory Cytoplasmic Male Sterility by Analysis RNA-seq and Small RNA-seq in Soybean

Cytoplasmic male sterility (CMS) is an important plant characteristic for exploiting heterosis to enhance crop traits during breeding. However, the CMS regulatory network remains unclear in plants, even though researchers have attempted to isolate genes associated with CMS. In this study, we performed high-throughput sequencing and degradome analyses to identify microRNAs (miRNAs) and their targets in a soybean CMS line (JLCMS9A) and its maintainer line (JLCMS9B). Additionally, the differentially expressed genes during reproductive development were identified using RNA-seq data. A total of 280 miRNAs matched soybean miRNA sequences in miRBase, including mature miRNAs and pre-miRNAs. Of the 280 miRNAs, 30, 23, and 21 belonged to the miR166, miR156, and miR171 families, respectively. Moreover, 410 novel low-abundant miRNAs were identified in the JLCMS9A and JLCMS9B flower buds. Furthermore, 303 and 462 target genes unique to JLCMS9A and JLCMS9B, respectively, as well as 782 common targets were predicted based on the degradome analysis. Target genes differentially expressed between the CMS line and the maintainer line were revealed by an RNA-seq analysis. Moreover, all target genes were annotated with diverse functions related to biological processes, cellular components, and molecular functions, including transcriptional regulation, the nucleus, meristem maintenance, meristem initiation, cell differentiation, auxin-activated signaling, plant ovule development, and anther development. Finally, a network was built based on the interactions. Analyses of the miRNA, degradome, and transcriptome datasets generated in this study provided a comprehensive overview of the reproductive development of a CMS soybean line. The data presented herein represent useful information for soybean hybrid breeding. Furthermore, the study results indicate that miRNAs might contribute to the soybean CMS regulatory network by modulating the expression of CMS-related genes. These findings lay the foundation for future studies on the molecular mechanisms underlying soybean CMS.

Transcriptome and MiRNAomics Analyses Identify Genes Associated with Cytoplasmic Male Sterility in Cotton (Gossypium hirsutum L.)

Cytoplasmic male sterility (CMS) is important for large-scale hybrid seed production. Rearrangements in the mitochondrial DNA (mtDNA) for the cotton (Gossypium hirsutum L.) CMS line J4A were responsible for pollen abortion. However, the expression patterns of nuclear genes associated with pollen abortion and the molecular basis of CMS for J4A are unknown, and were the objectives of this study by comparing J4A with the J4B maintainer line. Cytological evaluation of J4A anthers showed that microspore abortion occurs during meiosis preventing pollen development. Changes in enzyme activity of mitochondrial respiratory chain complex IV and mitochondrial respiratory chain complex V and the content of ribosomal protein and ATP during anther abortion were observed for J4A suggesting insufficient synthesis of ATP hindered pollen production. Additionally, levels of sucrose, starch, soluble sugar, and fructose were significantly altered in J4A during the meiosis stage, suggesting reduced sugar metabolism contributed to sterility. Transcriptome and miRNAomics analyses identified 4461 differentially expressed mRNAs (DEGs) and 26 differentially expressed microRNAs (DEMIs). Pathway enrichment analysis indicated that the DEMIs were associated with starch and sugar metabolism. Six deduced target gene regulatory pairs that may participate in CMS were identified, ghi-MIR7484-10/mitogen-activated protein kinase kinase 6 (MAPKK6), ghi-undef-156/agamous-like MADS-box protein AGL19 (AGL19), ghi-MIR171-1-22/SNF1-related protein kinase regulatory subunit gamma-1 and protein trichome birefringence-like 38, and ghi-MIR156-(8/36)/WRKY transcription factor 28 (WRKY28). Overall, a putative CMS mechanism involving mitochondrial dysfunction, the ghi-MIR7484-10/MAPKK6 network, and reduced glucose metabolism was suggested, and ghi-MIR7484-10/MAPKK6 may be related to abnormal microspore meiosis and induction of excessive sucrose accumulation in anthers.

Expression characteristics of pineal miRNAs at ovine different reproductive stages and the identification of miRNAs targeting the AANAT gene

Background

Many recent studies have shown that miRNAs play important roles in the regulation of animal reproduction, including seasonal reproduction. The pineal gland is a crucial hub in the regulation of seasonal reproduction. However, little is known about the expression characteristics of pineal miRNAs in different reproductive seasons (anestrus and breeding season). Therefore, the expression profiles and regulatory roles of ovine pineal miRNAs were investigated during different reproductive stages using Solexa sequencing technology and dual luciferase reporter assays.

Results

A total of 427 miRNAs were identified in the sheep pineal gland. Significant differences in miRNA expression were demonstrated between anestrus and the breeding season in terms of the frequency distributions of miRNA lengths, number of expressed miRNAs, and specifically and highly expressed miRNAs in each reproductive stage. KEGG analysis of the differentially expressed (DE) miRNAs between anestrus and the breeding season indicated that they are significantly enriched in pathways related to protein synthesis, secretion and uptake. Furthermore, transcriptome analysis revealed that many target genes of DE miRNAs in the ribosome pathway showed relatively low expression in the breeding season. On the other hand, analyses combining miRNA-gene expression data with target relationship validation in vitro implied that miR-89 may participate in the negative regulation of aralkylamine N-acetyltransferase (AANAT) mRNA expression by targeting its 3’UTR at a unique binding site.

Conclusions

Our results provide new insights into the expression characteristics of sheep pineal miRNAs at different reproductive stages and into the negative regulatory effects of pineal miRNAs on AANAT mRNA expression.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07536-y.

Identification of microRNAs and their gene targets in cytoplasmic male sterile and fertile maintainer lines of pigeonpea

Comparative analysis of genome-wide miRNAs and their gene targets between cytoplasmic male sterile (CMS) and fertile lines of pigeonpea suggests a possible role of miRNA-regulated pathways in reproductive development. Exploitation of hybrid vigor using CMS technology has delivered nearly 50% yield gain in pigeonpea. Among various sterility-inducing cytoplasms (A₁-A₉) reported so far in pigeonpea, A₂ and A₄ are the two major sources that facilitate hybrid seed production. Recent evidence suggests involvement of micro RNA in vast array of biological processes including plant reproductive development. In pigeonpea, information about the miRNAs is insufficient. In view of this, we sequenced six small RNA libraries of CMS line UPAS 120A and isogenic fertile line UPAS 120B using Illumina technology. Results revealed 316 miRNAs including 248 known and 68 novel types. A total of 637 gene targets were predicted for known miRNAs, while 324 genes were associated with novel miRNAs. Degradome analysis revealed 77 gene targets of predicted miRNAs, which included a variety of transcription factors playing key roles in plant reproduction such as F-box family proteins, apetala 2, auxin response factors, ethylene-responsive factors, homeodomain-leucine zipper proteins etc. Differential expression of both known and novel miRNAs implied roles for both conserved as well as species-specific players. We also obtained several miRNA families such as miR156, miR159, miR167 that are known to influence crucial aspects of plant fertility. Gene ontology and pathway level analyses of the target genes showed their possible implications for crucial events during male reproductive development such as tapetal degeneration, pollen wall formation, retrograde signaling etc. To the best of our knowledge, present study is first to combine deep sequencing of small RNA and degradome for elucidating the role of miRNAs in flower and male reproductive development in pigeonpea.

Integrated transcriptome, small RNA and degradome sequencing approaches proffer insights into chlorogenic acid biosynthesis in leafy sweet potato

Leafy sweet potato is rich in total phenolics (TP) which play key roles in health protection, the chlorogenic acid (CGA) constitutes the major components of phenolic compounds in leafy sweet potato. Unfortunately, the mechanism of CGA biosynthesis in leafy sweet potato is unclear. To dissect the mechanisms of CGA biosynthesis, we performed transcriptome, small RNA (sRNA) and degradome sequencing of one low-CGA content and one high-CGA content genotype at two stages. A total of 2,333 common differentially expressed genes (DEGs) were identified, and the enriched DEGs were related to photosynthesis, starch and sucrose metabolism and phenylpropanoid biosynthesis. The functional genes, such as CCR, CCoAOMT and HCT in the CGA biosynthetic pathway were down-regulated, indicating that the way to lignin was altered, and two possible CGA biosynthetic routes were hypothesized. A total of 38 DE miRNAs were identified, and 1,799 targets were predicated for 38 DE miRNAs by using in silico approaches. The target genes were enriched in lignin and phenylpropanoid catabolic processes. Transcription factors (TFs) such as apetala2/ethylene response factor (AP2/ERF) and Squamosa promoter binding protein-like (SPL) predicated in silico were validated by degradome sequencing. Association analysis of the DE miRNAs and transcriptome datasets identified that miR156 family negatively targeted AP2/ERF and SPL. Six mRNAs and six miRNAs were validated by qRT-PCR, and the results showed that the expression levels of the mRNAs and miRNAs were consistent with the sequencing data. This study established comprehensive functional genomic resources for the CGA biosynthesis, and provided insights into the molecular mechanisms involving in this process. The results also enabled the first perceptions of the regulatory roles of mRNAs and miRNAs, and offered candidate genes for leafy sweet potato improvements.

Knock-down of δ-aminolevulinic acid dehydratase via virus-induced gene silencing alters the microRNA biogenesis and causes stress-related reactions in citrus plants

The δ-aminolevulinic acid (δ-ALA) is an intermediate in the biosynthetic pathway of tetrapyrroles. Tetrapyrroles play vital roles in many biological processes such as photosynthesis, respiration, and light-sensing. ALA-dehydratase (ALAD) combines two molecules of δ-ALA to form porphobilinogen. In citrus, the silencing of ALAD caused discrete yellow spots and necrosis in leaves and stems. Additionally, it caused rapid death in developing new shoots. Herein, we hypothesize that the accumulation of δ-ALA results in severe stress and reduced meristem development. For that reason, we investigated the dynamic changes in the expression profiles of 23 microRNA (miRNA) identified through small RNA sequencing, from CTV-tALAD plants in comparison with healthy C. macrophylla and C. macrophylla infiltrated with CTV-wt. Furthermore, we reported the effect of ALAD silencing on the total phenolics, H₂O_2, and reactive oxygen species (ROS) levels, to examine the possibilities of miRNAs involving the regulation of these pathways. Our results showed that the total phenolics content, H₂O₂, and O₂^- levels were increased in CTV-tALAD plants. Moreover, 63 conserved miRNA members belonging to 23 different miRNA families were differentially expressed in CTV-tALAD plants compared to controls. The identified miRNAs are implicated in auxin biosynthesis and signaling, axillary shoot meristem formation and leaf morphology, starch metabolism, and oxidative stress. Collectively, our findings suggested that ALAD silencing initiates stress on citrus plants. As a result, CTV-tALAD plants exhibit reduced metabolic rate, growth, and development in order to cope with the stress that resulted from the accumulation of δ-ALA. This cascade of events led to leaf, stem, and meristem necrosis and failure of new shoot development.

The miR399-CsUBC24 Module Regulates Reproductive Development and Male Fertility in Citrus

MicroRNA399 (miR399) regulates phosphate homeostasis in plants by down-regulating the expression of PHOSPHATE2 (PHO2, or UBC24 encoding the ubiquitin-conjugating E2 enzyme). We previously identified CsmiR399a.1 in a small RNA sequencing screen of a male-sterile somatic cytoplasmic hybrid (or cybrid) of pummelo (Citrus grandis). Here, we report that miR399 affects reproductive development and male fertility in citrus. Down-regulation of CsmiR399a.1 using a short tandem target mimic (STTM) led to abnormal floral development, inhibition of anther dehiscence, and decreased pollen fertility. When grown in inorganic phosphate (Pi)-sufficient conditions, CsmiR399a.1-STTM plants had lower total phosphorus content in their leaves than the wild type and showed typical symptoms of Pi deficiency. In CsmiR399a.1-STTM plants, the expression of genes involved in starch metabolism and Pi homeostasis was significantly different than in the wild type. Thus, we conclude that miR399-STTM mimicked Pi deficiency, disturbed starch metabolism, and was responsible for pollen grain collapse in the transgenic lines. We identified CsUBC24, a citrus homolog of Arabidopsis (Arabidopsis thaliana) AtUBC24 (PHO2), as a target of CsmiR399a.1 that physically interacts with the floral development regulators SEPALLATA family (CsSEP1.1, CsSEP1.2, and CsSEP3) and the anther dehiscence regulator INDUCER OF CBF EXPRESSION1 (CsICE1). We hypothesize that CsUBC24 downregulates the CsSEPs, which disrupts the floral meristem identity regulatory network and leads to developmental abnormalities in flowers. By interacting with CsICE1, CsUBC24 disturbs stomate function on the anther surface, which inhibits anther dehiscence. These findings indicate that a miR399-based mechanism influences both reproductive development and male fertility in citrus.

Identification of microRNAs and their targets in inflorescences of an Ogura-type cytoplasmic male-sterile line and its maintainer fertile line of turnip (Brassica rapa ssp. rapifera) via high-throughput sequencing and degradome analysis

Cytoplasmic male sterility (CMS) is a widely used trait in angiosperms caused by perturbations in nucleus-mitochondrion interactions that suppress the production of functional pollen. MicroRNAs (miRNAs) are small non-coding RNAs that act as regulatory molecules of transcriptional or post-transcriptional gene silencing in plants. The discovery of miRNAs and their possible implications in CMS induction provides clues for the intricacies and complexity of this phenomenon. Previously, we characterized an Ogura-CMS line of turnip (Brassica rapa ssp. rapifera) that displays distinct impaired anther development with defective microspore production and premature tapetum degeneration. In the present study, high-throughput sequencing was employed for a genome-wide investigation of miRNAs. Six small RNA libraries of inflorescences collected from the Ogura-CMS line and its maintainer fertile (MF) line of turnip were constructed. A total of 120 pre-miRNAs corresponding to 89 mature miRNAs were identified, including 87 conversed miRNAs and 33 novel miRNAs. Among these miRNAs, the expression of 10 differentially expressed mature miRNAs originating from 12 pre-miRNAs was shown to have changed by more than two-fold between inflorescences of the Ogura-CMS line and inflorescences of the MF line, including 8 down- and 2 up-regulated miRNAs. The expression profiles of the differentially expressed miRNAs were confirmed by stem-loop quantitative real-time PCR. In addition, to identify the targets of the identified miRNAs, a degradome analysis was performed. A total of 22 targets of 25 miRNAs and 17 targets of 28 miRNAs were identified as being involved in the reproductive development for Ogura-CMS and MF lines of turnip, respectively. Negative correlations of expression patterns between partial miRNAs and their targets were detected. Some of these identified targets, such as squamosa promoter-binding-like transcription factor family proteins, auxin response factors and pentatricopeptide repeat-containing proteins, were previously reported to be involved in reproductive development in plants. Taken together, our results can help improve the understanding of miRNA-mediated regulatory pathways that might be involved in CMS occurrence in turnip.

Integrated Analysis of Small RNA, Transcriptome and Degradome Sequencing Provides New Insights into Floral Development and Abscission in Yellow Lupine (Lupinus luteus L.)

The floral development in an important legume crop yellow lupine (Lupinus luteus L., Taper cv.) is often affected by the abscission of flowers leading to significant economic losses. Small non-coding RNAs (sncRNAs), which have a proven effect on almost all developmental processes in other plants, might be of key players in a complex net of molecular interactions regulating flower development and abscission. This study represents the first comprehensive sncRNA identification and analysis of small RNA, transcriptome and degradome sequencing data in lupine flowers to elucidate their role in the regulation of lupine generative development. As shedding in lupine primarily concerns flowers formed at the upper part of the inflorescence, we analyzed samples from extreme parts of raceme separately and conducted an additional analysis of pedicels from abscising and non-abscising flowers where abscission zone forms. A total of 394 known and 28 novel miRNAs and 316 phased siRNAs were identified. In flowers at different stages of development 59 miRNAs displayed differential expression (DE) and 46 DE miRNAs were found while comparing the upper and lower flowers. Identified tasiR-ARFs were DE in developing flowers and were strongly expressed in flower pedicels. The DEmiR-targeted genes were preferentially enriched in the functional categories related to carbohydrate metabolism and plant hormone transduction pathways. This study not only contributes to the current understanding of how lupine flowers develop or undergo abscission but also holds potential for research aimed at crop improvement.

Small RNA profiling from meiotic and post-meiotic anthers reveals prospective miRNA-target modules for engineering male fertility in sorghum

Understanding male gametophyte development is essential to augment hybrid production in sorghum. Although small RNAs are known to critically influence anther/pollen development, their roles in sorghum reproduction have not been deciphered yet. Here, we report small RNA profiling and high-confidence annotation of microRNAs (miRNAs) from meiotic and post-meiotic anthers in sorghum. We identified 262 miRNAs (82 known and 180 novel), out of which 58 (35 known and 23 novel) exhibited differential expression between two stages. Out of 35 differentially expressed known miRNAs, 13 are known to regulate anther/pollen development in other plant species. We also demonstrated conserved spatiotemporal patterns of 21- and 24-nt phasiRNAs and their respective triggers, miR2118 and miR2275, in sorghum anthers as evidenced in other monocots. miRNA target identification yielded 5622 modules, of which 46 modules comprising 16 known and 8 novel miRNA families with 38 target genes are prospective candidates for engineering male fertility in grasses.

Uncovering anthocyanin biosynthesis related microRNAs and their target genes by small RNA and degradome sequencing in tuberous roots of sweetpotato

BACKGROUND

Compared with white-fleshed sweetpotato (WFSP), purple-fleshed sweetpotato (PFSP) is a desirable resource for functional food development because of the abundant anthocyanin accumulation in its tuberous roots. Some studies have shown that the expression regulation mediated by miRNA plays an important role in anthocyanin biosynthesis in plants. However, few miRNAs and their corresponding functions related to anthocyanin biosynthesis in tuberous roots of sweetpotato have been known.

RESULTS

In this study, small RNA (sRNA) and degradome libraries from the tuberous roots of WFSP (Xushu-18) and PFSP (Xuzishu-3) were constructed, respectively. Totally, 191 known and 33 novel miRNAs were identified by sRNA sequencing, and 180 target genes cleaved by 115 known ib-miRNAs and 5 novel ib-miRNAs were identified by degradome sequencing. Of these, 121 miRNAs were differently expressed between Xushu-18 and Xuzishu-3. Integrated analysis of sRNA, degradome sequencing, GO, KEGG and qRT-PCR revealed that 26 differentially expressed miRNAs and 36 corresponding targets were potentially involved in the anthocyanin biosynthesis. Of which, an inverse correlation between the expression of ib-miR156 and its target ibSPL in WFSP and PFSP was revealed by both qRT-PCR and sRNA sequencing. Subsequently, ib-miR156 was over-expressed in Arabidopsis. Interestingly, the ib-miR156 over-expressing plants showed suppressed abundance of SPL and a purplish phenotype. Concomitantly, upregulated expression of four anthocyanin pathway genes was detected in transgenic Arabidopsis plants. Finally, a putative ib-miRNA-target model involved in anthocyanin biosynthesis in sweetpotato was proposed.

CONCLUSIONS

The results represented a comprehensive expression profiling of miRNAs related to anthocyanin accumulation in sweetpotato and provided important clues for understanding the regulatory network of anthocyanin biosynthesis mediated by miRNA in tuberous crops.

Transcriptome analysis of miRNAs expression reveals novel insights into adventitious root formation in lotus (Nelumbo nucifera Gaertn.)

MicroRNA (miRNA)-regulated gene expression plays an important role in various plant metabolic processes. Although adventitious roots are critical to plant growth in lotus, the role of miRNA in AR formation remains unclear. Expression profiling of miRNAs was carried out during three different developmental stages of ARs in lotus: no induction of AR stage, initial stage of ARs, and maximum number of ARs. These data are referenced with the whole lotus genome as already identified through high-throughput tag-sequencing. 1.3 × 10⁷ tags were achieved, of which 11,035,798, 11,436,062, and 12,542,392 clean tags were obtained from each stage, respectively. miRNA analysis revealed that miRNAs were less than 10% among all small RNAs. In total, 310 miRNAs (240 up-regulated and 70 down-regulated miRNAs) exhibited expression changes from the no induction stage to the initial stage. Moreover, expression of 140 miRNAs was increased and that of 123 miRNAs was decreased between the initial and maximum AR stages, mostly by ~ - 4-4-fold. miRNAs involved in metabolic pathways differed between the initial stage/no induction stage and the maximum number stage/initial stage. Several miRNAs in the initial stage/no induction stage were related to plant hormone metabolism and pyruvate and MAPK pathways, while major miRNAs in the maximum number stage/initial stage were involved in carbohydrate metabolism. All differentially expressed miRNAs associated with AR formation from the initial stage to maximum stage were also analyzed. The expression of 16 miRNAs was determined using qRT-PCR. This work provides a general insight into miRNA regulation during AR formation in lotus.

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.

An integrated proteomic and metabolomic study to evaluate the effect of nucleus-cytoplasm interaction in a diploid citrus cybrid between sweet orange and lemon

Our results provide a comprehensive overview how the alloplasmic condition might lead to a significant improvement in citrus plant breeding, developing varieties more adaptable to a wide range of conditions. Citrus cybrids resulting from somatic hybridization hold great potential in plant improvement. They represent effective products resulting from the transfer of organelle-encoded traits into cultivated varieties. In these cases, the plant coordinated array of physiological, biochemical, and molecular functions remains the result of integration among different signals, which derive from the compartmentalized genomes of nucleus, plastids and mitochondria. To dissect the effects of genome rearrangement into cybrids, a multidisciplinary study was conducted on a diploid cybrid (C2N), resulting from a breeding program aimed to improve interesting agronomical traits for lemon, the parental cultivars 'Valencia' sweet orange (V) and 'femminello' lemon (F), and the corresponding somatic allotetraploid hybrid (V + F). In particular, a differential proteomic analysis, based on 2D-DIGE and MS procedures, was carried out on leaf proteomes of C2N, V, F and V + F, using the C2N proteome as pivotal condition. This investigation revealed differentially represented protein patterns that can be associated with genome rearrangement and cell compartment interplay. Interestingly, most of the up-regulated proteins in the cybrid are involved in crucial biological processes such as photosynthesis, energy production and stress tolerance response. The cybrid differential proteome pattern was concomitant with a general increase of leaf gas exchange and content of volatile organic compounds, highlighting a stimulation of specific pathways that can be related to observed plant performances. Our results contribute to a better understanding how the alloplasmic condition might lead to a substantial improvement in plant breeding, opening new opportunities to develop varieties more adaptable to a wide range of conditions.

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.

QTL mapping of male sterility and transmission pattern in progeny of Satsuma mandarin

Seedlessness is one of the important traits in citrus breeding. Male sterility derived from Satsuma mandarin (Citrus unshiu) has been used in Japanese citrus breeding programs to obtain seedless cultivars. The efficiency of seedless cultivar breeding would be improved by developing a selection marker linked to seedlessness. In this study, we performed QTL mapping in ‘Okitsu No. 46’ × ‘Okitsu No. 56’ (O46-O56) crosses for the number of pollen grains per anther (NPG) and apparent pollen fertility (APF), two traits used as an index of male sterility, and detected a candidate QTL for NPG (MS-P1) on linkage group 8 with a significant LOD score (7.31) and 47% of variance explained. The QTL for APF (MS-F1) was detected on linkage group 6 with a significant LOD score (5.71) and 63.6% of variance explained. The role of both MS-P1 in reducing NPG and MS-F1 in decreasing APF were confirmed with the ‘Okitsu No.46’ × ‘Kara’ (O46-K) cross. Pedigree analysis inferred that both MS-P1 and MS-F1 in ‘Okitsu No. 46’ were derived from kunenbo (Citrus nobilis) through hassaku (C. hassaku) and ‘Sweet Spring’. Cytoplasm analysis revealed that both male-sterile ‘Sweet Spring’ and ‘Okitsu No. 46’ have cytoplasm derived from Kishu (C. kinokuni hort. ex Tanaka), but the cytoplasm of male-sterile kunenbo and hassaku were derived from other varieties rather than Kishu. These results suggest that MS-P1 and MS-F1 primarily reduce the NPG and decrease APF, but their expression requires a cytoplasm derived from Kishu. These findings will improve our understanding of the molecular mechanism of male sterility in citrus and help to develop a DNA marker for seedless breeding in citrus.

Comparative transcript profiling reveals the mechanism of female sterility associated with seedless Ponkan mandarin (Citrus reticulata Blanco)

Seedlessness is a highly desirable trait in citrus varieties. Sterility is the key determination for seedlessness formation. However, the molecular basis for female sterility in seedless mandarin remains unclear. Thus, a seedless Ponkan mandarin (Citrus reticulata Blanco 'Lipeng No.2'), considered the bud mutation of normal seedy Ponkan, was collected to identify candidate genes involved in seedless variation. Suppression subtractive hybridization (SSH) screened 1091 uniESTs related to seedy and seedless Ponkan (727 singlets and 364 contigs), which mainly governed catalytic activity, transferase activity, and oxygen binding. By using RNA-Seq technology, 106 differentially expressed genes (DEGs) were captured, of which 74 were up-regulated and 32 were down-regulated. Gene Ontology and pathway analysis showed that six DEGs were enriched in the biosynthesis of secondary metabolite, whereas five DEGs were enriched in the signaling of plant hormones. The combined results of SSH and RNA-Seq indicated the importance of amino acid metabolism in seedless Ponkan. Our findings revealed that the mechanism of seedless Ponkan generation may be related to gene regulation, signal cascade, and hormone levels. This study provided a solid foundation for functional gene identification in seedless Ponkan and a good reference for relevant research on molecular mechanisms of female sterility in Ponkan mandarin.

Exploration of miRNAs and target genes of cytoplasmic male sterility line in cotton during flower bud development

Cytoplasmic male sterility (CMS) lines provide crucial material to harness heterosis for crop plants, which serves as an important strategy for hybrid seed production. However, the molecular mechanism remains obscure. Although microRNAs (miRNAs) play important roles in vegetative growth and reproductive growth, there are few reports on miRNAs regulating the development of male sterility in Upland cotton. In present study, 12 small RNA libraries were constructed and sequenced for two development stages of flower buds from a CMS line and its maintainer line. Based on the results, 256 novel miRNAs were allocated to 141 new miRNA families, and 77 known miRNAs belonging to 54 conserved miRNA families were identified as well. Comparative analysis revealed that 61 novel and 10 conserved miRNAs were differentially expressed. Further transcriptome analysis identified 232 target genes for these miRNAs, which participated in cellular developmental process, cell death, pollen germination, and sexual reproduction. In addition, expression patterns of typical miRNA and the negatively regulated target genes, such as PPR, ARF, AP2, and AFB, were verified by qRT-PCR in cotton flower buds. These targets were previously reported to be related to reproduction development and male sterility, suggesting that miRNAs might act as regulators of CMS occurrence. Some miRNAs displayed specific expression profiles in special developmental stages of CMS line and its fertile hybrid (F₁). Present study offers new information on miRNAs and their related target genes in exploiting CMS mechanism, and revealing the miRNA regulatory networks in Upland cotton.

Non-coding RNAs and plant male sterility: current knowledge and future prospects

Latest outcomes assign functional role to non-coding (nc) RNA molecules in regulatory networks that confer male sterility to plants. Male sterility in plants offers great opportunity for improving crop performance through application of hybrid technology. In this respect, cytoplasmic male sterility (CMS) and sterility induced by photoperiod (PGMS)/temperature (TGMS) have greatly facilitated development of high-yielding hybrids in crops. Participation of non-coding (nc) RNA molecules in plant reproductive development is increasingly becoming evident. Recent breakthroughs in rice definitively associate ncRNAs with PGMS and TGMS. In case of CMS, the exact mechanism through which the mitochondrial ORFs exert influence on the development of male gametophyte remains obscure in several crops. High-throughput sequencing has enabled genome-wide discovery and validation of these regulatory molecules and their target genes, describing their potential roles performed in relation to CMS. Discovery of ncRNA localized in plant mtDNA with its possible implication in CMS induction is intriguing in this respect. Still, conclusive evidences linking ncRNA with CMS phenotypes are currently unavailable, demanding complementing genetic approaches like transgenics to substantiate the preliminary findings. Here, we review the recent literature on the contribution of ncRNAs in conferring male sterility to plants, with an emphasis on microRNAs. Also, we present a perspective on improved understanding about ncRNA-mediated regulatory pathways that control male sterility in plants. A refined understanding of plant male sterility would strengthen crop hybrid industry to deliver hybrids with improved performance.

The Role of Non-Coding RNAs in Cytoplasmic Male Sterility in Flowering Plants

The interactions between mitochondria and nucleus substantially influence plant development, stress response and morphological features. The prominent example of a mitochondrial-nuclear interaction is cytoplasmic male sterility (CMS), when plants produce aborted anthers or inviable pollen. The genes responsible for CMS are located in mitochondrial genome, but their expression is controlled by nuclear genes, called fertility restorers. Recent explosion of high-throughput sequencing methods enabled to study transcriptomic alterations in the level of non-coding RNAs under CMS biogenesis. We summarize current knowledge of the role of nucleus encoded regulatory non-coding RNAs (long non-coding RNA, microRNA as well as small interfering RNA) in CMS. We also focus on the emerging data of non-coding RNAs encoded by mitochondrial genome and their possible involvement in mitochondrial-nuclear interactions and CMS development.

Distinct transcriptional and processing regulations control miR167a level in tomato during stress

Besides their definite role in plant developmental processes miR167 also serve as mediator of stress response. Although differential expression of miR167 occurs during stresses, the regulatory-mechanism of biogenesis remained elusive. Therefore, using tomato as the model plant we have explored the mechanism of regulation of miR167a expression during stresses. Fungus or virus infections and exposure to cold stress raised the level of miR167a expression. Whereas, salt, drought and heat treatments resulted in the downregulation, indicating different stresses activated alternative mechanisms for miR167a regulation. Interestingly, the relative expression level of precursors in control versus temperature stressed plants differed from the pattern observed in the mature miR167a expression, suggesting that both transcriptional and processing regulation were important for biogenesis. The promoter-regulatory sequence of the major isoform MIR167a harbours several development and stress-related regulatory sites. Accordingly, promoter assays using transient transformation and transgenic tobacco plants proved stress-dependent regulation of the promoter. Further analyses corroborated the role of tomato DREB2A protein in the transcriptional regulation during temperature stress. Finally, in vitro assays established the importance of processing factors in cold-stress dependent efficient processing of MIR167a precursors. These data confirm distinct role of transcriptional and processing machinery in stress-influenced regulation of tomato miR167a biogenesis.

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.

Comparative transcriptome analysis during early fruit development between three seedy citrus genotypes and their seedless mutants

Identification of genes with differential transcript abundance (GDTA) in seedless mutants may enhance understanding of seedless citrus development. Transcriptome analysis was conducted at three time points during early fruit development (Phase 1) of three seedy citrus genotypes: Fallglo (Bower citrus hybrid (Citrus reticulata×C. reticulata×C. paradisi)×Temple (C. reticulata×C. sinensis)), grapefruit (C. paradisi), Pineapple sweet orange (C. sinensis), and their seedless mutants. Seed abortion in seedless mutants was observed at 26 days post anthesis (Time point 2). Affymetrix transcriptomic analysis revealed 359 to 1077 probe sets with differential transcript abundance in the comparison of seedless versus seedy fruits for each citrus genotypes and time points. The GDTA identified by 18 microarray probe sets were validated by qPCR. Hierarchical clustering analysis revealed a range of GDTA associated with development, hormone and protein metabolism, all of which may reflect genes associated with seedless fruit development. There were 14, 9 and 12 genes found exhibiting similar abundance ratios in all three seedless versus seedy genotype comparisons at time point 1, 2 and 3, respectively. Among those genes were genes coding for an aspartic protease and a cysteine protease, which may play important roles in seedless fruit development. New insights into seedless citrus fruit development may contribute to biotech approaches to create seedless cultivars.