Late life disability is a highly devastating condition affecting 20% or more of persons aged 65 years and older in the USA; it is an important determinant of acute medical and long-term care costs which represent a growing burden on national economies. Disability is a multifactorial trait that contributes substantially to decline of health/wellbeing. Accordingly, gaining insights into the genetics of disability could help in identifying molecular mechanisms of this devastating condition and age-related processes contributing to a large fraction of specific geriatric conditions, concordantly with geroscience. We performed a genome-wide association study of disability in a sample of 24,068 subjects from five studies with 12,550 disabled individuals. We identified 30 promising disability-associated polymorphisms in 19 loci at p < 10-4; four of them attained suggestive significance, p < 10-5. In contrast, polygenic risk scores aggregating effects of minor alleles of independent SNPs that were adversely or beneficially associated with disability showed highly significant associations in meta-analysis, p = 3.13 × 10-45 and p = 5.60 × 10-23, respectively, and were replicated in each study. The analysis of genetic pathways, related diseases, and biological functions supported the connections of genes for the identified SNPs with disabling and age-related conditions primarily through oxidative/nitrosative stress, inflammatory response, and ciliary signaling. We identified musculoskeletal system development, maintenance, and regeneration as important components of gene functions. The beneficial and adverse gene sets may be differently implicated in the development of musculoskeletal-related disability with the beneficial set characterized, e.g., by regulation of chondrocyte proliferation and bone formation, and the adverse set by inflammation and bone loss.
Importance: Primary open-angle glaucoma presents with increased prevalence and a higher degree of clinical severity in populations of African ancestry compared with European or Asian ancestry. Despite this, individuals of African ancestry remain understudied in genomic research for blinding disorders. Objectives: To perform a genome-wide association study (GWAS) of African ancestry populations and evaluate potential mechanisms of pathogenesis for loci associated with primary open-angle glaucoma. Design, Settings, and Participants: A 2-stage GWAS with a discovery data set of 2320 individuals with primary open-angle glaucoma and 2121 control individuals without primary open-angle glaucoma. The validation stage included an additional 6937 affected individuals and 14 917 unaffected individuals using multicenter clinic- and population-based participant recruitment approaches. Study participants were recruited from Ghana, Nigeria, South Africa, the United States, Tanzania, Britain, Cameroon, Saudi Arabia, Brazil, the Democratic Republic of the Congo, Morocco, Peru, and Mali from 2003 to 2018. Individuals with primary open-angle glaucoma had open iridocorneal angles and displayed glaucomatous optic neuropathy with visual field defects. Elevated intraocular pressure was not included in the case definition. Control individuals had no elevated intraocular pressure and no signs of glaucoma. Exposures: Genetic variants associated with primary open-angle glaucoma. Main Outcomes and Measures: Presence of primary open-angle glaucoma. Genome-wide significance was defined as P < 5 × 10-8 in the discovery stage and in the meta-analysis of combined discovery and validation data. Results: A total of 2320 individuals with primary open-angle glaucoma (mean [interquartile range] age, 64.6 [56-74] years; 1055 [45.5%] women) and 2121 individuals without primary open-angle glaucoma (mean [interquartile range] age, 63.4 [55-71] years; 1025 [48.3%] women) were included in the discovery GWAS. The GWAS discovery meta-analysis demonstrated association of variants at amyloid-β A4 precursor protein-binding family B member 2 (APBB2; chromosome 4, rs59892895T>C) with primary open-angle glaucoma (odds ratio [OR], 1.32 [95% CI, 1.20-1.46]; P = 2 × 10-8). The association was validated in an analysis of an additional 6937 affected individuals and 14 917 unaffected individuals (OR, 1.15 [95% CI, 1.09-1.21]; P < .001). Each copy of the rs59892895*C risk allele was associated with increased risk of primary open-angle glaucoma when all data were included in a meta-analysis (OR, 1.19 [95% CI, 1.14-1.25]; P = 4 × 10-13). The rs59892895*C risk allele was present at appreciable frequency only in African ancestry populations. In contrast, the rs59892895*C risk allele had a frequency of less than 0.1% in individuals of European or Asian ancestry. Conclusions and Relevance: In this genome-wide association study, variants at the APBB2 locus demonstrated differential association with primary open-angle glaucoma by ancestry. If validated in additional populations this finding may have implications for risk assessment and therapeutic strategies.
Women of African ancestry have lower incidence of epithelial ovarian cancer (EOC) yet worse survival compared to women of European ancestry. We conducted a genome-wide association study in African ancestry women with 755 EOC cases, including 537 high-grade serous ovarian carcinomas (HGSOC) and 1,235 controls. We identified four novel loci with suggestive evidence of association with EOC (p < 1 × 10-6 ), including rs4525119 (intronic to AKR1C3), rs7643459 (intronic to LOC101927394), rs4286604 (12 kb 3' of UGT2A2) and rs142091544 (5 kb 5' of WWC1). For HGSOC, we identified six loci with suggestive evidence of association including rs37792 (132 kb 5' of follistatin [FST]), rs57403204 (81 kb 3' of MAGEC1), rs79079890 (LOC105376360 intronic), rs66459581 (5 kb 5' of PRPSAP1), rs116046250 (GABRG3 intronic) and rs192876988 (32 kb 3' of GK2). Among the identified variants, two are near genes known to regulate hormones and diseases of the ovary (AKR1C3 and FST), and two are linked to cancer (AKR1C3 and MAGEC1). In follow-up studies of the 10 identified variants, the GK2 region SNP, rs192876988, showed an inverse association with EOC in European ancestry women (p = 0.002), increased risk of ER positive breast cancer in African ancestry women (p = 0.027) and decreased expression of GK2 in HGSOC tissue from African ancestry women (p = 0.004). A European ancestry-derived polygenic risk score showed positive associations with EOC and HGSOC in women of African ancestry suggesting shared genetic architecture. Our investigation presents evidence of variants for EOC shared among European and African ancestry women and identifies novel EOC risk loci in women of African ancestry.
BACKGROUND: Genome-wide association studies (GWAS) in Parkinson's disease have increased the scope of biological knowledge about the disease over the past decade. We aimed to use the largest aggregate of GWAS data to identify novel risk loci and gain further insight into the causes of Parkinson's disease. METHODS: We did a meta-analysis of 17 datasets from Parkinson's disease GWAS available from European ancestry samples to nominate novel loci for disease risk. These datasets incorporated all available data. We then used these data to estimate heritable risk and develop predictive models of this heritability. We also used large gene expression and methylation resources to examine possible functional consequences as well as tissue, cell type, and biological pathway enrichments for the identified risk factors. Additionally, we examined shared genetic risk between Parkinson's disease and other phenotypes of interest via genetic correlations followed by Mendelian randomisation. FINDINGS: Between Oct 1, 2017, and Aug 9, 2018, we analysed 7·8 million single nucleotide polymorphisms in 37 688 cases, 18 618 UK Biobank proxy-cases (ie, individuals who do not have Parkinson's disease but have a first degree relative that does), and 1·4 million controls. We identified 90 independent genome-wide significant risk signals across 78 genomic regions, including 38 novel independent risk signals in 37 loci. These 90 variants explained 16-36% of the heritable risk of Parkinson's disease depending on prevalence. Integrating methylation and expression data within a Mendelian randomisation framework identified putatively associated genes at 70 risk signals underlying GWAS loci for follow-up functional studies. Tissue-specific expression enrichment analyses suggested Parkinson's disease loci were heavily brain-enriched, with specific neuronal cell types being implicated from single cell data. We found significant genetic correlations with brain volumes (false discovery rate-adjusted p=0·0035 for intracranial volume, p=0·024 for putamen volume), smoking status (p=0·024), and educational attainment (p=0·038). Mendelian randomisation between cognitive performance and Parkinson's disease risk showed a robust association (p=8·00 × 10-7). INTERPRETATION: These data provide the most comprehensive survey of genetic risk within Parkinson's disease to date, to the best of our knowledge, by revealing many additional Parkinson's disease risk loci, providing a biological context for these risk factors, and showing that a considerable genetic component of this disease remains unidentified. These associations derived from European ancestry datasets will need to be followed-up with more diverse data. FUNDING: The National Institute on Aging at the National Institutes of Health (USA), The Michael J Fox Foundation, and The Parkinson's Foundation (see appendix for full list of funding sources).
Parkinson's disease or first degree relation to individual with Parkinson's disease
Anxiety disorders are common, complex psychiatric disorders with twin heritabilities of 30-60%. We conducted a genome-wide association study of Lifetime Anxiety Disorder (ncase = 25 453, ncontrol = 58 113) and an additional analysis of Current Anxiety Symptoms (ncase = 19 012, ncontrol = 58 113). The liability scale common variant heritability estimate for Lifetime Anxiety Disorder was 26%, and for Current Anxiety Symptoms was 31%. Five novel genome-wide significant loci were identified including an intergenic region on chromosome 9 that has previously been associated with neuroticism, and a locus overlapping the BDNF receptor gene, NTRK2. Anxiety showed significant positive genetic correlations with depression and insomnia as well as coronary artery disease, mirroring findings from epidemiological studies. We conclude that common genetic variation accounts for a substantive proportion of the genetic architecture underlying anxiety.
Genomic studies in African populations provide unique opportunities to understand disease etiology, human diversity, and population history. In the largest study of its kind, comprising genome-wide data from 6,400 individuals and whole-genome sequences from 1,978 individuals from rural Uganda, we find evidence of geographically correlated fine-scale population substructure. Historically, the ancestry of modern Ugandans was best represented by a mixture of ancient East African pastoralists. We demonstrate the value of the largest sequence panel from Africa to date as an imputation resource. Examining 34 cardiometabolic traits, we show systematic differences in trait heritability between European and African populations, probably reflecting the differential impact of genes and environment. In a multi-trait pan-African GWAS of up to 14,126 individuals, we identify novel loci associated with anthropometric, hematological, lipid, and glycemic traits. We find that several functionally important signals are driven by Africa-specific variants, highlighting the value of studying diverse populations across the region.
OBJECTIVE: To study how genetics may play a role in determining risk of chemotherapy-related amenorrhea (CRA) in young women with breast cancer. DESIGN: Genome-wide association study. SETTING: Not applicable. PATIENT(S): Premenopausal women ≤45 years of age enrolled in one of these three trials were included if they had at least one menstrual case report form after chemotherapy ended and if they were of European ancestry. Forms during and up to 3 months after receipt of GnRH agonist were excluded. INTERVENTION(S): None. MAIN OUTCOME MEASURE(S): The association of single-nucleotide polymorphisms with post-chemotherapy menstruation adjusted for trial and arm, age, tamoxifen use, and nodal status. RESULT(S): The median age of the 1,168 women was 41 years (range 19-45). Among these, 457 (39%) never resumed menses after chemotherapy. Older age, tamoxifen use, and node-negative disease were associated with increased risk of CRA. Adjusting for these, rs147451859, in an intron of PPCDC (phosphopantothenoylcysteine decarboxylase), and rs17587029, located 5' upstream of RPS20P11 (ribosomal protein S20 pseudogene 11), were associated with post-chemotherapy menstruation. CONCLUSION(S): Genetic variation may contribute to risk of CRA. Better prediction of who will experience CRA may inform reproductive and treatment decision making in young women with cancer.
Adverse response to chemotherapy (amenorrhea) in breast cancer
Aortic calcification is an important independent predictor of future cardiovascular events. We performed a genome-wide association meta-analysis to determine SNPs associated with the extent of abdominal aortic calcification (n = 9,417) or descending thoracic aortic calcification (n = 8,422). Two genetic loci, HDAC9 and RAP1GAP, were associated with abdominal aortic calcification at a genome-wide level (P < 5.0 × 10-8). No SNPs were associated with thoracic aortic calcification at the genome-wide threshold. Increased expression of HDAC9 in human aortic smooth muscle cells promoted calcification and reduced contractility, while inhibition of HDAC9 in human aortic smooth muscle cells inhibited calcification and enhanced cell contractility. In matrix Gla protein-deficient mice, a model of human vascular calcification, mice lacking HDAC9 had a 40% reduction in aortic calcification and improved survival. This translational genomic study identifies the first genetic risk locus associated with calcification of the abdominal aorta and describes a previously unknown role for HDAC9 in the development of vascular calcification.
Creatine kinase (CK) and lactate dehydrogenase (LDH) are widely used markers of tissue damage. To search for sequence variants influencing serum levels of CK and LDH, 28.3 million sequence variants identified through whole-genome sequencing of 2,636 Icelanders were imputed into 63,159 and 98,585 people with CK and LDH measurements, respectively. Here we describe 13 variants associating with serum CK and 16 with LDH levels, including four that associate with both. Among those, 15 are non-synonymous variants and 12 have a minor allele frequency below 5%. We report sequence variants in genes encoding the enzymes being measured (CKM and LDHA), as well as in genes linked to muscular (ANO5) and immune/inflammatory function (CD163/CD163L1, CSF1, CFH, HLA-DQB1, LILRB5, NINJ1 and STAB1). A number of the genes are linked to the mononuclear/phagocyte system and clearance of enzymes from the serum. This highlights the variety in the sources of normal diversity in serum levels of enzymes.
Because metabolites are hypothesized to play key roles as markers and effectors of cardiometabolic diseases, recent studies have sought to annotate the genetic determinants of circulating metabolite levels. We report a genome-wide association study (GWAS) of 217 plasma metabolites, including >100 not measured in prior GWAS, in 2076 participants of the Framingham Heart Study (FHS). For the majority of analytes, we find that estimated heritability explains >20% of interindividual variation, and that variation attributable to heritable factors is greater than that attributable to clinical factors. Further, we identify 31 genetic loci associated with plasma metabolites, including 23 that have not previously been reported. Importantly, we include GWAS results for all surveyed metabolites and demonstrate how this information highlights a role for AGXT2 in cholesterol ester and triacylglycerol metabolism. Thus, our study outlines the relative contributions of inherited and clinical factors on the plasma metabolome and provides a resource for metabolism research.
Rationale: Idiopathic pulmonary fibrosis (IPF) is a complex lung disease characterized by scarring of the lung that is believed to result from an atypical response to injury of the epithelium. Genome-wide association studies have reported signals of association implicating multiple pathways including host defense, telomere maintenance, signaling, and cell-cell adhesion.Objectives: To improve our understanding of factors that increase IPF susceptibility by identifying previously unreported genetic associations.Methods: We conducted genome-wide analyses across three independent studies and meta-analyzed these results to generate the largest genome-wide association study of IPF to date (2,668 IPF cases and 8,591 controls). We performed replication in two independent studies (1,456 IPF cases and 11,874 controls) and functional analyses (including statistical fine-mapping, investigations into gene expression, and testing for enrichment of IPF susceptibility signals in regulatory regions) to determine putatively causal genes. Polygenic risk scores were used to assess the collective effect of variants not reported as associated with IPF.Measurements and Main Results: We identified and replicated three new genome-wide significant (P < 5 × 10-8) signals of association with IPF susceptibility (associated with altered gene expression of KIF15, MAD1L1, and DEPTOR) and confirmed associations at 11 previously reported loci. Polygenic risk score analyses showed that the combined effect of many thousands of as yet unreported IPF susceptibility variants contribute to IPF susceptibility.Conclusions: The observation that decreased DEPTOR expression associates with increased susceptibility to IPF supports recent studies demonstrating the importance of mTOR signaling in lung fibrosis. New signals of association implicating KIF15 and MAD1L1 suggest a possible role of mitotic spindle-assembly genes in IPF susceptibility.
A significantly higher proportion of patients with nonalcoholic steatohepatitis (NASH) who received obeticholic acid (OCA) had histological improvement relative to placebo in the FLINT (farnesoid X nuclear receptor ligand obeticholic acid for noncirrhotic, NASH treatment) trial. However, genetic predictors of response to OCA are unknown. We conducted a genome-wide association study (GWAS) in FLINT participants to identify variants associated with NASH resolution and fibrosis improvement. Genotyping was performed using the Omni2.5 content GWAS chip. To avoid false positives introduced by population stratification, we focused our GWAS on white participants. Six regions on chromosomes 1, 4, 6, 7, 15, and 17 had multiple single nucleotide polymorphisms (SNPs) with suggestive association (P < 1 × 10-4 ) with NASH resolution. A sentinel SNP, rs75508464, near CELA3B on chromosome 1 was associated with NASH resolution, improvement in the nonalcoholic fatty liver disease activity score, portal inflammation, and fibrosis. Among individuals carrying this allele, 83% achieved NASH resolution with OCA compared with only 33% with placebo. Eight regions on chromosomes 1, 2, 3, 11, 13, and 18 had multiple SNPs associated with fibrosis improvement; of these, rs12130403 near TDRD10 on chromosome 1 was also associated with improvement in NASH and portal inflammation, and rs4073431 near ANO3 on chromosome 11 was associated with NASH resolution and improvement in steatosis. Multiple SNPs on chromosome 11 had suggestive association with pruritus, with rs1379650 near ANO5 being the top SNP. Conclusion: We identified several variants that may be associated with histological improvement and pruritus in individuals with NASH receiving OCA. The rs75508464 variant near CELA3B may have the most significant effect on NASH resolution in those receiving OCA.
Fibrosis improvement in nonalcoholic steatohepatitis
BACKGROUND: Metabolism of inorganic arsenic (iAs) is subject to inter-individual variability, which is explained partly by genetic determinants. OBJECTIVES: We investigated the association of genetic variants with arsenic species and principal components of arsenic species in the Strong Heart Family Study (SHFS). METHODS: We examined variants previously associated with cardiometabolic traits (~ 200,000 from Illumina Cardio MetaboChip) or arsenic metabolism and toxicity (670) among 2,428 American Indian participants in the SHFS. Urine arsenic species were measured by high performance liquid chromatography-inductively coupled plasma mass spectrometry (HPLC-ICP-MS), and percent arsenic species [iAs, monomethylarsonate (MMA), and dimethylarsinate (DMA), divided by their sum × 100] were logit transformed. We created two orthogonal principal components that summarized iAs, MMA, and DMA and were also phenotypes for genetic analyses. Linear regression was performed for each phenotype, dependent on allele dosage of the variant. Models accounted for familial relatedness and were adjusted for age, sex, total arsenic levels, and population stratification. Single nucleotide polymorphism (SNP) associations were stratified by study site and were meta-analyzed. Bonferroni correction was used to account for multiple testing. RESULTS: Variants at 10q24 were statistically significant for all percent arsenic species and principal components of arsenic species. The index SNP for iAs%, MMA%, and DMA% (rs12768205) and for the principal components (rs3740394, rs3740393) were located near AS3MT, whose gene product catalyzes methylation of iAs to MMA and DMA. Among the candidate arsenic variant associations, functional SNPs in AS3MT and 10q24 were most significant (p < 9.33 × 10-5). CONCLUSIONS: This hypothesis-driven association study supports the role of common variants in arsenic metabolism, particularly AS3MT and 10q24. Citation: Balakrishnan P, Vaidya D, Franceschini N, Voruganti VS, Gribble MO, Haack K, Laston S, Umans JG, Francesconi KA, Goessler W, North KE, Lee E, Yracheta J, Best LG, MacCluer JW, Kent J Jr., Cole SA, Navas-Acien A. 2017. Association of cardiometabolic genes with arsenic metabolism biomarkers in American Indian communities: the Strong Heart Family Study (SHFS). Environ Health Perspect 125:15-22; http://dx.doi.org/10.1289/EHP251.
The triggering receptor expressed on myeloid cells 2 (TREM2) gene has been reported to increase the risk of Alzheimer's disease (AD). The soluble TREM2 protein (sTREM2) in cerebrospinal fluid (CSF) was also associated with AD. However, the role of sTREM2 in AD and its genetic modifiers remain unclear. We carried out a genome-wide association study for CSF sTREM2 levels using participants from the Alzheimer's Disease Neuroimaging Initiative and validated the significant association in an independent cohort from Chinese Alzheimer's Biomarker and LifestylE study. rs7232 in membrane spanning 4-domains A6A (MS4A6A) gene was associated with CSF sTREM2 levels at genome-wide significance (p = 1.42 × 10-15). The locus influences CSF sTREM2 levels especially in nondemented individuals. And the association was replicable in the validation cohort from Chinese Alzheimer's Biomarker and LifestylE study (p = 0.0106). Besides, the expressions of MS4A6A and TREM2 were correlated in brain regions (p < 2 × 10-16). The findings of our study suggest that the AD risk variant in the MS4A6A gene participates in the regulation of sTREM2.
BACKGROUND: Previous studies of radiological damage in rheumatoid arthritis (RA) have used candidate-gene approaches, or evaluated single genome-wide association studies (GWAS). We undertook the first meta-analysis of GWAS of RA radiological damage to: (1) identify novel genetic loci for this trait; and (2) test previously validated variants. METHODS: Seven GWAS (2,775 RA cases, of a range of ancestries) were combined in a meta-analysis. Radiological damage was assessed using modified Larsen scores, Sharp van Der Heijde scores, and erosive status. Single nucleotide polymophsim (SNP) associations with radiological damage were tested at a single time-point using regression models. Primary analyses included age and disease duration as covariates. Secondary analyses also included rheumatoid factor (RF). Meta-analyses were undertaken in trans-ethnic and European-only cases. RESULTS: In the trans-ethnic primary meta-analysis, one SNP (rs112112734) in close proximity to HLA-DRB1, and strong linkage disequilibrium with the shared-epitope, attained genome-wide significance (P = 4.2x10-8). In the secondary analysis (adjusting for RF) the association was less significant (P = 1.7x10-6). In both trans-ethnic primary and secondary meta-analyses 14 regions contained SNPs with associations reaching P<5x10-6; in the European primary and secondary analyses 13 and 10 regions contained SNPs reaching P<5x10-6, respectively. Of the previously validated SNPs for radiological progression, only rs660895 (tagging HLA-DRB1*04:01) attained significance (P = 1.6x10-5) and had a consistent direction of effect across GWAS. CONCLUSIONS: Our meta-analysis confirms the known association between the HLA-DRB1 shared epitope and RA radiological damage. The lack of replication of previously validated non-HLA markers highlights a requirement for further research to deliver clinically-useful prognostic genetic markers.
Infections and mental disorders are two of the major global disease burdens. While correlations between mental disorders and infections have been reported, the possible genetic links between them have not been assessed in large-scale studies. Moreover, the genetic basis of susceptibility to infection is largely unknown, as large-scale genome-wide association studies of susceptibility to infection have been lacking. We utilized a large Danish population-based sample (N = 65,534) linked to nationwide population-based registers to investigate the genetic architecture of susceptibility to infection (heritability estimation, polygenic risk analysis, and a genome-wide association study (GWAS)) and examined its association with mental disorders (comorbidity analysis and genetic correlation). We found strong links between having at least one psychiatric diagnosis and the occurrence of infection (P = 2.16 × 10-208, OR = 1.72). The SNP heritability of susceptibility to infection ranged from ~2 to ~7% in samples of differing psychiatric diagnosis statuses (suggesting the environment as a major contributor to susceptibility), and polygenic risk scores moderately but significantly explained infection status in an independent sample. We observed a genetic correlation of 0.496 (P = 2.17 × 10-17) between a diagnosis of infection and a psychiatric diagnosis. While our GWAS did not identify genome-wide significant associations, we found 90 suggestive (P ≤ 10-5) associations for susceptibility to infection. Our findings suggest a genetic component in susceptibility to infection and indicate that the occurrence of infections in individuals with mental illness may be in part genetically driven.
Elevated serum urate levels cause gout and correlate with cardiometabolic diseases via poorly understood mechanisms. We performed a trans-ancestry genome-wide association study of serum urate in 457,690 individuals, identifying 183 loci (147 previously unknown) that improve the prediction of gout in an independent cohort of 334,880 individuals. Serum urate showed significant genetic correlations with many cardiometabolic traits, with genetic causality analyses supporting a substantial role for pleiotropy. Enrichment analysis, fine-mapping of urate-associated loci and colocalization with gene expression in 47 tissues implicated the kidney and liver as the main target organs and prioritized potentially causal genes and variants, including the transcriptional master regulators in the liver and kidney, HNF1A and HNF4A. Experimental validation showed that HNF4A transactivated the promoter of ABCG2, encoding a major urate transporter, in kidney cells, and that HNF4A p.Thr139Ile is a functional variant. Transcriptional coregulation within and across organs may be a general mechanism underlying the observed pleiotropy between urate and cardiometabolic traits.
Bipolar disorder is a common, chronic psychiatric disorder. Despite high heritability, there is a paucity of identified genetic risk factors. Immune biomarkers are under more direct genetic influence than bipolar disorder. To explore the genetic associations with immune biomarker levels in cerebrospinal fluid (CSF) and blood serum which previously showed differences in bipolar disorder, we performed a study involving 291 individuals (184 bipolar disorder patients and 107 controls). The biomarkers assayed in both CSF and serum were: chitinase-3-like protein-1 (YKL-40), monocyte chemoattractant protein-1 (MCP-1), soluble cluster of differentiation (sCD14), tissue inhibitor of metalloproteinases-1 and 2 (TIMP-1 and TIMP-2). C-reactive protein (CRP) was only quantified in serum, and interleukin 8 (IL-8) measures were only available in CSF. Genome-wide association studies were conducted using PLINK for each of three genotyping waves and incorporated covariates for population substructure, age, sex, and body mass index (BMI). Results were combined by meta-analysis. Genome-wide significant associations were detected for all biomarkers except TIMP-1 and TIMP-2 in CSF. The strongest association in CSF was found for markers within the CNTNAP5 gene with YKL-40 (rs150248456, P = 2.84 × 10-10). The strongest association in serum was also for YKL-40 but localized to the FANCI gene (rs188263039, P = 5.80 × 10-26). This study revealed numerous biologically plausible genetic associations with immune biomarkers in CSF and blood serum. Importantly, the genetic variants regulating immune biomarker levels in CSF and blood serum differ. These results extend our knowledge of how biomarkers showing alterations in bipolar disorder are genetically regulated.
OBJECTIVE: A number of epidemiological studies have reported that decreased serum bilirubin, an endogenous antioxidant, is associated with cardiovascular disease. However, previous Mendelian randomization analyses conducted using a single sample have shown no evidence of association. Approach and Results: A 2-sample summary Mendelian randomization study was performed by obtaining exposure and outcome data from separate nonoverlapping samples. We utilized data from the KoGES (Korean Genome and Epidemiology Study; n=25 406) and KCPS-II (Korean Cancer Prevention Study-II; n=14 541) biobank for serum bilirubin and stroke, respectively. Using KoGES, a total of 1784 single nucleotide polymorphisms associated with serum bilirubin levels were discovered using a genome-wide significance threshold (P<5×10-8), of which 10 single nucleotide polymorphisms were identified as independent (R2<0.005) and adopted as genetic instruments. From KCPS-II, total and ischemic stroke cases were identified (n=1489 and n=686), with 12 366 acting as controls. Various 2-sample summary Mendelian randomization methods were employed, with Mendelian randomization estimates showing an inverse causal association between serum bilirubin levels and total stroke risk (odds ratio, 0.481 [95% CI, 0.234-0.988]; P=0.046). This association increased in magnitude when restricting the analysis to ischemic stroke cases (odds ratio, 0.302 [95% CI, 0.105-0.868]; P=0.026). CONCLUSIONS: Our findings provide evidence of significant causal relationship between high levels of bilirubin and decreased stroke risk in Korean population in agreement with observational approaches. This highlights the potential for bilirubin to serve as a therapeutic target for oxidative stress-related diseases such as stroke and suggests that previous findings were not a consequence of unmeasured confounding.
Motivated by observational studies that report associations between schizophrenia and traits, such as poor diet, increased body mass index and metabolic disease, we investigated the genetic contribution to dietary intake in a sample of 335,576 individuals from the UK Biobank study. A principal component analysis applied to diet question item responses generated two components: Diet Component 1 (DC1) represented a meat-related diet and Diet Component 2 (DC2) a fish and plant-related diet. Genome-wide association analysis identified 29 independent single-nucleotide polymorphisms (SNPs) associated with DC1 and 63 SNPs with DC2. Estimated from over 35,000 3rd-degree relative pairs that are unlikely to share close family environments, heritabilities for both DC1 and DC2 were 0.16 (standard error (s.e.) = 0.05). SNP-based heritability was 0.06 (s.e. = 0.003) for DC1 and 0.08 (s.e = 0.004) for DC2. We estimated significant genetic correlations between both DCs and schizophrenia, and several other traits. Mendelian randomisation analyses indicated a negative uni-directional relationship between liability to schizophrenia and tendency towards selecting a meat-based diet (which could be direct or via unidentified correlated variables), but a bi-directional relationship between liability to schizophrenia and tendency towards selecting a fish and plant-based diet consistent with genetic pleiotropy.