Effects of Heat Stress on Goat Production and Mitigating Strategies: A Review.

Goats, versatile creatures selectively bred for various purposes, have become pivotal in shaping the socioeconomic landscape, particularly in rural and economically challenged areas. Their remarkable ability to withstand and adapt to extreme heat has proven invaluable, allowing them to flourish and reproduce in even the harshest climates on Earth. Goat farming has emerged as a reliable and sustainable solution for securing food resources. However, despite its significance, the goat-producing industry has received less attention than other ruminants. Despite goats' inherent resilience to heat, their productivity and reproductive performance suffer under high ambient temperatures, leading to heat stress. This presents a significant challenge for goat production, necessitating a comprehensive multidisciplinary approach to mitigating the adverse effects of heat stress. This review aims to explore the diverse impacts of heat stress on goats and propose effective measures to address the sector's challenges. By understanding and addressing these issues, we can enhance the resilience and sustainability of goat farming, ensuring its continued contribution to food security and socioeconomic development.

PCB169 exposure aggravated the development of non-alcoholic fatty liver in high-fat diet-induced male C57BL/6 mice.

Polychlorinated biphenyls (PCBs) are lipophilic environmental toxicants. Epidemiological studies have established a link between PCBs and both metabolic syndrome and nonalcoholic fatty liver disease (NAFLD). Multiple studies have reported that exposure to both PCB156 and PCB126 among the 12 dioxin-like PCBs leads to the development of NAFLD. However, studies to elucidate whether PCB169 induces the development of NAFLD by constructing in vivo models have not been reported. Therefore, we evaluated the effects of exposure to PCB169 (5 mg/kg-bw) on hepatic lipid metabolism in C57BL/6 mice from control diet and high-fat diet cohorts. The results showed that PCB169 exposure reduced body weight and intraperitoneal fat mass in mice on the control diet, but the liver lipid levels were significantly increased, exacerbating NAFLD in mice on a high-fat diet. Through transcriptomics studies, it was found that PCB169 exposure induced significant up-regulation of Pparγ, Fasn, and Aacs genes involved in hepatic lipogenesis, as well as remarkable up-regulation of Hmgcr, Lss, and Sqle genes involved in cholesterol synthesis. Additionally, there was notable down-regulation of Pparα and Cpt1 genes involved in lipid ß-oxidation, leading to abnormal lipid accumulation in the liver. In addition, we found that PCB169 exposure significantly activated the Arachidonic acid metabolism, PPAR signaling pathway, Metabolism of xenobiotics by cytochrome P450, and Retinol metabolism pathways, and so on. Our study suggests that PCB169 can modify gene expression related to lipid metabolism, augument lipid accumulation in the liver, and further contribute to the development of NAFLD, thereby revealing the detrimental effects associated with PCB exposure on animal growth and metabolism.

Replacing alfalfa hay with amaranth hay: effects on production performance, rumen fermentation, nutrient digestibility and antioxidant ability in dairy cow.

OBJECTIVE: The aim of this research was to explore the effects of dietary substitution of alfalfa hay by amaranth hay on production performance, rumen fermentation, nutrient digestibility, serum biochemical parameters and antioxidant ability in dairy cows. METHODS: A total of 45 healthy Holstein cows with same parity and similar milk yield and body weight were randomly divided into 3 groups: control diet without amaranth hay (CON) or 50% and 100% alfalfa hay replaced by an equal amount of amaranth hay (dry matter basis, AH1 and AH2, respectively). All the cows were fed regularly 3 times a day at 06:30, 14:30, and 22:30 and had free access to water. The experiment lasted for 60 d. RESULTS: The dry matter intake of CON and AH1 groups was higher (p<0.05) than that of AH2 group. Compared with AH1 group, the milk yield of AH2 group was reduced (p<0.05). Moreover, dietary substitution of alfalfa hay by amaranth hay increased (p<0.05) milk fat, ammonia nitrogen and acetate concentrations. However, the crude protein digestibility of AH2 group was lower (p<0.05) than that of CON group, while an opposite tendency of serum urea nitrogen was found between two groups. The neutral detergent fiber digestibility of AH1 group was increased (p<0.05) when compared to AH2 group. Amaranth hay treatment increased (p<0.05) the serum concentration of glutathione peroxidase in dairy cows. Compared with CON group, the malonaldehyde activity of AH1 group was decreased (p<0.05). CONCLUSION: Dietary replacing alfalfa hay with amaranth hay (50% ratio) in dairy cows did not affect production performance but improved their antioxidant ability.

Grape Seed Extract as a Feed Additive Improves the Growth Performance, Ruminal Fermentation and Immunity of Weaned Beef Calves.

The purpose of this research was to evaluate effects of grape seed extract (Gse) supplementation on the growth performance; ruminal fermentation; nutrient digestibility; and serum biochemical, antioxidative, and immune parameters of weaned beef calves. A total of 30 Simmental crossbred male calves with similar age and body weight were randomly allocated to two groups: a control group with no Gse (CON) and a Gse supplementation group (GSE) (4 g/d Gse per animal). The results show that, compared with the CON group, the average daily gain significantly increased (p = 0.043) in the GSE group. The ruminal contents of microbial protein and butyrate in GSE group were higher (p < 0.05) than those in the CON group. Additionally, calves fed Gse displayed increased (p < 0.05) dry matter and neutral detergent fiber digestibility. Moreover, the serum concentrations of triglyceride, catalase, superoxide dismutase, immunoglobulin G and immunoglobulin M were higher (p < 0.05) in the GSE group than those in the CON group. However, opposite tendencies of non-esterified fatty acid, malondialdehyde, tumor necrosis factor-α and interleukin-6 were found between the two groups. Overall, the supplementation of Gse can improve ruminal fermentation, nutrient digestibility, antioxidant ability, and immunity, as well as promoting the healthy growth of weaned cross-breed beef calves.

An atlas of CNV maps in cattle, goat and sheep.

Copy number variation (CNV) is the most prevalent type of genetic structural variation that has been recognized as an important source of phenotypic variation in humans, animals and plants. However, the mechanisms underlying the evolution of CNVs and their function in natural or artificial selection remain unknown. Here, we generated CNV region (CNVR) datasets which were diverged or shared among cattle, goat, and sheep, including 886 individuals from 171 diverse populations. Using 9 environmental factors for genome-wide association study (GWAS), we identified a series of candidate CNVRs, including genes relating to immunity, tick resistance, multi-drug resistance, and muscle development. The number of CNVRs shared between species is significantly higher than expected (P<0.00001), and these CNVRs may be more persist than the single nucleotide polymorphisms (SNPs) shared between species. We also identified genomic regions under long-term balancing selection and uncovered the potential diversity of the selected CNVRs close to the important functional genes. This study provides the evidence that balancing selection might be more common in mammals than previously considered, and might play an important role in the daily activities of these ruminant species.

Base pair editing in goat: nonsense codon introgression into FGF5 results in longer hair.

The ability to alter single bases without homology directed repair (HDR) of double-strand breaks provides a potential solution for editing livestock genomes for economic traits, which are often multigenic. Progress toward multiplex editing in large animals has been hampered by the costly inefficiencies of HDR via microinjection of in vitro manipulated embryos. Here, we designed sgRNAs to induce nonsense codons (C-to-T transitions) at four target sites in caprine FGF5, which is a crucial regulator of hair length in mammals. Initial transfections of the third generation Base Editor (BE3) plasmid and four different sgRNAs into caprine fibroblasts were ineffective in altering FGF5. In contrast, all five progenies produced from microinjected single-cell embryos had alleles with a targeted nonsense mutation. The effectiveness of BE3 to make single base changes varied considerably based on sgRNA design. In addition, the rate of mosaicism differed between animals, target sites, and tissue type. The phenotypic effects on hair fiber were characterized by hematoxylin and eosin, immunofluorescence staining, and western blotting. Differences in morphology were detectable, even though mosaicism was probably affecting the levels of FGF5 expression. PCR amplicon and whole-genome resequencing analyses for off-target changes caused by BE3 were low at a genome-wide scale. This study provided the first evidence of base editing in large mammals produced from microinjected single-cell embryos. Our results support further optimization of BEs for introgressing complex human disease alleles into large animal models, to evaluate potential genetic improvement of complex health and production traits in a single generation.

CRISPR/Cas9-mediated VDR knockout plays an essential role in the growth of dermal papilla cells through enhanced relative genes.

BACKGROUND: Hair follicles in cashmere goats are divided into primary and secondary hair follicles (HFs). HF development, which determines the morphological structure, is regulated by a large number of vital genes; however, the key functional genes and their interaction networks are still unclear. Although the vitamin D receptor (VDR) is related to cashmere goat HF formation, its precise effects are largely unknown. In the present study, we verified the functions of key genes identified in previous studies using hair dermal papilla (DP) cells as an experimental model. Furthermore, we used CRISPR/Cas9 technology to modify the VDR in DP cells to dissect the molecular mechanism underlying HF formation in cashmere goats. RESULTS: The VDR expression levels in nine tissues of Shaanbei white cashmere goats differed significantly between embryonic day 60 (E60) and embryonic day 120 (E120). At E120, VDR expression was highest in the skin. At the newborn and E120 stages, the VDR protein was highly expressed in the root sheath and hair ball region of Shaanbei cashmere goats. We cloned the complete CDS of VDR in the Shaanbei white cashmere goat and constructed a VDR-deficient DP cell model by CRISPR/Cas9. Heterozygous and homozygous mutant DP cells were produced. The growth rate of mutant DP cells was significantly lower than that of wild-type DP cells (P < 0.05) and VDR mRNA levels in DP cells decreased significantly after VDR knockdown (P < 0.05). Further, the expression levels of VGF, Noggin, Lef1, and ß-catenin were significantly downregulated (P < 0.05). CONCLUSIONS: Our results indicated that VDR has a vital role in DP cells, and that its effects are mediated by Wnt and BMP4 signaling.

Allele-specific expression and alternative splicing in horse×donkey and cattle×yak hybrids.

Divergence of gene expression and alternative splicing is a crucial driving force in the evolution of species; to date, however the molecular mechanism remains unclear. Hybrids of closely related species provide a suitable model to analyze allele-specific expression (ASE) and allele-specific alternative splicing (ASS). Analysis of ASE and ASS can uncover the differences in cis-regulatory elements between closely related species, while eliminating interference of trans-regulatory elements. Here, we provide a detailed characterization of ASE and ASS from 19 and 10 transcriptome datasets across five tissues from reciprocal-cross hybrids of horse×donkey (mule/hinny) and cattle×yak (dzo), respectively. Results showed that 4.8%-8.7% and 10.8%-16.7% of genes exhibited ASE and ASS, respectively. Notably, lncRNAs and pseudogenes were more likely to show ASE than protein-coding genes. In addition, genes showing ASE and ASS in mule/hinny were found to be involved in the regulation of muscle strength, whereas those of dzo were involved in high-altitude adaptation. In conclusion, our study demonstrated that exploration of genes showing ASE and ASS in hybrids of closely related species is feasible for species evolution research.

Synchronous profiling and analysis of mRNAs and ncRNAs in the dermal papilla cells from cashmere goats.

BACKGROUND: Dermal papilla cells (DPCs), the "signaling center" of hair follicle (HF), delicately master continual growth of hair in mammals including cashmere, the fine fiber annually produced by secondary HF embedded in cashmere goat skins. Such unparalleled capacity bases on their exquisite character in instructing the cellular activity of hair-forming keratinocytes via secreting numerous molecular signals. Past studies suggested microRNA (miRNAs) and long non-coding RNAs (lncRNAs) play essential roles in a wide variety of biological process, including HF cycling. However, their roles and related molecular mechanisms in modulating DPCs secretory activities are still poorly understood. RESULTS: Here, we separately cultivated DPCs and their functionally and morphologically distinct dermal fibroblasts (DFs) from cashmere goat skins at anagen. With the advantage of high throughput RNA-seq, we synchronously identified 2540 lncRNAs and 536 miRNAs from two types of cellular samples at 4th passages. Compared with DFs, 1286 mRNAs, 18 lncRNAs, and 42 miRNAs were upregulated, while 1254 mRNAs, 53 lncRNAs and 44 miRNAs were downregulated in DPCs. Through overlapping with mice data, we ultimately defined 25 core signatures of DPCs, including HOXC8 and RSPO1, two crucial activators for hair follicle stem cells (HFSCs). Subsequently, we emphatically investigated the impacts of miRNAs and lncRNAs (cis- and trans- acting) on the genes, indicating that ncRNAs extensively exert negative and positive effects on their expressions. Furthermore, we screened lncRNAs acting as competing endogenous RNAs (ceRNAs) to sponge miRNAs and relief their repressive effects on targeted genes, and constructed related lncRNAs-miRNAs-HOXC8/RSPO1 interactive lines using bioinformatic tools. As a result, XR_310320.3-chi-miR-144-5p-HOXC8, XR_311077.2-novel_624-RSPO1 and others lines appeared, displaying that lncRNAs might serve as ceRNAs to indirectly adjust HFSCs status in hair growth. CONCLUSION: The present study provides an unprecedented inventory of lncRNAs, miRNAs and mRNAs in goat DPCs and DFs. We also exhibit some miRNAs and lncRNAs potentially participate in the modulation of HFSCs activation via delicately adjusting core signatures of DPCs. Our report shines new light on the latent roles and underlying molecular mechanisms of ncRNAs on hair growth.

Whole-genome bisulfite sequencing of goat skins identifies signatures associated with hair cycling.

BACKGROUND: Hair follicles (HFs), upon development, undergo repetitive cycles of growth (anagen), regression (catagen), and rest (telogen). The transition between the stages is determined by multiple molecular signals, including DNA methylation, which plays important roles in mammalian cellular identity and is essential for the development of HFs. Secondary hair follicles (SHFs) in cashmere goat exhibit classic cyclic hair development, and little has been done on a genome-wide scale to examine potentially methylated genes involved in the hair cyclic transition. RESULTS: Genome-wide DNA methylation profiles between skin tissues sampled during the anagen and telogen stages in cashmere goats were investigated using whole-genome bisulfite sequencing (WGBS). The methylation status was observed to be higher in the skin samples with HFs in the telogen than those in the anagen stage. A total of 1311 differentially methylated regions (DMRs) were identified between the two groups, which contained 493 fully annotated DMR-related genes (DMGs) (269 Hyper- DMGs and 224 Hypo-DMGs). Furthermore, a significant over-representation of the functional categories for DMGs related to immune response and intercellular crosstalk during hair cycling was observed. By integrating DNA methylation and mRNA expression data, we revealed that four genes (FMN1, PCOLCE, SPTLC3, and COL5A1) are crucial factors for elucidating epigenetic mechanisms contributing to the telogen-to-anagen transition. CONCLUSION: Our study provided systematic methylome maps pertaining to the hair cycling stages (anagen vs telogen) at a single-base resolution, and revealed stage-specific methylation loci during cashmere growth or quiescence. Furthermore, we identified epigenetically regulated genes that are potentially involved in HF development and growth in cashmere goats, and likely in other mammal species.

Effects of all-trans retinoic acid on goat dermal papilla cells cultured in vitro

Background: All-trans retinoic acid (ATRA), a vitamin A-derived active metabolite, exerts important functions in hair biology. Previous studies indicated that excess ATRA hampered hair follicle morphogenesis and cyclic regeneration in adulthood, but other studies stated that ATRA promoted hair growth. Dermal papilla (DP), a cluster of specialized fibroblasts, plays pivotal roles in controlling development and regeneration of hair follicle. Several lines of evidence indicated that DP might be the target cells of ATRA in the hair follicle. To confirm this hypothesis, the present study was performed to explore the biological effects of ATRA on goat dermal papilla cells (DPCs) and clarify the roles of ATRA in hair biology. Results: Our experimental results indicated that key signaling transducers of ATRA were dynamically expressed in distinct stages of goat cashmere growth cycle, and high-dose ATRA treatment (10-5 M) significantly impaired the viability of goat DPCs and lowered the ratio of proliferating cells. Otherwise, goat DPCs were stimulated to enter apoptosis and their cell cycle progression was severely blocked by ATRA. Moreover, the expression of fibroblast growth factor 7 (Fgf7), one of the potent hair growth stimulators secreted by DPCs, was transcriptionally repressed following ATRA treatment. Conclusion: DPCs are the targets of ATRA in the hair follicle, and ATRA negatively regulates hair growth by the targeted suppression of cell viability and growth factor expression of goat DPCs. Through these observations, we offer a new mechanistic insight into the roles of ATRA in hair biology.

Integrative analysis reveals ncRNA-mediated molecular regulatory network driving secondary hair follicle regression in cashmere goats.

BACKGROUND: Cashmere is a keratinized product derived from the secondary hair follicles (SHFs) of cashmere goat skins. The cashmere fiber stops growing following the transition from the actively proliferating anagen stage to the apoptosis-driven catagen stage. However, little is known regarding the molecular mechanisms responsible for the occurrence of apoptosis in SHFs, especially as pertains to the role of non-coding RNAs (ncRNAs) and their interactions with other molecules. Hair follicle (HF) degeneration is caused by localized apoptosis in the skin, while anti-apoptosis pathways may coexist in adjacent HFs. Thus, elucidating the molecular interactions responsible for apoptosis and anti-apoptosis in the skin will provide insights into HF regression. RESULTS: We used multiple-omics approaches to systematically identify long non-coding RNAs (lncRNAs), microRNAs (miRNAs) and mRNAs expressed in cashmere goat skins in two crucial phases (catagen vs. anagen) of HF growth. Skin samples were collected from three cashmere goats at the anagen (September) and catagen (February) stages, and six lncRNA libraries and six miRNA libraries were constructed for further analysis. We identified 1122 known and 403 novel lncRNAs in the goat skins, 173 of which were differentially expressed between the anagen and catagen stages. We further identified 3500 gene-encoding transcripts that were differentially expressed between these two phases. We also identified 411 known miRNAs and 307 novel miRNAs, including 72 differentially expressed miRNAs. We further investigated the target genes of lncRNAs via both cis- and trans-regulation during HF growth. Our data suggest that lncRNAs and miRNAs act synergistically in the HF growth transition, and the catagen inducer factors (TGFß1 and BDNF) were regulated by miR-873 and lnc108635596 in the lncRNA-miRNA-mRNA networks. CONCLUSION: This study enriches the repertoire of ncRNAs in goats and other mammals, and contributes to a better understanding of the molecular mechanisms of ncRNAs involved in the regulation of HF growth and regression in goats and other hair-producing species.

Comparative proteomic analyses using iTRAQ-labeling provides insights into fiber diversity in sheep and goats.

The structural component of wool and hair fibers, such as keratin-associated proteins (KAPs), has been well described, but the genetic determinants of fiber diameter are largely unknown. Here, we have used an iTRAQ-based proteomic approach to investigate differences in protein abundance among 18 samples from sheep and goats across a diverse range of fibers. We identified proteins with different abundance and are associated with variation in fiber features. Proteins with different abundance are mainly keratin or keratin-associated proteins (KRTAP11-1, KRT6A, KRT38), or are related to hair growth (DSC2, DSG3, EEF2, CALML5, TCHH, SELENBP1) and fatty acid synthesis (FABP4, FABP5). RNA-seq further confirmed the functional importance of the DSC2 gene in the determination of woolly phenotype in goat fibers. This comprehensive analysis of fibers from major fiber-producing animals is the first to provide a list of candidate proteins that are involved in fiber formation. This list will be valuable asset for future studies into the molecular mechanisms that underlie fiber diversity. BIOLOGICAL SIGNIFICANCE: Proteins are the basis for animal-derived hair fibers, yet proteins conferring fiber structure and characteristics in sheep and goats are largely elusive. By examining 27 fibers samples representing 9 fiber types from sheep and goats through the iTRAQ approach, we show a list of differentially abundant proteins that are important to hair structural component, or genes related to hair growth and fatty acid synthesis. RNA-seq further validated the DSC2 gene is key to the woolly/straight hair phenotype in goats.

Efficient generation of goats with defined point mutation (I397V) in GDF9 through CRISPR/Cas9.

The recent emergence of the clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated (Cas) 9 system has attracted significant attention for its potential to improve traits of agricultural importance. However, most applications in livestock species to date have depended on aberrant DNA repair to generate frameshifting indels. Whether this genomic engineering technique involving homology-dependent repair (HDR) can be used to introduce defined point mutations has been less explored. Previously, we reported a G→A point mutation (g.231A>G, p.Val397Ile) in the growth differentiation factor 9 (GDF9) gene that has a large effect on the litter size of cashmere goats. In the present study we report that by co-injecting synthesised RNAs and single-stranded oligo deoxynucleotide (ssODN) donor sequences into goat zygotes, we successfully introduced defined point mutations resulting in single amino acid substitutions in the proteins as expected. The efficiency of this precise single-nucleotide substitution in newborn kids was as high as 24% (4/17), indicating that ssODN-directed HDR via zygote injection is efficient at introducing point mutations in the goat genome. The findings of the present study further highlight the complex genome modifications facilitated by the CRISPR/Cas9 system, which is able to introduce defined point mutations. This represents a significant development for the improvement of reproduction traits in goats, as well as for validating the roles of specific nucleotides in functional genetic elements in large animals.

Comparative transcriptome analysis reveals potentially novel roles of Homeobox genes in adipose deposition in fat-tailed sheep.

Adipose tissues are phenotypically, metabolically and functionally heterogeneous based on the sites of their deposition. Undesirable fat deposits in the body are often detrimental to animal and human health. To unravel the potential underlying mechanisms governing accumulation of adipose tissues in various regions of the body, i.e., subcutaneous (SAT), visceral (VAT) and tail (TAT), we profiled transcriptomes from Tan sheep, a Chinese indigenous breed with notable fat tail using RNA-seq. Upon comparison, we identified a total of 1,058 differentially expressed genes (DEGs) between the three adipose types (218, 324, and 795 in SAT/VAT, SAT/TAT, and VAT/TAT, respectively), from which several known key players were identified that are involved in lipid metabolic process, Wnt signals, Vitamin A metabolism, and transcriptional regulation of adipocyte differentiation. We also found that many elevated genes in VAT were notably enriched for key biological processes such as cytokine secretion, signaling molecule interaction and immune systems. Several developmental genes including HOXC11, HOXC12 and HOXC13, and adipose-expressed genes in the tail region, such as HOTAIR_2, HOTAIR_3 and SP9 were specially highlighted, indicating their strong associations with tail fat development in fat-tailed sheep. Our results provide new insight into exploring the specific fat deposition in tail, also contribute to the understanding of differences between adipose depots.

Integrating miRNA and mRNA Expression Profiling Uncovers miRNAs Underlying Fat Deposition in Sheep.

MicroRNAs (miRNAs) are endogenous, noncoding RNAs that regulate various biological processes including adipogenesis and fat metabolism. Here, we adopted a deep sequencing approach to determine the identity and abundance of miRNAs involved in fat deposition in adipose tissues from fat-tailed (Kazakhstan sheep, KS) and thin-tailed (Tibetan sheep, TS) sheep breeds. By comparing HiSeq data of these two breeds, 539 miRNAs were shared in both breeds, whereas 179 and 97 miRNAs were uniquely expressed in KS and TS, respectively. We also identified 35 miRNAs that are considered to be putative novel miRNAs. The integration of miRNA-mRNA analysis revealed that miRNA-associated targets were mainly involved in the gene ontology (GO) biological processes concerning cellular process and metabolic process, and miRNAs play critical roles in fat deposition through their ability to regulate fundamental pathways. These pathways included the MAPK signaling pathway, FoxO and Wnt signaling pathway, and focal adhesion. Taken together, our results define miRNA expression signatures that may contribute to fat deposition and lipid metabolism in sheep.

Effect of miR-125b on dermal papilla cells of goat secondary hair follicle

Background: MicroRNAs (miRNAs) are endogenous noncoding RNAs that regulate various biological processes. miR-125b is a miRNA that has been reported to be critical for hair follicle (HF) morphogenesis and development. We identified that the expression of miR-125b varies during an individual hair cycle (anagen, catagen, and telogen) in the skin of cashmere goats. We constructed a gain model (by overexpressing miR-125b) and a loss model (by inhibiting endogenous miR-125b) based on dermal papilla cells (DPCs) to further investigate the role of miR-125b in HF cycle. In addition, we used a dual-luciferase system to highlight the predicated target genes of miR-125b. Results: We found that miR-125b affects the expression of FGF5, IGF-1, SHH, TNF-α, MSX2, LEF-1, FGF7, NOGGIN, BMP2, BMP4, TGF-ß1, and ß-catenin. The dual-luciferase assay further validated a direct interaction between miR-125b and FGF5 and TNF-α. Conclusion: miR-125b affects the expression levels of genes related to hair cycle and may also play a critical role in regulating the periodic development of HF.

Whole-genome sequencing of eight goat populations for the detection of selection signatures underlying production and adaptive traits.

The goat (Capra hircus) is one of the first farm animals that have undergone domestication and extensive natural and artificial selection by adapting to various environments, which in turn has resulted in its high level of phenotypic diversity. Here, we generated medium-coverage (9-13×) sequences from eight domesticated goat breeds, representing morphologically or geographically specific populations, to identify genomic regions representing selection signatures. We discovered ~10 million single nucleotide polymorphisms (SNPs) for each breed. By combining two approaches, ZHp and di values, we identified 22 genomic regions that may have contributed to the phenotypes in coat color patterns, body size, cashmere traits, as well as high altitude adaptation in goat populations. Candidate genes underlying strong selection signatures including coloration (ASIP, KITLG, HTT, GNA11, and OSTM1), body size (TBX15, DGCR8, CDC25A, and RDH16), cashmere traits (LHX2, FGF9, and WNT2), and hypoxia adaptation (CDK2, SOCS2, NOXA1, and ENPEP) were identified. We also identified candidate functional SNPs within selected genes that may be important for each trait. Our results demonstrated the potential of using sequence data in identifying genomic regions that are responsible for agriculturally significant phenotypes in goats, which in turn can be used in the selection of goat breeds for environmental adaptation and domestication.

Multiplex gene editing via CRISPR/Cas9 exhibits desirable muscle hypertrophy without detectable off-target effects in sheep.

The CRISPR/Cas9 system provides a flexible approach for genome engineering of genetic loci. Here, we successfully achieved precise gene targeting in sheep by co-injecting one-cell-stage embryos with Cas9 mRNA and RNA guides targeting three genes (MSTN, ASIP, and BCO2). We carefully examined the sgRNAs:Cas9-mediated targeting effects in injected embryos, somatic tissues, as well as gonads via cloning and sequencing. The targeting efficiencies in these three genes were within the range of 27-33% in generated lambs, and that of simultaneously targeting the three genes was 5.6%, which demonstrated that micro-injection of zygotes is an efficient approach for generating gene-modified sheep. Interestingly, we observed that disruption of the MSTN gene resulted in the desired muscle hypertrophy that is characterized by enlarged myofibers, thereby providing the first detailed evidence supporting that gene modifications had occurred at both the genetic and morphological levels. In addition, prescreening for the off-target effect of sgRNAs was performed on fibroblasts before microinjection, to ensure that no detectable off-target mutations from founder animals existed. Our findings suggested that the CRISPR/Cas9 method can be exploited as a powerful tool for livestock improvement by simultaneously targeting multiple genes that are responsible for economically significant traits.

Comparative Transcriptome Analysis of Fetal Skin Reveals Key Genes Related to Hair Follicle Morphogenesis in Cashmere Goats.

Cashmere goat skin contains two types of hair follicles (HF): primary hair follicles (PHF) and secondary hair follicles (SHF). Although multiple genetic determinants associated with HF formation have been identified, the molecules that determine the independent morphogenesis of HF in cashmere goats remain elusive. The growth and development of SHF directly influence the quantity and quality of cashmere production. Here, we report the transcriptome profiling analysis of nine skin samples from cashmere goats using 60- and 120-day-old embryos (E60 and E120, respectively), as well as newborns (NB), through RNA-sequencing (RNA-seq). HF morphological changes indicated that PHF were initiated at E60, with maturation from E120, while differentiation of SHF was identified at E120 until formation of cashmere occurred after birth (NB). The RNA-sequencing analysis generated over 20.6 million clean reads from each mRNA library. The number of differentially expressed genes (DEGs) in E60 vs. E120, E120 vs. NB, and E60 vs. NB were 1,024, 0 and 1,801, respectively, indicating that no significant differences were found at transcriptomic levels between E120 and NB. Key genes including B4GALT4, TNC, a-integrin, and FGFR1, were up-regulated and expressed in HF initiation from E60 to E120, while regulatory genes such as GPRC5D, PAD3, HOXC13, PRR9, VSIG8, LRRC15, LHX2, MSX-2, and FOXN1 were up-regulated and expressed in HF keratinisation and hair shaft differentiation from E120 and NB to E60. Several genes belonging to the KRT and KRTAP gene families were detected throughout the three HF developmental stages. The transcriptional trajectory analyses of all DEGs indicated that immune privilege, glycosaminoglycan biosynthesis, extracellular matrix receptor interaction, and growth factor receptors all played dominant roles in the epithelial-mesenchymal interface and HF formation. We found that the Wnt, transforming growth factor-beta/bone morphogenetic protein, and Notch family members played vital roles in HF differentiation and maturation. The DEGs we found could be attributed to the generation and development of HF, and thus will be critically important for improving the quantity and quality of fleece production in animals for fibres.