Genome-wide selection signatures analysis of litter size in Dazu black goats using single-nucleotide polymorphism

Genome-wide selection reveals candidate genes associated with multiple teats in Hu sheep

Increasing the number of teats in sheep helps to improve the survival rate of sheep lambs after birth. In order to analyze the candidate genes related to the formation of multiple teats in Hu sheep, the present study was conducted to investigate the genetic pattern of multiple teats in Hu sheep. In this study, based on genome-wide data from 157 Hu sheep, Fst, xp-EHH, Pi and iHS signaling were performed, and the top 5% signal regions of each analyzed result were annotated based on the Oar_v4.0 for sheep. The results show that a total of 142 SNP loci were selected. We found that PTPRG, TMEM117 and LRP1B genes were closely associated with polypodium formation in Hu sheep, in addition, among the candidate genes related to polypodium we found genes such as TMEM117, SLC25A21 and NCKAP5 related to milk traits. The present study screened out candidate genes for the formation of multiple teats at the genomic level in Hu sheep.

Population structure and selection signal analysis of indigenous sheep from the southern edge of the Taklamakan Desert

Analyzing the genetic diversity and selection characteristics of sheep (Ovis aries) holds significant value in understanding their environmental adaptability, enhancing breeding efficiency, and achieving effective conservation and rational utilization of genetic resources. In this study, we utilized Illumina Ovine SNP 50 K BeadChip data from four indigenous sheep breeds from the southern margin of the Taklamakan Desert (Duolang sheep: n = 36, Hetian sheep: n = 74, Kunlun sheep: n = 27, Qira black sheep: n = 178) and three foreign meat sheep breeds (Poll Dorset sheep: n = 105, Suffolk sheep: n = 153, Texel sheep: n = 150) to investigate the population structure, genetic diversity, and genomic signals of positive selection within the indigenous sheep. According to the Principal component analysis (PCA), the Neighbor-Joining tree (NJ tree), and Admixture, we revealed distinct clustering patterns of these seven sheep breeds based on their geographical distribution. Then used Cross Population Extended Haplotype Homozygosity (XP-EHH), Fixation Index (FST), and Integrated Haplotype Score (iHS), we identified a collective set of 32 overlapping genes under positive selection across four indigenous sheep breeds. These genes are associated with wool follicle development and wool traits, desert environmental adaptability, disease resistance, reproduction, and high-altitude adaptability. This study reveals the population structure and genomic selection characteristics in the extreme desert environments of native sheep breeds from the southern edge of the Taklimakan Desert, providing new insights into the conservation and sustainable use of indigenous sheep genetic resources in extreme environments. Additionally, these findings offer valuable genetic resources for sheep and other mammals to adapt to global climate change.

Structural variant landscapes reveal convergent signatures of evolution in sheep and goats

BACKGROUND

Sheep and goats have undergone domestication and improvement to produce similar phenotypes, which have been greatly impacted by structural variants (SVs). Here, we report a high-quality chromosome-level reference genome of Asiatic mouflon, and implement a comprehensive analysis of SVs in 897 genomes of worldwide wild and domestic populations of sheep and goats to reveal genetic signatures underlying convergent evolution.

RESULTS

We characterize the SV landscapes in terms of genetic diversity, chromosomal distribution and their links with genes, QTLs and transposable elements, and examine their impacts on regulatory elements. We identify several novel SVs and annotate corresponding genes (e.g., BMPR1B, BMPR2, RALYL, COL21A1, and LRP1B) associated with important production traits such as fertility, meat and milk production, and wool/hair fineness. We detect signatures of selection involving the parallel evolution of orthologous SV-associated genes during domestication, local environmental adaptation, and improvement. In particular, we find that fecundity traits experienced convergent selection targeting the gene BMPR1B, with the DEL00067921 deletion explaining ~10.4% of the phenotypic variation observed in goats.

CONCLUSIONS

Our results provide new insights into the convergent evolution of SVs and serve as a rich resource for the future improvement of sheep, goats, and related livestock.

Whole-genome resequencing reveals genomic variation and dynamics in Ethiopian indigenous goats

Ethiopia has about 52 million indigenous goats with marked phenotypic variability, which is the outcome of natural and artificial selection. Here, we obtained whole-genome sequence data of three Ethiopian indigenous goat populations (Arab, Fellata, and Oromo) from northwestern Ethiopia and analyzed their genome-wide genetic diversity, population structure, and signatures of selection. We included genotype data from four other Ethiopian goat populations (Abergelle, Keffa, Gumuz, and Woyto-Guji) and goats from Asia; Europe; and eastern, southern, western, and northern Africa to investigate the genetic predisposition of the three Ethiopian populations and performed comparative genomic analysis. Genetic diversity analysis showed that Fellata goats exhibited the lowest heterozygosity values (Ho = 0.288 ± 0.005 and He = 0.334 ± 0.0001). The highest values were observed in Arab goats (Ho = 0.310 ± 0.010 and He = 0.347 ± 4.35e-05). A higher inbreeding coefficient (FROH = 0.137 ± 0.016) was recorded for Fellata goats than the 0.105 ± 0.030 recorded for Arab and the 0.112 ± 0.034 recorded for Oromo goats. This indicates that the Fellata goat population should be prioritized in future conservation activities. The three goat populations showed the majority (∼63%) of runs of homozygosity in the shorter (100-150 Kb) length category, illustrating ancient inbreeding and/or small founder effects. Population relationship and structure analysis separated the Ethiopian indigenous goats into two distinct genetic clusters lacking phylogeographic structure. Arab, Fellata, Oromo, Abergelle, and Keffa represented one genetic cluster. Gumuz and Woyto-Guji formed a separate cluster and shared a common genetic background with the Kenyan Boran goat. Genome-wide selection signature analysis identified nine strongest regions spanning 163 genes influencing adaptation to arid and semi-arid environments (HOXC12, HOXC13, HOXC4, HOXC6, and HOXC9, MAPK8IP2), immune response (IL18, TYK2, ICAM3, ADGRG1, and ADGRG3), and production and reproduction (RARG and DNMT1). Our results provide insights into a thorough understanding of genetic architecture underlying selection signatures in Ethiopian indigenous goats in a semi-arid tropical environment and deliver valuable information for goat genetic improvement, conservation strategy, genome-wide association study, and marker-assisted breeding.

Genomic Inbreeding and Runs of Homozygosity Analysis of Cashmere Goat

Cashmere goats are valuable genetic resources which are famous worldwide for their high-quality fiber. Runs of homozygosity (ROHs) have been identified as an efficient tool to assess inbreeding level and identify related genes under selection. However, there is limited research on ROHs in cashmere goats. Therefore, we investigated the ROH pattern, assessed genomic inbreeding levels and examined the candidate genes associated with the cashmere trait using whole-genome resequencing data from 123 goats. Herein, the Inner Mongolia cashmere goat presented the lowest inbreeding coefficient of 0.0263. In total, we identified 57,224 ROHs. Seventy-four ROH islands containing 50 genes were detected. Certain identified genes were related to meat, fiber and milk production (FGF1, PTPRM, RERE, GRID2, RARA); fertility (BIRC6, ECE2, CDH23, PAK1); disease or cold resistance and adaptability (PDCD1LG2, SVIL, PRDM16, RFX4, SH3BP2); and body size and growth (TMEM63C, SYN3, SDC1, STRBP, SMG6). 135 consensus ROHs were identified, and we found candidate genes (FGF5, DVL3, NRAS, KIT) were associated with fiber length or color. These findings enhance our comprehension of inbreeding levels in cashmere goats and the genetic foundations of traits influenced by selective breeding. This research contributes significantly to the future breeding, reservation and use of cashmere goats and other goat breeds.

Single-cell transcriptomic atlas of goat ovarian aging

BACKGROUND

The ovaries are one of the first organs that undergo degenerative changes earlier in the aging process, and ovarian aging is shown by a decrease in the number and quality of oocytes. However, little is known about the molecular mechanisms of female age-related fertility decline in different types of ovarian cells during aging, especially in goats. Therefore, the aim of this study was to reveal the mechanisms driving ovarian aging in goats at single-cell resolution.

RESULTS

For the first time, we surveyed the single-cell transcriptomic landscape of over 27,000 ovarian cells from newborn, young and aging goats, and identified nine ovarian cell types with distinct gene-expression signatures. Functional enrichment analysis showed that ovarian cell types were involved in their own unique biological processes, such as Wnt beta-catenin signalling was enriched in germ cells, whereas ovarian steroidogenesis was enriched in granulosa cells (GCs). Further analysis showed that ovarian aging was linked to GCs-specific changes in the antioxidant system, oxidative phosphorylation, and apoptosis. Subsequently, we identified a series of dynamic genes, such as AMH, CRABP2, THBS1 and TIMP1, which determined the fate of GCs. Additionally, FOXO1, SOX4, and HIF1A were identified as significant regulons that instructed the differentiation of GCs in a distinct manner during ovarian aging.

CONCLUSIONS

This study revealed a comprehensive aging-associated transcriptomic atlas characterizing the cell type-specific mechanisms during ovarian aging at the single-cell level and offers new diagnostic biomarkers and potential therapeutic targets for age-related goat ovarian diseases.

Genome-Wide Association Study of the Reproductive Traits of the Dazu Black Goat (Capra hircus) Using Whole-Genome Resequencing

Reproductive traits are the basic economic traits of goats and important indicators in goat breeding. In this study, Dazu black goats (DBGs; n = 150), an important Chinese local goat breed with excellent reproductive performance, were used to screen for important variation loci and genes of reproductive traits. Through genome-wide association studies (GWAS), 18 SNPs were found to be associated with kidding traits (average litter size, average litter size in the first three parity, and average litter size in the first six parity), and 10 SNPs were associated with udder traits (udder depth, teat diameter, teat length, and supernumerary teat). After gene annotation of the associated SNPs and in combination with relevant references, the candidate genes, namely ATP1A1, LRRC4C, SPCS2, XRRA1, CELF4, NTM, TMEM45B, ATE1, and FGFR2, were associated with udder traits, while the ENSCHIG00000017110, SLC9A8, GLRB, GRIA2, GASK1B, and ENSCHIG00000026285 genes were associated with litter size. These SNPs and candidate genes can provide useful biological information for improvement of the reproductive traits of goats.

Uterus proliferative period ceRNA network of Yunshang black goat reveals candidate genes on different kidding number trait

Pregnancy loss that occurs in the uterus is an important and widespread problem in humans and farm animals and is also a key factor affecting the fecundity of livestock. Understanding the differences in the fecundity of goats may be helpful in guiding the breeding of goats with high fecundity. In this study, we performed RNA sequencing and bioinformatics analysis to study the uterus of Yunshang black goats with high and low fecundity in the proliferative period. We identified mRNAs, long non-coding RNAs (lncRNAs), and microRNAs (miRNAs) by analyzing the uterine transcriptomes. The target genes of the identified miRNAs and lncRNAs were predicted, and miRNA-mRNA interaction and competitive endogenous RNA (ceRNA) networks were constructed. By comparisons between low- and high-fecundity groups, we identified 1,674 differentially expressed mRNAs (914 were upregulated, and 760 were downregulated), 288 differentially expressed lncRNAs (149 were upregulated, and 139 were downregulated), and 17 differentially expressed miRNAs (4 were upregulated, and 13 were downregulated). In addition, 49 miRNA-mRNA pairs and 45 miRNA-lncRNA pairs were predicted in the interaction networks. We successfully constructed a ceRNA interaction network with 108 edges that contained 19 miRNAs, 11 mRNAs, and 73 lncRNAs. Five candidate genes (PLEKHA7, FAT2, FN1, SYK, and ITPR2) that were annotated as cell adhesion or calcium membrane channel protein were identified. Our results provide the overall expression profiles of mRNAs, lncRNAs, and miRNAs in the goat uterus during the proliferative period and are a valuable reference for studies into the mechanisms associated with the high fecundity, which may be helpful to guide goat to reduce pregnancy loss.

Single nucleotide polymorphisms in the 3′ UTR of follistatin-like 4 and scavenger receptor class B member 1 are associated with Dazu black goat litter size

The untranslated regions (UTR) of genes play crucial roles in regulating gene expression at the post-transcriptional level such as affecting mRNA stabilization. In this study, 26 single nucleotide polymorphisms (SNPs) and one deletion located in UTR were genotyped from 186 Dazu black goats via SNaPshot, and the correlation between genotype and litter size was analyzed. Results indicated that two SNP loci, SNP_chr17-20182525 and SNP_chr7-65652612, which are located at the 3′UTR of scavenger receptor class B member 1 and follistatin-like 4, are significantly (p<0.05) correlated with the litter size of first parity goats. SNP_chr7-65652612 is also significantly associated with the total litter size of first and second parity offspring (p<0.05). In conclusion, SNP_chr7-65652612 and SNP_chr17-20182525 have correlation with the litter size of Dazu black goat and they are potential genetic markers for litter size breeding.

Sequencing Reveals Population Structure and Selection Signatures for Reproductive Traits in Yunnan Semi-Fine Wool Sheep (Ovis aries)

Yunnan semi-fine wool sheep are among the most important cultivated sheep breeds in China. However, their population structure, genetic characteristics and traits of interest are poorly studied. In this study, we systematically studied the population characteristics and selection signatures of 40 Yunnan semi-fine wool sheep using SNPs obtained from whole-genome resequencing data. A total of 1393 Gb of clean data were acquired. The mapping rate against the reference genome was 91.23% on average (86.01%–92.26%), and the average sequence depth was 9.51X. After filtering, 28,593,198 SNPs and 4,725,259 indels with high quality were obtained. The heterozygosity rate, inbreeding coefficient and effective population size of the sheep were calculated to preliminarily explore their genetic characteristics. The average heterozygosity rate was 0.264, the average inbreeding coefficient was 0.0099, and the effective population size estimated from the heterozygote excess (HE) was 242.9. Based on the Tajima’s D and integrated haplotype score (iHS) approaches, 562 windows and 11,356 core SNPs showed selection signatures in the Yunnan semi-fine wool sheep population. After genome annotation and gene enrichment analysis, we found traces of early domestication in sensory organs, behavioural activity and the nervous system as well as adaptive changes in reproductive and wool traits under selection in this population. Some selected genes related to litter size, including FSHR, BMPR1B and OXT, were identified as being under selection. Specific missense mutations of the FSHR gene that differed from the reference genome were also identified in the population, and we found some SNP variations that may affect litter size. Our findings provide a theoretical basis for the conservation and utilization of Yunnan semi-fine wool sheep. Furthermore, our results reveal some changes common to sheep after domestication and provide a new opportunity to investigate the genetic variation influencing fecundity within a population evolving under artificial selection.

Genome-Wide Association Study of Body Conformation Traits by Whole Genome Sequencing in Dazu Black Goats

Simple Summary

Body conformation traits are economically important in the goat meat industry. Good growth performance in goats, including an accelerated growth rate, can improve carcass weight and meat yield. The identification of genetic variants associated with these traits provides a basis for the genetic improvement of growth performance. In this study, we measured six body conformation traits, including body height, body length, cannon circumference, chest depth, chest width, and heart girth. By a genome-wide association study of a Chinese meat goat breed, 53 significant single nucleotide polymorphisms and 42 candidate genes associated with these traits were detected. These findings improve our understanding of the genetic basis of body conformation traits in goats.

Abstract

Identifying associations between genetic markers and economic traits has practical benefits for the meat goat industry. To better understand the genomic regions and biological pathways contributing to body conformation traits of meat goats, a genome-wide association study was performed using Dazu black goats (DBGs), a Chinese indigenous goat breed. In particular, 150 DBGs were genotyped by whole-genome sequencing, and six body conformation traits, including body height (BH), body length (BL), cannon circumference (CC), chest depth (CD), chest width (CW), and heart girth (HG), were examined. In total, 53 potential SNPs were associated with these body conformation traits. A bioinformatics analysis was performed to evaluate the genes located close to the significant SNPs. Finally, 42 candidate genes (e.g., PSTPIP2, C7orf57, CCL19, FGF9, SGCG, FIGN, and SIPA1L) were identified as components of the genetic architecture underlying body conformation traits. Our results provide useful biological information for the improvement of growth performance and have practical applications for genomic selection in goats.

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.

Combined approaches to reveal genes associated with litter size in Yunshang black goats

Intensive artificial selection has been imposed in Yunshang black goats, the first black specialist mutton goat breed in China, with a breeding object of improving reproductive performance, which has contributed to reshaping of the genome including the characterization of SNP, ROH and haplotype. However, variation in reproductive ability exists in the present population. A WGS was implemented in two subpopulations (polytocous group, PG, and monotocous group, MG) with evident differences of litter size. Following the mapping to reference genome, and SNP calling and pruning, three approaches - GWAS, ROH analysis and detection of signatures of selection - were employed to unveil candidate genes responsible for litter size. Consequently, 12 candidate genes containing OSBPL8 with the minimum P-value were uncovered by GWAS. Differences were observed in the pattern of ROH between two subpopulations that shared similar low inbreeding coefficients. Two ROH hotspots and 12 corresponding genes emerged from ROH pool association analysis. Based on the nS_L statistic, 15 and 61 promising genes were disclosed under selection for MG and PG respectively. Of them, some promising genes participate in ovarian function (PPP2R5C, CDC25A, ESR1, RPS26 and SERPINBs), seasonal reproduction (DIO3, BTG1 and CRYM) and metabolism (OSBPL8, SLC39A5 and SERPINBs). Our study pinpointed some novel promising genes influencing litter size, provided a comprehensive insight into genetic makeup of litter size and might facilitate selective breeding in goats.

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, F_ST = 0.148, V_ST = 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.

Goat CMTM2: mRNA expression profiles of different alternative spliced variants and associations analyses with growth traits

CKLF like MARVEL transmembrane domain containing 2 (CMTM2) plays crucial roles in spermiogenesis, skeletogenous, growth, and development through PI3K/Akt and other pathways. The purpose of this study was to explore the expression profile and variation of different spliced CMTM2 gene in Shaanbei white cashmere goats, as well as to find the relationships between a CMTM2 promoter region 14 bp genetic variant and growth traits in 1366 Shaanbei white cashmere goats. In this study, we identified alternative CMTM2 splicing and detected the effects of the spliced variants on mRNA expression levels in tissues. Meanwhile, an unreported spliced variant of CMTM2 in goat was identified using in CDS cloning and RT-PCR, namely, CMTM2-AS2. Compared with the normal transcript (CMTM2-AS1), the novel variant had the higher expression level in muscle and liver tissues, indicating that it plays an effective role in growth traits. Furthermore, a 14 bp deletion was detected within CMTM2 promoter region, and the different genotypes were significantly associated with growth traits (e.g., body length, circumference of cannon bone) in the large group of 1366 individuals in Shaanbei white cashmere goats. We found that the body length of the individuals with II (n = 571) genotype had better phenotypes than those with DD (n = 118) and ID (n = 650) genotypes. These results have direct guiding significance for goat breeding in the future and provide a new idea for studying the characteristics and functions of CMTM2 gene in goats.

Myostatin (MSTN) Gene Indel Variation and Its Associations with Body Traits in Shaanbei White Cashmere Goat

Myostatin (MSTN) gene, also known as growth differentiation factor 8 (GDF8), is a member of the transforming growth factor-beta super-family and plays a negative role in muscle development. It acts as key points during pre- and post-natal life of amniotes that ultimately determine the overall muscle mass of animals. There are several studies that concentrate on the effect of a 5 bp insertion/deletion (indel) within the 5' untranslated region (5' UTR) of goat MSTN gene in goats. However, almost all sample sizes were below 150 individuals. Only in Boer goats, the sample sizes reached 482. Hence, whether the 5 bp indel was still associated with the growth traits of goats in large sample sizes which were more reliable is not clear. To find an effective and dependable DNA marker for goat rearing, we first enlarged the sample sizes (n = 1074, Shaanbei White Cashmere goat) which would enhance the robustness of the analysis and did the association analyses between the 5 bp indel and growth traits. Results uncovered that the 5 bp indel was significantly related to body height, height at hip cross, and chest width index (p < 0.05). In addition, individuals with DD genotype had a superior growing performance than those with the ID genotype. These findings suggested that the 5 bp indel in MSTN gene are significantly associated with growth traits and the specific genotype might be promising for maker-assisted selection (MAS) of goats.