About the Authors:
Wen Gao
Affiliation: College of Animal Science and Technology, Yangzhou University, Yangzhou, China
Wei Sun
* E-mail: [email protected]
Affiliation: College of Animal Science and Technology, Yangzhou University, Yangzhou, China
Jinfeng Yin
Affiliation: College of Animal Science and Technology, Yangzhou University, Yangzhou, China
Xiaoyang Lv
Affiliation: College of Animal Science and Technology, Yangzhou University, Yangzhou, China
Jianjun Bao
Affiliation: College of Animal Science and Technology, Yangzhou University, Yangzhou, China
Jiarui Yu
Affiliation: College of Animal Science and Technology, Yangzhou University, Yangzhou, China
Lihong Wang
Affiliation: College of Animal Science and Technology, Yangzhou University, Yangzhou, China
Chengyan Jin
Affiliation: College of Animal Science and Technology, Yangzhou University, Yangzhou, China
Liang Hu
Affiliation: College of Animal Science and Technology, Yangzhou University, Yangzhou, China
Introduction
Hu sheep is a unique breed of white lambskin sheep that has been domesticated in China for about 800 years. Hu sheep have several excellent characters such as good adaptability, high prolificacy and reproduction, early maturation, and high lactation. Furthermore, the value of Hu sheep lambskin in the international market is relatively high and this is mainly due to its superior quality, spotless white color, silky smooth texture, and peculiar flower patterns that are similar to water waves. Based on the width of the flower pattern, Hu sheep lambskin can be divided into three wave types, namely, large, medium, and small, with the small waves showing the best quality [1]. Several factors determine the type of lambskin pattern, which include fineness, density, and curvature of wool. The hair follicle is an accessory organ of the skin that has complex morphology and structure. Hair follicles control the growth of hair and can be divided into primary follicle and secondary follicle, which respectively grow myelinated hair and unmyelinated hair [2]. Therefore, studying the inheritance pattern of hair follicles may reveal the regulatory mechanism of flower pattern formation in Hu sheep lambskin at the molecular level, which in turn may be applied in breeding schemes for the improvement of Hu sheep lambskin quality.
Current research studies on miRNAs mainly focus on muscle, liver, ovary and cancer tissues, whereas most of studies on sheep hair follicles concentrate on its growth cycle and only a few studies centered on miRNA expression and hair follicle development in Hu sheep lambskin. He et al. [3] extracted the total RNA from alpaca skin, prepared miRNA chips and screened miRNAs in the ear and back skin by using Affymetrix miRNA chip interspecies hybridization. Andl et al. [4] reported that Dicer mRNA and miRNA were expressed in the skin tissue of mice, indicating that miRNA plays an important role in the development of skin and hair follicle. Zhu et al. [5] found that the key gene of chromogenesis, Mitf, could be inhibited by transducing the lentiviral vector of miR-25 into the melanophores of alpaca skin, prompting them to hypothesize that miR-25 inhibits Mitf expression, and miR-25 possibly participates in the regulation of hair color. Mardaryev et al. [6] found that the expression of miR-31 significantly increased during skin and hair follicle growth. Furthermore, they found that miR-31 plays a significant role in regulating hair growth and hair follicle development in mice through bone morphogenetic protein (BMP) and Wingless-Int (WNT) signaling pathways. In addition, Tang et al. [7] found that miR-31 is highly expressed during hair follicle growth and downregulated in telogen.
To date, most investigations have focused on polyembryony and meat performance instead of lambskin quality, thereby hindering further improvements in breeding Hu sheep with higher quality lambskin and in protecting germplasm resources. A number of histological assessment of Hu sheep lambskin have been performed; however, our understanding of the specific molecular mechanism underlying the growth and development of hair follicles and the formation of different flower patterns in Hu sheep lambskin is limited. In the present study, high-throughput sequencing and bioinformatics analysis were initially employed to identify differentially expressed miRNAs that influence the growth and development of hair follicles as well as the formation of specific flower patterns in Hu sheep lambskin. By combining histological observation and the micro-observation technology, the correlation between the expression of 14 miRNAs in different wave patterns and the histological properties of hair follicles was analyzed. Our results showed that seven candidate miRNAs, namely, miRNA-143, miRNA-10a, let-7i, NW_004080184.1_6326, NW_004080165.1_8572, NW_004080181.1_3961, and NW_004080190.1_13733, are involved in the development of Hu sheep lambskin hair follicles, which in turn may facilitate in the elucidation of the molecular mechanism underlying its growth and development in Hu sheep lambskin.
Materials and methods
Ethics statement
This study was conducted in strict compliance with the recommendations of the Guide for the Care and Use of Laboratory Animals of Jiangsu Province and of the Animal Care and Use Committee of the Chinese Ministry of Agriculture. The government of Jiangsu Province (Permit Number 45) and the Ministry of Agriculture of China (Permit Number 39) approved the protocol performed in this study. All efforts were made to minimize animal suffering.
Experimental animals
Six pairs of full-sib Hu sheep were selected at birth at the Suzhou stud farm in China. Each pair consisted of one individual with predominantly large-wave wool, one with medium-wave wool, and one with predominantly small-wave wool. To extract RNA from the hair follicles, about 1 cm of the hair root was excised and placed into a microtube that contained with 1 mL of TRIzol® Reagent (Invitrogen, Carlsbad, CA, USA) and surrounded by drikold. The quantity of hair root was collected up to a third of the volume of the microtube. Each individual was locally anaesthetized by subcutaneous injection of 0.1 mL of 2% procaine (Changzhou Sunchem Pharmaceutical Chemical Material Co., Ltd, Changzhou, China) prior to the removal of the back skin tissue at an approximate size of 1.5 cm2. Then, the back skin tissue was flattened on the cardboard, which was then immersed in 4% paraformaldehyde solution (Sinopharm Chemical Reagent Co., Ltd, Shanghai, China) in a centrifuge tube.
Total RNA extraction and small RNA library preparation
Total RNA was extracted from the hair roots of six pairs of full-sib individuals using a buffer containing ß-mercaptoethanol and guanidine that was provided in an RNAiso Plus kit in accordance with the manufacturer’s instructions (Takara Biotechnology Dalian, Co. Ltd., China). RNA was eluted with 40 uL of RNase-free water, and the quantity and purity of total RNA were measured by using a NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies, Wilmington, DE, U.S.). Samples with purity (A260/A280) of 1.8–2.0 were used in the subsequent experiment. Total RNA was purified using an RNAclean Kit (TIANGEN BIOTECH Co., Ltd, Beijing, China). The samples for small RNA sequencing were prepared in accordance with the instructions that were provided by the manufacturer of the Illumina sequencing kit.
Small RNA sequencing and processing of reads
Small RNA libraries were sequenced on an Illumina HiSeq 2500 system with 50-base pair (bp) single-end reads according to the manufacturer’s protocol. Then, a series of treatments, including quality inspection of bases and selection of length, as well as sequencing fragments of reliable quality were performed, which was followed by estimation of varieties and quantities by using a distribution statistic. The length range of small RNAs is 18–30 nt and that of miRNAs is 21–22 nt, the peak of the length distribution can help to determine small RNA species. All small sequencing data have been uploaded to the Sequence Read Archive database (SRR3099013, SRR3099014, and SRR3099015).
Validation of small RNA sequencing by real-time PCR
Differentially expressed miRNAs, namely, miRNA-10a, miR-199a-3p, let-7i, NW_004080190.1_13733, and NW_004080166.1_9139 were randomly selected for real-time PCR, and the U6 gene was used as reference (S1 Table). The SYBR Green I method was used for quantitative testing to verify the reliability of the sequencing results. A standard curve was established using a cDNA gradient dilution, and three replicates were prepared for each sample. The relative expression of the target gene was calculated by using the following formula: Relative expression = 2-ΔΔCt, ΔΔCt = (Ctarget gene—Chousekeeping gene)(large waves, medium waves)—(Ctarget gene—Chousekeeping gene)small waves.
Paraffin section fabrication
The lambskin tissues were treated according to the conventional paraffin section method. After dehydration, transparent, embedding, sectioning, and HE staining, the structure of the hair follicles was observed under a microscope and representative images were captured. Three different fields of view of each section were selected. The number and diameter of the primary and secondary hair follicles were respectively determined. Each hair follicle feature was compared among the large, medium, and small waves by using ANOVA.
Quantitation of 14 miRNAs by real-time PCR
Total RNA was extracted by using a buffer that contained β-mercaptoethanol and guanidine, which was provided in an RNAiso Plus Kit, following the manufacturer’s instructions (Takara Biotechnology Dalian, Co., Ltd.). RNA was eluted with 40 μL of RNase-free water. The concentration of RNA was measured by using a NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA). Samples with purity (A260/A280) of > 1.8 were used. A total of 2 μg of RNA from each sample was transcribed into cDNA by using a TaKaRa reserve transcription kit (TaKaRa Biotechnology Dalian, Co., Ltd.) according to manufacturer’s instructions. The primers used in real-time PCR are listed in S1 Table.
The U6 gene was employed as reference, and the relative expression level of the target gene was calculated by using 2-ΔΔCt method. According to a TIANGEN miRcute miRNA fluorescence quantitative detection kit, the optional reaction system was a 20-μL reaction containing 10 μL of 2× miRcute miRNA Premix, 0.4 μL of the forward primer, 0.4 μL of the reverse primer, 1 μL of the first-strand cDNA of the miRNA, and 8.2 μL of ddH2O. The PCR conditions included 40 cycles of 2 min at 94°C, 20 s at 94°C, and 34 s at 60°C.
Statistical analysis
SPSS17.0 was used to calculate the Ct values and standard errors among repeat samples, and the difference in relative gene expression levels was analyzed by using the 2-ΔΔCt method. The relative expression of the target gene was calculated by using the following formula: Relative expression = 2-ΔΔCt, ΔΔC t = (Ctarget gene—Chousekeeping gene)(large waves, medium waves)−(Ctarget gene—Chousekeeping gene)(small waves). The expression levels of miRNA in different wave patterns were compared using ANOVA.
Results
miRNA sequencing and assembly
Three kinds of waves (large, medium, and small) were selected for analysis. The miRNA library was established using a Illumina Hiseq 2500 system, and sequencing generated 7,964,087, 7,167,843, and 7,534,406 miRNAs from large, medium, and small wave wool, respectively. Detailed results of the sequencing and assembly are shown in Table 1.
[Figure omitted. See PDF.]
Table 1. Summary of Illumina sequencing and miRNA assembly.
https://doi.org/10.1371/journal.pone.0176532.t001
Differential expression analysis of miRNAs
A total of 522 differentially expressed miRNAs were screened using high-throughput sequencing. Of these miRNAs, 384 were co-expressed in the three patterns, whereas 9, 26, and 40 miRNAs were expressed in small waves, medium waves, and large waves, respectively (Fig 1). Statistical analysis was employed to identify significant differentially expressed miRNAs in the two kinds of patterns. The results showed that there were 36 differentially expressed miRNAs, including 29 novel miRNAs. Among these miRNAs, 10, 10, and 5 miRNAs were differentially expressed in small and large, small and medium, and medium and large waves, respectively, whereas no differentially expressed miRNAs in three kinds of waves were detected.
[Figure omitted. See PDF.]
Fig 1. The distribution of miRNA in small waves, medium waves and large waves.
https://doi.org/10.1371/journal.pone.0176532.g001
Predictive analysis of target genes of differentially expressed miRNAs
The target genes of 522 differentially expressed miRNAs were predicted by using the miRanda predict software to annotate their biological functions. The result showed that there were 1,008, 2,062, and 4,169 target genes in small and large, medium and large, and small and medium waves, respectively (S2–S4 Tables).
Functional analysis of target genes
Gene ontology (GO) and pathway analyses were performed for these target genes. GO analysis identified 306 GO terms in small and large waves (Table 2). There were 217 GO terms in small and medium waves (Table 3). Finally, there were 371 GO terms in medium and large waves (Table 4). These were mainly involved in cell process, metabolism, ion transport, growth, and development. Pathway analysis showed that the target genes were related to 53, 61, and 65 signal pathways in small and large, medium and large, small waves and medium waves, respectively (Tables 5–7), including metabolic pathways, ion transport, cell proliferation and differentiation, and cell adhesion.
[Figure omitted. See PDF.]
Table 2. The top 10 enriched GO terms of target genes for differentially expressed miRNAs between small and large waves.
https://doi.org/10.1371/journal.pone.0176532.t002
[Figure omitted. See PDF.]
Table 3. The top 10 enriched GO terms of target genes for differentially expressed miRNAs between small and medium waves.
https://doi.org/10.1371/journal.pone.0176532.t003
[Figure omitted. See PDF.]
Table 4. The top 10 enriched GO terms of target genes for differentially expressed miRNAs between medium and large waves.
https://doi.org/10.1371/journal.pone.0176532.t004
[Figure omitted. See PDF.]
Table 5. The top 10 enriched KEGG pathway of target genes for differentially expressed miRNAs between small and large waves.
https://doi.org/10.1371/journal.pone.0176532.t005
[Figure omitted. See PDF.]
Table 6. The top 10 enriched KEGG pathway of target genes for differentially expressed miRNAs between small and medium waves.
https://doi.org/10.1371/journal.pone.0176532.t006
[Figure omitted. See PDF.]
Table 7. The top 10 enriched KEGG pathway of target genes for differentially expressed miRNAs between medium and large waves.
https://doi.org/10.1371/journal.pone.0176532.t007
Validation of small RNA sequencing by RT-PCR
To confirm the reliability of sequencing results, five miRNAs were selected for RT-PCR verification. The sequencing results were basically identical to that of RT-PCR of differentially expressed miRNAs (Table 8).
[Figure omitted. See PDF.]
Table 8. The comparison on the results of sequencing and RT-PCR for differentially expressed miRNAs.
https://doi.org/10.1371/journal.pone.0176532.t008
Histological properties of different patterns of lambskin hair follicles
Cross-sectional images of skin tissue sections of large, medium, and small wave lambskins of Hu sheep is shown in Fig 2. The statistical results of each indicator of hair follicle are shown in Table 9.
[Figure omitted. See PDF.]
Fig 2. Transverse section images of hair follicles in large, medium and small waves.
https://doi.org/10.1371/journal.pone.0176532.g002
[Figure omitted. See PDF.]
Table 9. Histological properties of hair follicles in large, medium and small waves.
https://doi.org/10.1371/journal.pone.0176532.t009
The diameter of the primary follicles was larger than that of secondary follicles regardless of large, medium, or small waves. The distribution of hair follicles is mainly in the form of hair follicle groups, centering on primary hair follicles, and secondary hair follicles are located adjacent to the primary hair follicles (Fig 2).
In the same field of view, no significant differences in primary follicles among the large, medium and small wave lambskins were observed (P>0.05); the number of secondary follicles in the small waves was significantly higher than that in large and medium waves (P<0.01), whereas the number of secondary follicles in the medium wave was similar to that in the small waves (P>0.05). Consequently, the total number of hair follicles in the small waves was higher than that in the large and medium waves in the same unit area. The diameters of the primary and secondary follicles in the large waves were significantly larger than those in the medium and small waves (P<0.01), whereas no significant difference was observed between the medium and small waves (P>0.05) (Table 9).
Expression of the 14 miRNAs in large, medium, and small waves in lambskin hair follicles
The level of expression of miRNA-143 in the large waves was significantly lower than that in the large waves (P<0.01). The expression level of miRNA-10a in the small waves was significantly higher than that in the large and medium waves (P<0.05) and that of let-7i was significantly higher in the medium waves compared to that in the large waves (P<0.05). In the differentially expressed novel miRNAs, the expression level of NW_004080184.1_6326 in medium waves was significantly lower than those in the large and small waves, respectively (P<0.05; P<0.01). The expression of NW_004080165.1_8572 in the medium waves was significantly different compared to that in the large and small waves (P<0.01). NW_004080181.1_3961 was expressed at a higher level in small waves than in medium waves (P<0.01). Highly significant differences in NW_004080190.1_13733 expression was observed between medium and large waves (P<0.01). No significant difference in the expression levels of the rest of the miRNAs were observed among the large, medium, and small waves (Table 10).
[Figure omitted. See PDF.]
Table 10. Relative expression quantities of the 14 miRNAs in large, medium and small waves.
https://doi.org/10.1371/journal.pone.0176532.t010
Correlation analysis between 14 miRNAs expression and histological properties of hair follicles
The association between the expression of 14 miRNAs and each histological properties of hair follicle was analyzed by using SPSS 17.0 (Tables 11–12), which showed that these were involved in the growth and development of hair follicles.
[Figure omitted. See PDF.]
Table 11. Correlation between 14 miRNAs expressions and diameters of hair follicle.
https://doi.org/10.1371/journal.pone.0176532.t011
[Figure omitted. See PDF.]
Table 12. Correlation between 14 miRNAs expressions and numbers of hair follicle.
https://doi.org/10.1371/journal.pone.0176532.t012
Discussion
MiRNAs function by interacting with its target genes. The present study determined that miRNAs could regulate multiple target genes and each target gene could be regulated by multiple miRNAs. The relationship between miRNAs and target genes plays an important role in the regulation of hair follicle growth and pattern formation; therefore, it is of great significance to explore the regulatory mechanisms underlying the growth and development of hair follicles in Hu sheep. Bioinformatics methods are currently employed to predict target genes based on the fact that miRNAs can combine with its target genes; however, some incomplete complementary sequence may also be paired, so there must be some false positives in the generated results.
Andl et al.[4] and Yi et al. [8] demonstrated the main role of miRNAs in mice was to maintain skin hair follicle morphogenesis and development in mice, and this discovery has paved the way for further investigations on the role of miRNAs in hair follicles. Currently, miRNA expression profiling in skin and hair or feather follicles has been employed to identify miRNAs that are related to hair follicle development and growth, particularly in sheep, goat, and duck [9–11]. Liu et al. [12] detected that 316 conserved miRNAs and 22 novel miRNAs in the skin of cashmere goats by using Illumina’s Solexa sequencing technology. Wu et al. [13] identified 205 conserved miRNAs and 9 novel miRNAs by RNA-Seq in goat hair follicles. In the present study, we found 522 differentially expressed miRNAs in the hair follicles of different patterns. These differences in results may various among breeds and reference genomes.
Hair follicle development and its postnatal regeneration are characterized by dramatic changes in its microanatomy and cellular activity, which are controlled by multiple signaling pathways, transcription factors, and epigenetic regulators, including miRNAs [14]. Some researches that are based on gene regulatory networks of hair follicles showed that the Wnt [15], transforming growth factor-β (TGF-β) [16], mitogen-activated protein kinase (MAPK) [17], sonic hedgehog (Shh) [18], BMP [19], Notch, and JAK-STAT [20] signaling pathways widely participate in hair follicle development, morphogenesis, and hair follicle cycling. The Wnt signaling pathway stimulates hair follicle stem cells to transformation from telogen to growth. Yang et al. [21] showed that miRNA-148b activates the Wnt/nt/d via miRNA-148bs that activate follicle stem cells to promote the transformation of hair follicles. In most situations, MAPK signaling pathway usually plays a major role as a performer, and TGF-β signaling pathway is a leader to activate other signaling pathways, and Shh, Notch, and JAK-STAT signaling pathways play a fine-tuning role. Pathway analysis has indicated that the target genes are involved with 53, 65, and 61 signal pathways in small and large, medium and large, small and medium waves, respectively, including the metabolic, cancer, MAPK signaling, and PI3K-Akt signaling pathways. These pathways involve cell proliferation and differentiation, cell growth, and development, and thus these differentially expressed miRNAs and target genes may participate in hair follicle development.
The present study screened 36 differentially expressed miRNAs, including 29 novel miRNAs and 7 known miRNAs, although no differentially expressed miRNAs in the three patterns were observed. In these miRNAs, miRNA-143 was upregulated in hair follicles and is differentially expressed in small and large waves. No studies on the function of miRNA-143 in hair follicles have been conducted, although the results of the present study suggest that it might play an important role in the growth and development of hair follicles. Previous studies have shown that miRNA-143 plays an important role various malignancies [22–23]. miRNA-143 has been reported to be downregulated in gastric, lung, cervical, and bladder cancer. The upregulation of miR-143 in these malignancies inhibits cell growth, which in turn indicates that it participates in cell differentiation, proliferation, apoptosis, and metabolism [24–28]. In the present study, we found that LIM-only protein 4 (LMO4) was the direct target gene of miRNA-143 and NW_004080184.1_632. Lu et al. [29] reported that LMO4 interacts with Smad proteins and increases the effect of TGFβ on epithelial cell growth.What’s more, it was found that LMO4 was upregulated in epithelial cells, particularly during active mesenchymal-epithelial interactions, including the epidermis, hair follicles, and mammary gland [30–33]. let-7c and let-7i belong to the let-7 family, which is an important tumor-suppressor miRNA, and both are downregulated. These regulate cell cycle-related genes and participate in cell apoptosis [34–35]. miR-199a-3p is involved in cell proliferation and development. In liver cancer, miR-199a-3p inhibits the growth of hepatoma cells [36], and in testicular cancer, it reduces tumor cell proliferation, migration, and infection [37]. In addition, miR-199a-3p also participates in cell differentiation of C2C12 myoblasts. miR-199a-3p is the most widely reported in AKT/mTOR signaling pathway, which plays an important role in cell signaling and regulates various biological processes such as cell proliferation, metabolism, growth, differentiation, and apoptosis. Shen et al. [38] showed that miR-199a-3p suppresses the AKT/mTOR signaling pathway and cell proliferation in vitro. It has also been found that miR-10a participates in cell proliferation and migration, whereas in contrast to miR-199a-3p, miR-10a promotes tumor cell proliferation, invasion, and migration [39–40]. Fibroblast growth factor-1 (FGF1), the target gene of NW_004080189.1_7961, is a member of the fibroblast growth factor (FGF) family. FGF1 is involved in the growth and development of various tissues and organs and can significantly promote angiogenesis and cell division. FGF1 also plays a significant role in cell proliferation, migration, and differentiation [41–42]. Previous studies have shown that MT4 promotes tumor growth and metastasis and influences cell proliferation. Pave et al. [43] reported that MT4 promotes cancer cell proliferation in vivo and in vitro by promoting an outside-in signaling through the EGFR pathway. Based on this information, the differentially miRNAs evaluated in the present study may be involved in cell proliferation, differentiation, and apoptosis, as well as the growth and development of hair follicles in Hu sheep lambskin. In addition, nine highly and differentially expressed novel miRNAs (NW_004080189.1_7961, NW_004080184.1_6535, NW_004080184.1_6326, NW_004080165.1_8572, NW_004080166.1_9139, NW_004080181.1_3961, NW_004080174.1_726, NW_004080166.1_10417, and NW_004080190.1_13733) were screened and may also play important roles in the growth and development of hair follicles of Hu sheep.
To study the effect of 14 differentially expressed miRNAs on different wave patterns in Hu sheep lambskin, association analysis was employed to examine the miRNA expression in different wave patterns and histological properties of hair follicle such as the diameter and number of primary and secondary follicles. Seven miRNAs were differentially expressed in different wave patterns, namely, miRNA-143, miRNA-10a, and let-7i, and four novel miRNAs, including NW_004080184.1_6326, NW_004080165.1_8572, NW_004080181.1_3961, and NW_004080190.1_13733. Correlation analysis indicated that these 7 miRNAs could be considered as candidate miRNAs for the growth and development of hair follicles in Hu sheep lambskin. In the present study, miRNA-143 was one of the upregulated miRNAs in Hu sheep hair follicular tissues. Additionally, the results showed that the expression level of miRNA-143 in small waves was higher than that in large and medium waves at the same period, and highly significant differences between large and small wave lambskin hair follicles were observed. Consequently, we speculated that miRNA-143 regulates sheep hair growth during hair follicle cell development. Research studies on miRNA -10a, let-7c, and let-7i are mainly related to cancer and their roles in animal hair follicle development are largely unknown. The role of miRNA-10a in tumors is to promote tumor cell proliferation and differentiation [44–45]. Conversely, as an important tumor suppressor, the miRNAs in let7 families are downregulated in tumor cells and are involved in apoptosis [34,35]. Members of the let7 family regulate the formation of skin melanoma [46]. The essential functions of members of the let7 family in goat and sheep skin have been previously investigated [47]. Expressed in both anagen and telogen of hair follicles in cashmere goat, miRNA-1et7a plays a major role in regulating hair follicle development [48]. In the present study, miRNA expression in medium and small waves were significantly different, which was similar to that of let-7i expression in large and medium waves. These observations suggest that the differential expression of miRNA-10a and let-7i in different pattern lambskin hair follicles influences the growth and development of Hu sheep hair follicles.
In novel miRNAs, NW_004080184.1_6326 was significantly upregulated (P<0.01) in small and large waves, whereas it was downregulated in medium waves. The expression level of NW_004080165.1_8572 in large and small waves was higher than that in medium waves (P<0.01). Significant differences in the expression level of NW_004080181.1_3961 were observed between medium and small waves, as well as in NW_004080190.1_13733 expression between large and medium waves. Based on these findings, we inferred that the differential expression of these novel miRNAs in different lambskin hair follicles plays an important role in the formation of lambskin patterns. Additionally, the expression of the rest of the seven miRNAs showed no significant differences among the large, medium, and small waves, and thus were not considered as candidate miRNAs. We preliminarily designated miRNA-143, miRNA-10a, let-7i, NW_004080184.1_6326, NW_004080165.1_8572, NW_004080181.1_3961, and NW_004080190.1_13733 as candidate miRNAs that regulate the development of Hu sheep hair follicles.
miRNAs in different cells form a comprehensive and multi-level network system through the interactions between signaling pathways and regulatory factors [49]. A large number of signaling pathways and regulatory molecules are involved in the steady-state maintenance of skin and hair follicle development, which include the Wnt and Akt signaling pathways, BMP, FGF, and Hox [50–51]. Studies have confirmed that miRNAs play important roles in the development of tissues and organs such as the brain, heart, liver, kidney, and skin [52–56]. Furthermore, miRNAs show specific expression profiles at different developmental stages and in various tissues. Recent research investigations have demonstrated that the expression profiles of miRNAs are tissue-specific. Previous studies have indicated that miRNAs are involved in the growth and development of the skin and hair follicles of various animals. Further research has also shown that miRNAs are differentially expressed at various stages of hair follicle growth. Moreover, miRNAs are important factors in the control of the growth and development of skin and hair follicles [4]. The present study showed that at the same miRNA expression level, the number of secondary follicles in small waves was higher than that in large waves, which was in agreement with our histological results. Combined with the study on the differential expression of miRNA-143 in large, medium, and small wave patterns, we inferred that miRNA -143 regulates Hu sheep hair follicle development.
Although no studies on the involvement of miRNA-10a and let-7i in sheep hair follicle development have been conducted to date, our research investigation has shown that at the same miRNA-10a expression level, the diameter of primary follicles in large waves was larger than that in small waves, and the diameter of secondary follicles in large waves was smaller than that in medium and small waves, and these findings were in agreement with our biopsy results. Furthermore, at the same let-7i expression level, the number of secondary follicles in large waves was lower than that in small waves, whereas the biopsy showed that these were similar. These discrepancies may be attributable to our small sample size as it is very difficult to finding full sibs of Hu sheep with different patterns. However, at the same let-7i expression level, the diameter of primary follicles in large waves was larger than that in medium waves, and these results was in agreement with our biopsy findings. The results of our correlation analysis basically agree with the findings of our corresponding biopsy examinations. In addition, in novel miRNAs, a correlation between the expression levels of NW_004080184.1_6326, NW_004080165.1_8572, NW_004080181.1_3961, NW_004080190.1_13733, miRNA-143, miRNA-10a, and let-7i and the diameter of secondary follicles in large and medium waves and the number of secondary follicles in small waves was observed, which coincided with our biopsy results. No significant differences in the expression levels of other miRNAs in the large, medium, and small waves were observed, although these showed correlations with related hair follicle indicators and thus were excluded as candidate miRNAs.
Conclusions
In this study, 36 miRNAs were found to be differentially expressed, including 13 miRNAs among small and large waves, 20 miRNAs among small and medium waves, 14 miRNAs among medium and large waves. Gene ontology analysis indicated that they are involved in cell proliferation, differentiation, apoptosis, development, metabolization, and the immune response. The expression pattern of miR-143, miR-10a, miR-199p-3a, let-7c, let-7i, NW_004080189.1_7961, NW_004080184.1_6535, NW_004080184.1_6326, NW_004080165.1_8572, NW_004080166.1_9139, NW_004080181.1_3961, NW_004080174.1_726, NW_004080166.1_10417 and NW_004080190.1_13733 showed that they may participate in the formation of hair follicles and different flower pattern. The correlation analysis between these 14 miRNAs expression and histological properties of hair follicle in large, medium and small waves showed that miRNA-143, miRNA-10a, let-7i, NW_004080184.1_6326, NW_004080165.1_8572, NW_004080181.1_3961 and NW_004080190.1_13733 can be considered as the important candidate miRNA to Hu sheep lambskin hair follicle development. Our findings not only expanded the database on sheep miRNAs but also could be useful for future functional studies of miRNA involved in hair follicle development and different flower pattern formation in Hu sheep.
Supporting information
[Figure omitted. See PDF.]
S1 Table. Relevant information of genes and primer sequences for Real Time PCR.
https://doi.org/10.1371/journal.pone.0176532.s001
(XLSX)
S2 Table. Predicted targets of the significant differentially expressed miRNA in small waves and large waves.
https://doi.org/10.1371/journal.pone.0176532.s002
(XLSX)
S3 Table. Predicted targets of the significant differentially expressed miRNA in medium waves and large waves.
https://doi.org/10.1371/journal.pone.0176532.s003
(XLSX)
S4 Table. Predicted targets of the significant differentially expressed miRNA in small waves and medium waves.
https://doi.org/10.1371/journal.pone.0176532.s004
(XLSX)
Acknowledgments
This work was performed by the projects of domesticated animals platform of the Ministry of Science and Technology of China, the graduate education innovation project of Jiangsu province of China (KYLX16_1402), the Project Funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions, the Project of Jiangsu Province Engineering Research Center of China (BM2012308), the key project of the national Spark Program of China (2012GA690003), the Subei Science and Technology program of Jiangsu Province of China (BN2014003, BN2015013, BN2015014), the Project of Jiangsu Province Agricultural Science and Technology Innovation Fund (CX(14)2073)and the Project of six peak of talents of Jiangsu Province of China.
Author Contributions
1. Conceptualization: WS WG JFY.
2. Data curation: XYL WG.
3. Formal analysis: WG JFY.
4. Funding acquisition: WS.
5. Investigation: JJB LHW JRY.
6. Supervision: CYJ LH.
7. Writing – original draft: WG JFY.
8. Writing – review & editing: WS WG.
Citation: Gao W, Sun W, Yin J, Lv X, Bao J, Yu J, et al. (2017) Screening candidate microRNAs (miRNAs) in different lambskin hair follicles in Hu sheep. PLoS ONE 12(5): e0176532. https://doi.org/10.1371/journal.pone.0176532
1. Chen RC. The identification and classification of Hu sheep lambskin. Journal of Economic Animal. 1980; (1):40–41.
2. Li CQ, Yin J, Zhang YJ, Guo ZC, Zhang WG, Gao AQ et al. Comparative study on skin and hair follicles cycling between inner Mongolia and Liaoning cashmere goats. Acta Veterinaria et Zootechnica Sinica. 2005; 36(7): 674–679
3. He XY, Hao HQ, Liu DD, Fan RW, Cao J, Zhu ZW, et al. Difference of microRNA expression in the ear and back skin of young Alpaca(Lama acos). Chinese Journal of Biochemistry and Molecular Biology. 2010; 26(11): 1016–1022.
4. 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. Curr Biol. 2006; 16(10): 1041–1049. pmid:16682203
5. Zhu Z, He J, Jia X, Jiang J, Bai R, Yu X, et al. MicroRNA-25 functions in regulation of pigmentation by targeting the transcription factor MITF in alpaca (Lama pacos) skin melanocytes. Domest Anim Endocrinol. 2010; 38(3): 200–209. pmid:20036482
6. Mardaryev AN, Ahmed MI, Vlahov NV, Fessing MY, Gill JH, Sharov AA, et al. Micro-RNA-31 controls hair cycle-associated changes in gene expression programs of the skin and hair follicle. The FASEB Journal. 2010; 24(10): 3869–3881. pmid:20522784
7. Tang XH, Liu GB, Du XY, Cao JH, Li XY, Luo Z, et al. Expression variation of microRNA-31 in wool follicle of Tibetan sheep(Ovis aries). Journal of Agricultural Biotechnology. 2012; 20(10): 1178–1183.
8. Yi R, O'Carroll D, Pasolli HA, Zhang Z, Dietrich FS, Tarakhovsky A. Morphogenesis in skin is governed by discrete sets of differentially expressed microRNAs. Nature Genetics. 2006; 38(3): 356–362. pmid:16462742
9. Zhang L, Nie Q, Su Y, Xian X, Luo W, Jia X. et al. MicroRNA profile analysis on duck feather follicle and skin with high-throughput sequencing technology. Gene. 2013; 519(1): 77–81. pmid:23384715
10. Liu G, Liu R, Li Q, Tang X, Yu M, Li X, et al. Identification of microRNAs in wool follicles during anagen, catagen, and telogenphases in Tibetan sheep. PLoS One. 2013; 8(10): e77801. pmid:24204975
11. Zhang WG, Wu JH, Li JQ, Midori Y. A subset of skin-expressed microRNAs with possible roles in goat and sheep hair growth based on expression profiling of mammalian microRNAs. Omics A Journal of Integrative Biology. 2007; 11(4), 385–396. pmid:18092910
12. Liu Z, Xiao H, Li H, Zhao Y, Lai S, Yu X, et al. Identification of conserved and novel microRNAs in cashmere goat skin by deep sequencing. PLoS One. 2012; 7(12): e50001. pmid:23236360
13. Wu ZY, Fu YH, Cao JH, Yu M, Tang XH, Zhao SH. Identification of differentially expressed miRNAs between white and black hair follicles by RNA-Sequencing in the Goat(Capra hircus). International Journal of Molecular Sciences. 2014; 15(6), 9531–9545. pmid:24879525
14. Andl T, Botchkareva NV. MicroRNAs (miRNAs) in the control of HF development and cycling: the next frontiers in hair research. Exp Dermatol. 2015; 24(11): 821–826. pmid:26121602
15. Andl T, Reddy ST, Gaddapara T, Millar SE. WNT signals are required for the initiation of hair follicle development. Dev Cell. 2002; 2(5):643–653. pmid:12015971
16. Foitzik K, Paus R, Doetschman T, Dotto GP. The TGF-beta2 isoform is both a required and sufficient inducer of murine hair follicle morphogenesis. Dev Biol. 1999; 212(2): 278–289. pmid:10433821
17. Headon DJ, Overbeek PA. Involvement of a novel Tnf receptor homologue in hair follicle induction. Nat Genet. 1999; 22(4), 370–374. pmid:10431242
18. Boras K, Hamel PA. Alx4 binding to LEF-1 regulates N-CAM promoter activity. J Biol Chem. 2002; 277(2): 1120–1127. pmid:11696550
19. Botchkare VA, Sharov AA. BMP signaling in the control of skin development and hair follicle growth. Differentiation. 2004; 72(9–10):512–526. pmid:15617562
20. Assa-Kunik E, Torres IL, Schejter ED, Johnston DS, Shilo BZ. Drosophila follicle cells are patterned by multiple levels of Notch signaling and antagonism between the Notch and JAK/STAT pathways. Development. 2002; 134(6), 1161–1169.
21. Yang WB, Li QQ, Su B, Mei Y. MicroRNA-148b promotes proliferation of hair follicle cells by targeting. Frontiers of Agricultural Science and Engineering. 2016; 3, 72–80.
22. Tazawa H, Kagawa S, Fujiwara T. MicroRNAs as potential target gene in cancer gene therapy of gastrointestinal tumors. Expert Opin Biol Ther. 2011; 11(2): 145–55. pmid:21219233
23. Suzuki HI, Yamagata K, Sugimoto K, Lwamoto T, Kato S, Miyazono K. Modulation of microRNA processing by p53. Nature. 2009; 460: 529–533. pmid:19626115
24. Wu BL, Xu LY, Du ZP, Liao LD, Zhang HF, Huang Q, et al. MiRNA profile in esophageal squamous cell carcinoma: downregulation of miR-143 and miR-145. World J Gastroenterol. 2011; 17(1): 79–88. pmid:21218087
25. Takagi T, Iio A, Nakagawa Y, Naoe T, Tanigawa N, Akao Y. Decreased expression of microRNA-143 and -145 in human gastric cancers. Oncology. 2009; 77(1): 12–21. pmid:19439999
26. Navarro A, Diaz T, Gallardo E. Viñolas N, Marrades RM, Gel B, et al. Prognostic implications of miR-16 expression levels in resected non-small-cell lung cancer. J Surg Oncol. 2011; 103(5): 411–415. pmid:21400525
27. Lin T, Dong W., Huang J., Pan Q., Fan X., Zhang C, et al. MicroRNA-143 as a tumor suppressor for bladder cancer. J Urol. 2009; 181(3): 1372–1380. pmid:19157460
28. Wang X, Tang S, Le SY, Lu R, Rader JS, Meyers C, et all. Aberrant expression of oncogenic and tumor-suppressive microRNAs in cervical cancer is required for cancer cell growth. PLoS One. 2008; 3(7), e2557. pmid:18596939
29. Lu Z, Lam KS, Wang N, Xu X, Cortes M, Andersen B. LMO4 can interact with Smads Proteins and modulate transforming growth factor-beta signaling in epithelial cells. Oncogene. 2006; 25(20): 2920–2930. pmid:16331278
30. Sum EY, O'Reilly LA, Jonas N, Lindeman GJ, Visvader JE. The LIM domain protein Lmo4 is highly expressed in proliferating mouse epithelial tissues. J Histochem Cytochem. 2005; 53(4), 475–86. pmid:15805422
31. Holik AZ, Filby CE, Pasquet J, Viitaniemi K, Ciciulla J, Sutherland KD et al. The LIM-domain only protein 4 contributes to lung epithelial cell proliferation but is not essential for tumor progression. Respir Res. 2015; 16(), 67–81. pmid:26048572
32. Sugihara TM, Bach I, Kioussi C, Rosenfeld MG, Andersen B. Mouse deformed epidermal autoregulatory factor 1 recruits a LIM domain factor, LMO-4, and CLIM coregulators. Proc Natl Acad Sci USA. 1998; 95(26), 15418–15423. pmid:9860983
33. Wang N, Kudryavtseva E, Chen IL, McCormick J, Sugihara TM, Ruiz R, et al. Expression of an engrailed-LMO4 fusion protein in mammary epithelial cells inhibits mammary gland development in mice. Oncogene. 2004; 23(8), 1507–1513. pmid:14676840
34. Takamizawa J, Konishi H, Yanagisawa K, Tomida S, Osada H, Endoh H, et al. Reduced expression of the let-7 microRNAs in human lung cancers in association with shortened postoperative survival. Cancer Res. 2004; 64(11), 3753–3756. pmid:15172979
35. Yanaihara N, Caplen N, Bowman E, Seike M, Kumamoto K, Yi M, et al. Unique microRNA molecular profiles in lung cancer diagnosis and prognosis. Cancer Cell. 2006; 9(3), 189–198. pmid:16530703
36. Hou J, Lin L, Zhou W, Wang Z, Ding G, Dong Q, et al. Identification of miRNomes in Human Liver and Hepatocellular Carcinoma Reveals miR-199a/b-3p as Therapeutic Target for Hepatocellular Carcinoma. Cancer Cell. 2011; 19(2), 232–243. pmid:21316602
37. Cheng HH, Davis AJ, Lee TL, Pang AL, Nagrani S, Rennert OM, et al. Methylation of an intronic region regulates miR-199a in testicular tumor malignancy. Oncogene. 2011; 30(31), 3404–3415. pmid:21383689
38. Shen L, Sun CM, Li YY, Li XT, Sun T, Liu CJ, et al. MicroRNA-199a-3p suppresses glioma cell proliferation by regulating the AKT/mTOR signaling pathway. 2015; Tumour Biol. 36(9), 6929–6938. pmid:25854175
39. Yan Y, Luo YC, Wan HY, Wang J, Zhang PP, Liu M, et al. MicroRNA-10a is involved in the metastatic process by regulating Eph tyrosine kinase receptor A4-mediated epithelial -mesenchymal transition and adhesion in hepatoma cells. Hepatology. 2013; 57(2), 667–677. pmid:22996586
40. Agirre X, Jimenez-Velasco A, San Jose-Eneriz E, Bandres E, Corden L, Aparicio O, et al. Down-regulation of has-miR-10a in chronic myeloid leukemia CD34+ cells increases USE2-mediated cell growth. Mol Cancer Res. 2008; 6(12): 1830–1840. pmid:19074828
41. Grose R, Dickson C. Fibroblast growth factor signaling in tumorigenesis. Cytokine Growth Factor Rev. 2005; 16(2):179–186. pmid:15863033
42. Presta M, Dell’Era P, Mitola S, Moroni E, Ronca R, Rusnati M. Fibroblast growth factor/fibroblast growth factor receptor system in angiogenesis. Cytokine Growth Factor Rev. 2005; 16(2), 159–178. pmid:15863032
43. Pave A, Truonq A, Yip C. et al. EGFR activation and signaling in cancer cells are enhanced by the membrane-bound metalloprotease MT4-MMP. Cancer Res. 2014; 74(23), 6758–70. pmid:25320013
44. Yi R, O'Carroll D, Pasolli HA, Zhang Z, Dietrich FS, Tarakhovsky A. Morphogenesis in skin is governed by discrete sets of differentially expressed microRNAs. Nature Genetics. 2006; 38(3), 356–362. pmid:16462742
45. Yi R, Fuchs E. MicroRNA-mediated control in the skin. Cell Death Differ. 2010; 17(2), 229–235. pmid:19609273
46. Sehultz J, Lorenz P, Gross G, Ibrahim S, Kunz M. MicroRNA let-7b targets important celll cycle molecules in malignant melanomacells and interferes with anchomge independent growth. Cell Res. 2008; 18(5), 549–557. pmid:18379589
47. Corey S, Frucht , Santos-Sacchi J, Navaratnam DS. MicroRNA181a plays a key role in hair cell regeneration in avian auditory epithelium. Neurosci Lett. 2011; 493(1–2), 44–48. pmid:21316421
48. Jiang Q. Differential Expression of miR-let7a in hair follicle cycle of Liaoning Cashmere goats and Identification of its targets. Master’s Thesis, Jilin University, Ji Ling, 2015.
49. Cai T, Liu ZH, Wang ZX, Zhao M, Ju HL, Li JQ. miRNA in regulation of skin and hair follicle development. Yi Chuan. 2013; 35(9): 1087–1094. pmid:24400482
50. Sennett R, Rendl M. Mesenchymal-epithelial interactions during hair follicle morphogenesis and cycling. Semin Cell Dev Biol. 2012; 23(), 917–927. pmid:22960356
51. Wu JH, Yan ZW, Hu SL, Zhang GW,Li JQ. Hoxc13 and the development of hair follicle. Yi Chuan. 2010; 32(7): 656–662. pmid:20650845
52. Polajeva J, Swartling FJ, Jiang YW, Singh U, Pietras K, Uhrbom L, et al. miRNA-21 is developmentally regulated in mouse brain and is co-expressed with SOX2 in glioma. BMC Cancer. 2012; 12: 378. pmid:22931209
53. Stauffer BL, Russell G, Nunley K, Miyamoto SD, Sucharov CC. miRNA expression in pediatric failing human heart. J Mol Cell Cardiol. 2013; 57: 43–46. pmid:23333438
54. Zhang YQ, Wu LP, Wang Y, Zhang MC, Li LM, Zhu DH,et al. Protective role of estrogen-induced miRNA-29 expression in carbon tetrachloride-induced mouse liver injury. J Biol Chem. 2012; 287(18): 14851–14862. pmid:22393047
55. Timoneda O, Balcells I, Nunez JI, Egea R, Vera G, Castello A, et al. miRNA expression profile analysis in kidney of different porcine breeds. PLoS One. 2013; 8(1): e55402. pmid:23372853
56. Botchkareva NV. MicroRNA/mRNA regulatory networks in the control of skin development and regeneration. Cell Cycle. 2012; 11(3), 468–474. pmid:22262186
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2017 Gao et al. This is an open access article distributed under the terms of the Creative Commons Attribution License: http://creativecommons.org/licenses/by/4.0/ (the “License”), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
Hu sheep lambskin is a unique white lambskin from China that exhibits three types of flower patterns, including small waves, medium waves, and large waves, with small waves considered the best quality. However, our understanding of the molecular mechanism underlying flower pattern formation in Hu sheep lambskin is limited. The aim of the present study was to further explore the relevance between candidate microRNAs (miRNAs) and developmental characteristics of hair follicles and screen miRNAs for later functional validation. Herein, we employed Illumina Hiseq 2500 to identify differentially expressed miRNAs in hair follicles of different flower patterns with small, medium, and large waves to construct a comprehensive sequence database on the mechanism of hair follicle development. Paraffin sections of lambskin tissue were prepared to assess the structure of different hair follicles. Expression levels of candidate miRNAs in different flower patterns were analyzed by relative quantitation using real-time PCR, combined with histological observation and micro-observation technologies, and the correlation between expression levels of candidate miRNAs and histological properties of hair follicles was analyzed by using SPSS 17.0. A total of 522 differentially expressed miRNAs were identified, and RNA-seq analysis detected 7,266 target genes in different groups of flower patterns. Gene ontological analysis indicated these target genes were mainly involved in cell proliferation, differentiation, growth, apoptosis, and ion transport, and 14 miRNAs, including miR-143, miR-10a, and let-7 were screened as candidate miRNAs in Hu sheep hair follicle growth and development. In the same field of vision, variance analysis showed that the number of secondary follicles in small waves was significantly larger than that in large and medium waves (P<0.01); the diameter of the primary and secondary follicles in large waves was respectively larger than those in medium and small waves (P<0.01). Combined with correlation analysis between miRNA expression and histological properties of hair follicles, highly significant differences in miRNA-143 expression levels between large and small waves were observed (P<0.01), and significant differences in the miRNA-10a expression levels between large and small waves (P<0.05) and in let-7i expression levels between large and medium waves were observed (P<0.05). Significant differences in the expression of novel miRNAs of NW_004080184.1_6326 between medium and large waves were detected (P<0.05), and highly significant differences between medium and small waves were observed (P<0.01). Highly significant differences in the expression level of NW_004080165.1_8572 between medium and large and small waves (P<0.01), in that of NW_004080181.1_3961 between medium and small waves (P<0.01), and in that of NW_004080190.1_13733 between medium and large waves were observed, whereas no significant differences in the other miRNAs among large, medium, and small waves were detected. Overall, the present study showed that miRNA-143, miRNA-10a, let-7i, NW_004080184.1_6326, NW_004080165.1_8572, NW_004080181.1_3961, and NW_004080190.1_13733 could be considered as important candidate genes, indicating these seven miRNAs may play significant roles in hair follicle growth and development in Hu sheep lambskin.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer