Whole-genome resequencing reveals molecular imprints of anthropogenic and natural selection in wild and domesticated sheep

De-Yin Zhang, Xiao-Xue Zhang, Fa-Di Li,3, Lv-Feng Yuan, Xiao-Long Li, Yu-Kun Zhang, Yuan Zhao, Li-Ming Zhao,Jiang-Hui Wang, Dan Xu, Jiang-Bo Cheng, Xiao-Bin Yang, Wen-Xin Li, Chang-Chun Lin, Bu-Bo Zhou,Wei-Min Wang,＊

1 State Key Laboratory of Grassland Agro-Ecosystems, Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural Affairs, Engineering Research Center of Grassland Industry, Ministry of Education; College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou, Gansu 730020, China

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

3 Engineering Laboratory of Sheep Breeding and Reproduction Biotechnology in Gansu Province, Minqin, Gansu 733300, China

4 Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Sciences, Lanzhou, Gansu 730046, China

ABSTRACT The abundance of domesticated sheep varieties and phenotypes is largely the result of long-term natural and artificial selection. However, there is limited information regarding the genetic mechanisms underlying phenotypic variation induced by the domestication and improvement of sheep. In this study, to explore genomic diversity and selective regions at the genome level, we sequenced the genomes of 100 sheep across 10 breeds and combined these results with publicly available genomic data from 225 individuals, including improved breeds, Chinese indigenous breeds,African indigenous breeds, and their Asian mouflon ancestor. Based on population structure, the domesticated sheep formed a monophyletic group,while the Chinese indigenous sheep showed a clear geographical distribution trend. Comparative genomic analysis of domestication identified several selective signatures, including IFI44 and IFI44L genes and PANK2 and RNF24 genes, associated with immune response and visual function.Population genomic analysis of improvement demonstrated that candidate genes of selected regions were mainly associated with pigmentation,energy metabolism, and growth development.Furthermore, the IFI44 and IFI44L genes showed a common selection signature in the genomes of 30 domesticated sheep breeds. The IFI44 c. 54413058 C＞G mutation was selected for genotyping and population genetic validation. Results showed that the IFI44 polymorphism was significantly associated with partial immune traits. Our findings identified the population genetic basis of domesticated sheep at the whole-genome level, providing theoretical insights into the molecular mechanism underlying breed characteristics and phenotypic changes during sheep domestication and improvement.

Keywords: Sheep; Whole-genome resequencing;Selection signature analysis; Immunity; IFI44 gene

INTRODUCTION

As one of the first domesticated herbivores (Chessa et al.,2009), sheep (Ovisaries) remain an essential source of meat,wool, and milk for humans (Jiang et al., 2014). Sheep were originally domesticated 8 000-11 000 years ago from wild sheep species found in the Fertile Crescent (Diamond, 2002;Zeder, 2008). After domestication, sheep distribution expanded to meet the needs of different human populations,thereby adapting to distinct climatic environments and forming varieties with diverse phenotypes (Scher, 2000). Thus,understanding the genetic differences and diversity of species and varieties is a focus of animal breeding research.

Publication of the sheep reference genome provides an opportunity to determine the genetic mechanisms involved in artificial and natural selection and the phenotypic differentiation of domesticated sheep from their wild ancestors. For example, Zhao et al. (2017) used population single-nucleotide polymorphisms (SNPs) to identify candidate genes associated with high fertility, coat color, tail type, and horn size and type in sheep. Based on pooled whole-genome resequencing data, Wang et al. (2019) identified several vision-associated genes with functional loci in Chinese indigenous sheep breeds. In addition, studies have identified candidate genes related to pigmentation, nervous system,sensory perception, litter size, tail fat deposition, immunity,wool fineness, and climatic adaptation using genome-wide analysis of selection signatures (Alberto et al., 2018; Cao et al., 2021; Chen et al., 2021; Hu et al., 2019; Li et al., 2020; Lv et al., 2022; Yang et al., 2016). Although previous research has focused on sheep domestication, few studies have explored the function and population genetic effects of candidate genes associated with phenotypic variation during domestication. Furthermore, only a handful of candidate loci have been identified in sheep compared to those associated with the domestication of dogs (Axelsson et al., 2013), pigs (Li et al., 2017), chicken (Wang et al., 2020), and cattle (Chen et al., 2018).

Here, we generated whole-genome sequencing data for 100 samples across 10 domesticated sheep breeds, combined with publicly available whole-genome resequencing data from 208 individuals representing 20 domesticated sheep breeds and 17 wild sheep (Asiatic mouflon,O.orientalis), to characterize population genetic structure and genomic diversity at the genome-wide level and to elucidate the genome-wide genetic mechanisms involved in the wild-todomesticated process. Population genetic effects of important candidate genes were validated using phenotypic data of a Hu sheep population. Overall, we aimed to identify important genomic regions and molecular markers selected in sheep during domestication and improvement.

MATERIALS AND METHODS

Ethics statement

All applicable international, national, and/or institutional guidelines for the care and use of animals were strictly followed. All animal collection protocols complied with the current laws of China. The collection of blood samples and experimental protocols were approved by the Animal Care Committee of Lanzhou University (Permit No. 2010-1), in compliance with the recommendations of the Regulations for the Administration of Affairs Concerning Experimental Animals of China.

Sample collection and sequencing

Using the jugular vein method, we collected 100 blood samples from 10 different sheep breeds, including East Friesian milk sheep (EF), Dorper sheep (DP), Texel sheep(TK), South African mutton merino sheep (NM), Black Suffolk sheep (BS), Australian white sheep (AW), Mongolian sheep(MG), Lanzhou large-tailed sheep (LL), Altay sheep (AL), and large-tailed Han sheep (LH) (Supplementary Table S1). For each breed, 10 blood samples were obtained from five healthy females and five healthy males. Genomic DNA was extracted from whole blood of each individual using an EasyPure Blood Genomic DNA Kit (TransGen Biotech, China) according to the manufacturer’s recommended protocols. The A260/280 ratio and agarose gel electrophoresis were used to assess DNA quality and integrity, respectively. Paired-end sequencing libraries for each individual were constructed (mean insert size of 500 bp) using the Illumina NovaSeq 6000 platform (Illumina,USA). In addition, the genomic data of 225 individuals were downloaded from the NCBI database, including 34 African indigenous sheep, 61 improved sheep breeds, 113 Chinese indigenous sheep, and 17 Asiatic mouflons representing their wild ancestor (Supplementary Table S2).

Sequencing data processing

The raw sequencing data downloaded from the NCBI Sequence Read Archive (SRA) were converted to fastq files using SRAToolkit (v2.9.2) (Kodama et al., 2012). All fastq files were then filtered using Trimmomatic (v0.36) to obtain clean reads. The following filter criteria were used to remove adapters and low-quality bases: reads with ＞10% unknown nucleotides (N), reads with ＞50% low-quality (Q-value＜5)bases, and reads with ＞10 nucleotides aligned to the adaptor sequence with up to two mismatches. The resulting clean reads were mapped to the sheep Oar_v4.0 reference genome using BWA (Burrows-Wheeler Aligner) software (v0.7.8) with default parameters (Li & Durbin, 2009). Mapping files were converted to BAM files and sorted using SAMtools (v1.12) (Li et al., 2009). Picard software (v2.26.2) was used to mark potential polymerase chain reaction (PCR) duplicates for subsequent variant calling. SNP calling was performed using the Bayesian method in the GATK package (v3.4.0). VCFtools(v0.1.14) was used to correct the GATK results, with highquality SNPs retained for further analysis based on the following criteria: (1) mean coverage depth ≥5, (2) missing rates ＜10% and minor allele frequencies (MAF) ≥5%, and (3)root mean squared (RMS) mapping quality ≥20. After filtering,all SNPs were functionally annotated based on the Oar_v.4.0 sheep reference genome using ANNOVAR (v2013-05-20)(Wang et al., 2010). From the genome annotation, the SNPs were classified as variations in intronic regions, upstream and downstream regions, splicing sites, and exonic regions(synonymous or non-synonymous SNPs), while mutations causing stop gain and stop loss were grouped as nonsynonymous SNPs.

Population genetic structure analysis

To clarify the genetic relationships of domesticated sheep from a genome-wide perspective, a neighbor-joining (NJ) tree was constructed for the 325 sheep based on the matrix of pairwise genetic distances from the autosomal SNP data using TreeBeST v.1.9.2 (Vilella et al., 2009) and visualized using FigTree v.1.4.3. Population genetic structure was inferred using ADMIXTURE (v1.3.0) with default parameters(Alexander et al., 2009); the number of predefined ancestral clusters ranged fromk=2 tok=15. Principal component analysis (PCA) was performed for the 325 individuals using GCTA software (v1.26.0), with the first and second principal components displayed using the R software (v3.5.2) (Yang et al., 2011).

Analysis of genome-wide selective sweep regions

We analyzed genome-wide selective sweeps during domestication and improvement based on fixation index (FST)analysis using VCFtools. TheFSTvalues were calculated using the sliding window approach, with 150 kb windows and 75 kb sliding steps according to previous study (Wang et al.,2019). For domestication analysis, we combined the 308 domesticated sheep (African indigenous breeds, improved breeds, and Chinese indigenous breeds) into a group and compared them with Asiatic mouflons (representing their wild ancestor). The selective sweeps of 30 domesticated sheep breeds were also detected by comparison with wild sheep. To establish the domestication and improvement process (i.e.,wild sheep to indigenous breeds to improved breeds), we tested genomic selection signals for wild, indigenous, and improved sheep breeds. The formula used to calculateFSTfollowed previous study (Wang et al., 2019). Windows with the top 5‰ ofFSTvalues were defined as the selective regions,and genes in the selective regions were identified as candidate genes.

Functional enrichment analysis of candidate genes

Kyoto Encyclopedia of Genes and Genomes (KEGG)enrichment pathways and Human Phenotype Ontology (HPO)terms were analyzed to explore the most relevant functions of the protein-coding genes of selective regions using KOBAS-i(Bu et al., 2021) and g:Profiler, respectively.P＜0.05 was used as the threshold for significantly enriched pathways and functions.

Genotyping of candidate SNP loci

To verify that the candidate genes selected from the 30 domesticated sheep breeds during the domestication process were associated with immune response, blood DNA from 904 individuals was extracted from our previously established Hu sheep population with accurate immune trait data records(Zhang et al., 2022). The candidate gene loci were genotyped using competitive allele-specific fluorescence resonance energy transfer (FRET)-based PCR (KBioscience competitive allele-specific PCR amplification of target sequences and endpoint fluorescence genotyping (KASPar™)) assays (LGC Genomics, UK) according to a previously published method(Zhang et al., 2021b). The primer pairs for SNPs designed for KASPar genotyping are listed in Supplementary Table S3.

Statistical analysis

A mixed linear model (MLM) in the lem4 package in R was used to test the SNP effect on hematological parameters. The animal model was defined as:

where Yijklis the observed value of the hematological parameter indices, μ is the population mean, Genotypeiis the effect of each genotype at the ithSNP locus, Batchlis the batch effect, Genotypeiand Batchlare fixed effects, and Fatherj, Motherk, and εijklare random residual effects.P＜0.05 was considered as the criterion for statistical significance.

RESULTS

High-density sheep genomic variation map

In total, 100 domesticated sheep were collected for wholegenome resequencing, which yielded 16.97 billion clean reads, with a mean depth of 9.57× and an average genome coverage of 96.6% (Supplementary Table S1). The new genomic sequence data were combined with publicly available genomic data of 225 individuals from 21 breeds, for a total of 325 individuals (Figure 1A; Supplementary Table S2). After variant calling and filtration, we identified 39 331 475 highquality SNPs for subsequent population analyses(Supplementary Table S4).

Population genetic structure

To understand the genetic relationships and divergences between all domesticated and wild sheep, we determined the population-based proportions of allele frequency changes for each geographic group of breeds compared with Asiatic mouflons (representing their wild ancestor). When the allele frequency changed by less than 25%, the proportions of African indigenous, Chinese indigenous, and improved sheep breeds were 80.17%, 78.71%, and 76.51%, respectively(Figure 1A). Data from 30 domesticated sheep breeds were also compared with data from wild sheep, which showed consistent results as the different groups mentioned above(Figure 1C; Supplementary Table S5). We also observed higher genetic differentiation (FST) between the Asiatic mouflons and improved breeds than between the Asiatic mouflons and indigenous (African and Chinese) breeds(Figure 1D). These results indicated that the differences between African indigenous breeds and wild sheep were the smallest, followed by Chinese indigenous breeds and improved breeds.

Figure 1 Genomic diversity and population genetics of wild and domestic sheep

Using Asiatic mouflons as an outgroup, we characterized the genetic relationships among all individuals using the NJ tree, population structure analysis, and PCA. The NJ results suggested that the 325 individuals could be classified into four main groups (Figure 1B), representing wild sheep, African indigenous breeds, improved breeds, and Chinese indigenous breeds, respectively. Population structure analysis yielded similar results (Figure 1E, F). Based on PCA, the Asian mouflons were clustered together, while domesticated sheep were classified into African indigenous breeds, Chinese indigenous breeds, and improved breeds. The Chinese indigenous breeds were further divided into three groups,which showed strong geographical distribution characteristics.These results are consistent with the geographical and morphological classification of Chinese indigenous sheep as Mongolian, Kazakh, and Tibetan sheep (Figure 1G-I).

Selection signals during domestication in sheep

To explore the genome-wide selection signatures influenced by domestication, we combined the 308 domesticated sheep(African indigenous breeds, improved breeds, and Chinese indigenous breeds) as a group and compared them with Asiatic mouflons (representing their wild ancestor). The genome-wideFSTvalue was calculated between the domesticated and wild sheep based on a 150 kb sliding window and 75 kb shift across the genome along the autosomes and sex chromosome. In total, 87 and two putatively selected genomic regions were identified on the autosomes and sex chromosome, respectively, with the top 5‰ of globalFSTvalues (spanning 100.915 Mb) accounting for 3.901% of the complete genome and harboring 328 genes(Figure 2A; Supplementary Table S6). We then carried out KEGG pathway and HPO category enrichment analyses of these genes influenced by domestication. KEGG analysis identified 32 significantly enriched pathways, 13 of which were associated with immune response, six of which were related to metabolic processes, and one of which was associated with longevity regulating pathway (Supplementary Table S7). HPO analysis of the candidate genes identified 23 significantly enriched HPO terms, including terms associated with visual function, such as iris coloboma, abnormal iris morphology, and abnormality of the eye (Supplementary Table S8).

Figure 2 Genome-wide detection and annotations of regions selected during domestication

To investigate how many selective regions were shared among breeds,FSTvalues were calculated for the following 30 comparisons: five African indigenous breeds, 15 Chinese indigenous breeds, and 10 improved breeds versus Asian mouflons, respectively (Figure 2B). In selective sweep analysis, windows with the top 5‰ ofFSTvalues for the autosomes and sex chromosome were defined as selective windows. Results identified 25 windows that overlapped in more than 25 domestic sheep breeds. Among them, two regions on chromosomes 1 and 13 (Chr1: 54.3-54.525 Mb and Chr13: 50.325-50.55 Mb) were identified from the above 30 pairwise comparisons (Table 1). These two regions contained genes associated with the immune response (IFI44(encoding interferon induced protein 44 andIFI44L(encoding interferon induced protein 44 like) and visual function (PANK2(encoding pantothenate kinase 2) andRNF24(encoding ring finger protein 24)). Comparisons among the genomic regions adjacent to the four genes (IFI44,IFI44L,PANK2, andRNF24)revealed high genetic differentiation (FST) between the Asiatic mouflons and different groups of domestic sheep (Figure 2C).Furthermore, genotype pattern analysis of SNP loci revealed that the genotypes were significantly different between the Asiatic mouflons and domestic sheep (Figure 2D, E),suggesting the occurrence of an intensive selective sweep in the two regions. These results indicated that genes related to immunity and sensory ability were strongly selected during early domestication.

Selective imprints during domestication and breeding

To further determine the significance of genomic divergence in sheep domestication and breeding history, we analyzed and compared the genomes of wild and indigenous groups, as well as indigenous and improved groups, to identify potential selection imprints that occurred during this process(Figure 3A). The globalFSTvalues for each comparison were calculated using 150 kb sliding windows with 75 kb steps across the genome. The top 5‰ of windows were defined as selective windows. After the adjacent selective windows were merged, a total of 90 domestication regions and 72 improvement regions were identified, comprising 327 and 326 candidate genes, respectively (Figure 3B, F; Supplementary Tables S9, S10). TheFSIP2(encoding fibrous sheathinteracting protein 2) gene, which is associated with reproduction traits on chromosome 2, was identified during domestication. The genotype pattern ofFSIP2showed significant differences between the Asiatic mouflons and indigenous sheep, and linkage disequilibrium analysis indicated strong linkage of SNPs in this region (Figure 3C-E).To explore the strongly selected regions of chromosomes 10 and 19 during the breeding process, we analyzed genomic architecture by calculating allele frequency at nonsynonymous SNPs in two genes (FGF9(encoding fibroblast growth factor 9) andMITF(encoding melanocyte inducing transcription factor)) (Figure 3G-I). Two variant alleles (c.35570188 A＞G, c. 35570233 T＞C) located in the exon region ofFGF9showed higher frequencies in the improved breeds,but lower frequencies in the indigenous breeds. Similarly, one variant allele (c. 31605530 C＞T) located in the second exon ofMITFdiffered between white sheep and non-white sheep.Annotation of candidate genes for domestication indicated that they were mainly related to immune processes (RAP1 signaling pathway, lysosome, and NOD-like receptor signaling pathway) and the thyroid hormone signaling pathway(Figure 4A; Supplementary Table S11). Furthermore, genes selected during breeding and improvement were mainly association with melanoma, pathways in cancer, notch signaling pathway, and starch and sucrose metabolism. The most significantly enriched pathway was melanoma(Figure 4B; Supplementary Table S12). Notably, seven candidate genomic regions overlapped in both domestication and improvement (Supplementary Table S13), suggesting that several important domestication loci may have undergone a second round of artificial selection for continued improvement of vital economic traits. Furthermore, candidate genes of the overlapping genomic regions were significantly enriched in cholinergic synapse, glycosylphosphatidylinositol (GPI)-anchor biosynthesis, and phototransduction pathways (Figure 4C;Supplementary Table S14).

Table 1 Information on overlapping genomic regions on autosomes and X chromosome of more than 25 sheep breeds

Figure 3 Genome-wide distribution of selective sweeps during domestication and breeding

SNP genotyping and association analysis

To further explore the effects of key genes in significant candidate regions on immunity traits, a novel mutation (c.54413058 C＞G) in theIFI44gene was genotyped using KASPar technology (Supplementary Figure S1, blue, green,and red dots indicate three different genotypes, respectively).We performed association analysis of hematological parameter indices based on our previously established Hu sheep population with accurate phenotypic data records(Zhang et al., 2022). Results indicated that theIFI44c.54413058 C＞G polymorphism was significantly associated with white blood cells (WBCs), neutrophils, and monocytes(P＜0.05). Furthermore, the WBC phenotype value in Hu lambs with the CC and GC genotypes was significantly higher than in those with the GG genotype, while the difference between CC and GC lambs was not significant (Table 2). These results indicated that CC was the dominant genotype associated with WBC count.

DISCUSSION

Understanding the molecular basis underlying genetic variation and phenotypic changes during animal domestication and subsequent selection could contribute to animal breeding and help determine how phenotypes are influenced by genotypic changes. Whole-genome sequencing has been widely applied to reveal the genomic variants under selection in domestic animals (e.g., horse, pig, cattle, and sheep) (Chen et al., 2021; Daetwyler et al., 2014; Liu et al., 2019b; Xu et al.,2020). Here, we collected whole-genome data from 325 sheep for comprehensive population genetic and selective signaling analysis at the genome level.

Elucidating population structure of sheep at genome level

Analysis of population genetic structure can enhance our understanding of the domestication process of certain species or breeds. Genome-level analysis of allele frequency changes showed that differences were the smallest between African indigenous breeds and wild sheep, followed by differences between Chinese indigenous breeds and improved breeds(Figure 1A, C). Based on population structure analysis, the 325 individuals could be grouped into two clusters, i.e.,domestic sheep and wild sheep, suggesting that all domestic sheep originated from a single domestication event(Figure 1G). Domestic sheep were further clustered into three groups, i.e., African indigenous, Chinese indigenous, and improved sheep breeds (Figure 1H). The Chinese indigenous breeds were also grouped into three clusters, i.e., Mongolian,Kazakh, and Tibetan sheep, which showed strong geographical distribution and characteristic morphological patterns. These results are in accordance with previous study using the Illumina Ovine SNP 50 K Bead Chip assay (Wei et al., 2015).

Figure 4 Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of candidate genes during domestication and breeding

Table 2 Association results between genotypes of ovine IFI44 gene and hematological parameters

Genomic regions related to domestication and improvement in sheep

First, for selection signature analysis, we focused on the whole domestication process in domestic sheep. Based on genome-wide comparisons of wild sheep and 308 domesticated sheep from 30 breeds, we identified several genes associated with immune response, metabolic processes, longevity-regulating pathway, and visual function(Supplementary Tables S7, S8). Vision plays a vital role in the survival and evolution of animals, such as predator avoidance,mate selection, and foraging (Yokoyama, 2002). Second, we detected selective sweeps in the 30 domesticated sheep breeds in comparison to wild sheep. Results identified two overlapping genomic regions for the 30 pairwise comparisons,which contained genes associated with immune response(IFI44andIFI44L) and visual function (PANK2andRNF24).TheIFI44andIFI44Lgenes belong to the type I interferoninducible gene family and are located on the same chromosome. These two genes play important roles in regulating autoimmune disorders, inflammation, and immune response, and inhibit respiratory syncytial virus infection(Busse et al., 2020). Previous studies have also reported that variants of thePANK2 gene are related to pantothenate kinase-associated neurodegeneration (e.g., optic atrophy),and genetic variants ofRNF24andPANK2are associated with optic disc morphology (Axenovich et al., 2011). Thus,these two genes may be associated with the evolution of vision during sheep domestication. Vision plays a vital role in animal survival, and many studies have demonstrated that visual acuity is weaker in domestic animals (e.g., chickens,dogs, and ducks) compared to their wild ancestors(Henderson et al., 2000; Peichl, 1992; Wang et al., 2016). We hypothesized that sensory ability and immunity may have been important targets of selection during domestication.

We also analyzed the potential selection signatures of wildto-indigenous sheep domestication and identified candidate genes associated with immune response (IFI44L,IFI44),reproductive traits (SPAG16(encoding sperm associated antigen 16) andFSIP2), and visual function (PDE6B(encoding phosphodiesterase 6B),FOXC1(encoding forkhead box C1), andGMDS(GDP-mannose 4,6-dehydratase)) during domestication.FSIP2is a protein-coding gene and plays an important role in spermatogenesis (Fang et al., 2021).Homozygous loss-of-function mutations inFSIP2can lead to male infertility (Liu et al., 2019a). In the present study,FSIP2was located in a significant selective region, and the genotype pattern differed between the wild and indigenous sheep breeds, suggesting thatFSIP2may play an important role in sheep fertility.

In the process of breeding and improvement, we identified several candidate genes associated with pigmentation,including several identified in previous research (e.g.,ASIP(encoding agouti signaling protein) andMITF) (Li et al., 2014).Mutations inMITFcan decrease pigmentation in dogs(Karlsson et al., 2007), pigs (Chen et al., 2016), ducks (Zhou et al., 2018), and quails (Minvielle et al., 2010).MITFalso encodes a protein in the melanoma pathway. We identified a SNP inMITFthat differed between white and non-white sheep. Moreover, two non-synonymous mutations ofFGF9were found at higher frequencies in the improved breeds than in the indigenous breeds.FGF9is a member of the fibroblast growth factor family, and is involved in multiple biological processes, such as cartilage development, cell growth, and embryonic development (Zhang et al., 2021a, 2021c). These results imply that important economic traits (e.g., coat color,body size) are the preferred targets during breeding improvement. In addition, seven genomic regions were continuously selected at the two stages, indicating that some candidate loci may have undergone a second round of artificial selection for continued improvement of important economics traits.

Functional exploration of candidate SNPs in IFI44 gene

To further verify the function of theIFI44gene, anIFI44SNP(c. 54413058 C＞G) was genotyped and subjected to association analysis in a Hu sheep population. Results indicated that the SNP atIFI44c. 54 413 058 C＞G was significantly associated with WBCs, neutrophils, and monocytes (Table 2). WBCs play an important role in the immune system. A feature of the inflammatory response is an increase in WBC count following bacterial or viral infection,and neutrophils, lymphocytes, and monocytes also play critical roles in innate and adaptive immunity (Parish, 2006). Thus, we concluded that the polymorphic sites inIFI44may be important molecular markers of overcoming immune deficiency in animals. Nevertheless, the immune function of theIFI44gene needs to be verified at the cellular and protein levels.

DATA AVAILABILITY

Raw data were deposited in the National Center for Biotechnology Information database under BioProjectID PRJNA795904 and PRJNA777695, in the Genome Sequence Archive under Accession No. CRA007173, and in the Science Data Bank under DOI: 10.577 60/sciencedb.01846.

SUPPLEMENTARY DATA

Supplementary data to this article can be found online.

COMPETING INTERESTS

The authors declare that they have no competing interests.

AUTHORS’ CONTRIBUTIONS

W.M.W. and F.D.L. designed the project. D.Y.Z., X.X.Z.,F.D.L., and W.M.W. contributed to blood sample collection.D.Y.Z., L.F.Y., W.M.W., X.L.L., Y.K.Z., and Y.Z. analyzed the data. D.Y.Z., X.X.Z., L.M.Z., J.H.W., D.X., J.B.C., X.B.Y.,W.X.L., C.C.L., and B.B.Z. contributed to data collection for the validation of experimental populations. D.Y.Z., X.X.Z.,X.L.L., and X.B.Y. participated in DNA extraction. D.Y.Z. wrote the paper. W.M.W., X.X.Z., and D.Y.Z. reviewed and edited the manuscript. All authors read and approved the final version of the manuscript.

ACKNOWLEDGMENTS

We would like to thank the staff at our laboratory for their ongoing assistance.