Genetic diversity and selection of Tibetan sheep breeds revealed by whole-genome resequencing

Objective This study aimed to elucidate the underlying gene regions responsible for productive, phenotypic or adaptive traits in different ecological types of Tibetan sheep and the discovery of important genes encoding valuable traits. Methods We used whole-genome resequencing to explore the genetic relationships, phylogenetic tree, and population genetic structure analysis. In addition, we identified 28 representative Tibetan sheep single-nucleotide polymorphisms (SNPs) and genomic selective sweep regions with different traits in Tibetan sheep by fixation index (Fst) and the nucleotide diversity (θπ) ratio. Results The genetic relationships analysis showed that each breed partitioned into its own clades and had close genetic relationships. We also identified many potential breed-specific selective sweep regions, including genes associated with hypoxic adaptability (MTOR, TRHDE, PDK1, PTPN9, TMTC2, SOX9, EPAS1, PDGFD, SOCS3, TGFBR3), coat color (MITF, MC1R, ERCC2, TCF25, ITCH, TYR, RALY, KIT), wool traits (COL4A2, ERC2, NOTCH2, ROCK1, FGF5, SOX9), and horn phenotypes (RXFP2). In particular, a horn-related gene, RXFP2, showed the four most significantly associated SNP loci (g. 29481646 A>G, g. 29469024 T>C, g. 29462010 C>T, g. 29461968 C>T) and haplotypes. Conclusion This finding demonstrates the potential for genetic markers in future molecular breeding programs to improve selection for horn phenotypes. The results will facilitate the understanding of the genetic basis of production and adaptive unique traits in Chinese indigenous Tibetan sheep taxa and offer a reference for the molecular breeding of Tibetan sheep

Defining Key Genes Regulating Morphogenesis of Apocrine Sweat Gland in Sheepskin

The apocrine sweat gland is a unique skin appendage in humans compared to mouse and chicken models. The absence of apocrine sweat glands in chicken and murine skin largely restrains further understanding of the complexity of human skin biology and skin diseases, like hircismus. Sheep may serve as an additional system for skin appendage investigation owing to the distributions and histological similarities between the apocrine sweat glands of sheep trunk skin and human armpit skin. To understand the molecular mechanisms underlying morphogenesis of apocrine sweat glands in sheepskin, transcriptome analyses were conducted to reveal 1631 differentially expressed genes that were mainly enriched in three functional groups (cellular component, molecular function and biological process), particularly in gland, epithelial, hair follicle and skin development. There were 7 Gene Ontology (GO) terms enriched in epithelial cell migration and morphogenesis of branching epithelium that were potentially correlated with the wool follicle peg elongation. An additional 5 GO terms were enriched in gland morphogenesis (20 genes), gland development (42 genes), salivary gland morphogenesis and development (8 genes), branching involved in salivary gland morphogenesis (6 genes) and mammary gland epithelial cell differentiation (4 genes). The enriched gland-related genes and two Kyoto Encyclopedia of Genes and Genomes pathway genes (WNT and TGF-β) were potentially involved in the induction of apocrine sweat glands. Genes named BMPR1A, BMP7, SMAD4, TGFB3, WIF1, and WNT10B were selected to validate transcript expression by qRT-PCR. Immunohistochemistry was performed to localize markers for hair follicle (SOX2), skin fibroblast (PDGFRB), stem cells (SOX9) and BMP signaling (SMAD5) in sheepskin. SOX2 and PDGFRB were absent in apocrine sweat glands. SOX9 and SMAD5 were both observed in precursor cells of apocrine sweat glands and later in gland ducts. These results combined with the upregulation of BMP signaling genes indicate that apocrine sweat glands were originated from outer root sheath of primary wool follicle and positively regulated by BMP signaling. This report established the primary network regulating early development of apocrine sweat glands in sheepskin and will facilitate the further understanding of histology and pathology of apocrine sweat glands in human and companion animal skin

Transcriptome Reveals Long Non-coding RNAs and mRNAs Involved in Primary Wool Follicle Induction in Carpet Sheep Fetal Skin

Murine primary hair follicle induction is driven by the communication between the mesenchyme and epithelium and mostly governed by signaling pathways including wingless-related integration site (WNT), ectodysplasin A receptor (EDAR), bone morphogenetic protein (BMP), and fibroblast growth factor (FGF), as observed in genetically modified mouse models. Sheep skin may serve as a valuable system for hair research owing to the co-existence of sweat glands with wool follicles in trunk skin and asynchronized wool follicle growth pattern similar to that of human head hair follicles. However, the mechanisms underlying wool follicle development remain largely unknown. To understand how long non-coding RNAs (lncRNAs) and mRNAs function in primary wool follicle induction in carpet wool sheep, we conducted high-throughput RNA sequencing and revealed globally altered lncRNAs (36 upregulated and 26 downregulated), mRNAs (228 elevated and 225 decreased), and 80 differentially expressed novel transcripts. Several key signals in WNT (WNT2B and WNT16), BMP (BMP3, BMP4, and BMP7), EDAR (EDAR and EDARADD), and FGF (FGFR2 and FGF20) pathways, and a series of lncRNAs, including XLOC_539599, XLOC_556463, XLOC_015081, XLOC_1285606, XLOC_297809, and XLOC_764219, were shown to be potentially important for primary wool follicle induction. GO and KEGG analyses of differentially expressed mRNAs and potential targets of altered lncRNAs were both significantly enriched in morphogenesis biological processes and transforming growth factor-β, Hedgehog, and PI3K-Akt signaling, as well as focal adhesion and extracellular matrix-receptor interactions. The prediction of mRNA-mRNA and lncRNA-mRNA interaction networks further revealed transcripts potentially involved in primary wool follicle induction. The expression patterns of mRNAs and lncRNAs of interest were validated by qRT-PCR. The localization of XLOC_297809 and XLOC_764219 both in placodes and dermal condensations was detected by in situ hybridization, indicating important roles of lncRNAs in primary wool follicle induction and skin development. This is the first report elucidating the gene network of lncRNAs and mRNAs associated with primary wool follicle early development in carpet wool sheep and will shed new light on selective wool sheep breeding

Data_Sheet_1_Transcriptome Reveals Long Non-coding RNAs and mRNAs Involved in Primary Wool Follicle Induction in Carpet Sheep Fetal Skin.XLSX

<p>Murine primary hair follicle induction is driven by the communication between the mesenchyme and epithelium and mostly governed by signaling pathways including wingless-related integration site (WNT), ectodysplasin A receptor (EDAR), bone morphogenetic protein (BMP), and fibroblast growth factor (FGF), as observed in genetically modified mouse models. Sheep skin may serve as a valuable system for hair research owing to the co-existence of sweat glands with wool follicles in trunk skin and asynchronized wool follicle growth pattern similar to that of human head hair follicles. However, the mechanisms underlying wool follicle development remain largely unknown. To understand how long non-coding RNAs (lncRNAs) and mRNAs function in primary wool follicle induction in carpet wool sheep, we conducted high-throughput RNA sequencing and revealed globally altered lncRNAs (36 upregulated and 26 downregulated), mRNAs (228 elevated and 225 decreased), and 80 differentially expressed novel transcripts. Several key signals in WNT (WNT2B and WNT16), BMP (BMP3, BMP4, and BMP7), EDAR (EDAR and EDARADD), and FGF (FGFR2 and FGF20) pathways, and a series of lncRNAs, including XLOC_539599, XLOC_556463, XLOC_015081, XLOC_1285606, XLOC_297809, and XLOC_764219, were shown to be potentially important for primary wool follicle induction. GO and KEGG analyses of differentially expressed mRNAs and potential targets of altered lncRNAs were both significantly enriched in morphogenesis biological processes and transforming growth factor-β, Hedgehog, and PI3K-Akt signaling, as well as focal adhesion and extracellular matrix-receptor interactions. The prediction of mRNA-mRNA and lncRNA-mRNA interaction networks further revealed transcripts potentially involved in primary wool follicle induction. The expression patterns of mRNAs and lncRNAs of interest were validated by qRT-PCR. The localization of XLOC_297809 and XLOC_764219 both in placodes and dermal condensations was detected by in situ hybridization, indicating important roles of lncRNAs in primary wool follicle induction and skin development. This is the first report elucidating the gene network of lncRNAs and mRNAs associated with primary wool follicle early development in carpet wool sheep and will shed new light on selective wool sheep breeding.</p

Image_1_Transcriptome Reveals Long Non-coding RNAs and mRNAs Involved in Primary Wool Follicle Induction in Carpet Sheep Fetal Skin.TIF

<p>Murine primary hair follicle induction is driven by the communication between the mesenchyme and epithelium and mostly governed by signaling pathways including wingless-related integration site (WNT), ectodysplasin A receptor (EDAR), bone morphogenetic protein (BMP), and fibroblast growth factor (FGF), as observed in genetically modified mouse models. Sheep skin may serve as a valuable system for hair research owing to the co-existence of sweat glands with wool follicles in trunk skin and asynchronized wool follicle growth pattern similar to that of human head hair follicles. However, the mechanisms underlying wool follicle development remain largely unknown. To understand how long non-coding RNAs (lncRNAs) and mRNAs function in primary wool follicle induction in carpet wool sheep, we conducted high-throughput RNA sequencing and revealed globally altered lncRNAs (36 upregulated and 26 downregulated), mRNAs (228 elevated and 225 decreased), and 80 differentially expressed novel transcripts. Several key signals in WNT (WNT2B and WNT16), BMP (BMP3, BMP4, and BMP7), EDAR (EDAR and EDARADD), and FGF (FGFR2 and FGF20) pathways, and a series of lncRNAs, including XLOC_539599, XLOC_556463, XLOC_015081, XLOC_1285606, XLOC_297809, and XLOC_764219, were shown to be potentially important for primary wool follicle induction. GO and KEGG analyses of differentially expressed mRNAs and potential targets of altered lncRNAs were both significantly enriched in morphogenesis biological processes and transforming growth factor-β, Hedgehog, and PI3K-Akt signaling, as well as focal adhesion and extracellular matrix-receptor interactions. The prediction of mRNA-mRNA and lncRNA-mRNA interaction networks further revealed transcripts potentially involved in primary wool follicle induction. The expression patterns of mRNAs and lncRNAs of interest were validated by qRT-PCR. The localization of XLOC_297809 and XLOC_764219 both in placodes and dermal condensations was detected by in situ hybridization, indicating important roles of lncRNAs in primary wool follicle induction and skin development. This is the first report elucidating the gene network of lncRNAs and mRNAs associated with primary wool follicle early development in carpet wool sheep and will shed new light on selective wool sheep breeding.</p

Data_Sheet_6_Transcriptome Reveals Long Non-coding RNAs and mRNAs Involved in Primary Wool Follicle Induction in Carpet Sheep Fetal Skin.XLSX

<p>Murine primary hair follicle induction is driven by the communication between the mesenchyme and epithelium and mostly governed by signaling pathways including wingless-related integration site (WNT), ectodysplasin A receptor (EDAR), bone morphogenetic protein (BMP), and fibroblast growth factor (FGF), as observed in genetically modified mouse models. Sheep skin may serve as a valuable system for hair research owing to the co-existence of sweat glands with wool follicles in trunk skin and asynchronized wool follicle growth pattern similar to that of human head hair follicles. However, the mechanisms underlying wool follicle development remain largely unknown. To understand how long non-coding RNAs (lncRNAs) and mRNAs function in primary wool follicle induction in carpet wool sheep, we conducted high-throughput RNA sequencing and revealed globally altered lncRNAs (36 upregulated and 26 downregulated), mRNAs (228 elevated and 225 decreased), and 80 differentially expressed novel transcripts. Several key signals in WNT (WNT2B and WNT16), BMP (BMP3, BMP4, and BMP7), EDAR (EDAR and EDARADD), and FGF (FGFR2 and FGF20) pathways, and a series of lncRNAs, including XLOC_539599, XLOC_556463, XLOC_015081, XLOC_1285606, XLOC_297809, and XLOC_764219, were shown to be potentially important for primary wool follicle induction. GO and KEGG analyses of differentially expressed mRNAs and potential targets of altered lncRNAs were both significantly enriched in morphogenesis biological processes and transforming growth factor-β, Hedgehog, and PI3K-Akt signaling, as well as focal adhesion and extracellular matrix-receptor interactions. The prediction of mRNA-mRNA and lncRNA-mRNA interaction networks further revealed transcripts potentially involved in primary wool follicle induction. The expression patterns of mRNAs and lncRNAs of interest were validated by qRT-PCR. The localization of XLOC_297809 and XLOC_764219 both in placodes and dermal condensations was detected by in situ hybridization, indicating important roles of lncRNAs in primary wool follicle induction and skin development. This is the first report elucidating the gene network of lncRNAs and mRNAs associated with primary wool follicle early development in carpet wool sheep and will shed new light on selective wool sheep breeding.</p

Data_Sheet_3_Transcriptome Reveals Long Non-coding RNAs and mRNAs Involved in Primary Wool Follicle Induction in Carpet Sheep Fetal Skin.XLSX

<p>Murine primary hair follicle induction is driven by the communication between the mesenchyme and epithelium and mostly governed by signaling pathways including wingless-related integration site (WNT), ectodysplasin A receptor (EDAR), bone morphogenetic protein (BMP), and fibroblast growth factor (FGF), as observed in genetically modified mouse models. Sheep skin may serve as a valuable system for hair research owing to the co-existence of sweat glands with wool follicles in trunk skin and asynchronized wool follicle growth pattern similar to that of human head hair follicles. However, the mechanisms underlying wool follicle development remain largely unknown. To understand how long non-coding RNAs (lncRNAs) and mRNAs function in primary wool follicle induction in carpet wool sheep, we conducted high-throughput RNA sequencing and revealed globally altered lncRNAs (36 upregulated and 26 downregulated), mRNAs (228 elevated and 225 decreased), and 80 differentially expressed novel transcripts. Several key signals in WNT (WNT2B and WNT16), BMP (BMP3, BMP4, and BMP7), EDAR (EDAR and EDARADD), and FGF (FGFR2 and FGF20) pathways, and a series of lncRNAs, including XLOC_539599, XLOC_556463, XLOC_015081, XLOC_1285606, XLOC_297809, and XLOC_764219, were shown to be potentially important for primary wool follicle induction. GO and KEGG analyses of differentially expressed mRNAs and potential targets of altered lncRNAs were both significantly enriched in morphogenesis biological processes and transforming growth factor-β, Hedgehog, and PI3K-Akt signaling, as well as focal adhesion and extracellular matrix-receptor interactions. The prediction of mRNA-mRNA and lncRNA-mRNA interaction networks further revealed transcripts potentially involved in primary wool follicle induction. The expression patterns of mRNAs and lncRNAs of interest were validated by qRT-PCR. The localization of XLOC_297809 and XLOC_764219 both in placodes and dermal condensations was detected by in situ hybridization, indicating important roles of lncRNAs in primary wool follicle induction and skin development. This is the first report elucidating the gene network of lncRNAs and mRNAs associated with primary wool follicle early development in carpet wool sheep and will shed new light on selective wool sheep breeding.</p

Data_Sheet_5_Transcriptome Reveals Long Non-coding RNAs and mRNAs Involved in Primary Wool Follicle Induction in Carpet Sheep Fetal Skin.XLSX

<p>Murine primary hair follicle induction is driven by the communication between the mesenchyme and epithelium and mostly governed by signaling pathways including wingless-related integration site (WNT), ectodysplasin A receptor (EDAR), bone morphogenetic protein (BMP), and fibroblast growth factor (FGF), as observed in genetically modified mouse models. Sheep skin may serve as a valuable system for hair research owing to the co-existence of sweat glands with wool follicles in trunk skin and asynchronized wool follicle growth pattern similar to that of human head hair follicles. However, the mechanisms underlying wool follicle development remain largely unknown. To understand how long non-coding RNAs (lncRNAs) and mRNAs function in primary wool follicle induction in carpet wool sheep, we conducted high-throughput RNA sequencing and revealed globally altered lncRNAs (36 upregulated and 26 downregulated), mRNAs (228 elevated and 225 decreased), and 80 differentially expressed novel transcripts. Several key signals in WNT (WNT2B and WNT16), BMP (BMP3, BMP4, and BMP7), EDAR (EDAR and EDARADD), and FGF (FGFR2 and FGF20) pathways, and a series of lncRNAs, including XLOC_539599, XLOC_556463, XLOC_015081, XLOC_1285606, XLOC_297809, and XLOC_764219, were shown to be potentially important for primary wool follicle induction. GO and KEGG analyses of differentially expressed mRNAs and potential targets of altered lncRNAs were both significantly enriched in morphogenesis biological processes and transforming growth factor-β, Hedgehog, and PI3K-Akt signaling, as well as focal adhesion and extracellular matrix-receptor interactions. The prediction of mRNA-mRNA and lncRNA-mRNA interaction networks further revealed transcripts potentially involved in primary wool follicle induction. The expression patterns of mRNAs and lncRNAs of interest were validated by qRT-PCR. The localization of XLOC_297809 and XLOC_764219 both in placodes and dermal condensations was detected by in situ hybridization, indicating important roles of lncRNAs in primary wool follicle induction and skin development. This is the first report elucidating the gene network of lncRNAs and mRNAs associated with primary wool follicle early development in carpet wool sheep and will shed new light on selective wool sheep breeding.</p

Data_Sheet_2_Transcriptome Reveals Long Non-coding RNAs and mRNAs Involved in Primary Wool Follicle Induction in Carpet Sheep Fetal Skin.XLSX

<p>Murine primary hair follicle induction is driven by the communication between the mesenchyme and epithelium and mostly governed by signaling pathways including wingless-related integration site (WNT), ectodysplasin A receptor (EDAR), bone morphogenetic protein (BMP), and fibroblast growth factor (FGF), as observed in genetically modified mouse models. Sheep skin may serve as a valuable system for hair research owing to the co-existence of sweat glands with wool follicles in trunk skin and asynchronized wool follicle growth pattern similar to that of human head hair follicles. However, the mechanisms underlying wool follicle development remain largely unknown. To understand how long non-coding RNAs (lncRNAs) and mRNAs function in primary wool follicle induction in carpet wool sheep, we conducted high-throughput RNA sequencing and revealed globally altered lncRNAs (36 upregulated and 26 downregulated), mRNAs (228 elevated and 225 decreased), and 80 differentially expressed novel transcripts. Several key signals in WNT (WNT2B and WNT16), BMP (BMP3, BMP4, and BMP7), EDAR (EDAR and EDARADD), and FGF (FGFR2 and FGF20) pathways, and a series of lncRNAs, including XLOC_539599, XLOC_556463, XLOC_015081, XLOC_1285606, XLOC_297809, and XLOC_764219, were shown to be potentially important for primary wool follicle induction. GO and KEGG analyses of differentially expressed mRNAs and potential targets of altered lncRNAs were both significantly enriched in morphogenesis biological processes and transforming growth factor-β, Hedgehog, and PI3K-Akt signaling, as well as focal adhesion and extracellular matrix-receptor interactions. The prediction of mRNA-mRNA and lncRNA-mRNA interaction networks further revealed transcripts potentially involved in primary wool follicle induction. The expression patterns of mRNAs and lncRNAs of interest were validated by qRT-PCR. The localization of XLOC_297809 and XLOC_764219 both in placodes and dermal condensations was detected by in situ hybridization, indicating important roles of lncRNAs in primary wool follicle induction and skin development. This is the first report elucidating the gene network of lncRNAs and mRNAs associated with primary wool follicle early development in carpet wool sheep and will shed new light on selective wool sheep breeding.</p

Data_Sheet_4_Transcriptome Reveals Long Non-coding RNAs and mRNAs Involved in Primary Wool Follicle Induction in Carpet Sheep Fetal Skin.XLSX

<p>Murine primary hair follicle induction is driven by the communication between the mesenchyme and epithelium and mostly governed by signaling pathways including wingless-related integration site (WNT), ectodysplasin A receptor (EDAR), bone morphogenetic protein (BMP), and fibroblast growth factor (FGF), as observed in genetically modified mouse models. Sheep skin may serve as a valuable system for hair research owing to the co-existence of sweat glands with wool follicles in trunk skin and asynchronized wool follicle growth pattern similar to that of human head hair follicles. However, the mechanisms underlying wool follicle development remain largely unknown. To understand how long non-coding RNAs (lncRNAs) and mRNAs function in primary wool follicle induction in carpet wool sheep, we conducted high-throughput RNA sequencing and revealed globally altered lncRNAs (36 upregulated and 26 downregulated), mRNAs (228 elevated and 225 decreased), and 80 differentially expressed novel transcripts. Several key signals in WNT (WNT2B and WNT16), BMP (BMP3, BMP4, and BMP7), EDAR (EDAR and EDARADD), and FGF (FGFR2 and FGF20) pathways, and a series of lncRNAs, including XLOC_539599, XLOC_556463, XLOC_015081, XLOC_1285606, XLOC_297809, and XLOC_764219, were shown to be potentially important for primary wool follicle induction. GO and KEGG analyses of differentially expressed mRNAs and potential targets of altered lncRNAs were both significantly enriched in morphogenesis biological processes and transforming growth factor-β, Hedgehog, and PI3K-Akt signaling, as well as focal adhesion and extracellular matrix-receptor interactions. The prediction of mRNA-mRNA and lncRNA-mRNA interaction networks further revealed transcripts potentially involved in primary wool follicle induction. The expression patterns of mRNAs and lncRNAs of interest were validated by qRT-PCR. The localization of XLOC_297809 and XLOC_764219 both in placodes and dermal condensations was detected by in situ hybridization, indicating important roles of lncRNAs in primary wool follicle induction and skin development. This is the first report elucidating the gene network of lncRNAs and mRNAs associated with primary wool follicle early development in carpet wool sheep and will shed new light on selective wool sheep breeding.</p