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.

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.

Identification of differentially expressed long non-coding RNAs and messenger RNAs involved with muscle development in Dazu black goats through RNA sequencing.

This study aimed to explore the genetic basis of muscle development in goats. The transcriptome dataset for differentially expressed lncRNAs (DELs) and differentially expressed genes (DEGs) of goat muscle at different developmental stages were obtained using RNA-Seq. A total of 447,806,481 and 587,559,465 clean reads in the longissimus dorsi muscle of Dazu black goats between 75d embryonic stage and 1d after birth were generated through Illumina paired-end sequencing, and their mapping rates were 89.82 and 90.99%, respectively. Moreover, 4517 DEGs and 648 DELs were identified, and 4784 lncRNA-mRNA targeting relationships were predicted. Gene function annotation results showed that 4101 DEGs were significantly enriched to 1098 GO terms, and 2014 DEGs were significantly enriched to 40 KEGG pathways, including many GO terms and pathways related to muscle development, such as cell differentiation and Wnt signaling pathway. Then, 10 DELs and 20 DEGs were randomly selected for RT-qPCR verification, and the agreement rate between the verification and RNA-Seq results was 90%, indicating the high reliability of the RNA-Seq data analysis. In conclusion, this study obtained several mRNAs and lncRNAs related to the muscle development of Dazu black goats and identified several targeted regulatory pairs of lncRNA-mRNA. This study may serve as a reference to understand the genetic basis and molecular mechanism of muscle development in goats.

Genome-Wide Selective Analysis of Boer Goat to Investigate the Dynamic Heredity Evolution under Different Stages.

Boer goats, as kemp in meat-type goats, are selected and bred from African indigenous goats under a long period of artificial selection. Their advantages in multiple economic traits, particularly their plump growth, have attracted worldwide attention. The current study displayed the genome-wide selection signature analyses of South African indigenous goat (AF), African Boer (BH), and Australian Boer (AS) to investigate the hereditary basis of artificial selection in different stages. Four methods (principal component analysis, nucleotide diversity, linkage disequilibrium decay, and neighbor-joining tree) implied the genomic diversity changes with different artificial selection intensities in Boer goats. In addition, the θπ, FST, and XP-CLR methods were used to search for the candidate signatures of positive selection in Boer goats. Consequently, 339 (BH vs. AF) and 295 (AS vs. BH) candidate genes were obtained from SNP data. Especially, 10 genes (e.g., BMPR1B, DNER, ITGAL, and KIT) under selection in both groups were identified. Functional annotation analysis revealed that these genes are potentially responsible for reproduction, metabolism, growth, and development. This study used genome-wide sequencing data to identify inheritance by artificial selection. The results of the current study are valuable for future molecular-assisted breeding and genetic improvement of goats.

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.

Hereditary Basis of Coat Color and Excellent Feed Conversion Rate of Red Angus Cattle by Next-Generation Sequencing Data.

Angus cattle have made remarkable contributions to the livestock industry worldwide as a commercial meat-type breed. Some evidence supported that Angus cattle with different coat colors have different feed-to-meat ratios, and the genetic basis of their coat color is inconclusive. Here, genome-wide association study was performed to investigate the genetic divergence of black and red Angus cattle with 63 public genome sequencing data. General linear model analysis was used to identify genomic regions with potential candidate variant/genes that contribute to coat color and feed conversion rate. Results showed that six single nucleotide polymorphisms (SNPs) and two insertion−deletions, which were annotated in five genes (ZCCHC14, ANKRD11, FANCA, MC1R, and LOC532875 [AFG3-like protein 1]), considerably diverged between black and red Angus cattle. The strongest associated loci, namely, missense mutation CHIR18_14705671 (c.296T > C) and frameshift mutation CHIR18_12999497 (c.310G>-), were located in MC1R. Three consecutive strongly associated SNPs were also identified and located in FANCA, which is widely involved in the Fanconi anemia pathway. Several SNPs of highly associated SNPs was notably enriched in ZCCHC14 and ANKRD11, which are related to myofiber growth and muscle development. This study provides a basis for the use of potential genetic markers to be used in future breeding programs to improve cattle selection in terms of coat color and meat phenotype. This study is also helpful to understand the hereditary basis of different coat colors and meat phenotypes. However, the putative candidate genes or markers identified in this study require further investigation to confirm their phenotypic causality and potential effective genetic relationships.

Genome-Wide Selection Sweep between Wild and Local Pigs from Europe for the Investigation of the Hereditary Characteristics of Domestication in Sus Scrofa.

The phenotypic characteristics of existing domestic pigs (DPs) greatly differ from those of wild boar (WB) populations thousands of years ago. After thousands of years of human domestication, WB and DP have diverged greatly in terms of genetics. Theoretically, worldwide local pigs have independent contributions from their local WBs at the beginning of Sus scrofa domestication. The investigation of the vicissitude of the heredity material between domestic populations and their wild ancestors will help in further understanding the domestication history of domestic animals. In the present study, we performed a genome-wide association scan (GWSA) and phylogeny estimation with a total of 1098 public European Illumina 60K single nucleotide polymorphism data, which included 650 local DPs and 448 WBs. The results revealed that the phylogenetic relationship of WBs corresponds to their geographical distribution and carries large divergence with DPs, and all WB breeds (e.g., HRWB, SBWB, and TIWB) presents a closely linkage with the middle WB (e.g., HRWB, and PLWB). In addition, 64 selected candidate genes (e.g., IDH2, PIP5K1B, SMARCA2, KIF5C, and TJP2) were identified from GWSA. A total of 63 known multiple biological functional pathways were annotated by 22 genes, and ubiquinone and other terpenoid-quinone biosynthesis pathways that belong to the metabolism of cofactors and vitamins were significantly enriched (p < 0.05). The most frequent (28.57%) pathways were classified under metabolism. We confirmed that the middle European WB has made an important genetic contribution to the entire European WB populations. A series of selected genes discovered from this study provides the scientific community with a deeper understanding of the heredity performance of metabolism and emotion and the real purpose behind domestication.

Identification of differentially expressed microRNAs in ovulatory and subordinate follicles in Dazu black goats.

How multiple ovulations happen in prolific goats is still unknown even though studies on ovarian physiology and folliculogenesis have made extensive progress. MicroRNAs (miRNAs) are a class of endogenous regulatory factors and regulate structural gene expression mainly at the post-transcriptional level. In this study, in the follicular phase, six large follicles were collected from prolific Dazu black goat and used to generate RNA libraries for RNA sequencing. Based on the litter size and average number of ovulatory follicles in Dazu black goats, the largest three follicles were sorted as ovulatory follicles, and the remaining as subordinate ones. In total, 418 known miRNAs and 110 novel miRNAs were found, and the expression of six randomly selected miRNAs was validated by quantitative PCR analysis. Nine miRNAs were differently expressed between the ovulatory and subordinate follicles (p < 0.01). Chi-miR-582-5p, novel-130, chi-miR-214-3p, and chi-miR-500-5p were upregulated in the ovulatory group, and chi-miR-383, chi-miR-130b-5p, chi-miR-92a-3p, chi-miR-125b-5p, and novel-9 were downregulated in the same group. Chi-miR-130b-5p and chi-miR-214-3p were predicted to target at LHR (XM_013967581.1), GDF9 (NM_001285708.1), BMP15 (NM_001285588.1), and CYP19A1 (XM_013967046.1). In conclusion, nine miRNAs were differently expressed between ovulatory and subordinate follicles, and chi-miR-130b-5p and chi-miR-214-3p were predicted to regulate the expression of genes involved in gonadotropin hormone signaling and oocyte-derived growth factors. Additional studies are needed to confirm these findings.

Genetic Basis of Follicle Development in Dazu Black Goat by Whole-Transcriptome Sequencing.

The follicle development (FD) is an important factor determining litter size in animals. Recent studies have found that noncoding RNAs (ncRNAs) play an important role in FD. In particular, the role of the regulatory mechanism of competing endogenous RNAs (ceRNAs) that drive FD has attracted increasing attention. Therefore, this study explored the genetic basis of goat FD by obtaining the complete follicular transcriptome of Dazu black goats at different developmental stages. Results revealed that 128 messenger RNAs (mRNAs), 4 long noncoding RNAs (lncRNAs), 49 microRNAs (miRNAs), and 290 circular RNAs (circRNAs) were significantly differentially expressed (DE) between large and small follicles. Moreover, DEmRNAs were enriched in many signaling pathways related to FD, as well as GO terms related to molecular binding and enzyme activity. Based on the analysis of the ceRNA network (CRN), 34 nodes (1 DElncRNAs, 10 DEcircRNAs, 14 DEmiRNAs, and 9 DEmRNAs) and 35 interactions (17 DEcircRNAs-DEmRNAs, 2 DElncRNAs-DEmiRNAs, and 16 DEmRNA-DEmiRNAs) implied that the CRN could be involved in the FD of goats. In conclusion, we described gene regulation by DERNAs and lncRNA/circRNA-miRNA-mRNA CRNs in the FD of goats. This study provided insights into the genetic basis of FD in precise transcriptional regulation.

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.

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.