A 1.1 Mb duplication CNV on chromosome 17 contributes to skeletal muscle development in Boer goats.

The Boer goat is one of the top meat breeds in modern animal husbandry and has attracted widespread attention for its unique growth performance. However, the genetic basis of muscle development in the Boer goat remains obscure. In this study, we identified specific structural variants in the Boer goat based on genome-wide selection signals and analyzed the basis of the molecular heredity of related candidate genes in muscle development. A total of 9 959 autosomal copy number variations (CNVs) were identified through selection signal analysis in 127 goat genomes. Specifically, we confirmed that the highest signal CNV (HSV) was a chromosomal arrangement containing an approximately 1.11 Mb (CHIR17: 60062304-61171840 bp) duplicated fragment inserted in reverse orientation and a 5 362 bp deleted region (CHIR17:60145940-60151302 bp) with overlapping genes (e.g., ARHGAP10, NR3C2, EDNRA, PRMT9, and TMEM184C). The homozygous duplicated HSV genotype (+/+) was found in 96% of Boer goats but was not detected in Eurasian goats and was only detected in 4% of indigenous African goats. The expression network of three candidate genes ( ARHGAP10, NR3C2, and EDNRA) regulating dose transcription was constructed by RNA sequencing. Results indicated that these genes were involved in the proliferation and differentiation of skeletal muscle satellite cells (SMSCs) and their overexpression significantly increased the expression of SAA3. The HSV of the Boer goat contributed to superior skeletal muscle growth via the dose effects of overlapping genes.

Genome-wide selective detection of Mile red-bone goat using next-generation sequencing technology.

The ecotype population of goats (Capra hircus) was created by long-term artificial selection and natural adaptation. Mile red-bone goat is an indigenous breed with visible red bones, and its special bone structure has received extensive attention. This study aimed to identify genetic variants and candidate genes associated with specific bone phenotypes using next-generation sequencing technology (NGS). The results revealed that 31,828,206 single nucleotide polymorphisms (SNPs) were obtained from 72 goats (20 Mile red-bone goats and 52 common goats) by NGS. A total of 100 candidate genes were identified on the basis top 1% window interaction from nucleotide diversity (π), π ratio (π A/π B), and pairwise fixation index (F ST). Exactly 77 known signaling pathways were enriched. Specifically, three coding genes (NMNAT2, LOC102172983, and PNLIP) were annotated in the vitamin metabolism signaling pathways, and NCF2 was annotated to the osteoclast (OC) differentiation pathway. Furthermore, 5862 reliable copy number variations (CNVs) were obtained, and 14 and 24 genes were annotated with the top 1‰ CNV based on F ST (>0.490) and V ST (>0.527), respectively. Several pathways related to bone development and metabolism of exogenous substances in vivo, including calcium signaling pathway, OC differentiation, and glycerophospholipid metabolism, were annotated. Specifically, six genes from 19 candidate CNVs, which were obtained by interaction of the top 1‰ CNVs with F ST and V ST, were annotated to mucin-type O-glycan biosynthesis and metabolic pathways. Briefly, the results implied that pseudopurpurin and specific genetic variants work together to contribute to the red-bone color and specific bone structure of Mile red-bone goat. This study is helpful to understanding the genetic basis of the unique bone phenotype of Mile red-bone goats.

Evolutionary relationship and population structure of domestic Bovidae animals based on MHC-linked and neutral autosomal microsatellite markers.

Major histocompatibility complex (MHC) genes are critical for disease resistance or susceptibility responsible for host-pathogen interactions determined mainly by extensive polymorphisms in the MHC genes. Here, we examined the diversity and phylogenetic pattern of MHC haplotypes reconstructed using three MHC-linked microsatellite markers in 55 populations of five Bovidae species and compared them with those based on neutral autosomal microsatellite markers (NAMs). Three-hundred-and-forty MHC haplotypes were identified in 1453 Bovidae individuals, suggesting significantly higher polymorphism and heterozygosity compared with those based on NAMs. The ambitious boundaries in population differentiation (phylogenetic network, pairwise FST and STRUCTURE analyses) within and between species assessed using the MHC haplotypes were different from those revealed by NAMs associated closely with speciation, geographical distribution, domestication and management histories. In addition, the mean FST was significantly correlated negatively with the number of observed alleles (NA), observed (HO) and expected (HE) heterozygosity and polymorphism information content (PIC) (P < 0.05) in the MHC haplotype dataset while there was no correction of the mean FST estimates (P> 0.05) between the MHC haplotype and NAMs datasets. Analysis of molecular variance (AMOVA) revealed a lower percentage of total variance (PTV) between species/groups based on the MHC-linked microsatellites than NAMs. Therefore, it was inferred that individuals within populations accumulated as many MHC variants as possible to increase their heterozygosity and thus the survival rate of their affiliated populations and species, which eventually reduced population differentiation and thereby complicated their classification and phylogenetic relationship inference. In summary, host-pathogen coevolution and heterozygote advantage, rather than demographic history, act as key driving forces shaping the MHC diversity within the populations and determining the interspecific MHC diversity.

Identification of a Goat Intersexuality-Associated Novel Variant Through Genome-Wide Resequencing and Hi-C.

Background: Polled intersex syndrome (PIS) leads to reproductive disorders in goats and exerts a heavy influence on goat breeding. Since 2001, the core variant of an 11.7 kb deletion at ~129 Mb on chromosome 1 (CHI1) has been widely used as a genetic diagnostic criterion. In 2020, a ~0.48 Mb insertion within the PIS deletion was identified by sequencing in XX intersex goats. However, the suitability of this variation for the diagnosis of intersex goats worldwide and its further molecular genetic mechanism need to be clarified. Results: The whole-genome selective sweep of intersex goats from China was performed with whole-genome next-generation sequencing technology for large sample populations and a case-control study on interbreeds. A series of candidate genes related to the goat intersexuality phenotype were found. We further confirmed that a ~0.48 Mb duplicated fragment (including ERG and KCNJ15) downstream of the ~20 Mb PIS region was reversely inserted into the PIS locus in intersex Chinese goats and was consistent with that in European Saanen and Valais black-necked goats. High-throughput chromosome conformation capture (Hi-C) technology was then used to compare the 3D structures of the PIS variant neighborhood in CHI1 between intersex and non-intersex goats. A newly found structure was validated as an intrachromosomal rearrangement. This inserted duplication changed the original spatial structure of goat CHI1 and caused the appearance of several specific loop structures in the adjacent ~20 kb downstream region of FOXL2. Conclusions: Results suggested that the novel complex PIS variant genome was sufficient as a broad-spectrum clinical diagnostic marker of XX intersexuality in goats from Europe and China. A series of private dense loop structures caused by segment insertion into the PIS deletion might affect the expression of FOXL2 or other neighboring novel candidate genes. However, these structures require further in-depth molecular biological experimental verification. In general, this study provided new insights for future research on the molecular genetic mechanism underlying female-to-male sex reversal in goats.