Diplotyper: diplotype-based association analysis

Causal mutations for health conditions with genetic bases can be identified through finding associations with haplotypes, a form of correlation known as linkage disequilibrium (LD) [1]. Investigating significant haplotype structure has become a routine research activity. The Haploview tool provides computation of LD and population haplotype patterns from genotype data [2]. The PLINK tool set [3] accomplishes diverse functions including a module carrying out Expectation-Maximization (EM) algorithm [4]. PLINK focuses on fast calculations with large datasets. WHAP was developed to perform haplotype-based association analysis in population and family samples using single nucleotide polymorphism (SNP) data [5]. An additional software tool was elaborated for carrying out haplotype association analysis in unrelated individuals [6]. To provide a detailed genome structure, a recloning system [7] was developed to obtain the sequences of 20 haplotypes from a chimpanzee and a gorilla, across human leukocyte antigen (HLA) genes.

Meanwhile, rare haplotypes have been investigated to identify their roles in influencing disease susceptibility. Experimental data showed that two rare haplotypes of parathyroid hormone-related peptide receptor type 1 and vitamin D receptor genes, with frequencies of 1.1% and 2.9%, respectively, were significantly associated with osteoporosis phenotypes (P = 4.2 × 10^-6 and P = 1.6 × 10^-4, respectively) [8]. Recently, haplotypes in the regulatory regions of the HLA-G gene were examined to recognize possible associations with the implantation outcome in couples undergoing assisted reproduction treatments (ART). The results revealed a complete absence of some haplotypes in couples undergoing ART [9].

Notably, Durrant et al. proposed a novel approach to investigate associations between diseases and haplotype clusters in a logistic regression framework through cladistic analysis of SNP haplotypes. Substantial increases in power over single-locus tests were demonstrated by the simulation study. Their empirical data showed that a haplotype cluster that consisted of two haplotypes had the strongest effect on Cystic Fibrosis (OR = 96.8) [10]. Luo et al. used a novel analysis, diplotype trend regression (DTR) analysis, to investigate associations between certain diplotypes of alcohol dehydrogenase and aldehyde dehydrogenase genes, and alcohol dependence. They demonstrated that DTR outperformed other conventional association methods [11]. Both articles indicated that our new algorithm might provide a synergistic effect through combining analyses based on both haplotype clusters and diplotypes.

Here, we propose a novel method to investigate associations between diplotypes and diseases. We define a haplotype cluster as a set of haplotypes. We also define a diplotype as a haplotype cluster pair, the definition of which is extended from a haplotype pair. The first step of our method uses the Haploview tool to generate all possible haplotypes. Second, all possible haplotype pairs (diplotypes) from SNP genotypes of all samples are generated by PLINK. Third, all possible haplotype clusters are generated by our clustering algorithm from the haplotypes produced in the first step. Fourth, the patterns of all possible diplotypes are generated from those haplotype clusters. Fifth, to calculate regression by PLINK, the diplotypes of the samples are transformed into AA, AB, or BB formats according to the diplotype patterns produced in the fourth step. Finally, PLINK was used with a regression model to obtain the association results. All of these procedures are performed automatically by the software we developed, named Diplotyper, which was implemented in Python 2.7.

We applied this method to an association study between high-density lipoprotein cholesterol (HDL-C) and the hepatic lipase (HL) gene. HL is involved in lipoprotein metabolism through its bridging function, which facilitates the interaction between lipoproteins and lipoprotein receptors, and its activity plays an important role in plasma lipoprotein metabolism and the atherosclerotic process [12]. HL plays an important role in both reverse cholesterol transport and non-cholesterol-dependent mechanisms involved with HDL [13, 14]. Changes in HL activity can be associated with alterations in lipoprotein composition, which may contribute to the development of atherosclerosis [12, 14]. Low HDL-C levels are risk factors for coronary heart and cardiovascular diseases [15, 16]. Extensive research has provided evidence that increasing HDL-C levels can reduce the risk of cardiovascular disease [17–20]. The risk of developing coronary heart disease can be determined by the levels of HDL-C [21]. In particular, a recent report, based on the Korea National Health and Nutrition Examination Surveys I, II, III and IV, shows a growing prevalence of dyslipidemia and hypertriglyceridemia in Korea during the last decade [22].

LIPC encodes hepatic triglyceride lipase, which is expressed mainly in the liver and is located on 15q21-q23, where it spans 171 kb and comprises nine exons and eight introns. Two SNPs [-514C > T (rs1800588) and -250G > A (rs2070895)] in the promoter region are in almost perfect LD (R² = 0.97) [23, 24] and both the promoter SNP (rs1800588) and the intronic SNP (rs261332) has strong LD (R² = 0.92) in HapMap CEU (Caucasians of European descent from Utah) database [25, 26]. The intronic SNPs rs261332 [26] and rs11858164 [27], and the promoter SNPs rs1532085 [28] and rs10468017 [29], were associated with HDL-C levels in genome-wide association studies. The promoter SNP rs1800588 showed an increase in HDL-C of 0.04 mmol/l in the CT group and 0.09 mmol/l in the TT group, compared with the CC carriers [30]. The promoter SNP rs2070895 showed a highly significant association with a 0.057 mmol/l increase in HDL-C per A allele (P = 8 × 10^-10) [24]. We investigated possible associations between LIPC SNPs or diplotypes and HDL-C levels in a Korean population consisting of 7,536 individuals.

The total number of tests of association in the additive model were 224 for HDL-C (data not shown). As shown in Table 4, out of all LIPC 38 intronic SNPs available to the research community through the KARE project, only three SNPs were selected, since these were the only ones with strong associations with HDL-C that were able to withstand Bonferroni correction (P < 2.23 × 10^-4, see Methods), on the basis of their diplotypes as well as their SNP genotypes. The SNP rs261332 had the strongest association (P value = 3.03 × 10^-12) with HDL-C levels. Figure 1 shows the LD block and the haplotypes for this block that exhibited frequencies of greater than the 1% threshold.

SNP	Minor; Major	Genotype frequency	Effect	Standard error	P value
rs11631342	G; A	29/923/6583	0.037	0.0083	5.66E-06
rs6494005	G; A	569/2891/4018	-0.023	0.0047	7.43E-07
rs261332	T; C	350/2617/4556	0.035	0.0049	3.03E-12

Associations of single SNPs with HDL-C levels.

Effect sizes (mmol/l) shown are β-coefficients and measured as additive effects, which correspond to the average change in HDL-C levels when one major allele is replaced with one minor allele.

Table 5 represents the results of the Diplotyper tool, which indicate very diverse associations between diplotypes and HDL-C levels. The top 10 results out of the 40 diplotypes (data not shown) in block 1, in terms of P value, are shown. The abbreviation b1 means block 1, H12 means the twelfth haplotype cluster, and OH12(*) represents all other haplotypes (including haplotypes below 1% frequency), except the H12 haplotype cluster; b1_H12 is a minor haplotype cluster and b1_OH12(*) a major haplotype cluster. The number of one homozygous diplotype (a b1_H12 pair) was 612, the number of a heterozygous diplotype (b1_H12 and b1_OH12(*)) was 3,064, and the number of another homozygous diplotype (a b1_OH12(*) pair) was 3,858. The strongest association with HDL-C levels showed a P value of 9.09 × 10^-19 and was much stronger than the P value of the single SNP rs261332. Figure 2 shows the genealogical tree visualized from HapStar tool [34]. The haplotypes (AAT or GAC) were associated with higher HDL-C levels and the haplotypes (AAC or AGC) were associated with lower HDL-C levels. The association of AAC / AAT or GAC (different in bold) with HDL-C levels was statistically significant, but the association of AAC/AGC (different in bold) with HDL-C levels was not statistically significant in (Table 6, Figure 1, 2).

Diplotypes	Minor; Major	Diplotype frequency	Effect	Standard error	P value
b1_H12(AAT-or-GAC)/b1_OH12(*)	b1_H12; b1_OH12	612/3064/3858	0.041	0.0044	9.09E-19
b1_H6(AAC-or-AGC)/b1_OH6(*)	b1_OH6; b1_H6	617/3064/3853	0.040	0.0044	1.04E-18
b1_H6(AAC-or-AGC)/b1_H12(AAT-or-GAC)	b1_H12; b1_H6	612/3059/3853	0.040	0.0044	1.28E-18
b1_H9(AAT)/b1_OH9(*)	b1_H9; b1_OH9	346/2616/4572	0.036	0.0049	2.45E-12
b1_H7(AAC-or-AGC-or-GAC)/b1_OH7(*)	b1_OH7; b1_H7	350/2618/4566	0.035	0.0049	2.73E-12
b1_H7(AAC-or-AGC-or-GAC)/b1_H9(AAT)	b1_H9; b1_H7	346/2612/4566	0.036	0.0049	2.86E-12
b1_H6(AAC-or-AGC)/b1_H9(AAT)	b1_H9; b1_H6	346/2375/3853	0.037	0.0052	4.58E-12
b1_H1(AAC)/b1_H12(AAT-or-GAC)	b1_H12; b1_H1	612/1933/1498	0.039	0.0057	3.72E-11
b1_H12(AAT-or-GAC)/b1_H13(AGC)	b1_H13; b1_H12	564/1126/612	-0.050	0.0075	5.64E-11
b1_H13(AGC)/b1_OH13(*)	b1_H13; b1_OH13	564/2922/4048	-0.023	0.0047	6.72E-07

Associations of diplotypes with HDL-C levels (the top 10 results in terms of P value).

Effect sizes (mmol/l) shown are β-coefficients and measured as additive effects, which correspond to the average change in HDL-C levels when one major haplotype cluster is replaced with one minor haplotype cluster.

HL plays a key role in the interconversion between two HDL subspecies and the indirect pathway of hepatic cholesterol uptake in reverse cholesterol transport with involvement of HDL [13, 14]. HDL-C levels are used clinically to evaluate the risk of developing cardiovascular disease [21]. The two SNPs (rs1800588 and rs2070895) in the promoter region are in almost perfect linkage disequilibrium (R² = 0.97) [24], and both the promoter SNP (rs1800588) and the intronic SNP (rs261332) show strong LD (R² = 0.92). The rs2070895 SNP showed a highly significant association, with a 0.057 mmol/l increase per A allele (P = 8 × 10^-10) [24]. In the present study, the rs261332 showed a significant association (P = 3.03 × 10^-12), with a 0.035 mmol/l increase per T allele. The rs11631342 also showed a significant association (P = 5.66 × 10^-6), with a 0.037 mmol/l increase per G allele. The diplotypes b1_H12(AAT-or-GAC)/b1_OH12(*) (Table 5) showed a strong association (P = 9.09 × 10^-19), with a 0.04 mmol/l increase per (AAT-or-GAC) haplotype cluster. Transporting cholesterol from peripheral tissues to the liver, HL interferes with the interconversions of the mature HDL particles and acts on the triglyceride-rich HDL, which is able to undergo hydrolysis to form small HDL particles, which are then transported to begin anew the process of cholesterol uptake [13]. Therefore, carriers of rs261332 T-allele and rs11631342 G-allele may have high HDL-C levels. The effect of these polymorphisms on HL synthesis may be an increased susceptibility to cardiovascular disease. The rs1800588 T-allele, rs2070895 A-allele, rs261332 T-allele, and rs11631342 G-allele may be of clinical relevance, and conferring protection against cardiovascular disease, as there is increasing evidence from population studies that increasing HDL-C levels reduces the risk of cardiovascular disease [17–20]. This finding requires replication in an independent population sample.

Significant haplotype structure has been routinely investigated to identify haplotypes carrying causative mutations. Many software tools like Haploview [2] and PLINK [3] provide statistical methods for haplotype association tests, most of which focus on only a single-haplotype. Notably, Durrant et al. showed through the simulation studies and empirical data that their association analysis based on haplotype clusters had increased power over single-locus or single-haplotype tests [10]. That is, they demonstrated the superiority of haplotype cluster-based association analysis, which extended the previous method based on the single-haplotype. Nevertheless, the haplotype cluster-based analysis can have the more limitation of association analysis compared with diplotype-based analysis as if allele-based analysis can have less diversity than genotype-based analysis. Meanwhile, Luo et al. demonstrated that DTR analysis based on diplotypes outperformed other, more common association methods [11]; however, they did not apply haplotype clusters to the DTR analysis.

We developed a novel method that can yield a synergistic effect by combining the positive aspects of analyses based on haplotype clusters and diplotypes. The Diplotyper software employs Haploview tool, utilizing LD block and haplotypes with a frequency threshold. The software also uses PLINK tool to generate all possible haplotype pairs in given genotypes of samples and perform association analysis using linear or logistic regression. In addition to employing these existing software tools, we designed a new procedure to cluster haplotypes of interest from a large set of haplotypes and implemented a function to generate the patterns of all possible diplotypes from these haplotype clusters. Finally, a function to transform diplotypes into PLINK formats was implemented. All of these processes were fully automated. We tested our method by conducting an association analysis between LIPC SNPs or diplotypes and HDL-C levels. The result showed that our approach can identify more precise and distinct signals compared with single-locus tests.

This study was supported by a grant from the Korea Health technology R&D Project, Ministry of Health & Welfare, Republic of Korea (A110749), grants from Korea Centers for Disease Control and Prevention (4845-301, 4851-302, 4851-307), and the Priority Centers Program of the National Research Foundation of Korea (NRF), which is funded by the Ministry of Education, Science, and Technology (No. 2009-0093821).

The publication costs for this article were funded by the first author.

This article has been published as part of BMC Medical Genomics Volume 6 Supplement 2, 2013: Selected articles from the Second Annual Translational Bioinformatics Conference (TBC 2012). The full contents of the supplement are available online at http://www.biomedcentral.com/bmcmedgenomics/supplements/6/S2.