Skip to main content

Exploration of miRNAs in Neolamarckia cadamba and their potential cross-kingdom functions

Abstract

Neolamarckia cadamba (Roxb.) Bosser (Rubiaceae) is a widespread plant with medicine value and has been used for treating various diseases, such as coughs, fevers, anemia, blood disorders and stomach pains. It had been reported that plant miRNAs might enter mammalian intestines and exert a special ingredient across the different species. However, the knowledge about miRNAs in N.cadamba is scarce. In this study, we first applied high-throughput sequencing to identify miRNAs in N.cadamba. Sequencing revealed a total of 11,285,965 raw reads in the small RNA library of N.cadamba leaves. By bioinformatics analysis, we identified a total of 192 miRNAs, including 167 conserved miRNAs and 25 novel miRNAs, which were distributed among 30 families. Next, we used miRanda to predict the targets of those miRNAs in human, 4030 target genes were predicted. Furthermore, Gene ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses indicated that the identified miRNAs majored in binding, vibrio cholerae infection and insulin signaling, which were likely to provide valuable references for further understanding of medical functions of N. cadamba.

Introduction

The N. cadamba is a member of the tribe Naucleeae in the family Rubiaceae and distributed widely in South Asia and the south of China [28]. According to present data, N. cadamba as a representative medicinal plants which is widely known as India Ayurvedic system of Medicine (IASM) [15] The studies of N. cadamba have mainly been focused on its medicinal values in the treatment of eye infections, antidiarrheal, skin diseases, dyspepsia and the stomatitis, coughs, fevers, anaemia, blood disorders and stomach pains [4, 26]. The identified major active components in bark and leaf are alkaloids, saponins, flavoniuds, terpenoids and chlorogenic acids [20]. Until now, scientific bioactivity identification researches have uncovered its antimalarial [30], anti-hepatotoxic [20], anti-inflammatory [9], antioxidant, wound healing and antimicrobial activities [31]. Narayan’s group found that the application of defatted methanol extract N. cadamba (MENC) on Ehrlich ascites carcinoma (EAC) played a constructive role on antitumor activity and antioxidant activity in vivo [13].

MicroRNAs (miRNAs) are a class small non-coding RNAs, which are 18–24 nucleotide in length and inhibit gene expression by mRNA cleavage or translation repression in the 3′ untranslated region (3’UTR) [8]. They are crucial on gene expressions, and a number of articles indicated that miRNAs not only execute biological functions inherent in their system, but can also be transferred from one species to another. Recent studies have shown that in addition to regulating intracellular levels, microRNAs also affect intercellular levels, even in interspecies [11, 23]. For instance, several years ago, plant microRNAs were discovered in sera and tissues of human and other animals [35]. Further studies showed these plant microRNAs were absorbed from food; one that was relatively high in serum, miR168a, directly targeted low-density lipoprotein receptor adaptor protein 1 (LDLRAP1) in liver cells and decreased the clearance of LDL from the blood. In 2015, Zhen et al. found plant miRNA miR2911 (Lonicera japonica, Honeysuckle) can directly target influenza A viruses to inhibited H5N1 and H7N9 viral replication [36]. Derived plant miR159 significant inhibited breast cancer cell proliferation by targeting transcription factor 7 (TCF7) [12]. An ever increasing number of researches indicated that plant-engineered miRNAs can be absorbed in the mammalian digestive tract and act as a bioactive constituen to target mammalian genes to mediating cross-kingdom regulations [24, 25, 34].

It is well known that a wide range of biological activities are in the N.cadamba And these substances play an important role, Due to a large number of studies showing that plant miRNAs can regulate biological activity across boundaries, miRNAs of N. cadamba may also play a certain role in addition to various roles of active ingredients, however, miRNAs of N.cadamba still remain unknown. Because the yield of huangliang wood leaves is large and easy to get, we choose leaves as raw materials, In the present study, we used high-through put sequencing and qPCR, and identified conserved and novel miRNAs in N.cadamba leaves, and further analyzed the miRNAs functions in N.cadamba via bioinformatics analysis. Next, we used miRanda to predict the targets of those miRNAs in human. This is the first reporte about miRNAs of N.cadamba and will Provide foundation for further understanding of active components in herbs.

Results

Summary of small RNA library data via deep sequencing

Total RNA was extracted and a small RNA library was constructed from N. cadamba leaves. As a result, a total of 11,285,965 raw reads were obtained. With filtering the adapters and low quality reads, 10,878,015 clean reads remained (Table. 1). The length distributions of the small RNAs were shown in Fig. 1 and majority of reads were 21 to 24 nts in length, the main size groups were 24 nt, The amounts of 24 nt small RNAs were approximately 39.42%. After initial processing, the high quality small RNAs reads were mapped to the N. camdaba leaves transcriptome genome using the bowtie software, and the number of total/unique sequences that paired with the transcriptome were 6,935,424/3,328,945, approximately half (48%) of total reads were mapped successfully. Only a small portion of small RNAs could be annotate (530,525 redundant reads, 28,031 unique reads), thus, the unannotated reads were collected for further prediction of novel miRNAs. The top 10 miRNAs occupied a percentage of 68.2%, top 20 miRNAs and top 50 miRNAs has a dominant percentage of 85.9 and 96.4%, respectively (Fig. 2a). Therefore, the miRNAs are highly concentrated in the top 50, as shown in (Fig. 2b), and the top 10 miRNAs were miR159a, miR166u, miR395b, miR395a, miR166a-3p, miR166h-3p, miR166e, novel-1, and novel-4.

Table 1 Summary of reads in small RNA libraries in N. cadamba
Fig. 1
figure 1

Length distribution of small RNAs from N. cadamba leaves. A total of 6,935,424 clean reads were obtained, ranging from 18 to 30 nt. 24 nt length reads were the most abundant

Fig. 2
figure 2

Distribution of miRNA reads and top 10 miRNAs. a The distribution of miRNA reads showed that top 10, top 20, top 50 and top 100 miRNAs accounted for 68.2, 85.9 and 96.4% of total reads. b Cumulative proportions of top 10 miRNAs. The miR-159a ranked first, accounting for 10.61% of total reads. And the miRNA166 family accounted for four of them

Identification and characterization of the conserved miRNAs

To identify the conserved miRNAs in N. cadamba, the small RNA sequences were compared to known plant mature miRNAs in miRBase database (version 21). One hundred sixty-seven known miRNAs that distributed in 30 families were identified. The most abundant conserved miRNAs were miR159 (63,739 reads) and miR166u (55,327 reads), followed by miR395b (48,862 reads) (Additional file 1: Supplemental Table S1). We analyzed miRNA families of N. cadamba and found miR-159 and miR-166 represented large families with 19 and 18 members, respectively.

Identified of potentially novel miRNAs in N.cadamba

As a result, we detected a total of 25 novel miRNAs, with reads varying from 2 to 24,598 reads (Additional file 1: Supplemental Table S1). The lengths of the novel miRNAs ranged from 20 to 24nts, and the precursors ranged from 45 to 279 nts. The most abundant novel miRNAs were novel-1 with 24,598 reads, followed by novel-4 and novel-8. Although the major novel miRNAs were lower expressed than conserved miRNAs, the species-specific functions they may played must not be overlooked.

Validation of the expression of the predicted miRNAs with RT-PCR

To further validate miRNAs in N.cadamba, we randomly selected seven miRNAs including five conserved miRNAs and two novel miRNAs for stem-loop qRT-PCR detection. The qRT-PCR results showed that seven miRNAs were detected successfully in samples (Fig. 3) and their abundance were also in consistent with Illumina sequencing results.

Fig. 3
figure 3

Quantitative RT-PCR analysis of mature miRNAs in N. cadamba leaves. The expression level of miR-157a was set as control and taken as 1, and the expression levels of other miRNAs were quantified relative to it. The values represent the mean and standard deviation of three independent experiments

Prediction of human gene targets for N.cadamba miRNAs

To investigate whether N.cadamba miRNAs might regulate human gene expression we used miRanda software to search homologous human miRNAs. Result showed 43 miRNAs potentially targeted 4030 genes in human. GO analysis assigned these putative targets into three main categories in terms of biology processes, cellular components, and molecular functions (Additional file 2: Supplemental Table S2). The most enriched GO terms in both species were “cell process” “single-organism cellular process” (BP), “cell”, “cell part” (CC), and “binding” (MF). The predicted target genes were mapped to the KEGG database and classified into 279 and 276 pathways of which 18 and 9 signaling pathways were significantly enriched, respectively (Table 2).

Table 2 Highly enriched KEGG pathways for putative human targets

As shown in Table 2, 18 of the 43 miRNAs might target about 17 related genes (ATP6V0D2, ATP6V0E1, ATP6V0A4, ATP6V1D, SEC61A2, SEC61A1, PRKCB, PRKCA, PRKACA, ADCY3, ADCY9, ARF1, KDELR3, KCNQ1, PDIA4, AGL, SLC12A2) (Fig. 4a) and are likely to be involved in the vibrio cholera infection in human. In addition, a number of 22(EXOC7, SOCS2, INPPL1, CBL, RPS6KB2, PPARGC1A, IRS1, RPTOR, PPP1R3D, G6PC, TSC1, GSK3B, GYS1, MAPK9, MAPK8, PRKACB, PRKAA2, SHC3, INSR, RAPGEF1, PYGB, FASN) (Fig. 4b) genes targeted by N.cadamba miRNAs were concerned with insulin signaling pathway, and 32 targeted genes may be involved in lipogenesis, glycolysis, anti-lipolysis and glycogenesis process, which contributed to anti-diabetic (Additional file 2: Supplemental Table S2).

Fig. 4
figure 4

MiRNAs targeting the genes (In this image, the orange circles show miRNA,the blue circles show target genes). a The 19 miRNAs were analyzed to target 17 vibrio cholera infection genes. b The 34 miRNAs were analyzed to target 22 genes which were concerned with insulin signaling pathway

MiR-166e directly targets FASN

We used miRanda software to identify targets of miR166e in human. A putative binding site conserved among various species was located in FAS (Fig. 5a). To explore this further, a partial FAS normal 3′-UTR sequence and a sequence in which the miR166e binding site is deleted were cloned into a luciferase reporter plasmid (Fig. 5a). HeLa cells were transfected with the reporter plasmids along with synthetic miR166e mimics or NC. The results showed that miR166e significantly reduced luciferase activity, and deletion of the miR166e binding site diminished this reduction (Fig. 5b). These results indicate that miR166e can target FAS.

Fig. 5
figure 5

MiR166e targets FAS. a Prediction of miR166e targets in the FAS 3′-UTR. b HeLa cells were transfected with the pmirGLO-FAS-WT/pmirGLO-FAS-DEL/ pmirGLO-FAS-MUT luciferase reporters and co-transfected with miR166e mimics or NC mimics. The relative luciferase activity was measured 24 h later, n = 6 in each group. Data with error bars mean ± SE *, P < 0.05

Discussion

The miRNAs derived from plants may regulate human genes and have attracted great attention. Zhang and collaborators found rice miRNAs (miR168a) can enter mammalian bloodstream. Later, the miR2911 from honeysuckle were found to target influenza viruses and protect the host. Lately, the miRNA162a can direct target amTOR and influence honeybee caste development [35,36,37]. The plant miR159 significantly inhibited breast cancer cell proliferation by targeting transcription factor 7 (TFC7) [12]. These articles indicated that plant miRNAs could play a vital role in regulating human health.

N.cadamba was used widely in Indian traditional formulations, and it has been a remedy in the treatments of diarrhea, bacterial infection, and diabetes. However, few studies on the breadth of miRNA has yet been reported for N.cadamba, therefore, the knowledge about molecular mechanisms for compounds produced by N.cadamba is still limited. A total of 167 conserved miRNAs and 25 novel miRNAs were identified. From what has been reported about the medicinal plant Moringa oldifera L, plant miRNAs could regulate cell cycle, apoptosis and metabolic in human [27]. Therefore, we predicted the potential targets of N.cadamba miRNAs in human by bioinformatics analysis. The 43 miRNAs of the N.cadamba were targeted to 4030 putative genes in human. GO and KEGG pathways analyses show that those miRNAs may participate in vibrio cholera infection processes.

In human, vibrio cholera infection pathway was highly enrich relative to other pathways. Vibrio cholera is the predominate causes of acute dehydrating diarrhea in developing countries, and the primary virulence factors are cholera toxin (CTX) and toxin-coregulated pilus (TCP), both of which are required for colonization of the human small intestine [14, 19]. Cholera toxin ADP-ribosylates adenylate cyclase, resulting in increased cAMP production and activation of protein kinase A [22, 32].. SLC12A2 can co-express β-cells, and specific co-expression transporters are known to regulate Cl、K+, and KCNQ1 channels, maintaining the Cl secretory response by recycling K+ entering the cell [3, 6]. These relative regulations by vibrio cholerae leads to severe and progressing dehydration and shock. In our study, 18 miRNAs were found involved in vibrio cholera infection. Adenylyl cyclase (AC) was the predicted target of miR166a-3p, miR166u, miR166g, miR166m, miR166h, miR319a, miR396a-5p, miR396g-3p, and miR396h, which can activate cAMP and further influence some key genes that cause diarrhea. The miR396a-5p and miR396h were predicted to target KCNQ1, which were important for Cl secretion in the colon. An increase in secretion of Cl, followed by the movement of water, was able to restore the balance of host that results from diarrhea. The SLC12A2, also known as the Na+-K+-Cl co-transporters, was identified as a possible target of miR396g-3p. We hypothesize that some miRNAs identified in N.cadamba leaves regulate vibrio cholera infection by suppressing the processes related to ion transport. However, their regulatory mechanisms needs further studies.

The dynamic regulation of liver glucose metabolism is vital for systemic carbohydrate homeostasis, and the disorder of insulin secretion could lead to diabetes. Previous studies have reported that various extracts of N.cadamba can reduce the blood glucose level of hyperglycemic mice by increasing the glucose uptake [1, 2, 5, 29]. Based on our bioinformatics prediction, PI3K/Akt pathway were targeted by N.cadamba miRNAs. The high levels of blood glucose will imbalance the insulin secretion and lead to insulin resistance (IR), which is possibly attributed to a defect in the insulin receptor/IRS-1/PI3K/Akt cascade. IRS-1 was found targeted by miR395b-3p, miR858b, novel-27 and novel-8 in human. IRS-1 tyrosine phosphorylation in response to insulin stimulation, increases the association of IRS-1 with PI3-kinase, which leads to activation of serine/threonine kinase protein B (Akt), and p-Akt inhibites the GSK-3β protein directly. The miR159, miR395, miR396 family and novel-1, 4, 8, 23, 28 miRNAs were predicted to target GSK-3β, which phosphorylates and inactivates glycogen synthase (GS) [33]. Therefore, the downstream activated gene encoding GS regulate glucose metabolism through affecting glucose transport and the synthesis of glycogen [17]. Those miRNAs may play a key role in glucose regulation network, and reduce the level of blood glucose. While the direct role of miRNAs in glucose metabolism need to further demonstrated.

Conclusion

We present the first identification of miRNAs in N.cadamba, and use bioinformatic analysis indicates N.cadamba miRNAs are possibly potential medicinal components involved in regulating mammal targets. Our findings display a foundation for further investigation for functions of N.cadamba miRNAs, especially in cross-kingdom manner.

Materials and methods

Plant material

The sample of N.cadamba was collected from the College of Forestry, South China Agricultural University. The fresh leaves were immediately frozen in liquid nitrogen and stored at − 80 °C.

RNA extraction

Frozen tissue was grounded to a fine powder in liquid nitrogen using a mortar and pestle. 100 mg of the powder per sample were transferred into individual RNase-free 1.5 mL tubes containing 600 mL of pre-warmed extraction buffer at 60 °C. The extraction buffer consisted of the following: 2% CTAB, 2% polyvinylpyrrolidone (PVP) K− 40, 100 mmol/L TrisHCl (pH 8.0), 25 mmol/L ethylenediaminetetraacetic acid (EDTA; pH 8.0), 2.0 mol/L NaCl, 2 g/L spermidine and 2% b-mercaptoethanol (added immediately before use). Lysis buffer without spermidine or b-mercaptoethanol was treated with 0.1% diethyl pyrocarbonate (DEPC) and autoclaved. 2 g/L spermidine was added and the mixture was stored at room temperature. The extracts were mixed by vortexing and incubated at 60 °C in a water bath for 10 min with vigorous shaking for several times. Then a third of 5 mol/LKAc (pH 4.8) was added. The supernatant was obtained by ultra-centrifugation at 9000 rpm for 10 min after an ice water bath for 30 min. An equal volume of chloroform/isoamyl alcohol (24,1) was added to the homogenate and was mixed completely by vortexing. The mixture was centrifuged at 12,000 r/min for 5 min at 4 °C. The supernatant was transferred to a new tube with 500 μL Phenol Water (pH 5.2) and the previous step was repeated. The supernatant was then transferred to a new tube containing an equal volume of isopropyl alcohol for 30 min. RNA is acquired by centrifugal at 4 °C at 12,000 rpm for 10 min. An equal volume of ethanol was then applied to wash RNA.

Small RNA library construction and deep sequencing

A total amount of 3 μg total RNA per sample was used as input material for the small RNA library. Sequencing libraries were generated using NEBNext® Multiplex Small RNA Library Prep Set for Illumina® (NEB, USA.) The quantity and purity of total RNA was measured by NanoDrop ND-1000 spectrophotometer (Nano Drop, USA) at 260/280 nm (ratio = 2.0). The integrity of total RNA was confirmed by Bioanalyzer 2100 and RNA 6000Nano LabChip Kit (Agilent, USA) with RIN number [7.0. Briefly, small RNAs of 18–35 nt in length were first isolated from the total RNA by size fractionation. Then, these small RNAs were ligated with 50-RNA and 30-RNA adapters, and subsequently, reverse transcription PCR was used to create cDNAs. The amplified cDNAs were purified and sequenced with Illumina HiSeq 2500 platform.

After deep sequencing, raw data were processed through Novogene Company’s Perl and Python scripts. In this step, clean data were obtained by removing the contaminating reads, sequences containing adapters, without insert tags and reads containing poly A or T or G or C. Sequences from 18 to 30 nt in length were selected for further analysis. Then, the retained reads were searched from NCBI, Rfam, and Repbase database to remove known classes of RNAs (mRNA mRNA, rRNA, tRNA, snRNA, snoRNA, and repeats), so that every unique small RNA is mapped to only one annotation. Since there is no published genome information of N.cadamba, the high-quality reads were mapped to Arabidopsis thaliana genome sequence using the Bowtie software [21]. The mappable small RNA tags were aligned to the miRNA precursor in the miRNA database (miRBase. 21.0; released in June, 2014) to obtain the known miRNA count. Finally, novel miRNAs were predicted by exploring the secondary structure, the dicer cleavage site, and the minimum free energy of the former unannotated small RNA tags which could be mapped to the reference sequence by available software miRDeep 2[16].

Identification of conserved and novel miRNAs

Mapped small RNA tags were used to look for known miRNAs. MiRBase21. 0 was used as reference, and modified software mirdeep2 and srna-tools-cli were used to obtain the potential miRNAs and draw the secondary structures. Custom scripts were used to obtain the miRNA counts as well as base bias on the first position of identified miRNA with certain length and on each position of all identified miRNA respectively. Those lower than 10 reads were discarded. The characteristics of hairpin structures of miRNA precursors can be used to predict novel miRNAs [7]. The available softwares miREvo and mirdeep2 were integrated to predict novel miRNAs by exploring the secondary structures, the Dicer cleavage sites and the minimum free energy (less than − 18 kcal/mol) of the small RNA tags unannotated in the former steps.

Human target gene prediction for the N.cadamba miRNAs

The miRNAs were used for human target prediction. MiRanda was employed to predict putative target genes, and then all the genes were mapped to the KEGG database to identify significantly enriched signal transduction pathways in target genes compared with the whole genome background.

Validation of miRNAs by stem-loop qRT-PCR

To further verify our identification results, 7 randomly chosen sequences (including 2 novel miRNAs as well as 5 conserved ones) were subjected to stem-loop qRT-PCR. The primers for PCR were listed in (Additional file 2: Supplemental Table S2). The stem-loop qRT-PCR steps as previously described [10]. The novel-49 miRNA was used as the internal control miRNA. The reverse primer for miRNAs was the Uni-miR qPCR Primer offered by the kit One Step PrimeScript miRNA cDNA Synthesis Kit (TaKaRa, Dalian). All reactions were performed in triplicate for each sample. The melting curve was generated to test the specificity of PCR products and avoid the false-positive peaks. No template control and no reverse transcription control were included in all reactions.

GO and KEGG enrichment analysis

Gene Ontology (GO) enrichment analysis was used on the target gene candidates of differentially expressed miRNAs (“target gene candidates” in the following). GOseq based Wallenius non-central hyper-geometric distribution, which could adjust for gene length bias, was implemented for GO enrichment analysis. KEGG [18] is a database resource for understanding high-level functions and utilities of biological systems, such as cell, organism and ecosystem, from molecular-level information, especially large-scale molecular datasets generated by genome sequencing and other high-throughput experimental technologies (http://genome.jp/kegg/). We used KOBAS (http://kobas.cbi.pku.edu.cn/home.do) software to test the statistical enrichment of the target gene candidates in KEGG pathways. The relative correlations among miRNAs and miRNA-target genes in N.cadamba Cytoscape networks were constructed by the Cytoscape 3.6.1.

Dual-luciferase reporter assay

Based on the predicted miRNA-mRNA binding sequences, normal sequences bearing an miR166e seed binding site or with the FAS 3′-UTR site deleted were generated by two complementary chemically synthesised primers (Sangon Biotech) as follows:

wt-FAS-3’UTR-sense (ACTGACTTGGAGACACCCTGGTCTGTGAGAGTCAGTGGAG)

mut-FAS-3’UTR-sense (ACTGACTTGGAGACACAAGTTGATGTGAAGAGTCAGTGGAG)

del-FAS-3’UTR-sense (ACTGACTTGGAGACACTGTGAAGAGTCAGTGGAG)

The complementary oligonucleotides were resuspended at a 1:1 ratio (1 μg/μL each) in annealing buffer (10 mM TRIS pH 7.5–8.0, 50 mM NaCl, 1 mM EDTA) and heated at 95 °C for 10 min to denature secondary structure. The temperature was then gradually lowered to room temperature. Annealed products were cloned into the pmirGLO vector (Promega) downstream from the firefly luciferase coding region (between XhoI and XbaI sites). HeLa cells were seeded in 96-well cell culture plates (3.5 × 104 cells per well) and cultured in RPMI 1640 (Life Technologies, Grand Island, NY, USA) with 10% fetal bovine serum (FBS). The next day, cells were transfected with recombinant pmirGLO-3’UTR vector (100 ng/well) mixed with their corresponding miR166e mimics or NC (3 pmol/well, RiboBio) for 6 h using Lipofectamine 2000 (Life Technologies). Cells were harvested at 24 h after transfection, and luciferase activity was detected by a dual luciferase reporter assay system (Promega) according to the manufacturer’s recommendations. Normalised firefly luciferase activity (firefly luciferase activity/Renilla luciferase activity) for each construct was compared with that of the pmirGLO vector.

Availability of data and materials

N. cadamba leaf raw sequence reads was deposited (https://www.ncbi.nlm.nih.gov/sra/SRP144809).

Abbreviations

AMPK:

AMP-activated protein kinase

CTAB:

Cetyltrimethyl ammonium bromide

FAS:

Fatty acid synthase

GO:

Gene ontology

GSK3β:

Glycogen synthase kinase 3 beta

HCE:

Houttuynia cordata Thunb extract

KEGG:

Kyoto Encyclopedia of Genes and Genomes

LDLRAP1:

Low-density lipoprotein receptor adaptor protein 1

miRNAs:

microRNAs

N.cadamba :

Neolamarckia cadamba

qRT-PCR:

quantitative real time RT-PCR

SREBP-1c:

Sterol regulatory element binding protein-1c

TCF7:

Transcription factor 7

References

  1. Acharyya S, Dash GK, Mondal S, Dash SK. Studies on glucose lowering efficacy of the Anthocephalus cadamba (Roxb.) Miq. roots. Int J Pharm Bio Sci. 2010;1:1–9.

  2. Ahmed F, Rahman S, Ahmed N, Hossain M, Biswas A, Sarkar S, Banna H, Khatun A, Chowdhury MH, Rahmatullah M. Evaluation of Neolamarckia cadamba (Roxb.) Bosser leaf extract on glucose tolerance in glucose-induced hyperglycemic mice. Afr J Tradit Complement Altern Med. 2011;8:79–81.

    PubMed  Google Scholar 

  3. Al-Hazza A, Linley J, Aziz Q, Hunter M, Sandle G. Upregulation of basolateral small conductance potassium channels (KCNQ1/KCNE3) in ulcerative colitis. Biochem Biophys Res Commun. 2016;470:473–8.

    Article  CAS  Google Scholar 

  4. Alam MA, Akter R, Subhan N, Rahman MM, Majumder MM, Nahar L, Sarker SD. Antidiarrhoeal property of the hydroethanolic extract of the flowering tops of Anthocephalus cadamba. Rev Bras Farmacognosia-Braz J Pharmacognosy. 2008;18:155–9.

    Google Scholar 

  5. Alam MA, Subhan N, Chowdhury SA, Awal MA, Mostofa M, Rashid MA, Hasan CM, Nahar L, Sarker SD. Anthocephalus cadamba (Roxb.) Miq., Rubiaceae, extract shows hypoglycemic effect and eases oxidative stress in alloxan-induced diabetic rats. Rev Bras. 2011;21:0–0.

    CAS  Google Scholar 

  6. Alshahrani S, Alvarez-Leefmans FJ, Di Fulvio M. Expression of the Slc12a1 gene in pancreatic beta-cells: molecular characterization and in silico analysis. Cell Physiol Biochem. 2012;30:95–112.

    Article  CAS  Google Scholar 

  7. Ambros V, Bartel B, Bartel DP, Burge CB, Carrington JC, Chen X, Dreyfuss G, Eddy SR, Griffiths-Jones S, Marshall M, Matzke M, Ruvkun G, Tuschl T. A uniform system for microRNA annotation. RNA. 2003;9:277–9.

    Article  CAS  Google Scholar 

  8. Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell. 2009;136:215–33.

    Article  CAS  Google Scholar 

  9. Chandrashekar KS, Prasanna KS, Abinash B. Anti-inflammatory effect of the methanol extract from Anthocephalus cadamba stem bark in animal models. Int J Plant Biol. 2010;1:6.

    Article  Google Scholar 

  10. Chen T, Xi QY, Ye RS, Cheng X, Qi QE, Wang SB, Shu G, Wang LN, Zhu XT, Jiang QY, Zhang YL. Exploration of microRNAs in porcine milk exosomes. BMC Genomics. 2014;15:100.

    Article  Google Scholar 

  11. Chen X, Liang H, Zhang J, Zen K, Zhang CY. Secreted microRNAs: a new form of intercellular communication. Trends Cell Biol. 2012;22:125–32.

    Article  CAS  Google Scholar 

  12. Chin AR, Fong MY, Somlo G, Wu J, Swiderski P, Wu X, Wang SE. Cross-kingdom inhibition of breast cancer growth by plant miR159. Cell Res. 2016;26:217–28.

    Article  CAS  Google Scholar 

  13. Dolai N, Islam A, Haldar PK. Methanolic extract of Anthocephalus cadamba induces apoptosis in Ehrlich ascites carcinoma cells in experimental mice. Indian J Pharmacol. 2016;48:445–9.

    Article  CAS  Google Scholar 

  14. Faruque SM, Albert MJ, Mekalanos JJ. Epidemiology, genetics, and ecology of toxigenic Vibrio cholerae. Microbiol Mol Biol Rev. 1998;62:1301–14.

    Article  CAS  Google Scholar 

  15. Fatima N, Ahmad MK, Ansari JA, Khan HJ, Rastogi N, Srivastava SK, Ahmad S, Ali Z. Antiproliferative and antioxidant studies of Anthocephalus cadamba (Roxb.) Miq. Bark. Indian J Pharm Sci. 2016;78:525–31.

    CAS  Google Scholar 

  16. Friedlander MR, Mackowiak SD, Li N, Chen W, Rajewsky N. miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic Acids Res. 2012;40:37–52.

    Article  Google Scholar 

  17. Gao S, Guo Q, Qin C, Shang R, Zhang Z. Sea buckthorn fruit oil extract alleviates insulin resistance through the PI3K/Akt signaling pathway in type 2 diabetes mellitus cells and rats. J Agric Food Chem. 2017;65:1328–36.

    Article  CAS  Google Scholar 

  18. Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, Katayama T, Kawashima S, Okuda S, Tokimatsu T, Yamanishi Y. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008;36:D480–4.

    Article  CAS  Google Scholar 

  19. Kaper JB, Morris JG Jr, Levine MM. Cholera. Clin Microbiol Rev. 1995;8:48–86.

    Article  CAS  Google Scholar 

  20. Kapil A, Koul IB, Suri OP. Antihepatotoxic effects of chlorogenic acid from Anthocephalus cadamba. Phytother Res. 1995;9:189–93.

    Article  CAS  Google Scholar 

  21. Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10:R25.

    Article  Google Scholar 

  22. Lencer WI, Tsai B. The intracellular voyage of cholera toxin: going retro. Trends Biochem Sci. 2003;28:639–45.

    Article  CAS  Google Scholar 

  23. Liang H, Zen K, Zhang J, Zhang CY, Chen X. New roles for microRNAs in cross-species communication. RNA Biol. 2013;10:367–70.

    Article  CAS  Google Scholar 

  24. Liang H, Zhang S, Fu Z, Wang Y, Wang N, Liu Y, Zhao C, Wu J, Hu Y, Zhang J, Chen X, Zen K, Zhang CY. Effective detection and quantification of dietetically absorbed plant microRNAs in human plasma. J Nutr Biochem. 2015;26:505–12.

    Article  CAS  Google Scholar 

  25. Mlotshwa S, Pruss GJ, MacArthur JL, Endres MW, Davis C, Hofseth LJ, Peña MM, Vance V. A novel chemopreventive strategy based on therapeutic microRNAs produced in plants. Cell Res. 2015;25:521–4.

    Article  CAS  Google Scholar 

  26. Pandey A, Negi PS. Traditional uses, phytochemistry and pharmacological properties of Neolamarckia cadamba: a review. J Ethnopharmacol. 2016;181:118–35.

    Article  CAS  Google Scholar 

  27. Pirro S, Zanella L, Kenzo M, Montesano C, Minutolo A, Potesta M, Sobze MS, Canini A, Cirilli M, Muleo R, Colizzi V, Galgani A. MicroRNA from Moringa oleifera: identification by high throughput sequencing and their potential contribution to plant medicinal value. PLoS One. 2016;11:e0149495.

    Article  Google Scholar 

  28. Handa SS, Borris RP, Cordell GA, Phillipson JD. NMR spectral analysis of cadambine from Anthocephalus chinensis. Nat Prod. 1983;46:325.

    Article  CAS  Google Scholar 

  29. Sanadhya I, Lobo V, Bhot M, Varghese J, Chandra N. Antidiaebetic activity of Anthocephalus Indicus A. Rich. fruits in alloxan induced diabetic rats. Int J Pharm Pharm Sci. 2013;5:519–23.

    Google Scholar 

  30. Schwikkard S, van Heerden FR. Antimalarial activity of plant metabolitesElectronic supplementary information (ESI) available: IC50 values of plant metabolites. Nat Prod Rep. 2002;19:675–92 See http://www.rsc.org/suppdata/np/b0/b008980j.

    Article  CAS  Google Scholar 

  31. Umachigi SP, Kumar GS, Jayaveera KN, Kishore DVK, Ashok KCK, Dhanapal R. Antimicrobial, wound healing and antioxidant activities of Anthocephalus cadamba. Tradit Complement Altern Med. 2007;4:481–7.

    Article  Google Scholar 

  32. Viswanathan VK, Hodges K, Hecht G. Enteric infection meets intestinal function: how bacterial pathogens cause diarrhoea. Nat Rev Microbiol. 2009;7:110–9.

    Article  CAS  Google Scholar 

  33. Wang LY, Wang Y, Xu DS, Ruan KF, Feng Y, Wang S. MDG-1, a polysaccharide from Ophiopogon japonicus exerts hypoglycemic effects through the PI3K/Akt pathway in a diabetic KKAy mouse model. J Ethnopharmacol. 2012;143:347–54.

    Article  CAS  Google Scholar 

  34. Yang J, Farmer LM, Agyekum AA, Elbaz-Younes I, Hirschi KD. Detection of an abundant plant-based small RNA in healthy consumers. PLoS One. 2015;10:e0137516.

    Article  Google Scholar 

  35. Zhang L, Hou D, Chen X, Li D, Zhu L, Zhang Y, Li J, Bian Z, Liang X, Cai X, Yin Y, Wang C, Zhang T, Zhu D, Zhang D, Xu J, Chen Q, Ba Y, Liu J, Wang Q, Chen J, Wang J, Wang M, Zhang Q, Zhang J, Zen K, Zhang CY. Exogenous plant MIR168a specifically targets mammalian LDLRAP1: evidence of cross-kingdom regulation by microRNA. Cell Res. 2012;22:107–26.

    Article  CAS  Google Scholar 

  36. Zhou Z, Li X, Liu J, Dong L, Chen Q, Liu J, Kong H, Zhang Q, Qi X, Hou D, Zhang L, Zhang G, Liu Y, Zhang Y, Li J, Wang J, Chen X, Wang H, Zhang J, Chen H, Zen K, Zhang CY. Honeysuckle-encoded atypical microRNA2911 directly targets influenza a viruses. Cell Res. 2015;25:39–49.

    Article  CAS  Google Scholar 

  37. Zhu K, Liu M, Fu Z, Zhou Z, Kong Y, Liang H, Lin Z, Luo J, Zheng H, Wan P, Zhang J, Zen K, Chen J, Hu F, Zhang CY, Ren J, Chen X. Plant microRNAs in larval food regulate honeybee caste development. PLoS Genet. 2017;13:e1006946.

    Article  Google Scholar 

Download references

Acknowledgements

None.

Funding

This work was supported by grants from the National Natural Science Foundation of China [grant numbers 31472163], and The Chinese National Key Scientific Project (2016YFD0500503).

Author information

Authors and Affiliations

Authors

Contributions

JH and JS carried out the bioinformatics and analysis, JH, YR and TC participated in drafted the manuscript. QX, JL and ML collected the sample and q-PCR validation. JW, BZ and JH were performed extract total RNA. YZ conceived of the study and participated in its design and coordination and helped to draft the manuscript. All authors read and approved the manuscript.

Corresponding authors

Correspondence to Jia-Jie Sun or Yong-Liang Zhang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Supplemental Table 1.

Statistics of reads by sequencing (Expression profile of conserved and novel miRNAs in N. cadamba).

Additional file 2: Supplemental Table 2.

Reverse transcript primers for miRNAs and Highly enriched KEGG pathways for putative human targets.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

He, JJ., Ye, R., Chen, T. et al. Exploration of miRNAs in Neolamarckia cadamba and their potential cross-kingdom functions. ExRNA 2, 6 (2020). https://doi.org/10.1186/s41544-019-0043-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s41544-019-0043-8

Keywords