Loci identified by a genome-wide association study of carotid artery stenosis in the eMERGE network.

Carotid artery atherosclerotic disease (CAAD) is a risk factor for stroke. We used a genome-wide association (GWAS) approach to discover genetic variants associated with CAAD in participants in the electronic Medical Records and Genomics (eMERGE) Network. We identified adult CAAD cases with unilateral or bilateral carotid artery stenosis and controls without evidence of stenosis from electronic health records at eight eMERGE sites. We performed GWAS with a model adjusting for age, sex, study site, and genetic principal components of ancestry. In eMERGE we found 1793 CAAD cases and 17,958 controls. Two loci reached genome-wide significance, on chr6 in LPA (rs10455872, odds ratio [OR] (95% confidence interval [CI]) = 1.50 (1.30-1.73), p = 2.1 × 10-8 ) and on chr7, an intergenic single nucleotide variant (SNV; rs6952610, OR (95% CI) = 1.25 (1.16-1.36), p = 4.3 × 10-8 ). The chr7 association remained significant in the presence of the LPA SNV as a covariate. The LPA SNV was also associated with coronary heart disease (CHD; 4199 cases and 11,679 controls) in this study (OR (95% CI) = 1.27 (1.13-1.43), p = 5 × 10-5 ) but the chr7 SNV was not (OR (95% CI) = 1.03 (0.97-1.09), p = .37). Both variants replicated in UK Biobank. Elevated lipoprotein(a) concentrations ([Lp(a)]) and LPA variants associated with elevated [Lp(a)] have previously been associated with CAAD and CHD, including rs10455872. With electronic health record phenotypes in eMERGE and UKB, we replicated a previously known association and identified a novel locus associated with CAAD.

Genome-Wide Association Study Identifies Two Novel Loci Associated with Female Stress and Urgency Urinary Incontinence.

PURPOSE: Genome-wide association studies have not identified replicable genetic risk loci for stress or urgency urinary incontinence. MATERIALS AND METHODS: We carried out a discovery stage, case control, genome-wide association study in 3 independent discovery cohorts of European women (8,979) for stress incontinence, urgency incontinence, and any incontinence phenotypes. We conducted replication in 6 additional studies of European ancestry (4,069). We collected bladder biopsies from women with incontinence (50) to further investigate bladder expression of implicated genes and pathways and used symptom questionnaires for phenotyping. We conducted meta-analyses using inverse variance fixed effects models and whole transcriptome analyses using Affymetrix® arrays with replication with TaqMan® polymerase chain reaction. RESULTS: In the discovery stage, we identified 16 single nucleotide polymorphisms genotyped or imputed at 5 loci that reached genome-wide significance (p <5×10-8). In replication, rs138724718 on chromosome 2 near the macrophage receptor with collagenous structure (MARCO) gene (replication p=0.003) was associated with stress incontinence. In addition, rs34998271 on chromosome 6 near the endothelin 1 (EDN1) gene (replication p=0.0008) was associated with urgency incontinence. In combined meta-analyses of discovery and replication cohorts, associations with genome-wide significance for these 2 single nucleotide polymorphisms were confirmed. Transcriptomics analyses showed differential expression of 7 of 19 genes in the endothelin pathway between stress and urgency incontinence (p <0.0001). CONCLUSIONS: We uncovered 2 new risk loci near the genes endothelin 1 (EDN1), associated with urgency incontinence, and macrophage receptor with collagenous structure (MARCO), associated with stress incontinence. These loci are biologically plausible given their roles in smooth muscle contraction and innate host defense, respectively.

Genetic variation in the SIM1 locus is associated with erectile dysfunction.

Erectile dysfunction affects millions of men worldwide. Twin studies support the role of genetic risk factors underlying erectile dysfunction, but no specific genetic variants have been identified. We conducted a large-scale genome-wide association study of erectile dysfunction in 36,649 men in the multiethnic Kaiser Permanente Northern California Genetic Epidemiology Research in Adult Health and Aging cohort. We also undertook replication analyses in 222,358 men from the UK Biobank. In the discovery cohort, we identified a single locus (rs17185536-T) on chromosome 6 near the single-minded family basic helix-loop-helix transcription factor 1 (SIM1) gene that was significantly associated with the risk of erectile dysfunction (odds ratio = 1.26, P = 3.4 × 10-25). The association replicated in the UK Biobank sample (odds ratio = 1.25, P = 6.8 × 10-14), and the effect is independent of known erectile dysfunction risk factors, including body mass index (BMI). The risk locus resides on the same topologically associating domain as SIM1 and interacts with the SIM1 promoter, and the rs17185536-T risk allele showed differential enhancer activity. SIM1 is part of the leptin-melanocortin system, which has an established role in body weight homeostasis and sexual function. Because the variants associated with erectile dysfunction are not associated with differences in BMI, our findings suggest a mechanism that is specific to sexual function.

Genome-wide meta-analysis of 158,000 individuals of European ancestry identifies three loci associated with chronic back pain.

Back pain is the #1 cause of years lived with disability worldwide, yet surprisingly little is known regarding the biology underlying this symptom. We conducted a genome-wide association study (GWAS) meta-analysis of chronic back pain (CBP). Adults of European ancestry were included from 15 cohorts in the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium, and from the UK Biobank interim data release. CBP cases were defined as those reporting back pain present for ≥3-6 months; non-cases were included as comparisons ("controls"). Each cohort conducted genotyping using commercially available arrays followed by imputation. GWAS used logistic regression models with additive genetic effects, adjusting for age, sex, study-specific covariates, and population substructure. The threshold for genome-wide significance in the fixed-effect inverse-variance weighted meta-analysis was p<5×10(-8). Suggestive (p<5×10(-7)) and genome-wide significant (p<5×10(-8)) variants were carried forward for replication or further investigation in the remaining UK Biobank participants not included in the discovery sample. The discovery sample comprised 158,025 individuals, including 29,531 CBP cases. A genome-wide significant association was found for the intronic variant rs12310519 in SOX5 (OR 1.08, p = 7.2×10(-10)). This was subsequently replicated in 283,752 UK Biobank participants not included in the discovery sample, including 50,915 cases (OR 1.06, p = 5.3×10(-11)), and exceeded genome-wide significance in joint meta-analysis (OR 1.07, p = 4.5×10(-19)). We found suggestive associations at three other loci in the discovery sample, two of which exceeded genome-wide significance in joint meta-analysis: an intergenic variant, rs7833174, located between CCDC26 and GSDMC (OR 1.05, p = 4.4×10(-13)), and an intronic variant, rs4384683, in DCC (OR 0.97, p = 2.4×10(-10)). In this first reported meta-analysis of GWAS for CBP, we identified and replicated a genetic locus associated with CBP (SOX5). We also identified 2 other loci that reached genome-wide significance in a 2-stage joint meta-analysis (CCDC26/GSDMC and DCC).

The eMERGE genotype set of 83,717 subjects imputed to ~40 million variants genome wide and association with the herpes zoster medical record phenotype.

The Electronic Medical Records and Genomics (eMERGE) network is a network of medical centers with electronic medical records linked to existing biorepository samples for genomic discovery and genomic medicine research. The network sought to unify the genetic results from 78 Illumina and Affymetrix genotype array batches from 12 contributing medical centers for joint association analysis of 83,717 human participants. In this report, we describe the imputation of eMERGE results and methods to create the unified imputed merged set of genome-wide variant genotype data. We imputed the data using the Michigan Imputation Server, which provides a missing single-nucleotide variant genotype imputation service using the minimac3 imputation algorithm with the Haplotype Reference Consortium genotype reference set. We describe the quality control and filtering steps used in the generation of this data set and suggest generalizable quality thresholds for imputation and phenotype association studies. To test the merged imputed genotype set, we replicated a previously reported chromosome 6 HLA-B herpes zoster (shingles) association and discovered a novel zoster-associated loci in an epigenetic binding site near the terminus of chromosome 3 (3p29).

Probing the Virtual Proteome to Identify Novel Disease Biomarkers.

BACKGROUND: Proteomic approaches allow measurement of thousands of proteins in a single specimen, which can accelerate biomarker discovery. However, applying these technologies to massive biobanks is not currently feasible because of the practical barriers and costs of implementing such assays at scale. To overcome these challenges, we used a "virtual proteomic" approach, linking genetically predicted protein levels to clinical diagnoses in >40 000 individuals. METHODS: We used genome-wide association data from the Framingham Heart Study (n=759) to construct genetic predictors for 1129 plasma protein levels. We validated the genetic predictors for 268 proteins and used them to compute predicted protein levels in 41 288 genotyped individuals in the Electronic Medical Records and Genomics (eMERGE) cohort. We tested associations for each predicted protein with 1128 clinical phenotypes. Lead associations were validated with directly measured protein levels and either low-density lipoprotein cholesterol or subclinical atherosclerosis in the MDCS (Malmö Diet and Cancer Study; n=651). RESULTS: In the virtual proteomic analysis in eMERGE, 55 proteins were associated with 89 distinct diagnoses at a false discovery rate q<0.1. Among these, 13 associations involved lipid (n=7) or atherosclerosis (n=6) phenotypes. We tested each association for validation in MDCS using directly measured protein levels. At Bonferroni-adjusted significance thresholds, levels of apolipoprotein E isoforms were associated with hyperlipidemia, and circulating C-type lectin domain family 1 member B and platelet-derived growth factor receptor-ß predicted subclinical atherosclerosis. Odds ratios for carotid atherosclerosis were 1.31 (95% CI, 1.08-1.58; P=0.006) per 1-SD increment in C-type lectin domain family 1 member B and 0.79 (0.66-0.94; P=0.008) per 1-SD increment in platelet-derived growth factor receptor-ß. CONCLUSIONS: We demonstrate a biomarker discovery paradigm to identify candidate biomarkers of cardiovascular and other diseases.

Longitudinal Patterns of Occurrence and Remission of Erectile Dysfunction in Men With Type 1 Diabetes.

BACKGROUND: Men with diabetes are at greater risk of erectile dysfunction (ED). AIM: To describe the natural history of ED in men with type 1 diabetes. METHODS: We examined up to 30 years of prospectively collected annual ED status and demographic and clinical variables from 600 male participants in the Diabetes Control and Complications Trial (DCCT; 1983-1993) and its follow-up study, the Epidemiology of Diabetes Interventions and Complications (1994-present; data in this study are through 2012). OUTCOMES: Yes vs no response to whether the participant had experienced impotence in the past year and whether he had used ED medication. RESULTS: Sixty-one percent of men reported ED at least once during the study. For some men, the initial report of ED was permanent. For others, potency returned and was lost multiple times. Visual display of the data showed four longitudinal ED phenotypes: never (38.7%), isolated (6.7%), intermittent (41.8%), and persistent (12.8%). Men who never reported ED or in only 1 isolated year were younger, had lower body mass index, and better glycemic control than men in the intermittent and persistent groups at DCCT baseline. In a multivariable logistic model comparing men at their first year reporting ED, men who were older had lower odds of remission and men who were in the conventional DCCT treatment group had higher odds of remission. CLINICAL TRANSLATION: If validated in other cohorts, such findings could be used to guide individualized interventions for patients with ED. STRENGTHS AND LIMITATIONS: This is the first examination of ED with repeated measures at an annual resolution, with up to 30 years of responses for each participant. However, the yes vs no response is a limitation because the real phenotype is not binary and the question can be interpreted differently depending on the participant. CONCLUSIONS: Age, glycemic control, and BMI were important longitudinal predictors of ED. We have described a more complex ED phenotype, with variation in remission patterns, which could offer insight into different mechanisms or opportunities for intervention. If validated in other cohorts, such findings could be used to establish more accurate prognostication of outcomes for patients with ED to guide individualized interventions. Palmer MR, Holt SK, Sarma AV, et al. Longitudinal Patterns of Occurrence and Remission of Erectile Dysfunction in Men With Type 1 Diabetes. J Sex Med 2017;14:1187-1194.

Mass spectrometry and next-generation sequencing reveal an abundant and rapidly evolving abalone sperm protein.

Abalone, a broadcast spawning marine mollusk, is an important model for molecular interactions and positive selection in fertilization, but the focus has previously been on only two sperm proteins, lysin and sp18. We used genomic and proteomic techniques to bring new insights to this model by characterizing the testis transcriptome and sperm proteome of the Red abalone Haliotis rufescens. One pair of homologous, testis-specific proteins contains a secretion signal and is small, abundant, and associated with the acrosome. Comparative analysis revealed that homologs are extremely divergent between species, and show strong evidence for positive selection. The acrosomal localization and rapid evolution of these proteins indicates that they play an important role in fertilization, and could be involved in the species-specificity of sperm-egg interactions in abalone. Our genomic and proteomic characterization of abalone fertilization resulted in the identification of interesting, novel peptides that have eluded detection in this important model system for 20 years.

Risk factors for orgasmic and concomitant erectile dysfunction in men with type 1 diabetes: a cross-sectional study.

In this study, we sought to determine the burden and characteristics of orgasmic dysfunction (OD) and concomitant erectile dysfunction (ED) in men with type 1 diabetes (T1D) enrolled in the Epidemiology of Diabetes Interventions and Complications (EDIC) study. In 2010, we assessed orgasmic and erectile function using the International Index of Erectile Function (IIEF). Sociodemographic, clinical, and diabetes characteristics were compared by OD status (OD only, OD and ED, no ED or OD). Age-adjusted associations between risk factors and OD status were examined. OD and ED information was available from 563 men. Eighty-three men (14.7%) reported OD of whom 21 reported OD only and 62 reported OD and ED. Age-adjusted odds ratios demonstrated that men who reported OD only had higher odds of depression, low sexual desire, and decreased alcohol use compared with men reporting no dysfunction. Men with OD concomitant with ED had greater odds of elevated hemoglobin A1C, peripheral and autonomic neuropathy, and nephropathy. Men reporting both dysfunctions were also more likely to report smoking, lower urinary tract symptoms, and had greater odds of androgen deficiency than men with no sexual dysfunction. Men with longstanding T1D suffer from an increased burden of OD. Psychogenic factors predominate in men reporting OD only while men who present with concomitant ED report increased burden of diabetes severity, characteristics previously observed with incident ED. ED may be the central impediment to sexual function in men with OD and ED. Longitudinal studies to characterize OD and ED experience over time are warranted.

CNV Association of Diverse Clinical Phenotypes from eMERGE reveals novel disease biology underlying cardiovascular disease.

BACKGROUND: Cardiovascular disease is the leading cause of death in the United States. Consequently, individuals who are genetically predisposed for high risk of cardiovascular disease would benefit most from prevention and early intervention approaches. Among common health risk factors affecting adult populations, we evaluated 23 cardiovascular disease-related traits, including BMI, glucose levels and lipid profiling to determine their associations with low-frequency recurrent copy number variations (CNV) (population frequency < 5%). RESULTS: We examined 10,619 unrelated subjects of European ancestry from the Electronic Medical Records and Genomics (eMERGE) Network who were genotyped with 657,366 markers genome-wide on the Illumina Infinium Quad 660 array. We performed CNV calling based on array marker intensity and evaluated data quality, ancestry stratification, and relatedness to ensure unbiased association discovery. Using a segment-based scoring approach, we assessed the association of all CNVs with each trait. In this large genome-wide analysis of low-frequency CNVs, we observed 11 novel genome-wide significant associations of low-frequency CNVs with major cardiovascular disease traits. CONCLUSION: In one of the largest genome-wide studies for low-frequency recurrent CNVs, we identified 11 loci associated with cardiovascular disease and related traits at the genome-wide significance level that may serve as biomarkers for prevention and early intervention studies in subjects who are at elevated risk. Our study further supports the role of low-frequency recurrent CNVs in the pathogenesis of common complex disease traits.

Insight into the genetic architecture of back pain and its risk factors from a study of 509,000 individuals.

Back pain (BP) is a common condition of major social importance and poorly understood pathogenesis. Combining data from the UK Biobank and CHARGE consortium cohorts allowed us to perform a very large genome-wide association study (total N = 509,070) and examine the genetic correlation and pleiotropy between BP and its clinical and psychosocial risk factors. We identified and replicated 3 BP-associated loci, including one novel region implicating SPOCK2/CHST3 genes. We provide evidence for pleiotropic effects of genetic factors underlying BP, height, and intervertebral disk problems. We also identified independent genetic correlations between BP and depression symptoms, neuroticism, sleep disturbance, overweight, and smoking. A significant enrichment for genes involved in the central nervous system and skeletal tissue development was observed. The study of pleiotropy and genetic correlations, supported by the pathway analysis, suggests at least 2 strong molecular axes of BP genesis, one related to structural/anatomical factors such as intervertebral disk problems and anthropometrics, and another related to the psychological component of pain perception and pain processing. These findings corroborate with the current biopsychosocial model as a paradigm for BP. Overall, the results demonstrate BP to have an extremely complex genetic architecture that overlaps with the genetic predisposition to its biopsychosocial risk factors. The work sheds light on pathways of relevance in the prevention and management of low BP.

A study paradigm integrating prospective epidemiologic cohorts and electronic health records to identify disease biomarkers.

Defining the full spectrum of human disease associated with a biomarker is necessary to advance the biomarker into clinical practice. We hypothesize that associating biomarker measurements with electronic health record (EHR) populations based on shared genetic architectures would establish the clinical epidemiology of the biomarker. We use Bayesian sparse linear mixed modeling to calculate SNP weightings for 53 biomarkers from the Atherosclerosis Risk in Communities study. We use the SNP weightings to computed predicted biomarker values in an EHR population and test associations with 1139 diagnoses. Here we report 116 associations meeting a Bonferroni level of significance. A false discovery rate (FDR)-based significance threshold reveals more known and undescribed associations across a broad range of biomarkers, including biometric measures, plasma proteins and metabolites, functional assays, and behaviors. We confirm an inverse association between LDL-cholesterol level and septicemia risk in an independent epidemiological cohort. This approach efficiently discovers biomarker-disease associations.