This preprint is under consideration at BMC Developmental Biology. Preprints are preliminary reports that have not undergone peer review. They should not be considered conclusive, used to inform clinical practice, or referenced by the media as validated information.
Learn more about In Review
Research article
The functions of ocu-miR-205 in regulating hair follicle development in Rex rabbits
Fuchang Li
Shandong Agricultural University
Corresponding Author
ORCiD: https://orcid.org/0000-0002-5670-8038
License:
This work is licensed under a Creative Commons Attribution 4.0 International License. Read Full License
Background: Hair follicles is an appendage from the vertebrate skin epithelium, and arise from the embryonic ectoderm andregenerate cyclically during adult life. Dermal papilla cells (DPCs) is the key dermal component of the hair follicle that directly regulates hair follicle development, growth and regeneration. Recent studies have reported that miRNA plays an important role in regulating hair follicle morphogenesis, proliferation, differentiation and apoptosis of hair follicle stem cells.
Results: The miRNAs expression profile of the DPCs from different hair density Rex rabbits shown that 240 differentially expressed of miRNAs were screened (log 2 (HD/LD)｜>1.00 and Q-value≤0.001). Among them, the expression of ocu-miR-205-5p in low hair densities DPCs was higher than that in high hair densities, and it is highly expressed in the skin tissue of Rex rabbits ( P <0.05). ocu-miR-205 could increase cell proliferation and cell apoptosis ratio, change cell cycle process ( P <0.05), affect the genes expression of PI3K/Akt, Wnt, Notch and BMP signaling pathways in DPCs and skin tissue of Rex rabbits, inhibit the protein phosphorylation level of CTNNB1, GSK-3β and the protein expression level of noggin (NOG), promote Akt phosphorylation level ( P <0.05). There was no significant change in primary follicle density ( P >0.05), but the secondary follicle density and total follicle density ( P <0.05) were changed after ocu-miR-205-5p interfered expression, and secondary/primary ratio (S/P) in ocu-miR-205-5p interfered expression group increased at 14 days after injection ( P <0.05).
Conclusion: ocu-miR-205 could promote the apoptosis of DPCs, affect PI3K/Akt, Wnt, Notch and BMP signaling pathways genes and proteins expression in DPCs and skin of Rex rabbits, promote the transformation of hair follicles from growth phase to regression and resting phase, and affect hair density of Rex rabbits.
Keywords
ocu-miR-205, Dermal papilla cell, Rex rabbit, Hair follicle density
Figure 1
Figure 2
Figure 3
Figure 4
Figure 5
Figure 6
Rex rabbits are typically used for fur production because their fur is short, fine, dense, and smooth and has important economic value [1]. The most important indicator for evaluating the fur quality of Rex rabbit is the density of hair follicles [2]. In recent years, methods to improve the hair follicle density of Rex rabbits have become the most important concern in rabbit production. Hair follicles are an appendage from the vertebrate epithelium in the skin, arise from the embryonic ectoderm, and regenerate cyclically during adulthood [3]. The process of hair follicle formation and differentiation involves at least 20 different cells and tissues [4], including the dermal papilla, hair matrix, inner root sheath, and outer root sheath, as well as different signalling pathways, such as the Wnt, Notch, bone morphogenetic protein (BMP), and fibroblast growth factor (FGF) pathways, among others [5-9]. Phosphatidylinositol 3'-kinase (PI3K) preferentially phosphorylates PIP2 to produce PIP3, and PIP3 is an important second messenger in cells that subsequently activates Akt, and thus it plays an important role in regulating the proliferation and apoptosis of hair follicle cells. The Wnt signalling pathway regulates epithelial morphogenesis, hair follicle development and cell differentiation. The PI3K/Akt signalling pathway inhibits the phosphorylation of β-catenin by phosphorylating glycogen synthase kinase 3β (GSK-3β) and activates the Wnt signalling pathway[5]. Dickkopf-related protein 1 (DKK1) inhibits the Wnt signalling pathway by inhibiting the phosphorylation of β-catenin and induces hair follicle regression [10]. When Notch receptor binds to a ligand, it activates hair follicle stem cells and then promotes the transition of hair follicles from the resting stage to growing stage [6]. The BMP signalling pathway is involved in embryonic skin appendage organ morphogenesis and postnatal hair follicle growth [7]. The BMP2 and BMP4 genes inhibit hair follicle development and are associated with maintaining hair follicles in the resting stage[8]. Noggin (NOG) is an inhibitor of the BMP signalling pathway, and its abnormal expression leads to follicular enlargement [7]. The Notch signalling pathway interacts with the BMP signalling pathway, and the BMP signalling pathway inhibits the Wnt signalling pathway by regulating β-catenin expression [6]. Dermal papilla cells (DPCs) is the key dermal component of the hair follicle, and it can directly regulate hair follicle development, growth and regeneration [11]. Besides, the characteristic of DPCs also affects the size, shape and cycling of hair follicle [12]. Signal exchange between DPCs and hair follicle stem cells at telophase is the key to initiating the next hair follicle cycle [13].
MicroRNAs (miRNAs) are small non-coding RNAs that have been plays an important role in embryogenesis, organ development, cell proliferation, metabolism, embryogenesis, behaviour and other biological processes [14-19]. Recently, miRNAs were reported to play important roles in regulating hair follicle morphogenesis and the proliferation, differentiation and apoptosis of hair follicle stem cells in mice, rats, goats and sheep [20-22]. Notably, miR-214 inhibits hair follicle growth and development by modulating the expression of regulatory factors in Wnt signalling pathway, such as β-catenin and lymphoid enhancer-binding factor 1 (Lef-1) [23]. Upon the overexpression of DKK1 in transgenic mice, the expression of miR-200b and miR-196a in epidermis decreases significantly, which is possibly mediated by potential target genes acting on the Wnt signalling pathway [24]. As shown in a previous study, miR-let-7b promotes alpaca hair growth by inhibiting the transcription of transforming growth factor β receptor 1 (TGFβR1) [25]. BMP4 negatively regulates the expression of miR-21, and miR-21 negatively regulates the expression of the BMP-dependent tumour suppressor genes Pten, Pdcd4, Timp3 and Tpm1 [26]. Hair follicle development and hair growth in mice are regulated by miR-31 through effects on the BMP and Wnt signalling pathways [20]. Additionally, miR-205 is a highly conserved miRNA that shares a similar expression pattern with miR-200 family [27]. It is one of the miRNAs expressed at the highest levels in the epidermis [28, 29], and it plays an essential role in promoting the neonatal expansion of skin stem cells during early development by modulating the PI3K pathway [30].
We isolated the DPCs from the skin of Rex rabbits and analysed the miRNA expression profiles of the DPCs from Rex rabbits with different hair densities to improve the fur quality of Rex rabbits. Among the miRNAs analysed, ocu-miR-205 was one of the miRNAs expressed at the highest levels in DPCs. We examined the function of ocu-miR-205 in hair follicle development.
DPCs show a complex miRNA expression pattern
The varieties of miRNAs in DPCs from 30-day-old Rex rabbits with lower and higher hair densities (Additional file 1) were studied by subjecting RNA samples with a high integrity and qualified quality (Additional file 2) to high-throughput small RNA sequencing using the BGISEQ-500 platform. 37930744, 39442011, 40965907, 38502653, 40622117 and 41149163 clean reads were obtained from the six samples (Additional file 3), and the majority of clean reads had a length of 23 nucleotides (Additional file 4). Comparison with known small RNA databases, the percentage of matching reads for each library was 91.91 %, 91.05 %, 92.46 %, 92.77 %, 90.13 %, and 91.61 %, respectively (Additional file 5). The results of the small RNA classification showed that miRNAs accounted for 80.80 %, 82.50 %, 81.60 %, 85.80 %, 76.50 %, and 81.80 %, respectively (Additional file 6). The base Q20 value of filtered data was greater than 90 % and the Q30 value was greater than 80 %. Additional files 7 and 8 show the quantity and quality distribution maps of the bases in each sample, respectively.
By screening differentially expressed genes (DEGs), 240 differentially expressed miRNAs were identified (｜log2(HD/LD)｜>1.00 and Q-value≤0.001; Fig. 1a and Additional file 9), including 122 miRNAs that were upregulated and 118 miRNAs that were downregulated (Fig. 1b). The annotation of the target genes of differentially expressed miRNAs was performed using Gene Ontology (GO) enrichment, and a total of 205,661 target genes were enriched in GO terms (Fig. 1c). Specific GO terms of the target genes mainly involved in the biological process (BP), cellular component (CC) and molecular function (MF) categories. A directed acyclic graph (DAG) of the enriched GO terms in the BP, CC and MF categories is provided in Additional file 10. Following, the target genes were also uploaded into the Kyoto Encyclopaedia of Genes and Genomes (KEGG) database to identify the pathways that were actively regulated by miRNAs in DPCs. Three hundred twenty-five pathways were predicted, and the top enriched pathways were the Wnt signalling pathway and Notch signalling pathway, among others (Fig. 1d).
The expression of ocu-miR-205
The structure of one of the differentially expressed miRNAs, ocu-miR-205, is shown in Fig. 2a. The expression of ocu-miR-205 differed in DPCs from Rex rabbits with different hair densities (｜log2(HD/LD)｜>1.00 and Q-value≤0.001; Fig. 2b) and was expressed at low levels in DPCs from Rex rabbits with a high hair density (HD) and at high levels in DPCs from Rex rabbits with a low hair density (LD). Quantitative fluorescence PCR results also confirmed the accuracy of the ocu-miR-205-5p sequencing results (Fig. 2c). Furthermore, the expression of ocu-miR-205-5p in different tissues from Rex rabbits differed (P<0.05). The expression in the skin tissue was higher than in the stomach, intestine, spleen, liver, heart, lung, kidney, muscle and fat (P<0.05; Fig. 2d), suggesting a tissue-specific expression pattern.
Effects of ocu-miR-205 on DPCs from Rex rabbits
HBAD-GFP, HBAD-ocu-miR-205-GFP, and HBAD-ocu-miR-205-5p-sponge-GFP were constructed and used to transfect DPCs at a multiplicity of infection (MOI) of 200. Cell proliferation, the cell cycle, cell apoptosis, and changes in the expression of ocu-miR-205-5p and hair follicle development-related genes and proteins were detected 48 h after transfection. The constructed adenoviruses successfully transfected DPCs (Fig. 3a), and ocu-miR-205-5p expression increased significantly after overexpression, but decreased significantly after silencing (P<0.05; Fig. 3b). Importantly, ocu-miR-205 increased cell proliferation and the cell apoptosis rate, in addition to altering the progression of the cell cycle (P<0.05; Table 1). Moreover, ocu-miR-205 inhibited the expression of the Inppl1, Frk and Phlda3 mRNAs, which are involved in the PI3K/Akt signalling pathway (P<0.05). The expression of ocu-miR-205 induced the expression of the DKK1 mRNA and inhibited the expression of the Wnt10b, CTNNB1 and GSK-3β genes in the Wnt signalling pathway (P<0.05). The expression of the Notch1, Jagged1, Hes1 and Hes5 genes in the Notch signalling pathway was suppressed by ocu-miR-205 (P<0.05), whereas the expression of the BMP2, BMP4 and TGF-β1 genes in the BMP signalling pathway was induced (P<0.05; Table 2). The expression of ocu-miR-205 inhibited the phosphorylation of the CTNNB1 and GSK-3β proteins, decreased the level of the noggin (NOG) protein, and increased the level of Akt phosphorylation (P<0.05; Fig. 4).
Effects of ocu-miR-205 on the skin tissue of Rex rabbits
One hundred 3-month-old Rex rabbits with similar body weights were randomly divided into 4 groups. After local shaving of the back skin, 50 μL of HBAD-GFP, HBAD-ocu-miR-205-GFP and ocu-miR-205-5p-sponge-GFP were injected intradermally into the rabbits in each group. 24 h after the injected, one randomly selected Rex rabbit from each group was euthanized, and frozen sections were prepared from the locally injected skin to confirm the expression of the miRNAs from the injected adenoviruses. Eight Rex rabbits were randomly selected from each group 7 days, 14 days and 21 days after transfection. After euthanasia, the injected skin was collected and divided into two parts. A part placed in freezing tube and the other part was fixed at a 4 % paraformaldehyde fixative solution. The frozen samples were used to detect the changes in the expression of ocu-miR-205-5p and hair follicle development-related genes and proteins. The fixed samples were used to prepare paraffin sections, and the hair density was counted after HE staining. The constructed adenoviruses successfully transfected the skin and hair follicles of Rex rabbits (Fig. 5a), and the expression of ocu-miR-205-5p increased significantly after overexpression, but decreased significantly after silencing (P<0.05; Fig. 5b). Fourteen days after the injection, ocu-miR-205 significantly altered the expression of genes involved in the PI3K/Akt, Wnt, Notch and BMP signalling pathways (P<0.05; Table 3). Moreover, ocu-miR-205 significantly altered the levels of the phosphorylated CTNNB1, GSK-3β, and Akt proteins and the level of the NOG protein (P<0.05; Fig. 6). A significant change in the primary follicle density was not observed (P>0.05), but the secondary follicle density and total follicle density (P<0.05) were altered after ocu-miR-205-5p interfered expression, and the secondary/primary ratio (S/P) in the ocu-miR-205-5p interfered expression group increased at 14 days after the injection (P<0.05; Table 4).
Rex rabbits mainly used for fur production, and the fur is the most important livestock products. The hair follicle cycle can typically be divided into anagen, catagen, and telogen phases[31]. The hair follicle of Rex rabbits contains primary and secondary hair follicles, and the different sizes and types of primary and secondary hair follicles is easy distinguished[32]. Besides, the hair follicle growth cycle is easy distinguish the differentiation of hair cycle phases [33]. Many factors affect the growth and development of hair follicles in Rex rabbits, such as genetic, nutritional, age, among others [34-37]. The most hair follicles in anagen phases when the Rex rabbits at 4-5 weeks[37]. Notably, miRNAs plays an important functions in regulating gene expressions[38-40]. Among them, miR-205 is abundantly expressed in the epidermis[41,42], and it plays an essential role by modulating the PI3K/Akt signalling pathway in promoting the neonatal expansion of skin stem cells during early development[43]. In vitro models, miR-205 promotes keratinocyte migration by targeting the lipid phosphatase SHIP2 and KIR4.1[44, 45]. There is a potential therapeutic application of DPCs in the treatment of alopecia, because DPCs of various origins induce the de novo formation of the hair follicle structure in both follicular and afollicular epidermis [46-49]. In the present study, ocu-miR-205-5p was expressed at significantly higher levels in DPCs from rabbits with low hair densities than in DPCs from rabbits with high hair densities. Moreover, ocu-miR-205 induced G0/G1 arrest in DPCs, which was further confirmed by the reduction in the population of DPCs in the G0/G1 phase and increase in the apoptotic rate after the transfection of the ocu-miR-205 inhibitor. This finding is consistent with previous studies [50]. Furthermore, overexpression of miR-205 significantly inhibited cell proliferation. Additionally, ocu-miR-205 inhibited the expression of related gene and proteins in the PI3K/Akt, Wnt, and Notch signalling pathways, and activated the BMP signalling pathway. Therefore, ocu-miR-205 plays an important role in regulating hair follicle development.
In conclusion, ocu-miR-205 promoted the apoptosis of DPCs, inhibited cell proliferation, modulated the expression of genes and proteins involved in the PI3K/Akt, Wnt, Notch and BMP signalling pathways in DPCs and the skin of Rex rabbits, promoted the transition of hair follicles from the growth phase to the regression and resting phase, and altered the hair density of Rex rabbits.
Animal and sample collection
The Rex rabbits used in this study were purchased from Taishan Rabbit Farm (Shandong, China). Ten adult female Rex rabbits with a high (>14000/cm2, HD) or low hair density (<10000/cm2, LD) were chosen and divided into 2 groups. An adult male rabbit was selected to mate with the female rabbits within each group. The first generation of offspring (F1 generation) was obtained and then inbred with F1 in the same nest (HD and HD, LD and LD inbreeding). Three rabbits each with high and low hair densities were selected from the second generation of offspring (F2 generation) 30 days after birth. Selected animals were electrically stunned (120 V, pulsed direct current, 50 Hz for 5 s) and euthanized by exsanguination of the carotid artery and before skin was harvested from the experimental rabbits. The DPCs were separated and cultured, methods according to references[51]. After the cells had established a monolayer (approximately 12 days), the total RNA was extracted, tested for quality.
Construction of the small RNA libraries, sequencing analysis, miRNA identification and prediction of new miRNAs in DPCs
The methods of construction small RNA libraries according to references[52]. Solexa sequencing by synthesis method using the BGISEQ-500 platform at Shenzhen Huada Biotech Co., Ltd. (Shenzhen, China) [53,54]. The raw reads produced from sequencing were filtered to remove low-quality reads, and the clean reads were analysed using BLAST with Bowtie-1.0.0 software, Rfam [55]. The clean reads were used for miRNA identification and compared with the mature miRNAs and pre-miRNAs from Oryctolagus cuniculus listed in miRBase 21.0 [56]. Subsequently, the miRDeep 2 software was used to predict the novel miRNAs by exploring the secondary structure[57], minimum free energy and dicer cleavage sites of the unannotated clean reads that mapped to the Oryctolagus cuniculus genome. After the identification of conserved miRNAs, the clean reads were aligned to the Oryctolagus cuniculus genome for identify new Oryctolagus cuniculus miRNAs.
In order to eliminate the effect of different sequencing quantities on quantitative accuracy, transcripts per million (TPMs) were calculated to standardize the expression levels of small RNA [58]. Based on the assumption that RNA sequencing is a random process, the expression of each transcript is presumed to exhibit a binomial distribution[59]. DEGseq calculated differential expression based on MA-plot [60, 61], and the P-value of each miRNA was corrected by performing multiple hypothesis tests using Q-values. When the difference in coincidence was more than two-fold and the Q-value was less than or equal to 0.001, the miRNAs were considered significantly differentially expressed (｜log2(HD/LD)｜>1.00 and Q-value≤0.001).
The miRNAs were aligned to the EST unigenes of Oryctolagus cuniculus and the target genes were predicted using the miRanda algorithm to obtain a better understanding of the potential functions of the significantly differentially expressed miRNAs in Rex rabbits with different hair densities [62]. An enrichment analysis of the predicted target genes was conducted with GO terms and KEGG pathways [63].
Construction and identification of adenovirus vectors overexpressing and silencing ocu-miR-205
HBAD-GFP (HANBIO adenovirus-green fluorescent protein; empty vector), HBAD-ocu-miR-205-GFP (overexpression), HBAD-ocu-miR-205-5p-sponge-GFP (silencing) adenoviruses were synthesized and constructed by Hanheng Biotechnology Co., Ltd. (Shanghai, China). The infective titres of HBAD-GFP, HBAD-ocu-miR-205-GFP and HBAD-ocu-miR-205-5p-sponge-GFP were 1.26*1010 PFU/mL, 1.58*1010 PFU/mL and 1.26*1010 PFU/mL, respectively.
Third-generation DPCs displaying good growth conditions were inoculated into a disposable 6-well plate. The cell density was approximately 1.0*105 cells/mL. Prior to the infection, the virus was subjected to 10-fold gradient dilution. Generally, the MOI (multiplicity of infection) was controlled in the range of 10-1000. HBAD-GFP, HBAD-ocu-miR-205-GFP and HBAD-ocu-miR-205-GFP were individually transfected into Rex rabbit DPCs at an MOI 200, and a negative control was established using cells undergoing normal culture.
Fifty microliters of the purified adenovirus were injected into the skin of each Rex rabbit with microinjector at a concentration of 5.0*108 - 1.0 * 109 virus particles per Rex rabbit after shaving the middle part of the back of 100 3-month-old Rex rabbits with similar body weights and good health. 24 h after transfection, one Rex rabbit from each group was randomly selected, euthanized, and frozen sections were prepared from the locally injected skin. The adenovirus-transfected skin was observed under a positive fluorescence microscope (Nikon ECLIPSE 80i, Japan).
Assessment of the proliferation, cell cycle and apoptosis of DPCs
DPCs were plated in a 96-well plate at a density of 104 cells/well, cultured in basal medium for 24 h, and then transfected with the indicated adenoviruses. The methods assessment of the proliferation according to references[51]. Third-generation DPCs were plated in a disposable 6-well plate at a density of 104 cells/mL, with 2 mL of the cell suspension plated in each well. After a 24-h incubation to allow cells to adhere, the culture medium was removed. After treatment and culture for a certain time, the cells were digested with a trypsin digestion solution lacking EDTA (Solarbio, Beijing, China). The methods assessment of the cell cycle according to references[51]. For apoptosis assessment, cells were centrifuged and collected into a 1.5 mL centrifugal tube, and then washed with PBS. After centrifugation, 500 μL of 10X Annexin V Binding Buffer was added to re-suspended the cells, followed by the labelling of F-actin. Cells were incubated with the FITC Annexin V and Propidium Iodide Staining Solution for 15 min at 4 °C and then analysed using flow cytometry. The percentages of early apoptotic cells (Q4), late apoptotic cells (Q2) and total apoptotic cells (Q2 + Q4) in each sample were calculated.
Total RNA extraction and Real-time PCR analysis
The methods of total RNA extraction and real-time PCR of mRNA was according to references[64]. All quantitative PCR primers were designed using Primer Premier 5 software (Additional file 11). For the quantitative RT-PCR of miRNAs, 1 μg of total RNA was reverse transcribed with Bulge-Loop miRNA-specific reverse transcription primers (RiboBio, China), and quantitative PCR was performed using Fast Start Universal SYBR Green Master Mix (Roche Diagnostics GmbH Mannheim, Germany) and Bulge-Loop primers (RiboBio, Guangzhou, China) on the 7500 Fast System 1.4 system with small nuclear RNA U6 as the normalisation control. The volume of each reaction was 20 µL, including 2 µL of cDNAs, 10 µL of SYBR Green Master (2X), 0.8 μL of the Bulge-LoopTM miRNA Forward Primer (5 µM), 0.8 µL of the Bulge-LoopTM Reverse Primer (5 µM), 0.4 µL of ROX Reference Dye II (50×) and 6.0 µL of ddH2O. PCR was performed under the following conditions: 10 min of template denaturation at 95 °C, followed by 40 cycles of 95 °C for 2 s, 60 °C for 20 s, and 70 °C for 10 s. Melting curves (70 °C-95 °C) for each sample were analysed after each run to confirm the specificity of amplification reactions. The relative expression levels of mRNAs and miRNAs were calculated using the arithmetic formula 2−△△Ct [65].
Western immunoblotting.
The methods of total protein extraction and SDS-PAGE were according to references[64]. The membranes were blocked with 5 % skimmed milk in phosphate-buffered saline (PBS; Solarbio, China) at 4 °C overnight and incubated with primary antibodies (tubulin AT819, Beyotime, China; phospho-CTNNB1-S552 pAb, Abcam, US; phospho-GSK3B-S9 pAb, Abcam, US; phospho-AKT1-S473 pAb, Abcam, US; or NOG polyclonal antibody, Abcam, US). The membranes were then rinsed with Tris-buffered saline containing Tween (TBST; Solarbio, China), and subjected to detection with a 1:3000 dilution of a horseradish peroxidase (HRP)-conjugated goat anti-mouse IgG antibody (Beyotime, China) at 37 °C for 1 h. Proteins were visualized using BeyoECL reagents (Beyotime, China). The intensity of the bands was quantified with a Pro Plus 6.0 Biological Image Analysis System. The levels of phospho-CTNNB1, phospho-GSK3B, phospho-AKT1 and NOG were normalized to the internal control beta-tubulin, and the relative expression levels were calculated.
Statistical analysis
All data were analysed with SAS software (SAS version 8e; SAS Institute, Cary, NC, USA). A one-way ANOVA was used to evaluate the differences in mean values among various groups. The data are presented as the means and R-MSE. P < 0.05 was regarded as statistically significant.
BMP: bone morphogenetic protein; BP: biological process; CC: cellular component; DAG: Directed Acyclic Graph; DEGs: differentially expressed genes; DPCs: dermal papilla cells; DKK1: Dickkopf-related protein 1; FGF: fibroblast growth factor; GAPDH: glyceraldehyde-3-phosphate dehydrogenase; GFP: green fluorescent protein; GO: Gene Ontology; GSK-3β: glycogen synthase kinase 3β; HD: high hair density; HE: haematoxylin and eosin; Lef-1: lymphoid enhancer-binding factor 1; LD: low hair density; miRNAs: microRNAs; mRNA: messenger RNA; NOG: noggin; S/P: secondary/primary ratio; SDS: sodium dodecyl sulfate; TGFβ: transforming growth factor β; TPM: transcripts per million; ocu-miR-205: Oryctolagus cuniculus microRNA 205; PAGE: polyacrylamide gel electrophoresis; PBS: phosphate-buffered saline; PI3K: phosphatidylinositol 3'-kinase.
Ethics approval and consent to participate
The experimental procedures were approved by the Committee of Ethics in Research of Shandong Agricultural University (SDAUA-2017-029) and performed in accordance with the Guidelines for Experimental Animals of the Ministry of Science and Technology (Beijing, China). The experimental procedures in this test compliant with the ARRIVE guidelines, in accordance with the U.K. Animals (Scientific Procedures) Act, 1986 and associated guidelines, EU Directive 2010/63/EU for animal experiments, or with the National Institutes of Health guide for the care and use of Laboratory animals (NIH Publications No. 8023, revised 1978). This article does not any studies with human participants performed by any of the authors.
Consent for publication
Not applicable.
Availability of data and materials
Small RNA sequencing data has been deposited in Sequence Read Archive of the National Center for Biotechnology Information under temporary submission ID: SUB7193686, and BioSample accessions: SAMN14476334, SAMN14476335, SAMN14476336, SAMN14476337, SAMN14476338, SAMN14476339.
Competing interests
The authors declare that they have no competing interests.
Funding
This study was supported by the Earmarked Fund for Modern Agro-industry Technology Research System (CARS-43-B-1), National Natural Science Foundation of China (31972594), and Funds of Shandong “Double Tops” Programme (SYL2017YSTD11). The funding bodies had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Authors’ contributions
GL, SL and LB isolation, culture and identification of DPCs, and construction of the small RNA libraries, sequencing analysis, miRNAs identification and prediction of new miRNAs. GL and HS isolated and cultured hair follicles, and performed cellular experiments. GL, HL and SG raised the rabbits and conducted the animal experiments. GL, WJ and FL revised the manuscript. YZ, WJ and FL designed this work and wrote this manuscript. All authors read and approved the final manuscript.
Acknowledgements
Not applicable.
Author details
1 College of Animal Science and Technology, Shandong Agricultural University, Tai’an 271018, P. R. China; 2 Shandong Provincial Key Laboratory of Animal Biotechnology and Disease Control and Prevention, Tai’an 271018, P. R. China; 3 Animal Husbandry and Veterinary Institute, Shandong Academy of Agricultural Sciences, Jinan 251000, P. R. China；4 Shandong Key Laboratory of Animal Disease Control and Breeding, Jinan 251000, P. R. China.
- Ruben D. Choosing a Rex rabbit. Small mammal breeds. 2014.
- Gu Z, Ren W, Huang R, Huang Y, Chen B. Study on density of rex rabbit. Chinese Journal of Rabbit Farming.1999; 4:18-21.
- Millar SE. Molecular mechanisms regulating hair follicle development. Journal of Investigative Dermatology. 2002; 118:216–225.
- Fuchs E. Scratching the surface of skin development. Nature. 2007;445:834-842.
- Lin CM, Yuan Y, Chen X, Li H, Cai B, Liu Y, et al. Expression of Wnt/β–catenin signaling, stem-cell markers and proliferating cell markers in rat whisker hair follicles. Journal of Molecular Histology. 2015;46(3):233-240.
- Demehri S, Kopan R. Notch signaling in bulge stem cells is not required for selection of hair follicle fate. Development. 2009; 136(6): 891-896.
- Kulessa H, Turk G, Hogan BL. Inhibition of Bmp signaling affects growth and differentiation in the anagen hair follicle. Embo Journal. 2000;19(24):6664-6674.
- Su R, Li J, Zhang W, Yin J, Zhao J, Chang Z. Expression of BMP2 in the skin and hair follicle from different stage in Inner Mongolia cashmere goat. Scientia Agricultura Sinica. 2008; 23:559-563.
- Foitzik K, Lindner G, Mueller-Roever S, Maurer M, Botchkareva N, Botchkarev V, et al. Control of murine hair follicle regression (catagen) by TGF-beta1 in vivo. The FASEB Journal. 2000;14:752-760.
- Greco V, Chen T, Rendl M, Schober M, Pasolli HA, Stokes N, et al. Two-step mechanism for stem cell activation during hair regeneration. Cell Stem Cell. 2009; 4(2):155–169.
- Driskell RR, Clavel C, Rendl M, Watt FM. Hair follicle dermal papilla cells at a glance. J. Cell Sci. 2011;124:1179-1182.
- Chi W, Wu E, Morgan BA. Dermal papilla cell number specifies hair size, shape and cycling and its reduction causes follicular decline. Development. 2013; 140:1676-1683.
- Woo C, Eleanor W, Bruce A. Morgan permal papilla cell number specifies hair size, shape and cycling and its reduction causes follicular decline. Development and stem cells. 2013; 140: 1676-1683.
- Yi R, Fuchs E. MicroRNAs and their roles in mammalian stem cells. J. Cell Sci. 2011; 124: 1775–1783.
- Carrington JC, Ambros V. Role of microRNAs in plant and animal development. Science 2003; 301:336–338.
- Xu P, Vernooy SY, Guo M, Hay BA. The Drosophila microRNA Mir-14 suppresses cell death and is required for normal fat metabolism. Curr. Biol. 2003; 13: 790–795.
- Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell 2004; 116:281–297.
- Teleman AA., Maitra S, Cohen SM. Drosophila lacking microRNA miR-278 are defective in energy homeostasis. Genes Dev. 2006; 20: 417–422.
- He L, Lim LP, Stanchina E, Xuan Z, Liang Y, Xue W, et al. A microRNA component of the p53 tumour suppressor network. Nature. 2007; 447:1130–1134.
- Mardaryev AN, Ahmed MI, Vlahov NV, Fessing MY, Gill J H, Sharov AA, et al. Micro-RNA-31 controls hair cycle-associated changes in gene expression programs of the skin and hair follicle. FASEB Journal. 2010; 24 (10), 3869 -3881.
- 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): 50001.
- Wang P, Hong W, Zhen J, Ren J, Li Z, Xu A. The changes of microRNA expression profles and tyrosinase related proteins in MITF knocked down melanocytes. Molecular BioSystems. 2012; 8: 2924-2931.
- Ahmed MI, Alam M, Emelianov VU. MicroRNA214 controls skin and hair follicle development by modulating the activity of the Wnt pathway. Journal of Cell Biology. 2014; 207(4), 549-567.
- Andl T, Murchison EP, Liu F, Zhang Y, Yunta-Gonzalez M, Tobias JW, et al. The miRNA-processing enzyme dicer is essential for the morphogenesis and maintenance of hair follicles. Current Biology. 2006;16(10), 1041–1049.
- Shen Y, Zhang Y, Liu N, Wang H, Xie J, Gao S, et al. Let-7b promotes alpaca hair growth via transcriptional repression of TGFβRI. Gene. 2016; 577:32–36.
- Ahmed MI, Mardaryev AN, Lewis CJ, Sharov AA, Botchkareva NV. MicroRNA-21 is an important downstream component of BMP signaling in epidermal keratinocytes. Journal of Cell Science 2011; 124(20), 3399–3404.
- Gregory PA, Bert AG, Paterson EL, Barry SC, Tsykin A, Farshid G, et al. The miR-200 family and miR-205 regulate epithelial to mesenchymal transition by targeting ZEB1 and SIP1. Nat. Cell Biol. 2008; 10, 593–601.
- Yi R, O\"Carroll D, Hilda AZ, Zhihong D, Fred ST, Alexander F, et al. Morphogenesis in skin is governed by discrete sets of differentially expressed microRNAs. Nat. Genet. 2006; 38, 356–362.
- Ryan DG, Oliveira-Fernandes M, Lavker RM. MicroRNAs of the mammalian eye display distinct and overlapping tissue specificity. Mol. Vis. 2006; 12, 1175–1184.
- Wang D, Zhang Z, O'Loughlin E, Wang L, Fan X, Lai E. et al. MicroRNA-205 controls neonatal expansion of skin stem cells by modulating the PI(3)K pathway. Nat. Cell Biol. 2013;15, 1153–1163.
- 3 Alonso L, Fuchs E. The hair cycle. J Cell Sci 2006;119,391-3.
- 3 Galbraith H. Fundamental hair follicle biology and fine fibre production in animals. Animal 2010;4,1490-509.
- 3 Zhu B, Xu T, Yuan J, Guo X, Liu D. Transcriptome sequencing reveals differences between primary and secondary hair follicle derived dermal papilla cells of the Cashmere goat (Capra hircus). PloS one 2013;8, e76282.
- Wang ZP, Zhang H, Yang H, Wang SZ, Rong EG, Pei WY, Li H, Wang N. Genome-wide association study for wool production traits in a Chinese Merino sheep population. Plos One 2014;9, e107101.
- Liu NH, Li K, Liu J, Yu M, Cheng W, De J, Liu S, Shi Y, Zhao JS. Differential expression of genes and proteins associated with wool follicle cycling. Mol Biol Rep 2014; 41:5343-5349.
- Meale SJ, Chaves AV, Ding S, Bush RD, McAllister TA. Effects of crude glycerin supplementation on wool production, feeding behavior, and body condition of Merino ewes. J Anim Sci 2013; 91, 878-885.
- Wu Z, Sun L, Liu G, Liu H, Liu H, Yu Z, et al. Hair follicle development and related gene and protein expression of skins in Rex rabbits during the first 8 weeks of life. Asian-Australas J Anim Sci. 2019; 32(4):477-484.
- Wu M, Sun Q, Guo X, Liu H. hMSCs possess the potential to differentiate intoc DP cells in vivo and in vitro. Cell Biol Int Rep. 2010; 19 (2), e00019.
- Lim LP, Glasner ME, Yekta S, Burge CB, Bartel DP. Vertebrate microRNA genes. Science. 2003; 299, 1540.
- Wienholds E, Plasterk RH. MicroRNA function in animal development. FEBS Lett. 2005; 579, 5911–5922.
- Yi R., O\"Carroll D, Hilda AZ, Zhihong D, Fred ST, Alexander F, et al. Morphogenesis in skin is governed by discrete sets of differentially expressed microRNAs. Nat. Genet. 2006; 38:356–362.
- Ryan DG, Oliveira-Fernandes M, Lavker RM. MicroRNAs of the mammalian eye display distinct and overlapping tissue specificity. Mol. Vis. 2006; 12: 1175–1184.
- Wang D, Li Q, Feng N, Cheng G, Guan Z, Wang Y, et al. MicroRNA-205 controls neonatal expansion of skin stem cells by modulating the PI3K pathway. Nat. Cell Biol., 2013; 15:1153–1163.
- Yu J, Peng H, Ruan Q, Fatima A, Getsios S, Lavker RM, et al. MicroRNA-205 promotes keratinocyte migration via the lipid phosphatase SHIP2. FASEB J. 2010; 24:3950–3959.
- Lin D, Halilovic A, Yue P, Bellner L, Wang K, Wang L, et al. Inhibition of miR-205 impairs the wound-healing process in human corneal epithelial cells by targeting KIR4.1 (KCNJ10). Invest. Ophthalmol. Vis. Sci. 2013; 54:6167–6178.
- Du T, Zamore PD. Beginning to understand microRNA function. Cell Res. 2007; 17:661–663.
- Kollar EJ. The induction of hair follicles by embryonic dermal papillae. J. Invest. Dermatol. 1970; 55:374-378.
- Oliver RF. The induction of hair follicle formation in the adult hooded rat by vibrissa dermal papillae. J. Embryol. Exp. Morphol. 1970; 23:219-236.
- Jahoda CA, Horne KA, Oliver RF. Induction of hair growth by implantation of cultured dermal papilla cells. Nature1984; 311:560-562.
- Zhang J, He X, Tong W, Johnson T, Wiedemann L, Mishina Y, et al. Bone morphogenetic protein signaling inhibits hair follicle anagen induction by restricting epithelial stem/progenitor cell activation and expansion. Stem Cells 2006; 24:2826-2839.
- Liu G, Bai L, Li S, Liu H, Zhu Y, Sun H, et al. Isolation, culture and growth characteristics of dermal papilla cells from Rex rabbits. Tissue and Cell. 2020; 65,101348.
- Zha W, Zhou L, Li S, Liu K, Yang G, Chen Z, et al. Characterization and comparative profiling of the small RNA transcriptomes in the Hemipteran insect Nilaparvata lugens. Gene. 2016; 595:83-91.
- Wang Z, Gerstein M, Snyder M. RNA-Seq: a revolutionary tool for Transcriptomics. Nature Reviews Genetics. 2009; 10(1):57-63.
- Mortazavi A, Williams BA, McCue K, Schaeffer L. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nature Methods. 2008;5(7): 621-628.
- Eric PN, Sarah WB, Alex B, Jennifer D, Ruth YE, Sean RE, et al. Rfam 12.0: updates to the RNA families database. Nucleic Acids Research. 2014;43: D130-137.
- Kozomara A, Griffiths-Jones S. miRBase: annotating high confidence microRNAs using deep sequencing data. NAR. 2014;42: D68-D73.
- Friedlander MR, Chen W, Adamidic C, Maaskoia J, Einspanier R, Knespel S, et al. Discovering microRNAs from deep sequencing data using miRDeep. Nature Biotechnology 2008;26:407-415.
- Peter AC, Yavuz A, Heleneb HT, Erno V, Rolf H, Renée X, Judith M, et al. Deep sequencing-based expression analysis shows major advances in robustness, resolution and inter-lab portability over five microarray platforms. Nucleic Acids Res. 2008;36(21): e141.
- Jiang H, Wong WH. Statistical inferences for isoform expression in RNASeq. Bioinformatics 2009;25:1026-1032.
- Yang YH, Dudoit S, Luu P, Lin DM, Peng V, Ngai J. Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. Nucleic Acids Res. 2002;30: e15.
- Wang L, Feng Z, Wang X, Wang X, Zhang X. DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics. 2010;26(1):136-138.
- John B, Enright AJ, Aravin AA, Tuschl T, Sander C, Marks DS. Human MicroRNA targets. Plos. Biology. 2005;3(7), e264.
- Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008;36: D480-484.
- Liu GY, Wu ZY, Zhu YL, Liu L, Li FC. Effects of dietary vitamin B 6 on the skeletal muscle protein metabolism of growing rabbits, Animal Production Science, 2017, 57, 2007–2015.
- LiLivak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2 -ΔΔCT method. Methods. 2001; 25: 402–408.
Table 1. Effects of ocu-miR-205 on proliferation, cell cycle and apoptosis of dermal papilla cells (DPCs; %)
|
Items |
Group |
R-MSE |
P-value |
|||
|
Control |
HBAD-GFP |
HBAD-ocu-miR-205-GFP |
HBAD-ocu-miR-205-5p-sponge-GFP |
|||
|
Proliferation of dermal papilla cells |
||||||
|
Optical density (OD) value |
0.56±0.01b |
0.55±0.01b |
0.61±0.01a |
0.34±0.01c |
0.0356 |
<0.0001 |
|
Cell cycle of dermal papilla cells |
||||||
|
Resting state/first gap (G0/G1) |
87.79±1.00a |
84.68±1.01b |
77.71±0.88c |
85.57±1.03ab |
2.7795 |
<0.0001 |
|
Synthesis (S) |
7.15±0.85b |
9.48±0.95b |
14.73±0.62a |
8.97±0.68b |
2.2232 |
<0.0001 |
|
Second gap/mitosis (G2/M) |
5.07±0.24b |
5.37±0.49b |
7.56±0.37a |
5.46±0.35b |
1.0523 |
0.0002 |
|
Apoptosis of dermal papilla cells |
||||||
|
Early apoptotic ratio (Q4) |
33.95±0.40b |
32.86±1.04b |
36.71±0.67a |
29.16±0.47c |
1.9470 |
<0.0001 |
|
Later apoptotic ratio (Q2) |
33.46± 0.57b |
32.15±1.33b |
36.11±0.73a |
31.40±0.57b |
2.4323 |
0.0032 |
|
Total apoptosis ratio (Q4+Q2) |
67.41±0.51b |
65.01±1.50b |
72.83±1.14a |
60.56±0.47 b |
2.8310 |
<0.0001 |
Note: Data shown are mean values ± s. d., and n = 8 per group. In the same row, values with same letter superscripts mean no significant difference (P>0.05), with different letter superscripts mean significant difference (P<0.05).
Table 2. Effects of ocu-miR-205 on gene expression of signal pathway of dermal papilla cells (DPCs)
|
Gene |
Group |
R-MSE |
P-value |
|||||||||||
|
Control |
HBAD-GFP |
HBAD-ocu-miR-205-GFP |
HBAD-ocu-miR-205-5p-sponge-GFP |
|||||||||||
|
PI3K/Akt signal pathway |
||||||||||||||
|
Inppl1 |
1.00±0.09b |
1.02±0.04b |
0.80±0.07b |
1.68±0.12a |
0.2353 |
<0.0001 |
||||||||
|
Inpp4b |
1.00± 0.09 |
0.97±0.16 |
0.92±0.08 |
1.21±0.16 |
0.3598 |
0.3966 |
||||||||
|
Frk |
1.00±0.34b |
0.90±0.07b |
0.50±0.13c |
2.63±0.47a |
0.8442 |
0.0001 |
||||||||
|
Phlda3 |
1.00± 0.04b |
0.89±0.06b |
0.69±0.05c |
1.18±0.05a |
0.1491 |
<0.0001 |
||||||||
|
Wnt signal pathway |
||||||||||||||
|
Wnt10b |
1.00±0.11b |
0.95±0.09b |
0.38±0.11c |
2.24±0.17a |
0.3519 |
<0.0001 |
|
|||||||
|
CTNNB1 |
1.00±0.09b |
1.50±0.13b |
1.48±0.19b |
3.16±0.27a |
0.5186 |
<0.0001 |
||||||||
|
GSK-3β |
1.00± 0.14b |
0.97±0.05b |
2.66±0.68a |
3.22±0.71a |
1.4118 |
0.0050 |
|
|||||||
|
DKK1 |
1.00±0.09 c |
2.90±0.13b |
6.09±0.84 a |
0.84±0.07 c |
1.2141 |
<0.0001 |
||||||||
|
Notch signal pathway |
||||||||||||||
|
Notch1 |
1.00±0.08b |
0.96±0.07b |
0.39±0.02c |
1.19±0.03a |
0.1604 |
<0.0001 |
||||||||
|
Jagged1 |
1.00±0.07 ab |
0.81±0.13 bc |
0.60±0.10c |
1.12±0.07 a |
0.2679 |
0.0032 |
||||||||
|
Hes1 |
1.00±0.12b |
1.19±0.16b |
0.90±0.05b |
2.40±0.33a |
0.5475 |
<0.0001 |
||||||||
|
Hes5 |
1.00± 0.16 b |
1.04±0.06 b |
0.18±0.04c |
1.38±0.04a |
0.2477 |
<0.0001 |
||||||||
|
BMP signal pathway |
||||||||||||||
|
BMP2 |
1.00±0.15b |
0.51±0.05bc |
1.69±0.34a |
0.41±0.05c |
0.5408 |
0.0002 |
||||||||
|
BMP4 |
1.00± 0.20 ab |
1.11±0.18 ab |
1.34±0.08a |
0.63±0.13b |
0.4405 |
0.0254 |
||||||||
|
TGF-β1 |
1.00±0.05 ab |
0.89±0.18 b |
1.33±0.13a |
0.71±0.07 b |
0.3304 |
0.0071 |
||||||||
Note: Data shown are mean values ± s. d., and n = 8 per group. In the same row, values with no letter superscripts or same letter superscripts mean no significant difference (P>0.05), with different letter superscripts mean significant difference(P<0.05).
Table 3. Effects of ocu-miR-205 on gene expression of signal pathway of Rex rabbits skin
|
Gene |
Group |
R-MSE |
P-value |
||||||
|
Control |
HBAD-GFP |
HBAD-ocu-miR-205-GFP |
HBAD-ocu-miR-205-5p-sponge-GFP |
||||||
|
PI3K/Akt signal pathway |
|||||||||
|
Inppl1 |
1.00±0.09b |
0.99±0.07b |
0.54±0.03c |
2.30±0.15a |
0.2731 |
<0.0001 |
|||
|
Inpp4b |
1.00± 0.15ab |
0.89± 0.07a |
0.09±0.01b |
1.05±0.02a |
0.2348 |
<0.0001 |
|||
|
Frk |
1.00±0.26a |
1.16±0.05b |
0.31±0.02c |
1.94±0.20a |
0.4666 |
<0.0001 |
|||
|
Phlda3 |
1.00± 0.11b |
1.18± 0.11a |
0.27±0.03b |
1.18±0.11a |
0.2500 |
<0.0001 |
|||
|
Wnt signal pathway |
|||||||||
|
Wnt10b |
1.00±0.14b |
0.86±0.04bc |
0.74±0.04c |
1.26±0.09a |
0.2406 |
0.0011 |
|||
|
CTNNB1 |
1.00±0.07b |
0.98±0.05b |
0.75±0.04c |
1.30±0.08a |
0.1789 |
<0.0001 |
|||
|
GSK-3β |
1.00± 0.06ba |
1.04±0.07b |
0.85±0.04b |
1.10±0.05a |
0.1496 |
0.0198 |
|||
|
DKK1 |
1.00±0.36a |
0.91±0.10a |
1.28±0.10a |
0.18±0.07b |
0.5625 |
0.0041 |
|||
|
Notch signal pathway |
|||||||||
|
Notch1 |
1.00±0.18b |
0.83±0.05bc |
0.54±0.05c |
1.95±0.22a |
0.4131 |
<0.0001 |
|||
|
Jagged1 |
1.00±0.13a |
1.02±0.05 a |
0.36±0.03c |
1.34±0.22 a |
0.3670 |
0.0001 |
|||
|
Hes1 |
1.00±0.09b |
1.03±0.06b |
0.80±0.03b |
1.49±0.18a |
0.3042 |
0.0009 |
|||
|
Hes5 |
1.00± 0.13 b |
1.00±0.05b |
0.36±0.03c |
1.43±0.20a |
0.3435 |
<0.0001 |
|||
|
BMP signal pathway |
|||||||||
|
BMP2 |
1.00±0.15b |
1.02±0.07b |
1.39±0.08a |
0.64±0.10c |
0.3051 |
0.0005 |
|||
|
BMP4 |
1.00± 0.15b |
0.87±0.04b |
1.58±0.14a |
0.73±0.05b |
0.3081 |
<0.0001 |
|||
|
TGF-β1 |
1.00±0.12 ab |
0.88±0.05 bc |
1.15±0.06a |
0.67±0.05b |
0.2222 |
0.0014 |
|||
Note: Data shown are mean values ± s. d., and n = 8 per group. In the same row, values with same letter superscripts mean no significant difference (P>0.05), with different letter superscripts mean significant difference (P<0.05).
Table 4. Effects of ocu-miR-205 on hair follicle density of Rex rabbits (Count/mm2)
|
Items |
Group |
R-MSE |
P-value |
|
|||||||||
|
Control |
HBAD-GFP |
HBAD-ocu-miR-205-GFP |
HBAD-ocu-miR-205-5p-sponge-GFP |
|
|||||||||
|
Transfection for 7 days |
|
|
|
|
|
|
|
||||||
|
Hair follicle density |
125.57±9.63bc |
144.25±11.47ba |
105.02±5.89c |
163.94±16.68a |
32.7817 |
0.0085 |
|
||||||
|
Primary hair follicle density |
10.08±0.85 |
9.31±0.85 |
7.80±0.63 |
9.79±0.77 |
2.2046 |
0.1902 |
|
||||||
|
Secondary hair follicle density |
115.49±9.48bc |
134.93±12.01ba |
97.22±5.68c |
154.15±16.65a |
32.9728 |
0.0113 |
|
||||||
|
Secondary/Primary (S/P) ratio |
11.92±1.17 |
16.02±2.67 |
12.92±1.06 |
16.38±2.04 |
5.2469 |
0.2536 |
|
||||||
|
Transfection for 14 days |
|
|
|
|
|
|
|
||||||
|
Hair follicle density |
130.54±12.22b |
136.62±8.51b |
119.30±7.24b |
203.34±14.79a |
31.4005 |
<0.0001 |
|
||||||
|
Primary hair follicle density |
10.86±0.54 |
9.96±0.50 |
11.92±0.74 |
10.11±0.70 |
1.7739 |
0.1307 |
|
||||||
|
Secondary hair follicle density |
119.68±12.03b |
126.66±8.28b |
107.37±7.10b |
213.32±14.62a |
30.9035 |
0.0001 |
|
||||||
|
Secondary/Primary (S/P) ratio |
11.08±0.96 b |
12.82±0.82 b |
9.20±0.78 b |
21.62±1.94 a |
3.4505 |
<0.0001 |
|
||||||
|
Transfection for 21 days |
|
|
|
|
|
|
|
||||||
|
Hair follicle density |
145.43±13.71b |
139.58±16.53b |
126.38±9.30b |
208.69±21.42a |
44.8651 |
0.0049 |
|||||||
|
Primary hair follicle density |
14.30±1.19 |
12.65±0.80 |
12.92±0.51 |
15.29±1.28 |
2.8084 |
0.2244 |
|||||||
|
Secondary hair follicle density |
131.13±13.21b |
126.94±16.57b |
113.46±9.27b |
193.40±21.31a |
44.4697 |
0.0060 |
|||||||
|
Secondary/Primary (S/P) ratio |
9.45±0.91 |
10.28±1.40 |
8.86±0.75 |
13.28±1.85 |
3.6808 |
0.1018 |
|||||||
Note: Data shown are mean values ± s. d., and n = 8 per group. In the same row, values with no letter superscripts or same letter superscripts mean no significant difference (P>0.05), with different letter superscripts mean significant difference (P<0.05).
Additional file 1:
Supplementary Figure 1. Skin crosscutting HE staining of Rex rabbits with different hair density（100 magnification）(a) Low hair density (b) High hair density
Supplementary Table 1. Statistics of hair follicle numbers of Rex rabbits with different hair density
Note: Data shown are mean values ± s. d., and n =3 per group. In the same row, values with with different letter superscripts mean significant difference(P<0.05).
Additional file 2:
Supplementary Table 2. RNA quality of six DPC samples
Supplementary Figure 2. RNA quality (a) Agarose electrophoresis of RNA in DPCs; (b) Quality test results of total RNA by Agilent 2100
Additional file 3:
Supplementary Table 3. Statistics of sequence data of each samples
Additional file 4:
Supplementary Figure 3. Length distribution map of Small RNA
The X-axis is the length of small RNA, and the Y-axis is the corresponding number of small RNA
Additional file 5:
Supplementary Table 4. Alignment clean tag in genome
Additional file 6:
Supplementary Figure 4. Statistical distribution map of small RNA types
In order to make each unique small RNA have a unique annotation, the annotation statistics of small RNA traverse the annotation according to the priority order of miRNA > piRNA > snoRNA > Rfam > other sRNA.
Additional file 7:
Supplementary Figure 5. The quantity distribution maps of clean tag base in each sample
The X-axis is the position of the base in read, and the Y-axis represents the proportion of the base.
Additional file 8:
Supplementary Figure 6. The quality distribution maps of the clean tag base in each sample
The X-axis is the position of base in read, and the Y-axis represents the base mass value.
Additional file 9:
Supplementary Table 5. LD-vs-HD_DEGseq. diffexpfilter
Additional file 10:
Supplementary Figure 7. Directed acyclic graph of GO enrichment
(a) LD-vs-HD_DEGseq. Biological_Process. Top GO；LD-vs-HD_DEGseq. Cellular_Component. Top GO；LD-vs-HD_DEGseq. Molecular_Function. Top GO.
The branches in the figure represent the inclusion relationship, and the function range defined from top to bottom is getting smaller and smaller. The box represents the top 5 GO terms of each classification enrichment degree, and through the inclusion relationship, the associated GO terms are displayed together. The name of the term and the P-value corrected by enrichment analysis are displayed on each node. The darker the color (red) indicates the smaller the P-value and the higher the enrichment degree.
Additional file 11:
Supplementary Table 6. Primer sequence information in experiments
This is a list of supplementary files associated with this preprint. Click to download.
- SupplementaryFigure7b.jpg
- SupplementaryFigure7a.jpg
- SupplementaryFigure7c.jpg
- ARRIVEchecklist.docx
- SupplementaryFigure5LD3.png
- SupplementaryFigure6HD3.png
- SupplementaryFigure6LD1.png
- SupplementaryFigure6HD2.png
- SupplementaryFigure6LD2.png
- SupplementaryFigure6LD3.png
- SupplementaryFigure4LD3.png
- SupplementaryFigure5HD1.png
- SupplementaryFigure5HD3.png
- SupplementaryFigure5HD2.png
- SupplementaryFigure5LD1.png
- SupplementaryFigure5LD2.png
- SupplementaryFigure6HD1.png
- SupplementaryFigure1b.jpg
- SupplementaryFigure4HD1.png
- SupplementaryFigure4LD1.png
- SupplementaryFigure4LD2.png
- SupplementaryFigure4HD2.png
- SupplementaryFigure4HD3.png
- SupplementaryFigure3HD3.png
- SupplementaryFigure2bLD1.png
- SupplementaryFigure1a.jpg
- SupplementaryFigure3LD1.png
- SupplementaryFigure3LD3.png
- SupplementaryFigure3LD2.png
- SupplementaryFigure2bLD2.png
- SupplementaryTable6.docx
- SupplementaryFigure3HD1.png
- SupplementaryFigure2bHD1.png
- SupplementaryFigure2bLD3.png
- SupplementaryFigure3HD2.png
- SupplementaryTable14.docx
- SupplementaryFigure2bHD3.png
- SupplementaryFigure2bHD2.png
- SupplementaryFigure2a.png
- SupplementaryTable5.LDvsHDDEGseq.diffexpfilter.xls
Badges
Prescreen
This manuscript passed our Prescreen assessment
Complete author information
Appropriate declaration statements
Potential risks to human health
Prescreen
Peer Review Timeline
CURRENT STATUS: Under Review
Version 6
Posted 13 Apr, 2020
No community comments so far
No comments provided
Editorial decision: Minor revision
On 06 Apr, 2020
Editor assigned
On 01 Apr, 2020
Submission checks complete
On 01 Apr, 2020
Editor invited
On 01 Apr, 2020
No comments provided
Submission checks complete
On 31 Mar, 2020
Editorial decision: Minor revision
On 30 Mar, 2020
No comments provided
Submission checks complete
On 25 Mar, 2020
Editorial decision: Minor revision
On 24 Mar, 2020
No comments provided
Editorial decision: Minor revision
On 18 Mar, 2020
Review #2 received
Received 12 Mar, 2020
Reviewer #2 agreed
On 09 Mar, 2020
Reviewers invited
Invitations sent on 09 Mar, 2020
Reviewer #1 agreed
On 09 Mar, 2020
Review #1 received
Received 09 Mar, 2020
Editor assigned
On 08 Mar, 2020
Submission checks complete
On 08 Mar, 2020
Editor invited
On 08 Mar, 2020
No comments provided
Editorial decision: Major revision
On 25 Feb, 2020
Review #2 received
Received 21 Feb, 2020
Reviewer #2 agreed
On 07 Feb, 2020
Reviewer #1 agreed
On 02 Feb, 2020
Review #1 received
Received 02 Feb, 2020
Reviewers invited
Invitations sent on 16 Oct, 2019
Editor invited
On 27 Aug, 2019
Editor assigned
On 21 Aug, 2019
Submission checks complete
On 21 Aug, 2019
Metrics
Comments: 0
PDF Downloads: ...
HTML Views: ...
Learn more about our company.
This preprint is under consideration at BMC Developmental Biology. Preprints are preliminary reports that have not undergone peer review. They should not be considered conclusive, used to inform clinical practice, or referenced by the media as validated information.
Learn more about In Review
Research article
The functions of ocu-miR-205 in regulating hair follicle development in Rex rabbits
Fuchang Li
Shandong Agricultural University
Corresponding Author
ORCiD: https://orcid.org/0000-0002-5670-8038
Badges
Prescreen
This manuscript passed our Prescreen assessment
Complete author information
Appropriate declaration statements
Potential risks to human health
Prescreen
Peer Review Timeline
CURRENT STATUS: Under Review
Version 6
Posted 13 Apr, 2020
No community comments so far
No comments provided
Editorial decision: Minor revision
On 06 Apr, 2020
Editor assigned
On 01 Apr, 2020
Submission checks complete
On 01 Apr, 2020
Editor invited
On 01 Apr, 2020
No comments provided
Submission checks complete
On 31 Mar, 2020
Editorial decision: Minor revision
On 30 Mar, 2020
No comments provided
Submission checks complete
On 25 Mar, 2020
Editorial decision: Minor revision
On 24 Mar, 2020
No comments provided
Editorial decision: Minor revision
On 18 Mar, 2020
Review #2 received
Received 12 Mar, 2020
Reviewer #2 agreed
On 09 Mar, 2020
Reviewers invited
Invitations sent on 09 Mar, 2020
Reviewer #1 agreed
On 09 Mar, 2020
Review #1 received
Received 09 Mar, 2020
Editor assigned
On 08 Mar, 2020
Submission checks complete
On 08 Mar, 2020
Editor invited
On 08 Mar, 2020
No comments provided
Editorial decision: Major revision
On 25 Feb, 2020
Review #2 received
Received 21 Feb, 2020
Reviewer #2 agreed
On 07 Feb, 2020
Reviewer #1 agreed
On 02 Feb, 2020
Review #1 received
Received 02 Feb, 2020
Reviewers invited
Invitations sent on 16 Oct, 2019
Editor invited
On 27 Aug, 2019
Editor assigned
On 21 Aug, 2019
Submission checks complete
On 21 Aug, 2019
Metrics
Comments: 0
PDF Downloads: ...
HTML Views: ...
License:
This work is licensed under a Creative Commons Attribution 4.0 International License. Read Full License
Background: Hair follicles is an appendage from the vertebrate skin epithelium, and arise from the embryonic ectoderm andregenerate cyclically during adult life. Dermal papilla cells (DPCs) is the key dermal component of the hair follicle that directly regulates hair follicle development, growth and regeneration. Recent studies have reported that miRNA plays an important role in regulating hair follicle morphogenesis, proliferation, differentiation and apoptosis of hair follicle stem cells.
Results: The miRNAs expression profile of the DPCs from different hair density Rex rabbits shown that 240 differentially expressed of miRNAs were screened (log 2 (HD/LD)｜>1.00 and Q-value≤0.001). Among them, the expression of ocu-miR-205-5p in low hair densities DPCs was higher than that in high hair densities, and it is highly expressed in the skin tissue of Rex rabbits ( P <0.05). ocu-miR-205 could increase cell proliferation and cell apoptosis ratio, change cell cycle process ( P <0.05), affect the genes expression of PI3K/Akt, Wnt, Notch and BMP signaling pathways in DPCs and skin tissue of Rex rabbits, inhibit the protein phosphorylation level of CTNNB1, GSK-3β and the protein expression level of noggin (NOG), promote Akt phosphorylation level ( P <0.05). There was no significant change in primary follicle density ( P >0.05), but the secondary follicle density and total follicle density ( P <0.05) were changed after ocu-miR-205-5p interfered expression, and secondary/primary ratio (S/P) in ocu-miR-205-5p interfered expression group increased at 14 days after injection ( P <0.05).
Conclusion: ocu-miR-205 could promote the apoptosis of DPCs, affect PI3K/Akt, Wnt, Notch and BMP signaling pathways genes and proteins expression in DPCs and skin of Rex rabbits, promote the transformation of hair follicles from growth phase to regression and resting phase, and affect hair density of Rex rabbits.
Figure 1
Figure 2
Figure 3
Figure 4
Figure 5
Figure 6
Rex rabbits are typically used for fur production because their fur is short, fine, dense, and smooth and has important economic value [1]. The most important indicator for evaluating the fur quality of Rex rabbit is the density of hair follicles [2]. In recent years, methods to improve the hair follicle density of Rex rabbits have become the most important concern in rabbit production. Hair follicles are an appendage from the vertebrate epithelium in the skin, arise from the embryonic ectoderm, and regenerate cyclically during adulthood [3]. The process of hair follicle formation and differentiation involves at least 20 different cells and tissues [4], including the dermal papilla, hair matrix, inner root sheath, and outer root sheath, as well as different signalling pathways, such as the Wnt, Notch, bone morphogenetic protein (BMP), and fibroblast growth factor (FGF) pathways, among others [5-9]. Phosphatidylinositol 3'-kinase (PI3K) preferentially phosphorylates PIP2 to produce PIP3, and PIP3 is an important second messenger in cells that subsequently activates Akt, and thus it plays an important role in regulating the proliferation and apoptosis of hair follicle cells. The Wnt signalling pathway regulates epithelial morphogenesis, hair follicle development and cell differentiation. The PI3K/Akt signalling pathway inhibits the phosphorylation of β-catenin by phosphorylating glycogen synthase kinase 3β (GSK-3β) and activates the Wnt signalling pathway[5]. Dickkopf-related protein 1 (DKK1) inhibits the Wnt signalling pathway by inhibiting the phosphorylation of β-catenin and induces hair follicle regression [10]. When Notch receptor binds to a ligand, it activates hair follicle stem cells and then promotes the transition of hair follicles from the resting stage to growing stage [6]. The BMP signalling pathway is involved in embryonic skin appendage organ morphogenesis and postnatal hair follicle growth [7]. The BMP2 and BMP4 genes inhibit hair follicle development and are associated with maintaining hair follicles in the resting stage[8]. Noggin (NOG) is an inhibitor of the BMP signalling pathway, and its abnormal expression leads to follicular enlargement [7]. The Notch signalling pathway interacts with the BMP signalling pathway, and the BMP signalling pathway inhibits the Wnt signalling pathway by regulating β-catenin expression [6]. Dermal papilla cells (DPCs) is the key dermal component of the hair follicle, and it can directly regulate hair follicle development, growth and regeneration [11]. Besides, the characteristic of DPCs also affects the size, shape and cycling of hair follicle [12]. Signal exchange between DPCs and hair follicle stem cells at telophase is the key to initiating the next hair follicle cycle [13].
MicroRNAs (miRNAs) are small non-coding RNAs that have been plays an important role in embryogenesis, organ development, cell proliferation, metabolism, embryogenesis, behaviour and other biological processes [14-19]. Recently, miRNAs were reported to play important roles in regulating hair follicle morphogenesis and the proliferation, differentiation and apoptosis of hair follicle stem cells in mice, rats, goats and sheep [20-22]. Notably, miR-214 inhibits hair follicle growth and development by modulating the expression of regulatory factors in Wnt signalling pathway, such as β-catenin and lymphoid enhancer-binding factor 1 (Lef-1) [23]. Upon the overexpression of DKK1 in transgenic mice, the expression of miR-200b and miR-196a in epidermis decreases significantly, which is possibly mediated by potential target genes acting on the Wnt signalling pathway [24]. As shown in a previous study, miR-let-7b promotes alpaca hair growth by inhibiting the transcription of transforming growth factor β receptor 1 (TGFβR1) [25]. BMP4 negatively regulates the expression of miR-21, and miR-21 negatively regulates the expression of the BMP-dependent tumour suppressor genes Pten, Pdcd4, Timp3 and Tpm1 [26]. Hair follicle development and hair growth in mice are regulated by miR-31 through effects on the BMP and Wnt signalling pathways [20]. Additionally, miR-205 is a highly conserved miRNA that shares a similar expression pattern with miR-200 family [27]. It is one of the miRNAs expressed at the highest levels in the epidermis [28, 29], and it plays an essential role in promoting the neonatal expansion of skin stem cells during early development by modulating the PI3K pathway [30].
We isolated the DPCs from the skin of Rex rabbits and analysed the miRNA expression profiles of the DPCs from Rex rabbits with different hair densities to improve the fur quality of Rex rabbits. Among the miRNAs analysed, ocu-miR-205 was one of the miRNAs expressed at the highest levels in DPCs. We examined the function of ocu-miR-205 in hair follicle development.
DPCs show a complex miRNA expression pattern
The varieties of miRNAs in DPCs from 30-day-old Rex rabbits with lower and higher hair densities (Additional file 1) were studied by subjecting RNA samples with a high integrity and qualified quality (Additional file 2) to high-throughput small RNA sequencing using the BGISEQ-500 platform. 37930744, 39442011, 40965907, 38502653, 40622117 and 41149163 clean reads were obtained from the six samples (Additional file 3), and the majority of clean reads had a length of 23 nucleotides (Additional file 4). Comparison with known small RNA databases, the percentage of matching reads for each library was 91.91 %, 91.05 %, 92.46 %, 92.77 %, 90.13 %, and 91.61 %, respectively (Additional file 5). The results of the small RNA classification showed that miRNAs accounted for 80.80 %, 82.50 %, 81.60 %, 85.80 %, 76.50 %, and 81.80 %, respectively (Additional file 6). The base Q20 value of filtered data was greater than 90 % and the Q30 value was greater than 80 %. Additional files 7 and 8 show the quantity and quality distribution maps of the bases in each sample, respectively.
By screening differentially expressed genes (DEGs), 240 differentially expressed miRNAs were identified (｜log2(HD/LD)｜>1.00 and Q-value≤0.001; Fig. 1a and Additional file 9), including 122 miRNAs that were upregulated and 118 miRNAs that were downregulated (Fig. 1b). The annotation of the target genes of differentially expressed miRNAs was performed using Gene Ontology (GO) enrichment, and a total of 205,661 target genes were enriched in GO terms (Fig. 1c). Specific GO terms of the target genes mainly involved in the biological process (BP), cellular component (CC) and molecular function (MF) categories. A directed acyclic graph (DAG) of the enriched GO terms in the BP, CC and MF categories is provided in Additional file 10. Following, the target genes were also uploaded into the Kyoto Encyclopaedia of Genes and Genomes (KEGG) database to identify the pathways that were actively regulated by miRNAs in DPCs. Three hundred twenty-five pathways were predicted, and the top enriched pathways were the Wnt signalling pathway and Notch signalling pathway, among others (Fig. 1d).
The expression of ocu-miR-205
The structure of one of the differentially expressed miRNAs, ocu-miR-205, is shown in Fig. 2a. The expression of ocu-miR-205 differed in DPCs from Rex rabbits with different hair densities (｜log2(HD/LD)｜>1.00 and Q-value≤0.001; Fig. 2b) and was expressed at low levels in DPCs from Rex rabbits with a high hair density (HD) and at high levels in DPCs from Rex rabbits with a low hair density (LD). Quantitative fluorescence PCR results also confirmed the accuracy of the ocu-miR-205-5p sequencing results (Fig. 2c). Furthermore, the expression of ocu-miR-205-5p in different tissues from Rex rabbits differed (P<0.05). The expression in the skin tissue was higher than in the stomach, intestine, spleen, liver, heart, lung, kidney, muscle and fat (P<0.05; Fig. 2d), suggesting a tissue-specific expression pattern.
Effects of ocu-miR-205 on DPCs from Rex rabbits
HBAD-GFP, HBAD-ocu-miR-205-GFP, and HBAD-ocu-miR-205-5p-sponge-GFP were constructed and used to transfect DPCs at a multiplicity of infection (MOI) of 200. Cell proliferation, the cell cycle, cell apoptosis, and changes in the expression of ocu-miR-205-5p and hair follicle development-related genes and proteins were detected 48 h after transfection. The constructed adenoviruses successfully transfected DPCs (Fig. 3a), and ocu-miR-205-5p expression increased significantly after overexpression, but decreased significantly after silencing (P<0.05; Fig. 3b). Importantly, ocu-miR-205 increased cell proliferation and the cell apoptosis rate, in addition to altering the progression of the cell cycle (P<0.05; Table 1). Moreover, ocu-miR-205 inhibited the expression of the Inppl1, Frk and Phlda3 mRNAs, which are involved in the PI3K/Akt signalling pathway (P<0.05). The expression of ocu-miR-205 induced the expression of the DKK1 mRNA and inhibited the expression of the Wnt10b, CTNNB1 and GSK-3β genes in the Wnt signalling pathway (P<0.05). The expression of the Notch1, Jagged1, Hes1 and Hes5 genes in the Notch signalling pathway was suppressed by ocu-miR-205 (P<0.05), whereas the expression of the BMP2, BMP4 and TGF-β1 genes in the BMP signalling pathway was induced (P<0.05; Table 2). The expression of ocu-miR-205 inhibited the phosphorylation of the CTNNB1 and GSK-3β proteins, decreased the level of the noggin (NOG) protein, and increased the level of Akt phosphorylation (P<0.05; Fig. 4).
Effects of ocu-miR-205 on the skin tissue of Rex rabbits
One hundred 3-month-old Rex rabbits with similar body weights were randomly divided into 4 groups. After local shaving of the back skin, 50 μL of HBAD-GFP, HBAD-ocu-miR-205-GFP and ocu-miR-205-5p-sponge-GFP were injected intradermally into the rabbits in each group. 24 h after the injected, one randomly selected Rex rabbit from each group was euthanized, and frozen sections were prepared from the locally injected skin to confirm the expression of the miRNAs from the injected adenoviruses. Eight Rex rabbits were randomly selected from each group 7 days, 14 days and 21 days after transfection. After euthanasia, the injected skin was collected and divided into two parts. A part placed in freezing tube and the other part was fixed at a 4 % paraformaldehyde fixative solution. The frozen samples were used to detect the changes in the expression of ocu-miR-205-5p and hair follicle development-related genes and proteins. The fixed samples were used to prepare paraffin sections, and the hair density was counted after HE staining. The constructed adenoviruses successfully transfected the skin and hair follicles of Rex rabbits (Fig. 5a), and the expression of ocu-miR-205-5p increased significantly after overexpression, but decreased significantly after silencing (P<0.05; Fig. 5b). Fourteen days after the injection, ocu-miR-205 significantly altered the expression of genes involved in the PI3K/Akt, Wnt, Notch and BMP signalling pathways (P<0.05; Table 3). Moreover, ocu-miR-205 significantly altered the levels of the phosphorylated CTNNB1, GSK-3β, and Akt proteins and the level of the NOG protein (P<0.05; Fig. 6). A significant change in the primary follicle density was not observed (P>0.05), but the secondary follicle density and total follicle density (P<0.05) were altered after ocu-miR-205-5p interfered expression, and the secondary/primary ratio (S/P) in the ocu-miR-205-5p interfered expression group increased at 14 days after the injection (P<0.05; Table 4).
Rex rabbits mainly used for fur production, and the fur is the most important livestock products. The hair follicle cycle can typically be divided into anagen, catagen, and telogen phases[31]. The hair follicle of Rex rabbits contains primary and secondary hair follicles, and the different sizes and types of primary and secondary hair follicles is easy distinguished[32]. Besides, the hair follicle growth cycle is easy distinguish the differentiation of hair cycle phases [33]. Many factors affect the growth and development of hair follicles in Rex rabbits, such as genetic, nutritional, age, among others [34-37]. The most hair follicles in anagen phases when the Rex rabbits at 4-5 weeks[37]. Notably, miRNAs plays an important functions in regulating gene expressions[38-40]. Among them, miR-205 is abundantly expressed in the epidermis[41,42], and it plays an essential role by modulating the PI3K/Akt signalling pathway in promoting the neonatal expansion of skin stem cells during early development[43]. In vitro models, miR-205 promotes keratinocyte migration by targeting the lipid phosphatase SHIP2 and KIR4.1[44, 45]. There is a potential therapeutic application of DPCs in the treatment of alopecia, because DPCs of various origins induce the de novo formation of the hair follicle structure in both follicular and afollicular epidermis [46-49]. In the present study, ocu-miR-205-5p was expressed at significantly higher levels in DPCs from rabbits with low hair densities than in DPCs from rabbits with high hair densities. Moreover, ocu-miR-205 induced G0/G1 arrest in DPCs, which was further confirmed by the reduction in the population of DPCs in the G0/G1 phase and increase in the apoptotic rate after the transfection of the ocu-miR-205 inhibitor. This finding is consistent with previous studies [50]. Furthermore, overexpression of miR-205 significantly inhibited cell proliferation. Additionally, ocu-miR-205 inhibited the expression of related gene and proteins in the PI3K/Akt, Wnt, and Notch signalling pathways, and activated the BMP signalling pathway. Therefore, ocu-miR-205 plays an important role in regulating hair follicle development.
In conclusion, ocu-miR-205 promoted the apoptosis of DPCs, inhibited cell proliferation, modulated the expression of genes and proteins involved in the PI3K/Akt, Wnt, Notch and BMP signalling pathways in DPCs and the skin of Rex rabbits, promoted the transition of hair follicles from the growth phase to the regression and resting phase, and altered the hair density of Rex rabbits.
Animal and sample collection
The Rex rabbits used in this study were purchased from Taishan Rabbit Farm (Shandong, China). Ten adult female Rex rabbits with a high (>14000/cm2, HD) or low hair density (<10000/cm2, LD) were chosen and divided into 2 groups. An adult male rabbit was selected to mate with the female rabbits within each group. The first generation of offspring (F1 generation) was obtained and then inbred with F1 in the same nest (HD and HD, LD and LD inbreeding). Three rabbits each with high and low hair densities were selected from the second generation of offspring (F2 generation) 30 days after birth. Selected animals were electrically stunned (120 V, pulsed direct current, 50 Hz for 5 s) and euthanized by exsanguination of the carotid artery and before skin was harvested from the experimental rabbits. The DPCs were separated and cultured, methods according to references[51]. After the cells had established a monolayer (approximately 12 days), the total RNA was extracted, tested for quality.
Construction of the small RNA libraries, sequencing analysis, miRNA identification and prediction of new miRNAs in DPCs
The methods of construction small RNA libraries according to references[52]. Solexa sequencing by synthesis method using the BGISEQ-500 platform at Shenzhen Huada Biotech Co., Ltd. (Shenzhen, China) [53,54]. The raw reads produced from sequencing were filtered to remove low-quality reads, and the clean reads were analysed using BLAST with Bowtie-1.0.0 software, Rfam [55]. The clean reads were used for miRNA identification and compared with the mature miRNAs and pre-miRNAs from Oryctolagus cuniculus listed in miRBase 21.0 [56]. Subsequently, the miRDeep 2 software was used to predict the novel miRNAs by exploring the secondary structure[57], minimum free energy and dicer cleavage sites of the unannotated clean reads that mapped to the Oryctolagus cuniculus genome. After the identification of conserved miRNAs, the clean reads were aligned to the Oryctolagus cuniculus genome for identify new Oryctolagus cuniculus miRNAs.
In order to eliminate the effect of different sequencing quantities on quantitative accuracy, transcripts per million (TPMs) were calculated to standardize the expression levels of small RNA [58]. Based on the assumption that RNA sequencing is a random process, the expression of each transcript is presumed to exhibit a binomial distribution[59]. DEGseq calculated differential expression based on MA-plot [60, 61], and the P-value of each miRNA was corrected by performing multiple hypothesis tests using Q-values. When the difference in coincidence was more than two-fold and the Q-value was less than or equal to 0.001, the miRNAs were considered significantly differentially expressed (｜log2(HD/LD)｜>1.00 and Q-value≤0.001).
The miRNAs were aligned to the EST unigenes of Oryctolagus cuniculus and the target genes were predicted using the miRanda algorithm to obtain a better understanding of the potential functions of the significantly differentially expressed miRNAs in Rex rabbits with different hair densities [62]. An enrichment analysis of the predicted target genes was conducted with GO terms and KEGG pathways [63].
Construction and identification of adenovirus vectors overexpressing and silencing ocu-miR-205
HBAD-GFP (HANBIO adenovirus-green fluorescent protein; empty vector), HBAD-ocu-miR-205-GFP (overexpression), HBAD-ocu-miR-205-5p-sponge-GFP (silencing) adenoviruses were synthesized and constructed by Hanheng Biotechnology Co., Ltd. (Shanghai, China). The infective titres of HBAD-GFP, HBAD-ocu-miR-205-GFP and HBAD-ocu-miR-205-5p-sponge-GFP were 1.26*1010 PFU/mL, 1.58*1010 PFU/mL and 1.26*1010 PFU/mL, respectively.
Third-generation DPCs displaying good growth conditions were inoculated into a disposable 6-well plate. The cell density was approximately 1.0*105 cells/mL. Prior to the infection, the virus was subjected to 10-fold gradient dilution. Generally, the MOI (multiplicity of infection) was controlled in the range of 10-1000. HBAD-GFP, HBAD-ocu-miR-205-GFP and HBAD-ocu-miR-205-GFP were individually transfected into Rex rabbit DPCs at an MOI 200, and a negative control was established using cells undergoing normal culture.
Fifty microliters of the purified adenovirus were injected into the skin of each Rex rabbit with microinjector at a concentration of 5.0*108 - 1.0 * 109 virus particles per Rex rabbit after shaving the middle part of the back of 100 3-month-old Rex rabbits with similar body weights and good health. 24 h after transfection, one Rex rabbit from each group was randomly selected, euthanized, and frozen sections were prepared from the locally injected skin. The adenovirus-transfected skin was observed under a positive fluorescence microscope (Nikon ECLIPSE 80i, Japan).
Assessment of the proliferation, cell cycle and apoptosis of DPCs
DPCs were plated in a 96-well plate at a density of 104 cells/well, cultured in basal medium for 24 h, and then transfected with the indicated adenoviruses. The methods assessment of the proliferation according to references[51]. Third-generation DPCs were plated in a disposable 6-well plate at a density of 104 cells/mL, with 2 mL of the cell suspension plated in each well. After a 24-h incubation to allow cells to adhere, the culture medium was removed. After treatment and culture for a certain time, the cells were digested with a trypsin digestion solution lacking EDTA (Solarbio, Beijing, China). The methods assessment of the cell cycle according to references[51]. For apoptosis assessment, cells were centrifuged and collected into a 1.5 mL centrifugal tube, and then washed with PBS. After centrifugation, 500 μL of 10X Annexin V Binding Buffer was added to re-suspended the cells, followed by the labelling of F-actin. Cells were incubated with the FITC Annexin V and Propidium Iodide Staining Solution for 15 min at 4 °C and then analysed using flow cytometry. The percentages of early apoptotic cells (Q4), late apoptotic cells (Q2) and total apoptotic cells (Q2 + Q4) in each sample were calculated.
Total RNA extraction and Real-time PCR analysis
The methods of total RNA extraction and real-time PCR of mRNA was according to references[64]. All quantitative PCR primers were designed using Primer Premier 5 software (Additional file 11). For the quantitative RT-PCR of miRNAs, 1 μg of total RNA was reverse transcribed with Bulge-Loop miRNA-specific reverse transcription primers (RiboBio, China), and quantitative PCR was performed using Fast Start Universal SYBR Green Master Mix (Roche Diagnostics GmbH Mannheim, Germany) and Bulge-Loop primers (RiboBio, Guangzhou, China) on the 7500 Fast System 1.4 system with small nuclear RNA U6 as the normalisation control. The volume of each reaction was 20 µL, including 2 µL of cDNAs, 10 µL of SYBR Green Master (2X), 0.8 μL of the Bulge-LoopTM miRNA Forward Primer (5 µM), 0.8 µL of the Bulge-LoopTM Reverse Primer (5 µM), 0.4 µL of ROX Reference Dye II (50×) and 6.0 µL of ddH2O. PCR was performed under the following conditions: 10 min of template denaturation at 95 °C, followed by 40 cycles of 95 °C for 2 s, 60 °C for 20 s, and 70 °C for 10 s. Melting curves (70 °C-95 °C) for each sample were analysed after each run to confirm the specificity of amplification reactions. The relative expression levels of mRNAs and miRNAs were calculated using the arithmetic formula 2−△△Ct [65].
Western immunoblotting.
The methods of total protein extraction and SDS-PAGE were according to references[64]. The membranes were blocked with 5 % skimmed milk in phosphate-buffered saline (PBS; Solarbio, China) at 4 °C overnight and incubated with primary antibodies (tubulin AT819, Beyotime, China; phospho-CTNNB1-S552 pAb, Abcam, US; phospho-GSK3B-S9 pAb, Abcam, US; phospho-AKT1-S473 pAb, Abcam, US; or NOG polyclonal antibody, Abcam, US). The membranes were then rinsed with Tris-buffered saline containing Tween (TBST; Solarbio, China), and subjected to detection with a 1:3000 dilution of a horseradish peroxidase (HRP)-conjugated goat anti-mouse IgG antibody (Beyotime, China) at 37 °C for 1 h. Proteins were visualized using BeyoECL reagents (Beyotime, China). The intensity of the bands was quantified with a Pro Plus 6.0 Biological Image Analysis System. The levels of phospho-CTNNB1, phospho-GSK3B, phospho-AKT1 and NOG were normalized to the internal control beta-tubulin, and the relative expression levels were calculated.
Statistical analysis
All data were analysed with SAS software (SAS version 8e; SAS Institute, Cary, NC, USA). A one-way ANOVA was used to evaluate the differences in mean values among various groups. The data are presented as the means and R-MSE. P < 0.05 was regarded as statistically significant.
BMP: bone morphogenetic protein; BP: biological process; CC: cellular component; DAG: Directed Acyclic Graph; DEGs: differentially expressed genes; DPCs: dermal papilla cells; DKK1: Dickkopf-related protein 1; FGF: fibroblast growth factor; GAPDH: glyceraldehyde-3-phosphate dehydrogenase; GFP: green fluorescent protein; GO: Gene Ontology; GSK-3β: glycogen synthase kinase 3β; HD: high hair density; HE: haematoxylin and eosin; Lef-1: lymphoid enhancer-binding factor 1; LD: low hair density; miRNAs: microRNAs; mRNA: messenger RNA; NOG: noggin; S/P: secondary/primary ratio; SDS: sodium dodecyl sulfate; TGFβ: transforming growth factor β; TPM: transcripts per million; ocu-miR-205: Oryctolagus cuniculus microRNA 205; PAGE: polyacrylamide gel electrophoresis; PBS: phosphate-buffered saline; PI3K: phosphatidylinositol 3'-kinase.
Ethics approval and consent to participate
The experimental procedures were approved by the Committee of Ethics in Research of Shandong Agricultural University (SDAUA-2017-029) and performed in accordance with the Guidelines for Experimental Animals of the Ministry of Science and Technology (Beijing, China). The experimental procedures in this test compliant with the ARRIVE guidelines, in accordance with the U.K. Animals (Scientific Procedures) Act, 1986 and associated guidelines, EU Directive 2010/63/EU for animal experiments, or with the National Institutes of Health guide for the care and use of Laboratory animals (NIH Publications No. 8023, revised 1978). This article does not any studies with human participants performed by any of the authors.
Consent for publication
Not applicable.
Availability of data and materials
Small RNA sequencing data has been deposited in Sequence Read Archive of the National Center for Biotechnology Information under temporary submission ID: SUB7193686, and BioSample accessions: SAMN14476334, SAMN14476335, SAMN14476336, SAMN14476337, SAMN14476338, SAMN14476339.
Competing interests
The authors declare that they have no competing interests.
Funding
This study was supported by the Earmarked Fund for Modern Agro-industry Technology Research System (CARS-43-B-1), National Natural Science Foundation of China (31972594), and Funds of Shandong “Double Tops” Programme (SYL2017YSTD11). The funding bodies had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Authors’ contributions
GL, SL and LB isolation, culture and identification of DPCs, and construction of the small RNA libraries, sequencing analysis, miRNAs identification and prediction of new miRNAs. GL and HS isolated and cultured hair follicles, and performed cellular experiments. GL, HL and SG raised the rabbits and conducted the animal experiments. GL, WJ and FL revised the manuscript. YZ, WJ and FL designed this work and wrote this manuscript. All authors read and approved the final manuscript.
Acknowledgements
Not applicable.
Author details
1 College of Animal Science and Technology, Shandong Agricultural University, Tai’an 271018, P. R. China; 2 Shandong Provincial Key Laboratory of Animal Biotechnology and Disease Control and Prevention, Tai’an 271018, P. R. China; 3 Animal Husbandry and Veterinary Institute, Shandong Academy of Agricultural Sciences, Jinan 251000, P. R. China；4 Shandong Key Laboratory of Animal Disease Control and Breeding, Jinan 251000, P. R. China.
- Ruben D. Choosing a Rex rabbit. Small mammal breeds. 2014.
- Gu Z, Ren W, Huang R, Huang Y, Chen B. Study on density of rex rabbit. Chinese Journal of Rabbit Farming.1999; 4:18-21.
- Millar SE. Molecular mechanisms regulating hair follicle development. Journal of Investigative Dermatology. 2002; 118:216–225.
- Fuchs E. Scratching the surface of skin development. Nature. 2007;445:834-842.
- Lin CM, Yuan Y, Chen X, Li H, Cai B, Liu Y, et al. Expression of Wnt/β–catenin signaling, stem-cell markers and proliferating cell markers in rat whisker hair follicles. Journal of Molecular Histology. 2015;46(3):233-240.
- Demehri S, Kopan R. Notch signaling in bulge stem cells is not required for selection of hair follicle fate. Development. 2009; 136(6): 891-896.
- Kulessa H, Turk G, Hogan BL. Inhibition of Bmp signaling affects growth and differentiation in the anagen hair follicle. Embo Journal. 2000;19(24):6664-6674.
- Su R, Li J, Zhang W, Yin J, Zhao J, Chang Z. Expression of BMP2 in the skin and hair follicle from different stage in Inner Mongolia cashmere goat. Scientia Agricultura Sinica. 2008; 23:559-563.
- Foitzik K, Lindner G, Mueller-Roever S, Maurer M, Botchkareva N, Botchkarev V, et al. Control of murine hair follicle regression (catagen) by TGF-beta1 in vivo. The FASEB Journal. 2000;14:752-760.
- Greco V, Chen T, Rendl M, Schober M, Pasolli HA, Stokes N, et al. Two-step mechanism for stem cell activation during hair regeneration. Cell Stem Cell. 2009; 4(2):155–169.
- Driskell RR, Clavel C, Rendl M, Watt FM. Hair follicle dermal papilla cells at a glance. J. Cell Sci. 2011;124:1179-1182.
- Chi W, Wu E, Morgan BA. Dermal papilla cell number specifies hair size, shape and cycling and its reduction causes follicular decline. Development. 2013; 140:1676-1683.
- Woo C, Eleanor W, Bruce A. Morgan permal papilla cell number specifies hair size, shape and cycling and its reduction causes follicular decline. Development and stem cells. 2013; 140: 1676-1683.
- Yi R, Fuchs E. MicroRNAs and their roles in mammalian stem cells. J. Cell Sci. 2011; 124: 1775–1783.
- Carrington JC, Ambros V. Role of microRNAs in plant and animal development. Science 2003; 301:336–338.
- Xu P, Vernooy SY, Guo M, Hay BA. The Drosophila microRNA Mir-14 suppresses cell death and is required for normal fat metabolism. Curr. Biol. 2003; 13: 790–795.
- Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell 2004; 116:281–297.
- Teleman AA., Maitra S, Cohen SM. Drosophila lacking microRNA miR-278 are defective in energy homeostasis. Genes Dev. 2006; 20: 417–422.
- He L, Lim LP, Stanchina E, Xuan Z, Liang Y, Xue W, et al. A microRNA component of the p53 tumour suppressor network. Nature. 2007; 447:1130–1134.
- Mardaryev AN, Ahmed MI, Vlahov NV, Fessing MY, Gill J H, Sharov AA, et al. Micro-RNA-31 controls hair cycle-associated changes in gene expression programs of the skin and hair follicle. FASEB Journal. 2010; 24 (10), 3869 -3881.
- 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): 50001.
- Wang P, Hong W, Zhen J, Ren J, Li Z, Xu A. The changes of microRNA expression profles and tyrosinase related proteins in MITF knocked down melanocytes. Molecular BioSystems. 2012; 8: 2924-2931.
- Ahmed MI, Alam M, Emelianov VU. MicroRNA214 controls skin and hair follicle development by modulating the activity of the Wnt pathway. Journal of Cell Biology. 2014; 207(4), 549-567.
- Andl T, Murchison EP, Liu F, Zhang Y, Yunta-Gonzalez M, Tobias JW, et al. The miRNA-processing enzyme dicer is essential for the morphogenesis and maintenance of hair follicles. Current Biology. 2006;16(10), 1041–1049.
- Shen Y, Zhang Y, Liu N, Wang H, Xie J, Gao S, et al. Let-7b promotes alpaca hair growth via transcriptional repression of TGFβRI. Gene. 2016; 577:32–36.
- Ahmed MI, Mardaryev AN, Lewis CJ, Sharov AA, Botchkareva NV. MicroRNA-21 is an important downstream component of BMP signaling in epidermal keratinocytes. Journal of Cell Science 2011; 124(20), 3399–3404.
- Gregory PA, Bert AG, Paterson EL, Barry SC, Tsykin A, Farshid G, et al. The miR-200 family and miR-205 regulate epithelial to mesenchymal transition by targeting ZEB1 and SIP1. Nat. Cell Biol. 2008; 10, 593–601.
- Yi R, O\"Carroll D, Hilda AZ, Zhihong D, Fred ST, Alexander F, et al. Morphogenesis in skin is governed by discrete sets of differentially expressed microRNAs. Nat. Genet. 2006; 38, 356–362.
- Ryan DG, Oliveira-Fernandes M, Lavker RM. MicroRNAs of the mammalian eye display distinct and overlapping tissue specificity. Mol. Vis. 2006; 12, 1175–1184.
- Wang D, Zhang Z, O'Loughlin E, Wang L, Fan X, Lai E. et al. MicroRNA-205 controls neonatal expansion of skin stem cells by modulating the PI(3)K pathway. Nat. Cell Biol. 2013;15, 1153–1163.
- 3 Alonso L, Fuchs E. The hair cycle. J Cell Sci 2006;119,391-3.
- 3 Galbraith H. Fundamental hair follicle biology and fine fibre production in animals. Animal 2010;4,1490-509.
- 3 Zhu B, Xu T, Yuan J, Guo X, Liu D. Transcriptome sequencing reveals differences between primary and secondary hair follicle derived dermal papilla cells of the Cashmere goat (Capra hircus). PloS one 2013;8, e76282.
- Wang ZP, Zhang H, Yang H, Wang SZ, Rong EG, Pei WY, Li H, Wang N. Genome-wide association study for wool production traits in a Chinese Merino sheep population. Plos One 2014;9, e107101.
- Liu NH, Li K, Liu J, Yu M, Cheng W, De J, Liu S, Shi Y, Zhao JS. Differential expression of genes and proteins associated with wool follicle cycling. Mol Biol Rep 2014; 41:5343-5349.
- Meale SJ, Chaves AV, Ding S, Bush RD, McAllister TA. Effects of crude glycerin supplementation on wool production, feeding behavior, and body condition of Merino ewes. J Anim Sci 2013; 91, 878-885.
- Wu Z, Sun L, Liu G, Liu H, Liu H, Yu Z, et al. Hair follicle development and related gene and protein expression of skins in Rex rabbits during the first 8 weeks of life. Asian-Australas J Anim Sci. 2019; 32(4):477-484.
- Wu M, Sun Q, Guo X, Liu H. hMSCs possess the potential to differentiate intoc DP cells in vivo and in vitro. Cell Biol Int Rep. 2010; 19 (2), e00019.
- Lim LP, Glasner ME, Yekta S, Burge CB, Bartel DP. Vertebrate microRNA genes. Science. 2003; 299, 1540.
- Wienholds E, Plasterk RH. MicroRNA function in animal development. FEBS Lett. 2005; 579, 5911–5922.
- Yi R., O\"Carroll D, Hilda AZ, Zhihong D, Fred ST, Alexander F, et al. Morphogenesis in skin is governed by discrete sets of differentially expressed microRNAs. Nat. Genet. 2006; 38:356–362.
- Ryan DG, Oliveira-Fernandes M, Lavker RM. MicroRNAs of the mammalian eye display distinct and overlapping tissue specificity. Mol. Vis. 2006; 12: 1175–1184.
- Wang D, Li Q, Feng N, Cheng G, Guan Z, Wang Y, et al. MicroRNA-205 controls neonatal expansion of skin stem cells by modulating the PI3K pathway. Nat. Cell Biol., 2013; 15:1153–1163.
- Yu J, Peng H, Ruan Q, Fatima A, Getsios S, Lavker RM, et al. MicroRNA-205 promotes keratinocyte migration via the lipid phosphatase SHIP2. FASEB J. 2010; 24:3950–3959.
- Lin D, Halilovic A, Yue P, Bellner L, Wang K, Wang L, et al. Inhibition of miR-205 impairs the wound-healing process in human corneal epithelial cells by targeting KIR4.1 (KCNJ10). Invest. Ophthalmol. Vis. Sci. 2013; 54:6167–6178.
- Du T, Zamore PD. Beginning to understand microRNA function. Cell Res. 2007; 17:661–663.
- Kollar EJ. The induction of hair follicles by embryonic dermal papillae. J. Invest. Dermatol. 1970; 55:374-378.
- Oliver RF. The induction of hair follicle formation in the adult hooded rat by vibrissa dermal papillae. J. Embryol. Exp. Morphol. 1970; 23:219-236.
- Jahoda CA, Horne KA, Oliver RF. Induction of hair growth by implantation of cultured dermal papilla cells. Nature1984; 311:560-562.
- Zhang J, He X, Tong W, Johnson T, Wiedemann L, Mishina Y, et al. Bone morphogenetic protein signaling inhibits hair follicle anagen induction by restricting epithelial stem/progenitor cell activation and expansion. Stem Cells 2006; 24:2826-2839.
- Liu G, Bai L, Li S, Liu H, Zhu Y, Sun H, et al. Isolation, culture and growth characteristics of dermal papilla cells from Rex rabbits. Tissue and Cell. 2020; 65,101348.
- Zha W, Zhou L, Li S, Liu K, Yang G, Chen Z, et al. Characterization and comparative profiling of the small RNA transcriptomes in the Hemipteran insect Nilaparvata lugens. Gene. 2016; 595:83-91.
- Wang Z, Gerstein M, Snyder M. RNA-Seq: a revolutionary tool for Transcriptomics. Nature Reviews Genetics. 2009; 10(1):57-63.
- Mortazavi A, Williams BA, McCue K, Schaeffer L. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nature Methods. 2008;5(7): 621-628.
- Eric PN, Sarah WB, Alex B, Jennifer D, Ruth YE, Sean RE, et al. Rfam 12.0: updates to the RNA families database. Nucleic Acids Research. 2014;43: D130-137.
- Kozomara A, Griffiths-Jones S. miRBase: annotating high confidence microRNAs using deep sequencing data. NAR. 2014;42: D68-D73.
- Friedlander MR, Chen W, Adamidic C, Maaskoia J, Einspanier R, Knespel S, et al. Discovering microRNAs from deep sequencing data using miRDeep. Nature Biotechnology 2008;26:407-415.
- Peter AC, Yavuz A, Heleneb HT, Erno V, Rolf H, Renée X, Judith M, et al. Deep sequencing-based expression analysis shows major advances in robustness, resolution and inter-lab portability over five microarray platforms. Nucleic Acids Res. 2008;36(21): e141.
- Jiang H, Wong WH. Statistical inferences for isoform expression in RNASeq. Bioinformatics 2009;25:1026-1032.
- Yang YH, Dudoit S, Luu P, Lin DM, Peng V, Ngai J. Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. Nucleic Acids Res. 2002;30: e15.
- Wang L, Feng Z, Wang X, Wang X, Zhang X. DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics. 2010;26(1):136-138.
- John B, Enright AJ, Aravin AA, Tuschl T, Sander C, Marks DS. Human MicroRNA targets. Plos. Biology. 2005;3(7), e264.
- Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008;36: D480-484.
- Liu GY, Wu ZY, Zhu YL, Liu L, Li FC. Effects of dietary vitamin B 6 on the skeletal muscle protein metabolism of growing rabbits, Animal Production Science, 2017, 57, 2007–2015.
- LiLivak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2 -ΔΔCT method. Methods. 2001; 25: 402–408.
Table 1. Effects of ocu-miR-205 on proliferation, cell cycle and apoptosis of dermal papilla cells (DPCs; %)
|
Items |
Group |
R-MSE |
P-value |
|||
|
Control |
HBAD-GFP |
HBAD-ocu-miR-205-GFP |
HBAD-ocu-miR-205-5p-sponge-GFP |
|||
|
Proliferation of dermal papilla cells |
||||||
|
Optical density (OD) value |
0.56±0.01b |
0.55±0.01b |
0.61±0.01a |
0.34±0.01c |
0.0356 |
<0.0001 |
|
Cell cycle of dermal papilla cells |
||||||
|
Resting state/first gap (G0/G1) |
87.79±1.00a |
84.68±1.01b |
77.71±0.88c |
85.57±1.03ab |
2.7795 |
<0.0001 |
|
Synthesis (S) |
7.15±0.85b |
9.48±0.95b |
14.73±0.62a |
8.97±0.68b |
2.2232 |
<0.0001 |
|
Second gap/mitosis (G2/M) |
5.07±0.24b |
5.37±0.49b |
7.56±0.37a |
5.46±0.35b |
1.0523 |
0.0002 |
|
Apoptosis of dermal papilla cells |
||||||
|
Early apoptotic ratio (Q4) |
33.95±0.40b |
32.86±1.04b |
36.71±0.67a |
29.16±0.47c |
1.9470 |
<0.0001 |
|
Later apoptotic ratio (Q2) |
33.46± 0.57b |
32.15±1.33b |
36.11±0.73a |
31.40±0.57b |
2.4323 |
0.0032 |
|
Total apoptosis ratio (Q4+Q2) |
67.41±0.51b |
65.01±1.50b |
72.83±1.14a |
60.56±0.47 b |
2.8310 |
<0.0001 |
Note: Data shown are mean values ± s. d., and n = 8 per group. In the same row, values with same letter superscripts mean no significant difference (P>0.05), with different letter superscripts mean significant difference (P<0.05).
Table 2. Effects of ocu-miR-205 on gene expression of signal pathway of dermal papilla cells (DPCs)
|
Gene |
Group |
R-MSE |
P-value |
|||||||||||
|
Control |
HBAD-GFP |
HBAD-ocu-miR-205-GFP |
HBAD-ocu-miR-205-5p-sponge-GFP |
|||||||||||
|
PI3K/Akt signal pathway |
||||||||||||||
|
Inppl1 |
1.00±0.09b |
1.02±0.04b |
0.80±0.07b |
1.68±0.12a |
0.2353 |
<0.0001 |
||||||||
|
Inpp4b |
1.00± 0.09 |
0.97±0.16 |
0.92±0.08 |
1.21±0.16 |
0.3598 |
0.3966 |
||||||||
|
Frk |
1.00±0.34b |
0.90±0.07b |
0.50±0.13c |
2.63±0.47a |
0.8442 |
0.0001 |
||||||||
|
Phlda3 |
1.00± 0.04b |
0.89±0.06b |
0.69±0.05c |
1.18±0.05a |
0.1491 |
<0.0001 |
||||||||
|
Wnt signal pathway |
||||||||||||||
|
Wnt10b |
1.00±0.11b |
0.95±0.09b |
0.38±0.11c |
2.24±0.17a |
0.3519 |
<0.0001 |
|
|||||||
|
CTNNB1 |
1.00±0.09b |
1.50±0.13b |
1.48±0.19b |
3.16±0.27a |
0.5186 |
<0.0001 |
||||||||
|
GSK-3β |
1.00± 0.14b |
0.97±0.05b |
2.66±0.68a |
3.22±0.71a |
1.4118 |
0.0050 |
|
|||||||
|
DKK1 |
1.00±0.09 c |
2.90±0.13b |
6.09±0.84 a |
0.84±0.07 c |
1.2141 |
<0.0001 |
||||||||
|
Notch signal pathway |
||||||||||||||
|
Notch1 |
1.00±0.08b |
0.96±0.07b |
0.39±0.02c |
1.19±0.03a |
0.1604 |
<0.0001 |
||||||||
|
Jagged1 |
1.00±0.07 ab |
0.81±0.13 bc |
0.60±0.10c |
1.12±0.07 a |
0.2679 |
0.0032 |
||||||||
|
Hes1 |
1.00±0.12b |
1.19±0.16b |
0.90±0.05b |
2.40±0.33a |
0.5475 |
<0.0001 |
||||||||
|
Hes5 |
1.00± 0.16 b |
1.04±0.06 b |
0.18±0.04c |
1.38±0.04a |
0.2477 |
<0.0001 |
||||||||
|
BMP signal pathway |
||||||||||||||
|
BMP2 |
1.00±0.15b |
0.51±0.05bc |
1.69±0.34a |
0.41±0.05c |
0.5408 |
0.0002 |
||||||||
|
BMP4 |
1.00± 0.20 ab |
1.11±0.18 ab |
1.34±0.08a |
0.63±0.13b |
0.4405 |
0.0254 |
||||||||
|
TGF-β1 |
1.00±0.05 ab |
0.89±0.18 b |
1.33±0.13a |
0.71±0.07 b |
0.3304 |
0.0071 |
||||||||
Note: Data shown are mean values ± s. d., and n = 8 per group. In the same row, values with no letter superscripts or same letter superscripts mean no significant difference (P>0.05), with different letter superscripts mean significant difference(P<0.05).
Table 3. Effects of ocu-miR-205 on gene expression of signal pathway of Rex rabbits skin
|
Gene |
Group |
R-MSE |
P-value |
||||||
|
Control |
HBAD-GFP |
HBAD-ocu-miR-205-GFP |
HBAD-ocu-miR-205-5p-sponge-GFP |
||||||
|
PI3K/Akt signal pathway |
|||||||||
|
Inppl1 |
1.00±0.09b |
0.99±0.07b |
0.54±0.03c |
2.30±0.15a |
0.2731 |
<0.0001 |
|||
|
Inpp4b |
1.00± 0.15ab |
0.89± 0.07a |
0.09±0.01b |
1.05±0.02a |
0.2348 |
<0.0001 |
|||
|
Frk |
1.00±0.26a |
1.16±0.05b |
0.31±0.02c |
1.94±0.20a |
0.4666 |
<0.0001 |
|||
|
Phlda3 |
1.00± 0.11b |
1.18± 0.11a |
0.27±0.03b |
1.18±0.11a |
0.2500 |
<0.0001 |
|||
|
Wnt signal pathway |
|||||||||
|
Wnt10b |
1.00±0.14b |
0.86±0.04bc |
0.74±0.04c |
1.26±0.09a |
0.2406 |
0.0011 |
|||
|
CTNNB1 |
1.00±0.07b |
0.98±0.05b |
0.75±0.04c |
1.30±0.08a |
0.1789 |
<0.0001 |
|||
|
GSK-3β |
1.00± 0.06ba |
1.04±0.07b |
0.85±0.04b |
1.10±0.05a |
0.1496 |
0.0198 |
|||
|
DKK1 |
1.00±0.36a |
0.91±0.10a |
1.28±0.10a |
0.18±0.07b |
0.5625 |
0.0041 |
|||
|
Notch signal pathway |
|||||||||
|
Notch1 |
1.00±0.18b |
0.83±0.05bc |
0.54±0.05c |
1.95±0.22a |
0.4131 |
<0.0001 |
|||
|
Jagged1 |
1.00±0.13a |
1.02±0.05 a |
0.36±0.03c |
1.34±0.22 a |
0.3670 |
0.0001 |
|||
|
Hes1 |
1.00±0.09b |
1.03±0.06b |
0.80±0.03b |
1.49±0.18a |
0.3042 |
0.0009 |
|||
|
Hes5 |
1.00± 0.13 b |
1.00±0.05b |
0.36±0.03c |
1.43±0.20a |
0.3435 |
<0.0001 |
|||
|
BMP signal pathway |
|||||||||
|
BMP2 |
1.00±0.15b |
1.02±0.07b |
1.39±0.08a |
0.64±0.10c |
0.3051 |
0.0005 |
|||
|
BMP4 |
1.00± 0.15b |
0.87±0.04b |
1.58±0.14a |
0.73±0.05b |
0.3081 |
<0.0001 |
|||
|
TGF-β1 |
1.00±0.12 ab |
0.88±0.05 bc |
1.15±0.06a |
0.67±0.05b |
0.2222 |
0.0014 |
|||
Note: Data shown are mean values ± s. d., and n = 8 per group. In the same row, values with same letter superscripts mean no significant difference (P>0.05), with different letter superscripts mean significant difference (P<0.05).
Table 4. Effects of ocu-miR-205 on hair follicle density of Rex rabbits (Count/mm2)
|
Items |
Group |
R-MSE |
P-value |
|
|||||||||
|
Control |
HBAD-GFP |
HBAD-ocu-miR-205-GFP |
HBAD-ocu-miR-205-5p-sponge-GFP |
|
|||||||||
|
Transfection for 7 days |
|
|
|
|
|
|
|
||||||
|
Hair follicle density |
125.57±9.63bc |
144.25±11.47ba |
105.02±5.89c |
163.94±16.68a |
32.7817 |
0.0085 |
|
||||||
|
Primary hair follicle density |
10.08±0.85 |
9.31±0.85 |
7.80±0.63 |
9.79±0.77 |
2.2046 |
0.1902 |
|
||||||
|
Secondary hair follicle density |
115.49±9.48bc |
134.93±12.01ba |
97.22±5.68c |
154.15±16.65a |
32.9728 |
0.0113 |
|
||||||
|
Secondary/Primary (S/P) ratio |
11.92±1.17 |
16.02±2.67 |
12.92±1.06 |
16.38±2.04 |
5.2469 |
0.2536 |
|
||||||
|
Transfection for 14 days |
|
|
|
|
|
|
|
||||||
|
Hair follicle density |
130.54±12.22b |
136.62±8.51b |
119.30±7.24b |
203.34±14.79a |
31.4005 |
<0.0001 |
|
||||||
|
Primary hair follicle density |
10.86±0.54 |
9.96±0.50 |
11.92±0.74 |
10.11±0.70 |
1.7739 |
0.1307 |
|
||||||
|
Secondary hair follicle density |
119.68±12.03b |
126.66±8.28b |
107.37±7.10b |
213.32±14.62a |
30.9035 |
0.0001 |
|
||||||
|
Secondary/Primary (S/P) ratio |
11.08±0.96 b |
12.82±0.82 b |
9.20±0.78 b |
21.62±1.94 a |
3.4505 |
<0.0001 |
|
||||||
|
Transfection for 21 days |
|
|
|
|
|
|
|
||||||
|
Hair follicle density |
145.43±13.71b |
139.58±16.53b |
126.38±9.30b |
208.69±21.42a |
44.8651 |
0.0049 |
|||||||
|
Primary hair follicle density |
14.30±1.19 |
12.65±0.80 |
12.92±0.51 |
15.29±1.28 |
2.8084 |
0.2244 |
|||||||
|
Secondary hair follicle density |
131.13±13.21b |
126.94±16.57b |
113.46±9.27b |
193.40±21.31a |
44.4697 |
0.0060 |
|||||||
|
Secondary/Primary (S/P) ratio |
9.45±0.91 |
10.28±1.40 |
8.86±0.75 |
13.28±1.85 |
3.6808 |
0.1018 |
|||||||
Note: Data shown are mean values ± s. d., and n = 8 per group. In the same row, values with no letter superscripts or same letter superscripts mean no significant difference (P>0.05), with different letter superscripts mean significant difference (P<0.05).
Additional file 1:
Supplementary Figure 1. Skin crosscutting HE staining of Rex rabbits with different hair density（100 magnification）(a) Low hair density (b) High hair density
Supplementary Table 1. Statistics of hair follicle numbers of Rex rabbits with different hair density
Note: Data shown are mean values ± s. d., and n =3 per group. In the same row, values with with different letter superscripts mean significant difference(P<0.05).
Additional file 2:
Supplementary Table 2. RNA quality of six DPC samples
Supplementary Figure 2. RNA quality (a) Agarose electrophoresis of RNA in DPCs; (b) Quality test results of total RNA by Agilent 2100
Additional file 3:
Supplementary Table 3. Statistics of sequence data of each samples
Additional file 4:
Supplementary Figure 3. Length distribution map of Small RNA
The X-axis is the length of small RNA, and the Y-axis is the corresponding number of small RNA
Additional file 5:
Supplementary Table 4. Alignment clean tag in genome
Additional file 6:
Supplementary Figure 4. Statistical distribution map of small RNA types
In order to make each unique small RNA have a unique annotation, the annotation statistics of small RNA traverse the annotation according to the priority order of miRNA > piRNA > snoRNA > Rfam > other sRNA.
Additional file 7:
Supplementary Figure 5. The quantity distribution maps of clean tag base in each sample
The X-axis is the position of the base in read, and the Y-axis represents the proportion of the base.
Additional file 8:
Supplementary Figure 6. The quality distribution maps of the clean tag base in each sample
The X-axis is the position of base in read, and the Y-axis represents the base mass value.
Additional file 9:
Supplementary Table 5. LD-vs-HD_DEGseq. diffexpfilter
Additional file 10:
Supplementary Figure 7. Directed acyclic graph of GO enrichment
(a) LD-vs-HD_DEGseq. Biological_Process. Top GO；LD-vs-HD_DEGseq. Cellular_Component. Top GO；LD-vs-HD_DEGseq. Molecular_Function. Top GO.
The branches in the figure represent the inclusion relationship, and the function range defined from top to bottom is getting smaller and smaller. The box represents the top 5 GO terms of each classification enrichment degree, and through the inclusion relationship, the associated GO terms are displayed together. The name of the term and the P-value corrected by enrichment analysis are displayed on each node. The darker the color (red) indicates the smaller the P-value and the higher the enrichment degree.
Additional file 11:
Supplementary Table 6. Primer sequence information in experiments