Pleiotropy among Common Genetic Loci Identified for Cardiometabolic Disorders and C-Reactive Protein

Pleiotropic genetic variants have independent effects on different phenotypes. C-reactive protein (CRP) is associated with several cardiometabolic phenotypes. Shared genetic backgrounds may partially underlie these associations. We conducted a genome-wide analysis to identify the shared genetic background of inflammation and cardiometabolic phenotypes using published genome-wide association studies (GWAS). We also evaluated whether the pleiotropic effects of such loci were biological or mediated in nature. First, we examined whether 283 common variants identified for 10 cardiometabolic phenotypes in GWAS are associated with CRP level. Second, we tested whether 18 variants identified for serum CRP are associated with 10 cardiometabolic phenotypes. We used a Bonferroni corrected p-value of 1.1×10^-04 (0.05/463) as a threshold of significance. We evaluated the independent pleiotropic effect on both phenotypes using individual level data from the Women Genome Health Study. Evaluating the genetic overlap between inflammation and cardiometabolic phenotypes, we found 13 pleiotropic regions. Additional analyses showed that 6 regions (APOC1, HNF1A, IL6R, PPP1R3B, HNF4A and IL1F10) appeared to have a pleiotropic effect on CRP independent of the effects on the cardiometabolic phenotypes. These included loci where individuals carrying the risk allele for CRP encounter higher lipid levels and risk of type 2 diabetes. In addition, 5 regions (GCKR, PABPC4, BCL7B, FTO and TMEM18) had an effect on CRP largely mediated through the cardiometabolic phenotypes. In conclusion, our results show genetic pleiotropy among inflammation and cardiometabolic phenotypes. In addition to reverse causation, our data suggests that pleiotropic genetic variants partially underlie the association between CRP and cardiometabolic phenotypes.     

An integrative functional genomics approach for discovering biomarkers in schizophrenia

Schizophrenia (SZ) is a complex disorder resulting from both genetic and environmental causes with a lifetime prevalence world-wide of 1%; however, there are no specific, sensitive and validated biomarkers for SZ. A general unifying hypothesis has been put forward that disease-associated single nucleotide polymorphisms (SNPs) from genome-wide association study (GWAS) are more likely to be associated with gene expression quantitative trait loci (eQTL). We will describe this hypothesis and review primary methodology with refinements for testing this paradigmatic approach in SZ. We will describe biomarker studies of SZ and testing enrichment of SNPs that are associated both with eQTLs and existing GWAS of SZ. SZ-associated SNPs that overlap with eQTLs can be placed into gene-gene expression, protein-protein and protein-DNA interaction networks. Further, those networks can be tested by reducing/silencing the gene expression levels of critical nodes. We present pilot data to support these methods of investigation such as the use of eQTLs to annotate GWASs of SZ, which could be applied to the field of biomarker discovery. Those networks that have association with SNP markers, especially cis-regulated expression, might lead to a more clear understanding of important candidate genes that predispose to disease and alter expression. This method has general application to many complex disorders.     

Validation of Type 2 Diabetes Risk Variants Identified by Genome-Wide Association Studies in Han Chinese Population: A Replication Study and Meta-Analysis

Background

Several genome-wide association studies (GWAS) involving European populations have successfully identified risk genetic variants associated with type 2 diabetes mellitus (T2DM). However, the effects conferred by these variants in Han Chinese population have not yet been fully elucidated.

Methods

We analyzed the effects of 24 risk genetic variants with reported associations from European GWAS in 3,040 Han Chinese subjects in Taiwan (including 1,520 T2DM cases and 1,520 controls). The discriminative power of the prediction models with and without genotype scores was compared. We further meta-analyzed the association of these variants with T2DM by pooling all candidate-gene association studies conducted in Han Chinese.

Results

Five risk variants in IGF2BP2 (rs4402960, rs1470579), CDKAL1 (rs10946398), SLC30A8 (rs13266634), and HHEX (rs1111875) genes were nominally associated with T2DM in our samples. The odds ratio was 2.22 (95% confidence interval, 1.81-2.73, P<0.0001) for subjects with the highest genetic score quartile (score>34) as compared with subjects with the lowest quartile (score<29). The incoporation of genotype score into the predictive model increased the C-statistics from 0.627 to 0.657 (P<0.0001). These estimates are very close to those observed in European populations. Gene-environment interaction analysis showed a significant interaction between rs13266634 in SLC30A8 gene and age on T2DM risk (P<0.0001). Further meta-analysis pooling 20 studies in Han Chinese confirmed the association of 10 genetic variants in IGF2BP2, CDKAL1, JAZF1, SCL30A8, HHEX, TCF7L2, EXT2, and FTO genes with T2DM. The effect sizes conferred by these risk variants in Han Chinese were similar to those observed in Europeans but the allele frequencies differ substantially between two populations.

Conclusion

We confirmed the association of 10 variants identified by European GWAS with T2DM in Han Chinese population. The incorporation of genotype scores into the prediction model led to a small but significant improvement in T2DM prediction.     

The association of C-reactive protein and CRP genotype with coronary heart disease: findings from five studies with 4,610 cases amongst 18,637 participants

Background

It is unclear whether C-reactive protein (CRP) is causally related to coronary heart disease (CHD). Genetic variants that are known to be associated with CRP levels can be used to provide causal inference of the effect of CRP on CHD. Our objective was to examine the association between CRP genetic variant +1444C>T (rs1130864) and CHD risk in the largest study to date of this association.

Methods and Results

We estimated the association of CRP genetic variant +1444C>T (rs1130864) with CRP levels and with CHD in five studies and then pooled these analyses (N = 18,637 participants amongst whom there were 4,610 cases). CRP was associated with potential confounding factors (socioeconomic position, physical activity, smoking and body mass) whereas genotype (rs1130864) was not associated with these confounders. The pooled odds ratio of CHD per doubling of circulating CRP level after adjustment for age and sex was 1.13 (95%CI: 1.06, 1.21), and after further adjustment for confounding factors it was 1.07 (95%CI: 1.02, 1.13). Genotype (rs1130864) was associated with circulating CRP; the pooled ratio of geometric means of CRP level among individuals with the TT genotype compared to those with the CT/CC genotype was 1.21 (95%CI: 1.15, 1.28) and the pooled ratio of geometric means of CRP level per additional T allele was 1.14 (95%CI: 1.11, 1.18), with no strong evidence in either analyses of between study heterogeneity (I² = 0%, p>0.9 for both analyses). There was no association of genotype (rs1130864) with CHD: pooled odds ratio 1.01 (95%CI: 0.88, 1.16) comparing individuals with TT genotype to those with CT/CC genotype and 0.96 (95%CI: 0.90, 1.03) per additional T allele (I²<7.5%, p>0.6 for both meta-analyses). An instrumental variables analysis (in which the proportion of CRP levels explained by rs1130864 was related to CHD) suggested that circulating CRP was not associated with CHD: the odds ratio for a doubling of CRP level was 1.04 (95%CI: 0.61, 1.80).

Conclusions

We found no association of a genetic variant, which is known to be related to CRP levels, (rs1130864) and having CHD. These findings do not support a causal association between circulating CRP and CHD risk, but very large, extended, genetic association studies would be required to rule this out.     

Physical,Behavioural and Genetic Predictors of Adult Hypertension: The Findings of the Kaunas Cardiovascular Risk Cohort Study

Background

The roots of adult hypertension go back to childhood. This study aimed to examine the independent effects of physical, behavioural and genetic factors identified in childhood and mid-adulthood for prediction of adult hypertension.

Methods

The study subjects were participants of the Kaunas Cardiovascular Risk Cohort study started in 1977 (n = 1082, age 12–13 years). In 2012, a total of 507 individuals (63.9% of eligible sample) participated in the 35-year follow-up survey. Health examination involved measurements of blood pressure (BP), anthropometric parameters, and interview about health behaviours. Subjects were genotyped for AGT (M235T), ACE (I/D, rs4340), ADM (rs7129220), and CACNB2 (rs12258967) genes polymorphisms. A genetic risk score was calculated as the sum of the number of risk alleles at each of four single nucleotide polymorphisms.

Results

AGT TT genotype male carriers had the highest mean values of systolic BP in childhood. In females, ADM genotype AA was associated with the highest values of systolic and diastolic BP, while CACNB2 genotype CC carriers had the highest values of diastolic BP in childhood. Systolic and diastolic BP in childhood, gain in BMI from childhood to adulthood, and risky alcohol consumption predicted hypertension in middle-aged men. In women, genetic risk score together with diastolic BP in childhood and gain in BMI were significant predictors of adult hypertension. The comparison of four nested logistic regression models showed that the prediction of hypertension improved significantly after the addition of BMI gain. Genetic risk score had a relatively weak effect on the improvement of the model performance in women.

Conclusions

BP in childhood and the gain in BMI from childhood to adulthood were significant predictors of adult hypertension in both genders. Genetic risk score in women and risky alcohol consumption in men were independently related with the risk of adult hypertension.     

The Calcitonin Receptor Gene Is a Candidate for Regulation of Susceptibility to Herpes simplex Type 1 Neuronal Infection Leading to Encephalitis in Rat

Herpes simplex encephalitis (HSE) is a fatal infection of the central nervous system (CNS) predominantly caused by Herpes simplex virus type 1. Factors regulating the susceptibility to HSE are still largely unknown. To identify host gene(s) regulating HSE susceptibility we performed a genome-wide linkage scan in an intercross between the susceptible DA and the resistant PVG rat. We found one major quantitative trait locus (QTL), Hse1, on rat chromosome 4 (confidence interval 24.3–31 Mb; LOD score 29.5) governing disease susceptibility. Fine mapping of Hse1 using recombinants, haplotype mapping and sequencing, as well as expression analysis of all genes in the interval identified the calcitonin receptor gene (Calcr) as the main candidate, which also is supported by functional studies. Thus, using unbiased genetic approach variability in Calcr was identified as potentially critical for infection and viral spread to the CNS and subsequent HSE development.     

A Genome-wide In Vitro Bacterial-Infection Screen Reveals Human Variation in the Host Response Associated with Inflammatory Disease

Recent progress in cataloguing common genetic variation has made possible genome-wide studies that are beginning to elucidate the causes and consequences of our genetic differences. Approaches that provide a mechanistic understanding of how genetic variants function to alter disease susceptibility and why they were substrates of natural selection would complement other approaches to human-genome analysis. Here we use a novel cell-based screen of bacterial infection to identify human variation in Salmonella-induced cell death. A loss-of-function allele of CARD8, a reported inhibitor of the proinflammatory protease caspase-1, was associated with increased cell death in vitro (p = 0.013). The validity of this association was demonstrated through overexpression of alternative alleles and RNA interference in cells of varying genotype. Comparison of mammalian CARD8 orthologs and examination of variation among different human populations suggest that the increase in infectious-disease burden associated with larger animal groups (i.e., herds and colonies), and possibly human population expansion, may have naturally selected for loss of CARD8. We also find that the loss-of-function CARD8 allele shows a modest association with an increased risk of systemic inflammatory response syndrome in a small study (p = 0.05). Therefore, a by-product of the selected benefit of loss of CARD8 could be increased inflammatory diseases. These results demonstrate the utility of genome-wide cell-based association screens with microbes in the identification of naturally selected variants that can impact human health.     

Genetic Variations in the Transforming Growth Factor Beta Pathway as Predictors of Bladder Cancer Risk

Bladder cancer is the fifth most common cancer in the United States, and identifying genetic markers that may predict susceptibility in high-risk population is always needed. The purpose of our study is to determine whether genetic variations in the transforming growth factor-beta (TGF-β) pathway are associated with bladder cancer risk. We identified 356 single-nucleotide polymorphisms (SNPs) in 37 key genes from this pathway and evaluated their association with cancer risk in 801 cases and 801 controls. Forty-one SNPs were significantly associated with cancer risk, and after adjusting for multiple comparisons, 9 remained significant (Q-value ≤0.1). Haplotype analysis further revealed three haplotypes within VEGFC and two haplotypes in EGFR were significantly associated with increased bladder cancer risk compared to the most common haplotype. Classification and regression tree analysis further revealed potential high-order gene-gene interactions, with VEGFC: rs3775194 being the initial split, which suggests that this variant is responsible for the most variation in risk. Individuals carrying the common genotype for VEGFC: rs3775194 and EGFR: rs7799627 and the variant genotype for VEGFR: rs4557213 had a 4.22-fold increase in risk, a much larger effect magnitude than that conferred by common genotype for VEGFR: rs4557213. Our study provides the first epidemiological evidence supporting a connection between TGF-β pathway variants and bladder cancer risk.     

Leukocyte Telomere Length-Related rs621559 and rs398652 Genetic Variants Influence Risk of HBV-Related Hepatocellular Carcinoma

Recent genome-wide association studies (GWAS) have identified eleven leukocyte telomere length (LTL)-related single nucleotide polymorphisms (SNPs). Since LTL has been associated with risk of many malignancies, LTL-related SNPs may contribute to cancer susceptibility. To test this hypothesis in hepatitis B virus (HBV)-related hepatocellular carcinoma (HCC), we genotyped these eleven LTL-related SNPs in a case-control set including 1186 HBV-related HCC cases, 508 chronic HBV carriers and 1308 healthy controls at the discovery stage. The associations of HCC risk with these SNPs were further confirmed in an independent case-control set. We found that 1p34.2 rs621559 and 14q21 rs398652 were significantly associated with HBV-related HCC risk (both P<0.005 after Bonferroni corrections). There was no significant difference of either rs621559 or rs398652 genotypes between chronic HBV carriers and healthy controls, demonstrating that the association was not due to predisposition to HBV infection. In the pooled analyses (1806 HBV-related HCC cases and 1954 controls), we observed a decreased HCC risk, 0.72-times, associated with the 1p34.2 rs621559 AA genotype compared to the GG genotype (P = 1.6×10⁻⁶). Additionally, there was an increased HCC risk, 1.27-fold, associated with the rs398652 GG genotype (P = 3.3×10⁻⁶). A statistical joint effect between the rs621559 GG and rs398652 GG genotypes may exist in elevating risk of HBV-related HCC. We show, for the first time, that rs398652 and rs621559 might be marker genetic variants for risk of HBV-related HCC in the Chinese population.     

HLA-DQB1*03 Confers Susceptibility to Chronic Hepatitis C in Japanese: A Genome-Wide Association Study

Hepatitis C virus (HCV) establishes a chronic infection in 70-80% of infected individuals. Many researchers have examined the effect of human leukocyte antigen (HLA) on viral persistence because of its critical role in the immune response against exposure to HCV, but almost all studies have proven to be inconclusive. To identify genetic risk factors for chronic HCV infection, we analyzed 458,207 single nucleotide polymorphisms (SNPs) in 481 chronic HCV patients and 2,963 controls in a Japanese cohort. Next, we performed a replication study with an independent panel of 4,358 cases and 1,114 controls. We further confirmed the association in 1,379 cases and 25,817 controls. In the GWAS phase, we found 17 SNPs that showed suggestive association (P < 1 × 10^-5). After the first replication study, we found one intronic SNP in the HLA-DQ locus associated with chronic HCV infection, and when we combined the two studies, the association reached the level of genome-wide significance. In the second replication study, we again confirmed the association (P _combined = 3.59 × 10⁻¹⁶, odds ratio [OR] = 0.79). Subsequent analysis revealed another SNP, rs1130380, with a stronger association (OR=0.72). This nucleotide substitution causes an amino acid substitution (R55P) in the HLA-DQB1 protein specific to the DQB1*03 allele, which is common worldwide. In addition, we confirmed an association with the previously reported IFNL3-IFNL4 locus and propose that the effect of DQB1*03 on HCV persistence might be affected by the IFNL4 polymorphism. Our findings suggest that a common amino acid substitution in HLA-DQB1 affects susceptibility to chronic infection with HCV in the Japanese population and may not be independent of the IFNL4 genotype.     

Genome-wide association analysis identifies variants associated with nonalcoholic fatty liver disease that have distinct effects on metabolic traits

Nonalcoholic fatty liver disease (NAFLD) clusters in families, but the only known common genetic variants influencing risk are near PNPLA3. We sought to identify additional genetic variants influencing NAFLD using genome-wide association (GWA) analysis of computed tomography (CT) measured hepatic steatosis, a non-invasive measure of NAFLD, in large population based samples. Using variance components methods, we show that CT hepatic steatosis is heritable (∼26%–27%) in family-based Amish, Family Heart, and Framingham Heart Studies (n = 880 to 3,070). By carrying out a fixed-effects meta-analysis of genome-wide association (GWA) results between CT hepatic steatosis and ∼2.4 million imputed or genotyped SNPs in 7,176 individuals from the Old Order Amish, Age, Gene/Environment Susceptibility-Reykjavik study (AGES), Family Heart, and Framingham Heart Studies, we identify variants associated at genome-wide significant levels (p<5×10⁻⁸) in or near PNPLA3, NCAN, and PPP1R3B. We genotype these and 42 other top CT hepatic steatosis-associated SNPs in 592 subjects with biopsy-proven NAFLD from the NASH Clinical Research Network (NASH CRN). In comparisons with 1,405 healthy controls from the Myocardial Genetics Consortium (MIGen), we observe significant associations with histologic NAFLD at variants in or near NCAN, GCKR, LYPLAL1, and PNPLA3, but not PPP1R3B. Variants at these five loci exhibit distinct patterns of association with serum lipids, as well as glycemic and anthropometric traits. We identify common genetic variants influencing CT–assessed steatosis and risk of NAFLD. Hepatic steatosis associated variants are not uniformly associated with NASH/fibrosis or result in abnormalities in serum lipids or glycemic and anthropometric traits, suggesting genetic heterogeneity in the pathways influencing these traits.     

Association Study of N-Methyl-D-Aspartate Receptor Subunit 2B (GRIN2B) Polymorphisms and Schizophrenia Symptoms in the Han Chinese Population

Schizophrenia (SZ) is a common and complex psychiatric disorder that has a significant genetic component. The glutamatergic system is the major excitatory neurotransmitter system in the central nervous system, and is mediated by N-methyl-D-aspartate (NMDA) receptors. Disturbances in this system have been hypothesized to play a major role in SZ pathogenesis. Several studies have revealed that the NMDA receptor subunit 2B (GRIN2B) potentially associates with SZ and its psychiatric symptoms. In this study, we performed a case–control study to identify polymorphisms of the GRIN2B gene that may confer susceptibility to SZ in the Han Chinese population. Thirty-four single nucleotide polymorphisms (SNPs) were genotyped in 528 paranoid SZ patients and 528 control subjects. A significant association was observed in allele and genotype between SZ and controls at rs2098469 (χ² = 8.425 and 4.994; p = 0.025 and 0.014, respectively). Significant associations were found in the allele at rs12319804 (χ² = 4.436; p = 0.035), as well as in the genotype at rs12820037 and rs7298664 between SZ and controls (χ² = 11.162 and 38.204; p = 0.003 and 4.27×10^-8, respectively). After applying the Bonferroni correction, rs7298664 still had significant genotype associations with SZ (p = 1.71×10^-7). In addition, rs2098469 genotype and allele frequencies, and 12820037 allele frequencies were nominally associated with SZ. Three haplotypes, CGA (rs10845849—rs12319804—rs10845851), CC (rs12582848—rs7952915), and AAGAC (rs2041986—rs11055665—rs7314376—rs7297101—rs2098469), had significant differences between SZ and controls (χ² = 4.324, 4.582, and 4.492; p = 0.037, 0.032, and 0.034, respectively). In addition, three SNPs, rs2098469, rs12820037, and rs7298664, were significantly associated with cognition factors PANSS subscores in SZ (F = 16.799, 7.112, and 13.357; p = 0.000, 0.017, and 0.000, respectively). In conclusion, our study provides novel evidence for an association between GRIN2B polymorphisms and SZ susceptibility and symptoms in the Han Chinese population.     

Genome-wide Association and Population Genetic Analysis of C-Reactive Protein in African American and Hispanic American Women

C-reactive protein (CRP) is a systemic inflammation marker that predicts future cardiovascular risk. CRP levels are higher in African Americans and Hispanic Americans than in European Americans, but the genetic determinants of CRP in these admixed United States minority populations are largely unknown. We performed genome-wide association studies (GWASs) of 8,280 African American (AA) and 3,548 Hispanic American (HA) postmenopausal women from the Women''s Health Initiative SNP Health Association Resource. We discovered and validated a CRP-associated variant of triggering receptors expressed by myeloid cells 2 (TREM2) in chromosomal region 6p21 (p = 10⁻¹⁰). The TREM2 variant associated with higher CRP is common in Africa but rare in other ancestral populations. In AA women, the CRP region in 1q23 contained a strong admixture association signal (p = 10⁻¹⁷), which appears to be related to several independent CRP-associated alleles; the strongest of these is present only in African ancestral populations and is associated with higher CRP. Of the other genomic loci previously associated with CRP through GWASs of European populations, most loci (LEPR, IL1RN, IL6R, GCKR, NLRP3, HNF1A, HNF4A, and APOC1) showed consistent patterns of association with CRP in AA and HA women. In summary, we have identified a common TREM2 variant associated with CRP in United States minority populations. The genetic architecture underlying the CRP phenotype in AA women is complex and involves genetic variants shared across populations, as well as variants specific to populations of African descent.     

Altered Cortico-Striatal Connectivity in Offspring of Schizophrenia Patients Relative to Offspring of Bipolar Patients and Controls

Schizophrenia (SZ) and bipolar disorder (BD) share clinical features, genetic risk factors and neuroimaging abnormalities. There is evidence of disrupted connectivity in resting state networks in patients with SZ and BD and their unaffected relatives. Resting state networks are known to undergo reorganization during youth coinciding with the period of increased incidence for both disorders. We therefore focused on characterizing resting state network connectivity in youth at familial risk for SZ or BD to identify alterations arising during this period. We measured resting-state functional connectivity in a sample of 106 youth, aged 7–19 years, comprising offspring of patients with SZ (N = 27), offspring of patients with BD (N = 39) and offspring of community control parents (N = 40). We used Independent Component Analysis to assess functional connectivity within the default mode, executive control, salience and basal ganglia networks and define their relationship to grey matter volume, clinical and cognitive measures. There was no difference in connectivity within any of the networks examined between offspring of patients with BD and offspring of community controls. In contrast, offspring of patients with SZ showed reduced connectivity within the left basal ganglia network compared to control offspring, and they showed a positive correlation between connectivity in this network and grey matter volume in the left caudate. Our findings suggest that dysconnectivity in the basal ganglia network is a robust correlate of familial risk for SZ and can be detected during childhood and adolescence.     

Intrahepatic Cholestasis of Pregnancy (ICP) in U.S. Latinas and Chileans: Clinical features,Ancestry Analysis,and Admixture Mapping

In the Americas, women with Indigenous American ancestry are at increased risk of intrahepatic cholestasis of pregnancy (ICP), relative to women of other ethnicities. We hypothesized that ancestry-related genetic factors contribute to this increased risk. We collected clinical and laboratory data, and performed biochemical assays on samples from U.S. Latinas and Chilean women, with and without ICP. The study sample included 198 women with ICP (90 from California, U.S., and 108 from Chile) and 174 pregnant control women (69 from California, U.S., and 105 from Chile). SNP genotyping was performed using Affymetrix arrays. We compared overall genetic ancestry between cases and controls, and used a genome-wide admixture mapping approach to screen for ICP susceptibility loci. We identified commonalities and differences in features of ICP between the 2 countries and determined that cases had a greater proportion of Indigenous American ancestry than did controls (p = 0.034). We performed admixture mapping, taking country of origin into account, and identified one locus for which Native American ancestry was associated with increased risk of ICP at a genome-wide level of significance (P = 3.1 x 10^-5, P_corrected = 0.035). This locus has an odds ratio of 4.48 (95% CI: 2.21-9.06) for 2 versus zero Indigenous American chromosomes. This locus lies on chromosome 2, with a 10 Mb 95% confidence interval which does not contain any previously identified hereditary ‘cholestasis genes.’ Our results indicate that genetic factors contribute to the risk of developing ICP in the Americas, and support the utility of clinical and genetic studies of ethnically mixed populations for increasing our understanding of ICP.     

Viral adaptation to host immune responses occurs in chronic hepatitis B virus (HBV) infection, and adaptation is greatest in HBV e antigen-negative disease

Hepatitis B virus (HBV)-specific T-cell responses are important in the natural history of HBV infection. The number of known HBV-specific T-cell epitopes is limited, and it is not clear whether viral evolution occurs in chronic HBV infection. We aimed to identify novel HBV T-cell epitopes by examining the relationship between HBV sequence variation and the human leukocyte antigen (HLA) type in a large prospective clinic-based cohort of Asian patients with chronic HBV infection recruited in Australia and China (n = 119). High-resolution 4-digit HLA class I and II typing and full-length HBV sequencing were undertaken for treatment-naïve individuals (52% with genotype B, 48% with genotype C, 63% HBV e antigen [HBeAg] positive). Statistically significant associations between HLA types and HBV sequence variation were identified (n = 49) at 41 sites in the HBV genome. Using prediction programs, we determined scores for binding between peptides containing these polymorphisms and associated HLA types. Among the regions that could be tested, HLA binding was predicted for 14/18 (78%). We identified several HLA-associated polymorphisms involving likely known anchor residues that resulted in altered predicted binding scores. Some HLA-associated polymorphisms fell within known T-cell epitopes with matching HLA restriction. Enhanced viral adaptation (defined as the presence of the relevant HLA and the escaped amino acid) was independently associated with HBeAg-negative disease (P = 0.003). Thus, HBV appears to be under immune pressure in chronic HBV infection, particularly in HBeAg-negative disease.     

Network-Based Multiple Sclerosis Pathway Analysis with GWAS Data from 15,000 Cases and 30,000 Controls

Multiple sclerosis (MS) is an inflammatory CNS disease with a substantial genetic component, originally mapped to only the human leukocyte antigen (HLA) region. In the last 5 years, a total of seven genome-wide association studies and one meta-analysis successfully identified 57 non-HLA susceptibility loci. Here, we merged nominal statistical evidence of association and physical evidence of interaction to conduct a protein-interaction-network-based pathway analysis (PINBPA) on two large genetic MS studies comprising a total of 15,317 cases and 29,529 controls. The distribution of nominally significant loci at the gene level matched the patterns of extended linkage disequilibrium in regions of interest. We found that products of genome-wide significantly associated genes are more likely to interact physically and belong to the same or related pathways. We next searched for subnetworks (modules) of genes (and their encoded proteins) enriched with nominally associated loci within each study and identified those modules in common between the two studies. We demonstrate that these modules are more likely to contain genes with bona fide susceptibility variants and, in addition, identify several high-confidence candidates (including BCL10, CD48, REL, TRAF3, and TEC). PINBPA is a powerful approach to gaining further insights into the biology of associated genes and to prioritizing candidates for subsequent genetic studies of complex traits.     

Novel associations for hypothyroidism include known autoimmune risk loci

Hypothyroidism is the most common thyroid disorder, affecting about 5% of the general population. Here we present the current largest genome-wide association study of hypothyroidism, in 3,736 cases and 35,546 controls. Hypothyroidism was assessed via web-based questionnaires. We identify five genome-wide significant associations, three of which are well known to be involved in a large spectrum of autoimmune diseases: rs6679677 near PTPN22, rs3184504 in SH2B3, and rs2517532 in the HLA class I region (p-values 2.8·10(-13), 2.6·10(-12), and 1.3·10(-8), respectively). We also report associations with rs4915077 near VAV3 (p-value 7.5·10(-10)) and rs925489 near FOXE1 (p value 2.4·10(-19)). VAV3 is involved in immune function, and FOXE1 and PTPN22 have previously been associated with hypothyroidism. Although the HLA class I region and SH2B3 have previously been linked with a number of autoimmune diseases, this is the first report of their association with thyroid disease. The VAV3 association is also novel. We also show suggestive evidence of association for hypothyroidism with a SNP in the HLA class II region (independent of the other HLA association) as well as SNPs in CAPZB, PDE8B, and CTLA4. CAPZB and PDE8B have been linked to TSH levels and CTLA4 to a variety of autoimmune diseases. These results suggest heterogeneity in the genetic etiology of hypothyroidism, implicating genes involved in both autoimmune disorders and thyroid function. Using a genetic risk profile score based on the top association from each of the five genome-wide significant regions in our study, the relative risk between the highest and lowest deciles of genetic risk is 2.0.     

Online-Datenbanken und systematische Metaanalysen komplex-genetischer Erkrankungen

Susceptibility to genetically complex disorders is determined by an unknown number of genetic determinants, and decades of intensive research have yielded hundreds of such potential susceptibility loci for Alzheimer’s disease (AD), Parkinson’s disease (PD), schizophrenia (SZ), and multiple sclerosis (MS). The results of genome-wide association studies are now adding to an already vast and complicated body of data. To facilitate the evaluation and interpretation of these findings, we have recently created online databases for genetic association studies in AD, PD, SZ, and MS. In addition to providing detailed summaries for each eligible study, the databases present the results of allele-based meta-analyses for all polymorphisms with sufficient genotype data. In this review, we discuss the background and implications of the database approach developed by our group, using current findings from the AD (AlzGene) and PD (PDGene) databases as examples.