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.

Limited clinical utility of a genetic risk score for the prediction of fracture risk in elderly subjects.

It is important to identify the patients at highest risk of fractures. A recent large-scale meta-analysis identified 63 autosomal single-nucleotide polymorphisms (SNPs) associated with bone mineral density (BMD), of which 16 were also associated with fracture risk. Based on these findings, two genetic risk scores (GRS63 and GRS16) were developed. Our aim was to determine the clinical usefulness of these GRSs for the prediction of BMD, BMD change, and fracture risk in elderly subjects. We studied two male (Osteoporotic Fractures in Men Study [MrOS] US, MrOS Sweden) and one female (Study of Osteoporotic Fractures [SOF]) large prospective cohorts of older subjects, looking at BMD, BMD change, and radiographically and/or medically confirmed incident fractures (8067 subjects, 2185 incident nonvertebral or vertebral fractures). GRS63 was associated with BMD (≅3% of the variation explained) but not with BMD change. Both GRS63 and GRS16 were associated with fractures. After BMD adjustment, the effect sizes for these associations were substantially reduced. Similar results were found using an unweighted GRS63 and an unweighted GRS16 compared with those found using the corresponding weighted risk scores. Only minor improvements in C-statistics (AUC) for fractures were found when the GRSs were added to a base model (age, weight, and height), and no significant improvements in C-statistics were found when they were added to a model further adjusted for BMD. Net reclassification improvements with the addition of the GRSs to a base model were modest and substantially attenuated in BMD-adjusted models. GRS63 is associated with BMD, but not BMD change, suggesting that the genetic determinants of BMD differ from those of BMD change. When BMD is known, the clinical utility of the two GRSs for fracture prediction is limited in elderly subjects.

The chromosome 3q25 genomic region is associated with measures of adiposity in newborns in a multi-ethnic genome-wide association study.

Newborns characterized as large and small for gestational age are at risk for increased mortality and morbidity during the first year of life as well as for obesity and dysglycemia as children and adults. The intrauterine environment and fetal genes contribute to the fetal size at birth. To define the genetic architecture underlying the newborn size, we performed a genome-wide association study (GWAS) in 4281 newborns in four ethnic groups from the Hyperglycemia and Adverse Pregnancy Outcome Study. We tested for association with newborn anthropometric traits (birth length, head circumference, birth weight, percent fat mass and sum of skinfolds) and newborn metabolic traits (cord glucose and C-peptide) under three models. Model 1 adjusted for field center, ancestry, neonatal gender, gestational age at delivery, parity, maternal age at oral glucose tolerance test (OGTT); Model 2 adjusted for Model 1 covariates, maternal body mass index (BMI) at OGTT, maternal height at OGTT, maternal mean arterial pressure at OGTT, maternal smoking and drinking; Model 3 adjusted for Model 2 covariates, maternal glucose and C-peptide at OGTT. Strong evidence for association was observed with measures of newborn adiposity (sum of skinfolds model 3 Z-score 7.356, P = 1.90×10⁻¹³, and to a lesser degree fat mass and birth weight) and a region on Chr3q25.31 mapping between CCNL and LEKR1. These findings were replicated in an independent cohort of 2296 newborns. This region has previously been shown to be associated with birth weight in Europeans. The current study suggests that association of this locus with birth weight is secondary to an effect on fat as opposed to lean body mass.

Genomic, transcriptomic, and lipidomic profiling highlights the role of inflammation in individuals with low high-density lipoprotein cholesterol.

OBJECTIVE: Low high-density lipoprotein cholesterol (HDL-C) is associated with cardiometabolic pathologies. In this study, we investigate the biological pathways and individual genes behind low HDL-C by integrating results from 3 high-throughput data sources: adipose tissue transcriptomics, HDL lipidomics, and dense marker genotypes from Finnish individuals with low or high HDL-C (n=450). APPROACH AND RESULTS: In the pathway analysis of genetic data, we demonstrate that genetic variants within inflammatory pathways were enriched among low HDL-C associated single-nucleotide polymorphisms, and the expression of these pathways upregulated in the adipose tissue of low HDL-C subjects. The lipidomic analysis highlighted the change in HDL particle quality toward putatively more inflammatory and less vasoprotective state in subjects with low HDL-C, as evidenced by their decreased antioxidative plasmalogen contents. We show that the focal point of these inflammatory pathways seems to be the HLA region with its low HDL-associated alleles also associating with more abundant local transcript levels in adipose tissue, increased plasma vascular cell adhesion molecule 1 (VCAM1) levels, and decreased HDL particle plasmalogen contents, markers of adipose tissue inflammation, vascular inflammation, and HDL antioxidative potential, respectively. In a population-based look-up of the inflammatory pathway single-nucleotide polymorphisms in a large Finnish cohorts (n=11 211), no association of the HLA region was detected for HDL-C as quantitative trait, but with extreme HDL-C phenotypes, implying the presence of low or high HDL genes in addition to the population-genomewide association studies-identified HDL genes. CONCLUSIONS: Our study highlights the role of inflammation with a genetic component in subjects with low HDL-C and identifies novel cis-expression quantitative trait loci (cis-eQTL) variants in HLA region to be associated with low HDL-C.

Incidental genetic findings in randomized clinical trials: recommendations from the Genomics and Randomized Trials Network (GARNET).

Recommendations and guidance on how to handle the return of genetic results to patients have offered limited insight into how to approach incidental genetic findings in the context of clinical trials. This paper provides the Genomics and Randomized Trials Network (GARNET) recommendations on incidental genetic findings in the context of clinical trials, and discusses the ethical and practical issues considered in formulating our recommendations. There are arguments in support of as well as against returning incidental genetic findings in clinical trials. For instance, reporting incidental findings in clinical trials may improve the investigator-participant relationship and the satisfaction of participation, but it may also blur the line between clinical care and research. The issues of whether and how to return incidental genetic findings, including the costs of doing so, should be considered when developing clinical trial protocols. Once decided, plans related to sharing individual results from the aim(s) of the trial, as well as incidental findings, should be discussed explicitly in the consent form. Institutional Review Boards (IRBs) and other study-specific governing bodies should be part of the decision as to if, when, and how to return incidental findings, including when plans in this regard are being reconsidered.

Genetic variation associated with circulating monocyte count in the eMERGE Network.

With white blood cell count emerging as an important risk factor for chronic inflammatory diseases, genetic associations of differential leukocyte types, specifically monocyte count, are providing novel candidate genes and pathways to further investigate. Circulating monocytes play a critical role in vascular diseases such as in the formation of atherosclerotic plaque. We performed a joint and ancestry-stratified genome-wide association analyses to identify variants specifically associated with monocyte count in 11 014 subjects in the electronic Medical Records and Genomics Network. In the joint and European ancestry samples, we identified novel associations in the chromosome 16 interferon regulatory factor 8 (IRF8) gene (P-value = 2.78×10(-16), ß = -0.22). Other monocyte associations include novel missense variants in the chemokine-binding protein 2 (CCBP2) gene (P-value = 1.88×10(-7), ß = 0.30) and a region of replication found in ribophorin I (RPN1) (P-value = 2.63×10(-16), ß = -0.23) on chromosome 3. The CCBP2 and RPN1 region is located near GATA binding protein2 gene that has been previously shown to be associated with coronary heart disease. On chromosome 9, we found a novel association in the prostaglandin reductase 1 gene (P-value = 2.29×10(-7), ß = 0.16), which is downstream from lysophosphatidic acid receptor 1. This region has previously been shown to be associated with monocyte count. We also replicated monocyte associations of genome-wide significance (P-value = 5.68×10(-17), ß = -0.23) at the integrin, alpha 4 gene on chromosome 2. The novel IRF8 results and further replications provide supporting evidence of genetic regions associated with monocyte count.

Linkage analysis identifies a locus for plasma von Willebrand factor undetected by genome-wide association.

The plasma glycoprotein von Willebrand factor (VWF) exhibits fivefold antigen level variation across the normal human population determined by both genetic and environmental factors. Low levels of VWF are associated with bleeding and elevated levels with increased risk for thrombosis, myocardial infarction, and stroke. To identify additional genetic determinants of VWF antigen levels and to minimize the impact of age and illness-related environmental factors, we performed genome-wide association analysis in two young and healthy cohorts (n = 1,152 and n = 2,310) and identified signals at ABO (P < 7.9E-139) and VWF (P < 5.5E-16), consistent with previous reports. Additionally, linkage analysis based on sibling structure within the cohorts, identified significant signals at chromosome 2q12-2p13 (LOD score 5.3) and at the ABO locus on chromosome 9q34 (LOD score 2.9) that explained 19.2% and 24.5% of the variance in VWF levels, respectively. Given its strong effect, the linkage region on chromosome 2 could harbor a potentially important determinant of bleeding and thrombosis risk. The absence of a chromosome 2 association signal in this or previous association studies suggests a causative gene harboring many genetic variants that are individually rare, but in aggregate common. These results raise the possibility that similar loci could explain a significant portion of the "missing heritability" for other complex genetic traits.

New susceptibility loci associated with kidney disease in type 1 diabetes.

Diabetic kidney disease, or diabetic nephropathy (DN), is a major complication of diabetes and the leading cause of end-stage renal disease (ESRD) that requires dialysis treatment or kidney transplantation. In addition to the decrease in the quality of life, DN accounts for a large proportion of the excess mortality associated with type 1 diabetes (T1D). Whereas the degree of glycemia plays a pivotal role in DN, a subset of individuals with poorly controlled T1D do not develop DN. Furthermore, strong familial aggregation supports genetic susceptibility to DN. However, the genes and the molecular mechanisms behind the disease remain poorly understood, and current therapeutic strategies rarely result in reversal of DN. In the GEnetics of Nephropathy: an International Effort (GENIE) consortium, we have undertaken a meta-analysis of genome-wide association studies (GWAS) of T1D DN comprising ~2.4 million single nucleotide polymorphisms (SNPs) imputed in 6,691 individuals. After additional genotyping of 41 top ranked SNPs representing 24 independent signals in 5,873 individuals, combined meta-analysis revealed association of two SNPs with ESRD: rs7583877 in the AFF3 gene (P = 1.2 × 10(-8)) and an intergenic SNP on chromosome 15q26 between the genes RGMA and MCTP2, rs12437854 (P = 2.0 × 10(-9)). Functional data suggest that AFF3 influences renal tubule fibrosis via the transforming growth factor-beta (TGF-ß1) pathway. The strongest association with DN as a primary phenotype was seen for an intronic SNP in the ERBB4 gene (rs7588550, P = 2.1 × 10(-7)), a gene with type 2 diabetes DN differential expression and in the same intron as a variant with cis-eQTL expression of ERBB4. All these detected associations represent new signals in the pathogenesis of DN.

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.

Return of individual research results from genome-wide association studies: experience of the Electronic Medical Records and Genomics (eMERGE) Network.

PURPOSE: Return of individual genetic results to research participants, including participants in archives and biorepositories, is receiving increased attention. However, few groups have deliberated on specific results or weighed deliberations against relevant local contextual factors. METHODS: The Electronic Medical Records and Genomics (eMERGE) Network, which includes five biorepositories conducting genome-wide association studies, convened a return of results oversight committee to identify potentially returnable results. Network-wide deliberations were then brought to local constituencies for final decision making. RESULTS: Defining results that should be considered for return required input from clinicians with relevant expertise and much deliberation. The return of results oversight committee identified two sex chromosomal anomalies, Klinefelter syndrome and Turner syndrome, as well as homozygosity for factor V Leiden, as findings that could warrant reporting. Views about returning findings of HFE gene mutations associated with hemochromatosis were mixed due to low penetrance. Review of electronic medical records suggested that most participants with detected abnormalities were unaware of these findings. Local considerations relevant to return varied and, to date, four sites have elected not to return findings (return was not possible at one site). CONCLUSION: The eMERGE experience reveals the complexity of return of results decision making and provides a potential deliberative model for adoption in other collaborative contexts.

Genetic variants associated with the white blood cell count in 13,923 subjects in the eMERGE Network.

White blood cell count (WBC) is unique among identified inflammatory predictors of chronic disease in that it is routinely measured in asymptomatic patients in the course of routine patient care. We led a genome-wide association analysis to identify variants associated with WBC levels in 13,923 subjects in the electronic Medical Records and Genomics (eMERGE) Network. We identified two regions of interest that were each unique to subjects of genetically determined ancestry to the African continent (AA) or to the European continent (EA). WBC varies among different ancestry groups. Despite being ancestry specific, these regions were identifiable in the combined analysis. In AA subjects, the region surrounding the Duffy antigen/chemokine receptor gene (DARC) on 1q21 exhibited significant association (p value = 6.71e-55). These results validate the previously reported association between WBC and of the regulatory variant rs2814778 in the promoter region, which causes the Duffy negative phenotype (Fy-/-). A second missense variant (rs12075) is responsible for the two principal antigens, Fya and Fyb of the Duffy blood group system. The two variants, consisting of four alleles, act in concert to produce five antigens and subsequent phenotypes. We were able to identify the marginal and novel interaction effects of these two variants on WBC. In the EA subjects, we identified significantly associated SNPs tagging three separate genes in the 17q21 region: (1) GSDMA, (2) MED24, and (3) PSMD3. Variants in this region have been reported to be associated with WBC, neutrophil count, and inflammatory diseases including asthma and Crohn's disease.

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.

Common variants near CAV1 and CAV2 are associated with primary open-angle glaucoma in Caucasians from the USA.

Primary open-angle glaucoma (POAG) is a genetically complex common disease characterized by progressive optic nerve degeneration that results in irreversible blindness. Recently, a genome-wide association study (GWAS) for POAG in an Icelandic population identified significant associations with single nucleotide polymorphisms (SNPs) between the CAV1 and CAV2 genes on chromosome 7q31. In this study, we confirm that the identified SNPs are associated with POAG in our Caucasian US population and that specific haplotypes located in the CAV1/CAV2 intergenic region are associated with the disease. We also present data suggesting that associations with several CAV1/CAV2 SNPs are significant mostly in women.

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.

Common genetic variation in multiple metabolic pathways influences susceptibility to low HDL-cholesterol and coronary heart disease.

A low level of HDL-C is the most common plasma lipid abnormality observed in men with established coronary heart disease (CHD). To identify allelic variants associated with susceptibility to low HDL-C and CHD, we examined 60 candidate genes with key roles in HDL metabolism, insulin resistance, and inflammation using samples from the Veterans Affairs HDL Intervention Trial (VA-HIT; cases, n = 699) and the Framingham Offspring Study (FOS; controls, n = 705). VA-HIT was designed to examine the benefits of HDL-raising with gemfibrozil in men with low HDL-C (≤40 mg/dl) and established CHD. After adjustment for multiple testing within each gene, single-nucleotide polymorphisms (SNP) significantly associated with case status were identified in the genes encoding LIPC (rs4775065, P < 0.0001); CETP (rs5882, P = 0.0002); RXRA (rs11185660, P = 0.0021); ABCA1 (rs2249891, P = 0.0126); ABCC6 (rs150468, P = 0.0206; rs212077, P = 0.0443); CUBN (rs7893395, P = 0.0246); APOA2 (rs3813627, P = 0.0324); SELP (rs732314, P = 0.0376); and APOC4 (rs10413089, P = 0.0425). Included among the novel findings of this study are the identification of susceptibility alleles for low HDL-C/CHD risk in the genes encoding CUBN and RXRA, and the observation that genetic variation in SELP may influence CHD risk through its effects on HDL.

Quality control and quality assurance in genotypic data for genome-wide association studies.

Genome-wide scans of nucleotide variation in human subjects are providing an increasing number of replicated associations with complex disease traits. Most of the variants detected have small effects and, collectively, they account for a small fraction of the total genetic variance. Very large sample sizes are required to identify and validate findings. In this situation, even small sources of systematic or random error can cause spurious results or obscure real effects. The need for careful attention to data quality has been appreciated for some time in this field, and a number of strategies for quality control and quality assurance (QC/QA) have been developed. Here we extend these methods and describe a system of QC/QA for genotypic data in genome-wide association studies (GWAS). This system includes some new approaches that (1) combine analysis of allelic probe intensities and called genotypes to distinguish gender misidentification from sex chromosome aberrations, (2) detect autosomal chromosome aberrations that may affect genotype calling accuracy, (3) infer DNA sample quality from relatedness and allelic intensities, (4) use duplicate concordance to infer SNP quality, (5) detect genotyping artifacts from dependence of Hardy-Weinberg equilibrium test P-values on allelic frequency, and (6) demonstrate sensitivity of principal components analysis to SNP selection. The methods are illustrated with examples from the "Gene Environment Association Studies" (GENEVA) program. The results suggest several recommendations for QC/QA in the design and execution of GWAS.

Genetic variants at 2q24 are associated with susceptibility to type 2 diabetes.

To identify type 2 diabetes (T2D) susceptibility loci, we conducted genome-wide association (GWA) scans in nested case-control samples from two prospective cohort studies, including 2591 patients and 3052 controls of European ancestry. Validation was performed in 11 independent GWA studies of 10,870 cases and 73,735 controls. We identified significantly associated variants near RBMS1 and ITGB6 genes at 2q24, best-represented by SNP rs7593730 (combined OR=0.90, 95% CI=0.86-0.93; P=3.7x10(-8)). The frequency of the risk-lowering allele T is 0.23. Variants in this region were nominally related to lower fasting glucose and HOMA-IR in the MAGIC consortium (P<0.05). These data suggest that the 2q24 locus may influence the T2D risk by affecting glucose metabolism and insulin resistance.

Genetic variants in ABO blood group region, plasma soluble E-selectin levels and risk of type 2 diabetes.

Blood soluble E-selectin (sE-selectin) levels have been related to various conditions such as type 2 diabetes. We performed a genome-wide association study among women of European ancestry from the Nurses' Health Study, and identified genome-wide significant associations between a cluster of markers at the ABO locus (9q34) and plasma sE-selectin concentration. The strongest association was with rs651007, which explained approximately 9.71% of the variation in sE-selectin concentrations. SNP rs651007 was also nominally associated with soluble intracellular cell adhesion molecule-1 (sICAM-1) (P = 0.026) and TNF-R2 levels (P = 0.018), independent of sE-selectin. In addition, the genetic-inferred ABO blood group genotypes were associated with sE-selectin concentrations (P = 3.55 x 10(-47)). Moreover, we found that the genetic-inferred blood group B was associated with a decreased risk (OR = 0.44, 0.27-0.70) of type 2 diabetes compared with blood group O, adjusting for sE-selectin, sICAM-1, TNF-R2 and other covariates. Our findings indicate that the genetic variants at ABO locus affect plasma sE-selectin levels and diabetes risk. The genetic associations with diabetes risk were independent of sE-selectin levels.