Estimates of genomic inbreeding and identification of candidate regions that differ between Chinese indigenous sheep breeds

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.

Genome-wide detection of runs of homozygosity and heterozygosity in Tunchang pigs

Tunchang pigs, mainly distributed throughout Hainan Province of China, are well-known for their superior meat quality, crude feed tolerance, and adaptability to high temperatures and humidity. Runs of homozygosity (ROH) can provide valuable information about the inbreeding coefficient in individuals and selection signals that may reveal candidate genes associated with key functional traits. Runs of heterozygosity (ROHet) are commonly associated with balance selection, which can help us understand the adaptive evolutionary history of domestic animals. In this study, we investigated ROHs and ROHets in 88 Tunchang pigs. We also compared the estimates of inbreeding coefficients in individuals calculated based on four methods. In summary, we detected a total of 16 ROH islands in our study, and 100 genes were found within ROH regions. These genes were correlated with economically important traits such as reproduction (e.g., SERPIND1, HIRA), meat quality (e.g., PI4KA, TBX1), immunity (e.g., ESS2, RANBP1), adaption to heat stress (TXNRD2 and DGCR8), and crude food tolerance (TRPM6). Moreover, we discovered 18 ROHet islands harbouring genes associated with reproduction (e.g., ARHGEF12, BMPR2), immune system (e.g., BRD4, DNMT3B). These findings may help us design effective breeding and conservation strategies for this unique breed.

Integrating genomics and transcriptomics to identify candidate genes for high-altitude adaptation and egg production in Nixi chicken

1. This study combined genome-wide selection signal analysis with RNA-sequencing to identify candidate genes associated with high altitude adaptation and egg production performance in Nixi chickens (NXC).2. Based on the whole-genome data from 20 NXC (♂:10; ♀:10), the population selection signal was analysed by sliding window analysis. The selected genes were screened by combination with the population differentiation statistic (FST). The sequence diversity statistic (θπ). RNA-seq was performed on the ovarian tissues of NXC (n = 6) and Lohmann laying hens (n = 6) to analyse the differentially expressed genes (DEGs) between the two groups. The functional enrichment analysis of the selected genes and differentially expressed genes was performed.3. There were 742 genes under strong positive selection and 509 differentially expressed genes screened in NXC. Integrated analysis of the genome and transcriptome revealing 26 overlapping genes. The candidate genes for adaptation to a high-altitude environment, as well as for egg production, disease resistance, vision and pigmentation in NXC were preliminarily screened.4. The results provided theoretical guidance for further research on the genetic resource protection and utilisation of NXC.

Analysis of 206 whole-genome resequencing reveals selection signatures associated with breed-specific traits in Hu sheep

As an invaluable Chinese sheep germplasm resource, Hu sheep are renowned for their high fertility and beautiful wavy lambskins. Their distinctive characteristics have evolved over time through a combination of artificial and natural selection. Identifying selection signatures in Hu sheep can provide a straightforward insight into the mechanism of selection and further uncover the candidate genes associated with breed-specific traits subject to selection. Here, we conducted whole-genome resequencing on 206 Hu sheep individuals, each with an approximate 6-fold depth of coverage. And then we employed three complementary approaches, including composite likelihood ratio, integrated haplotype homozygosity score and the detection of runs of homozygosity, to detect selection signatures. In total, 10 candidate genomic regions displaying selection signatures were simultaneously identified by multiple methods, spanning 88.54 Mb. After annotating, these genomic regions harbored collectively 92 unique genes. Interestingly, 32 candidate genes associated with reproduction were distributed in nine genomic regions detected. Out of them, two stood out as star candidates: BMPR1B and GNRH2, both of which have documented associations with fertility, and a HOXA gene cluster (HOXA1-5, HOXA9, HOXA10, HOXA11 and HOXA13) had also been linked to fertility. Additionally, we identified other genes that are related to hair follicle development (LAMTOR3, EEF1A2), ear size (HOXA1, KCNQ2), fat tail formation (HOXA10, HOXA11), growth and development (FAF1, CCNDBP1, GJB2, GJA3), fat deposition (ACOXL, JAZF1, HOXA3, HOXA4, HOXA5, EBF4), immune (UBR1, FASTKD5) and feed intake (DAPP1, RNF17, NPBWR2). Our results offer novel insights into the genetic mechanisms underlying the selection of breed-specific traits in Hu sheep and provide a reference for sheep genetic improvement programs.

Genome-wide detections for runs of homozygosity and selective signatures reveal novel candidate genes under domestication in chickens

BACKGROUND

Indigenous chickens were developed through a combination of natural and artificial selection; essentially, changes in genomes led to the formation of these modern breeds via admixture events. However, their confusing genetic backgrounds include a genomic footprint regulating complex traits, which is not conducive to modern animal breeding.

RESULTS

To better evaluate the candidate regions under domestication in indigenous chickens, we considered both runs of homozygosity (ROHs) and selective signatures in 13 indigenous chickens. The genomes of Silkie feather chickens presented the highest heterozygosity, whereas the highest inbreeding status and ROH number were found in Luhua chickens. Short ROH (< 1 Mb), were the principal type in all chickens. A total of 291 ROH islands were detected, and QTLdb mapping results indicated that body weight and carcass traits were the most important traits. An ROH on chromosome 2 covering VSTM2A gene was detected in 12 populations. Combined analysis with the Tajima's D index revealed that 18 genes (e.g., VSTM2A, BBOX1, and RYR2) were under selection and covered by ROH islands. Transcriptional analysis results showed that RYR2 and BBOX1 were specifically expressed in the heart and muscle tissue, respectively.

CONCLUSION

Based on genome-wide scanning for ROH and selective signatures, we evaluated the genomic characteristics and detected significant candidate genes covered by ROH islands and selective signatures. The findings in this study facilitated the understanding of genetic diversity and provided valuable insights for chicken breeding and conservation strategies.

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.

Genomic structure analysis and construction of DNA fingerprint for four sheep populations

The Yongdeng Qishan sheep (QS) is a sheep population found locally in China. To gain in-depth knowledge of its population characteristics, three control groups were chosen, comprising the Lanzhou fat-tailed sheep (LFT), TAN sheep (TAN), and Minxian black fur sheep (MBF), inhabiting the nearby environments. This study genotyped a total of 120 individuals from four sheep populations: QS, LFT, TAN, and MBF. Using Specific-Locus Amplified Fragment Sequencing, we conducted genetic diversity, population structure, and selective sweep analysis, and constructed the fingerprint of each population. In total, there were 782 535 single nucleotide polymorphism (SNP) variations identified, with most being situated within regions that are intergenic or intronic. The genetic diversity analysis revealed that the QS population exhibited lower genetic diversity compared to the other three populations. Consistent results were obtained from the principal component, phylogenetic tree, and population structure analysis, indicating significant genetic differences between QS and the other three populations. However, a certain degree of differentiation was observed within the QS population. The linkage disequilibrium (LD) patterns among the four populations showed clear distinctions, with the QS group demonstrating the most rapid LD decline. Kinship analysis supported the findings of population structure, dividing the 90 QS individuals into two subgroups consisting of 23 and 67 individuals. Selective sweep analysis identified a range of genes associated with reproduction, immunity, and adaptation to high-altitude hypoxia. These genes hold potential as candidate genes for marker-assisted selection breeding. Additionally, a total of 86 523 runs of homozygosity (ROHs) were detected, showing non-uniform distribution across chromosomes, with chromosome 1 having the highest coverage percentage and chromosome 26 the lowest. In the high-frequency ROH islands, 79 candidate genes were associated with biological processes such as reproduction and fat digestion and absorption. Furthermore, a DNA fingerprint was constructed for the four populations using 349 highly polymorphic SNPs. In summary, our research delves into the genetic diversity and population structure of QS population. The construction of DNA fingerprint profiles for each population can provide valuable references for the identification of sheep breeds both domestically and internationally.

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.

Genomic signatures of selection, local adaptation and production type characterisation of East Adriatic sheep breeds

BACKGROUND

The importance of sheep breeding in the Mediterranean part of the eastern Adriatic has a long tradition since its arrival during the Neolithic migrations. Sheep production system is extensive and generally carried out in traditional systems without intensive systematic breeding programmes for high uniform trait production (carcass, wool and milk yield). Therefore, eight indigenous Croatian sheep breeds from eastern Adriatic treated here as metapopulation (EAS), are generally considered as multipurpose breeds (milk, meat and wool), not specialised for a particular type of production, but known for their robustness and resistance to certain environmental conditions. Our objective was to identify genomic regions and genes that exhibit patterns of positive selection signatures, decipher their biological and productive functionality, and provide a "genomic" characterization of EAS adaptation and determine its production type.

RESULTS

We identified positive selection signatures in EAS using several methods based on reduced local variation, linkage disequilibrium and site frequency spectrum (eROHi, iHS, nSL and CLR). Our analyses identified numerous genomic regions and genes (e.g., desmosomal cadherin and desmoglein gene families) associated with environmental adaptation and economically important traits. Most candidate genes were related to meat/production and health/immune response traits, while some of the candidate genes discovered were important for domestication and evolutionary processes (e.g., HOXa gene family and FSIP2). These results were also confirmed by GO and QTL enrichment analysis.

CONCLUSIONS

Our results contribute to a better understanding of the unique adaptive genetic architecture of EAS and define its productive type, ultimately providing a new opportunity for future breeding programmes. At the same time, the numerous genes identified will improve our understanding of ruminant (sheep) robustness and resistance in the harsh and specific Mediterranean environment.

Genetic gain and inbreeding from simulation of different genomic mating schemes for pig improvement

BACKGROUND

Genomic selection involves choosing as parents those elite individuals with the higher genomic estimated breeding values (GEBV) to accelerate the speed of genetic improvement in domestic animals. But after multi-generation selection, the rate of inbreeding and the occurrence of homozygous harmful alleles might increase, which would reduce performance and genetic diversity. To mitigate the above problems, we can utilize genomic mating (GM) based upon optimal mate allocation to construct the best genotypic combinations in the next generation. In this study, we used stochastic simulation to investigate the impact of various factors on the efficiencies of GM to optimize pairing combinations after genomic selection of candidates in a pig population. These factors included: the algorithm used to derive inbreeding coefficients; the trait heritability (0.1, 0.3 or 0.5); the kind of GM scheme (focused average GEBV or inbreeding); the approach for computing the genomic relationship matrix (by SNP or runs of homozygosity (ROH)). The outcomes were compared to three traditional mating schemes (random, positive assortative or negative assortative matings). In addition, the performance of the GM approach was tested on real datasets obtained from a Large White pig breeding population.

RESULTS

Genomic mating outperforms other approaches in limiting the inbreeding accumulation for the same expected genetic gain. The use of ROH-based genealogical relatedness in GM achieved faster genetic gains than using relatedness based on individual SNPs. The G_ROH-based GM schemes with the maximum genetic gain resulted in 0.9%-2.6% higher rates of genetic gain ΔG, and 13%-83.3% lower ΔF than positive assortative mating regardless of heritability. The rates of inbreeding were always the fastest with positive assortative mating. Results from a purebred Large White pig population, confirmed that GM with ROH-based GRM was more efficient than traditional mating schemes.

CONCLUSION

Compared with traditional mating schemes, genomic mating can not only achieve sustainable genetic progress but also effectively control the rates of inbreeding accumulation in the population. Our findings demonstrated that breeders should consider using genomic mating for genetic improvement of pigs.

Estimates of genomic inbreeding and identification of candidate regions in Beijing-You chicken populations

Runs of homozygosity (ROHs) has become an effective method for analysing inbreeding in livestock populations. Moreover, ROHs is well-suited to detect signatures of selection via ROH islands. This study aimed to investigate the occurrence and distribution of ROHs, compare the genomic inbreeding coefficients and identify the genomic regions with high ROH frequencies in different Beijing-You chicken (BY) populations, including a random conservation population (BY_R), a pedigree conservation population (BY_P), and a commercial population obtained from the market (BY_C). Among them, BY_R in 2010 and 2019 were BY_R1 and BY_R2 respectively. A total of 27 916 ROHs were identified. The average number of ROHs per individual across the three BY populations ranged from 213 (BY_P) to 161 (BY_C), and the average length of ROHs ranged from 0.432 Mb (BY_R2) to 0.451 Mb (BY_P). The highest inbreeding coefficient calculated based on ROHs (F_ROH ) was 0.1 in BY_P, whereas the lowest F_ROH was 0.0743 in BY_C. In addition, the inbreeding coefficient of BY_R2 (F_ROH  = 0.0798) was higher than that of BY_R1 (F_ROH  = 0.0579). Furthermore, the highest proportion of long ROH fragments (>4 Mb) was observed in BY_P and BY_C. This study showed the top 10 ROH islands of each population, and these ROH islands harboured 53 genes, some of which were related to limb development, body size and immune response. These findings contribute to the understanding of genetic diversity and population demography, and might help improve breeding and conservation strategies for BY populations.

Genome-wide assessment of the population structure and genetic diversity of four Portuguese native sheep breeds

As the effects of global warming become increasingly complex and difficult to manage, the conservation and sustainable use of locally adapted sheep breeds are gaining ground. Portuguese native sheep breeds are important reservoirs of genetic diversity, highly adapted to harsh environments and reared in low input production systems. Genomic data that would describe the breeds in detail and accelerate the selection of more resilient animals to be able to cope with climatic challenges are still lacking. Here, we sequenced the genomes of 37 animals from four Portuguese native sheep breeds (Campaniça, Bordaleira Serra da Estrela, Merino Branco and Merino Preto) and 19 crossbred sheep to make inferences on their genomic diversity and population structure. Mean genomic diversities were very similar across these breeds (.30 ≤ H_o ≤ .34; .30 ≤ H_e ≤ .35; 1.7 × 10^-3 ≤ π ≤ 3.1 × 10^-3) and the levels of inbreeding were negligible (.005 ≤ F_IS ≤ .038). The Principal Components, Bayesian clustering and Treemix analyses split the Portuguese breeds in two main groups which are consistent with historical records: one comprising Campaniça and Serra da Estrela together with other European and transboundary dairy breeds; and another of the well-differentiated multi-purpose Merino and Merino-related breeds. Runs of homozygosity analyses yielded 1,690 ROH segments covering an average of 2.27 Gb across the genome in all individuals. The overall genome covered by ROH segments varied from 27,75 Mb in Serra da Estrela to 61,29 Mb in Campaniça. The phylogenetic analysis of sheep mitogenomes grouped the Portuguese native breeds within sub-haplogroup B1a along with two animals of the Akkaraman breed from Turkey. This result provides additional support to a direct influence of Southwest Asian sheep in local breeds from the Iberian Peninsula. Our study is a first step pertaining to the genomic characterization of Portuguese sheep breeds and the results emphasize the potential of genomic data as a valid tool to guide conservation efforts in locally adapted sheep breeds. In addition, the genomic data we generated can be used to identify markers for breed assignment and traceability of certified breed-products.

Genetic diversity of two native sheep breeds by genome-wide analysis of single nucleotide polymorphisms

Wallachian and Sumava sheep are autochthonous breeds that have undergone a significant bottleneck effect and subsequent restoration efforts. The first objective of this study was to evaluate the degree of genetic variability of both breeds and, therefore, the current management of the breeding. The second was to determine whether these two breeds still retain their genetic uniqueness in relation to each other and other breeds, despite regenerative interventions. Our data consisted of 48 individuals of Sumava and 37 individuals of Wallachian sheep. The comparison data contained 25 other breeds (primarily European) from the HapMap dataset generated by the International Sheep Genomics Consortium. When comparing all 27 breeds, the Czech breeds clustered with 15 other breeds and formed a single branch with them according to Nei's distances. At the same time, however, the clusters of both breeds were integral and easily distinguishable from the others when displayed with principal component analysis (PCA). Population substructure analysis did not show any common genetic ancestry of the Czech national breeds and breeds used for regeneration or, eventually, breeds whose ancestral population was used for regeneration. The average values of F_ST were higher in Wallachian sheep (F_ST = 0.14) than in Sumava sheep (F_ST = 0.08). The linkage disequilibrium (LD) extension per autosome was higher in Wallachian than in Sumava sheep. Consequently, the N_e estimates five generations ago were 68 for Sumava versus 34 for Wallachian sheep. Both native Czech breeds exhibit a wide range of inbreeding based on the excess of homozygosity (F_HOM) among individuals, from -0.04 to 0.16 in Sumava and from -0.13 to 0.12 in Wallachian. Average inbreeding based on runs of homozygosity was 0.21 in Sumava and 0.27 in Wallachian. Most detected runs of homozygosity (ROH) were less than 5 Mb long for both breeds. ROH segments longer than 15 Mb were absent in Wallachian sheep. Concerning putative selection signatures, a total of 471 candidate genes in Wallachian sheep within 11 hotspots and 653 genes within 13 hotspots in Sumava sheep were identified. Czech breeds appear to be well differentiated from each other and other European breeds. Their genetic diversity is low, especially in the case of the Wallachian breed. Sumava is not so threatened by low diversity but has a larger share of the non-native gene pool.

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.

Within- and between-Breed Selection Signatures in the Original and Improved Valachian Sheep

This study explored the genomic diversity and selection signatures in two Slovakian national breeds, the Original Valachian and the Improved Valachian sheep. As they are an important animal genetic resource within the country, but with decreasing population size, our aim is to identify potentially valuable genomic regions. A total of 97 sheep (18 male and 79 female) from the Original Valachian, and 69 sheep (25 male and 44 female) from the Improved Valachian populations were genotyped using the GeneSeek GGP Ovine 50 K chip. The inbreeding levels were assessed with runs of homozygosity (ROH). The selection signatures within breeds were identified based on the top 1% of most homozygous regions within the breed, the so-called ROH islands. The selection signatures between breeds were assessed based on variance in linkage disequilibrium. Overall, we have identified selection signatures with quantitative trait loci (QTL) and genes pointing towards all three production purposes of the Valachian sheep, milk, meat, and wool, including their quality characteristics. Another group with apparent large importance was the various traits related to health and resistance to parasites, which is well in line with the sturdy nature of this breed.

Genome-wide run of homozygosity analysis reveals candidate genomic regions associated with environmental adaptations of Tibetan native chickens

Background

In Tibet, the two most important breeds are Tibetan chicken and Lhasa white chicken, and the duo exhibit specific adaptations to the high altitude thereby supplying proteins for humans living in the plateau. These breeds are partly included in the conservation plans because they represent important chicken genetic resources. However, the genetic diversity of these chickens is rarely investigated. Based on whole-genome sequencing data of 113 chickens from 4 populations of Tibetan chicken including Shigatse (SH), Nyemo (NM), Dagze (DZ) and Nyingchi (LZ), as well as Lhasa white (LW) chicken breed, we investigated the genetic diversity of these chicken breeds by genetic differentiation, run of homozygosity (ROH), genomic inbreeding and selection signature analyses.

Results

Our results revealed high genetic diversity across the five chicken populations. The linkage disequilibrium decay was highest in LZ, while subtle genetic differentiation was found between LZ and other populations (Fst ranging from 0.05 to 0.10). Furthermore, the highest ROH-based inbreeding estimate (F_ROH) of 0.11 was observed in LZ. In other populations, the F_ROH ranged from 0.04 to 0.06. In total, 74, 111, 62, 42 and 54 ROH islands containing SNPs ranked top 1% for concurrency were identified in SH, NM, DZ, LZ and LW, respectively. Genes common to the ROH islands in the five populations included BDNF, CCDC34, LGR4, LIN7C, GLS, LOC101747789, MYO1B, STAT1 and STAT4. This suggested their essential roles in adaptation of the chickens. We also identified a common candidate genomic region harboring AMY2A, NTNG1 and VAV3 genes in all populations. These genes had been implicated in digestion, neurite growth and high-altitude adaptation.

Conclusions

High genetic diversity is observed in Tibetan native chickens. Inbreeding is more intense in the Nyingchi population which is also genetically distant from other chicken populations. Candidate genes in ROH islands are likely to be the drivers of adaptation to high altitude exhibited by the five Tibetan native chicken populations. Our findings contribute to the understanding of genetic diversity offer valuable insights for the genetic mechanism of adaptation, and provide veritable tools that can help in the design and implementation of breeding and conservation strategies for Tibetan native chickens.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08280-z.