A genome-wide association study reveals candidate genes for the supernumerary nipple phenotype in sheep (Ovis aries)

Genome-wide association analysis of nine reproduction and morphological traits in three goat breeds from Southern China

OBJECTIVE

This study aimed to investigate the significant single nucleotide polymorphisms (SNPs) and genes associated with nine reproduction and morphological traits in three breed populations of Chinese goats.

METHODS

The genome-wide association of nine reproduction and morphological traits (litter size, nipple number, wattle, skin color, coat color, black dorsal line, beard, beard length, and hind leg hair) were analyzed in three Chinese native goat breeds (n = 336) using an Illumina Goat SNP50 Beadchip.

RESULTS

A total of 17 genome-wide or chromosome-wide significant SNPs associated with one reproduction trait (litter size) and six morphological traits (wattle, coat color, black dorsal line, beard, beard length, and hind leg hair) were identified in three Chinese native goat breeds, and the candidate genes were annotated. The significant SNPs and corresponding putative candidate genes for each trait are as follows: two SNPs located on chromosomes 6 (CSN3) and 24 (TCF4) for litter size trait; two SNPs located on chromosome 9 (KATNA1) and 1 (UBASH3A) for wattle trait; three SNPs located on chromosome 26 (SORCS3), 24 (DYM), and 20 (PDE4D) for coat color trait; two SNPs located on chromosome 18 (TCF25) and 15 (CLMP) for black dorsal line trait; four SNPs located on chromosome 8, 2 (PAX3), 5 (PIK3C2G), and 28 (PLA2G12B and OIT3) for beard trait; one SNP located on chromosome 18 (KCNG4) for beard length trait; three SNPs located on chromosome 17 (GLRB and GRIA2), 28 (PGBD5), and 4 for hind leg hair trait. In contrast, there were no SNPs identified for nipple number and skin color.

CONCLUSION

The significant SNPs or genes identified in this study provided novel insights into the genetic mechanism underlying important reproduction and morphological traits of three local goat breeds in Southern China as well as further potential applications for breeding goats.

Design and characterization of a high-resolution multiple-SNP capture array by target sequencing for sheep

The efficiency of molecular breeding largely depends on inexpensive genotyping arrays. In this study, we aimed to develop an ovine high-resolution multiple-single-nucleotide polymorphism (SNP) capture array, based on genotyping by target sequencing (GBTS) system with capture-in-solution (liquid chip) technology. All the markers were from 40K captured regions, including genes located within selective sweep regions, breed-specific regions, quantitative trait loci (QTL), and the potential functional SNPs on the sheep genome. The results showed that a total of 210K high-quality SNPs were identified in the 40K regions, indicating a high average capture ratio (99.7%) for the target genomic regions. Using genotyped data (n = 317) from liquid chip technology, we further performed genome-wide association studies (GWAS) to detect the genetic loci affecting sheep hair types and teat number. A single significant association signal for hair types was identified on 6.7-7.1 Mb of chromosome 25. The IRF2BP2 gene (chr25: 7,067,974-7,071,785), which is located within this genomic region, has been previously known to be involved in hair/wool traits in sheep. The results further showed a new candidate region around 26.4 Mb of chromosome 13, between the ARHGAP21 and KIAA1217 genes, that was significantly related to teat number in sheep. The haplotype patterns of this region also showed differences in animals with 2, 3, or 4 teats. Advances in using the high-accuracy and low-cost liquid chip are expected to accelerate sheep genomic and breeding studies in the coming years.

Genetics of the phenotypic evolution in sheep: a molecular look at diversity-driving genes

BACKGROUND

After domestication, the evolution of phenotypically-varied sheep breeds has generated rich biodiversity. This wide phenotypic variation arises as a result of hidden genomic changes that range from a single nucleotide to several thousands of nucleotides. Thus, it is of interest and significance to reveal and understand the genomic changes underlying the phenotypic variation of sheep breeds in order to drive selection towards economically important traits.

REVIEW

Various traits contribute to the emergence of variation in sheep phenotypic characteristics, including coat color, horns, tail, wool, ears, udder, vertebrae, among others. The genes that determine most of these phenotypic traits have been investigated, which has generated knowledge regarding the genetic determinism of several agriculturally-relevant traits in sheep. In this review, we discuss the genomic knowledge that has emerged in the past few decades regarding the phenotypic traits in sheep, and our ultimate aim is to encourage its practical application in sheep breeding. In addition, in order to expand the current understanding of the sheep genome, we shed light on research gaps that require further investigation.

CONCLUSIONS

Although significant research efforts have been conducted in the past few decades, several aspects of the sheep genome remain unexplored. For the full utilization of the current knowledge of the sheep genome, a wide practical application is still required in order to boost sheep productive performance and contribute to the generation of improved sheep breeds. The accumulated knowledge on the sheep genome will help advance and strengthen sheep breeding programs to face future challenges in the sector, such as climate change, global human population growth, and the increasing demand for products of animal origin.

Breeding of African sheep reared under low-input/output smallholder production systems for trypanotolerance

Trypanosomiasis is a disease caused by unicellular protozoan parasites. Small ruminants succumb to trypanosomiasis in areas of high tsetse fly challenge, resulting in serious economic loss often to farmers in low-input smallholder systems. At present, trypanosomiasis is treated with trypanocidal drugs, but access to these can be limited, and increasing parasite resistance raises questions about their efficacy. The development of trypanotolerance in small ruminant flocks through targeted breeding strategies is considered a sustainable and economical option for controlling African trypanosomiasis. Recently, quantitative trait loci (QTLs) associated with trypanotolerance traits in sheep have been reported. The results of these studies form the basis for more studies to identify QTLs associated with trypanosomiasis resistance, particularly in African livestock species. For example, signatures of positive selection for trypanotolerance have been identified using genome-wide single-nucleotide polymorphism data. However, there are several challenges in performing genetic analyses using data from low-input smallholder systems, including a lack of recorded pedigree and production records and the need for large sample sizes when flock sizes are often fewer than 50 animals. Breeding strategies to improve trypanotolerance should also preserve existing genetic diversity as well as minimize excessive genetic introgression by trypanosusceptible breeds. This review discusses the possibilities of breeding for trypanosome tolerance/resistance in low-input/low-output small ruminant production systems. Potential challenges are outlined, and potential available genetic resources are described as a foundation for future work.

A study using single-locus and multi-locus genome-wide association study to identify genes associated with teat number in Hu sheep

The multiple teats trait is common in many species of mammals and is considered related to lactation ability in swine. However, in Hu sheep, related gene research is still relatively limited. In this study, a genome‐wide association study was used to identify genetic markers and genes related to the number of teats in the Hu sheep population, a native Chinese sheep breed. A single marker method and several multi‐locus methods were utilized. A total of 61 SNPs were found to be related to the number of teats. Among these, 11 SNPs and one SNP were consistently detected by two and three multi‐locus models respectively. Four SNPs were concordantly identified between the single marker and multi‐locus methods. We also performed quantitative real‐time PCR testing of these identified candidate genes, identifying three genes with significantly different expression. Our study suggested that the LHFP, DPYSL2, and TDP‐43 genes may be related to the number of teats in sheep. The combination of single and multi‐locus GWAS detected additional SNPs not found with only one model. Our results provide new and important insights into the genetic mechanisms of the mammalian multiparous teat phenotype. These findings may be useful for future breeding and understanding the genetics of sheep and other livestock.

Identification of Copy Number Variations and Genetic Diversity in Italian Insular Sheep Breeds

Copy number variants (CNVs) are one of the major contributors to genetic diversity and phenotypic variation in livestock. The aim of this work is to identify CNVs and perform, for the first time, a CNV-based population genetics analysis with five Italian sheep breeds (Barbaresca, Comisana, Pinzirita, Sarda, and Valle del Belìce). We identified 10,207 CNVs with an average length of 1.81 Mb. The breeds showed similar mean numbers of CNVs, ranging from 20 (Sarda) to 27 (Comisana). A total of 365 CNV regions (CNVRs) were determined. The length of the CNVRs varied among breeds from 2.4 Mb to 124.1 Mb. The highest number of shared CNVRs was between Comisana and Pinzirita, and only one CNVR was shared among all breeds. Our results indicated that segregating CNVs expresses a certain degree of diversity across all breeds. Despite the low/moderate genetic differentiation among breeds, the different approaches used to disclose the genetic relationship showed that the five breeds tend to cluster in distinct groups, similar to the previous studies based on single-nucleotide polymorphism markers. Gene enrichment was described for the 37 CNVRs selected, considering the top 10%. Out of 181 total genes, 67 were uncharacterized loci. Gene Ontology analysis showed that several of these genes are involved in lipid metabolism, immune response, and the olfactory pathway. Our results corroborated previous studies and showed that CNVs represent valuable molecular resources for providing useful information for separating the population and could be further used to explore the function and evolutionary aspect of sheep genome.

Variation and frequency of supernumerary teats, litter size, histological features and the fibroblast growth factor 2 (FGF-2) gene expression pattern in goats

Historically, female domestic goats carrying multiple kids are mostly unable to express sufficient nursing ability due to a limited number of functional teats. Therefore, the functional teat is an important component in prolific goat breeding, and plays a key role in the future health of their kids. With this motivation, we wanted to investigate the phenotypic features, litter size, histology of adult female mammary glands, and the gene expression profile of the fibroblast growth factor 2 (FGF-2) gene in goats. To illustrate this, the initial dataset of the current study consists of an electronic questionnaire that includes 697 individuals (548 does and 149 bucks) of five endemic and three exotic goats from 2015 to 2020 in different geographic areas of Iran, from 59 Markhoz (MARG), 50 Azari (AZAR), 73 Busheri (BUSH), 69 Sarbisheh (SARB), 165 Mahabadi (MOHA) indigenous goats and also exotic breeds, including 183 Saanen (SANN), 39 Alpine (ALPN), and 59 Boer (BORE) goats. The results of this study confirmed that MOHA goats (4.16%), BORE (4.43%) and SANN goat breeds (5.75%) have larger litter sizes. Interestingly, the evidence gathering when SNTs occurred showed that both the BUSH and BORE goat breeds had the highest frequency of SNTs. Moreover, under the same physiological and lactation conditions, there was no statistically significant difference in histological features between the three compared does class consist of the two teats, SNTs, and four functional teats. In addition, the results of the gene expression profile significantly highlight the FGF-2 gene pattern in two teat groups compared to other SNT groups (P < 0.01). In summary, this scenario can be used to generate further research and facts on responsible candidate genes, the variations in teat numbers in goats, examining both the incidence of SNT and litter size.

Whole-genome sequence analysis unveils different origins of European and Asiatic mouflon and domestication-related genes in sheep

The domestication and subsequent development of sheep are crucial events in the history of human civilization and the agricultural revolution. However, the impact of interspecific introgression on the genomic regions under domestication and subsequent selection remains unclear. Here, we analyze the whole genomes of domestic sheep and their wild relative species. We found introgression from wild sheep such as the snow sheep and its American relatives (bighorn and thinhorn sheep) into urial, Asiatic and European mouflons. We observed independent events of adaptive introgression from wild sheep into the Asiatic and European mouflons, as well as shared introgressed regions from both snow sheep and argali into Asiatic mouflon before or during the domestication process. We revealed European mouflons might arise through hybridization events between a now extinct sheep in Europe and feral domesticated sheep around 6000-5000 years BP. We also unveiled later introgressions from wild sheep to their sympatric domestic sheep after domestication. Several of the introgression events contain loci with candidate domestication genes (e.g., PAPPA2, NR6A1, SH3GL3, RFX3 and CAMK4), associated with morphological, immune, reproduction or production traits (wool/meat/milk). We also detected introgression events that introduced genes related to nervous response (NEURL1), neurogenesis (PRUNE2), hearing ability (USH2A), and placental viability (PAG11 and PAG3) into domestic sheep and their ancestral wild species from other wild species.

Genome-Wide Detection of Copy Number Variations and Their Association With Distinct Phenotypes in the World's Sheep

Copy number variations (CNVs) are a major source of structural variation in mammalian genomes. Here, we characterized the genome-wide CNV in 2059 sheep from 67 populations all over the world using the Ovine Infinium HD (600K) SNP BeadChip. We tested their associations with distinct phenotypic traits by conducting multiple independent genome-wide tests. In total, we detected 7547 unique CNVs and 18,152 CNV events in 1217 non-redundant CNV regions (CNVRs), covering 245 Mb (∼10%) of the whole sheep genome. We identified seven CNVRs with frequencies correlating to geographical origins and 107 CNVRs overlapping 53 known quantitative trait loci (QTLs). Gene ontology and pathway enrichment analyses of CNV-overlapping genes revealed their common involvement in energy metabolism, endocrine regulation, nervous system development, cell proliferation, immune, and reproduction. For the phenotypic traits, we detected significantly associated (adjusted P < 0.05) CNVRs harboring functional candidate genes, such as SBNO2 for polycerate; PPP1R11 and GABBR1 for tail weight; AKT1 for supernumerary nipple; CSRP1, WNT7B, HMX1, and FGFR3 for ear size; and NOS3 and FILIP1 in Wadi sheep; SNRPD3, KHDRBS2, and SDCCAG3 in Hu sheep; NOS3, BMP1, and SLC19A1 in Icelandic; CDK2 in Finnsheep; MICA in Romanov; and REEP4 in Texel sheep for litter size. These CNVs and associated genes are important markers for molecular breeding of sheep and other livestock species.

Identification of genes associated with litter size combining genomic approaches in Luzhong mutton sheep

Litter size is one of the most important reproductive traits of sheep, which has pronounced effects on the profit of husbandry enterprises and enthusiasm of breeders. Despite the importance of litter size, the underlying genetic mechanisms have not been entirely elucidated. Therefore, based on a high-density SNP chip, genome-wide comparative analysis was performed between two groups with different fecundity to reveal candidate genes linked to litter size via detection of homozygosity and selection signatures in Luzhong mutton sheep. Consequently, nine promising genes were identified from six runs of homozygosity islands, and functionally linked to reproduction (ACTL7A, ACTL7B, and ELP1), embryonic development (KLF5 and PIBF1), and cell cycle (DACH1, BORA, DIS3, and MZT1). A total of 128 genes were observed under selection, of which HECW1 and HTR1E were related to total lambs born, GABRG3, LRP1B, and MACROD2 to teat number, and AGBL1 to reproductive seasonality. Additionally, the presence of inbreeding depression implies the urgency of reasonable mating system to increase litter size in the present herd. These findings provide a comprehensive insight to the genetic makeup of litter size, and also contribute to implementation of marker-assisted selection in sheep.

Whole-genome resequencing of wild and domestic sheep identifies genes associated with morphological and agronomic traits

Understanding the genetic changes underlying phenotypic variation in sheep (Ovis aries) may facilitate our efforts towards further improvement. Here, we report the deep resequencing of 248 sheep including the wild ancestor (O. orientalis), landraces, and improved breeds. We explored the sheep variome and selection signatures. We detected genomic regions harboring genes associated with distinct morphological and agronomic traits, which may be past and potential future targets of domestication, breeding, and selection. Furthermore, we found non-synonymous mutations in a set of plausible candidate genes and significant differences in their allele frequency distributions across breeds. We identified PDGFD as a likely causal gene for fat deposition in the tails of sheep through transcriptome, RT-PCR, qPCR, and Western blot analyses. Our results provide insights into the demographic history of sheep and a valuable genomic resource for future genetic studies and improved genome-assisted breeding of sheep and other domestic animals.

Size of supernumerary teats in sheep correlates with complexity of the anatomy and microenvironment

Supernumerary nipples or teats (polythelia) are congenital accessory structures that may develop at any location along the milk line and have been implicated in the pathogenesis of mastitis. We describe the anatomy and histology of 27 spontaneously occurring supernumerary teats from 16 sheep, delineating two groups of teats - simple and anatomically complex - according to the complexity of the anatomy and microenvironment. Anatomically complex supernumerary teats exhibited significantly increased length and barrel diameter compared with simple supernumerary teats. A teat canal and/or teat cistern was present in anatomically complex teats, with smooth muscle fibres forming a variably well-organised encircling teat sphincter. Complex supernumerary teats also exhibited immune cell infiltrates similar to those of normal teats, including lymphoid follicle-like structures at the folds of the teat cistern-teat canal junction, and macrophages that infiltrated the peri-cisternal glandular tissue. One complex supernumerary teat exhibited teat end hyperkeratosis. These anatomical and histological features allow inference that supernumerary teats may be susceptible to bacterial ingress through the teat canal and we hypothesise that this may be more likely in those teats with less well-organised encircling smooth muscle. The teat cistern of anatomically complex teats may also constitute a focus of milk accumulation and thus a possible nidus for bacterial infection, potentially predisposing to mastitis. We suggest that size of the supernumerary teat, and relationship to the main teats, particularly in the case of 'cluster teats', should be considerations if surgical removal is contemplated.

Genome-wide association study between CNVs and milk production traits in Valle del Belice sheep

Copy number variation (CNV) is a major source of genomic structural variation. The aim of this study was to detect genomic CNV regions (CNVR) in Valle del Belice dairy sheep population and to identify those affecting milk production traits. The GO analysis identified possible candidate genes and pathways related to the selected traits. We identified CNVs in 416 individuals genotyped using the Illumina OvineSNP50 BeadChip array. The CNV association using a correlation-trend test model was examined with the Golden Helix SVS 8.7.0 tool. Significant CNVs were detected when their adjusted p-value was <0.01 after false discovery rate (FDR) correction. We identified 7,208 CNVs, which gave 365 CNVRs after aggregating overlapping CNVs. Thirty-one CNVRs were significantly associated with one or more traits included in the analysis. All CNVRs, except those on OAR19, overlapped with quantitative trait loci (QTL), even if they were not directly related to the traits of interest. A total of 222 genes were annotated within the significantly associated CNVRs, most of which played important roles in biological processes related to milk production and health-related traits. Identification of the genes in the CNVRs associated with the studied traits will provide the basis for further investigation of their role in the metabolic pathways related to milk production and health traits.

Genome-wide association study reveals novel genes for the ear size in sheep (Ovis aries)

Variations in ear size can be observed in livestock such as sheep; however, the genetic basis of variable ear size in sheep is still poorly understood. To investigate causative genes associated with ear size in sheep, a genome-wide association study was performed in 115 adult Duolang sheep with different-sized floppy ears using the Ovine Infinium HD BeadChip. We found 38 significant SNPs at the genome-wide or chromosome-wise 5% significance level after Bonferroni correction. The most significant association (P = 1.61 × 10^-6 ) was found at SNP rs402740419, located in the DCC gene, which plays a critical role in ear development. Also, we observed two additional significant SNPs, rs407891215 in PTPRD and rs407769095 in SOX5, both of which are functionally associated with ear developmental processes. Our results are useful for future sheep breeding and provide insights into the genetic basis of ear size development in sheep and other livestock.

Genome-Wide Association Analyses Highlight the Potential for Different Genetic Mechanisms for Litter Size Among Sheep Breeds

Reproduction is an important trait in sheep breeding as well as in other livestock. However, despite its importance the genetic mechanisms of litter size in domestic sheep (Ovis aries) are still poorly understood. To explore genetic mechanisms underlying the variation in litter size, we conducted multiple independent genome-wide association studies in five sheep breeds of high prolificacy (Wadi, Hu, Icelandic, Finnsheep, and Romanov) and one low prolificacy (Texel) using the Ovine Infinium HD BeadChip, respectively. We identified different sets of candidate genes associated with litter size in different breeds: BMPR1B, FBN1, and MMP2 in Wadi; GRIA2, SMAD1, and CTNNB1 in Hu; NCOA1 in Icelandic; INHBB, NF1, FLT1, PTGS2, and PLCB3 in Finnsheep; ESR2 in Romanov and ESR1, GHR, ETS1, MMP15, FLI1, and SPP1 in Texel. Further annotation of genes and bioinformatics analyses revealed that different biological pathways could be involved in the variation in litter size of females: hormone secretion (FSH and LH) in Wadi and Hu, placenta and embryonic lethality in Icelandic, folliculogenesis and LH signaling in Finnsheep, ovulation and preovulatory follicle maturation in Romanov, and estrogen and follicular growth in Texel. Taken together, our results provide new insights into the genetic mechanisms underlying the prolificacy trait in sheep and other mammals, suggesting targets for selection where the aim is to increase prolificacy in breeding projects.