Knockdown of miR-361-5p promotes the induced activation of SHF-stem cells through FOXM1 mediated Wnt/ß-catenin pathway in cashmere goats.

The inducing activation event of secondary hair follicle (SHF)-stem cells is considered a key biological process in the SHF regeneration, and the morphogenesis of cashmere fiber in cashmere goats. The miR-361-5p was essentially implicated in the induced activation of SHF-stem cells of cashmere goats, but its functional mechanisms are unclear. Here, we confirmed miR-361-5p was significantly downregulated in anagen SHF bugle of cashmere goats compared with that at telogen, and miR-361-5p expression was significantly lower in SHF-stem cells after activation than its counterpart before activation. Further, we found that miR-361-5p could negatively regulate the induced activation event of SHF-stem cells in cashmere goats. Mechanistically, through dual-luciferase reporter assays, miR-361-5p specifically bound with FOXM1 mRNA in SHF-stem cells of cashmere goats and negatively regulated the expression of FOXM1 gene. Also, through overexpression/knockdown analysis of FOXM1 gene, our results indicated that FOXM1 upregulated the expression of Wnt/ß-catenin pathway related genes in SHF-stem cells. Moreover, based on TOP/FOP-flash Wnt report assays, the knockdown of miR-361-5p promotes the Wnt/ß-catenin pathway activation through upregulating the FOXM1 expression in SHF-stem cells. Finally, we demonstrated that miR-361-5p negatively regulated the induced activation of SHF-stem cells through FOXM1 mediated Wnt/ß-catenin pathway in cashmere goats.

CircRNA-1967 participates in the differentiation of goat SHF-SCs into hair follicle lineage by sponging miR-93-3p to enhance LEF1 expression.

Circular RNAs (circRNAs), a novel class of non-coding RNAs, can interact with miRNAs through a sequence-driven sponge mechanism, thereby regulating the expression of their downstream target genes. CircRNA-1967 was found in secondary hair follicles (SHFs) of cashmere goats, but its functions are not clear. Here, we showed that both circRNA-1967 and its host gene BNC2 had significantly higher expression in SHF bulge at anagen than those at telogen of cashmere goats. Also, circRNA-1967 participates in the differentiation of SHF stem cells (SHF-SCs) into hair follicle lineage in cashmere goats. RNA pull-down assay verified that circRNA-1967 interacts with miR-93-3p. We also indicated that circRNA-1967 promoted LEF1 expression in SHF-SCs of cashmere goats. By dual-luciferase reporter analysis, we found that circRNA-1967 up-regulated LEF1 expression through the miR-93-3p-mediated pathway. The results from this study demonstrated that circRNA-1967 participated in the differentiation of goat SHF-SCs into hair follicle lineage by sponging miR-93-3p to enhance LEF1 expression. Our founding might constitute a novel pathway for revealing the potential mechanism of the differentiation of SHF-SCs into hair follicle lineage in cashmere goats. Also, these results provided a valuable basis for further enhancing the intrinsic regeneration of cashmere goat SHFs with the formation and growth of cashmere fibers.

Analysis of m6A methylation in skin tissues of different sex Liaoning cashmere goats.

N6-methyladenosine (m6A) is the most frequent internal modification of mRNA and lncRNA in eukaryotes. We used two high-throughput sequencing method, m6A-seq and RNA-seq to identify pivotal m6A-modified genes in cashmere fineness and fiber growth. 8062 m6A peaks were detected by m6A-seq, including 2157 upregulated and 6445 downregulated. Furthermore, by comparing m6A-modified genes of the male Liaoning Cashmere Goat (M-LCG) and female Liaoning Cashmere Goat (F-LCG) skin tissues, we get 862 differentially expressed m6A-modified genes. To identify differently expressed m6A genes associated with cashmere fineness, 11 genes were selected for validation using real time fluorescent quantitative PCR in M-LCG and F-LCG. This study provides an acadamic basis on the molecular regulation mechanism of m6A modification in cashmere growth process.

[Progress in goat genome studies].

The goat genome is the research basis for the protection and utilization of goat resources, which is important for breeding and improving goat breeds. At present, with the continuous improvement of goat reference genome, various important research progress in goat origin, evolution and adaptability has been achieved. In this review, we summarize the research progress in the goat genome in detail, encompassing goat genome structure, genome map (genetic, physical and comparative maps), goat high throughput sequencing and SNP chip development. We aim to provide a theoretical foundation for the development of goat genome selection.

Intensity modulation of a terahertz bandpass filter: utilizing image currents induced on MEMS reconfigurable metamaterials.

We experimentally demonstrated a tunable terahertz bandpass filter based on microelectromechanical systems (MEMS) reconfigurable metamaterials. The unit cell of the filter consists of two split-ring resonators (SRRs) and a movable bar. Initially, the movable bar situates at the center of the unit cell, and the filter has two passbands whose central frequencies locate at 0.65 and 0.96 THz. The intensity of the two passbands can be actively modulated by the movable bar, and a maximum modulation depth of 96% is achieved at 0.96 THz. The mechanism of tunability is investigated using the finite-integration time-domain method. The result shows that the image currents induced on the movable bar are opposite the resonance currents induced on the SRRs and, thus, weaken the oscillating intensity of the resonance currents. This scheme paves the way to dynamically control and switch the terahertz wave at some constant frequencies utilizing induced image currents.

[The transcriptome research progresses of skin hair follicle development].

Recently, transcriptome sequencing technology has achieved significant progresses in gene network regulation of important economic traits in animals. As the derivative of mammalian skin, hair follicle is capable of self-renew. Its proliferation and differentiation result in hair formation. Researches have revealed that many growth factors and receptors coordinate genes and environment, as well as play an extremely important role during hair growth. In this review, we summarize the progresses that transcriptome sequencing technologies have achieved in researches of hair follicle development and renegeration in a variety of species, such as humans, mice, goats. We aim to provide theoretical mechanisms for the artificial interference of villus growth cycle, and new ideas for therapeutic treatment of skin hair follicle- related diseases.

CircERCC6 Positively Regulates the Induced Activation of SHF Stem Cells in Cashmere Goats via the miR-412-3p/BNC2 Axis in an m6A-Dependent Manner.

The cashmere, a kind of nature protein fiber, is one of the main use of cashmere goats. The induced activation of secondary hair follicle (SHF) stem cells by the dermal papilla cell-derived signals is a key biological process for the morphogenesis and growth of cashmere fiber in cashmere goats. Previously, the circRNA-ERCC6 (circERCC6) was identified from cashmere goat SHFs; however, its biological significance is unclear in the SHF physiology process of cashmere goats. In this study, we found that circERCC6 exhibited significantly higher expression at anagen SHF bulge compared with the counterpart of telogen and harbored three m6A modified sites (named m6A-685, m6A-862, and m6A-995) through methylation immunoprecipitation using a real-time quantitative polymerase chain reaction (Me-RIP-qPCR) technique. The knockdown experiments of circERCC6 in SHF stem cells showed that circERCC6 positively regulates the induced activation of SHF stem cells in cashmere goats. Through a dual-luciferase reporter assay, we demonstrated that m6A-modified circERCC6 (m6A-circERCC6) sponged miR-412-3p to upregulate the expression of BNC2 mRNA in SHFstem cells. Through m6A-deficient mutant assay in circERCC6 knockdown SHF stem cells, we further showed that m6A modification within circERCC6 is required to mediate the miR-412-3p/BNC2 axis to finally promote the proper induced activation of SHF stem cells in cashmere goats.

N6-Methyladenosine modification (m6A) of circRNA-ZNF638 contributes to the induced activation of SHF stem cells through miR-361-5p/Wnt5a axis in cashmere goats.

OBJECTIVE: The objective of this study was to investigate the effects of N6-Methyladenosine modification-circRNA-zinc finger protein 638 (m6A-circRNA-ZNF638) on the induced activation of secondary hair follicle (SHF) stem cells with its potential mechanisms in cashmere goats. METHODS: The m6A modification of ZNF638 was analyzed using methylation immunoprecipitation with real-time quantitative polymerase chain reaction technique in SHF stem cells. The effects of circRNA-ZNF638 on the induced activation of SHF stem cells in m6A dependence were evaluated through the overexpression of circRNA-ZNF638/its m6Adeficient mutants in circRNA-ZNF638 knockdown SHF stem cells. The competitive binding of miR-361-5p to circRNA-ZNF638/Wnt5a 3'- untranslated region was analyzed through Dual-luciferase reporter assay. RESULTS: The m6A-circRNA-ZNF638 had significantly higher transcription at anagen SHF bulge of cashmere goats compared with that at telogen, as well as it positively regulated the induced activation of SHF-stem cells in cashmere goats. Mechanismly, m6A-circRNA-ZNF638 sponged miR-361-5p to heighten the transcriptional expression of Wnt5a gene in SHFstem cells. We further demonstrated that the internal m6A modification within circRNAZNF638 is required for mediating the miR-361-5p/Wnt5a pathway to regulate the induced activation of SHF stem cells through an introducing of m6A-deficient mutant of circRNAZNF638. CONCLUSION: The circRNA-ZNF638 contributes the proper induced activation of SHF-stem cells in cashmere goats in m6A-dependent manner through miR-361-5p/Wnt5a axis.

CircRNA-0100 positively regulates the differentiation of cashmere goat SHF-SCs into hair follicle lineage via sequestering miR-153-3p to heighten the KLF5 expression.

Circular RNAs (circRNAs) have stable structures, being a covalently closed loop without 5 ' and 3 ' free ends. They can function as "miRNA sponges" in regulating the expression of their target genes. It was thought that circRNAs are involved in the development of the secondary hair follicle (SHF) in cashmere goats. In our previous investigation, a new circRNA named circRNA-0100 was identified from the SHF of cashmere goats, but its function is unknown. In this work, we found that circRNA-0100 exhibited significantly higher expression at anagen SHF bulge than its counterpart at telogen in cashmere goats. Based on the use of both overexpression and siRNA interference assays, our data indicated that circRNA-0100 promoted the differentiation of cashmere goat SHF stem cells (SHF-SCs) into hair follicle lineage, which was evaluated by analyzing the transcriptional level changes of six indicator genes in SHF-SCs of cashmere goats. Using the RNA pull-down technique, we showed that circRNA-0100 served as "molecular sponges" of miR-153-3p in SHF-SCs. Through the use of dual-luciferase reporter assays, our data indicated that circRNA-0100 positively regulated the transcriptional expression of the KLF5 gene via the miR-153-3p-mediated pathway. Ultimately, we showed that circRNA-0100 promoted the differentiation of SHF-SCs into hair lineage, which might be achieved via sequestering miR-153-3p to heighten the KLF5 expression in SHF-SCs of cashmere goats. Our results provide novel scientific evidence for revealing the potential molecular regulatory mechanisms on the differentiation of SHF-SCs into hair lineage in cashmere goats.

Identification of the Key Genes Associated with Different Hair Types in the Inner Mongolia Cashmere Goat.

The Inner Mongolia cashmere goat is an excellent local breed in China. According to the characteristics of wool quilts, the Inner Mongolia cashmere goat can be divided into three types: a long-hair type (hair length of >22 cm), a short-hair type (hair length of ≤13 cm), and an intermediate type (hair length of >13 cm and ≤22 cm). It is found that hair length has a certain reference value for the indirect selection of other important economic traits of cashmere. In order to explore the molecular mechanisms and related regulatory genes of the different hair types, a weighted gene coexpression network analysis (WGCNA) was carried out on the gene expression data and phenotypic data of 12-month-old Inner Mongolia cashmere goats with a long-hair type (LHG) and a short-hair type (SHG) to explore the coexpression modules related to different coat types and nine candidate genes, and detect the relative expression of key candidate genes. The results showed that the WGCNA divided these genes into 19 coexpression modules and found that there was a strong correlation between one module and different hair types. The expression trends of this module's genes were different in the two hair types, with high expression in the LHG and low expression in the SHG. GO functions are mainly concentrated in cellular components, including intermediate filaments (GO:0005882), intermediate filament cytoskeletons (GO:0045111), and cytoskeletal parts (GO:0044430). The KEGG pathway is mainly enriched in arginine as well as proline metabolism (chx00330) and the MAPK signaling pathway (chx04010). The candidate genes of the different hair types, including the KRT39, KRT74, LOC100861184, LOC102177231, LOC102178767, LOC102179881, LOC106503203, LOC108638293, and LOC108638298 genes, were screened. Through qRT-PCR, it was found that there were significant differences in these candidate genes between the two hair types, and most of them had a significant positive correlation with hair length. It was preliminarily inferred that these candidate genes could regulate the different hair types of cashmere goats and provide molecular markers for hair growth.

Identification of Genes Related to Hair Follicle Cycle Development in Inner Mongolia Cashmere Goat by WGCNA.

Cashmere goat from Inner Mongolia is an excellent local breed in China, and the related cashmere product is a kind of precious textile raw material with high price. Cashmere is generated from secondary hair follicles, which has obvious annual periodicity and includes three different stages: anagen, catagen, and telogen. Therefore, we investigated skin transcriptome data for 12 months using weighted gene co-expression network analysis (WGCNA) to explore essential modules, pathways, and genes responsible for the periodic growth and development of secondary hair follicles. A total of 17 co-expression modules were discovered by WGCNA, and there is a strong correlation between steelblue module and month (0.65, p = 3E-09), anagen (0.52, p = 1E-05), telogen (-0.6, p = 8E-08). Gene expression was generally high during late anagen to catagen (June to December), while expression was downregulated from telogen to early anagen (January-May), which is similar to the growth rule of hair follicle cycle. KEGG pathway enrichment analyses of the genes of steelblue module indicated that genes are mainly enriched in Cell cycle, Wnt signaling pathway, p53 signaling pathway and other important signal pathways. These genes were also significantly enriched in GO functional annotation of the cell cycle, microtubule movement, microtubule binding, tubulin binding, and so on. Ten genes (WIF1, WNT11, BAMBI, FZD10, NKD1, LEF1, CCND3, E2F3, CDC6, and CDC25A) were selected from these modules, and further identified as candidate biomarkers to regulate periodic development of hair follicles using qRT-PCR. The Wnt signaling pathway and Cell cycle play an important role in the periodic development of hair follicles. Ten genes were identified as essential functional molecules related to periodic development of hair follicle. These findings laid a foundation for understanding molecular mechanisms in biological functions such as hair follicle development and hair growth in cashmere goats.

Identification and Molecular Analysis of m6A-circRNAs from Cashmere Goat Reveal Their Integrated Regulatory Network and Putative Functions in Secondary Hair Follicle during Anagen Stage.

N6-methyladenosine (m6A) is the most abundant modification in linear RNA molecules. Over the last few years, interestingly, many circRNA molecules are also found to have extensive m6A modification sites with temporal and spatial specific expression patterns. To date, however, little information is available concerning the expression profiling and functional regulatory characteristics of m6A modified circRNAs (m6A-circRNAs) in secondary hair follicles (SHFs) of cashmere goats. In this study, a total of fifteen m6A-circRNAs were identified and characterized in the skin tissue of cashmere goats. Of these, six m6A-circRNAs were revealed to have significantly higher expression in skin at anagen compared with those at telogen. The constructed ceRNA network indicated a complicated regulatory relationship of the six anagen up-regulated m6A-circRNAs through miRNA mediated pathways. Several signaling pathways implicated in the physiological processes of hair follicles were enriched based on the potential regulatory genes of the six anagen up-regulated m6A-circRNAs, such as TGF-beta, axon guidance, ribosome, and stem cell pluripotency regulatory pathways, suggesting the analyzed m6A-circRNAs might be essentially involved in SHF development and cashmere growth in cashmere goats. Further, we showed that four m6A-circRNAs had highly similar expression trends to their host genes in SHFs of cashmere goats including m6A-circRNA-ZNF638, -TULP4, -DNAJB6, and -CAT. However, the expression patterns of two m6A-circRNAs (m6A-circRNA-STAM2 and -CAAP1) were inconsistent with the linear RNAs from their host genes in the SHFs of cashmere goats. These results provide novel information for eluci-dating the biological function and regulatory characteristics of the m6A-circRNAs in SHF development and cashmere growth in goats.

The regulation mechanism of different hair types in inner Mongolia cashmere goat based on PI3K-AKT pathway and FGF21.

There is genetic diversity of hair types in the Inner Mongolia cashmere goat population. Previous studies have found that fibroblast growth factor 21 (FGF21) and PI3K-AKT signal pathways may be related to different hair types in Inner Mongolia cashmere goats. Therefore, the purpose of this study was to explore the effects of the PI3K-AKT signal pathway on different hair types, the expression of mRNA and protein expression sites of FGF21 in the hair follicles of cashmere goats with different hair types, so as to lay a foundation for understanding the molecular mechanism of different hair types and the role of skin hair follicle development. In this experiment, the skin tissues of long hair type (LHG) and short hair type (SHG) of Inner Mongolia cashmere goat were collected in three key periods of secondary hair follicle growth, namely, anagen (September), catagen (December), and telogen (March). The relative expression of FGF21 and PI3K-AKT signal pathway candidate gene mRNA in different periods and different hair types was detected by real-time fluorescence quantitative technique (qRT-PCR), and the expression site of FGF21 protein was located by immunohistochemical technique. Through qRT-PCR, it was found that the relative expression of FGF21, FGFR1, AKT3, BRCA1, PKN3, SPP1, and GNG4 was significantly different between LHG and SHG. The expression of FGF21 in the skin of LHG was significantly higher than that of SHG in the three periods. Through immunohistochemical test, it was found that FGF21 protein was mainly expressed in primary hair follicle connective tissue sheath, primary hair follicle outer root sheath, secondary hair follicle outer root sheath, and sebaceous glands. It was also found that the expression of LHG skin tissue in the outer root sheath of primary hair follicles was higher than that of SHG in three periods. In summary, it is suggested that the PI3K-AKT signal pathway may play an important role in the formation of different hair types in Inner Mongolia cashmere goats.

There is genetic diversity of hair types in Inner Mongolia cashmere goat population. The purpose of this study was to explore the effects of the PI3K-AKT signal pathway on different hair types, the expression of mRNA and protein expression sites of FGF21 in the hair follicles of cashmere goats with different hair types, so as to lay a foundation for understanding the molecular mechanism of different hair types. It was found that the relative expression of FGF21, FGFR1, AKT3, BRCA1, PKN3, SPP1, and GNG4 was significantly different between LHG and SHG. It was found that FGF21 protein was mainly expressed in primary hair follicle connective tissue sheath, primary hair follicle outer root sheath, secondary hair follicle outer root sheath, and sebaceous glands. It was also found that the expression of LHG skin tissue in the outer root sheath of primary hair follicles was higher than that of SHG in three periods. So, it is suggested that the PI3K-AKT signal pathway and FGF21 may play an important role in the formation of different hair types in Inner Mongolia cashmere goats.

Single-Cell Sequencing Reveals Differential Cell Types in Skin Tissues of Liaoning Cashmere Goats and Key Genes Related Potentially to the Fineness of Cashmere Fiber.

Cashmere fineness is one of the important factors determining cashmere quality; however, our understanding of the regulation of cashmere fineness at the cellular level is limited. Here, we used single-cell RNA sequencing and computational models to identify 13 skin cell types in Liaoning cashmere goats. We also analyzed the molecular changes in the development process by cell trajectory analysis and revealed the maturation process in the gene expression profile in Liaoning cashmere goats. Weighted gene co-expression network analysis explored hub genes in cell clusters related to cashmere formation. Secondary hair follicle dermal papilla cells (SDPCs) play an important role in the growth and density of cashmere. ACTA2, a marker gene of SDPCs, was selected for immunofluorescence (IF) and Western blot (WB) verification. Our results indicate that ACTA2 is mainly expressed in SDPCs, and WB results show different expression levels. COL1A1 is a highly expressed gene in SDPCs, which was verified by IF and WB. We then selected CXCL8 of SDPCs to verify and prove the differential expression in the coarse and fine types of Liaoning cashmere goats. Therefore, the CXCL8 gene may regulate cashmere fineness. These genes may be involved in regulating the fineness of cashmere in goat SDPCs; our research provides new insights into the mechanism of cashmere growth and fineness regulation by cells.

Draft Genome of the European Mouflon (Ovis orientalis musimon).

Mouflon (Ovis orientalis) with its huge and beautiful horns is considered as one of the ancestors of domesticated sheep. The European mouflon (Ovis orientalis musimon) is in the Asiatic mouflon (O. orientalis) clade. In order to provide novel genome information for mouflon, moreover promote genetic analysis of genus Ovis both domestic and wild, we propose to sequence the mouflon genome. We assembled the highly heterozygous mouflon genome based on Illumina HiSeq platform using the next-generation sequencing technology. Finally, the draft genome we accessed approximately 2.69 Gb (42.15% GC), while N50 sizes of contig and scaffold are 110.1 kb and 10.4 Mb, respectively. The contiguity of this assembly is obviously better than earlier versions. Further analyses predicted 20,814 protein-coding genes in the mouflon genome and 12,390 shared gene families among bovine species. It is estimated that the divergence time between O. orientalis musimon and Ovis aries was 7.6 million years ago. The draft mouflon genome assembly will provide data support and theoretical basis for various investigations of the genus Ovis species in future.

m6A Methylation Analysis of Differentially Expressed Genes in Skin Tissues of Coarse and Fine Type Liaoning Cashmere Goats.

N6-methyladenosine (m6A) is the most common internal modification in mRNAs of all higher eukaryotes. Here we perform two high-throughput sequencing methods, m6A-modified RNA immunoprecipitation sequence (MeRIP-seq) and RNA sequence (RNA-seq) to identify key genes with m6A modification in cashmere fiber growth. A total of 9,085 m6A sites were differentially RNA m6A methylated as reported from by MeRIP-seq, including 7,170 upregulated and 1,915 downregulated. In addition, by comparing m6A-modified genes between the fine-type Liaoning cashmere goat (FT-LCG) and coarse-type Liaoning Cashmere Goat (CT-LCG) skin samples, we obtain 1,170 differentially expressed genes. In order to identify the differently methylated genes related to cashmere fiber growth, 19 genes were selected to validate by performing qRT-PCR in FT-LCG and CT-LCG. In addition, GO enrichment analysis shows that differently methylated genes are mainly involved in keratin filament and intermediate filament. These findings provide a theoretical basis for future research on the function of m6A modification during the growth of cashmere fiber.

Transcriptomic analysis reveals critical genes for the hair follicle of Inner Mongolia cashmere goat from catagen to telogen.

There are two main types of hair follicle in Inner Mongolia Cashmere goats, the primary hair follicle (PHF) producing hair fibers and the secondary hair follicle (SHF) producing cashmere fibers. Of both fibers from cashmere-bearing goats in Aerbasi, Inner Mongolia, the timing of cyclical phases for the cashmere have been well clarified but hair fibers have been less noticeable. Herein, we evaluated transcriptome of PHF and SHF from the same three goats in Aerbasi at the catagen- and telogen phase of cashmere growth. We totally found 1977 DEGs between PHFs at the telogen and catagen phases of SHF, 1199 DEGs between telogen- and catagen SHF, 2629 DEGs between PHF at the catagen phase of SHF and catagen SHF, and 755 DEGs between PHF at the telogen phase of SHF and telogen SHF. By analyzing gene functions based on GO and KEGG database, we found that the DEGs have functions in muscle contraction and muscle filament sliding between catagen- and telogen SHF, indicating that arrector pilli muscles might play a role on the transition from catagen to telogen. Moreover, considering that the enriched GO and KEGG categories of the DEGs between PHF and SHF, we suggested that part of PHF might rest in their own anagen phase when SHF are at catagen, but PHF might enter into the telogen phase at SHF's telogen. Finally, we highly recommended the several potential genes acting as the regulators of the transition between growth phases including IL17RB and eight members of ZNF. These results provide insight into molecular mechanisms on the transition of SHF from catagen to telogen together with PHF's growth situation at SHF's catagen and telogen in Inner Mongolia Cashmere goats.

The complete mitochondrial DNA analyses of Inner Mongolia Cashmere goat type of Erlangshan species.

Inner Mongolia Cashmere goat is evolved from Mongolia goat in long-term breeding, as a source of high-quality cashmere, Inner Mongolia Cashmere goat type of Erlangshan is famous in the world. By using molecular biology techniques, we sequenced the complete mitochondrial DNA of Erlangshan Cashmere goat. The complete length of sequence is 16,640 bp, containing 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes and a control region (D-loop region). Our mitogenome will enlarge genomic information for further studies on evolution and domestication of Inner Mongolia Cashmere goat, enhance germplasm conservation and breeding programs of domestic goat.

Identification of selection signals by large-scale whole-genome resequencing of cashmere goats.

Inner Mongolia and Liaoning cashmere goats are two outstanding Chinese multipurpose breeds that adapt well to the semi-arid temperate grassland. These two breeds are characterized by their soft cashmere fibers, thus making them great models to identify genomic regions that are associated with cashmere fiber traits. Whole-genome sequencing of 70 cashmere goats produced more than 5.52 million single-nucleotide polymorphisms and 710,600 short insertions and deletions. Further analysis of these genetic variants showed some population-specific molecular markers for the two cashmere goat breeds that are otherwise phenotypically similar. By analyzing F ST and θπ outlier values, we identified 135 genomic regions that were associated with cashmere fiber traits within the cashmere goat populations. These selected genomic regions contained genes, which are potential involved in the production of cashmere fiber, such as FGF5, SGK3, IGFBP7, OXTR, and ROCK1. Gene ontology enrichment analysis of identified short insertions and deletions also showed enrichment in keratinocyte differentiation and epidermal cell differentiation. These findings demonstrate that this genomic resource will facilitate the breeding of cashmere goat and other Capra species in future.

Genome-wide Target Enrichment-aided Chip Design: a 66 K SNP Chip for Cashmere Goat.

Compared with the commercially available single nucleotide polymorphism (SNP) chip based on the Bead Chip technology, the solution hybrid selection (SHS)-based target enrichment SNP chip is not only design-flexible, but also cost-effective for genotype sequencing. In this study, we propose to design an animal SNP chip using the SHS-based target enrichment strategy for the first time. As an update to the international collaboration on goat research, a 66 K SNP chip for cashmere goat was created from the whole-genome sequencing data of 73 individuals. Verification of this 66 K SNP chip with the whole-genome sequencing data of 436 cashmere goats showed that the SNP call rates was between 95.3% and 99.8%. The average sequencing depth for target SNPs were 40X. The capture regions were shown to be 200 bp that flank target SNPs. This chip was further tested in a genome-wide association analysis of cashmere fineness (fiber diameter). Several top hit loci were found marginally associated with signaling pathways involved in hair growth. These results demonstrate that the 66 K SNP chip is a useful tool in the genomic analyses of cashmere goats. The successful chip design shows that the SHS-based target enrichment strategy could be applied to SNP chip design in other species.