Transcriptome profile at different physiological stages reveals potential mode for curly fleece in Chinese tan sheep

Xiaolong Kang, Gang Liu, Yufang Liu, Qinqin Xu, Ming Zhang, Meiying Fang, Xiaolong Kang, Gang Liu, Yufang Liu, Qinqin Xu, Ming Zhang, Meiying Fang

Abstract

Tan sheep (Ovis aries), a Chinese indigenous breed, has special curly fleece after birth, especially at one month old. However, this unique phenotype disappears gradually with age and the underlying reasons of trait evolvement are still unknown. In this study, skin transcriptome data was used to study this issue. In total 51,215 transcripts including described transcripts and transfrags were identified. Pathway analysis of the top 100 most highly expressed transcripts, which included TCHH and keratin gene family members, such as KRT25, KRT5, KRT71, KRT14 and others, showed pathways known to be relevant to hair/fleece development and function. Six hundred differentially expressed (DE) transcripts were detected at two different physiological ages (one-month-old with curly fleece and 48-month-old without curly fleece) and were categorized into three major functional groups: cellular component, molecular function, and biological process. The top six functional categories included cell, cell part, cellular process, binding, intracellular, metabolic process. The detected differentially expressed genes were particularly involved in signal, signal peptide, disulfide bond, glycoprotein and secreted terms, respectively. Further splicing isoform analysis showed that the metallothionein 3 isoform was up-regulated in Tan lamb skin, indicating that it may be related to the conformation of curly fleece in Chinese Tan lamb. The hair-related important differentially expressed genes (SPINK4, FGF21, ESRα, EphA3, NTNG1 and GPR110) were confirmed by qPCR analysis. We deduced that the differences existed in expressed transcripts, splice isoforms and GO categories between the two different physiological stages, which might constitute the major reasons for explaining the trait evolvement of curly fleece in Chinese Tan sheep. This study provides some clues for elucidating the molecular mechanism of fleece change with age in Chinese Tan sheep, as well as supplying some potential values for understanding human hair disorder and texture changes.

Conflict of interest statement

Competing Interests: The authors have declared that no competing interests exist.

Figures

Figure 1. The Chinese Tan sheep and…
Figure 1. The Chinese Tan sheep and its fleece.
A: The Chinese Tan lamb (L); B: The pelt of Tan lamb at one month old; C: The Chinese adult Tan sheep (A); D: The pelt of Tan adult sheep.
Figure 2. Expression of the top 20…
Figure 2. Expression of the top 20 most highly expressed genes in sheep skin.
The x-axis shows gene ID; y-axis shows gene expression level (FPKM).
Figure 3. Statistics of mainly alternative splicing…
Figure 3. Statistics of mainly alternative splicing events.
The first column shows the types of alternative transcript events; the second column shows the splicing graphs and the third to fifth column shows the number of AS events in A and L combined group, A group and L group, respectively.
Figure 4. Functional categorization of differentially expressed…
Figure 4. Functional categorization of differentially expressed genes based on known genes in the Uniprot database.
The x-axis shows the 2nd level term of Gene Ontology; y-axis shows percent of genes in DE genes.
Figure 5. The expression level of differently…
Figure 5. The expression level of differently expressed genes validated by qPCR.

References

    1. Cui ZJ, Zhang YL, Jiang Y, Chen GN, Xu Z (1962) Breeding report of Tan sheep. Part 1: The relationship of ecology and reproduction of Chinese Tan sheep. China Anim Husb Veter Med 4: 1–5.
    1. Lu TA, Fan T, Zhang SD, Zhang HW, Ma ZY (1986) Fleece color inheritance on Tan sheep. J Gansu Agri University 3: 14–19.
    1. Wei ZQ, Yu HC, Zhang ZH, Chen RX, Gong WH (2003) Comparison of hemoglobin polymorphism of big size strain, fur strain and their hybrids in Tan sheep. China Herbivore Sci 1: 57–58.
    1. Zhang R, Li AH, Yang LJ, Chen ZQ, Ren DX, et al. (2010) Correlation between KRT1.2 gene and properties of lamb fur qualities of Tan sheep in Ningxia. J Agri Sci 31: 27–30.
    1. Yang LJ, Li AH, Zhang R, Chen ZQ, Ren DX, et al. (2010) Correlation between KAP1.3 gene and fur quality characteristics in Ningxia Tan sheep. J Ningxia University (Nat Sci Edition) 31: 381–388.
    1. Cox LA, Glenn JP, Spradling KD, Nijland MJ, Garcia R, et al. (2012) A genome resource to address mechanisms of developmental programming: determination of the fetal sheep heart transcriptome. Physiology 590: 2873–2884.
    1. Ren H, Li L, Su H, Xu L, Wei C, et al. (2011) Histological and transcriptome-wide level characteristics of fetal myofiber hyperplasia during the second half of gestation in Texel and Ujumqin sheep. BMC Genomics 12: 411.
    1. Jäger M, Ott CE, Grünhagen J, Hecht J, Schell H, et al. (2011) Composite transcriptome assembly of RNA-seq data in a sheep model for delayed bone healing. BMC Genomics 12: 158.
    1. Geng R, Yuan C, Chen Y (2013) Exploring Differentially Expressed Genes by RNA-Seq in Cashmere Goat (Capra hircus) Skin during Hair Follicle Development and Cycling. PloS One 8: e62704.
    1. Xu T, Guo X, Wang H, Du X, Gao X, et al. (2013) De novo transcriptome assembly and differential gene expression profiling of three capra hircus skin types during anagen of the hair growth cycle. Inter J Genomics 2013: 269191.
    1. McDermott BM, Baucom JM, Hudspeth AJ (2007) Analysis and functional evaluation of the hair-cell transcriptome. PNAS 104: 11820–11825.
    1. Kurek D, Garinis GA, van Doorninck JH, van der Wees J, Grosveld FG (2007) Transcriptome and phenotypic analysis reveals Gata3-dependent signalling pathways in murine hair follicles. Development 134: 261–272.
    1. Okano J, Levy C, Lichti U, Sun HW, Yuspa SH, et al. (2012) Cutaneous retinoic acid levels determine hair follicle development and downgrowth. J Biological Chemistry 287: 39304–39315.
    1. Langmead B, Trapnell C, Pop M, Salzberg SL (2009) Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol 10: R25.
    1. Trapnell C, Pachter L, Salzberg SL (2009) TopHat: discovering splice junctions with RNA-Seq. Bioinformatics 25: 1105–1111.
    1. Li RQ, Yu C, Li YR, Lam TW, Yiu SM, et al. (2009) SOAP2: an improved ultrafast tool for short read alignment. Bioinformatics 25: 1966–1967.
    1. Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, et al. (2010) Transcript assembly and quantification by RNA-seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotech 28: 511–515.
    1. Roberts A, Pimentel H, Trapnell C, Pachter L (2011) Identification of novel transcripts in annotated genomes using RNA-seq. Bioinformatics 27: 2325–2329.
    1. Sammeth M, Foissac S, Guigó R (2008) A general definition and nomenclature for alternative splicing events. PLoS Comput Biol 4: e1000147.
    1. Anders S, Huber W (2010) Differential expression analysis for sequence count data. Genome Biol 11: R106.
    1. Benjamini Y, Hochberg Y (1995) Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Royal Statistical Soc Series B (Methodol) 57: 289–300.
    1. Young MD, Wakefield MJ, Smyth GK, Oshlack A (2010) Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol 11: R14.
    1. Mao XZ, Cai T, Olyarchuk JG, Wei LP (2005) Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics 21: 3787–3793.
    1. Rabbani P, Takeo M, Chou W, Myung P, Bosenberg M, et al. (2011) Coordinated activation of Wnt in epithelial and melanocyte stem cells initiates pigmented hair regeneration. Cell 145: 941–955.
    1. Oh HS, Smart RC (1996) An estrogen receptor pathway regulates the telogen-anagen hair follicle transition and influences epidermal cell proliferation. Proc Natl Acad Sci USA 93: 12525–12530.
    1. Shimomura Y, Wajid M, Ishii Y, Shapiro L, Petukhova L, et al. (2008) Disruption of P2RY5, an orphan G protein–coupled receptor, underlies autosomal recessive woolly hair. Nat Genet 40: 335–339.
    1. Midorikawa T, Chikazawa T, Yoshino T, Takada K, Arase S (2004) Different gene expression profile observed in dermal papilla cells related to androgenic alopecia by DNA macroarray analysis. J Dermatol Sci 36: 25–32.
    1. Yamada Y, Midorikawa T, Oura H, Yoshino T, Ohdera M, et al. (2008) Ephrin-A3 not only increases the density of hair follicles but also accelerates anagen development in neonatal mice. J Dermatol Sci 52: 178–185.
    1. Ahmad N, Mukhtar H (2004) Cytochrome P450: a target for drug development for skin diseases. J Invest Dermatol 123: 417–425.
    1. Schröder N, Sekhar A, Geffers I, Müller J, Dittrich–Breiholz O, et al. (2006) Identification of mouse genes with highly specific expression patterns in differentiated intestinal epithelium. Gastroenterology 130: 902–907.
    1. Fang Y, Liu WJ, Zhang FQ, Shao YG, Yu SG (2010) The Polymorphism of a Novel Mutation of KAP13. 1 Gene and its Associations with Cashmere Traits on Xinjiang Local Goat Breed in China. Asian J Anim Vet 5: 34–42.
    1. Saito H, Yasumoto KI, Takeda K, Takahashi K, Yamanmoto H, et al. (2003) Microphthalmia associated transcription factor in the Wnt signaling pathway. Pigment Cell Res. 16: 261–265.
    1. Guenzi E, Töpolt K, Cornali E, Lubeseder-Martellato C, Jörg A, et al. (2001) The helical domain of GBP-1 mediates the inhibition of endothelial cell proliferation by inflammatory cytokines. J EMBO 20: 5568–5577.
    1. Thibaut S, Bernard BA (2005) The biology of hair shape. Inter J Dermatol 44: 2–3.
    1. Bawden CS, McLaughlan C, Nesci A, Rogers G (2001) A unique type I keratin intermediate filament gene family is abundantly expressed in the inner root sheaths of sheep and human hair follicles. J Invest Dermatol 116: 157–166.
    1. Medland SE, Nyholt DR, Painter JN, McEvoy BP, McRae AF, et al. (2009) Common variants in the trichohyalin gene are associated with straight hair in Europeans. The American J Human Genetics 85: 750–755.
    1. Fietz MJ, McLaughlan CJ, Campbell MT, Rogers G (1993) Analysis of the sheep trichohyalin gene: potential structural and calcium-binding roles of trichohyalin in the hair follicle. J Cell Biol 121: 855–865.
    1. Kikkawa Y, Oyama A, Ishii R, Miura I, Amano T, et al. (2003) A small deletion hotspot in the type II keratin gene mK6irs1/Krt2-6g on mouse chromosome 15, a candidate for causing the wavy hair of the caracul (Ca) mutation. Genetics 165: 721–733.
    1. Cadieu E, Neff MW, Quignon P, Walsh K, Chase K, et al. (2009) Coat variation in the domestic dog is governed by variants in three genes. Science 326: 150–153.
    1. Nakatake Y, Hoshikawa M, Asaki T, Kassai Y, Itoh N (2001) Identification of a novel fibroblast growth factor, FGF-22, preferentially expressed in the inner root sheath of the hair follicle. Biochimica et Biophysica Acta 1517: 460–463.
    1. Kawano M, Komi-Kuramochi A, Asada M, Suzuki M, Oki J, et al. (2005) Comprehensive analysis of FGF and FGFR expression in skin: FGF18 is highly expressed in hair follicles and capable of inducing anagen from telogen stage hair follicles. J Invest Dermatol 124: 877–885.
    1. Wu XH, Fu Y, Yang DY, Zhang RH, Zheng WP, et al. (2012) Detailed transcriptome description of the neglected cestode taenia multiceps. PLoS One 7: e45830.
    1. Mikkola ML (2008) TNF superfamily in skin appendage development. Cytokine Growth Factor Rev 19: 219–230.
    1. Sharov AA, Schroeder M, Sharova TY, et al. (2011) Matrix metalloproteinase-9 is involved in the regulation of hair canal formation. J Invest Dermatol 131: 257–260.
    1. Bejsovec A (2000) Wnt signaling: an embarrassment of receptors. Current Biol 10: R919–R922.
    1. Wodarz A, Nusse R (1998) Mechanisms of Wnt signaling in development. Annu Rev Cell Dev Biol 14: 59–88.
    1. Glinka A, Wu W, Delius H, Monaghan AP, Blumenstock C, et al. (1998) Dickkopf-1 is a member of a new family of secreted proteins and functions in head induction. Nature 391: 357–362.
    1. Packman S, Palmiter RD, Karin M, O'Toole C (1987) Metallothionein messenger RNA regulation in the mottled mouse and Menkes kinky hair syndrome. J Clinical Invest 79: 1338–1342.

Source: PubMed

3
Iratkozz fel