Novel breast cancer susceptibility loci under linkage peaks identified in African ancestry consortia.

BACKGROUND: Expansion of genome-wide association studies across population groups is needed to improve our understanding of shared and unique genetic contributions to breast cancer. We performed association and replication studies guided by a priori linkage findings from African ancestry (AA) relative pairs. METHODS: We performed fixed-effect inverse-variance weighted meta-analysis under three significant AA breast cancer linkage peaks (3q26-27, 12q22-23, and 16q21-22) in 9241 AA cases and 10 193 AA controls. We examined associations with overall breast cancer as well as estrogen receptor (ER)-positive and negative subtypes (193,132 SNPs). We replicated associations in the African-ancestry Breast Cancer Genetic Consortium (AABCG). RESULTS: In AA women, we identified two associations on chr12q for overall breast cancer (rs1420647, OR = 1.15, p = 2.50×10-6; rs12322371, OR = 1.14, p = 3.15×10-6), and one for ER-negative breast cancer (rs77006600, OR = 1.67, p = 3.51×10-6). On chr3, we identified two associations with ER-negative disease (rs184090918, OR = 3.70, p = 1.23×10-5; rs76959804, OR = 3.57, p = 1.77×10-5) and on chr16q we identified an association with ER-negative disease (rs34147411, OR = 1.62, p = 8.82×10-6). In the replication study, the chr3 associations were significant and effect sizes were larger (rs184090918, OR: 6.66, 95% CI: 1.43, 31.01; rs76959804, OR: 5.24, 95% CI: 1.70, 16.16). CONCLUSION: The two chr3 SNPs are upstream to open chromatin ENSR00000710716, a regulatory feature that is actively regulated in mammary tissues, providing evidence that variants in this chr3 region may have a regulatory role in our target organ. Our study provides support for breast cancer variant discovery using prioritization based on linkage evidence.

Efficacy and Safety of Bepirovirsen in Chronic Hepatitis B Infection.

BACKGROUND: Bepirovirsen is an antisense oligonucleotide that targets all hepatitis B virus (HBV) messenger RNAs and acts to decrease levels of viral proteins. METHODS: We conducted a phase 2b, randomized, investigator-unblinded trial involving participants with chronic HBV infection who were receiving or not receiving nucleoside or nucleotide analogue (NA) therapy. Participants were randomly assigned (in a 3:3:3:1 ratio) to receive weekly subcutaneous injections of bepirovirsen at a dose of 300 mg for 24 weeks (group 1), bepirovirsen at a dose of 300 mg for 12 weeks then 150 mg for 12 weeks (group 2), bepirovirsen at a dose of 300 mg for 12 weeks then placebo for 12 weeks (group 3), or placebo for 12 weeks then bepirovirsen at a dose of 300 mg for 12 weeks (group 4). Groups 1, 2, and 3 received loading doses of bepirovirsen. The composite primary outcome was a hepatitis B surface antigen (HBsAg) level below the limit of detection and an HBV DNA level below the limit of quantification maintained for 24 weeks after the planned end of bepirovirsen treatment, without newly initiated antiviral medication. RESULTS: The intention-to-treat population comprised 457 participants (227 receiving NA therapy and 230 not receiving NA therapy). Among those receiving NA therapy, a primary-outcome event occurred in 6 participants (9%; 95% credible interval, 0 to 31) in group 1, in 6 (9%; 95% credible interval, 0 to 43) in group 2, in 2 (3%; 95% credible interval, 0 to 16) in group 3, and 0 (0%; post hoc credible interval, 0 to 8) in group 4. Among participants not receiving NA therapy, a primary-outcome event occurred in 7 participants (10%; 95% credible interval, 0 to 38), 4 (6%; 95% credible interval, 0 to 25), 1 (1%; post hoc credible interval, 0 to 6), and 0 (0%; post hoc credible interval, 0 to 8), respectively. During weeks 1 through 12, adverse events, including injection-site reactions, pyrexia, fatigue, and increased alanine aminotransferase levels, were more common with bepirovirsen (groups 1, 2, and 3) than with placebo (group 4). CONCLUSIONS: In this phase 2b trial, bepirovirsen at a dose of 300 mg per week for 24 weeks resulted in sustained HBsAg and HBV DNA loss in 9 to 10% of participants with chronic HBV infection. Larger and longer trials are required to assess the efficacy and safety of bepirovirsen. (Funded by GSK; B-Clear ClinicalTrials.gov number, NCT04449029.).

An Accidental Genetic Epidemiologist.

I briefly describe my early life and how, through a series of serendipitous events, I became a genetic epidemiologist. I discuss how the Elston-Stewart algorithm was discovered and its contribution to segregation, linkage, and association analysis. New linkage findings and paternity testing resulted from having a genotyping lab. The different meanings of interaction-statistical and biological-are clarified. The computer package S.A.G.E. (Statistical Analysis for Genetic Epidemiology), based on extensive method development over two decades, was conceived in 1986, flourished for 20 years, and is now freely available for use and further development. Finally, I describe methods to estimate and test hypotheses about familial correlations, and point out that the liability model often used to estimate disease heritability estimates the heritability of that liability, rather than of the disease itself, and so can be highly dependent on the assumed distribution of that liability.

Comparison of Pharmacokinetics of the GalNAc-Conjugated Antisense Oligonucleotide GSK3389404 in Participants with Chronic Hepatitis B Infection across the Asia-Pacific Region.

GSK3389404, an N-acetyl galactosamine-conjugated antisense oligonucleotide (ASO), was in clinical development for chronic hepatitis B (CHB) treatment. Few studies have examined ASOs in Asian participants. In this analysis, the plasma pharmacokinetics (PK) of GSK3389404 were characterized and compared in patients with CHB across the Asia-Pacific region (N = 64), including mainland China (n = 16), Hong Kong (n = 8), Japan (n = 21), South Korea (n = 12), Singapore (n = 4), and the Philippines (n = 3), from a phase 2a, multicenter, randomized, double-blind, placebo-controlled study (NCT03020745). Hepatitis B(e) antigen-positive and -negative patients (on or not on stable nucleos[t]ide regimens) received single (30 mg or 120 mg) or multiple (30 mg, 60 mg, or 120 mg weekly or 120 mg biweekly) subcutaneous GSK3389404 injections. The plasma concentrations were measured on day 1 in all cohorts as well as on days 29 and 57 in the multiple-dose cohorts. The GSK3389404 plasma PK were similar to those reported in a previous study in non-Asian healthy participants with a median time to peak concentration (tmax) of 1 to 4 h postdose, a mean half-life of 3 to 5 h across cohorts, and no accumulation following repeat dosing. The GSK3389404 plasma tmax and half-life values were dose-independent. The increase in the plasma peak concentration (Cmax) and the area under the concentration versus time curve (AUC) was dose-proportional from 60 to 120 mg and greater than dose-proportional from 30 to 60 or 120 mg. The GSK3389404 plasma concentration versus time profiles, half-life, tmax, Cmax, and AUC values were all comparable across the Asia-Pacific populations. Given the similarity of the PK among ASOs, this analysis suggests that the PK from any Asia-Pacific population may be used to guide ASO dose selection in the Asia-Pacific region.

Phase IIa, randomised, double-blind study of GSK3389404 in patients with chronic hepatitis B on stable nucleos(t)ide therapy.

BACKGROUND & AIMS: Bepirovirsen, an antisense oligonucleotide targeting pregenomic and mRNA transcripts of HBV, has been conjugated to N-acetyl galactosamine (GSK3389404) to enhance hepatocyte delivery. This dose-finding study was the first to assess GSK3389404 for chronic HBV infection. METHODS: This phase IIa, randomised, double-blind, placebo-controlled, 2-part study was conducted in 22 centres in Asia (NCT03020745). Pharmacokinetic findings from Part 1 informed Part 2 dosing. In Part 2, patients with chronic hepatitis B on nucleos(t)ide analogue therapy were randomised 11:2 to GSK3389404 (30, 60, 120 mg weekly or 120 mg bi-weekly) or placebo until Day 85. Coprimary endpoints included HBsAg response (≥1.5 log10 IU/ml reduction from baseline) rate, safety and pharmacokinetics. RESULTS: Parts 1 and 2 included 12 (9 GSK3389404, 3 placebo) and 66 patients (56 GSK3389404, 10 placebo), respectively. In Part 2, one patient each in the 60 mg weekly, 120 mg weekly and 120 mg bi-weekly arms achieved a HBsAg response. HBsAg reductions were dose-dependent (Day 85: mean 0.34 [60 mg weekly] to 0.75 log10 IU/ml [120 mg weekly]) and occurred in hepatitis B e antigen-positive and -negative patients. No patient achieved HBsAg seroclearance. 43/56 (77%) GSK3389404- and 9/10 (90%) placebo-treated patients reported adverse events. No deaths were reported. Alanine aminotransferase flares (>2x upper limit of normal) occurred in 2 GSK3389404-treated patients (120 mg weekly, 120 mg bi-weekly); both were associated with decreased HBsAg, but neither was considered a responder. GSK3389404 plasma concentrations peaked 2-4 hours post dose; mean plasma half-life was 3-5 hours. CONCLUSIONS: GSK3389404 showed an acceptable safety profile and target engagement, with dose-dependent reductions in HBsAg. However, no efficacious dosing regimen was identified. CLINICAL TRIAL NUMBER: NCT03020745. LAY SUMMARY: Hepatitis B virus (HBV) can result in chronic HBV infection, which may ultimately lead to chronic liver disease, primary liver cancer and death; HBV proteins may prevent the immune system from successfully controlling the virus. GSK3389404 is an investigational agent that targets HBV RNA, resulting in reduced viral protein production. This study assessed the safety of GSK3389404 and its ability to reduce the viral proteins in patients with chronic HBV infection. GSK3389404 showed dose-dependent reduction in hepatitis B surface antigen, with an acceptable safety profile. While no clear optimal dose was identified, the findings from this study may help in the development of improved treatment options for patients with chronic HBV infections.

Fisher's influence on me.

This is the 100th year anniversary of Fisher's 1918 paper "The correlation between relatives on the supposition of Mendelian inheritance" (Transactions of the Royal Society of Edinburgh 1918, 52 pp 899-438). Fisher's work has had a strong influence on today's genetic epidemiology and this brief autobiographical note highlights a few of the ways his influence on me has affected the field. Although I once took a course of lectures from Fisher, it was mainly his writings that influenced my statistical thinking. Not only did the concept of maximum likelihood appeal to me, but also the concepts of interclass and intraclass correlations, discriminant analysis, and transforming semiquantitative scores to minimize interactions-all topics I first learned about from the 11th edition of his book on Statistical Methods for Research Workers. This, together with a few serendipitous events that shaped my career, had a large influence on me and hence also on the field of genetic epidemiology.

Adjustment for covariates using summary statistics of genome-wide association studies.

Linear regression is a standard approach to identify genetic variants associated with continuous traits in genome-wide association studies (GWAS). In a standard epidemiology study, linear regression is often performed with adjustment for covariates to estimate the independent effect of a predictor variable or to improve statistical power by reducing residual variability. However, it is problematic to adjust for heritable covariates in genetic association analysis. Here, we propose a new method that utilizes summary statistics of the covariate from additional samples for reducing the residual variability and hence improves statistical power. Our simulation study showed that the proposed methodology can maintain a good control of Type I error and can achieve much higher power than a simple linear regression. The method is illustrated by an application to the GWAS results from the Genetic Investigation of Anthropometric Traits consortium.

ONETOOL for the analysis of family-based big data.

Motivation: Despite the need for separate tools to analyze family-based data, there are only a handful of tools optimized for family-based big data compared to the number of tools available for analyzing population-based data. Results: ONETOOL implements the properties of well-known existing family data analysis tools and recently developed methods in a computationally efficient manner, and so is suitable for analyzing the vast amount of variant data available from sequencing family members, providing a rich choice of analysis methods for big data on families. Availability and implementation: ONETOOL is freely available from http://healthstat.snu.ac.kr/software/onetool/. Supplementary information: Supplementary data are available at Bioinformatics online.

Applying family analyses to electronic health records to facilitate genetic research.

Motivation: Pedigree analysis is a longstanding and powerful approach to gain insight into the underlying genetic factors in human health, but identifying, recruiting and genotyping families can be difficult, time consuming and costly. Development of high throughput methods to identify families and foster downstream analyses are necessary. Results: This paper describes simple methods that allowed us to identify 173 368 family pedigrees with high probability using basic demographic data available in most electronic health records (EHRs). We further developed and validate a novel statistical method that uses EHR data to identify families more likely to have a major genetic component to their diseases risk. Lastly, we showed that incorporating EHR-linked family data into genetic association testing may provide added power for genetic mapping without additional recruitment or genotyping. The totality of these results suggests that EHR-linked families can enable classical genetic analyses in a high-throughput manner. Availability and implementation: Pseudocode is provided as supplementary information. Contact: HEBBRING.SCOTT@marshfieldresearch.org. Supplementary information: Supplementary data are available at Bioinformatics online.

How Consistent are Genetic Factors in Explaining Leisure-Time Physical Activity and Sport Participation? The Portuguese Healthy Families Study.

This study investigates how consistent genetic factors are, as measured by heritability estimates (h2), in the leisure-time physical activity index (LTPAI) and sport participation index (SPI) from early (10-14 yrs) to late adolescence (15-19 yrs). The sample comprises 12,385 subjects from 3,378 Portuguese nuclear families. Height and weight were measured and body mass index (BMI) was calculated, and the LTPAI and SPI were estimated by questionnaire. Socioeconomic status (SES) was assessed by parental occupation. Analyses were done using S.A.G.E. software. Our results showed that h2 estimates for the LTPAI and SPI in the two age groups (10-14 yrs and 15-19 yrs) were stable: for the LTPAI, h2 = 0.297 and 0.322, respectively; and for the SPI, h2 = 0.413 and 0.428, respectively. Sibling correlations and environmental correlations are higher in the younger age group for both the LTPAI and the SPI. Spousal correlations are higher in the younger age group for the LTPAI and lower for the SPI than the older group. Parent-offspring correlations are similar in both age groups for the LTPAI and SPI. In conclusion, the influence of genetic factors on physical activity and sport participation remains stable across age in adolescence. However, variation in sibling correlations - in particular, environmental correlations - was observed. These findings suggest that shared/non-shared environmental factors express different degrees of importance across age. Future intervention programs aiming to promote change in behaviors need to consider these results to bring about positive changes in physical activity and sport participation behaviors within the family setting.

Novel recurrently mutated genes in African American colon cancers.

We used whole-exome and targeted sequencing to characterize somatic mutations in 103 colorectal cancers (CRC) from African Americans, identifying 20 new genes as significantly mutated in CRC. Resequencing 129 Caucasian derived CRCs confirmed a 15-gene set as a preferential target for mutations in African American CRCs. Two predominant genes, ephrin type A receptor 6 (EPHA6) and folliculin (FLCN), with mutations exclusive to African American CRCs, are by genetic and biological criteria highly likely African American CRC driver genes. These previously unsuspected differences in the mutational landscapes of CRCs arising among individuals of different ethnicities have potential to impact on broader disparities in cancer behaviors.

Molecular Basis of Gene-Gene Interaction: Cyclic Cross-Regulation of Gene Expression and Post-GWAS Gene-Gene Interaction Involved in Atrial Fibrillation.

Atrial fibrillation (AF) is the most common cardiac arrhythmia at the clinic. Recent GWAS identified several variants associated with AF, but they account for <10% of heritability. Gene-gene interaction is assumed to account for a significant portion of missing heritability. Among GWAS loci for AF, only three were replicated in the Chinese Han population, including SNP rs2106261 (G/A substitution) in ZFHX3, rs2200733 (C/T substitution) near PITX2c, and rs3807989 (A/G substitution) in CAV1. Thus, we analyzed the interaction among these three AF loci. We demonstrated significant interaction between rs2106261 and rs2200733 in three independent populations and combined population with 2,020 cases/5,315 controls. Compared to non-risk genotype GGCC, two-locus risk genotype AATT showed the highest odds ratio in three independent populations and the combined population (OR=5.36 (95% CI 3.87-7.43), P=8.00×10-24). The OR of 5.36 for AATT was significantly higher than the combined OR of 3.31 for both GGTT and AACC, suggesting a synergistic interaction between rs2106261 and rs2200733. Relative excess risk due to interaction (RERI) analysis also revealed significant interaction between rs2106261 and rs2200733 when exposed two copies of risk alleles (RERI=2.87, P<1.00×10-4) or exposed to one additional copy of risk allele (RERI=1.29, P<1.00×10-4). The INTERSNP program identified significant genotypic interaction between rs2106261 and rs2200733 under an additive by additive model (OR=0.85, 95% CI: 0.74-0.97, P=0.02). Mechanistically, PITX2c negatively regulates expression of miR-1, which negatively regulates expression of ZFHX3, resulting in a positive regulation of ZFHX3 by PITX2c; ZFHX3 positively regulates expression of PITX2C, resulting in a cyclic loop of cross-regulation between ZFHX3 and PITX2c. Both ZFHX3 and PITX2c regulate expression of NPPA, TBX5 and NKX2.5. These results suggest that cyclic cross-regulation of gene expression is a molecular basis for gene-gene interactions involved in genetics of complex disease traits.

Genome-Wide Association and Trans-ethnic Meta-Analysis for Advanced Diabetic Kidney Disease: Family Investigation of Nephropathy and Diabetes (FIND).

Diabetic kidney disease (DKD) is the most common etiology of chronic kidney disease (CKD) in the industrialized world and accounts for much of the excess mortality in patients with diabetes mellitus. Approximately 45% of U.S. patients with incident end-stage kidney disease (ESKD) have DKD. Independent of glycemic control, DKD aggregates in families and has higher incidence rates in African, Mexican, and American Indian ancestral groups relative to European populations. The Family Investigation of Nephropathy and Diabetes (FIND) performed a genome-wide association study (GWAS) contrasting 6,197 unrelated individuals with advanced DKD with healthy and diabetic individuals lacking nephropathy of European American, African American, Mexican American, or American Indian ancestry. A large-scale replication and trans-ethnic meta-analysis included 7,539 additional European American, African American and American Indian DKD cases and non-nephropathy controls. Within ethnic group meta-analysis of discovery GWAS and replication set results identified genome-wide significant evidence for association between DKD and rs12523822 on chromosome 6q25.2 in American Indians (P = 5.74x10-9). The strongest signal of association in the trans-ethnic meta-analysis was with a SNP in strong linkage disequilibrium with rs12523822 (rs955333; P = 1.31x10-8), with directionally consistent results across ethnic groups. These 6q25.2 SNPs are located between the SCAF8 and CNKSR3 genes, a region with DKD relevant changes in gene expression and an eQTL with IPCEF1, a gene co-translated with CNKSR3. Several other SNPs demonstrated suggestive evidence of association with DKD, within and across populations. These data identify a novel DKD susceptibility locus with consistent directions of effect across diverse ancestral groups and provide insight into the genetic architecture of DKD.

On the association analysis of CNV data: a fast and robust family-based association method.

BACKGROUND: Copy number variation (CNV) is known to play an important role in the genetics of complex diseases and several methods have been proposed to detect association of CNV with phenotypes of interest. Statistical methods for CNV association analysis can be categorized into two different strategies. First, the copy number is estimated by maximum likelihood and association of the expected copy number with the phenotype is tested. Second, the observed probe intensity measurements can be directly used to detect association of CNV with the phenotypes of interest. RESULTS: For each strategy we provide a statistic that can be applied to extended families. The computational efficiency of the proposed methods enables genome-wide association analysis and we show with simulation studies that the proposed methods outperform other existing approaches. In particular, we found that the first strategy is always more efficient than the second strategy no matter whether copy numbers for each individual are well identified or not. With the proposed methods, we performed genome-wide CNV association analyses of hematological trait, hematocrit, on 521 Korean family samples. CONCLUSIONS: We found that statistical analysis with the expected copy number is more powerful than the statistic with the probe intensity measurements regardless of the accuracy of the estimation of copy numbers.

Family-Based Rare Variant Association Analysis: A Fast and Efficient Method of Multivariate Phenotype Association Analysis.

Family-based designs have been repeatedly shown to be powerful in detecting the significant rare variants associated with human diseases. Furthermore, human diseases are often defined by the outcomes of multiple phenotypes, and thus we expect multivariate family-based analyses may be very efficient in detecting associations with rare variants. However, few statistical methods implementing this strategy have been developed for family-based designs. In this report, we describe one such implementation: the multivariate family-based rare variant association tool (mFARVAT). mFARVAT is a quasi-likelihood-based score test for rare variant association analysis with multiple phenotypes, and tests both homogeneous and heterogeneous effects of each variant on multiple phenotypes. Simulation results show that the proposed method is generally robust and efficient for various disease models, and we identify some promising candidate genes associated with chronic obstructive pulmonary disease. The software of mFARVAT is freely available at http://healthstat.snu.ac.kr/software/mfarvat/, implemented in C++ and supported on Linux and MS Windows.

The association of the vanin-1 N131S variant with blood pressure is mediated by endoplasmic reticulum-associated degradation and loss of function.

High blood pressure (BP) is the most common cardiovascular risk factor worldwide and a major contributor to heart disease and stroke. We previously discovered a BP-associated missense SNP (single nucleotide polymorphism)-rs2272996-in the gene encoding vanin-1, a glycosylphosphatidylinositol (GPI)-anchored membrane pantetheinase. In the present study, we first replicated the association of rs2272996 and BP traits with a total sample size of nearly 30,000 individuals from the Continental Origins and Genetic Epidemiology Network (COGENT) of African Americans (P=0.01). This association was further validated using patient plasma samples; we observed that the N131S mutation is associated with significantly lower plasma vanin-1 protein levels. We observed that the N131S vanin-1 is subjected to rapid endoplasmic reticulum-associated degradation (ERAD) as the underlying mechanism for its reduction. Using HEK293 cells stably expressing vanin-1 variants, we showed that N131S vanin-1 was degraded significantly faster than wild type (WT) vanin-1. Consequently, there were only minimal quantities of variant vanin-1 present on the plasma membrane and greatly reduced pantetheinase activity. Application of MG-132, a proteasome inhibitor, resulted in accumulation of ubiquitinated variant protein. A further experiment demonstrated that atenolol and diltiazem, two current drugs for treating hypertension, reduce the vanin-1 protein level. Our study provides strong biological evidence for the association of the identified SNP with BP and suggests that vanin-1 misfolding and degradation are the underlying molecular mechanism.

Selecting SNPs informative for African, American Indian and European Ancestry: application to the Family Investigation of Nephropathy and Diabetes (FIND).

BACKGROUND: The presence of population structure in a sample may confound the search for important genetic loci associated with disease. Our four samples in the Family Investigation of Nephropathy and Diabetes (FIND), European Americans, Mexican Americans, African Americans, and American Indians are part of a genome- wide association study in which population structure might be particularly important. We therefore decided to study in detail one component of this, individual genetic ancestry (IGA). From SNPs present on the Affymetrix 6.0 Human SNP array, we identified 3 sets of ancestry informative markers (AIMs), each maximized for the information in one the three contrasts among ancestral populations: Europeans (HAPMAP, CEU), Africans (HAPMAP, YRI and LWK), and Native Americans (full heritage Pima Indians). We estimate IGA and present an algorithm for their standard errors, compare IGA to principal components, emphasize the importance of balancing information in the ancestry informative markers (AIMs), and test the association of IGA with diabetic nephropathy in the combined sample. RESULTS: A fixed parental allele maximum likelihood algorithm was applied to the FIND to estimate IGA in four samples: 869 American Indians; 1385 African Americans; 1451 Mexican Americans; and 826 European Americans. When the information in the AIMs is unbalanced, the estimates are incorrect with large error. Individual genetic admixture is highly correlated with principle components for capturing population structure. It takes ~700 SNPs to reduce the average standard error of individual admixture below 0.01. When the samples are combined, the resulting population structure creates associations between IGA and diabetic nephropathy. CONCLUSIONS: The identified set of AIMs, which include American Indian parental allele frequencies, may be particularly useful for estimating genetic admixture in populations from the Americas. Failure to balance information in maximum likelihood, poly-ancestry models creates biased estimates of individual admixture with large error. This also occurs when estimating IGA using the Bayesian clustering method as implemented in the program STRUCTURE. Odds ratios for the associations of IGA with disease are consistent with what is known about the incidence and prevalence of diabetic nephropathy in these populations.

Immunoglobulin G genotypes and the risk of schizophrenia.

Genes of the immune system are relevant to the etiology of schizophrenia. However, to our knowledge, no large-scale studies, using molecular methods, have been undertaken to investigate the role of highly polymorphic immunoglobulin GM (γ marker) genes in this disorder. In this investigation, we aimed to determine whether particular GM genotypes were associated with susceptibility to schizophrenia. Using a matched case-control study design, we analyzed DNA samples from 798 subjects-398 patients with schizophrenia and 400 controls-obtained from the U.S. National Institute of Mental Health Repository. GM alleles were determined by the TaqMan(®) genotyping assay. The GM 3/3; 23-/23- genotype was highly significantly associated with susceptibility to schizophrenia (p = 0.0002). Subjects with this genotype were over three times (OR 3.4; 95 % CI 1.7-6.7) as likely to develop schizophrenia as those without this genotype. Our results show that immunoglobulin GM genes are risk factors for the development of schizophrenia. Since GM alleles have been implicated in gluten sensitivity and in immunity to neurotropic viruses associated with cognitive impairment, the results presented here may help unify these two disparate areas of pathology affected in this disorder.