Chi-miR-30b-5p inhibits dermal papilla cells proliferation by targeting CaMKIIδ gene in cashmere goat
BMC Genomics volume 21, Article number: 430 (2020)
During goat embryonic morphogenesis and postnatal initiation of hair follicle (HF) regeneration, dermal papilla (DP) cells play a vital role in hair formation. Growing evidence shows that microRNAs (miRNAs) participate in HF development and DP cell proliferation. However, the molecular mechanisms have not been thoroughly investigated.
In this study, we utilized miRNA sequencing (miRNA-Seq) to identify differentially expressed miRNAs at different HF cycling stages (anagen and telogen). MiRNA-Seq has identified 411 annotated miRNAs and 130 novel miRNAs in which 29 miRNAs were up-regulated and 32 miRNAs were down-regulated in the anagen phase compared to the telogen phase. Target gene prediction and functional enrichment analysis indicated some major biological pathways related to hair cycling, such as Wnt signaling pathways, ECM-receptor interaction, VEGF signaling pathway, biosynthesis of amino acids, metabolic pathways, ribosome and oxidative phosphorylation. Also, we explored the function of chi-miR-30b-5p in regulating hair growth cycle. Similar to the HF cycling, DP cells were isolated from skin and used to investigate miRNA functions. The MTT and EdU assays showed that the viability and proliferation of DP cells were inhibited or promoted after the transfection of chi-miR-30b-5p mimic or inhibitor, respectively. Bioinformatics analysis revealed CaMKIIδ as a candidate target gene of chi-miR-30b-5p, and the dual-luciferase and western blot assay demonstrated that chi-miR-30b-5p bound to the 3’UTR of CaMKIIδ and further inhibited its translation.
Chi-miR-30b-5p was found to be highly expressed in the telogen than that in the anagen phase and could inhibit the proliferation of DP cells by targeting CaMKIIδ. Our study provides new information on the regulatory functions of miRNAs during HF development.
Cashmere is a vital raw material for the clothing industry with the characteristics of great resilience, high moisture absorption, and fine heat preservation . The Shanbei white cashmere goat is noted for its remarkable and fine fiber production features. During postnatal life, HF undergoes three periods of alternation: anagen (active growth), catagen (regression), and telogen (relative resting) . The transformation and growth of HF are regulated by various complex factors in the skin [3, 4]. Vitamin D receptor shows an essential role in anagen-catagen transformation . Some anagen inducing signals are identified by mutant mouse experiments, like Wnt/β-catenin, BMP antagonist (Noggin), and Shh [6, 7]. Furthermore, FGF5 shows an important role in the occurrence of catagen . A new hair germ generates after the telogen-to-anagen transition, and the transition between these two periods is crucial for hair regeneration. Therefore, the two stages are selected for miRNA-seq in this study.
MiRNAs are a class of small endogenous non-coding RNAs, roughly 18 ~ 25 nt in length. They can post-transcriptionally regulate gene expression by binding to the seed region (the 5′ end sequences of miRNA) and the cognate 3’UTR of mRNAs sequences to inhibit translation or induce mRNA decay [9, 10]. Many studies have indicated that miRNAs play crucial roles in skin development and regeneration [11,12,13]. It has identified and characterized miR-214, miR-21, miR-24 and miR-200 family in HF biology of human or mouse [14,15,16,17]. MiR-214 inhibits the expression of ctnnb1 (a key factor of Wnt signaling pathway) to regulate HF cycling and skin morphogenesis . MiR-24 overexpression mice display the altered HF structure and a marked defect with thinner hair cover through the repression of Tcf-3 (hair keratinocyte stemness regulator) . MiR-21 promotes cell proliferation and migration by preventing the inhibitory effects of BMP4 in primary keratinocyte cells . MiR-200 families regulate the orientation in hair germ and cell adhesion, leading to precise hair morphogenesis and cell fate specification . Previous study has shown that miR-30b-5p performs its biological function as a tumor inhibitor in human hepatocellular carcinoma . In gastric cancer, miR-30b-5p suppresses tumor growth and promotes apoptosis . Also, miR-30b-5p suppresses tumor proliferation in non-small cell lung and colorectal cancer . Yuan et al. observed that chi-miR-30b-5p expressed significantly lower in the anagen than that in the telogen of cashmere goat . However, whether chi-miR-30b-5p plays a similar function in HF remains unknown.
In order to explore the role of miRNAs in cashmere development and cycling, the differentially expressed miRNAs between the telogen and anagen were identified using miRNA-seq. Then we further examined the function of chi-miR-30b-5p in regulating hair cycle and DP cells proliferation, which would help to illustrate the regulation mechanism of miRNA on hair cycle.
Sequencing of miRNAs
Two small RNA (sRNA) libraries were constructed from the anagen and telogen phases respectively to identify the miRNAs involving cashmere HF cycle. MiRNA seq data showed the average of 13,503,395 clean reads were retained from raw data (Table 1). Generally speaking, the length interval of mammal sRNA is 18 ~ 35 nt. Therefore, the clean reads within the range were screened for subsequent analysis (Additional file 1). To evaluate the sequencing quality, the length distribution was aligned out from the four libraries (Fig. 1a). Most of the sequences were arranged between 18 ~ 25 nt, and the length of 22 nt, which is the typical length of miRNAs, had the highest proportion.
In order to analyze the expression and distribution of miRNAs, all of the filtered clean reads were mapped with the goat genome sequence. Totally 25,870,562 (50.08%) reads were aligned to the same strand in the reference sequence direction and 16,148,170 (31.26%) reads were mapped to the antisense strand (Additional file 1).
Expression analysis of annotated and novel miRNAs
First, the clean reads from above were mapped to the annotated miRNA sequence of Capra hircus in miRbase 20.0. A total of 21,154,428 reads were mapped to the annotated miRNA, which represented 50.35% of the total mapped reads (Additional file 2). As a result, 411 annotated mature miRNAs and 130 novel miRNAs were determined at the two stages of HFs cycle. At the same time, 259 annotated miRNA and 139 novel miRNA precursors were identified (Additional files 3 and 4).
The expression levels of novel and annotated miRNAs in each library were counted and normalized by TPM (Additional file 5). A total of 486 miRNAs were co-expressed at the two HF cycling stages, whereas 36 and 19 miRNAs were specifically expressed in telogen and anagen, respectively (Fig. 1b).
Based on negative binomial distribution, the differentially expressed miRNAs between the two cycling stages were analyzed, and p-adjust value < 0.05 was considered as significantly differentially expressed miRNAs. 58 annotated and 3 novel differentially expressed miRNAs were identified between the telogen and anagen. Compared with the telogen, 29 miRNAs were up-regulated and 32 miRNAs were down-regulated in the anagen (Fig. 2a, Additional file 6). The differentially expressed miRNAs in each library were clustered and visually viewed in a heat map (Fig. 2b).
MiRNA targets prediction and KEGG pathway analysis
Target genes of differentially expressed miRNAs were predicted by miRanda and TargetScan. The important genes involving the processes of HF cycle and development were identified. Figure 2c showed the regulatory relationships between miRNAs and genes.
Furthermore, the target genes were conjointly analyzed with the differentially expressed genes (DEGs) (q value ≤0.05) in the anagen and telogen in our previous study , resulting in 2347 identified target genes. KEGG (Kyoto encyclopedia of genes and genomes) pathway analysis showed that the 2347 target genes were annotated for 289 biological functions. Some traditional pathways associated with HF cycling were presented in the top 20 KEGG pathways, including ECM-receptor interaction and Wnt signaling pathway (Fig. 2d, Additional file 7). In addition, the results indicated that biosynthesis of amino acids, metabolic pathways, ribosome, and oxidative phosphorylation played an important role in HFs cycling.
Chi-miR-30b-5p suppressed the proliferation of DP cells
In order to further identify the miRNAs related with HF cycling, the differentially expressed miRNAs in our study were compared with Yuan et al’s results in the anagen and telogen . Six miRNAs (chi-miR-150, chi-miR-151-5p, chi-miR-196a, chi-miR-196b, chi-miR-30b-5p, chi-miR-9-5p) were identified (Fig. 3a). Chi-miR-30b-5p was selected for further study combining with its abundant expression from our sequencing data and its function in tumor cell proliferation. The RT-qPCR result showed that chi-miR-30b-5p was significantly lower in the anagen than that in the telogen, which was consistent with the sequencing result (Fig. 3b). To evaluate the function of chi-miR-30b-5p, its mimic and inhibitor were transfected into HEK293T cells, respectively (Fig. 4a and b). The transfection efficiency was detected by RT-qPCR, and the results revealed that the supplementation of 40, 80 and 120 nM chi-miR-30b-5p mimic all significantly increased the expression level of chi-miR-30b-5p compared with the mimic-NC group, respectively (Fig. 4a). In contrast, the supplementation of 80 and 160 nM inhibitors resulted in significant decrease of chi-miR-30b-5p expression as compared to inhibitor-NC group (Fig. 4b). Therefore, 40 nM mimic and 80 nM inhibitor were transfected into DP cells for the subsequent assays, respectively. The RT-qPCR analysis revealed that they both could significantly affect the expression of chi-miR-30b-5p in DP cells (Fig. 4c).
Further, MTT and EdU assays were performed to explore the cell proliferation after the transfection of 24 h with chi-miR-30b-5p mimic or inhibitor. Cell viability was significantly reduced in DP cells with chi-miR-30b-5p mimic compared with the control group (P < 0.01), whereas cell viability was significantly increased with the supplementation of chi-miR-30b-5p inhibitor compared to its control (P < 0.01) (Fig. 5a). Similarly, the cultured DP cells had significantly less mitotic activity (low percent of EdU positive cells) with chi-miR-30b-5p overexpression and more mitotic activity (high percent of EdU positive cells) with chi-miR-30b-5p suppression (Fig. 5b and c). Combined with the MTT result, it could be confirmed that chi-miR-30b-5p inhibited the proliferation of DP cells.
Chi-miR-30b-5p directly targets the 3’UTR of CaMKIIδ gene
To identify the direct targets of chi-miR-30b-5p in DP cells, the DEGs (p value < 0.05) associated with HF development in the anagen and telogen were conjointly analyzed with the potential targets of chi-miR-30b-5p predicted by TargetScan (Fig. 6a). Three potential targets of chi-miR-30b-5p (Fig. 6b), FZD3, CaMKIIδ, and FRS2 genes expressed significantly higher in the anagen than that in the telogen from RNA-seq profile . Dual-luciferase assay was performed to validate and explore whether chi-miR-30b-5p could bind to the 3’UTRs of FZD3, CaMKIIδ and FRS2 or not (Fig. 6c). The renilla luciferase activity was significantly reduced after the co-transfection of CaMKIIδ 3’UTR psi-check2 reporter with chi-miR-30b-5p mimic compared with the corresponding control, whereas it was significantly increased with the transfection of inhibitor (Fig. 6d). However, FZD3 and FRS2 genes had no significant changes (Fig. 6e and f). The results indicated that CaMKIIδ was a target gene of chi-miR-30b-5p other than FZD3 and FRS2.
Chi-miR-30b-5p inhibits the expression of CaMKIIδ
To confirm the effect of chi-miR-30b-5p on CaMKIIδ translation in DP cells, the CaMKIIδ protein levels were examined with the transfection of mimic and inhibitor by Western blot. The result indicated that the overexpression of chi-miR-30b-5p suppressed the expression of CaMKIIδ, whereas the inhibition of chi-miR -30b-5p promoted CaMKIIδ protein expression (Fig. 7a and b). However, the mRNA expression level of CaMKIIδ was not significantly changed with the overexpression or inhibition of chi-miR-30b-5p in DP cells (Fig. 7c), suggesting that chi-miR-30b-5p could inhibit the translation of CaMKIIδ gene. Thus, we concluded that chi-miR-30b-5p directly targeted the 3’UTR of CaMKIIδ to inhibit its translation in DP cells.
MiRNA can play an important role in developmental timing, metabolism, occurrence of disease, and cell proliferation by gene regulation in animals and plants [23,24,25,26]. High-throughput sequencing technology has been widely used in identifying the differential expressions of miRNAs in various species [27,28,29]. Cashmere is produced by the secondary HF, and its growth has obvious periodic changes: anagen (growth), catagen (cessation of growth) and telogen (resting) . One of the most important events in HF growth process is the transition from telogen to anagen phase, which called hair regrowth. It has known that miR-214, miR-21, miR-24 and miR-196a play important roles in human or mouse skin and HF development [14,15,16]. In current study, a total of 58 annotated miRNAs and 3 novel significantly differentially expressed miRNAs, including chi-miR-214-3p and chi-miR-196a, were identified between the anagen and telogen through miRNA-seq, indicating that these two miRNAs may participate in the regulation of follicular transition. The target genes of differentially expressed miRNAs were significantly enriched in Wnt, VEGF, ECM-receptor interaction, FGF and BMP Signaling pathways. The result was consistent with the studies in humans and mice hair follicle development [3, 31,32,33].
Several previous researches had reported the discovery and identification of miRNAs in goat skin by sequencing [21, 34]. Corresponding with Yuan et al’s result , chi-miR-30b-5p expressed significantly lower in the anagen than that in the telogen. Therefore, chi-miR-30b-5p was selected for further investigation. Zhu et al. found that miR-30b-5p could suppress cell proliferation and promote cell apoptosis in gastric cancer . DP cells, as an important cell model playing a vital role in hair formation [6, 35, 36], were used to investigate the function of chi-miR-30b-5p. In this study, we found that chi-miR-30b-5p inhibited DP cells proliferation through MTT and EdU assays. The cell viability and mitotic activity of DP cells were significantly reduced with the overexpression of chi-miR-30b-5p.
Primary function of miRNAs is to reduce mRNA decay or inhibit translation of their target genes . Reports show that miR-30b-5p inhibits human colorectal cancer by targeting KRAS, PIK3CD and BCL2  and modulates glioma cell proliferation by directly targeting MTDH . In order to illustrate the molecular mechanism of chi-miR-30b-5p on DP cells proliferation, the potential target genes were predicted and screened by TargetScan combining the RNA-seq data from skin tissues at the anagen and telogen of Shanbei white cashmere goat . CaMKIIδ was selected as a target gene of chi-miR-30b-5p. CaMKIIδ is one of the four isoforms of CaMKII (Ca2+/calmodulin-dependent protein kinase II), a serine/threonine protein kinase with a broad spectrum of substrates [38,39,40]. CaMKII is an important member of Wnt/Ca2+ pathway which activates Ca2+/calmodulin-dependent protein kinase II (CaMKII), and thus multiple downstream signal pathways . The dual-luciferase assay demonstrated that chi-miR-30b-5p directly targeted the 3’UTR of CaMKIIδ gene. And the overexpression of miR-30b-5p significantly decreased the CaMKIIδ expression at protein level, but not at mRNA level. Obviously, our data indicated that miR-30b-5p could inhibit the translation of CaMKIIδ gene. A previous study revealed that miR-30b-5p was down-regulated in cardiac hypertrophy, and the restoration of its function inhibited the expression of CaMKIIδ, thereby preventing cellular hypertrophy. Their findings suggested that miR-30b-5p might function as a hypertrophic suppressor . Consistent with their results, we found that chi-miR-30b-5p could inhibit cell proliferation by suppressing the expression of CaMKIIδ at protein level in DP cells.
In summary, we identified 61 differentially expressed miRNAs in anagen and telogen. Among them, we proved that chi-miR-30b-5p could inhibit DP cells proliferation by inhibiting CaMKIIδ translation (Fig. 8). Our study provides an additional insight into understanding the regulatory mechanism of hair follicle growth.
Tissues collected, RNA isolation and sequencing
Six Shanbei White cashmere goats were from Shanbei Cashmere Goats Engineering Technology Research Center (Shaanxi, China) and randomly separated into two groups. The skin tissues were collected at cashmere anagen and telogen, respectively. Acquisition of cashmere goat skin samples had been described in our previous article . After sample collection, all the animals were raised normally and recovered in two weeks with appropriate care. Total RNA was extracted from the obtained skin tissues using Trizol reagent (Invitrogen, USA). The concentration and quality of total RNA were detected using an Agilent Bioanalyzer 2100 and stored at − 80 °C. For miRNA-Seq, four miRNA pool libraries were constructed from anagen (n = 2) and telogen (n = 2). Total RNA (3 μg) from three individuals at the same stage in each pool was adopted as input material for the miRNA library. An Illumina Hiseq 2500 platform was used in the present study and 50 bp single-end reads were generated.
The acquisition of clean data (clean reads), annotated miRNA and novel miRNA prediction was defined in previous study . The clean reads with Q20, Q30 > 99% and the length between 18 ~ 35 nt were subjected to downstream analysis. The genome mapped software Bowtie (v-0.12.9) was used to alignment with the filtered sRNA clean reads .
Differential expression of miRNAs
Transcript per million (TPM): Normalization formula: Normalized expression = (actual miRNA count/total count of clean reads) *1000000 were used for assessing miRNA expression level analysis . Differential miRNA expression analysis between the anagen and telogen was executed using the DESeq2 R package. Hochberg and Benjamini method were used to standardize the p-values . Clustering analysis of differentially expressed miRNAs was achieved utilizing pheatmap R package.
MiRNA targets prediction and KEGG pathway analyses
MiRanda (http://www.microrna.org/microrna/home.do) and TargetScan (http://www.targetscan.org/vert_71/) were used to predict the potential targets of miRNAs. The miRNA-mRNA interaction network was established and visually presented using Cytoscape software (v3.6.1). Then the predicted targets were conjointly analyzed with the DEGs (q value ≤0.05) in anagen and telogen. KOBAS 3.0 (http://kobas.cbi.pku.edu.cn/index.php) was used to investigate the statistical enrichment of candidate targets in KEGG pathways.
Real-time quantitative PCR (RT-qPCR)
Extracted RNA from cells or tissues was converted to cDNA using the PrimeScript™ RT regent Kit with gDNA Eraser (Takara, Dalian, China). RT-qPCR was executed on a Bio-Rad CFX96 Touch™ Real Time PCR Detection System (Bio-Rad, USA) with TB Green™ Premix Ex TaqTM II (Takara, Dalian, China). RT-qPCR procedure was as follows: 95 °C for 1 min, then following 40 cycles of 95 °C for 10 s and appropriate annealing temperatures for 30 s. RT-qPCR primers were listed in Additional file 8 and designed by Primer Premier 6.
FZD3, CaMKIIδ gene involving Wnt signal pathway and FRS2 gene in FGF signal pathway were considered as the potential targets of chi-miR-30b-5p by TargetScan software. The FZD3–3’UTR, CaMKIIδ-3’UTR and FRS2–3’UTR sequences including chi-miR-30b-5p binding site were amplified by PCR, then the amplified fragments and psi-CHECK-2 dual-luciferase reporter vector (Promega, Madison, WI, USA) were digested by Xho I/Bsa I and Not I enzymes (NEB, New England) and subsequently ligated by T4 DNA ligase (NEB, New England). The constructed vectors were further verified by sequencing. The PCR primer information was listed in Additional file 8.
Cell culture and transfection
HEK293T cell line (GNHu17, Cell Bank of the Chinese Academy of Sciences, Shanghai, China) was cultured in high-glucose Dulbecco’s modified Eagle’s medium (DMEM; Hyclone, USA) supplemented with 1% double antibiotics (Penicillin and Streptomycin) (GM) and 10% fetal bovine serum (FBS; Hyclone, USA) in a 37 °C cell incubator with 5% CO2. The DP cells were isolated from the secondary HFs of cashmere goat and cultured in DMEM/F12 media as previous described . When the cells reached 80% confluence were transfected with chi-miR-30b-5p mimic, chi-miR-30b-5p inhibitor at an optimal concentration using Lipofectamine 2000 transfection reagent ((Invitrogen, Carlsbad, CA) and Opti-MEM (Gibco). The transfection efficiency was identified by RT-qPCR.
Cell proliferation assays
96-well plate was used to seed single cell pellets with 1 × 104 cells per well. The mimic (40 nM) and inhibitor (80 nM) of chi-miR-30b-5p were transfected into DP cells using Lipofectamine 2000. After the transfection of 24 h, the cultured medium was aspirated, then added 50 μL of MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide, 0.5 mg/mL, Sigma, St. Louis, MO) to each well, and the plates were further incubated for 4 h in a 37 °C incubator in the dark. Then, 150 μL of dimethyl sulfoxide (DMSO) (D8371, Solarbio) was added to each well to dissolve the formed formazan crystals. Cell viability was measured using optical density at 490 nm (OD490) absorbance value ratios using a SynergyH1 multi-detector microplate reader (BioTek, Winooski, VT).
EdU proliferation assay
DP cell proliferation was also evaluated using Cell-Light EdU DNA cell proliferation kit (RiboBio, Guangzhou, China). The DP cells were cultured for 24 h in a 96-well plate, then incubated with EdU medium for 2 h and the EdU positive cells were detected following the manufacturer’s protocol.
Dual-luciferase activity assay
The chi-miR-30b-5p mimic (40 nM), inhibitor (80 nM), and their corresponding negative controls were cotransfected into HEK293T cells with psiCHECK-2-FZD3–3’UTR, psiCHECK-2-CaMKIIδ-3’UTR, or psiCHECK-2-FRS2–3’UTR vector when the cells were about 80% confluence in 24-well plates by Lipofectamine 2000, respectively. The transfection reagent was changed with fresh growth medium (DMEM with 1% double antibiotics, 10% FBS) after the transfection of 6 ~ 8 h. Then the cells were rinsed with PBS one time and harvested following the manufacturer’s protocol. The luciferase activities were measured by dual luciferase reporter kit (Promega) and the Renilla Luciferase activity was normalized against Firefly Luciferase activity. The transfections were performed in triplicate.
The process of total proteins extraction and protocol of western blot were described in previous study . The antibody information was as follows: anti-CaMKIIδ (ab181052; Abcam, England), anti-GAPDH (60004–1-Ig; proteintech, USA), secondary antibodies of anti-immune rabbit IgG-HRP (Zhongshan-Bio, Beijing, China) for CaMKIIδ and anti-immune-mouse IgG-HRP (Zhongshan-Bio, Beijing, China) for GAPDH. All antibodies were diluted into 1:1000.
The fold-change of mRNA or miRNA expression was analyzed using the 2-△△CT method. β-actin or U6 was used as internal control to normalize the data. RT-qPCR data were presented as mean ± standard error of the mean (SEM). Histograms in this research were analyzed by GraphPad Prism 7.0 software (GraphPad Software, San Diego, CA, USA) with t-test. Statistical significance was considered as *P < 0.05, **P < 0.01, and ***P < 0.001. Image J was handled for gray value analysis of gels image.
Availability of data and materials
The miRNA-seq datasets supporting the conclusions of this paper are available in the NCBI SRA repository (https://www.ncbi.nlm.nih.gov/sra/). BioProject accession: PRJNA572591. BioSamples: SAMN12796315, SAMN12796316, SAMN12796317, SAMN12796318.
The lncRNA-seq datasets supporting the DEGs used in this study are available in the NBCI SRA repository. BioProject accession: PRJNA477237. BioSamples: SAMN09464004, SAMN09464005, SAMN09464006, SAMN09464007, SAMN09464008, SAMN09464009.
The miRNA-seq datasets supporting the differentially expressed miRNAs in Yuan et al’s results used the current research are available in Gene Expression Omnibus (GEO) (https://www.ncbi.nlm.nih.gov/gds/) with the accession number GSE47742.
Differentially expressed genes
Kyoto encyclopedia of genes and genomes
Real-time quantitative PCR
Transcript per million
Berger J, Buuveibaatar B, Mishra C. Globalization of the cashmere market and the decline of large mammals in Central Asia. Conservation biology : the journal of the Society for Conservation Biology. 2013;27(4):679–89.
Schneider MR, Schmidt-Ullrich R, Paus R. The hair follicle as a dynamic miniorgan. Current biology : CB. 2009;19(3):R132–42.
Tsai SY, Sennett R, Rezza A, Clavel C, Grisanti L, Zemla R, et al. Wnt/beta-catenin signaling in dermal condensates is required for hair follicle formation. Dev Biol. 2014;385(2):179–88.
Wang S, Luo Z, Zhang Y, Yuan D, Ge W, Wang X. The inconsistent regulation of HOXC13 on different keratins and the regulation mechanism on HOXC13 in cashmere goat (Capra hircus). 2018;19(1):630.
Bikle DD, Elalieh H, Chang S, Xie Z, Sundberg JP. Development and progression of alopecia in the vitamin D receptor null mouse. J Cell Physiol. 2006;207(2):340–53.
Stenn KS, Paus R. Controls of hair follicle cycling. Physiol Rev. 2001;81(1):449–94.
Paus R, Foitzik K. In search of the "hair cycle clock": a guided tour. Differentiation; research in biological diversity. 2004;72(9–10):489–511.
Hebert JM, Rosenquist T, Gotz J, Martin GR. FGF5 as a regulator of the hair growth cycle: evidence from targeted and spontaneous mutations. Cell. 1994;78(6):1017–25.
Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell. 2009;136(2):215–33.
Ambros V. MicroRNAs and developmental timing. Curr Opin Genet Dev. 2011;21(4):511–7.
Yi R, Poy MN, Stoffel M, Fuchs E. A skin microRNA promotes differentiation by repressing 'stemness'. Nature. 2008;452(7184):225–9.
Ning MS, Andl T. Control by a hair's breadth: the role of microRNAs in the skin. Cellular and molecular life sciences : CMLS. 2013;70(7):1149–69.
Botchkarev VA, Gdula MR, Mardaryev AN, Sharov AA, Fessing MY. Epigenetic regulation of gene expression in keratinocytes. The Journal of investigative dermatology. 2012;132(11):2505–21.
Ahmed MI, Alam M, Emelianov VU, Poterlowicz K, Patel A, Sharov AA, et al. MicroRNA-214 controls skin and hair follicle development by modulating the activity of the Wnt pathway. J Cell Biol. 2014;207(4):549–67.
Amelio I, Lena AM, Bonanno E, Melino G, Candi E. miR-24 affects hair follicle morphogenesis targeting Tcf-3. Cell death & disease. 2013;4:e922.
Ahmed MI, Mardaryev AN, Lewis CJ, Sharov AA, Botchkareva NV. MicroRNA-21 is an important downstream component of BMP signalling in epidermal keratinocytes. J Cell Sci. 2011;124(Pt 20):3399–404.
Hoefert JE, Bjerke GA. The microRNA-200 family coordinately regulates cell adhesion and proliferation in hair morphogenesis. 2018;217(6):2185–204.
Qin X, Chen J, Wu L, Liu Z. MiR-30b-5p acts as a tumor suppressor, repressing cell proliferation and cell cycle in human hepatocellular carcinoma. Biomedicine & pharmacotherapy = Biomedecine & pharmacotherapie. 2017;89:742–50.
Zhu ED, Li N, Li BS, Li W, Zhang WJ, Mao XH, et al. miR-30b, down-regulated in gastric cancer, promotes apoptosis and suppresses tumor growth by targeting plasminogen activator inhibitor-1. PloS one. 2014;9(8):e106049.
Liao WT, Ye YP, Zhang NJ, Li TT, Wang SY, Cui YM, et al. MicroRNA-30b functions as a tumour suppressor in human colorectal cancer by targeting KRAS, PIK3CD and BCL2. J Pathol. 2014;232(4):415–27.
Yuan C, Wang X, Geng R, He X, Qu L, Chen Y. Discovery of cashmere goat (Capra hircus) microRNAs in skin and hair follicles by Solexa sequencing. BMC Genomics. 2013;14:511.
Wang S, Ge W, Luo Z, Guo Y, Jiao B, Qu L, et al. Integrated analysis of coding genes and non-coding RNAs during hair follicle cycle of cashmere goat (Capra hircus). BMC Genomics. 2017;18(1):767.
Chen CZ, Li L, Lodish HF, Bartel DP. MicroRNAs modulate hematopoietic lineage differentiation. Science (New York, NY). 2004;303(5654):83–6.
Zhang D, Liu Z, Zheng N, Wu H, Zhang Z, Xu J. MiR-30b-5p modulates glioma cell proliferation by direct targeting MTDH. Saudi journal of biological sciences. 2018;25(5):947–52.
Liu Y, Liu R, Yang F, Cheng R, Chen X, Cui S, et al. miR-19a promotes colorectal cancer proliferation and migration by targeting TIA1. Molecular cancer. 2017;16(1):53.
Aryal B, Singh AK, Rotllan N, Price N, Fernandez-Hernando C. MicroRNAs and lipid metabolism. Curr Opin Lipidol. 2017;28(3):273–80.
Viollet C, Davis DA, Reczko M, Ziegelbauer JM, Pezzella F, Ragoussis J, et al. Next-generation sequencing analysis reveals differential expression profiles of MiRNA-mRNA target pairs in KSHV-infected cells. PLoS One. 2015;10(5):e0126439.
Ji Z, Liu Z, Chao T, Hou L, Fan R, He R, et al. Screening of miRNA profiles and construction of regulation networks in early and late lactation of dairy goat mammary glands. Sci Rep. 2017;7(1):11933.
Ivey KN, Srivastava D. microRNAs as Developmental Regulators. Cold Spring Harbor perspectives in biology. 2015;7(7):a008144.
Ibraheem M, Galbraith H, Scaife J, Ewen S. Growth of secondary hair follicles of the cashmere goat in vitro and their response to prolactin and melatonin. J Anat. 1994;185(Pt 1):135–42.
Kandyba E, Kobielak K. Wnt7b is an important intrinsic regulator of hair follicle stem cell homeostasis and hair follicle cycling. Stem cells (Dayton, Ohio). 2014;32(4):886–901.
Jacques BE, Dabdoub A, Kelley MW. Fgf signaling regulates development and transdifferentiation of hair cells and supporting cells in the basilar papilla. Hear Res. 2012;289(1–2):27–39.
Hou C, Miao Y, Wang J, Wang X, Chen CY, Hu ZQ. Collagenase IV plays an important role in regulating hair cycle by inducing VEGF, IGF-1, and TGF-β expression. Drug design, development and therapy. 2015;9:5373–83.
Liu Z, Xiao H, Li H, Zhao Y, Lai S, Yu X, et al. Identification of conserved and novel microRNAs in cashmere goat skin by deep sequencing. PLoS One. 2012;7(12):e50001.
Driskell RR, Clavel C, Rendl M, Watt FM. Hair follicle dermal papilla cells at a glance. J Cell Sci. 2011;124(Pt 8):1179–82.
Sennett R, Rendl M. Mesenchymal-epithelial interactions during hair follicle morphogenesis and cycling. Semin Cell Dev Biol. 2012;23(8):917–27.
Eulalio A, Huntzinger E, Izaurralde E. Getting to the root of miRNA-mediated gene silencing. Cell. 2008;132(1):9–14.
Braun AP, Schulman H. The multifunctional calcium/calmodulin-dependent protein kinase: from form to function. Annu Rev Physiol. 1995;57:417–45.
Means AR. Regulatory cascades involving calmodulin-dependent protein kinases. Molecular endocrinology (Baltimore, Md). 2000;14(1):4–13.
Tombes RM, Faison MO, Turbeville JM. Organization and evolution of multifunctional Ca(2+)/CaM-dependent protein kinase genes. Gene. 2003;322:17–31.
De A. Wnt/Ca2+ signaling pathway: a brief overview. Acta Biochim Biophys Sin. 2011;43(10):745–56.
He J, Jiang S, Li FL, Zhao XJ, Chu EF, Sun MN, et al. MicroRNA-30b-5p is involved in the regulation of cardiac hypertrophy by targeting CaMKIIdelta. Journal of investigative medicine : the official publication of the American Federation for Clinical Research. 2013;61(3):604–12.
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(3):R25.
Zhou L, Chen J, Li Z, Li X, Hu X, Huang Y, et al. Integrated profiling of microRNAs and mRNAs: microRNAs located on Xq27.3 associate with clear cell renal cell carcinoma. PloS one. 2010;5(12):e15224.
Storey JD. The positive false discovery rate: a Bayesian interpretation and the q -value. Ann Stat. 2003;31(6):2013–35.
Zhang X, Cheng Z, Wang L, Jiao B, Yang H, Wang X. MiR-21-3p centric regulatory network in dairy cow mammary epithelial cell proliferation. J Agric Food Chem. 2019;67(40):11137–47.
This research was supported by the National Natural Science Foundation of China (No.31772573 and No. 31972556), Key Project of Natural Science Basic Research Program of Shaanxi Province (No.2019JZ-07) and National Science and Technology Major Project of China (2018ZX08010-09B). The funding bodies played no role in the design of the study and collection, analysis, and interpretation of data, nor in writing the manuscript.
The experimental animals were fed according to the local cashmere goat standard of Shaanxi province (DB61/T583–2013, http://www.sxny.gov.cn/) with the coefficient less than 0.125. All the experimental procedures with goats used in the present study had been given prior approval by the Experimental Animal Manage Committee of Northwest A&F University (2011–31,101,684), which included written informed consent to use the animals from the owners of the animals.
Consent for publication
The authors declared no interest conflict.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
List of sRNA quantity and type after length screening (18-35 nt) and mapped to goat genome.
The information of mapped sRNA.
The information of mapped annotated miRNA.
The information of novel miRNAs.
The expression levels of annotated and novel miRNAs (normalized by TPM).
The differentially expressed miRNAs of annotated and novel miRNAs (analyzed by DESeq2).
The information of KEGG.
The information of primers.
About this article
Cite this article
Zhang, Y., Xia, S., Wang, T. et al. Chi-miR-30b-5p inhibits dermal papilla cells proliferation by targeting CaMKIIδ gene in cashmere goat. BMC Genomics 21, 430 (2020). https://doi.org/10.1186/s12864-020-06799-1