Skip to main content
Download PDF

Identification of copy number variation in Tibetan sheep using whole genome resequencing reveals evidence of genomic selection

BMC Genomics volume 24, Article number: 555 (2023) Cite this article

Abstract

Background

Copy number variation (CNV) is an important source of structural variation in the mammalian genome. CNV assays present a new method to explore the genomic diversity of environmental adaptations in animals and plants and genes associated with complex traits. In this study, the genome-wide CNV distribution characteristics of 20 Tibetan sheep from two breeds (10 Oula sheep and 10 Panou sheep) were analysed using whole-genome resequencing to investigate the variation in the genomic structure of Tibetan sheep during breeding.

Results

CNVs were detected using CNVnator, and the overlapping regions of CNVs between individual sheep were combined. Among them, a total of 60,429 CNV events were detected between the indigenous sheep breed (Oula) and the synthetic sheep breed (Panou). After merging the overlapping CNVs, 4927 CNV regions (CNVRs) were finally obtained. Of these, 4559 CNVRs were shared by two breeds, and there were 368 differential CNVRs. Deletion events have a higher percentage of occurrences than duplication events. Functional enrichment analysis showed that the shared CNVRs were significantly enriched in 163 GO terms and 62 KEGG pathways, which were mainly associated with organ development, neural regulation, immune regulation, digestion and metabolism. In addition, 140 QTLs overlapped with some of the CNVRs at more than 1 kb, such as average daily gain QTL, body weight QTL, and total lambs born QTL. Many of the CNV-overlapping genes such as PPP3CA, SSTR1 and FASN, overlap with the average daily weight gain and carcass weight QTL regions. Moreover, VST analysis showed that XIRP2, ABCB1, CA1, ASPA and EEF2 differed significantly between the synthetic breed and local sheep breed. The duplication of the ABCB1 gene may be closely related to adaptation to the plateau environment in Panou sheep, which deserves further study. Additionally, cluster analysis, based on all individuals, showed that the CNV clustering could be divided into two origins, indicating that some Tibetan sheep CNVs are likely to arise independently in different populations and contribute to population differences.

Conclusions

Collectively, we demonstrated the genome-wide distribution characteristics of CNVs in Panou sheep by whole genome resequencing. The results provides a valuable genetic variation resource and help to understand the genetic characteristics of Tibetan sheep. This study also provides useful information for the improvement and breeding of Tibetan sheep in the future.

Peer Review reports

Background

Copy number variants (CNVs) are a type of genomic polymorphism that may be an important component of phenotypic variation [1]. CNVs are defined as duplications or deletions of genomic segments that range in size from 50 base pairs (bp) to millions of base pairs (Mb) and vary between individuals or species [2]. CNVs have greater effects than single nucleotide polymorphisms (SNPs) and indels in altering gene expression, and they regulate gene expression through gene dosage and positional effects [3]. With the development of sequencing technology and the related advances in molecular biology, the genome has been increasingly understood, and a large number of genetic variants present in the genome have been discovered to have different effects on individual traits. CNVs are an important genetic component of the human genome for many diseases and a major force in human evolution [4]. In the study of human diseases, a large number of genetic diseases were caused by alterations in genome structure. It has been found that CNVs can not only cause major diseases, such as birth defects, but can also be widely distributed in the human genome to cause genetic diversity. CNVs are predicted to be associated with 15% of human genetic diseases; thus, CNVs play an important role in the generation of genetic variation and individual phenotypes in humans [5]. In addition, a wide range of phenotypes and functions in animals and plants are associated with CNVs [6,7,8]. White fur is more common in sheep, but white wool is easily dyed into other colours. CNVs affect the phenotype of pigmentation and coat colour in domestic animals. The expression level of the agouti-signalling protein (ASIP) gene was significantly increased in sheep wool capsules containing a copy of this gene, which was located directly downstream of the ITCH promoter [9, 10]. This modified ASIP duplication is associated with the typical white coat colour of domesticated sheep breeds. The strong selection signature from the ASIP and KIT loci in sheep provides conclusive evidence for the strong artificial selection of these duplicated alleles [11]. It was found that the duplication of the ASIP gene in goats can also promote white hair colour [9]. Similar to that found in sheep, the coat colour of cattle is mainly determined by the inheritance of different alleles of the KIT, ASIP, TYRP1 and MC1R genes [12]. Furthermore, studies have found strong evidence that many genes belonging to the innate and adaptive immune systems have variable copy numbers in vertebrate species, particularly in the major histocompatibility complex (MHC) genes [13].

It is well known that the evolutionary impact of livestock is due to the environment and selective pressure from tamers. Species diversity arises from the continuous selection of unique external phenotypes by people in the process of repeated breeding [14]. Thus, the external phenotype remains the most obvious result of artificial selection in our domestic species, although this phenotype often escapes simple classification at the molecular level. In recent years, research advances in genomics have resulted in new genotyping tools that allow breeders to identify specific genomic segments to be transferred from parents to offspring for a more precise manual selection of traits. In China, local sheep are divided into three ancestral groups based on their geographical distribution and their genetic origin: Mongolian sheep, Kazakh sheep and Tibetan sheep [15]. With the migration of human beings, Tibetan sheep gradually became the most widely distributed and common livestock on the Tibetan Plateau [16]. Variation in the sheep genome was the basis of widespread phenotypic differences among sheep breeds. Sheep breeds were bred both to adapt them to the local natural environment and to meet human needs. With the long-term presence of natural and artificial selection pressures, different types of domesticated sheep breeds have gradually developed from the wild sheep population [17]. Both the Panou (PO) sheep breed and the Oula (OL) sheep breed are Tibetan sheep, but their individual breeding histories are different. It was reported that Oula sheep were a local breed formed by long-term and continuous hybridization between wild Argali sheep and local Tibetan sheep [18, 19]. Compared with other Tibetan sheep breeds, Oula sheep had unique breed differentiation characteristics [20]. The Panou sheep is an improved breed that was formed after a long period of cross breeding with wild Argali sheep as the male parent and Oula sheep as the female parent using modern genomic polymerization breeding technology [21]. Little research has been reported on the genomic information of Oula sheep and Panou sheep, and there is a gap in research related to CNV in particular.

Currently, the use of whole genome sequencing (WGS) technology to identify CNVs in animal genomes has the advantages of identifying small fragments of CNVs, high coverage, and more accurate results than using aCGH (comparative genomic hybridization arrays) and SNP (single nucleotide polymorphism) arrays [22]. With the dramatic increase in sequencing capacity and the corresponding decrease in sequencing costs, whole genome sequencing data are becoming the primary source of CNV detection. This study aimed to understand the genomic characteristics of Oula sheep (local breed) and Panou sheep (cultivated breed) at the CNV level using a whole genome resequencing approach. This study will deepen the understanding of the physiology and genome characteristics of Panou sheep and provide some theoretical basis for future breeding of Tibetan sheep.

Materials and methods

Sample collection and DNA extraction

Blood samples from two populations were collected from the Maqu grassland in the area of Oula Town, Maqu County, Gannan Tibetan Autonomous Prefecture, Gansu Province, China (see Additional file 1: Table S1). The two sheep populations were distributed in different locations in the Maqu grassland, with an average number of more than 300 per population. We randomly selected 10 ewes from each population to collect their blood samples for subsequent genome sequencing. With anaesthesia, 10 ml of blood was collected from the jugular vein of each sheep with a blood collection needle. It was placed in a blood collection tube containing EDTA, stored in liquid nitrogen. The anaesthetic used was procaine, and the instructions for its use were strictly followed during the procedure. Genomic DNA for whole-genome sequencing was extracted from whole blood using a Genomic DNA Purification Kit (Thermo Fisher Scientific Inc., Waltham, MA, USA) [23].

Library construction and whole genome sequencing

Restriction endonuclease were used for biological digestion of the genomic DNA samples. Then, the repair of the terminus and dA-tail was accomplished by the Paired-End DNA Sample Prep Kit (Illumina Inc., San Diego, CA, USA), which was ligated to the paired-end DNA adaptors. Finally, a 300-bp insert was used for PCR amplification. Genomic libraries were constructed and sequenced using the Illumina HiSeq X Ten (Illumina Inc., San Diego, CA, USA) NGS platform at Genedenovo company (Guangzhou, China) to obtain 150 bp (PE150) paired-end read data. The genome sequencing coverage of all samples was 10×. Raw reads were processed to obtain high-quality clean reads following three stringent filtering standards: (1) reads with ≥ 10% unidentified nucleotides (N) were removed; (2) reads with > 50% bases having Phred quality scores of ≤ 20 were removed; and (3) reads aligned to the barcode adapter were removed [24].

Sequence alignment and variant calling

The Burrows-Wheeler Aligner [25] software was then used to map the clean reads to the Tibetan sheep reference genome CAU_O.aries_1.0 (https://www.ncbi.nlm.nih.gov/assembly/GCA_017524585.1) with the default parameters. High-quality reads were retained for variant calling by removing multiple aligned reads. Using the Markduplicates function of Picard-2.9.2 software to remove PCR duplicate reads [26]. Genome AnalysisToolkit-3.3.0 was used to realign reads to correct errors caused by InDels.

CNV and CNVR detection

CNVnator-v0.4 program [27] was used for CNV prediction relative to the CAU_O.aries_1.0 reference assembly, and the optimal bin size was unified to 300 bp for all samples. The CNV calls were filtered with p value < 0.01, q0 (zero mapping quality) < 0.5, and size > 1 kb. Estimated copy number per region and copy number between breeds were calculated using the ‘-genotype’ parameter of CNVnator. The ‘CNV_overlap.py’ script on GitHub (https://github.com/bjtrost/TCAG-WGS-CNV-workflow) was used to obtain CNVRs [28]. CNVR could only be used for downstream analyses if it appears in three or more individuals per breed to remove false-positive results [29].

CNVR annotation

All CNVRs were aligned to the sheep reference genome, and candidate genes were functionally annotated using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis [30, 31]. GO term categories included molecular function, biological process and cellular component, with a significance level of p value < 0.05 for each category. In addition, quantitative trait loci (QTL) of sheep were downloaded (https://www.animalgenome.org/cgi-bin/QTLdb/OA/index) and compared with the identified CNVRs. The Bedtools-v2.27.1 [32] ‘intersect’ command was used to detect the QTLs that overlapped with the identified CNVRs.

Breed differentiation

Three methods were used to detect breed differentiation based on variation in CNVRs: cluster analysis, VST statistic and ANOVA test. All sheep individuals were analysed by clustering using the pvcluster R program [33]. Pvclust can perform both clustering and branching point tests, making the clustering results more stable and reliable. In addition, it can provide more detailed statistical information, such as bootstrap confidence and p value of clustering branches, which helps to interpret the reliability and significance of clustering results [34, 35]. A scoring matrix coded as 0 or 1 was used, depending on the presence or absence of CNVRs in each individual. Approximate unbiased (AU) p-values and bootstrap probability (BP) values were calculated by unweighted pair group averaging method (UPGMA) based on the clustering edges shown in red and green, respectively [33]. The highly differentiated regions between the two sheep populations were detected using the VST statistic method [36]. Only protein-coding genes identified in CNVRs could be used for VST calculations to study the effects of CNVRs on each breed. We calculated the VST using VST = (VT-VS)/VT, where VT is the total variance in copy number between the two sheep breeds and VS is the mean variance within each species [37]. Differences in genes copy number between varieties were calculated using the ANOVA test. Tukey’s HSD test was selected for post-hoc analysis since the sample size was the same for each species (10 samples per species) [26].

Results

Identification of CNVs and CNVRs

The sequencing data were mapped to the CAU_O.aries_1.0 reference set, and the read mapping statistics are shown in Table S2 (see Additional file 2: Table S2). The average mapping rate was 99.70%, and the sequencing depth was 10.82×. An average of 1788 deletions and 1232 duplications were detected on the autosomes, with an average deletion size of 20.91 kb and an average duplication size of 19.66 kb (Table 1).

Table 1 Statistical summary of detected CNVs
Full size table

We finally obtained 4927 CNVRs in 20 samples. Of these, 3410 were deleted regions ranging in size from 1 kb to 1650.2 kb, and 1517 were duplicated regions ranging in size from 1.3 kb to 368 kb. The numbers of breed-specific CNVRs were 185 (OL) and 183 (PO), (Fig. 1A). Among all loci, OAR 11 had the largest number of CNVRs, with 294 deletions and 188 duplications (see Additional file 3: Table S3). Approximately 3137 deletions were between 1 kb and 50 kb in size, of which 1954 (62.28%) were less than 10 kb. A total of 1377 duplications were between 1.3 kb and 50 kb in size, and 735 (53.37%) were less than 10 kb (Fig. 1B). The Fig. 1C shows the density of genes and distribution of all CNVRs used for analysis. The number of CNVRs on chromosomes did not differ significantly between the two breeds. In general, the number of detected CNVRs was proportional to chromosome size (Fig. 1D).

Fig. 1
figure 1

General description of the identified copy number variation regions (CNVRs). (A) Venn diagram for overlapping CNVRs between Oula (OL) and Panou (PO) breeds; (B) distribution of CNVRs size by state; (C) distribution of CNVRs on chromosome ideogram according to their state (deletion and duplication). This was drawn using Rideogram R package. The color painted on the chromosome represents the gene density; (D) chromosomal distribution of CNVRs for two Tibetan sheep breeds and length of Tibetan sheep autosome.

Full size image

Annotation of CNVR

We downloaded the Sheep Quantitative Trait Locus Database (Sheep QTLdb, https://www.animalgenome.org/cgi-bin/QTLdb/OA/index) from the Animal QTLdb database to detect if any identified QTLs associated with CNVRs. A total of 520 QTLs that had been reported in OAR version 3.1 were used for homology comparison (see Additional file 4: Table S4). The results showed that 140 QTLs overlapped with some of the CNVRs by more than 1 kb. Production_QTLs were the most frequently detected, with a total of 26 QTLs overlapping with CNVRs. The most detected Production_QTL was average daily gain, followed by body weight and total lambs born (see Additional file 5: Table S5). A total of 2215 protein-coding genes were annotated in 4927 CNVRs. Approximately 137 CNVRs were annotated to only one gene. At position 3.33 Mb-3.54 Mb in the location of OAR24, up to 27 genes were annotated (see Additional file 6: Table S6).

Functional enrichment analysis of GO (Gene Ontology) terms and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways was performed by DAVID. Protein-coding genes were enriched in a total of 128 biological processes, 17 cellular components, and 18 molecular functions (Additional file 7: Table S7). In addition, 62 KEGG pathways were significantly enriched (p value < 0.05) (Additional file 8: Table S8). The GO terms and KEGG pathways were mainly involved in organ development, neural regulation, immune regulation and digestive metabolism. The GO terms included biological regulation (GO:0065007), regulation of biological processes (GO:0050789), regulation of cellular processes (GO:0050794), developmental processes (GO:0032502) and cellular communication (GO:0007154). Among the KEGG pathway categories, signal transduction (ko04014: Ras signalling pathway, ko04630: Jak-STAT signalling pathway, ko04310: Wnt signalling pathway), skin pigmentation (ko04916: melanogenesis), membrane transport (ko02010: ABC transporters), and digestive system (ko04974: protein digestion and absorption) were significantly enriched.

Breed differentiation of CNVR

In the two breeds, the hierarchical clustering map was drawn using the scoring matrix to determine whether CNVRs exist in each individual. This matrix indicates the presence or absence of CNVR in an individual and calculates the p value by multiscale bootstrap resampling. The results are shown together with the Approximately Unbiased (AU) p value and Bootstrap Probability (BP) value, and it was confirmed that they were significantly clustered within each breed. The clustering results show that there is some stratification between Panou and Oula sheep (Fig. 2). However, there were some samples that were not well clustered, such as OL-3, OL-8, PO-7, and PO-9. We speculate that OL-3, OL-8, PO-7 and PO-9 belong to heterogeneous individuals present within the group, and their characteristics may be more similar to another group. This phenomenon is quite normal and suggests that the characteristics of individuals are not entirely determined by the group to which they belong but may be influenced by other factors. These factors may include the individual’s genes, environment, etc. [38].

Fig. 2
figure 2

Hierarchical clustering analysis according to the presence or absence of CNVR.

Full size image

All CNVRs were compared to the sheep genome database, and VST analysis, ANOVA and Tukey’s test were performed for each CNVR-related gene to identify population-specific selection signatures. The Manhattan plot shows the results of VST, with chromosomes in the horizontal coordinates and VST values in the vertical coordinates. To observe genes with high differentiation between different breeds, the top 98th [26] percentile was taken as the threshold line (Fig. 3) and genes above this threshold were functionally annotated. The results showed that 41 genes exceeded the threshold in the two breeds (Additional file 9: Table S9). Subsequently, ANOVA and Tukey’s test were used. As a result, 27 genes showed statistically significant differences (p value < 0.01) among the varieties (Additional file 10: Table S10). Through functional analysis of the genes above the threshold detection, five candidate genes, XIRP2, ABCB1, CA1, ASPA and EEF2, related to growth and development, digestion and metabolism, environmental adaptation and reproduction were selected (Table 2).

Fig. 3
figure 3

Manhattan plot of genome-wide genes’ VST value. The red line represents the threshold line

Full size image
Table 2 CNVR-harboring genes showing high divergence in two population pairs
Full size table

Discussion

Sheep are raised for a variety of purposes, including meat, wool and skin [39]. Sheep can be bred to promote increased productivity, and various forms of variations have been identified in the genome. CNVs affect animal production traits in various ways, such as by deleting genes or altering gene dosage [40]. Among the CNVs identified in this study, “deletion” events were more common than “duplication” events, which is consistent with observations in studies of other species [41]. The results of cluster and VST analyses showed that there were copy number differences in many genes between the two breeds and that there was some divergence between the Panou and Oula sheep populations. However, the two breeds were relatively close in genetic evolutionary relationships, probably because the two breeds share a common wild lineage ancestor. In our previous PCA results on Oula and Panou sheep, principal components could clearly distinguish Panou sheep from Oula sheep. The minimum number of splits K = 1 in the results of the population structure analysis, indicating that the two breeds are not currently producing a larger population structure. When K = 2 corresponds to our maximum number of clusters, there are two different genetic lineages between the two breeds. Thus, the populations are in a state of low differentiation from each other. When the K value was greater than 2, only four individuals, OL-1, OL-5, PO-2, and PO-6, consistently maintained a separate genetic structural cluster, while all other individuals exhibited more complex genetic structural clusters. Thus, the current artificial and environmental selection pressure is causing the Oula and Panou sheep populations to gradually and continuously diverge in different directions [21].

Genome-wide CNV detection is also a strategy to identify potential key genes for traits of interest by mining genes within CNVRs in a specific experimental design. Even at low coverage, CNV detection using the WGS-based method exhibits better performance than the array-based method [42]. Therefore, we used the different CNVRs between the local breed (Oula sheep) and the cultivated breed (Panou sheep) for the identification of production performance candidate genes to provide a theoretical basis for the subsequent continuous selection of Tibetan sheep. However, due to the possibility of false-positives, we analysed CNVRs for at least three individuals in each breed. We identified 4927 CNVRs in the two breeds, including 3410 deletions and 1517 duplications, and there were a large number of identical CNVRs. This is probably due to their common evolutionary history prior to artificial selection. In this study, the number of CNVRs had a significant positive linear relationship with the length of chromosomes. This result is also consistent with the reports of others [43]. In addition, similar to the results of Huang et al., the CNV assay showed more deletion events than duplication events [44]. Therefore, we speculate that this may be a characteristic of Tibetan sheep. GO enrichment analysis showed that the detected CNVRs overlapping genes common to the two sheep breeds were significantly enriched in organ development, neural regulation, immune regulation and digestive metabolism. The functions related to immune regulation, metabolism and growth and development were also significantly enriched in yak [45] and other sheep breeds [46]. Like the yak, the Oula and Panou sheep live in high altitudes with harsh climates. The strong immune system and the metabolic capacity of the organisms facilitate their adaptation to cold environments with low oxygen at high altitudes and reduce the occurrence of diseases [47]. In the KEGG pathway analysis, it is noteworthy that CNVR-carrying genes common to the two breeds were significantly enriched in digestive metabolism, ABC transporter proteins, disease defence and melanogenesis. In high-altitude and cold regions, the level of nutrient digestion and metabolism is important for Tibetan sheep to maintain homeostasis of their growth and development. This result is similar to that of previous studies on cattle [48] and goats [49]. Studies in mammals have found that ABC transporter proteins can carry a variety of endogenous metabolites through lipid membranes, thereby promoting the absorption and utilization of these nutrients [50]. Thus, the functional genes related to feed digestion and utilization screened in CNVRs may be related to the adaptation of Tibetan sheep to the local living environment.

Many studies have shown that CNVRs contain QTLs associated with economically important traits in animals [51]. Therefore, the CNVRs detected in this study were compared with the QTLs reported in the sheep QTL database. It was found that QTLs affecting economic traits of genetic variation in livestock could be identified by screening the relevant genes contained in genomic CNVs [52]. After integrating the CNVs into the QTLs, 140 QTLs partially overlapped the CNVRs by more than 1 kb. Many overlapping CNV genes, such as PPP3CA, SSTR1 and FASN, are located in the average daily weight gain and carcass weight QTL regions. Wan et al. found that the PPP3CA gene was involved in the regulation of skeletal muscle growth and development by inhibiting the expression of fibroblast growth factor 23 (FGF23). FGF23 plays a key role in the control of skeletal muscle fibres and is expressed in various muscle tissues [53]. The PPP3CA gene was also expressed at high levels in the biceps femoris, longus and abdominal muscles of Tianfu goats, suggesting that the PPP3CA gene has an important regulatory role in muscle growth and development [54]. SSTR1 plays a more important role than other receptors in insulin and GH secretion, especially in the maintenance of basal levels of GH [55]. Li et al. found that the SSTR1 gene was closely associated with the growth and development of Hulun Buir sheep [56]. The fatty acid synthase (FASN) gene is used to improve fatty acid composition [57], and FASN is also considered to be an important gene affecting carcass traits in sheep [58, 59]. Therefore, in this study, the candidate genes screened by CNVs could be used as relevant molecular markers for the improvement of Tibetan sheep breeds.

Selective sweeps can reveal putative regions under environmental and artificial selection during local adaptation and domestication [60]. In the genome, calculating paired VST values can be used to screen for key CNVRs that differ significantly between populations [43]. We identified a total of 368 differential CNVRs between the two breeds and further identified 41 functional genes located within these CNVRs. Through annotation analysis, it was found that these genes were enriched in aspartic acid metabolism, kinase binding, ATPase-coupled transmembrane transporter activity regulation and other processes and participated in multiple metabolic pathways. In the present study, five CNVRs carrying the XIRP2, ABCB1, CA1, ASPA and EEF2 genes showed significant pairwise differentiation between Oula and Panou sheep. The XIRP2 gene has been reported to be associated with feed efficiency, average daily gain and residual feed intake in cattle [61]. Subsequently, Seifi Moroudi et al. also demonstrated that the XIRP2 gene was associated with growth and developmental biology pathways in sheep and cattle [62]. This function of the XIRP2 gene is important for Tibetan sheep living in a relatively pasture-deficient area. A previous study found that the ABCB1 gene was hypoxia-responsive and could be regulated by HIF-1α due to a hypoxia-responsive element on its promoter [63]. Moreover, in the gut of sheep, the relative expression of the ABCB1 gene in different compartments differed greatly, and there was a certain gradient, which indicates that it plays an important role in the digestive metabolic process in sheep [64]. Berton et al. found that the ABCB1 gene was also involved in the gastrointestinal immune response and inflammatory processes in sheep [65]. Sales et al. explored the relationship between the ABCB1 gene and folliculogenesis by vitrification and in vitro culture of sheep ovarian tissue, indicating that the ABCB1 gene may affect the development of ewe follicles [66]. Thus, it seems that the ABCB1 gene plays an important role in the adaptation of Tibetan sheep to the plateau environment. Carbonic anhydrase 1 (CA1) is an enzyme abundant in colonic epithelial cells and plays a role in the maintenance of cellular acid-base homeostasis [67]. CA1 is associated with many physiological and pathological processes, including the development of various digestive disorders, gluconeogenesis and adipogenesis. In addition, polymorphic loci in the CA1 gene could be used as potential genetic markers for improving feed efficiency in sheep [68]. The CA1 gene is important for Tibetan sheep living at high altitudes to maintain their growth and metabolic balance, enhance disease resistance and improve feed utilization. Melatonin secretion is closely related to oestrous in sheep, and the EEF2 gene induces melatonin secretion from the pituitary gland to participate in the oestrous of sheep [69]. Zhang et al. found that the relative expression of the ASPA gene in the hypothalamus changed significantly during oestrous and was associated with litter size in sheep [70]. During domestication, several genes associated with the physical characteristics of sheep were artificially selected in a targeted manner. Therefore, under these selection pressures, CNVs may accumulate in sheep breeds, thus forming the genetic basis of important economic characteristics. In this study, the Panou sheep is a cultivated breed with the local Oula Tibetan sheep as the hybrid female and its wild-line ancestor, the Argali sheep, as the hybrid male, after a long period of adaptation to the low oxygen environment of the plateau and artificial selection. Therefore, we believe that the CNVRs of Panou sheep carry five genes, XIRP2, ABCB1, CA1, ASPA and EEF2, which are also subject to strict selection pressure during breed formation. This may help them better adapt to the low oxygen and cold environment in high-altitude areas and the lack of growth conditions for pasture, which is in line with our initial breeding goals.

Conclusions

In this study, based on whole genome resequencing of 20 Tibetan sheep, we explored the genomic characteristics and structural variation of two breeds of Tibetan sheep and 368 CNVRs were obtained. By analysing the functions of the annotated genes in the CNVRs and the genes showing copy number differences among varieties, we found that the genetic variation caused by artificial selection was effective at the CNV level. Some CNVRs carry differential genes, such as XIRP2, ABCB1, CA1, ASPA and EEF2, that may be related to body weight, environmental adaptation and reproductive capacity. Our research results provide valuable resources of the genome-wide variation for Tibetan sheep and are useful for elucidating the genetic basis of plateau adaptation in Tibetan sheep. Furthermore, these results will contribute greatly to the future breeding and development of economic traits in Tibetan sheep.

Data Availability

The datasets generated in this paper can be found at Sequence Read Archive: PRJNA797957. Link address: https://www.ncbi.nlm.nih.gov/bioproject/PRJNA797957.

References

  1. Stranger B, Forrest M, Dunning M, Ingle C, Beazley C, Thorne N, Redon R, Bird C, de Grassi A, Lee C, et al. Relative impact of nucleotide and copy number variation on gene expression phenotypes. Sci (New York NY). 2007;315(5813):848–53.

    CAS  Google Scholar 

  2. Conrad D, Pinto D, Redon R, Feuk L, Gokcumen O, Zhang Y, Aerts J, Andrews T, Barnes C, Campbell P, et al. Origins and functional impact of copy number variation in the human genome. Nature. 2010;464(7289):704–12.

    CAS  PubMed  Google Scholar 

  3. Chiang C, Scott A, Davis J, Tsang E, Li X, Kim Y, Hadzic T, Damani F, Ganel L, Montgomery S, et al. The impact of structural variation on human gene expression. Nat Genet. 2017;49(5):692–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  4. Sudmant P, Mallick S, Nelson B, Hormozdiari F, Krumm N, Huddleston J, Coe B, Baker C, Nordenfelt S, Bamshad M, et al. Global diversity, population stratification, and selection of human copy-number variation. Sci (New York NY). 2015;349(6253):aab3761.

    Google Scholar 

  5. Luo X, Qin F, Cai G, Xiao F. Integrating genomic correlation structure improves copy number variations detection. Bioinf (Oxford England). 2021;37(3):312–7.

    CAS  Google Scholar 

  6. Norris B, Whan V. A gene duplication affecting expression of the ovine ASIP gene is responsible for white and black sheep. Genome Res. 2008;18(8):1282–93.

    CAS  PubMed  PubMed Central  Google Scholar 

  7. Wright D, Boije H, Meadows J, Bed’hom B, Gourichon D, Vieaud A, Tixier-Boichard M, Rubin C, Imsland F, Hallböök F, et al. Copy number variation in intron 1 of SOX5 causes the pea-comb phenotype in chickens. PLoS Genet. 2009;5(6):e1000512.

    PubMed  PubMed Central  Google Scholar 

  8. Zhou Z, Jiang Y, Wang Z, Gou Z, Lyu J, Li W, Yu Y, Shu L, Zhao Y, Ma Y, et al. Resequencing 302 wild and cultivated accessions identifies genes related to domestication and improvement in soybean. Nat Biotechnol. 2015;33(4):408–14.

    CAS  PubMed  Google Scholar 

  9. Fontanesi L, Beretti F, Riggio V, Gómez González E, Dall’Olio S, Davoli R, Russo V, Portolano B. Copy number variation and missense mutations of the agouti signaling protein (ASIP) gene in goat breeds with different coat colors. Cytogenet Genome Res. 2009;126(4):333–47.

    CAS  PubMed  Google Scholar 

  10. Fontanesi L, Beretti F, Martelli P, Colombo M, Dall’olio S, Occidente M, Portolano B, Casadio R, Matassino D, Russo V. A first comparative map of copy number variations in the sheep genome. Genomics. 2011;97(3):158–65.

    CAS  PubMed  Google Scholar 

  11. Kijas J, Lenstra J, Hayes B, Boitard S, Porto Neto L, San Cristobal M, Servin B, McCulloch R, Whan V, Gietzen K, et al. Genome-wide analysis of the world’s sheep breeds reveals high levels of historic mixture and strong recent selection. PLoS Biol. 2012;10(2):e1001258.

    CAS  PubMed  PubMed Central  Google Scholar 

  12. Trigo B, Utsunomiya A, Fortunato A, Milanesi M, Torrecilha R, Lamb H, Nguyen L, Ross E, Hayes B, Padula R, et al. Variants at the ASIP locus contribute to coat color darkening in Nellore cattle. Genet Selection Evolution: GSE. 2021;53(1):40–51.

    CAS  PubMed Central  Google Scholar 

  13. Balakrishnan C, Ekblom R, Völker M, Westerdahl H, Godinez R, Kotkiewicz H, Burt D, Graves T, Griffin D, Warren W, et al. Gene duplication and fragmentation in the zebra finch major histocompatibility complex. BMC Biol. 2010;8:29–47.

    PubMed  PubMed Central  Google Scholar 

  14. Bickhart D, Liu G. The challenges and importance of structural variation detection in livestock. Front Genet. 2014;5:37–50.

    PubMed  PubMed Central  Google Scholar 

  15. Wei C, Wang H, Liu G, Wu M, Cao J, Liu Z, Liu R, Zhao F, Zhang L, Lu J, et al. Genome-wide analysis reveals population structure and selection in chinese indigenous sheep breeds. BMC Genomics. 2015;16:194–205.

    PubMed  PubMed Central  Google Scholar 

  16. Liu J, Ding X, Zeng Y, Yue Y, Guo X, Guo T, Chu M, Wang F, Han J, Feng R, et al. Genetic diversity and phylogenetic evolution of Tibetan Sheep Based on mtDNA D-Loop sequences. PLoS ONE. 2016;11(7):e0159308.

    PubMed  PubMed Central  Google Scholar 

  17. Li X, Yang J, Shen M, Xie X, Liu G, Xu Y, Lv F, Yang H, Yang Y, Liu C, et al. Whole-genome resequencing of wild and domestic sheep identifies genes associated with morphological and agronomic traits. Nat Commun. 2020;11(1):2815–30.

    CAS  PubMed  PubMed Central  Google Scholar 

  18. Liu J-b, Guo J, Wang F, Yue Y-j, Zhang W-l. Feng R-l, Guo T-t, Yang B-h, Sun X-p: carcass and meat quality characteristics of Oula lambs in China. Small Ruminant Research. 2015;123(2–3):251–9.

    Google Scholar 

  19. Guo X, Liu J, Zeng Y, Ding X, Bao P, Yan P, Pei J. Study on complete mitochondrial genome of Oula Sheep (Ovis aries). Agricultural Sci & Technol. 2017;18(8):1365–6.

    Google Scholar 

  20. Zhang Y, Xue X, Liu Y, Abied A, Ding Y, Zhao S, Wang W, Ma L, Guo J, Guan W, et al. Genome-wide comparative analyses reveal selection signatures underlying adaptation and production in Tibetan and Poll Dorset sheep. Sci Rep. 2021;11(1):2466.

    CAS  PubMed  PubMed Central  Google Scholar 

  21. Shi H, Li T, Su M, Wang H, Li Q, Lang X, Ma Y. Whole genome sequencing revealed genetic diversity, population structure, and selective signature of Panou Tibetan sheep. BMC Genomics. 2023;24(1):50.

    CAS  PubMed  PubMed Central  Google Scholar 

  22. Escaramís G, Docampo E, Rabionet R. A decade of structural variants: description, history and methods to detect structural variation. Brief Funct Genomics. 2015;14(5):305–14.

    PubMed  Google Scholar 

  23. Abied A, Bagadi A, Bordbar F, Pu Y, Augustino S, Xue X, Xing F, Gebreselassie G, Mwacharo J, Ma Y, et al. Genomic diversity, Population structure, and signature of selection in five chinese native Sheep Breeds adapted to Extreme environments. Genes. 2020;11(5):494–509.

    CAS  PubMed  PubMed Central  Google Scholar 

  24. Ciani E, Lasagna E, D’Andrea M, Alloggio I, Marroni F, Ceccobelli S, Delgado Bermejo JV, Sarti FM, Kijas J, Lenstra JA, et al. Merino and Merino-derived sheep breeds: a genome-wide intercontinental study. Genet Selection Evol. 2015;47:64–74.

    Google Scholar 

  25. Li H, Durbin R. Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinf (Oxford England). 2010;26(5):589–95.

    Google Scholar 

  26. Seol D, Ko B, Kim B, Chai H, Lim D, Kim H. Identification of Copy Number Variation in Domestic Chicken using whole-genome sequencing reveals evidence of selection in the genome. Animals: An open Access Journal from MDPI. 2019;9(10):809–23.

    PubMed  Google Scholar 

  27. Abyzov A, Urban A, Snyder M, Gerstein M. CNVnator: an approach to discover, genotype, and characterize typical and atypical CNVs from family and population genome sequencing. Genome Res. 2011;21(6):974–84.

    CAS  PubMed  PubMed Central  Google Scholar 

  28. Trost B, Walker S, Wang Z, Thiruvahindrapuram B, MacDonald J, Sung W, Pereira S, Whitney J, Chan A, Pellecchia G, et al. A Comprehensive Workflow for Read depth-based identification of Copy-Number variation from whole-genome sequence data. Am J Hum Genet. 2018;102(1):142–55.

    CAS  PubMed  PubMed Central  Google Scholar 

  29. Pierce M, Dzama K, Muchadeyi F. Genetic diversity of seven cattle breeds inferred using Copy Number Variations. Front Genet. 2018;9:163–76.

    PubMed  PubMed Central  Google Scholar 

  30. Huang dW, Sherman B, Lempicki R. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37(1):1–13.

    Google Scholar 

  31. Huang dW, Sherman B, Lempicki R. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.

    Google Scholar 

  32. Quinlan A, Hall I. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinf (Oxford England). 2010;26(6):841–2.

    CAS  Google Scholar 

  33. Suzuki R, Shimodaira H. Pvclust: an R package for assessing the uncertainty in hierarchical clustering. Bioinf (Oxford England). 2006;22(12):1540–2.

    CAS  Google Scholar 

  34. Davidson R, Sullivant S. Distance-Based phylogenetic methods around a polytomy. IEEE/ACM Trans Comput Biol Bioinf. 2014;11(2):325–35.

    Google Scholar 

  35. Takezaki N. Tie trees generated by distance methods of phylogenetic reconstruction. Mol Biol Evol. 1998;15(6):727–37.

    CAS  PubMed  Google Scholar 

  36. Redon R, Ishikawa S, Fitch K, Feuk L, Perry G, Andrews T, Fiegler H, Shapero M, Carson A, Chen W, et al. Global variation in copy number in the human genome. Nature. 2006;444(7118):444–54.

    CAS  PubMed  PubMed Central  Google Scholar 

  37. Wu Q, Zhou Y, Wang Y, Zhang Y, Shen Y, Su Q, Gao G, Xu H, Zhou X, Liu B. Whole-genome sequencing reveals breed-differential CNVs between Tongcheng and large White pigs. Anim Genet. 2020;51(6):940–4.

    CAS  PubMed  Google Scholar 

  38. McHugh C, Brown L, Thornton T. Detecting heterogeneity in Population structure across the genome in Admixed populations. Genetics. 2016;204(1):43–56.

    PubMed  PubMed Central  Google Scholar 

  39. Pei WY, Yun J, Rong EG, Yang H, Wang ZP, Wang SZ, Hui LI, Wang N. Promoter characterization of Sheep Dlx3 Gene and Association of promoter polymorphisms with wool quality traits in chinese Merino. Scientia Agricultura Sinica. 2013;46(3):614–22.

    CAS  Google Scholar 

  40. Ma, Liu XX, Pan JF, He LNYH. Genome-wide detection of copy number variation in chinese indigenous sheep using an ovine high-density 600 K SNP array. Sci Rep Uk. 2017;7:912–21.

    Google Scholar 

  41. Stafuzza N, Silva R, Fragomeni B, Masuda Y, Huang Y, Gray K, Lourenco D. A genome-wide single nucleotide polymorphism and copy number variation analysis for number of piglets born alive. BMC Genomics. 2019;20(1):321–31.

    PubMed  PubMed Central  Google Scholar 

  42. Zhou B, Ho S, Zhang X, Pattni R, Haraksingh R, Urban A. Whole-genome sequencing analysis of CNV using low-coverage and paired-end strategies is efficient and outperforms array-based CNV analysis. J Med Genet. 2018;55(11):735–43.

    CAS  PubMed  Google Scholar 

  43. Yuan C, Lu Z, Guo T, Yue Y, Wang X, Wang T, Zhang Y, Hou F, Niu C, Sun X, et al. A global analysis of CNVs in chinese indigenous fine-wool sheep populations using whole-genome resequencing. BMC Genomics. 2021;22(1):78–87.

    CAS  PubMed  PubMed Central  Google Scholar 

  44. Huang Y, Li Y, Wang X, Yu J, Cai Y, Zheng Z, Li R, Zhang S, Chen N, Asadollahpour Nanaei H, et al. An atlas of CNV maps in cattle, goat and sheep. Sci China Life Sci. 2021;64(10):1747–64.

    CAS  PubMed  Google Scholar 

  45. Wang H, Chai Z, Hu D, Ji Q, Xin J, Zhang C, Zhong J. A global analysis of CNVs in diverse yak populations using whole-genome resequencing. BMC Genomics. 2019;20(1):61–72.

    PubMed  PubMed Central  Google Scholar 

  46. Toremurat Z, Ibrahim EE, Huang Y-Z, Lan X, Pi L, Chaogetu B, Hu L, Chen H. Copy number variations of TOP2B gene are associated with growth traits in chinese sheep breeds. Animal Biotechnol. 2022;33(1):85–9.

    CAS  Google Scholar 

  47. Qiu Q, Zhang G, Ma T, Qian W, Wang J, Ye Z, Cao C, Hu Q, Kim J, Larkin D, et al. The yak genome and adaptation to life at high altitude. Nat Genet. 2012;44(8):946–9.

    CAS  PubMed  Google Scholar 

  48. Liu G, Brown T, Hebert D, Cardone M, Hou Y, Choudhary R, Shaffer J, Amazu C, Connor E, Ventura M, et al. Initial analysis of copy number variations in cattle selected for resistance or susceptibility to intestinal nematodes. Mammalian Genome: Official Journal of the International Mammalian Genome Society. 2011;22:111–21.

    PubMed  Google Scholar 

  49. Guan D, Martínez A, Castelló A, Landi V, Luigi-Sierra M, Fernández-Álvarez J, Cabrera B, Delgado J, Such X, Jordana J, et al. A genome-wide analysis of copy number variation in Murciano-Granadina goats. Genet Selection Evolution: GSE. 2020;52(1):44–53.

    CAS  PubMed Central  Google Scholar 

  50. Dron J, Wang J, Berberich A, Iacocca M, Cao H, Yang P, Knoll J, Tremblay K, Brisson D, Netzer C, et al. ABCA1Large-scale deletions of the gene in patients with hypoalphalipoproteinemia. J Lipid Res. 2018;59(8):1529–35.

    CAS  PubMed  PubMed Central  Google Scholar 

  51. Cao X, Huang Y, Ma Y, Cheng J, Qu Z, Ma Y, Bai Y, Tian F, Lin F, Ma Y, et al. Integrating CNVs into meta-QTL identified GBP4 as positional candidate for adult cattle stature. Funct Integr Genom. 2018;18(5):559–67.

    CAS  Google Scholar 

  52. Hu Z, Park C, Reecy J. Bringing the animal QTLdb and CorrDB into the future: meeting new challenges and providing updated services. Nucleic Acids Res. 2022;50:D956–61.

    CAS  PubMed  Google Scholar 

  53. Wan L, Ma J, Xu G, Wang D, Wang N. Molecular cloning, structural analysis and tissue expression of protein phosphatase 3 catalytic subunit alpha isoform (PPP3CA) gene in Tianfu goat muscle. Int J Mol Sci. 2014;15(2):2346–58.

    PubMed  PubMed Central  Google Scholar 

  54. Bär L, Großmann C, Gekle M, Föller M. Calcineurin inhibitors regulate fibroblast growth factor 23 (FGF23) synthesis. Naunyn Schmiedebergs Arch Pharmacol. 2017;390(11):1117–23.

    PubMed  Google Scholar 

  55. Kailey B, van de Bunt M, Cheley S, Johnson P, MacDonald P, Gloyn A, Rorsman P, Braun M. SSTR2 is the functionally dominant somatostatin receptor in human pancreatic β- and α-cells. Am J Physiol Endocrinol Metabolism. 2012;303(9):E1107–1116.

    CAS  Google Scholar 

  56. Ding N, Tian D, Li X, Zhang Z, Tian F, Liu S, Han B, Liu D, Zhao K. Genetic polymorphisms of IGF1 and IGF1R genes and their Effects on Growth Traits in Hulun Buir Sheep. Genes. 2022;13(4):666–79.

    CAS  PubMed  PubMed Central  Google Scholar 

  57. Zhang Y, Fan X, Qiu L, Zhu W, Huang L, Miao Y. Liver X receptor α promotes milk fat synthesis in buffalo mammary epithelial cells by regulating the expression of FASN. J Dairy Sci. 2021;104(12):12980–93.

    CAS  PubMed  Google Scholar 

  58. Wen Y, Li S, Bao G, Wang J, Liu X, Hu J, Zhao F, Zhao Z, Shi B, Luo Y. Comparative transcriptome analysis reveals the mechanism Associated with dynamic changes in Meat Quality of the Longissimus Thoracis muscle in Tibetan Sheep at different growth stages. Front Veterinary Sci. 2022;9:926725–39.

    Google Scholar 

  59. Al-Thuwaini TM, Kareem ZA. Novel missense variant L46Q of fatty acid synthase gene and fatty acids content in Awassi sheep. Acta Scientiarum Animal Sciences. 2022;44:e56273–3.

    Google Scholar 

  60. Chen Q, Qu K, Ma Z, Zhan J, Zhang F, Shen J, Ning Q, Jia P, Zhang J, Chen N, et al. Genome-wide Association Study identifies genomic loci Associated with Neurotransmitter concentration in cattle. Front Genet. 2020;11:139–48.

    CAS  PubMed  PubMed Central  Google Scholar 

  61. Connor E, Kahl S, Elsasser T, Parker J, Li R, Van Tassell C, Baldwin R, Barao S. Enhanced mitochondrial complex gene function and reduced liver size may mediate improved feed efficiency of beef cattle during compensatory growth. Funct Integr Genom. 2010;10(1):39–51.

    CAS  Google Scholar 

  62. Seifi Moroudi R, Ansari Mahyari S, Vaez Torshizi R, Lanjanian H, Masoudi-Nejad A. Identification of new genes and quantitative trait locis associated with growth curve parameters in F2 chicken population using genome-wide association study. Anim Genet. 2021;52(2):171–84.

    CAS  PubMed  Google Scholar 

  63. Comerford K, Wallace T, Karhausen J, Louis N, Montalto M, Colgan S. Hypoxia-inducible factor-1-dependent regulation of the multidrug resistance (MDR1) gene. Cancer Res. 2002;62(12):3387–94.

    CAS  PubMed  Google Scholar 

  64. Ballent M, Wilkens MR, Mate L, Muscher AS, Virkel G, Sallovitz J, Schroeder B, Lanusse C, Lifschitz A. P-glycoprotein in sheep liver and small intestine: gene expression and transport efflux activity. J Vet Pharmacol Ther. 2013;36(6):576–82.

    CAS  PubMed  Google Scholar 

  65. Berton M, da Silva R, Banchero G, Mourão G, Ferraz J, Schenkel F, Baldi F. Genomic integration to identify molecular biomarkers associated with indicator traits of gastrointestinal nematode resistance in sheep. J Anim Breed Genetics = Zeitschrift fur Tierzuchtung und Zuchtungsbiologie. 2022;139(5):502–16.

    CAS  PubMed  Google Scholar 

  66. Sales A, Duarte A, Rocha R, Brito I, Locatelli Y, Alves B, Alves K, Figueiredo J, Rodrigues A. Transcriptional downregulation of ABC transporters is related to follicular degeneration after vitrification and in vitro culture of ovine ovarian tissue. Theriogenology. 2022;177:127–32.

    CAS  PubMed  Google Scholar 

  67. Viikilä P, Kivelä A, Mustonen H, Koskensalo S, Waheed A, Sly W, Pastorek J, Pastorekova S, Parkkila S, Haglund C. Carbonic anhydrase enzymes II, VII, IX and XII in colorectal carcinomas. World J Gastroenterol. 2016;22(36):8168–77.

    PubMed  PubMed Central  Google Scholar 

  68. Zhang D, Zhang X, Li F, Yuan L, Zhang Y, Li X, Zhao Y, Song Q, Li G, Wang W. Polymorphisms in ovine ME1 and CA1 genes and their association with feed efficiency in Hu sheep. J Anim Breed Genetics = Zeitschrift fur Tierzuchtung und Zuchtungsbiologie. 2021;138(5):589–99.

    CAS  PubMed  Google Scholar 

  69. Chen S, Guo X, He X, Di R, Zhang X, Zhang J, Wang X, Chu M. FecBInsight into Pituitary lncRNA and mRNA at two estrous stages in small tail Han Sheep with different genotypes. Front Endocrinol. 2021;12:789564.

    Google Scholar 

  70. Zhang Z, Tang J, He X, Di R, Zhang X, Zhang J, Hu W, Chu M. Identification and characterization of hypothalamic alternative splicing events and variants in Ovine Fecundity-Related genes. Animals: An open Access Journal from MDPI. 2020;10(11):2111–26.

    PubMed  Google Scholar 

Download references

Acknowledgements

Thanks to all participants for their advice and support of this study. The authors also thank the reviewers for constructive comments.

Funding

This project was financed by the Education science and technology innovation project of Gansu Province (GSSYLXM-02), and National Key R&D Program of China (2021YFD1100502), and Purchase contract project of the Ministry of Agriculture and Rural Affairs (19221204), and Discipline Team Project of Gansu Agricultural University (NO: GAU-XKTD-2022-20), and Project of Research Condition Construction and Achievement Transformation of Gansu Academy of Agricultural Sciences (Modern Biological Breeding, 2021GAAS01).

Author information

Authors and Affiliations

  1. College of Animal Science and Technology, Gansu Agricultural University, Lanzhou, 730070, China

    Huibin Shi, Taotao Li, Manchun Su, Huihui Wang, Qiao Li & Youji Ma

  2. Gansu Key Laboratory of Animal Generational Physiology and Reproductive Regulation, Lanzhou, 730070, China

    Huibin Shi, Taotao Li, Manchun Su, Huihui Wang, Qiao Li & Youji Ma

  3. College of Animal Science & Technology, Henan University of Animal Husbandry and Economy, Zhengzhou, 450046, China

    Huibin Shi

  4. Institute of Animal & Pasture Science and Green Agriculture, Gansu Academy of Agricultural Science, Lanzhou, 730070, China

    Xia Lang

Authors
  1. Huibin Shi
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Taotao Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Manchun Su
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Huihui Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Qiao Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Xia Lang
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Youji Ma
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

Y.J.M. and T.T.L. led the project and designed and conceived the study. H.B.S performed data analysis, interpreted results, prepared and developed the manuscript. M.C.S. and Q.L. revised the manuscript. X.L. and H.H.W. performed sampling. All authors contributed and approved the final manuscript.

Corresponding author

Correspondence to Youji Ma.

Ethics declarations

Competing interests

The authors declare no competing interests.

Ethics approval and consent to participate

All experimental procedures in the present study were performed in accordance with the Ministry of Science and Technology of the People’s Republic of China (Approval No. 2006 − 398). All experimental animals in this study were approved by the relevant animal ethics committee (Animal Care Committee of Gansu Agricultural University, GSAU-AEW-2020-0057). All the experimental procedures were performed in accordance with the relevant animal ethics and welfare guidelines and regulations of the Animal Care Committee of Gansu Agricultural University (GSAU-AEW-2020-0057). We have complied with ARRIVE at submission.

Consent for publication

Not applicable.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Shi, H., Li, T., Su, M. et al. Identification of copy number variation in Tibetan sheep using whole genome resequencing reveals evidence of genomic selection. BMC Genomics 24, 555 (2023). https://doi.org/10.1186/s12864-023-09672-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-023-09672-z

Keywords

BMC Genomics

ISSN: 1471-2164

Contact us