Runs of homozygosity in Swiss goats reveal genetic changes associated with domestication and modern selection

Identification of consensus homozygous regions and their associations with growth and feed efficiency traits in American mink

The recent chromosome-based genome assembly and the newly developed 70K single nucleotide polymorphism (SNP) array for American mink (Neogale vison) facilitate the identification of genetic variants underlying complex traits in this species. The objective of this study was to evaluate the association between consensus runs of homozygosity (ROH) with growth and feed efficiency traits in American mink. A subsample of two mink populations (n = 2,986) were genotyped using the Affymetrix Mink 70K SNP array. The identified ROH segments were included simultaneously, concatenated into consensus regions, and the ROH-based association studies were carried out with linear mixed models considering a genomic relationship matrix for 11 growth and feed efficiency traits implemented in ASReml-R version 4. In total, 298,313 ROH were identified across all individuals, with an average length and coverage of 4.16 Mb and 414.8 Mb, respectively. After merging ROH segments, 196 consensus ROH regions were detected and used for genome-wide ROH-based association analysis. Thirteen consensus ROH regions were significantly (P < 0.01) associated with growth and feed efficiency traits. Several candidate genes within the significant regions are known for their involvement in growth and body size development, including MEF2A, ADAMTS17, POU3F2, and TYRO3. In addition, we found ten consensus ROH regions, defined as ROH islands, with frequencies over 80% of the population. These islands harbored 12 annotated genes, some of which were related to immune system processes such as DTX3L, PARP9, PARP14, CD86, and HCLS1. This is the first study to explore the associations between homozygous regions with growth and feed efficiency traits in American mink. Our findings shed the light on the effects of homozygosity in the mink genome on growth and feed efficiency traits, that can be utilized in developing a sustainable breeding program for mink.

Genetic diversity and selection signatures in Hainan black goats revealed by whole-genome sequencing data

Understanding the genetic characteristics of indigenous goat breeds is crucial for their conservation and breeding efforts. Hainan black goats, as a native breed of south China's tropical island province of Hainan, possess distinctive traits such as black hair, a moderate growth rate, good meat quality, and small body size. However, they exhibit exceptional resilience to rough feeding conditions, possess high-quality meat, and show remarkable resistance to stress and heat. In this study, we resequenced the whole genome of Hainan black goats to study the economic traits and genetic basis of these goats, we leveraged whole-genome sequencing data from 33 Hainan black goats to analyze single nucleotide polymorphism (SNP) density, Runs of homozygosity (ROH), Integrated Haplotype Score (iHS), effective population size (Ne), Nucleotide diversity Analysis (Pi) and selection characteristics. Our findings revealed that Hainan black goats harbor a substantial degree of genetic variation, with a total of 23 608 983 SNPs identified. Analysis of ROHs identified 53 710 segments, predominantly composed of short fragments, with inbreeding events mainly occurring in ancient ancestors, the estimates of inbreeding based on ROH in Hainan black goats typically exhibit moderate values ranging from 0.107 to 0.186. This is primarily attributed to significant declines in the effective population size over recent generations. Moreover, we identified 921 candidate genes within the intersection candidate region of ROH and iHS. Several of these genes are associated with crucial traits such as immunity (PTPRC, HYAL1, HYAL2, HYAL3, CENPE and PKN1), heat tolerance (GNG2, MAPK8, CAPN2, SLC1A1 and LEPR), meat quality (ACOX1, SSTR1, CAMK2B, PPP2CA and PGM1), cashmere production (AKT4, CHRM2, OXTR, AKT3, HMCN1 and CDK19), and stress resistance (TLR2, IFI44, ENPP1, STK3 and NFATC1). The presence of these genes may be attributed to the genetic adaptation of Hainan black goats to local climate conditions. The insights gained from this study provide valuable references and a solid foundation for the preservation, breeding, and utilization of Hainan black goats and their valuable genetic resources.

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.

Differential patterns in runs of homozygosity in two mice lines under divergent selection for environmental variability for birth weight

Runs of homozygosity (ROH) are defined as long continuous homozygous stretches in the genome which are assumed to originate from a common ancestor. It has been demonstrated that divergent selection for variability in mice is possible and that low variability in birth weight is associated with robustness. To analyse ROH patterns and ROH-based genomic inbreeding, two mouse lines that were divergently selected for birth weight variability for 26 generations were used, with: 752 individuals for the high variability line (H-Line), 766 individuals for the low variability line (L-Line) and 74 individuals as a reference population. Individuals were genotyped using the high density Affymetrix Mouse Diversity Genotyping Array. ROH were identified using both the sliding windows (SW) and the consecutive runs (CR) methods. Inbreeding coefficients were calculated based on pedigree (FPED ) information, on ROH identified using the SW method (FROHSW ) and on ROH identified using the CR method (FROHCR ). Differences in genomic inbreeding were not consistent across generations and these parameters did not show clear differences between lines. Correlations between FPED and FROH were high, particularly for FROHSW . Moreover, correlations between FROHSW and FPED were even higher when ROH were identified with no restrictions in the number of heterozygotes per ROH. The comparison of FROH estimates between either of the selected lines were based on significant differences at the chromosome level, mainly in chromosomes 3, 4, 6, 8, 11, 15 and 19. ROH-based inbreeding estimates that were computed using longer homozygous segments had a higher relationship with FPED . Differences in robustness between lines were not attributable to a higher homozygosis in the L-Line, but maybe to the different distribution of ROH at the chromosome level between lines. The analysis identified a set of genomic regions for future research to establish the genomic basis of robustness.

Genome-wide runs of homozygosity signatures in diverse Indian goat breeds

The present study analyzed ROH and consensus ROH regions in 102 animals of eleven diverse Indian goat (Capra hircus) breeds using whole genome sequencing. A total of 51,705 ROH and 21,271 consensus regions were identified. The mean number of ROH per animal was highest in the meat breed, Jharkhand Black (2693) and lowest in the pashmina breed, Changthangi (60). The average length of ROH (ALROH) was maximum in Kanniadu (974.11 Kb) and minimum in Tellicherry (146.98 Kb). Long ROH is typically associated with more recent inbreeding, whereas short ROH is connected to more ancient inbreeding. The overall ROH-based genomic inbreeding (FROH) was highest for Jharkhand Black (0.602) followed by Kanniadu (0.120) and Sangamneri (0.108) among all breeds. FROH of Jharkhand Black was higher than Kanniadu  up to 5 Mb ROH length category. However, in > 20 Mb ROH length category, Kanniadu (0.98)  exhibited significantly higher FROH than Jharkhand Black (0.46). This  implies that Kanniadu had higher levels of recent inbreeding than Jharkhand Black. Despite this, due to the presence of both recent and ancient inbreeding, Jharkhand Black  demonstrated higher overall FROH compared to Kanniadu. ROH patterns revealed dual purpose (meat and dairy) and pashmina breeds as less consanguineous while recent inbreeding was apparent in meat breeds. Analysis of ROH consensus regions identified selection sweeps in key genes governing intramuscular fat deposition, meat tenderisation, lean meat production and carcass weight (CDK4, ALOX15, CASP9, PRDM16, DVL1) in meat breeds; milk fat percentage and mammary gland development (POLD1, NOTCH2, ARHGAP35) in dual purpose (meat and dairy) breeds; while cold adaptation and hair follicle development (APOBEC1, DNAJC3, F2RL1, FGF9) in pashmina breed. MAPK, RAS, BMP and Wnt signaling pathways associated with hair follicle morphogenesis in Changthangi were also identified. PCA analysis based on ROH consensus regions revealed that meat breeds are more diverse than other goat breeds/populations.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-03921-y.

Runs of Homozygosity Detection and Selection Signature Analysis for Local Goat Breeds in Yunnan, China

Runs of Homozygosity (ROH) are continuous homozygous DNA segments in diploid genomes, which have been used to estimate the genetic diversity, inbreeding levels, and genes associated with specific traits in livestock. In this study, we analyzed the resequencing data from 10 local goat breeds in Yunnan province of China and five additional goat populations obtained from a public database. The ROH analysis revealed 21,029 ROH segments across the 15 populations, with an average length of 1.27 Mb, a pattern of ROH, and the assessment of the inbreeding coefficient indicating genetic diversity and varying levels of inbreeding. iHS (integrated haplotype score) was used to analyze high-frequency Single-Nucleotide Polymorphisms (SNPs) in ROH regions, specific genes related to economic traits such as coat color and weight variation. These candidate genes include OCA2 (OCA2 melanosomal transmembrane protein) and MLPH (melanophilin) associated with coat color, EPHA6 (EPH receptor A6) involved in litter size, CDKAL1 (CDK5 regulatory subunit associated protein 1 like 1) and POMC (proopiomelanocortin) linked to weight variation and some putative genes associated with high-altitude adaptability and immune. This study uncovers genetic diversity and inbreeding levels within local goat breeds in Yunnan province, China. The identification of specific genes associated with economic traits and adaptability provides actionable insights for utilization and conservation efforts.

Genome-wide detection of positive and balancing signatures of selection shared by four domesticated rainbow trout populations (Oncorhynchus mykiss)

BACKGROUND

Evolutionary processes leave footprints along the genome over time. Highly homozygous regions may correspond to positive selection of favorable alleles, while maintenance of heterozygous regions may be due to balancing selection phenomena. We analyzed data from 176 fish from four disconnected domestic rainbow trout populations that were genotyped using a high-density Axiom Trout genotyping 665K single nucleotide polymorphism array, including 20 from the US and 156 from three French lines. Using methods based on runs of homozygosity and extended haplotype homozygosity, we detected signatures of selection in these four populations.

RESULTS

Nine genomic regions that included 253 genes were identified as being under positive selection in all four populations Most were located on chromosome 2 but also on chromosomes 12, 15, 16, and 20. In addition, four heterozygous regions that contain 29 genes that are putatively under balancing selection were also shared by the four populations. These were located on chromosomes 10, 13, and 19. Regardless of the homozygous or heterozygous nature of the regions, in each region, we detected several genes that are highly conserved among vertebrates due to their critical roles in cellular and nuclear organization, embryonic development, or immunity. We identified new candidate genes involved in rainbow trout fitness, as well as 17 genes that were previously identified to be under positive selection, 10 of which in other fishes (auts2, atp1b3, zp4, znf135, igf-1α, brd2, col9a2, mrap2, pbx1, and emilin-3).

CONCLUSIONS

Using material from disconnected populations of different origins allowed us to draw a genome-wide map of signatures of positive selection that are shared between these rainbow trout populations, and to identify several regions that are putatively under balancing selection. These results provide a valuable resource for future investigations of the dynamics of genetic diversity and genome evolution during domestication.

Whole-Genome Resequencing Reveals Genetic Diversity and Wool Trait-Related Genes in Liangshan Semi-Fine-Wool Sheep

Understanding the genetic makeup of local sheep breeds is essential for their scientific conservation and sustainable utilization. The Liangshan semi-fine-wool sheep (LSS), a Chinese semi-fine-wool breed renowned for its soft wool, was analyzed using whole-genome sequencing data including 35 LSS, 84 sheep from other domestic breeds, and 20 Asiatic mouflons. We investigated the genetic composition of LSS by conducting analyses of the population structure, runs of homozygosity, genomic inbreeding coefficients, and selection signature. Our findings indicated that LSS shares greater genetic similarity with Border Leicester and Romney sheep than with Tibetan (TIB), Yunnan (YNS), and Chinese Merino sheep. Genomic analysis indicated low to moderate inbreeding coefficients, ranging from 0.014 to 0.154. In identifying selection signals across the LSS genome, we pinpointed 195 candidate regions housing 74 annotated genes (e.g., IRF2BP2, BVES, and ALOX5). We also found the overlaps between the candidate regions and several known quantitative trait loci related to wool traits, such as the wool staple length and wool fiber diameter. A selective sweep region, marked by the highest value of cross-population extended haplotype homozygosity, encompassed IRF2BP2-an influential candidate gene affecting fleece fiber traits. Furthermore, notable differences in genotype frequency at a mutation site (c.1051 + 46T > C, Chr25: 6,784,190 bp) within IRF2BP2 were observed between LSS and TIB and YNS sheep (Fisher's exact test, p < 2.2 × 10-16). Taken together, these findings offer insights crucial for the conservation and breeding enhancement of LSS.

Genome-Wide Association Study of Milk Composition in Karachai Goats

This study is first to perform a genome-wide association study (GWAS) to investigate the milk quality traits in Karachai goats. The objective of the study was to identify candidate genes associated with milk composition traits based on the identification and subsequent analysis of all possible SNPs, both genome-wide (high-confidence) and suggestive (subthreshold significance). To estimate the milk components, 22 traits were determined, including several types of fatty acids. DNA was extracted from ear tissue or blood samples. A total of 167 Karachai goats were genotyped using an Illumina GoatSNP53K BeadChip panel (Illumina Inc., San Diego, CA, USA). Overall, we identified 167 highly significant and subthreshold SNPs associated with the milk components of Karachai goats. A total of 10 SNPs were located within protein-coding genes and 33 SNPs in close proximity to them (±0.2 Mb). The largest number of genome-wide significant SNPs was found on chromosomes 2 and 8 and some of them were associated with several traits. The greatest number of genome-wide significant SNPs was identified for crude protein and lactose (6), and the smallest number-only 1 SNP-for freezing point depression. No SNPs were identified for monounsaturated and polyunsaturated fatty acids. Functional annotation of all 43 SNPs allowed us to identify 66 significant candidate genes on chromosomes 1, 2, 3, 4, 5, 8, 10, 13, 16, 18, 21, 23, 25, 26, and 27. We considered these genes potential DNA markers of the fatty acid composition of Karachai goat milk. Also, we found 12 genes that had a polygenic effect: most of them were simultaneously associated with the dry matter content and fatty acids (METTL, SLC1A 8, PHACTR1, FMO2, ECI1, PGP, ABCA3, AMDHD2). Our results suggest that the genes identified in our study affecting the milk components in Karachai goats differed from those identified in other breeds of dairy goats.

Genomic selection pressure discovery using site-frequency spectrum and reduced local variability statistics in Pakistani Dera-Din-Panah goat

BACKGROUND

Population geneticists have long sought to comprehend various selection traces accumulated in the goat genome due to natural or human driven artificial selection through breeding practices, which led the wild animals to domestication, so understanding evolutionary process may helpful to utilize the full genetic potential of goat genome.

METHODS AND RESULTS

As a step forward to pinpoint the selection signals in Pakistani Dera-Din-Panah (DDP) goat, whole-genome pooled sequencing (n = 12) was performed, and 618,236,192 clean paired-end reads were mapped against ARS1 reference goat assembly. Five different selection signature statistics were applied using four site-frequency spectrum (SFS) methods (Tajima's D ([Formula: see text]), Fay and Wu's H ([Formula: see text]), Zeng's E ([Formula: see text]), [Formula: see text]) and one reduced local variability approach named pooled heterozygosity ([Formula: see text]). The under-selection regions were annotated with significant threshold values of [Formula: see text]≥4.7, [Formula: see text]≥6, [Formula: see text]≥2.5, Pool-HMM ≥ 12, and [Formula: see text]≥5 that resulted in accumulative 364 candidate gene hits. The highest genomic selection signals were observed on Chr. 4, 6, 10, 12, 15, 16, 18, 20, and 27 and harbor ADAMTS6, CWC27, RELN, MYCBP2, FGF14, STIM1, CFAP74, GNB1, CALML6, TMEM52, FAM149A, NADK, MMP23B, OPN3, FH, MFHAS1, KLKB1, RRM1, KMO, SPEF2, F11, KIT, KMO, ERI1, ATP8B4, and RHOG genes. Next, the validation of our captured genomic hits was also performed by more than one applied statistics which harbor meat production, immunity, and reproduction associated genes to strengthen our hypothesis of under-selection traits in this Pakistani goat breed. Furthermore, common candidate genes captured by more than one statistical method were subjected to gene ontology and KEGG pathway analysis to get insights of particular biological processes associated with this goat breed.

CONCLUSION

Current perception of genomic architecture of DDP goat provides a better understanding to improve its genetic potential and other economically important traits of medium to large body size, milk, and fiber production by updating the genomic insight driven breeding strategies to boost the livestock and agriculture-based economy of the country.

Identifying Candidate Genes for Litter Size and Three Morphological Traits in Youzhou Dark Goats Based on Genome-Wide SNP Markers

This study aimed to reveal the potential genetic basis for litter size, coat colour, black middorsal stripe and skin colour by combining genome-wide association analysis (GWAS) and selection signature analysis and ROH detection within the Youzhou dark (YZD) goat population (n = 206) using the Illumina GoatSNP54 BeadChip. In the GWAS, we identified one SNP (snp54094-scaffold824-899720) on chromosome 11 for litter size, two SNPs on chromosome 26 (snp11508-scaffold142-1990450, SORCS3) and chromosome 12 (snp55048-scaffold842-324525, LOC102187779) for coat colour and one SNP on chromosome 18 (snp56013-scaffold873-22716, TCF25) for the black middorsal stripe. In contrast, no SNPs were identified for skin colour. In selection signature analysis, 295 significant iHS genomic regions with a mean |iHS| score > 2.66, containing selection signatures encompassing 232 candidate genes were detected. In particular, 43 GO terms and one KEGG pathway were significantly enriched in the selected genes, which may contribute to the excellent environmental adaptability and characteristic trait formation during the domestication of YZD goats. In ROH detection, we identified 4446 ROH segments and 282 consensus ROH regions, among which nine common genes overlapped with those detected using the iHS method. Some known candidate genes for economic traits such as reproduction (TSHR, ANGPT4, CENPF, PIBF1, DACH1, DIS3, CHST1, COL4A1, PRKD1 and DNMT3B) and development and growth (TNPO2, IFT80, UCP2, UCP3, GHRHR, SIM1, CCM2L, CTNNA3 and CTNNA1) were revealed by iHS and ROH detection. Overall, this study is limited by the small population size, which affects the results of GWAS to a certain extent. Nevertheless, our findings could provide the first overview of the genetic mechanism underlying these important traits and provide novel insights into the future conservation and utilisation of Chinese goat germplasm resources.

Molecular Cytogenetics in Domestic Bovids: A Review

The discovery of the Robertsonian translocation (rob) involving cattle chromosomes 1 and 29 and the demonstration of its deleterious effects on fertility focused the interest of many scientific groups on using chromosome banding techniques to reveal chromosome abnormalities and verify their effects on fertility in domestic animals. At the same time, comparative banding studies among various species of domestic or wild animals were found useful for delineating chromosome evolution among species. The advent of molecular cytogenetics, particularly the use of fluorescence in situ hybridization (FISH), has allowed a deeper investigation of the chromosomes of domestic animals through: (a) the physical mapping of specific DNA sequences on chromosome regions; (b) the use of specific chromosome markers for the identification of the chromosomes or chromosome regions involved in chromosome abnormalities, especially when poor banding patterns are produced; (c) better anchoring of radiation hybrid and genetic maps to specific chromosome regions; (d) better comparisons of related and unrelated species by comparative FISH mapping and/or Zoo-FISH techniques; (e) the study of meiotic segregation, especially by sperm-FISH, in some chromosome abnormalities; (f) better demonstration of conserved or lost DNA sequences in chromosome abnormalities; (g) the use of informatic and genomic reconstructions, in addition to CGH arrays, to predict conserved or lost chromosome regions in related species; and (h) the study of some chromosome abnormalities and genomic stability using PCR applications. This review summarizes the most important applications of molecular cytogenetics in domestic bovids, with an emphasis on FISH mapping applications.

Genetic Diversity and Selection Signatures in Jianchang Black Goats Revealed by Whole-Genome Sequencing Data

Understanding the genetic composition of indigenous goats is essential to promote the scientific conservation and sustainable utilization of these breeds. The Jianchang Black (JC) goat, a Chinese native breed, is solid black and exhibits crude feed tolerance, but is characterized by a low growth rate and small body size. Based on the whole-genome sequencing data for 30 JC, 41 Jintang Black (JT), and 40 Yunshang Black (YS) goats, and 21 Bezoar ibexes, here, we investigated the genetic composition of JC goats by conducting analyses of the population structure, runs of homozygosity (ROH), genomic inbreeding, and selection signature. Our results revealed that JT and YS showed a close genetic relationship with a non-negligible amount of gene flows but were genetically distant from JC, apart from Bezoars. An average of 2039 ROHs were present in the autosomal genome per individual. The ROH-based inbreeding estimates in JC goats generally showed moderate values ranging from 0.134 to 0.264, mainly due to rapid declines in the effective population size during recent generations. The annotated genes (e.g., IL2, IL7, and KIT) overlapping with ROH islands were significantly enriched in immune-related biological processes. Further, we found 61 genes (e.g., STIM1, MYO9A, and KHDRBS2) under positive selection in JC goats via three complementary approaches, which may underly genetic adaptations to local environmental conditions. Our findings provided references for the conservation and sustainable utilization of JC goats.

Whole-Genome Sequencing Reveals the Genomic Characteristics and Selection Signatures of Hainan Black Goat

Goats have become one of the most adaptive and important livestock species distributed in developing countries in recent years. The Hainan Black goat is a native goat breed of the Hainan region that is generally well-liked by the local population and is thus raised in large numbers. However, the genomic diversity and selective signals of the Hainan Black goat have not been clearly elucidated yet. Therefore, in this study, we performed whole-genome resequencing of 16 Hainan Black goats and compared the results with those of 71 goats of 6 other breeds from different geographic regions. Principal component analysis (PCA) and phylogenetic analysis identified seven lineages for all goats. Hainan Black goats showed the most similarity with Leizhou goats and the least similarity with Boer goats. Selective sweep analysis identified candidate genes associated with various functions, including immune resistance to disease (TNFAIP2 (TNF alpha induced protein 2) and EXOC3L4 (exocyst complex component 3 like 4)), melanin biosynthetic process (CDH15 (cadherin 15), ASIP (agouti signaling protein), and PARD3 (par-3 family cell polarity regulator)), and light sensitivity (CNGB3 (cyclic nucleotide gated channel subunit beta 3) and CNBD1 (cyclic nucleotide binding domain containing 1)), underlying strong selection signatures in Hainan Black goats. The melanin biosynthetic process, circadian entrainment, regulation of cyclic adenosine 3,5-monophosphate (cAMP)-mediated signaling, and the Rap-1 signaling pathway were significantly enriched in Hainan Black and Alashan Cashmere goats. This result may be important for understanding each trait. Selection signature analysis revealed candidate single nucleotide polymorphisms (SNPs) and genes correlated with the traits of Hainan Black goats. Collectively, our results provide valuable insights into the genetic basis of specific traits correlated with the Hainan island climate, artificial selection in certain local goat breeds, and the importance of protecting breed resources.