Genetic variation in the human leukocyte antigen region confers susceptibility to Clostridioides difficile infection

Clostridioides difficile (C. diff.) infection (CDI) is a leading cause of hospital acquired diarrhea in North America and Europe and a major cause of morbidity and mortality. Known risk factors do not fully explain CDI susceptibility, and genetic susceptibility is suggested by the fact that some patients with colons that are colonized with C. diff. do not develop any infection while others develop severe or recurrent infections. To identify common genetic variants associated with CDI, we performed a genome-wide association analysis in 19,861 participants (1349 cases; 18,512 controls) from the Electronic Medical Records and Genomics (eMERGE) Network. Using logistic regression, we found strong evidence for genetic variation in the DRB locus of the MHC (HLA) II region that predisposes individuals to CDI (P > 1.0 × 10-14; OR 1.56). Altered transcriptional regulation in the HLA region may play a role in conferring susceptibility to this opportunistic enteric pathogen.

Genome-wide trans-ancestry meta-analysis provides insight into the genetic architecture of type 2 diabetes susceptibility

To further understanding of the genetic basis of type 2 diabetes (T2D) susceptibility, we aggregated published meta-analyses of genome-wide association studies (GWAS) including 26,488 cases and 83,964 controls of European, East Asian, South Asian, and Mexican and Mexican American ancestry. We observed significant excess in directional consistency of T2D risk alleles across ancestry groups, even at SNPs demonstrating only weak evidence of association. By following up the strongest signals of association from the trans-ethnic meta-analysis in an additional 21,491 cases and 55,647 controls of European ancestry, we identified seven novel T2D susceptibility loci. Furthermore, we observed considerable improvements in fine-mapping resolution of common variant association signals at several T2D susceptibility loci. These observations highlight the benefits of trans-ethnic GWAS for the discovery and characterisation of complex trait loci, and emphasize an exciting opportunity to extend insight into the genetic architecture and pathogenesis of human diseases across populations of diverse ancestry.

Genome-wide association study of obsessive-compulsive disorder

Obsessive-compulsive disorder (OCD) is a common, debilitating neuropsychiatric illness with complex genetic etiology. The International OCD Foundation Genetics Collaborative (IOCDF-GC) is a multi-national collaboration established to discover the genetic variation predisposing to OCD. A set of individuals affected with DSM-IV OCD, a subset of their parents, and unselected controls, were genotyped with several different Illumina SNP microarrays. After extensive data cleaning, 1465 cases, 5557 ancestry-matched controls and 400 complete trios remained, with a common set of 469,410 autosomal and 9657 X-chromosome single nucleotide polymorphisms (SNPs). Ancestry-stratified case-control association analyses were conducted for three genetically-defined subpopulations and combined in two meta-analyses, with and without the trio-based analysis. In the case-control analysis, the lowest two P-values were located within DLGAP1 (P=2.49 × 10(-6) and P=3.44 × 10(-6)), a member of the neuronal postsynaptic density complex. In the trio analysis, rs6131295, near BTBD3, exceeded the genome-wide significance threshold with a P-value=3.84 × 10(-8). However, when trios were meta-analyzed with the case-control samples, the P-value for this variant was 3.62 × 10(-5), losing genome-wide significance. Although no SNPs were identified to be associated with OCD at a genome-wide significant level in the combined trio-case-control sample, a significant enrichment of methylation QTLs (P<0.001) and frontal lobe expression quantitative trait loci (eQTLs) (P=0.001) was observed within the top-ranked SNPs (P<0.01) from the trio-case-control analysis, suggesting these top signals may have a broad role in gene expression in the brain, and possibly in the etiology of OCD.

Large-scale association analysis provides insights into the genetic architecture and pathophysiology of type 2 diabetes

To extend understanding of the genetic architecture and molecular basis of type 2 diabetes (T2D), we conducted a meta-analysis of genetic variants on the Metabochip, including 34,840 cases and 114,981 controls, overwhelmingly of European descent. We identified ten previously unreported T2D susceptibility loci, including two showing sex-differentiated association. Genome-wide analyses of these data are consistent with a long tail of additional common variant loci explaining much of the variation in susceptibility to T2D. Exploration of the enlarged set of susceptibility loci implicates several processes, including CREBBP-related transcription, adipocytokine signaling and cell cycle regulation, in diabetes pathogenesis.

Evaluating genetic risk for prostate cancer among Japanese and Latinos

BACKGROUND

There have been few genome-wide association studies (GWAS) of prostate cancer among diverse populations. To search for novel prostate cancer risk variants, we conducted GWAS of prostate cancer in Japanese and Latinos. In addition, we tested prostate cancer risk variants and developed genetic risk models of prostate cancer for Japanese and Latinos.

METHODS

Our first-stage GWAS of prostate cancer included Japanese (cases/controls = 1,033/1,042) and Latino (cases/controls = 1,043/1,057) from the Multiethnic Cohort (MEC). Significant associations from stage I (P < 1.0 × 10(-4)) were examined in silico in GWAS of prostate cancer (stage II) in Japanese (cases/controls = 1,583/3,386) and Europeans (cases/controls = 1,854/1,894).

RESULTS

No novel stage I single-nucleotide polymorphism (SNP) outside of known risk regions reached genome-wide significance. For Japanese, in stage I, the most notable putative novel association was seen with 10 SNPs (P ≤ 8.0 × 10(-6)) at chromosome 2q33; however, this was not replicated in stage II. For Latinos, the most significant association was observed with rs17023900 at the known 3p12 risk locus (stage I: OR = 1.45; P = 7.01 × 10(-5) and stage II: OR = 1.58; P = 3.05 × 10(-7)). The majority of the established risk variants for prostate cancer, 79% and 88%, were positively associated with prostate cancer in Japanese and Latinos (stage I), respectively. The cumulative effects of these variants significantly influence prostate cancer risk (OR per allele = 1.10; P = 2.71 × 10(-25) and OR = 1.07; P = 1.02 × 10(-16) for Japanese and Latinos, respectively).

CONCLUSION AND IMPACT

Our GWAS of prostate cancer did not identify novel genome-wide significant variants. However, our findings show that established risk variants for prostate cancer significantly contribute to risk among Japanese and Latinos.

Pitfalls of merging GWAS data: lessons learned in the eMERGE network and quality control procedures to maintain high data quality

Genome-wide association studies (GWAS) are a useful approach in the study of the genetic components of complex phenotypes. Aside from large cohorts, GWAS have generally been limited to the study of one or a few diseases or traits. The emergence of biobanks linked to electronic medical records (EMRs) allows the efficient reuse of genetic data to yield meaningful genotype-phenotype associations for multiple phenotypes or traits. Phase I of the electronic MEdical Records and GEnomics (eMERGE-I) Network is a National Human Genome Research Institute-supported consortium composed of five sites to perform various genetic association studies using DNA repositories and EMR systems. Each eMERGE site has developed EMR-based algorithms to comprise a core set of 14 phenotypes for extraction of study samples from each site's DNA repository. Each eMERGE site selected samples for a specific phenotype, and these samples were genotyped at either the Broad Institute or at the Center for Inherited Disease Research using the Illumina Infinium BeadChip technology. In all, approximately 17,000 samples from across the five sites were genotyped. A unified quality control (QC) pipeline was developed by the eMERGE Genomics Working Group and used to ensure thorough cleaning of the data. This process includes examination of sample and marker quality and various batch effects. Upon completion of the genotyping and QC analyses for each site's primary study, eMERGE Coordinating Center merged the datasets from all five sites. This larger merged dataset reentered the established eMERGE QC pipeline. Based on lessons learned during the process, additional analyses and QC checkpoints were added to the pipeline to ensure proper merging. Here, we explore the challenges associated with combining datasets from different genotyping centers and describe the expansion to eMERGE QC pipeline for merged datasets. These additional steps will be useful as the eMERGE project expands to include additional sites in eMERGE-II, and also serve as a starting point for investigators merging multiple genotype datasets accessible through the National Center for Biotechnology Information in the database of Genotypes and Phenotypes. Our experience demonstrates that merging multiple datasets after additional QC can be an efficient use of genotype data despite new challenges that appear in the process.

Quality control procedures for genome-wide association studies

Genome-wide association studies (GWAS) are being conducted at an unprecedented rate in population-based cohorts and have increased our understanding of the pathophysiology of complex disease. Regardless of context, the practical utility of this information will ultimately depend upon the quality of the original data. Quality control (QC) procedures for GWAS are computationally intensive, operationally challenging, and constantly evolving. Here we enumerate some of the challenges in QC of GWAS data and describe the approaches that the electronic MEdical Records and Genomics (eMERGE) network is using for quality assurance in GWAS data, thereby minimizing potential bias and error in GWAS results. We discuss common issues associated with QC of GWAS data, including data file formats, software packages for data manipulation and analysis, sex chromosome anomalies, sample identity, sample relatedness, population substructure, batch effects, and marker quality. We propose best practices and discuss areas of ongoing and future research.

Variants in interferon-alpha pathway genes and response to pegylated interferon-Alpha2a plus ribavirin for treatment of chronic hepatitis C virus infection in the hepatitis C antiviral long-term treatment against cirrhosis trial

Combination treatment with pegylated-interferon-alpha (PEG IFN-alpha) and ribavirin, the current recommended therapy for chronic hepatitis C virus (HCV) infection, results in a sustained virological response (SVR) in only about half of patients. Because genes involved in the interferon-alpha pathway may affect antiviral responses, we analyzed the relationship between variants in these genes and SVR among participants in the Hepatitis C Antiviral Long-Term treatment Against Cirrhosis (HALT-C) trial. Patients had advanced chronic hepatitis C that had previously failed to respond to interferon-based treatment. Participants were treated with peginterferon-alpha2a and ribavirin during the trial. Subjects with undetectable HCV RNA at week 72 were considered to have had an SVR. Subjects with detectable HCV RNA at week 20 were considered nonresponders. We used TaqMan assays to genotype 56 polymorphisms found in 13 genes in the interferon-alpha pathway. This analysis compares genotypes for participants with an SVR to nonresponders. The primary analysis was restricted to European American participants because a priori statistical power was low among the small number (n = 131) of African American patients. We used logistic regression to control the effect of other variables that are associated with treatment response. Among 581 European American patients, SVR was associated with IFNAR1 IVS1-22G (adjusted odds ratio, 0.57; P = 0.02); IFNAR2 Ex2-33C (adjusted odds ratio, 2.09; P = 0.02); JAK1 IVS22+112T (adjusted odds ratio, 1.66; P = 0.04); and ADAR Ex9+14A (adjusted odds ratio, 1.67; P = 0.03). For the TYK2-2256A promoter region variant, a borderline association was present among European American participants (OR, 1.51; P = 0.05) and a strong relationship among African American patients; all 10 with SVR who were genotyped for TYK2 -2256 carried the A variant compared with 68 of 120 (57%) nonresponders (P = 0.006).

CONCLUSION

Genetic polymorphisms in the interferon-alpha pathway may affect responses to antiviral therapy of chronic hepatitis C.