miRNAs and lncRNAs in reproductive development

Whole-transcriptome analyses identify key differentially expressed mRNAs, lncRNAs, and miRNAs associated with male sterility in watermelon

Male sterility is a valuable trait for watermelon breeding, as watermelon hybrids exhibit obvious heterosis. However, the underlying regulatory mechanism is still largely unknown, especially regarding the related non-coding genes. In the present study, approximately 1035 differentially expressed genes (DEGs), as well as 80 DE-lncRNAs and 10 DE-miRNAs, were identified, with the overwhelming majority down-regulated in male-sterile floral buds. Enrichment analyses revealed that the general phenylpropanoid pathway as well as its related metabolisms was predicted to be altered in a mutant compared to its fertile progenitor. Meanwhile, the conserved genetic pathway DYT1-TDF1-AMS-MS188-MS1, as well as the causal gene ClAMT1 for the male-sterile mutant Se18, was substantially disrupted during male reproductive development. In addition, some targets of the key regulators AMS and MS188 in tapetum development were also down-regulated at a transcriptional level, such as ABCG26 (Cla004479), ACOS5 (Cla022956), CYP703A2 (Cla021151), PKSA (Cla021099), and TKPR1 (Cla002563). Considering lncRNAs may act as functional endogenous target mimics of miRNAs, competitive endogenous RNA networks were subsequently constructed, with the most complex one containing three DE-miRNAs, two DE-lncRNAs, and 21 DEGs. Collectively, these findings not only contribute to a better understanding of genetic regulatory networks underlying male sterility in watermelon, but also provide valuable candidates for future research.

A combined transcriptome - miRNAome approach revealed that a kinesin gene is differentially targeted by a novel miRNA in an apomictic genotype of Eragrostis curvula

Weeping lovegrass (Eragrostis curvula [Shrad.] Nees) is a perennial grass typically established in semi-arid regions, with good adaptability to dry conditions and sandy soils. This polymorphic complex includes both sexual and apomictic cytotypes, with different ploidy levels (2x-8x). Diploids are known to be sexual, while most polyploids are facultative apomicts, and full apomicts have also been reported. Plant breeding studies throughout the years have focused on achieving the introgression of apomixis into species of agricultural relevance, but, given the complexity of the trait, a deeper understanding of the molecular basis of regulatory mechanisms of apomixis is still required. Apomixis is thought to be associated with silencing or disruption of the sexual pathway, and studies have shown it is influenced by epigenetic mechanisms. In a previous study, we explored the role of miRNA-mRNA interactions using two contrasting E. curvula phenotypes. Here, the sexual OTA-S, the facultative Don Walter and the obligate apomictic Tanganyika cDNA and sRNA libraries were inquired, searching for miRNA discovery and miRNA expression regulation of genes related to the reproductive mode. This allowed for the characterization of seven miRNAs and the validation of their miRNA-target interactions. Interestingly, a kinesin gene was found to be repressed in the apomictic cultivar Tanganyika, targeted by a novel miRNA that was found to be overexpressed in this genotype, suggestive of an involvement in the reproductive mode expression. Our work provided additional evidence of the contribution of the epigenetic regulation of the apomictic pathway.

Unrevealing functional candidate genes for bovine fertility through RNA sequencing meta-analysis and regulatory elements networks of co-expressed genes and lncRNAs

The high genetic heterogeneity and environmental effects of subfertility in livestock species make the elucidation of the genetic mechanisms associated with reproductive efficiency a difficult task. Network and co-expression network meta-analyses were applied alongside genetic variant calling and long non-coding RNA (lncRNA) characterization to identify functionally relevant target genes and regulatory subnetworks associated with fertility in dairy cattle. In total, 505 lncRNAs (441 previously annotated in the bovine reference genome ARS-UCD 1.2 and 64 novel lncRNAs) were identified. Seven differentially expressed genes between high-fertile (HF) and sub-fertile (SF) Holstein cows were identified in the network meta-analysis (CA5A, ENSBTAG00000051149, ENSBTAG00000003272, DEFB7, DIO2, TRPV3, and COL4A4). Additionally, seven functional candidate differentially co-expressed (DcoExp) modules with a differential regulatory pattern (|z-score|>2) were identified between HF and SF cows. The functional candidate genes and DcoExp modules identified were associated with fertility relevant processes such as the regulation of embryonic implantation and proliferation, interaction and molecule transfer between the fetus and the cow, and the immune system. These results help to better understand the genetic mechanisms associated with reproductive efficiency in dairy cattle through the identification of potential biomarkers and genetic variants associated with differentially expressed regulatory gene and lncRNAs regulatory element networks.

Reproductive tissue-specific translatome of a rice thermo-sensitive genic male sterile line

Translational regulation, especially tissue- or cell type-specific gene regulation, plays essential roles in plant growth and development. Thermo-sensitive genic male sterile (TGMS) lines have been widely used for hybrid breeding in rice (Oryza sativa). However, little is known about translational regulation during reproductive stage in TGMS rice. Here, we used translating ribosome affinity purification (TRAP) combined with RNA sequencing to investigate the reproductive tissue-specific translatome of TGMS rice expressing FLAG-tagged ribosomal protein L18 (RPL18) from the germline-specific promoter MEIOSIS ARRESTED AT LEPTOTENE1 (MEL1). Differentially expressed genes at the transcriptional and translational levels were enriched in pollen and anther-related formation and development processes. These contained a number of genes reported to be involved in tapetum programmed cell death (PCD) and lipid metabolism during pollen development and anther dehiscence in rice, including several encoding transcription factors and key enzymes, as well as several long non-coding RNAs (lncRNAs) that potentially affect tapetum and pollen-related genes in male sterility. This study represents the first comprehensive reproductive tissue-specific characterization of the translatome in TGMS rice. These results contribute to our understanding of the molecular basis of sterility in TGMS rice and will facilitate further genetic manipulation of TGMS rice in two-line breeding systems.

Comprehensive analysis of formin gene family highlights candidate genes related to pollen cytoskeleton and male fertility in wheat (Triticum aestivum L.)

Background

Formin, a highly conserved multi-domain protein, interacts with microfilaments and microtubules. Although specifically expressed formin genes in anthers are potentially significant in research on male sterility and hybrid wheat breeding, similar reports in wheat, especially in thermo-sensitive genic male sterile (TGMS) wheat, remain elusive.

Results

Herein, we systematically characterized the formin genes in TGMS wheat line BS366 named TaFormins (TaFHs) and predicted their functions in inducing stress response. In total, 25 TaFH genes were uncovered, majorly localized in 2A, 2B, and 2D chromosomes. According to the neighbor-joining (NJ) method, all TaFH proteins from wheat and other plants clustered in 6 sub-groups (A-F). The modeled 3D structures of TaFH1-A/B, TaFH2-A/B, TaFH3-A/B and TaFH3-B/D were validated. And different numbers of stress and hormone-responsive regulatory elements in their 1500 base pair promoter regions were contained in the TaFH genes copies. TaFHs had specific temporal and spatial expression characteristics, whereby TaFH1, TaFH4, and TaFH5 were expressed highly in the stamen of BS366. Besides, the accumulation of TaFHs was remarkably lower in a low-temperature sterile condition (Nanyang) than fertile condition (Beijing), particularly at the early stamen development stage. The pollen cytoskeleton of BS366 was abnormal in the three stages under sterile and fertile environments. Furthermore, under different stress levels, TaFHs expression could be induced by drought, salt, abscisic acid (ABA), salicylic acid (SA), methyl jasmonate (MeJA), indole-3-acetic acid (IAA), polyethylene glycol (PEG), and low temperature. Some miRNAs, including miR167, miR1120, and miR172, interacts with TaFH genes; thus, we constructed an interaction network between microRNAs, TaFHs, phytohormone responses, and distribution of cytoskeleton to reveal the regulatory association between upstream genes of TaFH family members and sterile.

Conclusions

Collectively, this comprehensive analysis provides novel insights into TaFHs and miRNA resources for wheat breeding. These findings are, therefore, valuable in understanding the mechanism of TGMS fertility conversion in wheat.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07878-7.

Genome-wide analysis of long noncoding RNAs affecting floral bud dormancy in pears in response to cold stress

The versatile role of long noncoding RNAs (lncRNAs) in plant growth and development has been established, but a systematic identification and analysis of lncRNAs in the pear has not been reported. Bud dormancy is a crucial and complicated protective mechanism for plants in winter. The roles of lncRNAs in the dormancy process remain largely unclear. In this study, we induced pear floral buds to enter into different dormant statuses by simulating four different chilling accumulation conditions. Then, a time series of RNA-seq analysis was performed and we identified 7594 lncRNAs in Pyrus pyrifolia (Burm. F.) Nakai that have not been identified. The sequence and expression of the lncRNAs were confirmed by PCR analysis. In total, 6253 lncRNAs were predicted to target protein-coding genes including 692 cis-regulated pairs (596 lncRNAs) and 13,158 trans-regulated pairs (6181 lncRNAs). Gene Ontology analysis revealed that most of lncRNAs' target genes were involved in catalytic activity, metabolic processes and cellular processes. In the trend analysis, 124 long-term cold response lncRNAs and 80 short-term cold response lncRNAs were predicted. Regarding the lncRNA-miRNA regulatory networks, 59 lncRNAs were identified as potential precursors for miRNA members of 20 families, 586 lncRNAs were targets of 261 pear miRNAs and 53 lncRNAs were endogenous target mimics for 26 miRNAs. In addition, three cold response lncRNAs, two miRNAs and their target genes were selected for expression confirmed. The trend of their expression was consistent with the predicted relationships among them and suggested possible roles of lncRNAs in ABA metabolic pathway. Our findings not only suggest the potential roles of lncRNAs in regulating the dormancy of pear floral buds but also provide new insights into the lncRNA-miRNA-mRNA regulatory network in plants.

Systematic identification and characterization of circular RNAs involved in flag leaf senescence of rice

Circular RNAs (circRNAs) identification, expression profiles, and construction of circRNA-parental gene relationships and circRNA-miRNA-mRNA ceRNA networks indicate that circRNAs are involved in flag leaf senescence of rice. Circular RNAs (circRNAs) are a class of 3'-5' head-to-tail covalently closed non-coding RNAs which have been proved to play important roles in various biological processes. However, no systematic identification of circRNAs associated with leaf senescence in rice has been studied. In this study, a genome-wide high-throughput sequencing analysis was performed using rice flag leaves developing from normal to senescence. Here, a total of 6612 circRNAs were identified, among which, 113 circRNAs were differentially expressed (DE) during the leaf senescence process. Moreover, 4601 (69.59%) circRNAs were derived from the exons or introns of their parental genes, while 2110 (71%) of the parental genes produced only one circRNA. The sequence alignment analysis showed that hundreds of rice circRNAs were conserved among different plant species. Gene Ontology (GO) enrichment analysis revealed that parental genes of DE circRNAs were enriched in many biological processes closely related to leaf senescence. Through weighted gene co-expression network analysis (WGCNA), six continuously down-expressed circRNAs, 18 continuously up-expressed circRNAs and 15 turn-point high-expressed circRNAs were considered to be highly associated with leaf senescence. Additionally, a total of 17 senescence-associated circRNAs were predicted to have parental genes, in which, regulations of three circRNAs to their parental genes were validated by qRT-PCR. The competing endogenous RNA (ceRNA) networks were also constructed. And a total of 11 senescence-associated circRNAs were predicted to act as miRNA sponges to regulate mRNAs, in which, regulation of two circRNAs to eight mRNAs was validated by qRT-PCR. It is discussed that senescence-associated circRNAs were involved in flag leaf senescence probably through mediating their parental genes and ceRNA networks, to participate in several well-studied senescence-associated processes, mainly including the processes of transcription, translation, and posttranslational modification (especially protein glycosylation), oxidation-reduction process, involvement of senescence-associated genes, hormone signaling pathway, proteolysis, and DNA damage repair. This study not only showed the systematic identification of circRNAs involved in leaf senescence of rice, but also laid a foundation for functional research on candidate circRNAs.

Uncloaking lncRNA-meditated gene expression as a potential regulator of CMS in cotton (Gossypium hirsutum L.)

Cytoplasmic male sterility is a well-proven mechanism for cotton hybrid production. Long non-coding RNAs belong to a class of transcriptional regulators that function in multiple biological processes. The cDNA libraries from the flower buds of the cotton CGMS, it's restorer (Rf) and maintainer lines were sequenced using high throughput NGS technique. A total of 1531 lncRNAs showed significant differential expression patterns between these three lines. Functional analysis of the co-expression network of lncRNA-mRNA using gene ontology vouchsafes that, lncRNAs play a crucial role in cytoplasmic male sterility and fertility restoration through pollen development, INO80 complex, development of anther wall tapetum, chromatin remodeling, and histone modification. Additionally, 94 lncRNAs were identified as putative precursors of 49 miRNAs. qRT-PCR affirms the concordance of expression pattern to RNA-seq data. These findings divulge the lncRNA driven miRNA-mediated regulation of gene expression profiling superintended for a better understanding of the CMS mechanisms of cotton.

Elucidating the Functional Role of Predicted miRNAs in Post- Transcriptional Gene Regulation Along with Symbiosis in Medicago truncatula

Background:

microRNAs are small non-coding RNAs which inhibit translational and post-transcriptional processes whereas long non-coding RNAs are found to regulate both transcriptional and post-transcriptional gene expression. Medicago truncatula is a well-known model plant for studying legume biology and is also used as a forage crop. In spite of its importance in nitrogen fixation and soil fertility improvement, little information is available about Medicago non-coding RNAs that play important role in symbiosis.

Objective:

In this study we have tried to understand the role of Medicago ncRNAs in symbiosis and regulation of transcription factors.

Methods:

We have identified novel miRNAs by computational methods considering various parameters like length, MFEI, AU content, SSR signatures and tried to establish an interaction model with their targets obtained through psRNATarget server.

Results:

149 novel miRNAs are predicted along with their 770 target proteins. We have also shown that 51 of these novel miRNAs are targeting 282 lncRNAs.

Conclusion:

In this study role of Medicago miRNAs in the regulation of various transcription factors are elucidated. Knowledge gained from this study will have a positive impact on the nitrogen fixing ability of this important model plant, which in turn will improve the soil fertility.

Variation in plastid genomes in the gynodioecious species Silene vulgaris

BACKGROUND

Gynodioecious species exist in two sexes - male-sterile females and hermaphrodites. Male sterility in higher plants often results from mitonuclear interaction between the CMS (cytoplasmic male sterility) gene(s) encoded by mitochondrial genome and by nuclear-encoded restorer genes. Mitochondrial and nuclear-encoded transcriptomes in females and hermaphrodites are intensively studied, but little is known about sex-specific gene expression in plastids. We have compared plastid transcriptomes between females and hermaphrodites in two haplotypes of a gynodioecious species Silene vulgaris with known CMS candidate genes.

RESULTS

We generated complete plastid genome sequences from five haplotypes S. vulgaris including the haplotypes KRA and KOV, for which complete mitochondrial genome sequences were already published. We constructed a phylogenetic tree based on plastid sequences of S. vulgaris. Whereas lowland S. vulgaris haplotypes including KRA and KOV clustered together, the accessions from high European mountains diverged early in the phylogram. S. vulgaris belongs among Silene species with slowly evolving plastid genomes, but we still detected 212 substitutions and 112 indels between two accessions of this species. We estimated elevated Ka/Ks in the ndhF gene, which may reflect the adaptation of S. vulgaris to high altitudes, or relaxed selection. We compared depth of coverage and editing rates between female and hermaphrodite plastid transcriptomes and found no significant differences between the two sexes. We identified 51 unique C to U editing sites in the plastid genomes of S. vulgaris, 38 of them in protein coding regions, 2 in introns, and 11 in intergenic regions. The editing site in the psbZ gene was edited only in one of two plastid genomes under study.

CONCLUSIONS

We revealed no significant differences between the sexes in plastid transcriptomes of two haplotypes of S. vulgaris. It suggests that gene expression of plastid genes is not affected by CMS in flower buds of S. vulgaris, although both sexes may still differ in plastid gene expression in specific tissues. We revealed the difference between the plastid transcriptomes of two S. vulgaris haplotypes in editing rate and in the coverage of several antisense transcripts. Our results document the variation in plastid genomes and transcriptomes in S. vulgaris.

miR-374a/Myc axis modulates iron overload-induced production of ROS and the activation of hepatic stellate cells via TGF-β1 and IL-6

The transformation of hepatic stellate cells (HSCs) to activated myofibroblasts plays a critical role in the progression of hepatic fibrosis, while iron-catalyzed production of free radical, including reaction and active oxygen (ROS), and activation and transformation of HSC into a myofibroblasts has been regarded as a major mechanism. In the present study, we attempted to investigate the mechanism of iron overload in hepatic fibrosis from the perspective of regulating HSC activation via oxidative stress and miR-374a/Myc axis. FAC stimulation significantly increased ROS production and TGF-β1 and IL-6 release dose-dependently in hepatocytes. miR-374a could target Myc, a co-transcription factor of both TGF-β1 and IL-6, to negatively regulate Myc expression; FAC stimulation significantly suppressed miR-374a expression, whereas the suppressive effect of FAC stimulation on miR-374a expression could be reversed by ROS inhibitor NAC, indicating that miR-374a could be modulated by iron overload-induced ROS. Via targeting Myc, miR-374a overexpression significantly reduced FAC-induced increases in TGF-β1 and IL-6 levels within L02 cells, whereas the effects of miR-374a overexpression were significantly attenuated via Myc overexpression. Finally, miR-374a overexpression attenuated FAC-induced activity of HSCs by decreasing α-SMA and Collagen I levels whereas Myc overexpression enhanced FAC-induced activity of HSCs by increasing α-SMA and Collagen I levels; the effects of miR-374a overexpression could also be significantly reversed by Myc overexpression, indicating that miR-374a suppresses the activation of HSCs by inhibiting Myc to reduce FAC-induced increases in TGF-β1 and IL-6 release. In conclusion, we demonstrate a novel mechanism of miR-374a/Myc axis modulating iron overload-induced production of ROS and the activation of HSCs via TGF-β1 and IL-6.

Network analysis of noncoding RNAs in pepper provides insights into fruit ripening control

Pepper is an important vegetable worldwide and is a model plant for nonclimacteric fleshy fruit ripening. Drastic visual changes and internal biochemical alterations are involved in fruit coloration, flavor, texture, aroma, and palatability to animals during the pepper fruit ripening process. To explore the regulation of bell pepper fruit ripening by noncoding RNAs (ncRNAs), we examined their expression profiles; 43 microRNAs (miRNAs), 125 circular RNAs (circRNAs), 366 long noncoding RNAs (lncRNAs), and 3266 messenger RNAs (mRNAs) were differentially expressed (DE) in mature green and red ripe fruit. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses revealed that the targets of the DE ncRNAs and DE mRNAs included several kinds of transcription factors (TFs) (ERF, bHLH, WRKY, MYB, NAC, bZIP, and ARF), enzymes involved in cell wall metabolism (beta-galactosidase, beta-glucosidase, beta-amylase, chitinase, pectate lyase (PL), pectinesterase (PE) and polygalacturonase (PG)), enzymes involved in fruit color accumulation (bifunctional 15-cis-phytoene synthase, 9-cis-epoxycarotenoid dioxygenase, beta-carotene hydroxylase and carotene epsilon-monooxygenase), enzymes associated with fruit flavor and aroma (glutamate-1-semialdehyde 2,1-aminomutase, anthocyanin 5-aromatic acyltransferase, and eugenol synthase 1) and enzymes involved in the production of ethylene (ET) (ACO1/ACO4) as well as other plant hormones such as abscisic acid (ABA), auxin (IAA), and gibberellic acid (GA). Based on accumulation profiles, a network of ncRNAs and mRNAs associated with bell pepper fruit ripening was developed that provides a foundation for further developing a more refined understanding of the molecular biology of fruit ripening.

The multiple roles of OsmiR535 in modulating plant height, panicle branching and grain shape

The miR156/miR529-SPL module acts a vital role in regulating plant growth and development. Though miR535 shows very high sequence identity to miR156 and miR529, it is still unknown whether miR535 could control plant growth and development. In this study, we performed the evolutionary analyses of miR535s in land plants and found that miR535s were less conserved than miR156s during evolution. In rice, miR535 expressed at a very low level during the vegetative growth but highly accumulated in young panicles, which is similar with OsmiR529, but opposite to OsmiR156. Expectedly, OsmiR535 overexpression in rice reduced plant height by decreasing the 1^st and 2^nd internode length. Furthermore, OsmiR535 overexpression imposed great influence in panicle architecture, such as more but shorter panicles, and fewer primary/secondary panicle branches. Moreover, OsmiR535 overexpression increased the grain length, but did not affect grain width. Through quantitative real-time PCR analyses, we further revealed that OsmiR535 overexpression repressed the expression of OsSPL7/12/16, as well as the OsSPLs downstream panicle related genes, including OsPIN1B, OsDEP1, OsLOG and OsSLR1. Taken together, our findings suggest that OsmiR535 multiply modulates plant height, panicle architecture and grain shape possibly by regulating OsSPLs genes in rice.

Identification and characterization of circular RNAs during the sea buckthorn fruit development

As a rising star of noncoding RNA, circular RNAs (circRNAs) have a covalently closed loop structure, which formed by 3'-5' ligation during splicing. A few circRNAs were identified and thought to be transcriptional noise due to cognitive defect over the past 40 years. Recently, with the development of high-throughput RNA sequencing techniques and specific algorithms for circRNA detection and quantification, plenty of potential circRNAs were identified in many species which play important roles in various biological processes. However, researches on circRNAs in fruit ripening process were lacking. Here, we totally identified 2616 circRNAs in sea buckthorn fruit development process, which uniformly distributed in sea buckthorn chromosome. Among them, 1721 (65.8%) circRNAs were arising from the exons of their host genes, 252 circRNAs were identified as the differentially expressed circRNAs (DEcircRNAs) between three different development stages, and 181 (71.8%) DEcircRNAs had sequence similarity with 235 identified circRNAs from five know plant species. Functional annotation revealed that host genes of DEcircRNAs were predicted to be involved in carotenoid biosynthesis, lipid synthesis and plant hormone signal transduction. Additionally, 53 DEcircRNAs were predicted as the corresponding nine miRNAs sponges in sea buckthorn. Divergent reverse-transcription PCR and RT-qPCR were used for validate the differential expression and back-splicing sites of six DEcircRNAs. These results revealed the role of circRNAs in sea buckthorn fruit ripening process and promoted the noncoding RNA researches in plants.

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.

Long Noncoding RNA CAT104 Promotes Cell Viability, Migration, and Invasion in Gastric Carcinoma Cells Through Activation of MicroRNA-381-Inhibiting Zinc Finger E-box-Binding Homeobox 1 (ZEB1) Expression

Gastric carcinoma (GC) remains the second leading cause of cancer-related deaths worldwide. Good biomarkers are of paramount importance for GC therapy. This study aimed to assess the role of long noncoding RNA (lncRNA) CAT104 in GC. We found that CAT104 was highly expressed in human GC NCI-N87, SGC7901, BGC823, BGC803, and AGS cells. Suppression of CAT104 decreased NCI-N87 cell viability, migration, and invasion, but promoted apoptosis. CAT104 knockdown enhanced the expression of microRNA-381 (miR-381) expression in NCI-N87 cells. miR-381 participated in the regulatory effects of CAT104 on NCI-N87 cell viability, migration, invasion, and apoptosis. Zinc finger E-box-binding homeobox 1 (ZEB1) was identified as a direct target of miR-381. Overexpression of ZEB1 reversed the miR-381 mimic-induced cell viability, migration, and invasion inhibition. Suppression of ZEB1 reversed the miR-381 inhibitor-induced activation of the c-Jun N-terminal kinase (JNK) pathway and Wnt/β-catenin signaling pathways in NCI-N87 cells. In conclusion, CAT104 might function as an oncogenic factor in GC cells via regulating the expression of miR-381 and ZEB1.

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.

Structure, target-specificity and expression of PN_LNC_N13, a long non-coding RNA differentially expressed in apomictic and sexual Paspalum notatum

KEY MESSAGE: ncRNA PN_LNC_N13 shows contrasting expression in reproductive organs of sexual and apomictic Paspalum notatum genotypes. Apomictic plants set genetically maternal seeds whose embryos derive by parthenogenesis from unreduced egg cells, giving rise to clonal offspring. Several Paspalum notatum apomixis related genes were identified in prior work by comparative transcriptome analyses. Here, one of these candidates (namely N13) was characterized. N13 belongs to a Paspalum gene family including 30-60 members, of which at least eight are expressed at moderate levels in florets. The sequences of these genes show no functional ORFs, but include segments of different protein coding genes. Particularly, N13 shows partial identity to maize gene BT068773 (RESPONSE REGULATOR 6). Secondary structure predictions as well as mature miRNA and target cleavage detection suggested that N13 is not a miRNA precursor. Moreover, N13 family members produce abundant 24-nucleotide small RNAs along extensive parts of their sequences. Surveys in the GREENC and CANTATA databases indicated similarity with plant long non-coding RNAs (lncRNAs) involved in splicing regulation; consequently, N13 was renamed as PN_LNC_N13. The Paspalum BT068773 predicted ortholog (N13TAR) originates floral transcript variants shorter than the canonical maize isoform and with possible structural differences between the apomictic and sexual types. PN_LNC_N13 is expressed only in apomictic plants and displays quantitative representation variation across reproductive developmental stages. However, PN_LNC_N13-like homologs and/or its derived sRNAs showed overall a higher representation in ovules of sexual plants at late premeiosis. Our results suggest the existence of a whole family of N13-like lncRNAs possibly involved in splicing regulation, with some members characterized by differential activity across reproductive types.

An Integrated Analysis of miRNAs and Methylated Genes Encoding mRNAs and lncRNAs in Sheep Breeds with Different Fecundity

In our previous study, we investigated the regulatory relationship between lncRNAs, miRNA, and mRNAs in an effort to shed light onto the regulatory mechanisms involved in sheep fecundity. As an extension of this study, here, we aimed to identify potential regulators of sheep fecundity using a genome-wide analysis of miRNAs and the methylated genes encoding mRNAs and lncRNAs in the ovaries of Dorset sheep (low fecundity) and Small Tail Han ewes (high fecundity) with the genotype BB (Han BB) and the genotype ++ (Han ++) by performing RNA-Seq and MeDIP-Seq analyses. Methylated coding-non-coding gene co-expression networks for Han and Dorset sheep were constructed using the methylated genes encoding the differentially expressed mRNAs and lncRNAs identified in this study. In the Han BB vs. Dorset comparison, the lncRNAs TTC26 and MYH15 had the largest degree. Similarly, the lncRNA NYAP1 had the largest degree in the Han ++ vs. Dorset comparison. None of the methylated genes encoding lncRNAs were co-expressed with the methylated genes encoding mRNAs in the Han BB vs. Han ++ comparison. The methylated genes encoding lncRNAs identified here may play a vital regulatory role in sheep breeding. Our results suggest that miRNAs might play a key role in sheep prolificacy by regulating target genes related to thyroid hormone synthesis, and methylated genes encoding lncRNAs associated with tight junctions might contribute to the high breeding rate in Han sheep. These findings may contribute to a deeper understanding of sheep prolificacy.

Identification and expression analysis of microRNAs during ovule development in rice (Oryza sativa) by deep sequencing

MicroRNA (miRNA) expression profiles during rice ovule development revealed the possible miRNA-mediated regulation between ovule sporophytic tissue and female gametophyte and the involvement of miRNAs in programmed cell death. MiRNAs are 20-24-nucleotide small RNAs that play key roles in the regulation of many growth and developmental processes in plants. Rice ovule development comprises a series of biological events, which are regulated by complex molecular mechanisms. To gain insight into miRNA-mediated regulation of rice ovule development, Illumina sequencing was used to examine the expression of miRNAs from the megaspore mother cell meiosis stage to the fertilized ovule stage. Based on the sequencing data, 486 known and 204 novel miRNAs were identified during rice ovule development. Moreover, 56, 65 and 11 differentially expressed miRNAs between adjacent developmental stages were identified. By analyzing transcriptome and degradome data, we identified 41, 65 and 12 coherent target genes for the differentially expressed miRNAs in ovule development. We found that changes in the expression of plant hormone-related miRNAs may play important roles in embryo sac development, providing evidence for cross-talk communication between sporophytic tissue and the female gametophyte. Additionally, we revealed that miRNAs may be involved in programmed cell death after fertilization. Finally, we constructed miRNA-mediated regulatory networks that are active during rice ovule development.

Vernalisation mediated LncRNA-like gene expression in Beta vulgaris

Sugar beet (Beta vulgaris L.) cannot form reproductive shoots during the first year of their life cycle. Flowering only occurs if plants are vernalised and are subsequently exposed to long days. However, the vernalisation mechanism remains poorly understood in sugar beet. Three putative lncRNAs associated with vernalisation (AGL15X1, AGL15X2 and CAULIFLOWER A) were investigated and the hypothesis that their expression occurred in response to vernalisation was experimentally tested. The regulation mechanisms of BvRAV1-like, lncRNA-like genes, BvFT1 and BvFT2 were also examined. The BvRAV1-like gene associated with vernalisation in sugar beet was validated for the first time. Our data confirmed the hypothesis that AGLX2 was the first candidate lncRNA of sugar beet and the BvRAV1-like gene was expressed in response to vernalisation. BvRAV1-like and AGLX2 genes might be coordinated with BvFT2 to promote reproductive growth by repressing BvFT1 during cold exposure followed by long day conditions. A new complementary flowering model of sugar beet was proposed. Our findings opened up new possibility for future studies and further illuminated the molecular mechanism of vernalisation in sugar beet.

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.