Multiomic meta-analysis suggests a correlation between steroid hormone-related genes and litter size in goats.

Litter size is a key indicator of production performance in livestock. However, its genetic basis in goats remains poorly understood. In this work, a genome-wide selection sweep analysis (GWSA) on 100 published goat genomes with different litter rates was performed for the first time to identify candidate genes related to kidding rate. This analysis was combined with the public RNA-sequencing data of ovary tissues (follicular phase) from high- and low-yielding goats. A total of 2278 genes were identified by GWSA. Most of these genes were enriched in signaling pathways related to ovarian follicle development and hormone secretion. Moreover, 208 differentially expressed genes between groups were obtained from the ovaries of goats with different litter sizes. These genes were substantially enriched in the cholesterol and steroid synthesis signaling pathways. Meanwhile, the weighted gene co-expression network was used to perform modular analysis of differentially expressed genes. The results showed that seven modules were reconstructed, of which one module showed a very strong correlation with litter size (r = -0.51 and p-value <0.001). There were 51 genes in this module, and 39 hub genes were screened by Pearson's correlation coefficient between core genes > 0.4, correlation coefficient between module members > 0.80 and intra-module connectivity ≥5. Finally, based on the results of GWSA and hub gene Venn analysis, seven key genes (ACSS2, HECW2, KDR, LHCGR, NAMPT, PTGFR and TFPI) were found to be associated with steroid synthesis and follicle growth development. This work contributes to understanding of the genetic basis of goat litter size and provides theoretical support for goat molecular breeding.

Candidate circRNAs related to skeletal muscle development in Dazu black goats.

Circular RNA (CircRNA), as a classical noncoding RNA, has been proven to regulate skeletal muscle development (SMD). However, the molecular genetic basis of circRNA regulation in muscle cells remains unclear. In this study, the expression patterns of circRNAs in the longissimus dorsi muscle at embryonic day 75 and postnatal day 1 in DBGs were investigated to identify the key circRNAs that play an important role in SMD in goats. A total of 140 significantly and differentially expressed circRNAs (DEcircRNAs) were identified among the groups at different developmental stages. Among the 116 host genes (HGs) of DEcircRNAs, 76 were significantly and differentially expressed, which was confirmed by previous RNA_seq data. Furthermore, the expression pattern of 10 DEcircRNAs with RT-qPCR was verified, which showed 80% concordance rate with that of RNA_seq datasets. Moreover, the authenticity of seven randomly selected DEcircRNAs was verified by PCR Sanger sequencing. Based on the functional annotation results, among the 76 significantly and differentially expressed HGs, 74 were enriched in 845 GO terms, whereas 35 were annotated to 85 KEGG pathways. The results of this study could provide a comprehensive understanding of the genetic basis of circRNAs involved in SMD and muscle growth.

Combined Proteomic and Metabolomic Analysis Reveals Comprehensive Regulation of Somatostatin DNA Vaccine in Goats.

Somatostatin (SS) plays crucial regulatory roles in animal growth and reproduction by affecting the synthesis and secretion of growth hormone (GH). However, the mechanism by which SS regulates growth and development in goats is still unclear. In order to investigate the regulatory networks of the hypothalamus and pituitary in goats affected by SS DNA vaccines, in this study, we used a previously established oral attenuated Salmonella typhimurium SS DNA vaccine, X9241 (ptCS/2SS-asd), to treat wethers. We analyzed the protein changes in hypothalamic and pituitary tissues using a TMT-based proteomics approach. Additionally, we examined the metabolic profiles of the serum of control and immunized wethers through untargeted metabolomics using liquid chromatography-mass spectrometry (LC-MS). Key signaling pathways were identified based on differentially expressed metabolites (DEMs) and differentially expressed proteins (DEPs). Furthermore, the effect of critical DEPs on signaling pathways was confirmed through Western blotting (WB) experiments, which elucidated the mechanism of active SS immunization in wethers. A proteomics analysis revealed that the expression of 58 proteins in the hypothalamus and 124 in the pituitary gland was significantly altered following SS vaccine treatment (fold change > 1.2 or < 0.83, p < 0.05). In the hypothalamus, many DEPs were associated with gene ontology (GO) terms related to neuronal signaling. In contrast, most DEPs were associated with metabolic pathways. In the pituitary gland, the DEPs were largely related to immune and nutrient metabolism functions, with significant enrichment in KEGG pathways, particularly those involving the metabolic pathway, sphingolipid signaling, and the cGMP-PKG signaling pathway. A metabolomic analysis further showed that active SS immunization in wethers led to significant alterations in seven serum metabolites. Notably, the sphingolipid signaling pathway, secondary bile acid synthesis, sphingolipid metabolism, and lysine synthesis were significantly disrupted. SS vaccines induced marked changes in hypothalamic-pituitary proteins in wethers, facilitating alterations in their growth processes. This study not only provides insights into the mechanism of the SS gene in regulating GH secretion in wethers but also establishes a basis for hormone immunoregulation technology to enhance livestock production performance.

Comparative analysis of the genetic diversity of the neutral microsatellite loci and second exon of the goat MHC-DQB1 gene.

This study compared and analyzed the genetic diversity and population structure of exon 2 of the DQB1 gene and 13 autosomal neutral microsatellite markers from 14 Chinese goat breeds to explore the potential evolutionary mechanism of the major histocompatibility complex (MHC). A total of 287 haplotypes were constructed from MHC-DQB1 exon 2 from 14 populations, and 82 nucleotide polymorphic sites (SNPs, 31.78%) and 172 heterozygous individuals (79.12%) were identified. The FST values of the microsatellites and MHC-DQB ranged between 0.01831-0.26907 and 0.00892-0.38871, respectively. Furthermore, 14 goat populations showed rich genetic diversity in the microsatellite loci and MHC-DQB1 exon 2. However, the population structure and phylogenetic relationship represented by the two markers were different. Positive selection and Tajima's D test results showed the occurrence of a diversified selection mechanism, which was primarily based on a positive and balancing selection in goat DQB. This study also found that the DQB sequences of bovines exhibited trans-species polymorphism (TSP) among species and families. In brief, this study indicated that positive and balancing selection played a major role in maintaining the genetic diversity of DQB, and TSP of MHC in bovines was common, which enhanced the understanding of the MHC evolution.

Geographical contrasts of Y-chromosomal haplogroups from wild and domestic goats reveal ancient migrations and recent introgressions.

By their paternal transmission, Y-chromosomal haplotypes are sensitive markers of population history and male-mediated introgression. Previous studies identified biallelic single-nucleotide variants in the SRY, ZFY and DDX3Y genes, which in domestic goats identified four major Y-chromosomal haplotypes, Y1A, Y1B, Y2A and Y2B, with a marked geographical partitioning. Here, we extracted goat Y-chromosomal variants from whole-genome sequences of 386 domestic goats (75 breeds) and seven wild goat species, which were generated by the VarGoats goat genome project. Phylogenetic analyses indicated domestic haplogroups corresponding to Y1B, Y2A and Y2B, respectively, whereas Y1A is split into Y1AA and Y1AB. All five haplogroups were detected in 26 ancient DNA samples from southeast Europe or Asia. Haplotypes from present-day bezoars are not shared with domestic goats and are attached to deep nodes of the trees and networks. Haplogroup distributions for 186 domestic breeds indicate ancient paternal population bottlenecks and expansions during migrations into northern Europe, eastern and southern Asia, and Africa south of the Sahara. In addition, sharing of haplogroups indicates male-mediated introgressions, most notably an early gene flow from Asian goats into Madagascar and the crossbreeding that in the 19th century resulted in the popular Boer and Anglo-Nubian breeds. More recent introgressions are those from European goats into the native Korean goat population and from Boer goat into Uganda, Kenya, Tanzania, Malawi and Zimbabwe. This study illustrates the power of the Y-chromosomal variants for reconstructing the history of domestic species with a wide geographical range.

A deletion variant within the FGF5 gene in goats is associated with gene expression levels and cashmere growth.

The FGF5 gene has been associated with the regulation of fibre length in mammals, including cashmere goats. A deletion variant at ~14 kb downstream of the FGF5 gene showed significant divergence between cashmere and non-cashmere goats in previous studies. In this study, we designed specific primers to genotype the deletion variant. The results of gel electrophoresis and Sanger sequencing revealed that a 507-bp deletion mutation is located at 95 454 685-95 455 191 of chromosome 6 in goats. Genotyping data from a large panel of 288 goats showed that the deletion at the FGF5 gene locus appeared to be associated with cashmere length. The deletion variant was close to fixation (frequency 0.97) in cashmere goats. Furthermore, electrophoretic mobility shift assays for evaluating DNA-protein interaction and mRNA expression levels of FGF5 suggested that the deletion variant may serve as a cis-acting element by specifically binding transcription factors to mediate quantitative changes in FGF5 mRNA expression. Our study illustrates how a structural mutation of the FGF5 gene has contributed to the cashmere growth phenotype in domestic goats. The deletion mutation within the FGF5 gene could potentially serve as a molecular marker of cashmere growth in cashmere goat breeding.

Genetic basis investigation of wattle phenotype in goat using genome-wide sequence data.

In domestic goats, wattles often appear in even numbers, mostly on the neck and a few under the ear. Goat wattle is composed of ectopic cartilage tissue covered by skin and was reported as a dominant inheritance. Thirty-eight goats from two Southwest Chinese breeds were studied to elucidate the genetic basis of wattle phenotype in goat. Their genomes were sequenced for wide-genome selective sweep analysis (WGSA) and a genome-wide association study (GWAS). The WGSA results revealed 500 candidate genes identified by fixation index and π ratio and 261 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways enriched with 195 genes and 38 significantly enriched KEGG items. In particular, three chondrogenesis-related pathways (Wnt, Hippo and MAPK signaling pathways) were found. Among the 500 genes, 474 were enriched to 2855 Gene Ontology items, and four (BMP2, BMP4, RARA and MSX1) were annotated in the regulation and development of chondrogenesis. Four chondrogenesis-related genes (GREM1, NEDD4, ATG7 and ITGA1) were identified from 519 single-nucleotide polymorphisms (SNPs) with a GWAS above the threshold. Six and 11 SNPs on chromosome 10 are located on GREM1 and NEDD4 respectively, and the highest numbers of SNPs on chromosomes 20 and 22 are located on ITGA1 and ATG7 respectively. All of these genes are related to cartilage development. This study identified a series of genes related to chondroplasia by GWAS and WGSA and presented the possibility that wattle inheritance may be influenced by multiple genes. This work provides a new theoretical understanding of the hereditary basis of wattle phenotype.

Investigating genetic diversity and population phylogeny of five Chongqing local chicken populations autosomal using microsatellites.

The genetic diversity and population structures of five Chongqing local chicken populations were investigated using by 24 microsatellite markers. Results revealed that the mean number of alleles (NA) ranged from 7.08 (Daninghe chicken, DN) to 8.46 (Nanchuan chicken, NC). The highest observed heterozygosity (HO) and expected heterozygosity (HE) were observed in DN (HO = 0.7252; HE = 0.7409) and the lowest HO and HE were observed in XS (Xiushan native chicken [XS], HO = 0.5910 and HE = 0.6697). The inbreeding coefficient (FIS) within population ranged from 0.022 (DN) to 0.119 (XS). Among the 24 microsatellite markers, four loci (MCW0111, MCW0016, ADL0278, and MCW0104) deviated from the Hardy-Weinberg equilibrium in all the studied populations. The results of population polygenetic analysis based on Nei's genetic distance and STRUCTURE software showed that the clustering of the five populations was incomplete consistent with geographical distribution. Moreover, a large number of gene flows were widespread among different populations, suggesting that genetic material exchanges occurred due to human activities and migration which was also verified by PCoA. In summary, this study preliminarily showed that Chongqing local chicken populations had rich genetic diversity and remarkable genetic divergence, but still high risk in conversion. These findings would be useful to the management of conservation strategies and the utilization of local chicken populations in further.

Genome-wide selection signal analysis of Australian Boer goat reveals artificial selection imprinting on candidate genes related to muscle development.

As one of the best-known commercial goat breeds in the world, Boer goat has undergone long-term artificial selection for nearly 100 years, and its excellent growth rate and meat production performance have attracted considerable worldwide attention. Herein, we used single nucleotide polymorphisms (SNPs) called from the whole-genome sequencing data of 46 Australian Boer goats to detect polymorphisms and identify genomic regions related to muscle development in comparison with those of 81 non-specialized meat goat individuals from Europe, Africa, and Asia. A total of 13 795 202 SNPs were identified, and the whole-genome selective signal screen with a π ratio of nucleotide diversity (πcase /πcontrol ) and pairwise fixation index (FST ) was analyzed. Finally, we identified 1741 candidate selective windows based on the top 5% threshold of both parameters; here, 449 candidate genes were only found in 727 of these regions. A total of 433 genes out of the 449 genes obtained were annotated to 2729 gene ontology terms, of which 51 were directly linked to muscle development (e.g., muscle organ development, muscle cell differentiation) by 30 candidate genes (e.g., JAK2, KCNQ1, PDE5A, PDLIM5, TBX5). In addition, 246 signaling pathways were annotated by 178 genes, and two pathways related to muscle contraction, including vascular smooth muscle contraction (ADCY7, PRKCB, PLA2G4E, ROCK2) and cardiac muscle contraction (CACNA2D3, CASQ2, COX6B1), were identified. The results could improve the current understanding of the genetic effects of artificial selection on the muscle development of goat. More importantly, this study provides valuable candidate genes for future breeding of goats.

SNP screening of the MSTN gene and correlation analysis between genetic polymorphisms and growth traits in Dazu black goat.

The aim of this study was to detect SNPs in myostatin (MSTN) gene of four goat breeds, and analyze the correlation of these markers on body measurement traits in the Dazu black goat breed. In total, twenty polymorphic sites were found in one hundred forty-eight individuals, and all SNPs were distributed in introns 1 and 2, except g. 425 C > T, which was found in the regulatory region. Three SNPs (g. 2732 C > T, g. 2752 G > A and g. 4552 A > C) were polymorphic in all four breeds. None of the tag SNPs (g. 425 C > T, g. 1583 A > G, 2732 C > T, g. 4552 A > C and g. 5167 C > T) were significantly correlated with body measurement traits (p > 0.05) in the Dazu black goat. However, individuals with genotype combination 3 (GtC 3) of tag SNPs had higher birth weight and weaning weight than individuals with the other genotype combinations (p < 0.05). Moreover, the genotype combination 4 (GtC 4) was significantly associated with body length and height at the age of 2 months (p < 0.05), and genotype combination 13 (GtC 13) was significantly correlated with body height at 6 months (p < 0.05). Briefly, the combined tag SNP markers might be useful for goat marker-associated selective breeding.

Immunization with oral KISS1 DNA vaccine inhibits testicular Leydig cell proliferation mainly via the hypothalamic-pituitary-testicular axis and apoptosis-related genes in goats.

This study aimed to analyze the effect and mechanism of immunization of oral KISS1 DNA vaccine on the proliferation of goat testicular Leydig cells. Ten 8-week-old male goats were randomly divided into KISS1 DNA vaccine and control groups for immunization (five goats each group). These goats were sacrificed at 8 weeks after primary immunization, and the tissue samples of hypothalamus, pituitary, and testis and Leydig cell samples were collected for RT-PCR and CCK8 assay. Immunization with the oral KISS1 DNA vaccine effectively inhibited the proliferation of Leydig cells, the expression of hypothalamus KISS1, GPR54, and GnRH mRNA, pituitary GnRHR and LH mRNA, testicular LHR mRNA, and apoptosis-inhibitory gene Bcl-2 mRNA in Leydig cells. By contrast, the immunization enhanced the mRNA expression of apoptosis-promoting gene Bax and Clusterin in Leydig cells. These findings indicate that immunization with the oral KISS1 DNA vaccine can inhibit the proliferation of goat testicular Leydig cells mainly via the hypothalamic-pituitary-testicular axis and apoptosis-related genes.

Comparative and selection sweep analysis of CNV was associated to litter size in Dazu black goats.

This study aims to identify the relative Copy number variation (CNV) associated with the litter size of Dazu black goats based on the unpublished CNV analytical results of our previously published sequencing data, in which the litter-size groups were classified into extreme low- and high-yield groups. Firstly, to compare the existence of valuable CNV in Dazu black goats with different fertility levels with mixed pools. We obtained 4992 and 4888 CNVs from the HY and LY, which overlapping 1461 genes, and classified on the original CNV type. Three genes [LOC108633278, PPP1R12A, and YIPF4] were observed in the intersection between the HY deletion and the LY duplication groups. Secondly, on individuals level, we identified a novel candidate CNV (Chr1_50215501, FST = 0.148, VST = 0.347) from 214 autosomal credible CNVs to be significant with litter size in the Dazu black goat, which located in the CBLB gene. This finding indicates the CBLB gene may affect the litter size of the Dazu black goats through structural variations, and Chr1_50215501 can be an effective genetic marker for marker-assisted selection breeding, and this study was also helps understand the molecular mechanism related to the goat litter size.

Three SNPs within exons of INHA and ACVR2B genes are significantly associated with litter size in Dazu black goats.

The aim of this study was to analyse the association between single-nucleotide polymorphisms within INHA and ACVR2B and litter size in Dazu black goats. In total, twenty-two SNPs were genotyped in 190 individuals by SNaPshot and resequencing. The results showed that three SNPs (SNP_1, SNP_12 and SNP_13 in this study) were detected to have significant additive genetic effect on the recorded goat litter size (p < .05). The SNP_1 (NC_030809.1), a non-synonymous substitution of G for T at chr2-g. 28314990 in the exon 2 of INHA gene (NM_001285606.1), resulted in homozygote 2 (HOM2) contributed 0.25 and heterozygote (HET) contributed 0.12 larger litter than homozygote 1 (HOM1). Meanwhile, SNP_12 (Chr22-g. 11721225 A > T) and SNP_13 (Chr22-g. 11721227 A > C) (NC_030829.1) simultaneously mutated at the first and third position of a triplet AAA (lysine, K) in the exon 4 of ACVR2B gene (XM_018066623.1) had estimated genetic effects of HOM1 (0.00) and HOM2 (0.03) larger than HET (-0.12). In conclusion, one SNPs (chr2-g. 28314990 T > G) within the exon 2 of INHA and two SNPs (Chr22-g. 11721225 A > T and Chr22-g. 11721227 A > C) in the exon 4 of ACVR2B gene were highly recommended as candidate markers of litter size in Dazu black goats. A large-scale association study to assess the impact of these variants on litter size is still necessary.

Selection signatures of litter size in Dazu black goats based on a whole genome sequencing mixed pools strategy.

Litter size is considered to be the most important index for estimating domestic animal productivity. Due to its complexity, the molecular mechanism of litter size has not been elucidated, and it has restricted the use of marker-assisted selection to create high-yield populations in goats. A genome-wide selective sweep analysis was performed with 31 Dazu black goats to identify significant genomic regions and candidate genes related to litter size by a mixed pools strategy. A total of 96 candidate genes were identified, including NR6A1, STK3, IGF2BP2, AR, HMGA2, NPTX1, ANKRD17, DPYD, CLRB, PPP3CA, PLCB1, STK3 and HMGA2, using mixed pool analysis with ZHp and FST. We classified these candidate genes based on the functional classification and annotation of signaling pathways. According to the GO and KEGG analysis results, a total of 43 GO terms and 108 pathways were annotated from these genes. In particular, some novel candidate genes were enriched in reproduction-related pathways, such as the estrogen signaling pathway and oocyte meiosis. These findings provide insight into the influences of coding genes on the fecundity traits of goats.

Current genetic diversity in eight local Chinese sheep populations.

China has numerous local domestic sheep breeds. In this study, the genetic diversity of eight sheep populations was estimated using 17 microsatellites. Knowledge of such diversity provides novel insight into the degree of breed protection needed and the prediction of hybrid advantage. In total, 17 microsatellites were genotyped in 186 individuals from eight populations. The mean number of alleles (± SD) ranged from 3.71 ± 1.36 in Zhaotong sheep to 11.94 ± 3.58 in small-tailed Han sheep. The observed heterozygote frequency (± SD) within a population ranged from 0.482 ± 0.025 in Zhaotong sheep to 0.664 ± 0.023 in Tibetan sheep. In addition, using pairwise difference (FST) analysis, the highest within-population diversity was observed in Tibetan sheep (πX = 12.8098) and small-tailed Han (πX = 12.67873), and the lowest diversity was observed in Zhaotong sheep (πX = 7.90337). The results for genetic divergence between populations indicated that the populations were significantly different (P < 0.05) based on the average number of pairwise differences between populations (πXY) and the corrected average pairwise differences. Both phylogenetic networks and structure analysis showed that these eight populations were separated into three clusters in accordance with their geographical habitat, except Tibetan and Hu sheep. In short, we genotyped eight local Chinese sheep populations using 17 microsatellites, and the results indicated that their current genetic diversity is decreasing and that new conservation strategies are needed. In addition, significant genetic differences between populations could be used in cross breeding.

Screening for signatures of selection of Tianzhu white yak using genome-wide re-sequencing.

The Tianzhu white yak, a domestic yak indigenous to the Qilian Mountains, migrated inland from the Qinghai-Tibet Plateau. Specific ecological and long-term artificial selection influenced the evolution of its pure white coat and physiological characteristics. Therefore, it is not only a natural population that represents a genomic selective region of environmental adaptability but is also an animal model for studying the pigmentation of the yak coat. A total of 24 261 829 variants, including 22 445 252 SNPs, were obtained from 29 yaks by genome-wide re-sequencing. According to the results of a selective sweep analysis of Tianzhu white yak in comparison to Tibetan yaks, nine candidate genes under selection in Tianzhu white yak were identified by combining π, Tajima's D, πA/πB and FST statistics, with threshold standards of 5%. These genes include PDCD1, NUP210, ABCG8, NEU4, LOC102287650, D2HGDH, COL4A1, RTP5 and HDAC11. Five of the nine genes were classified into 12 molecular signaling pathways, and most of these signaling pathways are involved in environmental information processing, organismal systems and metabolism. A majority of these genes has not been implicated in previous studies of yak coat color and high-altitude animals. Our findings are helpful not only for explaining the molecular mechanism of yak coat pigmentation but also for exploring the genetic changes in Tianzhu white yak due to environmental adaptation.

Whole-genome analysis identifying candidate genes of altitude adaptive ecological thresholds in yak populations.

The domestic yak (Bos grunniens) is an iconic symbol of animal husbandry on the Qinghai-Tibet Plateau. Long-term domestication and natural selection have led to a wide distribution of yak, forming many ecological populations to adapt to the local ecological environment. High altitude is closely related to oxygen density, and it is an important environmental ecological factor for biological survival and livestock production. The aim of the present study was to perform a preliminary analysis to identify the candidate genes of altitude distribution adapted ecological thresholds in yak using next-generation sequence technology. A total of 15,762,829 SNPs were obtained from 29 yaks with high- and low-altitude distribution by genome-wide sequencing. According to the results of the selective sweep analysis with FST and ZHp, 21 candidate genes were identified. 14 genes (serine/threonine protein kinase TNNI3K, TEN1, DYM, ITPR1, ZC4H2, KNTC1, ADGRB3, CLYBL, TANGO6, ASCC3, KLHL3, PDE4D, DEPDC1B and AGBL4) were grouped into 32 Gene Ontology terms, and four genes (RPS6KA6, ITPR1, GNAO1 and PDE4D) annotated in 35 pathways, including seven environmental information processing and one environmental adaptation. Therefore, the novel candidate genes found in the current study do not only support new theories about high-altitude adaptation, but also further explain the molecular mechanisms of altitude adaptation threshold in yaks.

WITHDRAWN:Phylogenetic relationship and variation among four strains of Small East African goats of Tanzania based on mitochondrial DNA.

Ahead of Print article withdrawn by publisher.