Genome-Wide Runs of Homozygosity, Effective Population Size, and Detection of Positive Selection Signatures in Six Chinese Goat Breeds

Whole-genome resequencing reveals new mutations in candidate genes for Beichuan-white goat prolificacya

In this study, we evaluated the copy number variation in the genomes of two groups of Beichuan-white goat populations with large differences in litter size by FST method, and identified 1739 genes and 485 missense mutations in the genes subject to positive selection. Through functional enrichment, ITGAV, LRP4, CDH23, TPRN, RYR2 and CELSR1 genes, involved in embryonic morphogenesis, were essential for litter size trait, which received intensive attention. In addition, some mutation sites of these genes have been proposed (ITGAV: c.38C > T; TPRN: c.133A > T, c.1192A > G, c.1250A > C; CELSR1: c.7640T > C), whose allele frequencies were significantly changed in the high fecundity goat group. Besides, we found that new mutations at these sites altered the hydrophilicity and 3D structure of the protein. Candidate genes related to litter size in this study and their missense mutation sites were identified. These candidate genes are helpful to understand the genetic mechanism of fecundity in Beichuan white goat, and have important significance for future goat breeding.

Whole-genome resequencing reveals diversity and selective signals in the Wuxue goat

Animal genetic resources are crucial for ensuring global food security. However, in recent years, a noticeable decline in the genetic diversity of livestock has occurred worldwide. This decline is pronounced in developing countries, where the management of these resources is insufficient. In the current study, we performed whole genome sequencing for 20 Wuxue (WX) and five Guizhou White (GW) goats. Additionally, we utilized the published genomes of 131 samples representing five different goat breeds from various regions in China. We investigated and compared the genetic diversity and selection signatures of WX goats. Whole genome sequencing analysis of the WX and GW populations yielded 120 425 063 SNPs, which resided primarily in intergenic and intron regions. Population genetic structure revealed that WX exhibited genetic resemblance to GW, Chengdu Brown, and Jintang Black and significant differentiation from the other goat breeds. In addition, three methods (nucleotide diversity, linkage disequilibrium decay, and runs of homozygosity) showed moderate genetic diversity in WX goats. We used nucleotide diversity and composite likelihood ratio methods to identify within-breed signatures of positive selection in WX goats. A total of 369 genes were identified using both detection methods, including genes related to reproduction (GRID2, ZNF276, TCF25, and SPIRE2), growth (HMGA2 and GJA3), and immunity (IRF3 and SRSF3). Overall, this study explored the adaptability of WX goats, shedding light on their genetic richness and potential to thrive in challenges posed by climatic changes and diseases. Further investigations are warranted to harness these insights to enhance more efficient and sustainable goat breeding initiatives.

Genomic inbreeding analysis reveals resilience and genetic diversity in Indian yak populations

The yak (Bos grunniens), renowned for its adaptability to extreme cold and hypoxic conditions, stands as a remarkable domestic animal crucial for sustaining livelihoods in harsh climates. We conducted a comprehensive analysis of the whole genome sequence data from three distinct Indian yak populations: Arunachali yak (n = 10), Himachali yak (n = 10), and Ladakhi yak (n = 10). The genomic data for Indian yaks were meticulously generated by our laboratory and compared with their Chinese counterpart, the Jinchuan yak (n = 8), for a more nuanced understanding. Our investigation revealed a total of 37,437 runs of homozygosity (ROH) segments in 34 animals representing four distinct yak populations. The Jinchuan yak population exhibited the highest proportion, constituting 80.8 % of total ROHs, predominantly as small segments (<0.1 Mb), accounting for 63 % of the overall ROHs. Further analysis uncovered a significantly higher degree of inbreeding in Chinese yaks compared to their Indian counterparts. The Indian yak populations, in contrast, demonstrated relatively lower and consistent levels of inbreeding. Moreover, we identified ROH hotspots that covered at least 60 % of individuals in our study, indicating their pivotal role in environmental adaptation. A total of five hotspot regions were detected, housing genes such as ENSBGRG00000015023 (WNT2), YIPF4, SPAST, TLN2, and DSG4. These genes are associated with traits including hair follicle initiation, nutrient stress response, microtubule assembly, development of cardiac muscle, hair follicle, and coat color. This observation strongly suggests that there is substantial selection acting on these genes, emphasizing their important role in environmental adaptation among yak populations.

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.

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.

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.

The coefficients of inbreeding revealed by ROH study among inbred individuals belonging to each type of the first cousin marriage: A preliminary report from North India

BACKGROUND

Genome-wide runs of homozygosity (ROH) are appropriate to estimate genomic inbreeding, determine population history, unravel the genetic architecture of complex traits and disorders.

OBJECTIVE

The study sought to investigate and compare the actual proportion of homozygosity or autozygosity in the genomes of progeny of four subtypes of first cousin mating in humans, using both pedigree and genomic measures for autosomes and sex chromosomes.

METHODS

For this purpose, Illumina Global Screening Array-24 v1.0 BeadChip followed by cyto-ROH analysis through Illumina Genome Studio was used to characterise the homozygosity in five participants from North Indian state (Uttar Pradesh). PLINK v.1.9 software was used to estimate the genomic inbreeding coefficients viz. ROH-based inbreeding estimate (F_ROH) and homozygous loci-based inbreeding estimate (F_HOM).

RESULTS

A total of 133 ROH segments were detected with maximum number and genomic coverage in Matrilateral Parallel (MP) type and minimum in outbred individual. ROH pattern revealed that MP type has a higher degree of homozygosity than other subtypes. The comparison of F_ROH, F_HOM, and pedigree-based inbreeding estimate (F_PED) showed some difference in theoretical and realised proportion of homozygosity for sex-chromosomal loci but not for autosome for each type of consanguinity.

CONCLUSIONS

This is the very first study to compare and estimate the pattern of homozygosity among the kindreds of first cousin unions. However, a greater number of individuals from each type of marriage is required for statistical inference of no difference between theoretical and realized homozygosity among different degrees of inbreeding prevalent in humans worldwide.

Screening Indels from the whole genome to identify the candidates and their association with economic traits in several goat breeds

In the present study, the re-sequencing of our previous whole-genome sequencing (WGS) for selected individuals of Dazu-black goat (DBG) and Inner-Mongolia Cashmere goat (IMCG) breeds was used to detect and compare the differentiation in Indels depending on the reference genome of goat. Then, three selected candidate Indels rs668795676, rs657996810, and rs669452874 of the three genes SUFU, SYCP2L and GLIPR1L1, respectively, have been chosen, based on the results of prior GWAS across the genome, and examined for their association with body weight and dimensions (body height, body length, heart girth, chest width, cannon circumference, and chest depth) by kompetitive allele specific PCR assay for 342 goats from the three studied goat breeds (DBG, n = 203; ♂99, ♀104), IMCG (n = 65; 15♂, 50♀), and Hechuan white goat (HWG, n = 74; 34♂, 40♀) breeds. The analysis of 192.747 Gb WGS revealed an average 334,151 Indels in the whole genome of DBG and IMCG breeds. Chromosome 1 had a maximum number of mutations (Indels) of 58,497 and 55,527 for IMCG and DBG, respectively, while chromosome 25 had the least number of mutations of 15,680 and 16,103 for IMCG and DBG, respectively. The majority of Indels were either Ins or Del of short fragments of 1-5 bp, which covered 79.06 and 71.78% of the total number of Indels mutations in IMCG and DBG, respectively. Comparing the differences of Indels between the studied goat breeds revealed 100 and 110 unique Indels for IMCG and DBG, respectively. The Indels loci in the intron region were unique for both studied goat breeds which were related to 30 and 38 candidate genes in IMCG and DBG, respectively, including SUFU, SYCP2L, and GLIPR1L1 genes. Concerning rs669452874 locus, body height and body length of Del/Del genotype in DBG were significantly higher (P < 0.05) than that of Ins/Del genotype, while body height and body length of Del/Del genotype in IMCG were significantly higher (P < 0.01) than those of Ins/Ins genotype, whereas body height and body length and heart girth of Del/Del genotype in HWG were significantly higher (P < 0.01) than those of the Ins/Del and Ins/Ins genotypes. Thus, Del/Del genotype of rs669452874 locus can be used as a candidate molecular marker related to the body dimensions in the studied goat breeds.

Discovering novel clues of natural selection on four worldwide goat breeds

In goat breeds, the domestication followed by artificial selection for economically important traits have shaped genetic variation within populations, leading to the fixation of specific alleles for specific traits. This led to the formation and evolution of many different breeds specialised and raised for a particular purpose. However, and despite the intensity of artificial selection, natural selection continues acting, possibly leaving a more diluted contribution over time, whose traces may be more difficult to capture. In order to explore selection footprints as response of environmental adaptation, we analysed a total of 993 goats from four transboundary goats breeds (Angora, Boer, Nubian and Saanen) genotyped with the SNP chip 50 K using outlier detection, runs of homozygosity and haplotype-based detection methods. Our results showed that all methods identified footprints on chromosome 6 (from 30 to 49 Mb) for two specific populations of Nubian goats sampled in Egypt. In Angora and Saanen breeds, we detected two selective sweeps using HapFLK, on chromosome 21 (from 52 to 55 Mb) and chromosome 25 (from 1 to 5 Mb) respectively. The analysis of runs of homozygosity showed some hotspots in all breeds. The overall investigation of the selected regions detected combining the different approaches and the gene ontology exploration revealed both novel and well-known loci related to adaptation, especially for heat stress. Our findings can help to better understand the balance between the two selective pressures in commercial goat breeds providing new insights on the molecular mechanisms of adaptation.

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.

Runs of homozygosity reveal past bottlenecks and contemporary inbreeding across diverging populations of an island-colonizing bird

Genomes retain evidence of the demographic history and evolutionary forces that have shaped populations and drive speciation. Across island systems, contemporary patterns of genetic diversity reflect population demography, including colonization events, bottlenecks, gene flow and genetic drift. Here, we investigate genome-wide diversity and the distribution of runs of homozygosity (ROH) using whole-genome resequencing of individuals (>22× coverage) from six populations across three archipelagos of Berthelot's pipit (Anthus berthelotii)-a passerine that has recently undergone island speciation. We show the most dramatic reduction in diversity occurs between the mainland sister species (the tawny pipit) and Berthelot's pipit and is lowest in the populations that have experienced sequential bottlenecks (i.e., the Madeiran and Selvagens populations). Pairwise sequential Markovian coalescent (PSMC) analyses estimated that Berthelot's pipit diverged from its sister species ~2 million years ago, with the Madeiran archipelago founded 50,000 years ago, and the Selvagens colonized 8000 years ago. We identify many long ROH (>1 Mb) in these most recently colonized populations. Population expansion within the last 100 years may have eroded long ROH in the Madeiran archipelago, resulting in a prevalence of short ROH (<1 Mb). However, the extensive long and short ROH detected in the Selvagens suggest strong recent inbreeding and bottleneck effects, with as much as 38% of the autosomes consisting of ROH >250 kb. These findings highlight the importance of demographic history, as well as selection and genetic drift, in shaping contemporary patterns of genomic diversity across diverging populations.

Genomic Regions Associated with Wool, Growth and Reproduction Traits in Uruguayan Merino Sheep

The aim of this study was to identify genomic regions and genes associated with the fiber diameter (FD), clean fleece weight (CFW), live weight (LW), body condition score (BCS), pregnancy rate (PR) and lambing potential (LP) of Uruguayan Merino sheep. Phenotypic records of approximately 2000 mixed-age ewes were obtained from a Merino nucleus flock. Genome-wide association studies were performed utilizing single-step Bayesian analysis. For wool traits, a total of 35 genomic windows surpassed the significance threshold (PVE ≥ 0.25%). The proportion of the total additive genetic variance explained by those windows was 4.85 and 9.06% for FD and CFW, respectively. There were 42 windows significantly associated with LWM, which collectively explained 43.2% of the additive genetic variance. For BCS, 22 relevant windows accounted for more than 40% of the additive genetic variance, whereas for the reproduction traits, 53 genomic windows (24 and 29 for PR and LP, respectively) reached the suggestive threshold of 0.25% of the PVE. Within the top 10 windows for each trait, we identified several genes showing potential associations with the wool (e.g., IGF-1, TGFB2R, PRKCA), live weight (e.g., CAST, LAP3, MED28, HERC6), body condition score (e.g., CDH10, TMC2, SIRPA, CPXM1) or reproduction traits (e.g., ADCY1, LEPR, GHR, LPAR2) of the mixed-age ewes.

Genomic insights into the conservation and population genetics of two Chinese native goat breeds

Chinese goats are an important group of goats worldwide. However, there are few studies on the conservation priority, genetic relationship, and potential gene flow between Chinese and global goat breeds. Here, we genotyped 239 goats from conservation populations of the Chinese Guangfeng and Ganxi breeds using the GoatSNP50 BeadChip. The conservation priority, population structure, selection signatures and introgression of these goats were analyzed in the context of 36 global goat breeds. First, we showed that Guangfeng and Ganxi goats had the largest effective population sizes across the global breeds 13 generations ago. Nevertheless, Ganxi goats have recently experienced a high degree of inbreeding, resulting in their conservation priority based on total gene and allelic diversities being lower than that of most other Chinese breeds (including Guangfeng goats). Population structure and admixture analyses showed that an average of 18% of Guangfeng genomic components were introgressed from Boer goats approximately 18-yr ago. Next, we reconstructed the subfamily structure of the core populations of Guangfeng and Ganxi goats, and proposed reasonable conservation strategies for inbreeding management. Moreover, a list of candidate genes under selection for fertility, immunity, growth, and meat quality were detected in Guangfeng and Ganxi goats. Finally, we identified some genes related to body development and reproduction, which were introgressed from Boer goats and may be beneficial for improving performance and productivity of Guangfeng goats. In conclusion, this study not only provides new insights into the conservation and utilization of Guangfeng and Ganxi goats but also enriches our understanding of artificial introgression from exotic goats into Chinese local goats.

Challenges in quantifying genome erosion for conservation

Massive defaunation and high extinction rates have become characteristic of the Anthropocene. Genetic effects of population decline can lead populations into an extinction vortex, where declining populations show lower genetic fitness, in turn leading to lower populations still. The lower genetic fitness in a declining population due to a shrinking gene pool is known as genetic erosion. Three different types of genetic erosion are highlighted in this review: overall homozygosity, genetic load and runs of homozygosity (ROH), which are indicative of inbreeding. The ability to quantify genetic erosion could be a very helpful tool for conservationists, as it can provide them with an objective, quantifiable measure to use in the assessment of species at risk of extinction. The link between conservation status and genetic erosion should become more apparent. Currently, no clear correlation can be observed between the current conservation status and genetic erosion. However, the high quantities of genetic erosion in wild populations, especially in those species dealing with habitat fragmentation and habitat decline, may be early signs of deteriorating populations. Whole genome sequencing data is the way forward to quantify genetic erosion. Extra screening steps for genetic load and hybridization can be included, since they could potentially have great impact on population fitness. This way, the information yielded from genetic sequence data can provide conservationists with an objective genetic method in the assessment of species at risk of extinction. However, the great complexity of genome erosion quantification asks for consensus and bridging science and its applications, which remains challenging.

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.

Integrated multi-omic data reveal the potential molecular mechanisms of the nutrition and flavor in Liancheng white duck meat

The Liancheng white (LW) duck is one of the most valued Chinese indigenous poultry breeds. Its meat is rich in nutrients and has distinct flavors, but the molecular mechanisms behind them are unknown. To address this issue, we measured and compared multi-omic data (genome, transcriptome, and metabolome) of breast meat from LW ducks and the Mianyang Shelduck (MS) ducks. We found that the LW duck has distinct breed-specific genetic features, including numerous mutant genes with differential expressions associated with amino acid metabolism and transport activities. The metabolome driven by genetic materials was also seen to differ between the two breeds. For example, several amino acids that are beneficial for human health, such as L-Arginine, L-Ornithine, and L-lysine, were found in considerably higher concentrations in LW muscle than in MS duck muscle (p < 0.05). SLC7A6, a mutant gene, was substantially upregulated in the LW group (p < 0.05), which may lead to excessive L-arginine and L-ornithine accumulation in LW duck meat through transport regulation. Further, guanosine monophosphate (GMP), an umami-tasting molecule, was considerably higher in LW muscle (p < 0.05), while L-Aspartic acid was significantly abundant in MS duck meat (p < 0.05), showing that the LW duck has a different umami formation. Overall, this study contributed to our understanding of the molecular mechanisms driving the enriched nutrients and distinct umami of LW duck meat, which will provide a useful reference for duck breeding.

Trajectory of livestock genomics in South Asia: A comprehensive review

Livestock plays a central role in sustaining human livelihood in South Asia. There are numerous and distinct livestock species in South Asian countries. Several of them have experienced genetic development in recent years due to the application of genomic technologies and effective breeding programs. This review discusses genomic studies on cattle, buffalo, sheep, goat, pig, horse, camel, yak, mithun, and poultry. The frontiers covered in this review are genetic diversity, admixture studies, selection signature research, QTL discovery, genome-wide association studies (GWAS), and genomic selection. The review concludes with recommendations for South Asian livestock systems to increasingly leverage genomic technologies, based on the lessons learned from the numerous case studies. This paper aims to present a comprehensive analysis of the dichotomy in the South Asian livestock sector and argues that a realistic approach to genomics in livestock can ensure long-term genetic advancements.

Conservation status and historical relatedness of South African communal indigenous goat populations using a genome-wide single-nucleotide polymorphism marker

Indigenous goats form the majority of populations in smallholder, low input, low output production systems and are considered an important genetic resource due to their adaptability to different production environments and support of communal farming. Effective population size (N_e), inbreeding levels, and the runs of homozygosity (ROHs) are effective tools for exploring the genetic diversity and understanding the demographic history in efforts to support breeding strategies to use and conserve genetic resources. Across populations, the current N_e of Gauteng was the lowest at 371 animals, while the historical N_e across populations suggests that the ancestor N_e has decreased by 53.86%, 44.58%, 42.16%, and 41.16% in Free State (FS), North West (NW), Limpopo (LP), and Gauteng (GP), respectively, over the last 971 generations. Genomic inbreeding levels related to ancient kinship (F_ROH > 5 Mb) were highest in FS (0.08 ± 0.09) and lowest in the Eastern Cape (EC) (0.02 ± 0.02). A total of 871 ROH island regions which include important environmental adaptation and hermo-tolerance genes such as IL10RB, IL23A, FGF9, IGF1, EGR1, MTOR, and MAPK3 were identified (occurring in over 20% of the samples) in FS (n = 37), GP (n = 42), and NW (n = 2) populations only. The mean length of ROH across populations was 7.76 Mb and ranged from 1.61 Mb in KwaZulu-Natal (KZN) to 98.05 Mb (GP and NW). The distribution of ROH according to their size showed that the majority (n = 1949) of the detected ROH were > 5 Mb in length compared to the other categories. Assuming two hypothetical ancestral populations, the populations from KZN and LP are revealed, supporting PC 1. The genomes of KZN and LP share a common origin but have substantial admixture from the EC and NW populations. The findings revealed that the occurrence of high N_e and autozygosity varied largely across breeds in communal indigenous goat populations at recent and ancient events when a genome-wide single-nucleotide polymorphism (SNP) marker was used. The use of Illumina goat SNP50K BeadChip shows that there was a migration route of communal indigenous goat populations from the northern part (LP) of South Africa to the eastern areas of the KZN that confirmed their historical relatedness and coincides with the migration periods of the Bantu nation.

Simultaneous testing of rule- and model-based approaches for runs of homozygosity detection opens up a window into genomic footprints of selection in pigs

BACKGROUND

Past selection events left footprints in the genome of domestic animals, which can be traced back by stretches of homozygous genotypes, designated as runs of homozygosity (ROHs). The analysis of common ROH regions within groups or populations displaying potential signatures of selection requires high-quality SNP data as well as carefully adjusted ROH-defining parameters. In this study, we used a simultaneous testing of rule- and model-based approaches to perform strategic ROH calling in genomic data from different pig populations to detect genomic regions under selection for specific phenotypes.

RESULTS

Our ROH analysis using a rule-based approach offered by PLINK, as well as a model-based approach run by RZooRoH demonstrated a high efficiency of both methods. It underlined the importance of providing a high-quality SNP set as input as well as adjusting parameters based on dataset and population for ROH calling. Particularly, ROHs ≤ 20 kb were called in a high frequency by both tools, but to some extent covered different gene sets in subsequent analysis of ROH regions common for investigated pig groups. Phenotype associated ROH analysis resulted in regions under potential selection characterizing heritage pig breeds, known to harbour a long-established breeding history. In particular, the selection focus on fitness-related traits was underlined by various ROHs harbouring disease resistance or tolerance-associated genes. Moreover, we identified potential selection signatures associated with ear morphology, which confirmed known candidate genes as well as uncovered a missense mutation in the ABCA6 gene potentially supporting ear cartilage formation.

CONCLUSIONS

The results of this study highlight the strengths and unique features of rule- and model-based approaches as well as demonstrate their potential for ROH analysis in animal populations. We provide a workflow for ROH detection, evaluating the major steps from filtering for high-quality SNP sets to intersecting ROH regions. Formula-based estimations defining ROHs for rule-based method show its limits, particularly for efficient detection of smaller ROHs. Moreover, we emphasize the role of ROH detection for the identification of potential footprints of selection in pigs, displaying their breed-specific characteristics or favourable phenotypes.

Genome-Wide Estimates of Runs of Homozygosity, Heterozygosity, and Genetic Load in Two Chinese Indigenous Goat Breeds

Runs of homozygosity (ROH) and heterozygosity (ROHet) are windows into population demographic history and adaptive evolution. Numerous studies have shown that deleterious mutations are enriched in the ROH of humans, pigs, cattle, and chickens. However, the relationship of deleterious variants to ROH and the pattern of ROHet in goats have been largely understudied. Here, 240 Guangfeng and Ganxi goats from Jiangxi Province, China, were genotyped using the Illumina GoatSNP50 BeadChip and genome-wide ROH, ROHet, and genetic load analyses were performed in the context of 32 global goat breeds. The classes with the highest percentage of ROH and ROHet were 0.5–2 Mb and 0.5–1 Mb, respectively. The results of inbreeding coefficients (based on SNP and ROH) and ROHet measurements showed that Guangfeng goats had higher genetic variability than most Chinese goats, while Ganxi goats had a high degree of inbreeding, even exceeding that of commercial goat breeds. Next, the predicted damaging homozygotes were more enriched in long ROHs, especially in Guangfeng goats. Therefore, we suggest that information on damaging alleles should also be incorporated into the design of breeding and conservation programs. A list of genes related to fecundity, growth, and environmental adaptation were identified in the ROH hotspots of two Jiangxi goats. A sense-related ROH hotspot (chromosome 12: 50.55–50.81 Mb) was shared across global goat breeds and may have undergone selection prior to goat domestication. Furthermore, an identical ROHet hotspot (chromosome 1: 132.21–132.54 Mb) containing two genes associated with embryonic development (STAG1 and PCCB) was detected in domestic goat breeds worldwide. Tajima’s D and BetaScan2 statistics indicated that this region may be caused by long-term balancing selection. These findings not only provide guidance for the design of conservation strategies for Jiangxi goat breeds but also enrich our understanding of the adaptive evolution of goats.

Selection Signature Analyses Revealed Genes Associated With Adaptation, Production, and Reproduction in Selected Goat Breeds in Kenya

Artificial and natural selection in livestock is expected to leave unique footprints on their genomes. Goat breeds in Kenya have evolved for survival, breeding, and production in various harsh ecological areas, and their genomes are likely to have acquired unique alleles for adaptation to such diverse production environments and other traits of economic importance. To investigate signals of selection for some selected goat breeds in Kenya, Alpine (n = 29), Galla (n = 12), Saanen (n = 24), and Toggenburg (n = 31) were considered. A total of 53,347 single-nucleotide polymorphisms (SNPs) generated using the Illumina GoatSNP50 BeadChip were analyzed. After quality control, 47,663 autosomal single-nucleotide polymorphisms remained for downstream analyses. Several complementary approaches were applied for the following analyses: integrated Haplotype Score (iHS), cross-population-extended haplotype homozygosity (XP-EHH), hapFLK, and FLK. A total of 404 top genomic regions were identified across all the four breeds, based on the four complementary analyses. Out of the 16 identified putative selection signature regions by the intersection of multiple-selective signal analyses, most of the putative regions were found to overlap significantly with the iHS and XP-EHH analyses on chromosomes 3, 4, 10, 15, 22, and 26. These regions were enriched with some genes involved in pathways associated directly or indirectly with environmental adaptation regulating immune responses (e.g., HYAL1 and HYAL3), milk production (e.g., LEPR and PDE4B), and adaptability (e.g., MST1 and PCK). The results revealed few intersect between breeds in genomic selection signature regions. In general, this did not present the typical classic selection signatures as predicted due to the complex nature of the traits. The results support that some various selection pressures (e.g., environmental challenges, artificial selection, and genome admixture challenges) have molded the genome of goat breeds in Kenya. Therefore, the research provides new knowledge on the conservation and utilization of these goat genetic resources in Kenya. In-depth research is needed to detect precise genes connected with adaptation and production in goat breeds in Kenya.

The individual and global impact of copy-number variants on complex human traits

The impact of copy-number variations (CNVs) on complex human traits remains understudied. We called CNVs in 331,522 UK Biobank participants and performed genome-wide association studies (GWASs) between the copy number of CNV-proxy probes and 57 continuous traits, revealing 131 signals spanning 47 phenotypes. Our analysis recapitulated well-known associations (e.g., 1q21 and height), revealed the pleiotropy of recurrent CNVs (e.g., 26 and 16 traits for 16p11.2-BP4-BP5 and 22q11.21, respectively), and suggested gene functionalities (e.g., MARF1 in female reproduction). Forty-eight CNV signals (38%) overlapped with single-nucleotide polymorphism (SNP)-GWASs signals for the same trait. For instance, deletion of PDZK1, which encodes a urate transporter scaffold protein, decreased serum urate levels, while deletion of RHD, which encodes the Rhesus blood group D antigen, associated with hematological traits. Other signals overlapped Mendelian disorder regions, suggesting variable expressivity and broad impact of these loci, as illustrated by signals mapping to Rotor syndrome (SLCO1B1/3), renal cysts and diabetes syndrome (HNF1B), or Charcot-Marie-Tooth (PMP22) loci. Total CNV burden negatively impacted 35 traits, leading to increased adiposity, liver/kidney damage, and decreased intelligence and physical capacity. Thirty traits remained burden associated after correcting for CNV-GWAS signals, pointing to a polygenic CNV architecture. The burden negatively correlated with socio-economic indicators, parental lifespan, and age (survivorship proxy), suggesting a contribution to decreased longevity. Together, our results showcase how studying CNVs can expand biological insights, emphasizing the critical role of this mutational class in shaping human traits and arguing in favor of a continuum between Mendelian and complex diseases.

Runs of Homozygosity Uncover Potential Functional-Altering Mutation Associated With Body Weight and Length in Two Duroc Pig Lines

Runs of homozygosity (ROH) are widely used to investigate genetic diversity, demographic history, and positive selection signatures of livestock. Commercial breeds provide excellent materials to reveal the landscape of ROH shaped during the intense selection process. Here, we used the GeneSeek Porcine 50K single-nucleotide polymorphism (SNP) Chip data of 3,770 American Duroc (AD) and 2,096 Canadian Duroc (CD) pigs to analyze the genome-wide ROH. First, we showed that AD had a moderate genetic differentiation with CD pigs, and AD had more abundant genetic diversity and significantly lower level of inbreeding than CD pigs. In addition, sows had larger levels of homozygosity than boars in AD pigs. These differences may be caused by differences in the selective intensity. Next, ROH hotspots revealed that many candidate genes are putatively under selection for growth, sperm, and muscle development in two lines. Population-specific ROHs inferred that AD pigs may have a special selection for female reproduction, while CD pigs may have a special selection for immunity. Moreover, in the overlapping ROH hotspots of two Duroc populations, we observed a missense mutation (rs81216249) located in the growth and fat deposition-related supergene (ARSB-DMGDH-BHMT) region. The derived allele of this variant originated from European pigs and was nearly fixed in Duroc pigs. Further selective sweep and association analyses indicated that this supergene was subjected to strong selection and probably contributed to the improvement of body weight and length in Duroc pigs. These findings will enhance our understanding of ROH patterns in different Duroc lines and provide promising trait-related genes and a functional-altering marker that can be used for genetic improvement of pigs.

Genetic Diversity, Population Structure and Phylogeny of Indigenous Goats of Mongolia Revealed by SNP Genotyping

Simple Summary

We performed SNP genotyping of indigenous goats of Mongolia to explore their demographic history in the global context and to estimate their genetic risks. Recently, these risks have become a subject of concern due to recent climatic disasters and uncontrolled massive breeding. Various clustering methods demonstrated close genetic relations among Mongolian, Russian, Chinese, and West Asian breeds. Mongolian goats themselves exhibited low to moderate estimates of genetic differentiation. We identified genetic features highlighting the distinct origin and breeding history of Mongolian goat breeds, as well as traces of artificial selection and geographic isolation. However, none of them met formal criteria to be considered as endangered.

Abstract

Mongolian goats are of great interest for studying ancient migration routes and domestication, and also represent a good model of adaptability to harsh environments. Recent climatic disasters and uncontrolled massive breeding endangered the valuable genetic resources of Mongolian goats and raised the question of their conservation status. Meanwhile, Mongolian goats have never been studied on genomic scale. We used Illumina Goat SNP50 to estimate genetic risks in five Mongolian goat breeds (Buural, Ulgii Red, Gobi GS, Erchim, Dorgon) and explored phylogenic relationships among these populations and in the context of other breeds. Various clustering methods showed that Mongolian goats grouped with other Asian breeds and were especially close to some neighboring Russian and Chinese breeds. The Buural breed showed the lowest estimates of inbreeding and exhibited the shortest genetic distances within the other Mongolian breeds, especially, to Ulgii Red and Gobi GS. These three breeds formed a single core group, being weakly differentiated from each other. Among them, Gobi GS displayed obvious signs of inbreeding probably resulted from artificial selection pressure. Dorgon and especially Erchim goats stand apart from the other Mongolian breeds according to various types of analyses, and bear unique features pointing to different breeding histories or distinct origins of these breeds. All populations showed strong decline in effective population size. However, none of them met formal criteria to be considered as endangered breeds. The SNP data obtained in this study improved the knowledge of Mongolian goat breeds and could be used in future management decisions in order to preserve their genetic diversity.

Genome-wide selection signatures detection in Shanghai Holstein cattle population identified genes related to adaption, health and reproduction traits

Background

Over several decades, a wide range of natural and artificial selection events in response to subtropical environments, intensive pasture and intensive feedlot systems have greatly changed the customary behaviour, appearance, and important economic traits of Shanghai Holstein cattle. In particular, the longevity of the Shanghai Holstein cattle population is generally short, approximately the 2nd to 3rd lactation. In this study, two complementary approaches, integrated haplotype score (iHS) and runs of homozygosity (ROH), were applied for the detection of selection signatures within the genome using genotyping by genome-reduced sequence data from 1092 cows.

Results

In total, 101 significant iHS genomic regions containing selection signatures encompassing a total of 256 candidate genes were detected. There were 27 significant |iHS| genomic regions with a mean |iHS| score > 2. The average number of ROH per individual was 42.15 ± 25.47, with an average size of 2.95 Mb. The length of 78 % of the detected ROH was within the range of 1–2 MB and 2–4 MB, and 99 % were shorter than 8 Mb. A total of 168 genes were detected in 18 ROH islands (top 1 %) across 16 autosomes, in which each SNP showed a percentage of occurrence > 30 %. There were 160 and 167 genes associated with the 52 candidate regions within health-related QTL intervals and 59 candidate regions within reproduction-related QTL intervals, respectively. Annotation of the regions harbouring clustered |iHS| signals and candidate regions for ROH revealed a panel of interesting candidate genes associated with adaptation and economic traits, such as IL22RA1, CALHM3, ITGA9, NDUFB3, RGS3, SOD2, SNRPA1, ST3GAL4, ALAD, EXOSC10, and MASP2. In a further step, a total of 1472 SNPs in 256 genes were matched with 352 cis-eQTLs in 21 tissues and 27 trans-eQTLs in 6 tissues. For SNPs located in candidate regions for ROH, a total of 108 cis-eQTLs in 13 tissues and 4 trans-eQTLs were found for 1092 SNPs. Eighty-one eGenes were significantly expressed in at least one tissue relevant to a trait (P value < 0.05) and matched the 256 genes detected by iHS. For the 168 significant genes detected by ROH, 47 gene-tissue pairs were significantly associated with at least one of the 37 traits.

Conclusions

We provide a comprehensive overview of selection signatures in Shanghai Holstein cattle genomes by combining iHS and ROH. Our study provides a list of genes associated with immunity, reproduction and adaptation. For functional annotation, the cGTEx resource was used to interpret SNP-trait associations. The results may facilitate the identification of genes relevant to important economic traits and can help us better understand the biological processes and mechanisms affected by strong ongoing natural or artificial selection in livestock populations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08042-x.

Genetic structure analysis in several populations of cattle using SNP genotypes

Parameters such as effective population size (Ne), runs of homozygosity (ROH), and inbreeding based on ROH (FROH) can give new insight into the level of genetic diversity for the population under selection. This research aimed to measure the extent of linkage disequilibrium (LD), effective population size (Ne), Haplotype Block Structure, and runs of homozygosity (ROHs) in several populations of cattle using SNP genotypes. In this study, that the average r² decreased with the increasing distance of SNP pairs. A general decrease in Ne can be seen for all four populations, indicating a loss of genetic diversity. The Iranian Holstein had the lowest level of genomic inbreeding at an ROH of 1, 5, 10 Mb, while the French Holstein had the highest. The maximum number of ROH is seen at a distance of less than 1 Mb, and the lowest number of ROH is seen at a distance of 10 Mb. The number of ROH decreases with increasing distance due to the increased recombination rate. This is a concern as an increase in inbreeding leads to a reduction in the effective population size, which was also evident in the study populations.

Genome-Wide Association Study Identifies Candidate Genes Associated with Feet and Leg Conformation Traits in Chinese Holstein Cattle

Simple Summary

Feet and leg problems are among the major reasons for dairy cows leaving the herd, as well as having direct association with production and reproduction efficiency, health (e.g., claw disorders and lameness) and welfare. Hence, understanding the genetic architecture underlying feet and conformation traits in dairy cattle offers new opportunities toward the genetic improvement and long-term selection. Through a genome-wide association study on Chinese Holstein cattle, we identified several candidate genes associated with feet and leg conformation traits. These results could provide useful information about the molecular breeding basis of feet and leg traits, thus improving the longevity and productivity of dairy cattle.

Abstract

Feet and leg conformation traits are considered one of the most important economical traits in dairy cattle and have a great impact on the profitability of milk production. Therefore, identifying the single nucleotide polymorphisms (SNPs), genes and pathways analysis associated with these traits might contribute to the genomic selection and long-term plan selection for dairy cattle. We conducted genome-wide association studies (GWASs) using the fixed and random model circulating probability unification (FarmCPU) method to identify SNPs associated with bone quality, heel depth, rear leg side view and rear leg rear view of Chinese Holstein cows. Phenotypic measurements were collected from 1000 individuals of Chinese Holstein cattle and the GeneSeek Genomic Profiler Bovine 100 K SNP chip was utilized for individual genotyping. After quality control, 984 individual cows and 84,906 SNPs remained for GWAS work; as a result, we identified 20 significant SNPs after Bonferroni correction. Several candidate genes were identified within distances of 200 kb upstream or downstream to the significant SNPs, including ADIPOR2, INPP4A, DNMT3A, ALDH1A2, PCDH7, XKR4 and CADPS. Further bioinformatics analyses showed 34 gene ontology terms and two signaling pathways were significantly enriched (p ≤ 0.05). Many terms and pathways are related to biological quality, metabolism and development processes; these identified SNPs and genes could provide useful information about the genetic architecture of feet and leg traits, thus improving the longevity and productivity of Chinese Holstein dairy cattle.

SNP-Based Genotyping Provides Insight Into the West Asian Origin of Russian Local Goats

Specific local environmental and sociocultural conditions have led to the creation of various goat populations in Russia. National goat diversity includes breeds that have been selected for down and mohair production traits as well as versatile local breeds for which pastoralism is the main management system. Effective preservation and breeding programs for local goat breeds are missing due to the lack of DNA-based data. In this work, we analyzed the genetic diversity and population structure of Russian local goats, including Altai Mountain, Altai White Downy, Dagestan Downy, Dagestan Local, Karachaev, Orenburg, and Soviet Mohair goats, which were genotyped with the Illumina Goat SNP50 BeadChip. In addition, we addressed genetic relationships between local and global goat populations obtained from the AdaptMap project. Russian goats showed a high level of genetic diversity. Although a decrease in historical effective population sizes was revealed, the recent effective population sizes estimated for three generations ago were larger than 100 in all studied populations. The mean runs of homozygosity (ROH) lengths ranged from 79.42 to 183.94 Mb, and the average ROH number varied from 18 to 41. Short ROH segments (<2 Mb) were predominant in all breeds, while the longest ROH class (>16 Mb) was the least frequent. Principal component analysis, Neighbor-Net graph, and Admixture clustering revealed several patterns in Russian local goats. First, a separation of the Karachaev breed from other populations was observed. Moreover, genetic connections between the Orenburg and Altai Mountain breeds were suggested and the Dagestan breeds were found to be admixed with the Soviet Mohair breed. Neighbor-Net analysis and clustering of local and global breeds demonstrated the close genetic relations between Russian local and Turkish breeds that probably resulted from past admixture events through postdomestication routes. Our findings contribute to the understanding of the genetic relationships of goats originating in West Asia and Eurasia and may be used to design breeding programs for local goats to ensure their effective conservation and proper management.

Microbiome analysis reveals the alterations in gut microbiota in different intestinal segments of Yimeng black goats

Mounting evidence revealed the importance of gut microbiota in host metabolism, immunity and physiology, and health. Yimeng black goats (YBGs) mainly distributed in Shandong province of China, displayed a complicated intestinal microecosystem, but studies of its gut microbiota are still insufficient to report. Therefore, this study was performed with an objective to characterize the intestinal microbial community structure and diversity in the small intestine (duodenum, jejunum and ileum) and cecum of YBGs and investigated the variability of gut microbiota of different intestinal segments. A total of 12 intestinal samples were collected from YBGs for high-throughput sequencing analysis based on V3-V4 variable region of 16S rRNA genes. Our results revealed alterations in gut microbial composition with obvious differences in relative abundance between the different intestinal segments. Additionally, small intestine including duodenum, jejunum and ileum not only displayed higher species abundance and diversity than cecum but also showed a significant difference among the main components of gut microbiota based on the analytical results of alpha and beta diversities. At the phylum level, Firmicutes and Proteobacteria were the most preponderant phyla in all the samples regardless of intestinal sites. Moreover, the microbiota in small intestine was significantly different from cecum, which were characterized by the higher relative abundance of Butyrivibrio_2, Megasphaera, Halomonas, Delftia, Hydrogenophaga, Limnobacter, Pseudoxanthomonas, Novosphingobium, Janibacter and Erythrobacter, whereas the levels of Butyricicoccus, unidentified_Lachnospiraceae, Fusicatenibacter, Akkermansia, Ruminococcaceae_NK4A214_group and Lactobacillus were lower. Overall, this study first characterized the profile of gut microbiota composition in different intestinal sites and provide better insight into intestinal microbial community structure and diversity of YBGs.

Runs of homozygosity analysis of South African sheep breeds from various production systems investigated using OvineSNP50k data

Background

Population history, production system and within-breed selection pressure impacts the genome architecture resulting in reduced genetic diversity and increased frequency of runs of homozygosity islands. This study tested the hypothesis that production systems geared towards specific traits of importance or natural or artificial selection pressures influenced the occurrence and distribution of runs of homozygosity (ROH) in the South African sheep population. The Illumina OvineSNP50 BeadChip was used to genotype 400 sheep belonging to 13 breeds from South Africa representing mutton, pelt and mutton and wool dual-purpose breeds, including indigenous non-descript breeds that are reared by smallholder farmers. To get more insight into the autozygosity and distribution of ROH islands of South African breeds relative to global populations, 623 genotypes of sheep from worldwide populations were included in the analysis. Runs of homozygosity were computed at cut-offs of 1–6 Mb, 6–12 Mb, 12–24 Mb, 24–48 Mb and > 48 Mb, using the R package detectRUNS. The Golden Helix SVS program was used to investigate the ROH islands.

Results

A total of 121,399 ROH with mean number of ROH per animal per breed ranging from 800 (African White Dorper) to 15,097 (Australian Poll Dorset) were obtained. Analysis of the distribution of ROH according to their size showed that, for all breeds, the majority of the detected ROH were in the short (1–6 Mb) category (88.2%). Most animals had no ROH > 48 Mb. Of the South African breeds, the Nguni and the Blackhead Persian displayed high ROH based inbreeding (F_ROH) of 0.31 ± 0.05 and 0.31 ± 0.04, respectively. Highest incidence of common runs per SNP across breeds was observed on chromosome 10 with over 250 incidences of common ROHs. Mean proportion of SNPs per breed per ROH island ranged from 0.02 ± 0.15 (island ROH224 on chromosome 23) to 0.13 ± 0.29 (island ROH175 on chromosome 15). Seventeen (17) of the islands had SNPs observed in single populations (unique ROH islands). The MacArthur Merino (MCM) population had five unique ROH islands followed by Blackhead Persian and Nguni with three each whilst the South African Mutton Merino, SA Merino, White Vital Swakara, Karakul, Dorset Horn and Chinese Merino each had one unique ROH island. Genes within ROH islands were associated with predominantly metabolic and immune response traits and predomestic selection for traits such as presence or absence of horns.

Conclusions

Overall, the frequency and patterns of distribution of ROH observed in this study corresponds to the breed history and implied selection pressures exposed to the sheep populations under study.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07314-2.

Genome-Wide Analysis Revealed Homozygosity and Demographic History of Five Chinese Sheep Breeds Adapted to Different Environments

Homozygosity of long sequence genotypes are a result of parents transmitting identical haplotypes, which can be used to estimate their auto-zygosity. Therefore, we used high-density SNP Chip data to characterize the auto-zygosity of each breed according to the occurrence and distribution of runs of homozygosity (ROH). Subsequently, we identified the genomic regions with high runs of homozygosity frequencies within individuals of each breed. We selected 96 sheep samples from five local Chinese sheep breeds belonging to different geographical locations. We identified 3046 ROHs within the study breed individuals, among which the longer segments (>1–5 Mb) were dominant. On average, ROH segments covered about 12% of the genomes; the coverage rate of OAR20 was the lowest and that of OAR2 was the highest. The distribution analysis of runs of homozygosity showed that the detected ROH mainly distributed between >26 and 28 Mb. The Hetian and Hu sheep showed the lowest ROH distribution. The estimation of homozygosity level reflects the history of modern and ancient inbreeding, which may affect the genomes of Chinese indigenous sheep breeds and indicate that some animals have experienced recent self-pollination events (Yabuyi, Karakul and Wadi). In these sheep breeds, the genomic regions were assumed to be under selection signatures frequently in line with long ROH. These regions included candidate genes associated with disease resistance traits (5S_rRNA), the innate and adaptive immune response (HERC2 and CYFIP1), digestion and metabolism (CENPJ), growth (SPP1), body size and developments (GJB2 and GJA3). This study highlighted new insights into the ROH patterns and provides a basis for future breeding and conservation strategies of Chinese sheep breeds.

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.

Conservation Assessment of the State Goat Farms by Using SNP Genotyping Data

Conservation of genetic resources is of great concern globally to maintain genetic diversity for sustainable food security. Comprehensive identification of the breed composition, estimation of inbreeding and effective population size are essential for the effective management of farm animal genetic resources and to prevent the animals from genetic erosion. The Zhongwei male (ZWM), Arbas Cashmere male (ACM) and Jining Grey male (JGM) goats are conserved in three different state goat farms in China but their family information, level of inbreeding and effective population size are unknown. We investigated the genomic relationship, inbreeding coefficient and effective population size in these three breeds from three state goat farms using the Illumina goat SNP50 BeadChip. Genomic relationships and phylogenetic analysis revealed that the breeds are clearly separated and formed separate clusters based on their genetic relationship. We obtained a high proportion of informative SNPs, ranging from 91.8% in the Arbas Cashmere male to 96.2% in the Jining Grey male goat breeds with an average mean of 96.8%. Inbreeding, as measured by F_ROH, ranged from 1.79% in ZWM to 8.62% in ACM goat populations. High F_ROH values, elevated genomic coverage of very long ROH (>30 Mb) and severe decline in effective population size were recorded in ACM goat farm. The existence of a high correlation between F_HOM and F_ROH indicates that F_ROH can be used as an alternative to inbreeding estimates in the absence of pedigree records. The N_e estimates 13 generations ago were 166, 69 and 79 for ZWM, ACM and JGM goat farm, respectively indicating that these goat breeds were strongly affected by selection pressure or genetic drift. This study provides insight into the genomic relationship, levels of inbreeding and effective population size in the studied goat populations conserved in the state goat farms which will be valuable in prioritizing populations for conservation and for developing suitable management practices for further genetic improvement of these Chinese male goats.

Linkage Disequilibrium, Effective Population Size and Genomic Inbreeding Rates in American Mink Using Genotyping-by-Sequencing Data

Knowledge of linkage disequilibrium (LD) patterns is necessary to determine the minimum density of markers required for genomic studies and to infer historical changes as well as inbreeding events in the populations. In this study, we used genotyping-by-sequencing (GBS) approach to detect single nucleotide polymorphisms (SNPs) across American mink genome and further to estimate LD, effective population size (Ne), and inbreeding rates based on excess of homozygosity (F_HOM) and runs of homozygosity (ROH). A GBS assay was constructed based on the sequencing of ApeKI-digested libraries from 285 American mink using Illumina HiSeq Sequencer. Data of 13,321 SNPs located on 46 scaffolds was used to perform LD analysis. The average LD (r² ± SD) between adjacent SNPs was 0.30 ± 0.35 over all scaffolds with an average distance of 51 kb between markers. The average r² < 0.2 was observed at inter-marker distances of >40 kb, suggesting that at least 60,000 informative SNPs would be required for genomic selection in American mink. The Ne was estimated to be 116 at five generations ago. In addition, the most rapid decline of population size was observed between 100 and 200 generations ago. Our results showed that short extensions of homozygous genotypes (500 kb to 1 Mb) were abundant across the genome and accounted for 33% of all ROH identified. The average inbreeding coefficient based on ROH longer than 1 Mb was 0.132 ± 0.042. The estimations of F_HOM ranged from −0.44 to 0.34 among different samples with an average of 0.15 over all individuals. This study provided useful insights to determine the density of SNP panel providing enough statistical power and accuracy in genomic studies of American mink. Moreover, these results confirmed that GBS approach can be considered as a useful tool for genomic studies in American mink.