Skip to main content

Genome-wide association studies for hematological traits in Chinese Sutai pigs

Abstract

Background

It has been shown that hematological traits are strongly associated with the metabolism and the immune system in domestic pig. However, little is known about the genetic architecture of hematological traits. To identify quantitative trait loci (QTL) controlling hematological traits, we performed single marker Genome-wide association studies (GWAS) and haplotype analysis for 15 hematological traits in 495 Chinese Sutai pigs.

Results

We identified 161 significant SNPs including 44 genome-wide significant SNPs associated with 11 hematological traits by single marker GWAS. Most of them were located on SSC2. Meanwhile, we detected 499 significant SNPs containing 154 genome-wide significant SNPs associated with 9 hematological traits by haplotype analysis. Most of the identified loci were located on SSC7 and SSC9.

Conclusions

We detected 4 SNPs with pleiotropic effects on SSC2 by single marker GWAS and (or) on SSC7 by haplotype analysis. Furthermore, through checking the gene functional annotations, positions and their expression variation, we finally selected 7 genes as potential candidates. Specially, we found that three genes (TRIM58, TRIM26 and TRIM21) of them originated from the same gene family and executed similar function of innate and adaptive immune. The findings will contribute to dissection the immune gene network, further identification of causative mutations underlying the identified QTLs and providing insights into the molecular basis of hematological trait in domestic pig.

Background

Hematological cells play essential roles in the immune responding process of disease resistance [1, 2]. Hematological cells are composed of three components, including leukocytes (white blood cells, WBCs), erythrocytes (red blood cells, RBCs) and platelets [3]. The major functions of leukocytes are innate and adaptive immunity and defending subject from pathogens [4, 5]. White blood cell count is a strong indicator of infectious and inflammatory diseases, such as leukaemia and lymphoma. Erythrocytes execute a range of functions such as transporting oxygen, carbon dioxide and killing pathogens [6, 7]. RBCs disorders indicate the increasing risk of anemia, polycythemia, hypertension and heart failure. Platelets play important roles in hemostasis, the initiation of wound repair and can be the strong effector cells of innate immune response [811]. Accompanying with reduction of platelet count, idiopathic thrombocytopenic purpura (ITP), often an idiopathic immune thrombocytopenia, may result in lower gastrointestinal bleeding or other internal bleeding in human [12]. Simply speaking, they are routinely measured as important clinical indicators to diagnose and monitor hematologic diseases and ascertain overall patient health condition.

The domestic pig is being increasingly exploited as an ideal model animal in human genetic diseases due to the high similarity with human physiological characteristics [13]. Therefore, discovering new loci for hematological traits and revealing their genetic mechanisms in domestic pig are conducive to the human blood disease. However, little is known about the association between genetic variation and hematological traits [1417]. To our knowledge, 239 genome-wide significant quantitative trait loci (QTL) have been identified so far, which only explained a small fraction of the genetic variance (http://www.animalgenome.org/cgi-bin/QTLdb/SS/index) [18]. In these identified QTLs, the confidence intervals are generally large (> 20 cM) [19] and harbor thousands of functional genes, thereby hampering the identification of plausible candidate genes. Compared with traditional QTL mapping strategies, single marker GWAS [20, 21] take advantage of linkage disequilibrium using high-density molecular markers rather than the linkage information using low-density markers in the intercross populations. Therefore, single marker GWAS could efficiently narrow down confidence interval of detected QTL and pick up the most associated markers for trait of interest. On the other hand, if the causative mutation is ancient, the LD between markers and mutated loci is too small to be captured with current marker density. Haplotype integrates linkage and linkage disequilibrium information [22] together, it is considered with the ability of overcoming the shortcoming in linkage and (or) single marker GWAS. Theoretically, haplotype analysis could acquire more accurate positions and shorter confidence intervals compared with separately performing linkage analysis or linkage disequilibrium analysis.

In this study, we conducted single marker GWAS and haplotype analysis of 15 hematological traits in Chinese Sutai population. The main purpose of the study is to reveal new loci associated with hematological trait and discover potential causative genes combining with biological and bioinformatics annotation. Furthermore, our result may also provide insights into the molecular basis of hematological trait in human.

Results

Phenotype statistics and SNP characteristics after quality control

The means, standard errors and coefficient of variation (C.V) of the phenotypic observations of the 15 hematological traits in the current experimental population are presented in Table 1. The C.V ranges from 3.73 to 38.71 as the minimum and maximum value for MCHC and PLT, respectively.

Table 1 Descriptive statistics of 15 hematological traits in the Sutai population

After quality control, none of the individuals had a genotyping call rate < 95%, resulting in 495 individuals remained for the association analyses. Additionally, 3610 SNPs with call rates < 90%, 16242 SNPs with minor allele frequency (MAF) < 0.05, 64 SNPs severely departing from Hardy Weinberg Equilibrium (HWE) (P-value < 10-5) and 149 makers exhibiting Mendelian inconsistency were excluded, remained a total of 44650 SNPs. We also removed 4864 SNPs, including unmapped SNPs or SNPs on the sex chromosome. Finally, a total of 39786 SNPs were remained for further analyses.

Erythrocyte traits

Single marker GWAS: In total, 141 significant SNPs (including 40 genome-wide and 101 suggestive SNPs) were identified for 8 erythrocyte traits: 5 for HCT, 1 for HGB, 40 for MCH, 22 for MCHC, 56 for MCV, 4 for RBC and 13 for RDW (Table 2 and Additional file 1: Table S1). All these 141 SNPs were located on SSC1, 2, 3, 4, 6,13 and 16; most of them located on SSC2 and SSC6 (Figure 1). No significant SNPs were detected for RDW-SD (Additional file 2: Figure S1). Eighty-three of the identified SNPs were located within 39 annotated genes, and 58 markers were located in region of 65 to 473458 bp away from the nearest annotated gene. In the 141 SNPs, 40 SNPs associated with at least two traits. They were mainly located on SSC2 and 13. And the most significant SNP (ss478944677) was associated with three erythrocyte traits: MCV (P-value =3.00 × 1011), MCH (P-value =1.10 × 109) and RDW (P-value =1.86 × 106).

Table 2 Description of lead SNPs showing significant association with hematological traits by GWAS
Figure 1
figure1

Manhattan plots for the single marker analysis of hematological traits surpass genome-wide significant threshold. log10(1/P-value) values are shown for all SNPs that passed quality control. The numbers indicate the chromosomes in the genome. The solid line and dotted line denotes the Bonferroni-corrected genome-wide and suggestive significant threshold, respectively. SNPs surpassing the genome-wide threshold are highlighted in pink and SNPs reaching the suggestive threshold in green. MCH: mean corpuscular hemoglobin; MCV: mean corpuscular volume; P-LCR: platelet-large cell ratio; RDW: red blood cell volume distribution width.

Haplotype analysis: Totally, 498 significant SNPs (including 154 genome-wide and 344 suggestive SNPs) were identified for 8 erythrocyte traits: 192 for HCT, 60 for MCH, 68 for MCV, 165 for RBC and 13 for RDW-SD (Table 3 and Additional file 3: Table S2). These significant SNPs were located on SSC1, 2, 4, 5, 7, 8, 9, 11, 12, 14 and 15 and most of them were located on SSC7 and 9 (Figure 2). No significant SNPs were detected in association with HGB, MCH and RDW (Additional file 4: Figure S2). The top SNP ss107842725 located in ENSSSCG00000001232 gene on SSC7 was associated with HCT, RBC and MCV. Furthermore, 38 of 154 genome-wide significant SNPs were located within the regions of 24 annotated genes and the others were located in the regions of the nearest known genes with the distance from 62 to 757213 bp.

Table 3 Description of lead SNPs showing significant association with hematological traits by LDLA
Figure 2
figure2

Manhattan plots for the haplotype analysis of hematological traits surpass genome-wide significant threshold. log10(1/P-value) values are shown for all SNPs that passed quality control. The numbers indicate the chromosomes in the genome. The solid line and dotted line denotes the Bonferroni-corrected genome-wide and suggestive significant threshold, respectively. SNPs surpassing the genome-wide threshold are highlighted in pink and SNPs reaching the suggestive threshold in green. HCT: hematocrit; RBC: red blood cell count.

White blood cell counts

Single marker GWAS: Analysis of white blood cell counts revealed two significant loci on SSC2 by single marker GWAS. The most significant SNP ss107857076 (P-value = 6.03 × 10-6) associated with WBC was located at 105499649 bp on SSC2 with a distance of 95033 bp away from ENSSSCG00000030166 gene. The remaining SNP ss131195511 was located at 101149437 bp on SSC2 and 277118 bp away from gene GPR98 (G protein-coupled receptor 98).

Haplotype analysis: One significant locus associated with WBC was identified by haplotype analysis. The SNP ss131152863 was located at 289943447 bp on SSC1 and 157600 bp away from TLR4 (toll-like receptor 4) gene.

Platelet traits

Single marker GWAS: Eighteen SNPs significantly associated with two platelet traits were detected by single marker GWAS: 13 for P-LCR and 5 for MPV. They were located on SSC2 and distributed within a 10.7 Mb region (54474152–65200938 bp). Both P-LCR and MPV shared the same top SNP of ss107886044 which was located in an annotated gene TRIM58 (tripartite motif containing 58) at 105499649 bp.

Haplotype analysis: No significant SNP was detected by haplotype analysis.

Discussion

The Sutai pigs were generated by intercross of Meishan (Erhualian) female and Duroc male for about 25 generations. Their genome was composed of a mosaic of small pieces of haplotype segments derived from both breed. As a result, their LD block was much smaller than classic QTL mapping populations [23]. Sutai pigs included two kinds of LD: LD between breeds created by intercross and LD within each breed created in the ancestor history, and they hence become very good experimental population for QTL mapping and single marker GWAS analysis.

Comparison with previous studies

By performing single marker GWAS and haplotype analyses, we identified 651 SNPs associated with the 15 hematologic traits. Of these SNPs, 253 located within known genes’ region, 265 located near to the annotated gene and 133 weren’t mapped to the current assembled genome (Sus Scrofa Build 10.2, http://asia.ensembl.org/index.html). So far, several papers reported single marker GWAS result for hematological traits in pig. Zhang et al. using similar association strategy revealed 185 genome-wide significant SNPs for 18 hematological traits in 1020 white Duroc x Erhualian F2 intercross [24]. Most of the identified significant SNPs were located on SSC8. Luo et al.[25] detected 62 genome-wide significant and 3 chromosome-wide significant SNPs associated with erythrocyte traits on a Large White × Chinese Min F2 intercross and most of them also retained on SSC8. Both of them pinpoint that KIT (v-kit Hardy-Zuckerman 4 feline sarcoma viral oncogene homolog) gene as the potential candidate. In our study, we didn’t detect any signal associated with erythrocyte traits in this region. KIT is essential for coat color while all individuals’ in our study is black. Hence there was no variation at KIT gene and of course without association signal. None significant SNP in Luo et al. and Zhang et al. was overlapped with our study. The reasons for the inconsistence by similar analysis strategy could be monomorphic at the causative locus, the population heterogeneity and the complex genetic background. These results also hint that hematological traits was a complex trait which affected by multiple genes. Wang et al.[26] identified 111 significant SNPs for 18 hematological traits after injected classical fever vaccine in 2 Western breeds and one Chinese synthetic breed by similar single marker association study. Their mapping result might include both QTL affecting immune responses and QTL affecting base hematological traits. Herein we found 9 SNPs on SSC6 were identical with the results of our present study, while none functional gene was posited in that region.

Comparison between single marker GWAS and haplotype analysis results

In this study, we performed both single marker GWAS and haplotype analysis to explore potential causal gene(s) for hematological traits in Chinese Sutai pigs. Only 9 SNPs located on SSC2 were overlapped by both analyses for MCH. The basic principle of single marker GWAS was to compare phenotypic differences grouped by alleles. If the marker density was not high enough, the significant SNPs may lose because of low LD between markers and causative mutation. However, haplotype will surmount this disadvantage. Druet and Georges [27] have fully descripted haplotype analysis, which took advantage of recent and ancestral recombination events simultaneously. In here, we used haplotype analysis and identified 490 SNPs located on SSC1, 2, 4, 5, 7, 8, 9, 11, 12, 14 and 15 which can’t be detected by single marker GWAS. However, one drawback of haplotype analysis is the reduction of detection power, because its degree of freedom is generally bigger than single marker analysis. Zhang et al.[28] also pinpoint this phenomenon due to the degree of freedom. However, the balance between increasing LD and decreasing power by the degree of freedom is hard to weight. Moreover, LD across whole genome is inhomogeneous – there are high LD in some regions and low LD in other regions. In this case, we recommend performing both single marker and haplotype analyses strategies to capture more associated SNPs. We obtained 141 significant SNPs by single marker analysis and 498 SNPs by haplotype analyses for 8 erythrocyte traits. In together, 651 significant were identified associating with hematological traits, which was more than any one analyses strategy.

Possible pleiotropic QTLs

The patterns of Manhattan plots of MCH, MCV and P-LCR were similar, and they shared a common region ranging from 54.47 Mb to 55.24 Mb containing three SNPs (ss131191392, ss478944677 and ss131085967) on SSC2. MCH and MCV are parameters reflecting average weight of hemoglobin per RBC and average volume of RBC, respectively. By analysis of the correlation among the hematological traits (Additional file 5: Table S3), high correlation between the two traits was observed (r = 0.804, P-value < 1.0 × 10-16). This result implied the QTL on SSC2 might be pleiotropic. Marker ss107842725, located at 24777963 bp on SSC7, was the top SNP associated with HCT, MCH and RBC. The Manhattan plot also explored very similar patterns for the three phenotypes. HCT, MCH and RBC mainly measure fluctuation of red blood cell and they may segregate dependently. Our results indicated that pleiotropic QTL was common on hematological traits. In clinical diagnosis, the three parameters (HCT, MCH and RBC) could be integrated together for more precisely diagnose.

Potential candidate functional genes

In total, we identified 161 significant SNPs on 7 different chromosomes associated with hematological traits by single marker GWAS (Additional file 1: Table S1). Among these SNPs, 25 SNPs were found within 14 annotated genes from 52.14 to 90.17 Mb. Through checking these annotated gene functions, we eventually selected four genes as potential candidate genes. The four genes, TRIM58, CPAMD8 (C3 and PZP-like, alpha-2-macroglobulin domain containing 8), ABCA7 (ATP-binding cassette, sub-family A (ABC1), member 7) and JAK3 (Janus kinase 3), were functionally associated with hematological related cells or immune function.

The SNP ss107886044 located in TRIM58 gene explained 15.43% (Table 2) of phenotypic variants of P-LCR. Christopher et al. regarded TRIM58 as an E3 ubiquitin ligase that regulated terminal erythroid cell cycles and enucleation [29]. Moreover, the TRIM58 protein was involved in pathogen-recognition [30] and the regulation of innate immune responses [31]. Therefore, the TRIM58 gene could be regarded as a strong candidate gene controlling P-LCR. In addition to TRIM58, the SNP marker (ss131190955) within the CPAMD8 gene also showed high association signal with a P-value of 1.33 × 10-10. The CPAMD8 gene was highly conserved, which may have similar function like other members of the C3/α2M family and also be involved in innate immunity [3234]. A promising gene ABCA7 was a member of the super family of ATP-binding cassette (ABC) transporters, expression of which was induced during vitro differentiation of human monocytes into macrophages. Besides, ABCA7 mRNA was detected predominantly in myelo-lymphatic tissues with highest expression in peripheral leukocytes [35, 36]. JAK3 was predominantly expressed in hematopoietic cells, such as NK cells, T cells and B cells [37] and transduced a signal in response to its activation. Furthermore, mutations which abrogated JAK3 might cause an autosomal SCID (severe combined immunodeficiency disease) [38].

By haplotype analysis, we identified 154 genome-wide significantly loci mainly SSC7 and SSC9. Among them, 50 significant SNPs for HCT were found within 34 annotated genes and 4 significant SNPs for RBC were found within 4 annotated genes. In these annotated genes, three genes TRIM26 on SSC7, TRIM21 (tripartite motif containing 21) and NUP98 (nucleoporin 98 kDa) on SSC9 were picked up as potential candidates by checking their gene functions. These genes were functionally associated with hematological related cells or immune function.

The TRIM26, encoding a member of the tripartite motif (TRIM) family, was located within the SLA region [39]. Lee et al. also speculated that the TRIM26 gene played essential roles in the human immune response because of its predicted protein function [40]. In addition to TRIM26, the TRIM21 gene also belonged to the tripartite motif (TRIM) family. It was an E3 ubiquitin ligase for IFN regulatory factor IRF3 and IRF8 with the function of innate and adaptive immunity [41]. Yang et al. demonstrated that TRIM21 interacts with PIN1 mediates the ubiquitination and degradation of IRF3 during virus infection [42]. Besides, it was reported that TRIM21 may regulate T-cell activation or proliferation, since overexpression of TRIM21 had been shown to increase IL-2 production in CD28-stimulated Jurkat T cells [43]. Therefore we could regard the TRIM21 gene, involved in both physiological immune responses and pathological autoimmune processes [44], as a strong candidate gene. The NUP98 fusion proteins had been shown to inhibit differentiation of hematopoietic precursors and to increase self-renewal of hematopoietic stem or progenitor cells [45]. The NUP98 gene also was known to be fused to at least 28 different partner genes in patients with hematopoietic malignancies, including acute myeloid leukemia, chronic myeloid leukemia in blast crisis, myelodysplastic syndrome, acute lymphoblastic leukemia, and bilineage/biphenotypic leukemia.

In all identified genes, we specially pointed out three genes (TRIM58, TRIM26 and TRIM21), which belonged to the same gene family. The three genes executed similar function of innate and adaptive immune and communicated together in the immune network system. Our result revealed a series of key driver genes in the immune network system.

Conclusions

In summary, we identified 651 SNPs, some of which were pleiotropic. Such as three SNPs on SSC2 associated with MCV, MCH and P-LCR and ss107842725 on SSC7 associated with HCT, MCH and RBC. What’s more, we selected 7 genes as potential candidates based on their functional annotations, positions and reported expression variation. Especially, three strong candidate genes (TRIM58, TRIM26 and TRIM21) may be the key driver genes in the immune network system. These findings will conduct further studies to examine the identified SNPs in other diverse population and pursue functional validation for identification of the causal mutation.

Methods

Ethics statement

All procedures involving animals followed the guidelines for the care and use of experimental animals approved by the State Council of the People’s Republic of China.

Study populations and phenotype measurement

The Sutai population comprised with 436 offspring of 4 boars and 55 sows. Each boar mated with 13 to 15 sows to make the family structure in balance. There were three batches of piglets which were almost born in three different months (April, June and July, 2011) at Sutai Pig Breeding Center in Suzhou city. At the age of 2–3 month, then the piglets were transferred to a farm in Nanchang city. All Sutai piglets were castrated and weaned at 18 days and 28 days after birth, respectively. They were fed with same diet (formulated according to age) under a standardized feeding and management regimen, and given free access to water. At 240 ± 6 days of age, a total of 436 Sutai offspring including 206 gilts and 230 barrows were slaughtered at a commercial abattoir.

Blood samples of 5 ml were immediately collected from each animal when it was slaughtered and directly injected into eppendorf tubes containing 30 ul of 20% EDTA in polybutadiene-styrene. A standard set of hematological data were recorded using a CD1700 whole blood analyzer (Abbott, USA) in 24 h postmortem at the First Affiliated Hospital of NanChang University, China. Fifteen hematological parameters including 8 baseline erythroid traits (hematocrit (HCT), hemoglobin (HGB), mean corpuscular hemoglobin (MCH), mean corpuscular hemoglobin concentration (MCHC), mean corpuscular volume (MCV), red blood cell count (RBC), red blood cell volume distribution width-SD (RDW-SD), and red blood cell volume distribution width (RDW)), 3 leukocyte traits (lymphocyte count (LYM), lymphocyte count percentage (LYMA), and white blood cell count (WBC)), and 4 platelet traits (platelet distribution width (PDW), platelet count (PLT), platelet-large cell ratio (P-LCR) and mean platelet volume (MPV)) were used for performing single marker GWAS. The correlations between 15 hematological parameters were performed by R psych package (http://personality-project.org/r/psych.manual.pdf).

Genotyping and quality control

Genomic DNA was extracted from ear tissues using a standard phenol/chloroform method [46]. All DNA samples were qualified and standardized into a final concentration of 50 ng/ul. A total of 436 offspring and their 59 parents in the Sutai pedigree were genotyped for the Porcine SNP60 Beadchips on an iScan System (Illumina, USA) following the manufacturer’s protocol. Quality control was carried out using PLINK (version 1.07) [47] and executed to exclude SNPs with parameter of call rate < 90%, minor allele frequency (MAF) < 5%, severely departed from HWE (P-value < 10-5) and Mendelian inconsistency rate > 10%. Moreover, individuals with missing genotypes > 10% or Mendelian errors > 5% were discarded for further analysis.

Statistical analyses

The genome-wide and suggestive significance thresholds in the two association strategies were determined by the Bonferroni correction, in which the conventional P-value was divided by the number of tests performed [48]. A SNP was considered to have genome-wide significance at P-value < 0.05/N and have suggestive significance at P-value < 1/N, where N is the number of SNPs tested in the analyses. The corresponding thresholds were set as 1.26 × 10-6 (0.05/39786) and 2.51 × 10-5 (1/39786) in this study.

Single marker GWAS

The linear tendency of allelic and phenotypic traits was tested using a general linear mixed model for each SNP [4951]. The model included a random polygenic effect and the variance-covariance matrix was proportionate to genome-wide identity-by-state [52]. The model was described as following: Y = Xb + Sα + Zu + e, where Y is the vector of phenotypes, b is the estimator of fixed effects including sex and batch, α is the SNP substitution effect and u is the random additive genetic effect following multinomial distribution u ~ N (0, Gσα2), in here G is the genomic similarity matrix that was constructed based on SNP markers as described in Eding et al.[53], and σα2 is the polygenetic additive variance. X, S and Z are the incidence matrices for b, α and u, respectively. e is a vector of residual errors with a distribution of N (0, Iσe2). The single marker GWAS were conducted by GenABEL packages in the R software [54, 55].

Haplotype analysis

The haplotypes were constructed following Druet & Georges, by using a Hidden Markov Model via PHASEBOOK [27] that assumes the existence of a predetermined number of ancestral haplotype states (K = 20) from which all haplotypes in the population are derived [56]. The statistical model used for the haplotype analysis was identical to that of single marker GWAS except that a haplotype effect was fitted instead of a SNP effect [57]. The haplotype was followed approximately the approach of Meuwissen and Goddard [22, 31, 58, 59], except that haplotypes were assumed to be completely uncorrelated, instead of fitting a more differentiating identity by descent (IBD) matrix G.

Phenotypic variants analysis

The fraction of phenotype variances explained by detected SNP was computed by following formula:

Var% = M S reduce 1 - M S full M S reduce × 100

Where MS full , MS reduce1 and MS reduce were the mean square (MS) in the linear models including three effects (mean, sex and SNP), including two effects (mean and sex) and only including mean, respectively.

Authors’ information

Feng Zhang and Zhiyan Zhang are co-first authors.

References

  1. 1.

    Muller M, Brem G: Disease resistance in farm animals. Experientia. 1991, 47 (9): 923-934. 10.1007/BF01929883.

    PubMed  CAS  Article  Google Scholar 

  2. 2.

    Oddgeirsson O, Simpson SP, Morgan AL, Ross DS, Spooner RL: Relationship between the bovine major histocompatibility complex (BoLA), erythrocyte markers and susceptibility to mastitis in Icelandic cattle. Anim Genet. 1988, 19 (1): 11-16.

    PubMed  CAS  Article  Google Scholar 

  3. 3.

    Tullis JL: Separation and purification of leukocytes and platelets. Blood. 1952, 7 (9): 891-896.

    PubMed  CAS  Google Scholar 

  4. 4.

    Beutler B: Innate immunity: an overview. Mol Immunol. 2004, 40 (12): 845-859. 10.1016/j.molimm.2003.10.005.

    PubMed  CAS  Article  Google Scholar 

  5. 5.

    van de Vosse E, van Dissel JT, Ottenhoff TH: Genetic deficiencies of innate immune signalling in human infectious disease. Lancet Infect Dis. 2009, 9 (11): 688-698. 10.1016/S1473-3099(09)70255-5.

    PubMed  Article  Google Scholar 

  6. 6.

    Arosa FA, Pereira CF, Fonseca AM: Red blood cells as modulators of T cell growth and survival. Curr Pharm Des. 2004, 10 (2): 191-201. 10.2174/1381612043453432.

    PubMed  CAS  Article  Google Scholar 

  7. 7.

    Nikinmaa M: Oxygen and carbon dioxide transport in vertebrate erythrocytes: an evolutionary change in the role of membrane transport. J Exp Biol. 1997, 200 (Pt 2): 369-380.

    PubMed  CAS  Google Scholar 

  8. 8.

    Elzey BD, Sprague DL, Ratliff TL: The emerging role of platelets in adaptive immunity. Cell Immunol. 2005, 238 (1): 1-9. 10.1016/j.cellimm.2005.12.005.

    PubMed  CAS  Article  Google Scholar 

  9. 9.

    Klinger MH, Jelkmann W: Review: role of blood platelets in infection and inflammation. J Interferon Cytokine Res. 2002, 22 (9): 913-922. 10.1089/10799900260286623.

    PubMed  CAS  Article  Google Scholar 

  10. 10.

    Smyth SS, McEver RP, Weyrich AS, Morrell CN, Hoffman MR, Arepally GM, French PA, Dauerman HL, Becker RC, Platelet Colloquium P: Platelet functions beyond hemostasis. J Thromb Haemost. 2009, 7 (11): 1759-1766. 10.1111/j.1538-7836.2009.03586.x.

    PubMed  CAS  Article  Google Scholar 

  11. 11.

    Yeaman MR: Platelets in defense against bacterial pathogens. Cell Mol Life Sci. 2010, 67 (4): 525-544. 10.1007/s00018-009-0210-4.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  12. 12.

    Idiopathic thrombocytopenic purpura. [http://www.mindheal.org/itp.html]

  13. 13.

    Swindle M, Makin A, Herron A, Clubb F, Frazier K: Swine as models in biomedical research and toxicology testing. Vet Pathol. 2012, 49 (2): 344-356. 10.1177/0300985811402846.

    PubMed  CAS  Article  Google Scholar 

  14. 14.

    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. J Immunol. 1998, 161 (2): 829-835.

    PubMed  CAS  Google Scholar 

  15. 15.

    Reiner G, Fischer R, Hepp S, Berge T, Kohler F, Willems H: Quantitative trait loci for red blood cell traits in swine. Anim Genet. 2007, 38 (5): 447-452. 10.1111/j.1365-2052.2007.01629.x.

    PubMed  CAS  Article  Google Scholar 

  16. 16.

    Reiner G, Fischer R, Hepp S, Berge T, Kohler F, Willems H: Quantitative trait loci for white blood cell numbers in swine. Anim Genet. 2008, 39 (2): 163-168. 10.1111/j.1365-2052.2008.01700.x.

    PubMed  CAS  Article  Google Scholar 

  17. 17.

    Wattrang E, Almqvist M, Johansson A, Fossum C, Wallgren P, Pielberg G, Andersson L, Edfors-Lilja I: Confirmation of QTL on porcine chromosomes 1 and 8 influencing leukocyte numbers, haematological parameters and leukocyte function. Anim Genet. 2005, 36 (4): 337-345. 10.1111/j.1365-2052.2005.01315.x.

    PubMed  CAS  Article  Google Scholar 

  18. 18.

    Hu ZL, Park CA, Wu XL, Reecy JM: Animal QTLdb: an improved database tool for livestock animal QTL/association data dissemination in the post-genome era. Nucleic Acids Res. 2013, 41 (Database issue): D871-D879.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  19. 19.

    Pearson TA, Manolio TA: How to interpret a genome-wide association study. JAMA. 2008, 299 (11): 1335-1344. 10.1001/jama.299.11.1335.

    PubMed  CAS  Article  Google Scholar 

  20. 20.

    Hu Z, Xu S: PROC QTL-a SAS procedure for mapping quantitative trait loci. Int J Plant Genomics. 2009, 2009: 141234-

    PubMed  PubMed Central  Article  Google Scholar 

  21. 21.

    Terwilliger JD: A powerful likelihood method for the analysis of linkage disequilibrium between trait loci and one or more polymorphic marker loci. Am J Hum Genet. 1995, 56 (3): 777-

    PubMed  CAS  PubMed Central  Google Scholar 

  22. 22.

    Meuwissen TH, Karlsen A, Lien S, Olsaker I, Goddard ME: Fine mapping of a quantitative trait locus for twinning rate using combined linkage and linkage disequilibrium mapping. Genetics. 2002, 161 (1): 373-379.

    PubMed  CAS  PubMed Central  Google Scholar 

  23. 23.

    Zhuang Q, Liang W, Cong S, Chen G, Shi Q: New breed of China Lean type pig - Sutai pig. Modernizing Agric. 2007, 12: 018-

    Google Scholar 

  24. 24.

    Zhang Z, Hong Y, Gao J, Xiao S, Ma J, Zhang W, Ren J, Huang L: Genome-wide association study reveals constant and specific loci for hematological traits at three time stages in a white Duroc x Erhualian F2 resource population. PLoS One. 2013, 8 (5): e63665-10.1371/journal.pone.0063665.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  25. 25.

    Luo W, Chen S, Cheng D, Wang L, Li Y, Ma X, Song X, Liu X, Li W, Liang J: Genome-wide association study of porcine hematological parameters in a large white × Minzhu F2 resource population. Int J Biol Sci. 2012, 8 (6): 870-

    PubMed  PubMed Central  Article  Google Scholar 

  26. 26.

    Wang JY, Luo YR, Fu WX, Lu X, Zhou JP, Ding XD, Liu JF, Zhang Q: Genome-wide association studies for hematological traits in swine. Anim Genet. 2013, 44 (1): 34-43. 10.1111/j.1365-2052.2012.02366.x.

    PubMed  Article  Google Scholar 

  27. 27.

    Druet T, Georges M: A hidden markov model combining linkage and linkage disequilibrium information for haplotype reconstruction and quantitative trait locus fine mapping. Genetics. 2010, 184 (3): 789-798. 10.1534/genetics.109.108431.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  28. 28.

    Zhang Z, Guillaume F, Sartelet A, Charlier C, Georges M, Farnir F, Druet T: Ancestral haplotype-based association mapping with generalized linear mixed models accounting for stratification. Bioinformatics. 2012, 28 (19): 2467-2473. 10.1093/bioinformatics/bts348.

    PubMed  CAS  Article  Google Scholar 

  29. 29.

    Trim58 is a putative E3 ubiquitin ligase that functions in late stage erythropoiesis. [https://ash.confex.com/ash/2012/webprogram/Paper51089.html]

  30. 30.

    Ozato K, Shin DM, Chang TH, Morse HC: TRIM family proteins and their emerging roles in innate immunity. Nat Rev Immunol. 2008, 8 (11): 849-860. 10.1038/nri2413.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  31. 31.

    Kawai T, Akira S: Regulation of innate immune signalling pathways by the tripartite motif (TRIM) family proteins. EMBO Mol Med. 2011, 3 (9): 513-527. 10.1002/emmm.201100160.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  32. 32.

    Armstrong PB, Quigley JP: Alpha2-macroglobulin: an evolutionarily conserved arm of the innate immune system. Dev Comp Immunol. 1999, 23 (4–5): 375-390.

    PubMed  CAS  Article  Google Scholar 

  33. 33.

    Studd J, Blainey J, Bailey D: A study of serum protein changes in late pregnancy and identification of the pregnancy zone protein using antigen antibody crossed immunoelectrophoresis. BJOG. 1970, 77 (1): 42-51. 10.1111/j.1471-0528.1970.tb03406.x.

    CAS  Article  Google Scholar 

  34. 34.

    Volanakis JE: The role of complement in innate and adaptive immunity. Curr Top Microbiol. 2002, 266: 41-56.

    CAS  Google Scholar 

  35. 35.

    Iwamoto N, Abe-Dohmae S, Sato R, Yokoyama S: ABCA7 expression is regulated by cellular cholesterol through the SREBP2 pathway and associated with phagocytosis. J Lipid Res. 2006, 47 (9): 1915-1927. 10.1194/jlr.M600127-JLR200.

    PubMed  CAS  Article  Google Scholar 

  36. 36.

    Kaminski WE, Orsó E, Diederich W, Klucken J, Drobnik W, Schmitz G: Identification of a novel human sterol-sensitive ATP-binding cassette transporter (ABCA7). Biochem Biophys Res Commun. 2000, 273 (2): 532-538. 10.1006/bbrc.2000.2954.

    PubMed  CAS  Article  Google Scholar 

  37. 37.

    Leonard WJ, O’Shea JJ: JAKS AND STATS: biological implications. Ann Rev Immunol. 1998, 16 (1): 293-322. 10.1146/annurev.immunol.16.1.293.

    CAS  Article  Google Scholar 

  38. 38.

    Russell SM, Tayebi N, Nakajima H, Riedy MC, Roberts JL, Aman MJ, Migone TS, Noguchi M, Markert ML, Buckley RH, O'Shea JJ, Leonard WJ: Mutation of Jak3 in a patient with SCID: essential role of Jak3 in lymphoid development. Science. 1995, 270 (5237): 797-800. 10.1126/science.270.5237.797.

    PubMed  CAS  Article  Google Scholar 

  39. 39.

    de Jong S, van Eijk KR, Zeegers DW, Strengman E, Janson E, Veldink JH, van den Berg LH, Cahn W, Kahn RS, Boks MP: Expression QTL analysis of top loci from GWAS meta-analysis highlights additional schizophrenia candidate genes. Eur J Hum Genet. 2012, 20 (9): 1004-1008. 10.1038/ejhg.2012.38.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  40. 40.

    Lee JS, Bae JS, Kim JH, Kim JY, Park TJ, Pasaje CF, Park BL, Cheong HS, Jang AS, Uh ST, Park CS, Shin HD: Association study between TRIM26 polymorphisms and risk of aspirin-exacerbated respiratory disease. Int J Mol Med. 2012, 29 (5): 927-933.

    PubMed  CAS  Google Scholar 

  41. 41.

    Yoshimi R, Chang TH, Wang H, Atsumi T, Morse HC, Ozato K: Gene disruption study reveals a nonredundant role for TRIM21/Ro52 in NF-kappaB-dependent cytokine expression in fibroblasts. J Immunol. 2009, 182 (12): 7527-7538. 10.4049/jimmunol.0804121.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  42. 42.

    Yang K, Shi H-X, Liu X-Y, Shan Y-F, Wei B, Chen S, Wang C: TRIM21 is essential to sustain IFN regulatory factor 3 activation during antiviral response. J Immunol. 2009, 182 (6): 3782-3792. 10.4049/jimmunol.0803126.

    PubMed  CAS  Article  Google Scholar 

  43. 43.

    Ishii T, Ohnuma K, Murakami A, Takasawa N, Yamochi T, Iwata S, Uchiyama M, Dang NH, Tanaka H, Morimoto C: SS-A/Ro52, an autoantigen involved in CD28-mediated IL-2 production. J Immunol. 2003, 170 (7): 3653-3661.

    PubMed  CAS  Article  Google Scholar 

  44. 44.

    Yoshimi R, Ishigatsubo Y, Ozato K: Autoantigen TRIM21/Ro52 as a possible target for treatment of systemic lupus erythematosus. Int J Rheum. 2012, 2012: 718237-

    Google Scholar 

  45. 45.

    Gough SM, Slape CI, Aplan PD: NUP98 gene fusions and hematopoietic malignancies: common themes and new biologic insights. Blood. 2011, 118 (24): 6247-6257. 10.1182/blood-2011-07-328880.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  46. 46.

    Moore D: Preparation and analysis of DNA. Short Protocols in Molecular Biology, Volume 1. Edited by: Ausubel FM, Brent R, Kingston RE, Moore DD, Seidman J, Smith JA, Struhl K. 1995, New York: John Wiley & Sons Inc, 24: 68, 3

    Google Scholar 

  47. 47.

    Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ, Sham PC: PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007, 81 (3): 559-575. 10.1086/519795.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  48. 48.

    Yang Q, Cui J, Chazaro I, Cupples LA, Demissie S: Power and type I error rate of false discovery rate approaches in genome-wide association studies. BMC Genet. 2005, 6 (1): S134-10.1186/1471-2164-6-134.

    PubMed  PubMed Central  Article  Google Scholar 

  49. 49.

    Bradbury PJ, Zhang Z, Kroon DE, Casstevens TM, Ramdoss Y, Buckler ES: TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics. 2007, 23 (19): 2633-2635. 10.1093/bioinformatics/btm308.

    PubMed  CAS  Article  Google Scholar 

  50. 50.

    Breslow NE, Clayton DG: Approximate inference in generalized linear mixed models. J Am Stat Assoc. 1993, 88 (421): 9-25.

    Google Scholar 

  51. 51.

    Yu J, Pressoir G, Briggs WH, Bi IV, Yamasaki M, Doebley JF, McMullen MD, Gaut BS, Nielsen DM, Holland JB: A unified mixed-model method for association mapping that accounts for multiple levels of relatedness. Nat Genet. 2005, 38 (2): 203-208.

    PubMed  Article  Google Scholar 

  52. 52.

    Hayes BJ, Goddard ME: Technical note: prediction of breeding values using marker-derived relationship matrices. J Anim Sci. 2008, 86 (9): 2089-2092. 10.2527/jas.2007-0733.

    PubMed  CAS  Article  Google Scholar 

  53. 53.

    Eding H: Marker‒based estimates of between and within population kinships for the conservation of genetic diversity. J Anim Breed Genet. 2001, 118 (3): 141-159. 10.1046/j.1439-0388.2001.00290.x.

    CAS  Article  Google Scholar 

  54. 54.

    Aulchenko YS, de Koning D-J, Haley C: Genomewide rapid association using mixed model and regression: a fast and simple method for genomewide pedigree-based quantitative trait loci association analysis. Genetics. 2007, 177 (1): 577-585. 10.1534/genetics.107.075614.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  55. 55.

    Aulchenko YS, Ripke S, Isaacs A, van Duijn CM: GenABEL: an R library for genome-wide association analysis. Bioinformatics. 2007, 23 (10): 1294-1296. 10.1093/bioinformatics/btm108.

    PubMed  CAS  Article  Google Scholar 

  56. 56.

    Sartelet A, Druet T, Michaux C, Fasquelle C, Geron S, Tamma N, Zhang Z, Coppieters W, Georges M, Charlier C: A splice site variant in the bovine RNF11 gene compromises growth and regulation of the inflammatory response. PLoS Genet. 2012, 8 (3): e1002581-10.1371/journal.pgen.1002581.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  57. 57.

    Grindflek E, Lien S, Hamland H, Hansen MH, Kent M, van Son M, Meuwissen TH: Large scale genome-wide association and LDLA mapping study identifies QTLs for boar taint and related sex steroids. BMC Genomics. 2011, 12: 362-10.1186/1471-2164-12-362.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  58. 58.

    Goddard M: Mapping multiple QTL by combined linkage disequilibrium/linkage analysis in outbred populations. 7th WCGALP, Montpellier, France, August. 2002, INRA, 1-4. Session 21: 2002

    Google Scholar 

  59. 59.

    Meuwissen TH, Goddard ME: Prediction of identity by descent probabilities from marker-haplotypes. Genet Sel Evol. 2001, 33 (6): 605-634. 10.1186/1297-9686-33-6-605.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

Download references

Acknowledgements

This study was supported by National Natural Science Foundation of China (31200926) and National 973 Program (2012CB722502).

Author information

Affiliations

Authors

Corresponding author

Correspondence to Lusheng Huang.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

LH conceived and led the coordination of the study. ZZ and FZ led the data analysis and the manuscript preparation. XY, WZ, HC and YH contributed to blood collection and slaughter. FZ, WZ, HC and YH directed the genotyping work and recording hematological data. XY, ZZ and FZ interpreted the results and contributed to edit the manuscript. All authors read and approved the final manuscript.

Feng Zhang, Zhiyan Zhang contributed equally to this work.

Electronic supplementary material

Additional file 1: Table S1: Description of all identified SNPs showing significant association with hematological traits by single marker GWAS. (DOCX 30 KB)

Additional file 2: Figure S1: Manhattan plots for the single marker analysis of hematological traits surpass suggestive significant threshold. log10 (1/P-value) values are shown for all SNPs that passed quality control. The solid line and dotted line denotes the Bonferroni-corrected genome-wide and suggestive significant threshold respectively. SNPs reaching the suggestive threshold are highlighted in green. HCT: hematocrit; HGB: hemoglobin; MCHC: mean corpuscular hemoglobin content; RBC: red blood cell; WBC: white blood cell count; MPV: mean platelet volume. (TIFF 2 MB)

Additional file 3: Table S2: Description of all identified SNPs showing significant association with hematological traits by haplotype analysis. (DOCX 44 KB)

Additional file 4: Figure S2: Manhattan plots for the haplotype analysis of hematological traits surpass suggestive significant threshold. log10 (1/P-value) values are shown for all SNPs that passed quality control. The solid line and dotted line denotes the Bonferroni-corrected genome-wide and suggestive significant threshold respectively. SNPs reaching the suggestive threshold are highlighted in green. MCH: mean corpuscular hemoglobin; MCV: mean corpuscular volume; RDW-SD: red blood cell volume distribution width-SD; WBC: white blood cell count. (TIFF 1 MB)

Additional file 5: Table S3: Description of correlation and P-value among the 15 hematological traits. (DOCX 23 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

Reprints and Permissions

About this article

Cite this article

Zhang, F., Zhang, Z., Yan, X. et al. Genome-wide association studies for hematological traits in Chinese Sutai pigs. BMC Genet 15, 41 (2014). https://doi.org/10.1186/1471-2156-15-41

Download citation

Keywords

  • Single marker GWAS
  • Haplotype analysis
  • Hematological traits
  • Pig