The dynamic impact approach (DIA) represents an alternative to overrepresentation analysis (ORA) for functional analysis of time-course experiments or those involving multiple treatments. The DIA can be used to estimate the biological impact of the differentially expressed genes (DEGs) associated with particular biological functions, for example, as represented by the Kyoto encyclopedia of genes and genomes (KEGG) annotations. However, the DIA does not take into account the correlated dependence structure of the KEGG pathway hierarchy. We have developed herein a path analysis model (KEGG-PATH) to subdivide the total effect of each KEGG pathway into the direct effect and indirect effect by taking into account not only each KEGG pathway itself, but also the correlation with its related pathways. In addition, this work also attempts to preliminarily estimate the impact direction of each KEGG pathway by a gradient analysis method from principal component analysis (PCA). As a result, the advantage of the KEGG-PATH model is demonstrated through the functional analysis of the bovine mammary transcriptome during lactation.

Recent advances in the study of the CRISPR/Cas9 system have provided a precise and versatile approach for genome editing in various species. However, the applicability and efficiency of this method in large animal models, such as the goat, have not been extensively studied. Here, by co-injection of one-cell stage embryos with Cas9 mRNA and sgRNAs targeting two functional genes (MSTN and FGF5), we successfully produced gene-modified goats with either one or both genes disrupted. The targeting efficiency of MSTN and FGF5 in cultured primary fibroblasts was as high as 60%, while the efficiency of disrupting MSTN and FGF5 in 98 tested animals was 15% and 21% respectively and 10% for double gene modifications. The on- and off-target mutations of the target genes in fibroblasts, as well as in somatic tissues and testis of founder and dead animals, were carefully analyzed. The results showed that simultaneous editing of several sites was achieved in large animals, demonstrating that the CRISPR/Cas9 system has the potential to become a robust and efficient gene engineering tool in farm animals and therefore will be critically important and applicable for breeding.

Goat domestication was critical for agriculture and civilization, but its underlying genetic changes and selection regimes remain unclear. Here, we analyze the genomes of worldwide domestic goats, wild caprid species, and historical remains, providing evidence of an ancient introgression event from a West Caucasian tur-like species to the ancestor of domestic goats. One introgressed locus with a strong signature of selection harbors the MUC6 gene, which encodes a gastrointestinally secreted mucin. Experiments revealed that the nearly fixed introgressed haplotype confers enhanced immune resistance to gastrointestinal pathogens. Another locus with a strong signal of selection may be related to behavior. The selected alleles at these two loci emerged in domestic goats at least 7200 and 8100 years ago, respectively, and increased to high frequencies concurrent with the expansion of the ubiquitous modern mitochondrial haplogroup A. Tracking these archaeologically cryptic evolutionary transformations provides new insights into the mechanisms of animal domestication.

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, ZH_p 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.

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.

Bacterial communities in gastrointestinal tracts (GIT) play an important role in animal health and performance. Despite its importance, little information is available on the establishment of microbial populations in the goat GIT or on changes occurring during early development. Therefore, this study investigated the bacterial community dynamics of the rumen, duodenum, jejunum, ileum, cecum, and colon in 15 goats at five developmental stages (0, 14, 28, 42, and 56 days old) by using 16S rDNA sequencing and quantitative real-time PCR technology. 940 genera were found to belong to 44 phyla distributed along the GIT. As a whole, the microbial richness and diversity showed a clear increasing trend as the kids aged and alpha diversity differed significantly among GIT compartments mainly occurring at middle day ages (14 and 28 days). Principal coordinate analysis indicated that the bacterial community displayed distinct temporal and spatial specificity along the GIT in preweaning goats. As kids aged, the phylum Firmicutes was replaced by Bacteroidetes in rumen, whereas Proteobacteria in the large intestine was displaced by Firmicutes. The phylum Proteobacteria was mainly present in the small intestine in older animals. In the rumen, taxa, such as Bacillus and Lactococcus decreased and Prevotella, Treponema, Ruminococcus, and unclassified Prevotellaceae increased with the age of kids. Furthermore, a lower proportion of taxa, such as Lactobacillus and Bacteroides was observed with higher abundances of both Christensenellaceae_R_7 and Ruminococcus in duodenum and jejunum in older animals. In the large intestine, the microbiota displayed taxonomic dynamics with increases of Ruminococcaceae UCG 005, unclassified Lachnospiraceae, Barnesiella, and Blautia as kids aged. Predicted pathway analysis suggested that genes involved in amino acid metabolism, and translation were abundant in both rumen and duodenum, while genes involved in membrane transport and carbohydrate metabolism were enriched in the large intestine. These results indicate that both the microbial colonization process and potential function exert a temporal-spatial specificity throughout the GIT of goats. This study provides new insight into the temporal dynamics of GIT microbiota development during preweaning and will aid to develop strategies for improving animal health and downstream production.