Assessing departure from Hardy-Weinberg equilibrium in the presence of disease association

Association between rs11614913 Polymorphism of The MiR-196-a2 Gene and Colorectal Cancer in The Presence of Departure from Hardy-Weinberg Equilibrium

Objective

Colorectal cancer (CRC) is the fourth most common and the second most lethal cancer worldwide. CRC mortality is increasing in Iran. In the current study, we aimed to investigate association between rs11614913 polymorphism of the miR-196-a2 gene and CRC.

Materials and Methods

In this case-control study, we assessed association of the rs11614913 polymorphism in 194 patients with CRC (case) and 286 healthy individuals (control). The expectation-maximization (EM) algorithm method was used to adjust deviation from Hardy-Weinberg equilibrium (HWE).

Results

There was no significant difference between genotypic frequencies of rs11614913 polymorphism in the control and case groups. Genotypic frequencies differed in the adjusted and unadjusted deviations from the HWE. Analysis of unadjusted and adjusted independent variables showed that age, sex, alcohol consumption, and drug use were statistically significant.

Conclusion

Our findings showed that rs11614913 polymorphism was not associated with CRC risk. Deviation from HWE affected the results. It is recommended to perform further studies to establish HWE. Ignoring the equilibrium can cause in consistencies in the results of studies.

Robust, flexible, and scalable tests for Hardy-Weinberg equilibrium across diverse ancestries

Traditional Hardy-Weinberg equilibrium (HWE) tests (the χ2 test and the exact test) have long been used as a metric for evaluating genotype quality, as technical artifacts leading to incorrect genotype calls often can be identified as deviations from HWE. However, in data sets composed of individuals from diverse ancestries, HWE can be violated even without genotyping error, complicating the use of HWE testing to assess genotype data quality. In this manuscript, we present the Robust Unified Test for HWE (RUTH) to test for HWE while accounting for population structure and genotype uncertainty, and to evaluate the impact of population heterogeneity and genotype uncertainty on the standard HWE tests and alternative methods using simulated and real sequence data sets. Our results demonstrate that ignoring population structure or genotype uncertainty in HWE tests can inflate false-positive rates by many orders of magnitude. Our evaluations demonstrate different tradeoffs between false positives and statistical power across the methods, with RUTH consistently among the best across all evaluations. RUTH is implemented as a practical and scalable software tool to rapidly perform HWE tests across millions of markers and hundreds of thousands of individuals while supporting standard VCF/BCF formats. RUTH is publicly available at https://www.github.com/statgen/ruth.

Association of PICALM Gene Polymorphisms with Alzheimer's Disease: Evidence from an Updated Meta-Analysis

BACKGROUND

Previous studies have examined the roles of three polymorphisms (rs3851179, rs541458, and rs592297) of the PICALM gene in susceptibility to Alzheimer's disease (AD) with inconclusive findings.

OBJECTIVE

We performed a meta-analysis to explore whether these three polymorphisms in the PICALM gene were associated with susceptibility to AD.

METHODS

Bibliographical searches were conducted in the PubMed, Embase, Web of Science, and China National Knowledge Infrastructure (CNKI) databases. Summary Odds Ratios (ORs) with 95% Confidence Intervals (CIs) were used to assess the strength of association in a random effects model. Potential sources of heterogeneity were identified by subgroup and meta-regression analyses.

RESULTS

Twenty studies (9,017 cases and 15,448 controls) on rs3851179, 12 studies (8,077 cases and 12,022 controls) on rs541458, and 4 studies (2,106 cases and 2,234 controls) on rs592297 were considered eligible for meta-analyses. For both rs3851179 and rs541458, the overall ORs were significant under all genetic models with mild heterogeneity. Compared with G carriers, A carriers of rs3851179 were associated with a decreased risk of AD (OR = 0.88; 95% CI 0.84, 0.91, P for Z-test <0.001, I2 = 0.0%). Compared with T carriers, C carriers of rs541458 were inversely associated with AD risk (OR = 0.86; 95% CI 0.81, 0.92, P for Z-test <0.001, I2 = 39.5%). No association was observed for rs592297. Subgroup and meta-regression analyses indicated that the protective effect of the rs541458 C allele was observed only among Caucasians, not among Asians (P for interaction: 0.021~<0.001).

CONCLUSION

rs3851179 and rs541458 appear to be associated with decreased AD risk. The null associations for rs592297 with AD risk need further confirmation with a larger number of participants.

A Bayesian analysis for investigating the association between rs13266634 polymorphism in SLC30A8 gene and type 2 diabetes

Purpose

The purpose of this study was to evaluate the association between rs13266634 polymorphism in SLC30A8 gene and type 2 diabetes in Iranian population, and also to provide a way for adjusting the deviation from the Hardy-Weinberg equilibrium.

Methods

This was a case-control study, the patients were selected from the East Azerbaijan province, Iran. In this study, 125 patients with type 2 diabetes (cases) and 125 healthy individuals (controls) were studied. Genotype and allele frequencies were determined in both groups, and the deviations from the Hardy-Weinberg equilibrium were assessed using Bayesian analysis.

Results

A statistically significant association was observed between rs13266634 polymorphism in SLC30A8 gene and type 2 diabetes. In genotype assessing, data analysis showed that, TT genotype play a role in diabetes type 2 risk (P = 0.001). Subjects with TT genotype had a lower risk of diabetes compared to those with CC and CT genotypes. Also, there was no significant relationship between this polymorphism and type 2 diabetes mellitus in the absence of Hardy-Weinberg equilibrium.

Conclusion

Our findings show that, rs13266634 polymorphism was associated with the type 2 diabetes risk in the population of Eastern Azerbaijan province; however, the low number of TT homozygous genotypes affected the precision of the results. Also, the deviation from HWE affected the results. It is recommended to perform further studies to establish Hardy-Weinberg equilibrium. The inconsistency in the results may be due to the ignorance of this equilibrium.

Susceptibility to leishmaniasis is affected by host SLC11A1 gene polymorphisms: a systematic review and meta-analysis

Leishmaniases are cutaneous, mucocutaneous, and visceral diseases affecting humans and domesticated animals mostly in the tropical and subtropical areas of the planet. Host genetics have been widely investigated for their role in developing various infectious diseases. The SLC11A1 gene has been reported to play a role in neutrophil function and is associated with susceptibility to infectious and inflammatory diseases such as tuberculosis or rheumatoid arthritis. In the present meta-analysis, we investigate the genetic association of SLC11A1 polymorphisms with susceptibility to leishmaniasis. Genotypes and other risk-related data were collected from 13 case-control and family-based studies (after literature search). Conventional random-effects meta-analysis was performed using STATA 13. To pool case-control and family-based data, the weighted Stouffer's method was also applied. Eight polymorphisms were investigated: rs2276631, rs3731865, rs3731864, rs17221959, rs201565523, rs2279015, rs17235409, and rs17235416. We found that rs17235409 (D543N) and rs17235416 (1729 + 55del4) are significantly associated with a risk for cutaneous leishmaniasis (CL), whereas rs17221959, rs2279015, and rs17235409 are associated with visceral leishmaniasis (VL). Our results suggest that polymorphisms in SLC11A1 affect susceptibility to CL and VL. These findings open new pathways in understanding macrophage response to Leishmania infection and the genetic factors predisposing to symptomatic CL or VL that can lead to the usage of predictive biomarkers in populations at risk.

rs2267531, a promoter SNP within glypican-3 gene in the X chromosome, is associated with hepatocellular carcinoma in Egyptians

Hepatocellular carcinoma (HCC) is a major health concern in Egypt owing to the high prevalence of hepatitis C virus (HCV) infection. HCC incidence is characterized by obvious male predominance, yet the molecular mechanisms behind this gender bias are still unidentified. Functional variations in X-linked genes have more impact on males than females. Glypican-3 (GPC3) gene, located in the Xq26 region, has lately emerged as being potentially implicated in hepatocellular carcinogenesis. The current study was designed to examine the association of -784 G/C single nucleotide polymorphism (SNP) in GPC3 promoter region (rs2267531) with HCC susceptibility in male and female Egyptian HCV patients. Our results revealed a significant association between GPC3 and HCC risk in both males and females, evidenced by higher C allele and CC/C genotype frequencies in HCC patients when compared to controls. However, no such association was found when comparing HCV patients to controls. Moreover, GPC3 gene and protein expression levels were significantly higher in CC/C than in GG/G genotype carriers in males and females. The CC/C genotype exhibited a significant shorter overall survival than GG/G genotype in HCC patients. In conclusion, GPC3 rs2267531 on the X chromosome is significantly associated with HCC, but not with HCV infection, in the Egyptian population.

The genetic component of human longevity: New insights from the analysis of pathway-based SNP-SNP interactions

In human longevity studies, single nucleotide polymorphism (SNP) analysis identified a large number of genetic variants with small effects, yet not easily replicable in different populations. New insights may come from the combined analysis of different SNPs, especially when grouped by metabolic pathway. We applied this approach to study the joint effect on longevity of SNPs belonging to three candidate pathways, the insulin/insulin‐like growth factor signalling (IIS), DNA repair and pro/antioxidant. We analysed data from 1,058 tagging SNPs in 140 genes, collected in 1825 subjects (1,089 unrelated nonagenarians from the Danish 1905 Birth Cohort Study and 736 Danish controls aged 46–55 years) for evaluating synergic interactions by SNPsyn. Synergies were further tested by the multidimensional reduction (MDR) approach, both intra‐ and interpathways. The best combinations (FDR<0.0001) resulted those encompassing IGF1R‐rs12437963 and PTPN1‐rs6067484, TP53‐rs2078486 and ERCC2‐rs50871, TXNRD1‐rs17202060 and TP53‐rs2078486, the latter two supporting a central role of TP53 in mediating the concerted activation of the DNA repair and pro‐antioxidant pathways in human longevity. Results were consistently replicated with both approaches, as well as a significant effect on longevity was found for the GHSR gene, which also interacts with partners belonging to both IIS and DNA repair pathways (PAPPA,PTPN1,PARK7, MRE11A). The combination GHSR‐MREA11, positively associated with longevity by MDR, was further found influencing longitudinal survival in nonagenarian females (p = .026). Results here presented highlight the validity of SNP‐SNP interactions analyses for investigating the genetics of human longevity, confirming previously identified markers but also pointing to novel genes as central nodes of additional networks involved in human longevity.

The interaction of miR-34b/c polymorphisms and negative life events increases susceptibility to major depressive disorder in Han Chinese population

BACKGROUND

Previous studies have shown that microRNAs(miRNAs) are involved in the pathogenesis of MDD; in particular, miR-34b/c has been implicated in MDD risk and found to exert antidepressant effects. However, the effects of miR-34b/c polymorphisms on MDD risk have not been investigated.

METHODS

In this study, we evaluated the effect of miR-34b/c gene polymorphisms and their interaction with negative life events in relation to MDD, using data from 381 Han Chinese patients with MDD and 291 healthy volunteers. Allelic, genotypic, haplotypic, and gene-environment associations were analyzed using UNPHASED and SPSS software.

RESULTS

After discarding data with extremely severe negative life events in our study population, we found an association between rs4938723, rs2187473 polymorphisms and MDD in the dominant models (TC/CC vs. TT, OR=1.45, P=0.027; TC/CC vs. TT, OR=3.32, P=0.030). In haplotype analysis, the C-G haplotype (rs4938723/rs28757623) showed the strongest association with MDD (OR=1.95, P=0.026). Additionally, we found significant gene-environment combination rs4938723 C allele, rs28757623 G allele and high level of negative life events (C-G-HN) was significantly associated with MDD (OR, 3.85; 95% CI, 1.62-9.13). In addition, the combination of (C-C-HN) is of significance (OR, 2.99; 95% CI, 1.36-6.60), indicating that the rs28757623 C allele may contribute to the risk of MDD as well.

LIMITATIONS

The sample size was small and the role of miR-34b/c polymorphisms for MDD should be assessed using independent samples from other ethnic populations.

CONCLUSIONS

Our results suggest that miR-34b/c is a susceptibility factor for MDD stratified by negative life events and that rs4938723 is a significant association locus for gene-environment interaction in relation to MDD risk.

A Case-Control Study of the Association between Polymorphisms in the Fibrinogen Alpha Chain Gene and Schizophrenia

Our previous studies using the mass spectrum analysis provided evidence that fibrinopeptide A (FPA) could be a potential biomarker for schizophrenia diagnosis. We sought further to demonstrate that variants in the fibrinogen alpha chain gene (FGA) coded FPA might confer vulnerability to schizophrenia. 1,145 patients with schizophrenia and 1,016 healthy volunteers from the Han population in Northeast China were recruited. The association of three tag single nucleotide polymorphisms (SNPs) (rs2070011 in the 5'UTR, rs2070016 in intron 4, and rs2070022 in the 3'UTR) in FGA and schizophrenia was examined using a case-control study design. Genotypic distributions of these three SNPs were not found to be significantly different between cases and controls (rs2070011: χ² = 1.28, P = 0.528; rs2070016: χ² = 4.11, P = 0.128; rs2070022: χ² = 1.23, P = 0.541). There were also no significant differences in SNP allelic frequencies between cases and controls (all P > 0.05). Additionally, the frequency of haplotypes consisting of alleles of these three SNPs was not significantly different between cases and healthy control subjects (global χ² = 9.27, P = 0.159). Our study did not show a significant association of FGA SNPs with schizophrenia. Future studies may need to test more FGA SNPs in a larger sample to identify those SNPs with a minor or moderate effect on schizophrenia.

Testing Departure from Hardy-Weinberg Proportions

The Hardy-Weinberg principle, one of the most important principles in population genetics, was originally developed for the study of allele frequency changes in a population over generations. It is now, however, widely used in studies of human diseases to detect inbreeding, population stratification, and genotyping errors. For assessment of deviation from Hardy-Weinberg proportions in data, the most popular approaches include the asymptotic Pearson's chi-squared goodness-of-fit test and the exact test. Pearson's chi-squared goodness-of-fit test is simple and straightforward, but is very sensitive to a small sample size or rare allele frequency. The exact test of Hardy-Weinberg proportions is preferable in these situations. The exact test can be performed through complete enumeration of heterozygote genotypes or on the basis of the Markov chain Monte Carlo procedure. In this chapter, we describe the Hardy-Weinberg principle and the commonly used Hardy-Weinberg proportion tests and their applications, and we demonstrate how the chi-squared test and exact test of Hardy-Weinberg proportions can be performed step-by-step using the popular software programs SAS, R, and PLINK, which have been widely used in genetic association studies, along with numerical examples. We also discuss approaches for testing Hardy-Weinberg proportions in case-control study designs that are better than traditional approaches for testing Hardy-Weinberg proportions in controls only. Finally, we note that deviation from the Hardy-Weinberg proportions in affected individuals can provide evidence for an association between genetic variants and diseases.

Ocular toxoplasmosis: susceptibility in respect to the genes encoding the KIR receptors and their HLA class I ligands

The objective of this study was to investigate the influence of the genes encoding the KIR receptors and their HLA ligands in the susceptibility of ocular toxoplasmosis. A total of 297 patients serologically-diagnosed with toxoplasmosis were selected and stratified according to the presence (n = 148) or absence (n = 149) of ocular scars/lesions due to toxoplasmosis. The group of patients with scars/lesions was further subdivided into two groups according to the type of ocular manifestation observed: primary (n = 120) or recurrent (n = 28). Genotyping was performed by PCR-SSOP. Statistical analyses were conducted using the Chi-square test, and odds ratio with a 95% confidence interval was also calculated to evaluate the risk association. The activating KIR3DS1 gene was associated with increased susceptibility for ocular toxoplasmosis. The activating KIR together with their HLA ligands (KIR3DS1-Bw4-80Ile and KIR2DS1⁺/C2⁺⁺ KIR3DS1⁺/Bw4-80Ile⁺) were associated with increased susceptibility for ocular toxoplasmosis and its clinical manifestations. KIR-HLA inhibitory pairs -KIR2DL3/2DL3-C1/C1 and KIR2DL3/2DL3-C1- were associated with decreased susceptibility for ocular toxoplasmosis and its clinical forms, while the KIR3DS1⁻/KIR3DL1⁺/Bw4-80Ile⁺ combination was associated as a protective factor against the development of ocular toxoplasmosis and, in particular, against recurrent manifestations. Our data demonstrate that activating and inhibitory KIR genes may influence the development of ocular toxoplasmosis.

The interaction of combined effects of the BDNF and PRKCG genes and negative life events in major depressive disorder

Major depressive disorder (MDD) is a mental disorder that results from complex interplay between multiple and partially overlapping sets of susceptibility genes and environmental factors. The brain derived neurotrophic factor (BDNF) and Protein kinase C gamma type (PRKCG) are logical candidate genes in MDD. Among diverse environmental factors, negative life events have been suggested to exert a crucial impact on brain development. In the present study, we hypothesized that interactions between genetic variants in BDNF and PRKCG and negative life events may play an important role in the development of MDD. We recruited a total of 406 patients with MDD and 391 age- and gender-matched control subjects. Gene-environment interactions were analyzed using generalized multifactor dimensionality reduction (GMDR). Under a dominant model, we observed a significant three-way interaction among BDNF rs6265, PRKCG rs3745406, and negative life events. The gene-environment combination of PRKCG rs3745406 C allele, BDNF rs6265 G allele and high level of negative life events (C-G-HN) was significantly associated with MDD (OR, 5.97; 95% CI, 2.71-13.15). To our knowledge, this is the first report of evidence that the BDNF-PRKCG interaction may modify the relationship between negative life events and MDD in the Chinese population.

Genetic Susceptibility to Cardiac and Digestive Clinical Forms of Chronic Chagas Disease: Involvement of the CCR5 59029 A/G Polymorphism

The clinical manifestations of chronic Chagas disease include the cardiac form of the disease and the digestive form. Not all the factors that act in the variable clinical course of this disease are known. This study investigated whether the CCR5Δ32 (rs333) and CCR5 59029 A/G (promoter region--rs1799987) polymorphisms of the CCR5 gene are associated with different clinical forms of chronic Chagas disease and with the severity of left ventricular systolic dysfunction in patients with chronic Chagas heart disease (CCHD). The antibodies anti-T. cruzi were identified by ELISA. PCR and PCR-RFLP were used to identify the CCR5Δ32 and CCR5 59029 A/G polymorphisms. The chi-square test was used to compare variables between groups. There was a higher frequency of the AA genotype in patients with CCHD compared with patients with the digestive form of the disease and the control group. The results also showed a high frequency of the AG genotype in patients with the digestive form of the disease compared to the other groups. The results of this study show that the CCR5Δ32 polymorphism does not seem to influence the different clinical manifestations of Chagas disease but there is involvement of the CCR5 59029 A/G polymorphism in susceptibility to the different forms of chronic Chagas disease. Besides, these polymorphisms do not influence left ventricular systolic dysfunction in patients with CCHD.

Polymorphism analysis in estrogen receptors alpha and beta genes and their association with infertile population in Pakistan

Studies on polymorphism of estrogen receptor (ESR) alpha and beta genes have been mostly implicated in infertility, but the results have been controversial due to lack of comprehensive data. The present study focused on association of ESR genes with both male and female infertility. In ESRα, PvuII (rs2234693) and XbaI (rs9340799) were studied while in ESRβ gene, risk of infertility was determined for silent G/A RsaI (rs1256049) polymorphism. Total 124 subjects (74 cases and 50 controls) were part of this study having primary infertility. Restriction fragment length polymorphism (RFLP) was performed with PvuII, XbaI and RsaI to determine polymorphism. Correlation between age and follicle stimulating hormone (FSH) of cases and controls was determined and no association was found between infertility and FSH hormone. Heterozygous AG genotype of XbaI polymorphism (P= 2.505e-06) and heterozygous TC genotype (P= 0.00003) in PvuII polymorphism were strongly associated with risk of infertility. In ESRβ gene, there was lack of polymorphism for RsaI in our population as all subjects were homozygous (GG). Haplotype frequencies showed that XbaI and PvuII polymorphisms are in strong linkage disequilibrium. This study shows that in our population XbaI and PvuII polymorphisms of ESRα are associated with risk of infertility.

Overrepresentation of the COL3A1 AA genotype in Polish skiers with anterior cruciate ligament injury

Although various intrinsic and extrinsic risk factors for anterior cruciate ligament (ACL) rupture have been identified, the exact aetiology of the injury is not yet fully understood. Type III collagen is an important factor in the repair of connective tissue, and certain gene polymorphisms may impair the tensile strength. The aim of this study was to examine the association of the COL3A1 rs1800255 polymorphism with ACL rupture in Polish male recreational skiers. A total of 321 male Polish recreational skiers were recruited for this study; 138 had surgically diagnosed primary ACL ruptures (ACL-injured group) and 183 were apparently healthy male skiers (control group – CON) who had no self-reported history of ligament or tendon injury. Both groups had a comparable level of exposure to ACL injury. Genomic DNA was extracted from the oral epithelial cells. All samples were genotyped on a real-time polymerase chain reaction instrument. The genotype distribution in the ACL-injured group was significantly different than in CON (respectively: AA=10.1 vs 2.2%, AG=22.5 vs 36.1, GG=67.4 vs 61.8%; p=0.0087). The AA vs AG+GG genotype of COL3A1 (odds ratio (OR)=5.05; 95% confidence interval (CI), 1.62-15.71, p=0.003) was significantly overrepresented in the ACL-injured group compared with CON. The frequency of the A allele was higher in the ACL-injured group (21.4%) compared with CON (20.2%), but the difference was not statistically significant (p=0.72). This study revealed an association between the COL3A1 rs1800255 polymorphism and ACL ruptures in Polish skiers.

Association between AKT1 gene polymorphisms and depressive symptoms in the Chinese Han population with major depressive disorder

For this study, 461 Chinese Han patients with depressive disorder were recruited. The AKT1 genotype and allele distribution were determined by PCR amplification and direct sequencing. UNPHASED software was used to analyze associations between the 17-item Hamilton Depression Rating Scale, total score, four factors and the AKT1 rs2494746 and rs3001371 polymorphisms. The results indicate that there is a significant association between suicidal ideation and anxiety symptoms in depressed patients and the rs2494746 polymorphism. The other AKT1 polymorphism, rs3001371, was significantly associated with work and activities. Patients with the rs3001371-A allele had a significantly more severe illness compared to patients with the rs3001371-G allele. Thus, AKT1 polymorphisms appear to be associated with depression severity, anxiety symptoms, work and activities, and suicide attempts in patients with depressive disorder.

Association of vWA and TPOX polymorphisms with venous thrombosis in Mexican mestizos

Objective. Venous thromboembolism (VTE) is a multifactorial disorder and, worldwide, the most important cause of morbidity and mortality. Genetic factors play a critical role in its aetiology. Microsatellites are the most important source of human genetic variation having more phenotypic effect than many single nucleotide polymorphisms. Hence, we evaluate a possible relationship between VTE and the genetic variants in von Willebrand factor, human alpha fibrinogen, and human thyroid peroxidase microsatellites to identify possible diagnostic markers. Methods. Genotypes were obtained from 177 patients with VTE and 531 nonrelated individuals using validated genotyping methods. The allelic frequencies were compared; Bayesian methods were used to correct population stratification to avoid spurious associations. Results. The vWA-18, TPOX-9, and TPOX-12 alleles were significantly associated with VTE. Moreover, subjects bearing the combination vWA-18/TPOX-12 loci exhibited doubled risk for VTE (95% CI = 1.02–3.64), whereas the combination vWA-18/TPOX-9 showed an OR = 10 (95% CI = 4.93–21.49). Conclusions. The vWA and TPOX microsatellites are good candidate biomarkers in venous thromboembolism diseases and could help to elucidate their origins. Additionally, these polymorphisms could become useful markers for genetic studies of VTE in the Mexican population; however, further studies should be done owing that this data only show preliminary evidence.

Statistical inference for Hardy-Weinberg proportions in the presence of missing genotype information

In genetic association studies, tests for Hardy-Weinberg proportions are often employed as a quality control checking procedure. Missing genotypes are typically discarded prior to testing. In this paper we show that inference for Hardy-Weinberg proportions can be biased when missing values are discarded. We propose to use multiple imputation of missing values in order to improve inference for Hardy-Weinberg proportions. For imputation we employ a multinomial logit model that uses information from allele intensities and/or neighbouring markers. Analysis of an empirical data set of single nucleotide polymorphisms possibly related to colon cancer reveals that missing genotypes are not missing completely at random. Deviation from Hardy-Weinberg proportions is mostly due to a lack of heterozygotes. Inbreeding coefficients estimated by multiple imputation of the missings are typically lowered with respect to inbreeding coefficients estimated by discarding the missings. Accounting for missings by multiple imputation qualitatively changed the results of 10 to 17% of the statistical tests performed. Estimates of inbreeding coefficients obtained by multiple imputation showed high correlation with estimates obtained by single imputation using an external reference panel. Our conclusion is that imputation of missing data leads to improved statistical inference for Hardy-Weinberg proportions.

A shrinkage method for testing the Hardy-Weinberg equilibrium in case-control studies

Testing for the Hardy-Weinberg equilibrium (HWE) is often used as an initial step for checking the quality of genotyping. When testing the HWE for case-control data, the impact of a potential genetic association between the marker and the disease must be controlled for otherwise the results may be biased. Li and Li [2008] proposed a likelihood ratio test (LRT) that accounts for this potential genetic association and it is more powerful than the commonly used control-only χ² test. However, the LRT is not efficient when the marker is independent of the disease, and also requires numerical optimization to calculate the test statistic. In this article, we propose a novel shrinkage test for assessing the HWE. The proposed shrinkage test yields higher statistical power than the LRT when the marker is independent of or weakly associated with the disease, and converges to the LRT when the marker is strongly associated with the disease. In addition, the proposed shrinkage test has a closed form and can be easily used to test the HWE for large datasets that result from genome-wide association studies. We compare the performance of the shrinkage test with existing methods using simulation studies, and apply the shrinkage test to a genome-wide association dataset for Alzheimer's disease.

A multistage genetic association study identifies breast cancer risk loci at 10q25 and 16q24

BACKGROUND

Heritable risk for breast cancer includes an increasing number of common, low effect risk variants. We conducted a multistage genetic association study in a series of independent epidemiologic breast cancer study populations to identify novel breast cancer risk variants.

METHODS

We tested 1,162 SNPs of greatest nominal significance from stage I of the Cancer Genetic Markers of Susceptibility breast cancer study (CGEMS; 1,145 cases, 1,142 controls) for evidence of replicated association with breast cancer in the Nashville Breast Cohort (NBC; 599 cases, 1,161 controls), the Collaborative Breast Cancer Study (CBCS; 1,552 cases, 1,185 controls), and BioVU Breast Cancer Study (BioVU; 1,172 cases, 1,172 controls).

RESULTS

Among these SNPs, a series of validated breast cancer risk variants yielded expected associations in the study populations. In addition, we observed two previously unreported loci that were significantly associated with breast cancer risk in the CGEMS, NBC, and CBCS study populations and had a consistent, although not statistically significant, risk effect in the BioVU study population. These were rs1626678 at 10q25.3 near ENO4 and KIAA1598 (meta-analysis age-adjusted OR = 1.13 [1.07-1.20], P = 5.6 × 10(-5)), and rs8046508 at 16q23.1 in the eighth intron of WWOX (meta-analysis age-adjusted OR = 1.20 [1.10-1.31], P = 3.5 × 10(-5)).

CONCLUSIONS

Our data supports the association of two novel loci, at 10q25.3 and 16q23.1, with risk of breast cancer.

IMPACT

The expanding compendium of known breast cancer genetic risk variants holds increasing power for clinical risk prediction models of breast cancer, improving upon the Gail model.

Testing departure from Hardy-Weinberg proportions

The Hardy-Weinberg principle, one of the most important principles in population genetics, was originally developed for the study of allele frequency changes in a population over generations. It is now, however, widely used in studies of human diseases to detect inbreeding, populations stratification, and genotyping errors. For assessment of deviation from the Hardy-Weinberg proportions in data, the most popular approaches include the asymptotic Pearson's chi-square goodness-of-fit test and the exact test. The Pearson's chi-square goodness-of-fit test is simple and straightforward, but it is very sensitive to small sample size or rare allele frequency. The exact test of Hardy-Weinberg proportions is preferable in these situations. The exact test can be performed through complete enumeration of heterozygote genotypes or on the basis of the Markov chain Monte Carlo procedure. In this chapter, we describe the Hardy-Weinberg principle and the commonly used Hardy-Weinberg proportions tests and their applications, and we demonstrate how the chi-square test and exact test of Hardy-Weinberg proportions can be performed step-by-step using the popular software programs SAS, R, and PLINK, which have been widely used in genetic association studies, along with numerical examples. We also discuss recent approaches for testing Hardy-Weinberg proportions in case-control study designs that are better than traditional approaches for testing Hardy-Weinberg proportions in controls only. Finally, we note that deviation from the Hardy-Weinberg proportions in affected individuals can provide evidence for an association between genetic variants and diseases.

A flexible likelihood framework for detecting associations with secondary phenotypes in genetic studies using selected samples: application to sequence data

For most complex trait association studies using next-generation sequencing, in addition to the primary phenotype of interest, many clinically important secondary traits are also available, which can be analyzed to map susceptibility genes. Owing to high sequencing costs, most studies use selected samples, and the sampling mechanisms of these studies can be complicated. When the primary and secondary traits are correlated, analyses of secondary phenotypes can cause spurious associations in selected samples and existing methods are inadequate to adjust for them. To address this problem, a likelihood-based method, MULTI-TRAIT-ASSOCIATION (MTA) was developed. MTA is flexible and can be applied to any study with known sampling mechanisms. It also allows efficient inferences of genetic parameters. To investigate the power of MTA and different study designs, extensive simulations were performed under rigorous population genetic and phenotypic models. It is demonstrated that there are great benefits for analyzing secondary phenotypes in selected samples. In particular, using case-control samples and samples with extreme primary phenotypes can be more powerful than analyzing random samples of equivalent size. One major challenge for sequence-based association studies is that most data sets are not of sufficient size to be adequately powered. By applying MTA, data sets ascertained under distinct mechanisms or targeted at different primary traits can be jointly analyzed to map common phenotypes and greatly increase power. The combined analysis can be performed using freely available data sets from public repositories, for example, dbGaP. In conclusion, MTA will have an important role in dissecting the etiology of complex traits.

Testing Hardy-Weinberg proportions in a frequency-matched case-control genetic association study

In case-control genetic association studies, cases are subjects with the disease and controls are subjects without the disease. At the time of case-control data collection, information about secondary phenotypes is also collected. In addition to studies of primary diseases, there has been some interest in studying genetic variants associated with secondary phenotypes. In genetic association studies, the deviation from Hardy-Weinberg proportion (HWP) of each genetic marker is assessed as an initial quality check to identify questionable genotypes. Generally, HWP tests are performed based on the controls for the primary disease or secondary phenotype. However, when the disease or phenotype of interest is common, the controls do not represent the general population. Therefore, using only controls for testing HWP can result in a highly inflated type I error rate for the disease- and/or phenotype-associated variants. Recently, two approaches, the likelihood ratio test (LRT) approach and the mixture HWP (mHWP) exact test were proposed for testing HWP in samples from case-control studies. Here, we show that these two approaches result in inflated type I error rates and could lead to the removal from further analysis of potential causal genetic variants associated with the primary disease and/or secondary phenotype when the study of primary disease is frequency-matched on the secondary phenotype. Therefore, we proposed alternative approaches, which extend the LRT and mHWP approaches, for assessing HWP that account for frequency matching. The goal was to maintain more (possible causative) single-nucleotide polymorphisms in the sample for further analysis. Our simulation results showed that both extended approaches could control type I error probabilities. We also applied the proposed approaches to test HWP for SNPs from a genome-wide association study of lung cancer that was frequency-matched on smoking status and found that the proposed approaches can keep more genetic variants for association studies.

Approximate and exact tests of Hardy-Weinberg equilibrium using uncertain genotypes

Testing for Hardy-Weinberg equilibrium (HWE) is commonly used as a quality control filter in genome-wide scans for markers with experimentally determined genotypes. In contrast, for markers with imputed genotypes, there are post-imputation metrics of quality that can be used as screens but there are no formal tests of deviation from HWE. Similarly, there are no formal tests of deviation from HWE for probabilistic genotypes that are generated by sequencing projects. Here, I describe generalizations of the approximate χ(2) and exact tests of HWE for use with uncertain genotypes. The tests fully account for the probabilities of all possible genotypes at a marker for each individual. By computer simulation, the approximate and exact tests are shown to maintain valid control of the type I error rate. Calculations of the loss of power as the uncertainty in genotypes increases are illustrated. The tests are compatible with chip-based genotypes for single-nucleotide polymorphisms and copy number polymorphisms, imputed genotypes, and probabilistic assignments of genotype from variable-coverage sequence data.

A test of Hardy-Weinberg equilibrium in structured populations

Testing for Hardy-Weinberg equilibrium (HWE) is used routinely as an important initial step for genotype data quality checking. Departure from HWE can be caused by many factors, such as genotyping errors, population stratification, and disease association, if we use affected individuals only. In a structured population, even if a marker is in HWE in each subpopulation, data may show departure from HWE if allele frequencies are different in different subpopulations and such a departure can be misinterpreted as a potential problem in genotyping quality, resulting in false exclusion from future analysis. In this article, we propose a new HWE test, a test for HWE in structured populations (HWES) that can assess departure from HWE and take into account of population stratification at the same time. Our proposed test can distinguish departure from HWE caused by population stratification and departure from HWE caused by other factors. We use simulation studies as well as applications to real data sets to evaluate the performance of the proposed test. Results show that, for a wide range of population structures, our proposed test has correct type I error rates while the traditional χ(2) test will lead to false-positive results. In homogenous populations, our proposed test has comparable power with the traditional χ(2) test.

Quantitative trait analysis suggests polymorphisms of estrogen-related genes regulate human sperm concentrations and motility

BACKGROUND

Human spermatogenesis is regulated by complex networks, and estrogens are recognized as one of the significant regulators of spermatogenesis. We tested the associations between variants of estrogen-related genes and semen parameters.

METHODS

We performed genotyping for genetic variants of estrogen-related genes and quantitative trait analysis of fertile and infertile men with well-characterized reproductive phenotypes. Men with known semen parameters (n= 677) were enrolled, including 210 fertile men and 467 infertile men. A total of 17 genetic markers from 10 genes, including 2 estrogen receptors (ER-α, ER-β), 7 estrogen synthesizing/metabolizing genes (CYP19A1, HSD17B1, CYP1A1, CYP1B1, COMT, GSTM1, GSTT1) and 1 transport gene (SHBG) were genotyped. Sperm concentration, motility and morphology were taken as quantitative traits to correlate with genetic variants in the estrogen-related genes.

RESULTS

Five genes (rs1801132 and rs2228480 of the ER-α gene, rs1256049 and rs4986938 of the ER-β gene, rs605059 of the HSD17B1 gene, rs1799941 of the SHBG gene and rs1048943 and rs4646903 of the CYP1A1 gene) were found to be significantly associated with sperm concentration (P< 0.01), while five genes (rs1801132 of the ER-a gene, rs1256049 of the ER-β gene, rs1048943 of the CYP1A1 gene, rs605059 of the HSD17B1 gene and rs1799941 along with rs6259 of the SHBG gene) were associated with sperm motility (P< 0.01). None of the estrogen-related genes were associated with sperm morphology. With an increasing number of risk alleles, sperm concentration and motility tended to deteriorate and show a loci-dosage effect.

CONCLUSIONS

Quantitative trait analysis based on a limited number of genetic markers suggests that estrogen-related genes mainly regulate sperm concentration and motility.

How to link call rate and p-values for Hardy-Weinberg equilibrium as measures of genome-wide SNP data quality

We study the link between two quality measures of SNP (single nucleotide polymorphism) data in genome-wide association (GWA) studies, that is, per SNP call rates (CR) and p-values for testing Hardy-Weinberg equilibrium (HWE). The aim is to improve these measures by applying methods based on realized randomized p-values, the false discovery rate and estimates for the proportion of false hypotheses. While exact non-randomized conditional p-values for testing HWE cannot be recommended for estimating the proportion of false hypotheses, their realized randomized counterparts should be used. P-values corresponding to the asymptotic unconditional chi-square test lead to reasonable estimates only if SNPs with low minor allele frequency are excluded. We provide an algorithm to compute the probability that SNPs violate HWE given the observed CR, which yields an improved measure of data quality. The proposed methods are applied to SNP data from the KORA (Cooperative Health Research in the Region of Augsburg, Southern Germany) 500 K project, a GWA study in a population-based sample genotyped by Affymetrix GeneChip 500 K arrays using the calling algorithm BRLMM 1.4.0. We show that all SNPs with CR = 100 per cent are nearly in perfect HWE which militates in favor of the population to meet the conditions required for HWE at least for these SNPs. Moreover, we show that the proportion of SNPs not being in HWE increases with decreasing CR. We conclude that using a single threshold for judging HWE p-values without taking the CR into account is problematic. Instead we recommend a stratified analysis with respect to CR.

Genetic association between single nucleotide polymorphisms in the paraoxonase 1 (PON1) gene and small-for-gestational-age birth in related and unrelated subjects

Paraoxonase 1 (PON1) protects against oxidative modification of low density lipoproteins. The PON1 gene has 3 polymorphisms considered strong determinants of PON1 levels: Q192R and L55M in the coding region and C-108T in the promoter region. PON1 levels are also influenced by smoking. The authors hypothesized that PON1 variants could increase the risk of vascular thrombosis, leading in turn to placental insufficiency and small-for-gestational-age (SGA) birth. The author compared PON1 variants and haplotypes between 448 newborn SGA cases and 431 newborn controls from Montréal, Québec, Canada (1998-2000) and studied the effects of interaction with maternal smoking. Transmission of the variants in case-parent trios was used as validation of the case-control results; the authors combined case-control and family data to analyze the associations of variants with SGA birth. In the case-control analysis, the T allele from C-108T increased the risk of SGA birth (additive odds ratio = 1.30, 95% confidence interval (CI): 1.06, 1.59), and the TRL haplotype (T from C-108T, R from Q192R, and L from L55M) was associated with an odds ratio of 1.51 (95% CI: 1.07, 2.15). Among smokers, the CRL haplotype was protective (odds ratio = 0.48, 95% CI: 0.28, 0.82). Case-parent trio results were compatible with case-control results.

Using both cases and controls for testing hardy-weinberg proportions in a genetic association study

OBJECTIVES

Assessment of the Hardy-Weinberg proportion (HWP) in controls has been widely used as a quality control measure in case-control association studies. However, when the disease being studied is common, controls might not represent the general population, which could result in inaccurate HWP test results. Such results could lead investigators to discard important single-nucleotide polymorphisms (SNPs) that could potentially be causal. In this paper, we showed the inappropriateness of the HWP test in controls and proposed a mixture HWP (mHWP) exact test using a mixture sample that mimics the general population.

METHODS

The mHWP exact test estimates HWP in a mixture sample that is a combination of both cases and controls proportional to the prevalence of disease. We implemented a re-sampling procedure to construct mixture samples and then obtained the empirical p value of HWP in the general population. Simulation studies were performed to investigate the performance of the proposed mHWP exact test. The method was also applied to a genetic association study of obesity.

RESULTS

The results showed that the mHWP exact test is more likely than either the traditional HWP method in controls or the likelihood-based approach to keep causal SNPs for further analysis when the disease is more common.

CONCLUSION

The mHWP exact test using a mixture sample is a better HWP test for case-control genetic association studies than the traditional HWP in controls or the likelihood-based approach, and it will improve our ability to keep causal SNPs in the case-control genetic association studies.

Associations of polymorphisms of eight muscle- or metabolism-related genes with performance in Mount Olympus marathon runners

Athletic endurance performance is probably partly under genetic control, but genetic association studies have yielded inconclusive results. The objective of the present study was to evaluate the association of polymorphisms in eight muscle- or metabolism-related genes with endurance performance in participants of the Olympus Marathon running race. We recruited 438 athletes who participated in the 2007 and 2008 annual running events of the Olympus Marathon: a 43.8-km race with an ascent from sea level to 2,690-m altitude and then a descent to 300 m. Phenotypes of interest were the competitive event time at the specific Olympus Marathon where the athlete was enrolled, the fastest reported timing ever achieved in an Olympus Marathon, and how many kilometers per week the athlete ran during the previous year. Eleven polymorphisms in alpha(3)-actinin (ACTN3), AMP deaminase-1 (AMPD1), bradykinin B(2) receptor (BDKRB2), beta(2)-adrenergic receptor (ADRB2), peroxisome proliferator-activated receptor (PPAR)-gamma coactivator-1 alpha (PPARGC1A), PPAR-alpha (PPARA), PPAR-delta (PPARD), and apoliprotein E (APOE) were evaluated. Hardy-Weinberg equilibrium testing on the overall cohort of male athletes showed a significant deviation for BDKRB2 rs1799722 (P = 0.018; P = 0.006 when limited to 316 habitual male runners) with an excess of the TT genotype. Across all athletes, no associations showed nominal statistical significance for any of the three phenotypes, and the same was true when analyses were limited to men (n = 417). When limited to 316 male athletes who identified running as their preferred sport, ADRB2 rs1042713 had nominally significant associations with faster times for the minor (A) allele for the fastest time ever (P = 0.01). The direction of effect was identical as previously postulated only for BDKRB2 rs1799722 and ADRB2 rs1042713, indicating consistency. BDKRB2 rs1799722 and ADRB2 rs1042713 have some support for being implicated in endurance performance among habitual runners and require further investigation.

Interleukin-2 receptor beta gene C627T polymorphism in Korean women with endometriosis: a case-control study

BACKGROUND

The purpose of this study was to investigate the potential association of the C627T polymorphism in the interleukin-2 receptor beta gene (IL-2R beta) with the risk of endometriosis in Korean women.

METHODS

Two hundred and thirty-seven women with surgically or histologically diagnosed endometriosis of stages III and IV were recruited for this study, and 164 patients with no evidence of endometriosis diagnosed by laparoscopy or laparotomy served as controls. The C627T polymorphism of the IL-2R beta was assessed using the TaqMan allelic discrimination assay. Chi2 analysis was used to examine any differences in genotype distributions and allele frequencies of the IL-2R beta C627T polymorphism between the endometriosis cases and the controls.

RESULTS

There was no statistically significant difference in the frequency of the IL-2R beta C627T polymorphism between the endometriosis patients and the controls (28.7% C/C, 48.1% C/T and 23.2% T/T versus 29.3, 44.5 and 26.2%, respectively, P = 0.72) or in the T allele frequencies (47.3 versus 48.5%, respectively, P = 0.73). Even when the endometriosis cases were subdivided into stages III and IV, no statistically significant differences in genotype distributions or allele frequencies were observed among the three groups.

CONCLUSIONS

Contrary to the recent data reported in a Taiwanese population, our results suggest that the C627T polymorphism of the IL-2R beta gene may not be associated with the risk of endometriosis in the Korean population.

Testing Hardy-Weinberg equilibrium using family data from complex surveys

Genetic data collected during the second phase of the Third National Health and Nutrition Examination Survey (NHANES III) enable us to investigate the association of a wide variety of health factors with regard to genetic variation. The classic question when looking into the genetic variations in a population is whether the population is in the state of Hardy-Weinberg Equilibrium (HWE). Our objective was to develop test procedures using family data from complex surveys such as NHANES III. We developed six Pearson chi(2) based tests for a diallelic locus of autosomal genes. The finite sample properties of the proposed test procedures were evaluated via Monte Carlo simulation studies and the Rao-Scott first order corrected test was recommended. Test procedures were applied to three loci from NHANES III genetic databases, i.e., ADRB2, TGFB1, and VDR. HWE was shown to hold at 0.05 level for all three loci when only families with genotypic information available for two parents and for one or more children were used in the analysis.

Recovering unused information in genome-wide association studies: the benefit of analyzing SNPs out of Hardy-Weinberg equilibrium

Although the rapid advancements in high throughput genotyping technology have made genome-wide association studies possible, these studies remain an expensive undertaking, especially when considering the large sample sizes necessary to find the small to moderate effect sizes that define complex diseases. It is therefore prudent to utilize all possible information contained in a genome-wide scan. We propose a straightforward analytical approach that tests often unused SNP data without sacrificing statistical validity. We simulate genotype miscalls under a variety of models consistent with observed miscall rates and test for departures from HWE using the standard Pearson's chi(2)-test. We find that true disease susceptibility loci subjected to various patterns of genotype miscalls can be largely out of HWE and, thus, be candidates for removal before association testing. These loci, we demonstrate, can maintain sufficient statistical power even under extreme error models. We additionally show that random miscalls of null SNPs, independent of the phenotype, do not induce bias in case-control or cohort studies, and we suggest that a significant HWE test should not prevent a SNP from being tested when conducting genome-wide association studies in these scenarios. However, association findings for SNPs that are out of HWE must be treated more carefully than 'regular' findings, for example, by re-genotyping the SNP in the same study using a different genotyping technology.

Associations among GABRG1, level of response to alcohol, and drinking behaviors

BACKGROUND

Recent studies of the genetics of alcoholism have focused on a cluster of genes encoding for gamma-aminobutyric acid (GABA(A)) receptor subunits, which is thought to play a role in the expression of addiction phenotypes. This study examined allelic associations between 2 single nucleotide polymorphisms (SNPs) of the GABRG1 gene (rs1391166 and rs1497571) and alcohol phenotypes, namely level of response to alcohol, alcohol use patterns, and alcohol-related problems.

METHOD

Participants were non-treatment-seeking seeking hazardous drinkers (n = 124) who provided DNA samples, participated in a face-to-face interview for level of response to alcohol, and completed a series of drinking and individual differences measures.

RESULTS

Analyses revealed that a SNP of the GABRG1 gene (rs1497571) was associated with level of response to alcohol and drinking patterns in this subclinical sample. Follow-up mediational analyses were also conducted to examine putative mechanisms underlying these associations.

DISCUSSION

These findings replicate and extend recent research suggesting that genetic variation at the GABRG1 locus may underlie the expression of alcohol phenotypes, including level of response to alcohol.

Genetic associations with schizophrenia: meta-analyses of 12 candidate genes

Genetic association studies on schizophrenia (SZ) have been repeatedly performed over the last two decades, resulting in a consensus that results are generally inconsistent. This consensus has begun to change as a result of meta-analyses (e.g., [Glatt, S.J. and Jonsson, E.G., 2006. The Cys allele of the DRD2 Ser311Cys polymorphism has a dominant effect on risk for schizophrenia: evidence from fixed- and random-effects meta-analyses. Am. J. Med. Genet. B. Neuropsychiatr. Genet. 141, 149-154.]). The SchizophreniaGene database (http://www.schizophreniaforum.org/res/sczgene/default.asp) has been a leader in meta-analyses of SZ association data, by dynamically and comprehensively cataloging all public genetic association studies, and preparing meta-analyses of case-control data. There are 19 "top" candidate genes from these analyses (access on December 20, 2007), showing the highest effect sizes and nominally significant associations of at least one variant in the meta-analyses of all ethnic samples or of samples of Caucasian ancestry. We selected 40 polymorphisms in 12 selected "top" genes for additional meta-analyses, which had at least one familial association data. We found gene-wide (correction for the number of meta-analyses for each gene) significant allelic association evidence for seven genes in the combined samples. The odds ratios (ORs) of the associated minor risk alleles range from 1.072 to 1.121, for DRD4, MTHFR, PPP3CC and TP53. For protective allele associations, the ORs are between 0.842 and 0.886, for DAO, IL1B, and SLC6A4. In population-based sub-analyses, we found significant results in four genes in Asians (ORs between 1.084 and 1.309 for DRD4, GABRB2, PPP3CC, and TP53), and one gene in European (OR of 0.888 for SLC6A4). The association of rs1816072 of GABRB2 with SZ in Asians was significant (adjusted P=0.048 after correction for 80 tests). No significant heterogeneity between case-control and family-based study designs was detected in 35 out of 40 polymorphisms. Our results further support eight potential SZ candidate genes and suggest that family data can reasonably be included in the meta-analysis of genetic associations.