Milk production is the most important economic trait of dairy goats and a key indicator for genetic improvement and breeding. However, milk yield is a complex phenotypic trait, and its genetic mechanisms are still not fully understood. This study focuses on dairy goats and non-dairy goats. By analyzing the population structure of these two groups, we found that there is a significant genetic distance between the populations of dairy goats and non-dairy goats. Using SNP and Indel analyses to identify selection signals, we identified several genes associated with milk production traits, including MPP7, PRPF6, DNAJC5, TPD52L2, HNF4G, LAMA3, FAM13A, and EPHA5. Through longitudinal GWAS of the milk production traits of 298 dairy goats, we discovered additional genes such as TRNAS-GGA-102, TTC39C, LAMA3, ANKRD29, NPC1, C24H18orf8, LOC108633789, RIOK3, TMEM241, CABLES1, LOC108633781, and RBBP8. Transcriptome sequencing of breast tissues at different lactation stages reveals dynamic LAMA3 expression changes. Three non-synonymous mutations in LAMA3 are identified, with the TT genotype at one site correlating significantly with average milk production in dairy goats. Our study discovered new genetic markers for improving dairy goat genetics and provided valuable insights into the genetic mechanisms underlying complex traits.

Lipids are intimately related to the unique flavor and nutritional values of goat milk. MicroRNAs (miRNA) participate in the regulation of various biological functions, including the synthesis and degradation of lipids. Several studies have shown that miR-103 is involved in the regulation of lipid metabolism, however, the molecular mechanism by which miR-103 regulates lipid metabolism in goat mammary gland is poorly understood. In this study, miR-103 was knocked out in goat mammary epithelial cells (GMECs) by CRISPR/Cas9, and the accumulation of lipid droplets, triglycerides, and cholesterol in the cells was suppressed subsequently. Overexpression or knockdown of miR-103-5p and miR-103-3p in GMECs revealed that it was miR-103-5p that promoted lipid accumulation but not miR-103-3p. In addition, Pantothenate Kinase 3 (PANK3), the host gene of miR-103, and Phospholipid Scramblase 4 (PLSCR4) were identified as the target genes of miR-103-5p by dual fluorescein and miRNA pulldown. Furthermore, we identified that cellular lipid levels were negatively regulated by PANK3 and PLSCR4. Lastly, in miR-103 knockout GMECs, the knockdown of PANK and PLSCR4 rescued the lipid accumulation. These findings suggest that miR-103-5p promotes lipid accumulation by targeting PLSCR4 and the host gene PANK3 in GMECs, providing new insights for the regulation of goat milk lipids via miRNAs.

The milk flavor can be attributed to the presence of numerous flavor molecules and precursors. In this study, we employed widely targeted metabolomic analysis techniques to analyze the metabolic profiles of various milk samples obtained from goats, sheep, dairy cows, and buffaloes. A total of 631 metabolites were identified in the milk samples, which were further categorized into 16 distinct classes. Principal component analysis (PCA) suggested that the metabolite profiles of samples from the same species exhibit clustering, while separated patterns of metabolite profiles are observed across goat, sheep, cow, and buffalo species. The differential metabolites between the groups of each species were screened based on fold change and variable importance in projection (VIP) values. Five core differential metabolites were subsequently identified, including 3-(3-hydroxyphenyl)-3-hydroxypropanoic acid, inosine 5′-triphosphate, methylcysteine, N-cinnamylglycine, and small peptide (L-tyrosine–L-aspartate). Through multiple comparisons, we also screened biomarkers of each type of milk. Our metabolomic data showed significant inter-species differences in the composition and concentration of some compounds, such as organic acids, amino acids, sugars, nucleotides, and their derivatives, which may affect the overall flavor properties of the milk sample. These findings provided insights into the molecular basis underlying inter-species variations in milk flavor.

BackgroundBreeding polled goats is a welfare-friendly approach for horn removal in comparison to invasive methods. To gain a comprehensive understanding of the genetic basis underlying polledness in goats, we conducted whole-genome sequencing of 106 Xinong Saanen dairy goats, including 33 horned individuals, 70 polled individuals, and 3 polled intersexuality syndrome (PIS) individuals.MethodsThe present study employed a genome-wide association study (GWAS) and linkage disequilibrium (LD) analysis to precisely map the genetic locus underlying the polled phenotype in goats.ResultsThe analysis conducted in our study revealed a total of 320 genome-wide significant single nucleotide polymorphisms (SNPs) associated with the horned/polled phenotype in goats. These SNPs exhibited two distinct peaks on chromosome 1, spanning from 128,817,052 to 133,005,441 bp and from 150,336,143 to 150,808,639 bp. The present study identified three genome-wide significant SNPs, namely Chr1:129789816, Chr1:129791507, and Chr1:129791577, as potential markers of PIS-affected goats. The results of our LD analysis suggested a potential association between MRPS22 and infertile intersex individuals, as well as a potential association between ERG and the polled trait in goats.ConclusionWe have successfully identified three marker SNPs closely linked to PIS, as well as several candidate genes associated with the polled trait in goats. These results may contribute to the development of SNP chips for early prediction of PIS in goats, thereby facilitating breeding programs aimed at producing fertile herds with polled traits.

Goat milk exhibits a robust and distinctive “goaty” flavor. However, the underlying genetic basis of goaty flavor remains elusive and requires further elucidation at the genomic level. Through comparative genomics analysis, we identified divergent signatures of certain proteins in goat, sheep, and cow. MMUT has undergone a goat-specific mutation in the B12 binding domain. We observed the goat FASN exhibits nonsynonymous mutations in the acyltransferase domain. Structural variations in these key proteins may enhance the capacity for synthesizing goaty flavor compounds in goat. Integrated omics analysis revealed the catabolism of branched-chain amino acids contributed to the goat milk flavor. Furthermore, we uncovered a regulatory mechanism in which the transcription factor ZNF281 suppresses the expression of the ECHDC1 gene may play a pivotal role in the accumulation of flavor substances in goat milk. These findings provide insights into the genetic basis underlying the formation of goaty flavor in goat milk.Statement of significanceBranched-chain fatty acids (BCFAs) play a crucial role in generating the distinctive “goaty” flavor of goat milk. Whether there is an underlying genetic basis associated with goaty flavor is unknown. To begin deciphering mechanisms of goat milk flavor development, we collected transcriptomic data from mammary tissue of goat, sheep, cow, and buffalo at peak lactation for cross-species transcriptome analysis and downloaded nine publicly available genomes for comparative genomic analysis. Our data indicate that the catabolic pathway of branched-chain amino acids (BCAAs) is under positive selection in the goat genome, and most genes involved in this pathway exhibit significantly higher expression levels in goat mammary tissue compared to other species, which contributes to the development of flavor in goat milk. Furthermore, we have elucidated the regulatory mechanism by which the transcription factor ZNF281 suppresses ECHDC1 gene expression, thereby exerting an important influence on the accumulation of flavor compounds in goat milk. These findings provide insights into the genetic mechanisms underlying flavor formation in goat milk and suggest further research to manipulate the flavor of animal products.

Understanding the differences in genetic variation between local Chinese dairy goat breeds and imported breeds can help germplasm innovation and molecular breeding. However, the research is limited in this area. In this study, whole-genome resequencing data from 134 individuals of both local and imported dairy goat breeds were analyzed, and their differences in genomic genetic variation, genetic diversity, and population structure were subsequently identified. We also screened candidate genes associated with important traits of dairy goats such as milk production (STK3, GHR, PRELID3B), reproduction (ATP5E), growth and development (CTSZ, GHR), and immune function (CTSZ, NELFCD). Furthermore, we examined allele frequency distributions for the genes of interest and found significant differences between the two populations. This study provides valuable resources for the study of genetic diversity in dairy goats and lays the foundation for the selective breeding of dairy goats in the future.

MicroRNA-26 (miR-26a and miR-26b) plays a critical role in lipid metabolism, but its endogenous regulatory mechanism in fatty acid metabolism is not clear in goat mammary epithelial cells (GMECs). GMECs with the simultaneous knockout of miR-26a and miR-26b were obtained using the CRISPR/Cas9 system with four sgRNAs. In knockout GMECs, the contents of triglyceride, cholesterol, lipid droplets, and unsaturated fatty acid (UFA) were significantly reduced, and the expression of genes related to fatty acid metabolism was decreased, but the expression level of miR-26 target insulin-induced gene 1 (INSIG1) was significantly increased. Interestingly, the content of UFA in miR-26a and miR-26b simultaneous knockout GMECs was significantly lower than that in wild-type GMECs and miR-26a- and miR-26b-alone knockout cells. After decreasing INSIG1 expression in knockout cells, the contents of triglycerides, cholesterol, lipid droplets, and UFAs were restored, respectively. Our studies demonstrate that the knockout of miR-26a/b suppressed fatty acid desaturation by upregulating the target INSIG1. This provides reference methods and data for studying the functions of miRNA families and using miRNAs to regulate mammary fatty acid synthesis.

Goat milk is rich in fat and protein, thus, has high nutritional values and benefits human health. However, goaty flavour is a major concern that interferes with consumer acceptability of goat milk and the 4-alkyl-branched-chain fatty acids (vBCFAs) are the major substances relevant to the goaty flavour in goat milk. Previous research reported that the acyl-coenzyme A synthetases (ACSs) play a key role in the activation of fatty acids, which is a prerequisite for fatty acids entering anabolic and catabolic processes and highly involved in the regulation of vBCFAs metabolism. Although ACS genes have been identified in humans and mice, they have not been systematically characterized in goats. In this research, we performed genome-wide characterization of the ACS genes in goats, identifying that a total of 25 ACS genes (without ACSM2A) were obtained in the Capra hircus and each ACS protein contained the conserved AMP-binding domain. Phylogenetic analysis showed that out of the 25 genes, 21 belonged to the ACSS, ACSM, ACSL, ACSVL, and ACSBG subfamilies. However, AACS, AASDH, ACSF, and ACSF3 genes were not classified in the common evolutionary branch and belonged to the ACS superfamily. The genes in the same clade had similar conserved structures, motifs and protein domains. The expression analysis showed that the majority of ACS genes were expressed in multi tissues. The comparative analysis of expression patterns in non-lactation and lactation mammary glands of goat, sheep and cow indicated that ACSS2 and ACSF3 genes may participate in the formation mechanisms of goaty flavour in goat milk. In conclusion, current research provides important genomic resources and expression information for ACSs in goats, which will support further research on investigating the formation mechanisms of the goaty flavour in goat milk.