Detection of quantitative trait loci affecting haematological traits in swine via genome scanning
BMC Genetics volume 11, Article number: 56 (2010)
Haematological traits, which consist of mainly three components: leukocyte traits, erythrocyte traits and platelet traits, play extremely important role in animal immune function and disease resistance. But knowledge of the genetic background controlling variability of these traits is very limited, especially in swine.
In the present study, 18 haematological traits (7 leukocyte traits, 7 erythrocyte traits and 4 platelet traits) were measured in a pig resource population consisting of 368 purebred piglets of three breeds (Landrace, Large White and Songliao Black Pig), after inoculation with the swine fever vaccine when the pigs were 21 days old. A whole-genome scan of QTL for these traits was performed using 206 microsatellite markers covering all 18 autosomes and the X chromosome. Using variance component analysis based on a linear mixed model and the false discovery rate (FDR) test, 35 QTL with FDR < 0.10 were identified: 3 for the leukocyte traits, 28 for the erythrocyte traits, and 4 for the platelet traits. Of the 35 QTL, 25 were significant at FDR < 0.05 level, including 9 significant at FDR < 0.01 level.
Very few QTL were previously identified for hematological traits of pigs and never in purebred populations. Most of the QTL detected here, in particular the QTL for the platelet traits, have not been reported before. Our results lay important foundation for identifying the causal genes underlying the hematological trait variations in pigs.
In pork-producing community, various infectious diseases caused by viral or bacterial pathogens badly restrict the efficiency of pork industries, and also seriously affect human health in the whole world. To most infectious diseases, the resistance of an individual is resulted from both innate immunity and acquired immunity. The capacity of innate immunity or acquired immunity is more or less controlled by genes [1–4]. The resistance to some diseases is due to a single gene [5, 6], but for most diseases it is due to multiple genes , among which a few genes (QTL or major genes) have relative large effects and thus play important role in disease resistance. Therefore, detection of these genes or QTL is important for disease-resistant breeding.
Haematological traits, which consist of three components: leukocyte traits, erythrocyte traits, and platelet traits, are important indicators of immune capacity of animals [8–10]. The influence of genetics on baseline peripheral blood cell parameters has been firmly established. In human, the heritability estimates for red blood cell, white blood cell, and platelet numbers derived from a twin study  were 0.42, 0.62, and 0.57, respectively, indicating a large genetic component underlies these traits. Study in baboons  confirmed that genetic influences account for a significant portion of the variance in these blood parameters as well as in mean corpuscular volume (MCV) and mean platelet volume (MPV). In swine, significant differences in a panel of haematological traits between Meishan and Large White pigs were observed by Clapperton et al. , which may have implications in the resistance to infection by a broad range of pathogens and subsequent disease effects in these breeds. However, knowledge of the genetic background controlling the variability of haematological traits, is very limited, especially in swine.
With the continuously increasing of DNA molecular markers and deep development of pig genome research, a large number of QTL or major genes controlling quantitative traits in pigs have been found and located in the relevant regions of chromosomes [13–15]. However, compared to the production traits, QTL mapping for haematological traits has been carried out only in a few research groups. Edfors-Lilja et al. [16, 17] first reported QTL on SSC1 and 8 for stress-mediated alterations in some leukocyte traits in swine, which were further confirmed later in another mapping study by the same group . Reiner et al.  reported nine genome-wide significant and 29 putative QTL distributed on 15 chromosomes controlling baseline leukocyte traits and leukocyte traits at distinct stages of the disease model. QTL for erythroid traits were also reported by the same group . 43 QTL controlling baseline erythroid traits and erythroid traits at distinct stages of the disease model were identified on 16 chromosomes, of which twelve were significant at the genome-wide level. In a very recent study , a total of 101 QTL, including 46 genome-wide significant QTL, regulating baseline erythroid traits at different age stages were found on all pig chromosomes except for SSC15 and SSC18. For platelet traits, there are no reports of QTL mapping studies in swine up to now.
The immune response to "stress" is important indicator of disease resistance capacity of animals. In the present study, we report the results of QTL mapping for the stress induced immune response of pigs which were challenged with swine fever vaccine. Unlike the previous studies mentioned above, where F2 animals derived from cross of two breeds was used as mapping population, we used purebred animals with half- and full-sib families as mapping population. 18 haematological traits, including seven leukocyte traits, seven erythrocyte traits, and four platelet traits, were measured after they were inoculated with the swine fever vaccine.
The means and standard errors of the phenotypic values of the 18 haematological traits in the mapping population are presented in Table 1. These results are generally consistent with previous results [19–21]. Significant differences (P < 0.05) between measurements before (on day 20) and after (on day 35) inoculation with the swine fever vaccine were observed for all traits except GR%, LY% and MO%. Differences between the three breeds are significant for most of the leukocyte traits and two of the four platelet traits, but not significant for all erythrocyte traits except MCHC. No significant differences between sexes were observed for all traits except RBC and MCHC.
QTL for the 18 haematological traits
The QTL detected with FDR < 0.10 for the 18 haematological traits are presented in Table 2. For the leukocyte traits, only three QTL, including one with FDR < 0.05, were detected for GR, MO% and WBC, respectively. For the erythrocyte traits, 28 QTL, including 21 with FDR < 0.05, were detected: ten for MCHC, seven for RDW, five for MCV, three for MCH, two for HGB, and one for HCT. For the platelet traits, four QTL, including three with FDR < 0.05, were detected: two for PCT and two for PDW. The positions of all of the detected QTL are shown in Fig. 1.
So far, there are only limited reports on QTL mapping for haematological traits in pigs [16–21]. All these studies used F2 animals as mapping population and only one category of haematological traits (leukocyte traits or erythrocyte traits) were involved. In the present study, we used purebred animals as mapping population and all three categories of haematological traits were analyzed. There is large inconsistency between the QTL for leukocyte and erythrocyte traits reported in the present and the previous studies. For the leukocyte traits, no QTL are in agreement with those reported in the previous studies. In particular, for most traits, the QTL found here and in the previous studies were mapped to different chromosomes. For example, for the trait WBC, we mapped two significant QTL on SSCX, but the QTL for this trait were found on SSC1, SSC8, SSC10, and SSC12 in the previous studies [13–16]. For the erythrocyte traits, only a few QTL are in agreement with those reported in their studies. For HGB and MCHC, we mapped a QTL at 31 cM on SSC2 and a QTL at 19 cM on SSC16, respectively, which overlap those reported by Reiner et al. . For RDW, we mapped a QTL at 20 cM on SSC1, which overlaps that reported by Zou et al. . Large inconsistency also exists between the results from Reiner et al.  and those from Zou et al. . Both studies involved in mapping QTL for erythrocyte traits in swine and employed the regression approach and a F2 design. 43 and 101 QTL were found by Reiner et al.  and Zou et al. , respectively. However, only four are in accordance with each other. A major possible reason for these inconsistencies is that the measurements for the traits are different in different studies. In the study of Edfors-Lilja et al. , the leukocyte traits were measured immediately before and after the pigs were transported from the farrowing pen to a finishing unit and mixed with non-littermates when they were about three months old. In the study of Zou et al. , the erythrocyte traits were measured when the pigs were 18, 46, and 240 days old without any stress induction. In the study of Reiner et al. [19, 20], the leukocyte and erythrocyte traits were measured on the day (when the pigs were approximately 100 days old) of and shortly before infection with 50 000 sporocysts of Sarcocytis miescheriana of and on days 14, 28, and 34 after infection. While in the present study, all haematological traits were measured one day before and 14 days after the inoculation with the swine fever vaccine when the pigs were 21 days old. The immune system of pigs is very sensitive to different external stresses. The genetic mechanism of immune responses to different external stress and in different physiological stage is very complex and currently unknown. Therefore, the phenotypes for the same haematological parameters measured in different studies may represent different traits, although they may be highly correlated. Another possible important reason for the discrepancies is that different mapping populations were used in different studies, in which the QTL segregation status may be different. In all cited studies, QTL were detected in crossbred populations while they were detected in purebred populations in the present study.
QTL mapping for platelet traits in swine has not been reported yet. However, QTL for platelet traits had been studied in mice [22, 23]. Among the QTL found in the two studies, a QTL for PLT was mapped to 52 Mb  and 51 Mb  on MMU7, respectively. This region is syntenic with the region on SSC2, where we found a putative QTL (P < 0.05, data not shown) for PLT in swine. In addition, Chueng et al.  mapped a QTL for PLT in mice to 48 Mb on MMU17, which is syntenic with the region of 71 cM on SSC7, where we also found a putative QTL (P < 0.05, data not shown) for PLT in swine.
Similar to the results of previous studies [20, 21], some QTL located at the same or nearby positions were found to have pleiotropic effects on two to three traits in the present study. These include the QTL located at 38 cM on SSC2 (affecting PDW and PCT), the QTL located at 124 cM on SSC9 (affecting MCH and MCHC), the QTL located at 8 and 9 cM on SSC12 (affecting HGB and HCT), and the QTL located at 51 cM on SSC14 (affecting MCH, MCHC, and RDW). These pleiotropic QTL contribute largely to the correlations between traits. For example, we can imagine a high genetic correlation between traits MCH and MCHC, and indeed, two QTL, one at 124 cM on SSC9 and the other at 51 cM on SSC14, were found having significant effects on both traits.
The candidate genes related to leucocyte traits have been discussed in detail by Wattrang et al. . However, none of these genes are located in the QTL regions identified in this study. For the erythrocyte traits, three candidate genes were reported, including the porcine stem cell growth factor receptor (KIT) gene , the erythropoietin (EPO) gene , and the erythropoietin receptor (EPOR) gene . The KIT gene is located at p12-p21on SSC8 and is essential for normal hematopoiesis, melanogenesis, and germ cell development . The EPO gene is located at p16-p15 on SSC3 and plays important role in controlling the level of circulating erythrocyte mass by promoting erythroid differentiation and initiating haemoglobin synthesis . The EPOR gene is located at q1.2-q2.1 on SSC2 and is critical for erythroid proliferation and differentiation . In the present study, no QTL was found around the position of the KIT gene, a putative QTL (P < 0.01, data not shown) affecting HGB was mapped between markers SW256 and S0141 on SSC2 which is close to that of the EPOR gene, and a QTL affecting MCV was mapped between markers SW2618 and SW836 on SSC3 which is near the position of the EPO gene.
Very few QTL were previously identified for hematological traits of pigs and never in purebred populations. In this study, 35 QTL with FDR < 0.10 were identified for 18 haematological traits: 3 for the leukocyte traits, 28 for the erythrocyte traits, and 4 for the platelet traits. Most of the QTL detected here, in particular the QTL for the platelet traits, have not been reported before. Our results lay important foundation for identifying the causal genes underlying the hematological trait variations in pigs.
The mapping population consisted of 368 purebred piglets of three breeds, including 5 Landrace boar families (15 sows and 87 piglets), 7 Large White boar families (33 sows and 190 piglets), and 4 Songliao Black Pig boar families (15 sows and 90 piglets), respectively. In each litter, 5~7 piglets were randomly selected. Songliao Black Pig is a breed derived from cross of Landrace, Duroc and Ming Pig (a local Chinese pig breed), so it has close genetic link with Landrace. All pigs were raised under standard indoor conditions at the experimental farm of the Institute of Animal Sciences, Chinese Academy of Agricultural Sciences.
Measurement of phenotypic traits
The piglets were inoculated with the swine fever vaccine on day 21 after birth. Blood samples were collected from each piglet one day before the inoculation (day 20) and two weeks after the inoculation (day 35), respectively. The samples were directly injected into eppendorf tubes containing 60 μl of 20% EDTA in phosphate-buffered saline (PBS). 18 hematologic traits (Table 1) were measured immediately after sample collection using a TEK-II mini automatic hemocyte analysor (Jiangxi Tekang Science and Technology Company, China) with swine specific parameter configuration at Beijing Xiyuan Hospital, China.
Genomic DNA was extracted from ear tissues of all piglets and their parents (445 animals in total) using a standard phenol/chloroform method. Microsatellite markers were initially selected from the USDA-MARC porcine reference map http://www.marc.usda.gov/ and were then analyzed for their polymorphisms using DNA from the boars. A final set of 206 informative markers distributed on 19 chromosomes (18 autosomes and chromosome X) was genotyped for all individuals. The total length covered by these markers is 2261.7 centi-Morgans (cM) with an average distance of 11.6 cM between adjacent markers according to the MARC map. The number of markers and their mean polymorphic information content (PIC) on each chromosome are shown in Table 3. For the QTL analysis, the linkage maps of the selected markers were derived directly from the MARC map.
where y is a vector of the phenotypic values for one of the haematological traits measured on day 35, a is a vector of fixed effects including breed (three), sex (two) and sampling season (three), c is a vector of the values for that haematological trait measured on day 20, b is the regression coefficient, u is a vector of polygenic effects, v is a vector of QTL allelic effects, e is a vector of residuals, X, Z, and T are incidence matrices for a, u, and v, respectively, A is the additive genetic relationship matrix among all individuals, is the additive polygenetic variance, Q is the IBD probability matrix among QTL alleles, is the QTL allelic variance, I is an unit matrix, and is the residual variance.
The QTL was scanned from left to right at 1 cM intervals along each chromosome. Restricted maximum likelihood (REML) was used to estimate the three variance components in the model and the likelihood ratio (LR) was calculated as test statistic for each particular location on the chromosome
where LMAX | H0 and LMAX | HA are restricted maximum likelihood functions corresponding to the null hypothesis (there is no QTL) and the alternative hypothesis (there is a QTL), respectively.
It is generally regarded that the LR follows approximately a χ2 distribution with degrees of freedom between one and two . To be more conservative, we assumed here it follows a χ2 distribution with two degrees of freedom. In addition, in order to avoid the increase of false positive rate caused by multiple tests, the FDR (false discovery rate) control approach [32, 33] was adopted to determine the significance thresholds based on the P values from the χ2 distribution. The positions with FDR < 0.05 were considered harboring significant QTL. Let m be the total number of tests involved in the analysis and P1 ≤ P2 ≤⋯≤ P m be the ordered observed P-values for the m tests, the FDR for P i is calculated as
Since the QTL analysis was performed at 1 cM intervals along each chromosome, each genome scan is equivalent to 2,261 tests given the assumed genome length. As these tests are not independent, Weller et al  recommended as a cautious choice considering the number of markers genotyped. This number was 206 in the present study, leading to the total number of tests m = 18 (traits) × 206 (markers) = 3,708.
The LOD-drop off method  was used to estimate the confidence intervals of the QTL positions.
Since the variance component method requires normal distribution of the phenotypes, we tested the normality of all traits. It turned out that most of the traits were not normally distributed. The Box-Cox transformation was performed for the traits which severely depart from normal distribution before the QTL analysis was carried out.
Clapperton M, Bishop SC, Glass EJ: Innate immune traits differ between Meishan and Large White pigs. Vet Immunol Immunopathol. 2005, 104 (3-4): 131-144. 10.1016/j.vetimm.2004.10.009.
Vincent AL, Thacker BJ, Halbur PG, Rothschild MF, Thacker EL: In vitro susceptibility of macrophages to porcine reproductive and respiratory syndrome virus varies between genetically diverse lines of pigs. Viral Immunology. 2005, 18 (3): 506-512. 10.1089/vim.2005.18.506.
Vincent AL, Thacker BJ, Halbur PG, Rothschild MF, Thacker EL: An investigation of susceptibility to porcine reproductive and respiratory syndrome virus between two genetically diverse commercial lines of pigs. Journal of Animal Science. 2006, 84 (1): 49-57.
Wimmers K, Kumar KG, Schellander K, Ponsuksili S: Porcine IL 12A and IL 12B gene mapping, variation and evidence of association with lytic complement and blood leucocyte proliferation traits. International Journal of Immunogenetics. 2008, 35 (1): 75-85.
Fujii J, Otsu K, Zorzato F, de Leon S, Khanna VK, Weiler JE, O'Brien PJ, MacLennan DH: Identification of a mutation in porcine ryanodine receptor associated with malignant hyperthermia. Science. 1991, 253 (5018): 448-451. 10.1126/science.1862346.
Jansen JH, Høgåsen K, Grøndahl AM: Porcine membranoproliferative glomerulonephritis type II: an autosomal recessive deficiency of factor H. the Veterinary Record. 1995, 137 (10): 240-244.
Kadarmideen HN: Biochemical, ECF18R, and RYR1 gene polymorphisms and their associations with osteochondral diseases and production traits in pigs. Biochemical Genetics. 2008, 46 (1-2): 41-53. 10.1007/s10528-007-9127-5.
Ford I, Douglas CW: The role of platelets in infective endocarditis. Platelets. 1997, 8 (5): 285-294. 10.1080/09537109777159.
Oddgeirsson O, Simpson SP, Morgan ALG, Ross DS, Spooner RL: Relationship between the bovine major histocompatibility complex (BoLA), erythrocyte markers and susceptibility to mastitis in Icelandic cattle. Animal Genetics. 1988, 19 (1): 11-16.
Müller M, Brem G: Disease resistance in farm animals. Experientia. 1991, 47 (9): 923-934. 10.1007/BF01929883.
Garner C, Tatu T, Reittie JE, Littlewood T, Darley J, Cervino S, Farrall M, Kelly P, Spector TD, Thein L: Genetic influences on F cells and other hematologic variables: a twin heritability study. Blood. 2000, 95 (1): 342-346.
Mahaney MC, Brugnara C, Lease LR, Platt OS: Genetic influences on peripheral blood counts: a study in baboons. Blood. 2005, 106 (4): 1210-1214. 10.1182/blood-2004-12-4863.
Hu ZL, Dracheva S, Jang W, Maglott D, Bastiaansen J, Rothschild MF, Reecy JM: A QTL resource and comparison tool for pigs: PigQTLDB. Mammalian Genome. 2005, 16 (10): 792-800. 10.1007/s00335-005-0060-9.
Hu ZL, Fritz ER, Reecy JM: AnimalQTLdb: a livestock QTL database tool set for positional QTL information mining and beyond. Nucleic Acids Research. 2007, 35 (Suppl 1): D604-609. 10.1093/nar/gkl946.
Rothschild MF, Hu ZL, Jiang ZH: Advances in QTL Mapping in Pigs. International Journal of Biological Sciences. 2007, 3 (3): 192-197.
Edfors-Lilja I, Wattrang E, Marklund L, Moller M, Andersson-Eklund L, Andersson L, Fossum C: Mapping quantitative trait loci for immune capacity in the pig. The Journal of Immunology. 1998, 161 (2): 829-835.
Edfors-Lilja I, Wattrang E, Anderesson L, Fossum C: Mapping quantitative trait loci for stress induced alterations in porcine leukocyte numbers and functions. Animal Genetics. 2000, 31 (3): 186-193. 10.1046/j.1365-2052.2000.00628.x.
Wattrang E, Almgvist M, Johansson A, Fossum C, Wallgren P, Pielberg G, Anderesson L, Edfors-Lilja I: Confirmation of QTL on porcine chromosomes and 8 influencing leukocyte numbers, haematological parameters and leukocyte function. Animal Genetics. 2005, 36 (4): 337-345. 10.1111/j.1365-2052.2005.01315.x.
Reiner G, Fischer R, Hepp S, Berge T, Köhler F, Willems H: Quantitative trait loci for white blood cell numbers in swine. Animal Genetics. 2008, 39 (2): 163-168. 10.1111/j.1365-2052.2008.01700.x.
Reiner G, Fischer R, Hepp S, Berge T, Köhler F, Willems H: Quantitative trait loci for red blood cell traits in swine. Animal Genetics. 2007, 38 (5): 447-452. 10.1111/j.1365-2052.2007.01629.x.
Zou ZZ, Ren J, Yan XM, Huang X, Yang SJ, Zhang ZY, Yang B, Li WB, Huang LS: Quantitative trait loci for porcine baseline erythroid traits at three growth ages in a White Duroc × Erhualian F2 resource population. Mammalian Genome. 2008, 19 (9): 640-646. 10.1007/s00335-008-9142-9.
Peters LL, Zhang W, Lambert AJ, Brugnara C, Churchill GA, Platt OS: Quantitative trait loci for baseline white blood cell count, platelet count and mean platelet volume. Mammalian Genome. 2005, 16 (10): 749-763. 10.1007/s00335-005-0063-6.
Cheung CC, Martin IC, Zenger KR, Donald JA, Thomson PC, Moran C, Buckey MF: Quantitative traits loci for steady platelet count in mice. Mammalian Genome. 2004, 15 (10): 784-797. 10.1007/s00335-004-2408-y.
Sakurai M, Zhou JH, Ohtaki M, Itoh T, Murakami Y, Yasue H: Assignment of c-KIT gene to swine chromosome 8p12-p21 by fluorescence in situ hybridization. Mammalian Genome. 1996, 7 (5): 397-10.1007/BF03035316.
Wan-sheng L, Ingrid H, Ingemar G, Bhanu PC: Mapping of the porcine erythropoietin gene to chromosome 3p16-pl5 and ordering of four related subclones by fiber-FISH and DNA-combing. Hereditas. 1998, 128 (1): 77-81. 10.1111/j.1601-5223.1998.00077.x.
Fahrenkrug SC, Campbell EM, Vallet JL, Rohrer GA: Physical assignment of the porcine erythropoietin receptor gene to SSC2. Animal Genetics. 2000, 31 (1): 69-70. 10.1111/j.1365-2052.2000.579-2.x.
Russell ES: Hereditary anemias of the mouse: a review for geneticists. Advances in Genetics. 1979, 20: 357-459. full_text.
Goldwasser E: Erythropoietin and the differentiation of red blood cells. Federation Proceedings. 1975, 34 (13): 2285-2292.
Grignola FE, Hoeschele I, Tier B: Mapping quantitative trait loci in outcross populations via residual maximum likelihood. I. Methodology. Genetics Selection Evolution. 1996, 28 (6): 479-490. 10.1186/1297-9686-28-6-479.
Grignola F, Hoeschele I, Zhang Q, Thaller G: Mapping quantitative trait loci in outcross populations via residual maximum likelihood.II. A simulation study. Genetics Selection Evolution. 1996, 28 (6): 491-504. 10.1186/1297-9686-28-6-491.
Zhang Q, Boichard D, Hoeschele I, Ernst C, Eggen A, Murkve B, Pfister-Genskow M, Witte LA, Grignola FE, Uimari P, Thaller G, Bishop MD: Mapping quantitative trait Loci for milk production and health of dairy cattle in a large outbred pedigree. Genetics. 1998, 149 (4): 1959-1973.
Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society. 1995, 57 (1): 289-300.
Weller JI, Song JZ, Heyen DW, Lewin HA, Ron M: A new approach to the problem of multiple comparisons in the genetic dissection of complex traits. Genetics. 1998, 150 (4): 1699-1706.
Lander ES, Botstein D: Mapping mendelian factors underlying quantitative traits using RFLP linkage maps. Genetics. 1989, 121 (1): 185-199.
This study was supported by the National Key Basic Research Program (Grant NO.2006CB102104) and the National Natural Science Foundation of China (Grant No. U0631005). We thank Dr. Lixian Wang (Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China) for help with sample collection. We also thank the peer reviewers for their valuable comments and suggestions.
YFG and XL are the major executive persons of all jobs of this study, including collection of the phenotypes, marker genotyping, statistical analysis, and drafting this manuscript. ZPW and JZ assisted in the statistical analysis. FH, YRL, SQC, CMQ, SL, XYN and XTQ assisted in phenotype collecting and marker genotyping. QZ planed and supervised the whole study. All authors read and approved the final manuscript.
Yuan-Fang Gong, Xin Lu contributed equally to this work.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Gong, YF., Lu, X., Wang, ZP. et al. Detection of quantitative trait loci affecting haematological traits in swine via genome scanning. BMC Genet 11, 56 (2010). https://doi.org/10.1186/1471-2156-11-56
- False Discovery Rate
- Mapping Population
- Mean Platelet Volume
- Mean Corpuscular Volume