A meta-analysis of genome-wide studies of resilience in the German population.

Resilience is the capacity to adapt to stressful life events. As such, this trait is associated with physical and mental functions and conditions. Here, we aimed to identify the genetic factors contributing to shape resilience. We performed variant- and gene-based meta-analyses of genome-wide association studies from six German cohorts (N = 15822) using the 11-item version of the Resilience Scale (RS-11) as outcome measure. Variant- and gene-level results were combined to explore the biological context using network analysis. In addition, we conducted tests of correlation between RS-11 and the polygenic scores (PGSs) for 12 personality and mental health traits in one of these cohorts (PROCAM-2, N = 3879). The variant-based analysis found no signals associated with resilience at the genome-wide level (p < 5 × 10-8), but suggested five genomic loci (p < 1 × 10-5). The gene-based analysis identified three genes (ROBO1, CIB3 and LYPD4) associated with resilience at genome-wide level (p < 2.48 × 10-6) and 32 potential candidates (p < 1 × 10-4). Network analysis revealed enrichment of biological pathways related to neuronal proliferation and differentiation, synaptic organization, immune responses and vascular homeostasis. We also found significant correlations (FDR < 0.05) between RS-11 and the PGSs for neuroticism and general happiness. Overall, our observations suggest low heritability of resilience. Large, international efforts will be required to uncover the genetic factors that contribute to shape trait resilience. Nevertheless, as the largest investigation of the genetics of resilience in general population to date, our study already offers valuable insights into the biology potentially underlying resilience and resilience's relationship with other personality traits and mental health.

Thyroid function, pernicious anemia and erythropoiesis: a two-sample Mendelian randomization study.

Autoimmune thyroid disease (AITD) and pernicious anemia (PA) often coexist, but the directionality is unknown. In a two-sample Mendelian randomization (MR) analysis, using summary statistics from large genome-wide association studies (GWASs) in Europeans (N = 49 269-755 406), we examined the genetic associations between thyroid function, PA and markers of erythropoiesis. We performed inverse variance weighted random-effects MR, several sensitivity MR analyses, and bidirectional MR and MR Steiger for directionality. AITD and PA were associated bidirectionally (P ≤ 8 × 10-6). Neither euthyroid thyroid stimulating hormone (TSH) nor free thyroxine (FT4) were causally associated with PA. One standard deviation (SD) increase in euthyroid FT4 regulated by genetic variants in deiodinases 1 and 2 genes (DIO1/DIO2), corresponding to low-normal free triiodothyronine (FT3) levels, was causally associated with a pernicious/macrocytic anemia pattern, i.e. decreased erythrocyte counts (rank-based inverse normal transformed ß = -0,064 [95% confidence interval: -0,085, -0,044], P = 8 × 10-10) and hemoglobin (-0.028 [-0.051, -0.005], P = 0.02) and increased mean corpuscular hemoglobin (0.058 [0.025, 0.091], P = 5 × 10-4) and mean corpuscular volume levels (0.075 [0.052, 0.098], P = 1 × 10-8). Meanwhile, subclinical hyperthyroidism mirrored that pattern. AITD was causally associated with increased erythrocyte distribution width (P = 0.007) and decreased reticulocyte counts (P ≤ 0.02), whereas high-normal FT4 regulated by DIO1/DIO2 variants was causally associated with decreased bilirubin (-0.039 (-0.064, -0.013), P = 0.003). In conclusion, the bidirectional association between AITD and PA suggests a shared heritability for these two autoimmune diseases. AITD was causally associated with impaired erythropoiesis and not autoimmune hemolysis. Additionally, in euthyroid individuals, local regulation of thyroid hormones by deiodinases likely plays a role in erythropoiesis.

Alpha-1 Antitrypsin Augmentation and the Liver Phenotype of Adults With Alpha-1 Antitrypsin Deficiency (Genotype Pi∗ZZ).

BACKGROUND & AIMS: α1-Antitrypsin (AAT) is a major protease inhibitor produced by hepatocytes. The most relevant AAT mutation giving rise to AAT deficiency (AATD), the 'Pi∗Z' variant, causes harmful AAT protein accumulation in the liver, shortage of AAT in the systemic circulation, and thereby predisposes to liver and lung injury. Although intravenous AAT augmentation constitutes an established treatment of AATD-associated lung disease, its impact on the liver is unknown. METHODS: Liver-related parameters were assessed in a multinational cohort of 760 adults with severe AATD (Pi∗ZZ genotype) and available liver phenotyping, of whom 344 received augmentation therapy and 416 did not. Liver fibrosis was evaluated noninvasively via the serum test AST-to-platelet ratio index and via transient elastography-based liver stiffness measurement. Histologic parameters were compared in 15 Pi∗ZZ adults with and 35 without augmentation. RESULTS: Compared with nonaugmented subjects, augmented Pi∗ZZ individuals displayed lower serum liver enzyme levels (AST 71% vs 75% upper limit of normal, P < .001; bilirubin 49% vs 58% upper limit of normal, P = .019) and lower surrogate markers of fibrosis (AST-to-platelet ratio index 0.34 vs 0.38, P < .001; liver stiffness measurement 6.5 vs 7.2 kPa, P = .005). Among biopsied participants, augmented individuals had less pronounced liver fibrosis and less inflammatory foci but no differences in AAT accumulation were noted. CONCLUSIONS: The first evaluation of AAT augmentation on the Pi∗ZZ-related liver disease indicates liver safety of a widely used treatment for AATD-associated lung disease. Prospective studies are needed to confirm the beneficial effects and to demonstrate the potential efficacy of exogenous AAT in patients with Pi∗ZZ-associated liver disease.

Epigenetic and integrative cross-omics analyses of cerebral white matter hyperintensities on MRI.

Cerebral white matter hyperintensities on MRI are markers of cerebral small vessel disease, a major risk factor for dementia and stroke. Despite the successful identification of multiple genetic variants associated with this highly heritable condition, its genetic architecture remains incompletely understood. More specifically, the role of DNA methylation has received little attention. We investigated the association between white matter hyperintensity burden and DNA methylation in blood at ∼450 000 cytosine-phosphate-guanine (CpG) sites in 9732 middle-aged to older adults from 14 community-based studies. Single CpG and region-based association analyses were carried out. Functional annotation and integrative cross-omics analyses were performed to identify novel genes underlying the relationship between DNA methylation and white matter hyperintensities. We identified 12 single CpG and 46 region-based DNA methylation associations with white matter hyperintensity burden. Our top discovery single CpG, cg24202936 (P = 7.6 × 10-8), was associated with F2 expression in blood (P = 6.4 × 10-5) and co-localized with FOLH1 expression in brain (posterior probability = 0.75). Our top differentially methylated regions were in PRMT1 and in CCDC144NL-AS1, which were also represented in single CpG associations (cg17417856 and cg06809326, respectively). Through Mendelian randomization analyses cg06809326 was putatively associated with white matter hyperintensity burden (P = 0.03) and expression of CCDC144NL-AS1 possibly mediated this association. Differentially methylated region analysis, joint epigenetic association analysis and multi-omics co-localization analysis consistently identified a role of DNA methylation near SH3PXD2A, a locus previously identified in genome-wide association studies of white matter hyperintensities. Gene set enrichment analyses revealed functions of the identified DNA methylation loci in the blood-brain barrier and in the immune response. Integrative cross-omics analysis identified 19 key regulatory genes in two networks related to extracellular matrix organization, and lipid and lipoprotein metabolism. A drug-repositioning analysis indicated antihyperlipidaemic agents, more specifically peroxisome proliferator-activated receptor-alpha, as possible target drugs for white matter hyperintensities. Our epigenome-wide association study and integrative cross-omics analyses implicate novel genes influencing white matter hyperintensity burden, which converged on pathways related to the immune response and to a compromised blood-brain barrier possibly due to disrupted cell-cell and cell-extracellular matrix interactions. The results also suggest that antihyperlipidaemic therapy may contribute to lowering risk for white matter hyperintensities possibly through protection against blood-brain barrier disruption.

Methylation Patterns of the FKBP5 Gene in Association with Childhood Maltreatment and Depressive Disorders.

Childhood maltreatment is an important risk factor for adult depression and has been associated with changes in the hypothalamic pituitary adrenal (HPA) axis, including cortisol secretion and methylation of the FKBP5 gene. Furthermore, associations between depression and HPA changes have been reported. This study investigated the associations of whole-blood FKBP5 mRNA levels, serum cortisol levels, childhood maltreatment, and depressive symptoms with the whole-blood methylation status (assessed via target bisulfite sequencing) of 105 CpGs at the FKBP5 locus using data from the general population-based Study of Health in Pomerania (SHIP) (N = 203). Both direct and interaction effects with the rs1360780 single-nucleotide polymorphism were investigated. Nominally significant associations of main effects on methylation of a single CpG site were observed at intron 3, intron 7, and the 3'-end of the gene. Additionally, methylation at two clusters at the 3'-end and intron 7 were nominally associated with childhood maltreatment × rs1360780 and depressive symptoms × rs1360780, respectively. The results add to the understanding of molecular mechanisms underlying the emergence of depression and could aid the development of personalised depression therapy and drug development.

Genome-wide association study of circulating interleukin 6 levels identifies novel loci.

Interleukin 6 (IL-6) is a multifunctional cytokine with both pro- and anti-inflammatory properties with a heritability estimate of up to 61%. The circulating levels of IL-6 in blood have been associated with an increased risk of complex disease pathogenesis. We conducted a two-staged, discovery and replication meta genome-wide association study (GWAS) of circulating serum IL-6 levels comprising up to 67 428 (ndiscovery = 52 654 and nreplication = 14 774) individuals of European ancestry. The inverse variance fixed effects based discovery meta-analysis, followed by replication led to the identification of two independent loci, IL1F10/IL1RN rs6734238 on chromosome (Chr) 2q14, (Pcombined = 1.8 × 10-11), HLA-DRB1/DRB5 rs660895 on Chr6p21 (Pcombined = 1.5 × 10-10) in the combined meta-analyses of all samples. We also replicated the IL6R rs4537545 locus on Chr1q21 (Pcombined = 1.2 × 10-122). Our study identifies novel loci for circulating IL-6 levels uncovering new immunological and inflammatory pathways that may influence IL-6 pathobiology.

Genetic Variants in WNT16 and PKD2L1 Locus Affect Heel Ultrasound Bone Stiffness: Analyses from the General Population and Patients Evaluated for Osteoporosis.

Osteoporosis, a complex chronic disease with increasing prevalence, is characterised by reduced bone mineral density (BMD) and increased fracture risk. The high heritability of BMD suggests substantial impact of the individual genetic disposition on bone phenotypes and the development of osteoporosis. In the past years, genome-wide association studies (GWAS) identified hundreds of genetic variants associated with BMD or osteoporosis. Here, we analysed 1103 single nucleotide polymorphisms (SNPs), previously identified as associated with estimated BMD (eBMD) in the UK Biobank. We assessed whether these SNPs are related to heel stiffness index obtained by quantitative ultrasound in 5665 adult participants of the Study of Health in Pomerania (SHIP). We confirmed 45 significant associations after correction for multiple testing. Next, we analysed six selected SNPs in 631 patients evaluated for osteoporosis [rs2707518 (CPED1/WNT16), rs3779381 (WNT16), rs115242848 (LOC101927709/EN1), rs10239787 (JAZF1), rs603424 (PKD2L1) and rs6968704 (JAZF1)]. Differences in minor allele frequencies (MAF) of rs2707518 and rs3779381 between SHIP participants (higher MAF) and patients evaluated for osteoporosis (lower MAF) indicated a protective effect of the minor allele on bone integrity. In contrast, differences in MAF of rs603424 indicated a harmful effect. Co-localisation analyses indicated that the rs603424 effect may be mediated via stearoyl-CoA desaturase (SCD) expression, an enzyme highly expressed in adipose tissue with a crucial role in lipogenesis. Taken together, our results support the role of the WNT16 pathway in the regulation of bone properties and indicate a novel causal role of SCD expression in adipose tissue on bone integrity.

Genome-wide interaction study with major depression identifies novel variants associated with cognitive function.

Major Depressive Disorder (MDD) often is associated with significant cognitive dysfunction. We conducted a meta-analysis of genome-wide interaction of MDD and cognitive function using data from four large European cohorts in a total of 3510 MDD cases and 6057 controls. In addition, we conducted analyses using polygenic risk scores (PRS) based on data from the Psychiatric Genomics Consortium (PGC) on the traits of MDD, Bipolar disorder (BD), Schizophrenia (SCZ), and mood instability (MIN). Functional exploration contained gene expression analyses and Ingenuity Pathway Analysis (IPA®). We identified a set of significantly interacting single nucleotide polymorphisms (SNPs) between MDD and the genome-wide association study (GWAS) of cognitive domains of executive function, processing speed, and global cognition. Several of these SNPs are located in genes expressed in brain, with important roles such as neuronal development (REST), oligodendrocyte maturation (TNFRSF21), and myelination (ARFGEF1). IPA® identified a set of core genes from our dataset that mapped to a wide range of canonical pathways and biological functions (MPO, FOXO1, PDE3A, TSLP, NLRP9, ADAMTS5, ROBO1, REST). Furthermore, IPA® identified upstream regulator molecules and causal networks impacting on the expression of dataset genes, providing a genetic basis for further clinical exploration (vitamin D receptor, beta-estradiol, tadalafil). PRS of MIN and meta-PRS of MDD, MIN and SCZ were significantly associated with all cognitive domains. Our results suggest several genes involved in physiological processes for the development and maintenance of cognition in MDD, as well as potential novel therapeutic agents that could be explored in patients with MDD associated cognitive dysfunction.

Allometric body shape indices, type 2 diabetes and kidney function: A two-sample Mendelian randomization study.

AIM: To examine the association between body mass index (BMI)-independent allometric body shape indices and kidney function. MATERIALS AND METHODS: We performed a two-sample Mendelian randomization (MR) analysis, using summary statistics from UK Biobank, CKDGen and DIAGRAM. BMI-independent allometric body shape indices were: A Body Shape Index (ABSI), Waist-Hip Index (WHI) and Hip Index (HI). Kidney function outcomes were: urinary albumin-to-creatinine ratio (UACR), estimated glomerular filtration rate and blood urea nitrogen. Furthermore, we investigated type 2 diabetes (T2D) as a potential mediator on the pathway to albuminuria. The main analysis was inverse variance-weighted random-effects MR in participants of European ancestry. We also performed several sensitivity MR analyses. RESULTS: A 1-standard deviation (SD) increase in genetically predicted ABSI and WHI levels was associated with higher UACR (ß = 0.039 [95% confidence interval: 0.016, 0.063] log [UACR], P = 0.001 for ABSI, and ß = 0.028 [0.012, 0.044] log [UACR], P = 6 x 10-4 for WHI) in women, but not in men. Meanwhile, a 1-SD increase in genetically predicted HI was associated with lower UACR in women (ß = -0.021 [-0.041, 0.000] log [UACR], P = 0.05) and in men (ß = -0.026 [-0.058, 0.005] log [UACR], P = 0.10). Corresponding estimates in individuals with diabetes were substantially augmented. Risk of T2D increased for genetically high ABSI and WHI in women (P < 6 x 10-19 ) only, but decreased for genetically high HI in both sexes (P < 9 x 10-3 ). No other associations were observed. CONCLUSIONS: Genetically high HI was associated with decreased risk of albuminuria, mediated through decreased T2D risk in both sexes. Opposite associations applied to genetically high ABSI and WHI in women only.

Genetic loci and prioritization of genes for kidney function decline derived from a meta-analysis of 62 longitudinal genome-wide association studies.

Estimated glomerular filtration rate (eGFR) reflects kidney function. Progressive eGFR-decline can lead to kidney failure, necessitating dialysis or transplantation. Hundreds of loci from genome-wide association studies (GWAS) for eGFR help explain population cross section variability. Since the contribution of these or other loci to eGFR-decline remains largely unknown, we derived GWAS for annual eGFR-decline and meta-analyzed 62 longitudinal studies with eGFR assessed twice over time in all 343,339 individuals and in high-risk groups. We also explored different covariate adjustment. Twelve genome-wide significant independent variants for eGFR-decline unadjusted or adjusted for eGFR-baseline (11 novel, one known for this phenotype), including nine variants robustly associated across models were identified. All loci for eGFR-decline were known for cross-sectional eGFR and thus distinguished a subgroup of eGFR loci. Seven of the nine variants showed variant-by-age interaction on eGFR cross section (further about 350,000 individuals), which linked genetic associations for eGFR-decline with age-dependency of genetic cross-section associations. Clinically important were two to four-fold greater genetic effects on eGFR-decline in high-risk subgroups. Five variants associated also with chronic kidney disease progression mapped to genes with functional in-silico evidence (UMOD, SPATA7, GALNTL5, TPPP). An unfavorable versus favorable nine-variant genetic profile showed increased risk odds ratios of 1.35 for kidney failure (95% confidence intervals 1.03-1.77) and 1.27 for acute kidney injury (95% confidence intervals 1.08-1.50) in over 2000 cases each, with matched controls). Thus, we provide a large data resource, genetic loci, and prioritized genes for kidney function decline, which help inform drug development pipelines revealing important insights into the age-dependency of kidney function genetics.

Assessment of significance of conditionally independent GWAS signals.

MOTIVATION: Multiple independently associated SNPs within a linkage disequilibrium region are a common phenomenon. Conditional analysis has been successful in identifying secondary signals. While conditional association tests are limited to specific genomic regions, they are benchmarked with genome-wide scale criterion, a conservative strategy. Within the weighted hypothesis testing framework, we developed a 'quasi-adaptive' method that uses the pairwise correlation (r2) and physical distance (d) from the index association to construct priority functions G =G(r2, d), which assign an SNP-specific α-threshold to each SNP. Family-wise error rate (FWER) and power of the approach were evaluated via simulations based on real GWAS data. We compared a series of different G-functions. RESULTS: Simulations under the null hypothesis on 1,100 primary SNPs confirmed appropriate empirical FWER for all G-functions. A G-function with optimal r2 = 0.3 between index and secondary SNP which down-weighted SNPs at higher distance step-wise-strong and gave more emphasis on d than on r2 had overall best power. It also gave the best results in application to the real datasets. As a proof of concept, 'quasi-adaptive' method was applied to GWAS on free thyroxine (FT4), inflammatory bowel disease (IBD) and human height. Application of the algorithm revealed 5 secondary signals in our example GWAS on FT4, 5 secondary signals in case of the IBD and 19 secondary signals on human height, that would have gone undetected with the established genome-wide threshold (α=5×10-8). AVAILABILITY AND IMPLEMENTATION: https://github.com/sghasemi64/Secondary-Signal. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Obesity and Kidney Function: A Two-Sample Mendelian Randomization Study.

BACKGROUND: Obesity and type 2 diabetes (T2D) are correlated risk factors for chronic kidney disease (CKD). METHODS: Using summary data from GIANT (Genetic Investigation of Anthropometric Traits), DIAGRAM (DIAbetes Genetics Replication And Meta-analysis), and CKDGen (CKD Genetics), we examined causality and directionality of the association between obesity and kidney function. Bidirectional 2-sample Mendelian randomization (MR) estimated the total causal effects of body mass index (BMI) and waist-to-hip ratio (WHR) on kidney function, and vice versa. Effects of adverse obesity and T2D were examined by stratifying BMI variants by their association with WHR and T2D. Multivariable MR estimated the direct causal effects of BMI and WHR on kidney function. The inverse variance weighted random-effects MR for Europeans was the main analysis, accompanied by several sensitivity MR analyses. RESULTS: One standard deviation (SD ≈ 4.8 kg/m2) genetically higher BMI was associated with decreased estimated glomerular filtration rate (eGFR) [ß=-0.032 (95% confidence intervals: -0.036, -0.027) log[eGFR], P = 1 × 10-43], increased blood urea nitrogen (BUN) [ß = 0.010 (0.005, 0.015) log[BUN], P = 3 × 10-6], increased urinary albumin-to-creatinine ratio [ß = 0.199 (0.067, 0.332) log[urinary albumin-to-creatinine ratio (UACR)], P = 0.003] in individuals with diabetes, and increased risk of microalbuminuria [odds ratios (OR) = 1.15 [1.04-1.28], P = 0.009] and CKD [1.13 (1.07-1.19), P = 3 × 10-6]. Corresponding estimates for WHR and for trans-ethnic populations were overall similar. The associations were driven by adverse obesity, and for microalbuminuria additionally by T2D. While genetically high BMI, unlike WHR, was directly associated with eGFR, BUN, and CKD, the pathway to albuminuria was likely through T2D. Genetically predicted kidney function was not associated with BMI or WHR. CONCLUSIONS: Genetically high BMI is associated with impaired kidney function, driven by adverse obesity, and for albuminuria additionally by T2D.

Plasma Proteomics of Renal Function: A Transethnic Meta-Analysis and Mendelian Randomization Study.

BACKGROUND: Studies on the relationship between renal function and the human plasma proteome have identified several potential biomarkers. However, investigations have been conducted largely in European populations, and causality of the associations between plasma proteins and kidney function has never been addressed. METHODS: A cross-sectional study of 993 plasma proteins among 2882 participants in four studies of European and admixed ancestries (KORA, INTERVAL, HUNT, QMDiab) identified transethnic associations between eGFR/CKD and proteomic biomarkers. For the replicated associations, two-sample bidirectional Mendelian randomization (MR) was used to investigate potential causal relationships. Publicly available datasets and transcriptomic data from independent studies were used to examine the association between gene expression in kidney tissue and eGFR. RESULTS: In total, 57 plasma proteins were associated with eGFR, including one novel protein. Of these, 23 were additionally associated with CKD. The strongest inferred causal effect was the positive effect of eGFR on testican-2, in line with the known biological role of this protein and the expression of its protein-coding gene (SPOCK2) in renal tissue. We also observed suggestive evidence of an effect of melanoma inhibitory activity (MIA), carbonic anhydrase III, and cystatin-M on eGFR. CONCLUSIONS: In a discovery-replication setting, we identified 57 proteins transethnically associated with eGFR. The revealed causal relationships are an important stepping stone in establishing testican-2 as a clinically relevant physiological marker of kidney disease progression, and point to additional proteins warranting further investigation.

TREML2 Gene Expression and Its Missense Variant rs3747742 Associate with White Matter Hyperintensity Volume and Alzheimer's Disease-Related Brain Atrophy in the General Population.

Although the common pathology of Alzheimer's disease (AD) and white matter hyperintensities (WMH) is disputed, the gene TREML2 has been implicated in both conditions: its whole-blood gene expression was associated with WMH volume and its missense variant rs3747742 with AD risk. We re-examined those associations within one comprehensive dataset of the general population, additionally searched for cross-relations and illuminated the role of the apolipoprotein E (APOE) ε4 status in the associations. For our linear regression and linear mixed effect models, we used 1949 participants from the Study of Health in Pomerania (Germany). AD was assessed using a continuous pre-symptomatic MRI-based score evaluating a participant's AD-related brain atrophy. In our study, increased whole-blood TREML2 gene expression was significantly associated with reduced WMH volume but not with the AD score. Conversely, rs3747742-C was significantly associated with a reduced AD score but not with WMH volume. The APOE status did not influence the associations. In sum, TREML2 robustly associated with WMH volume and AD-related brain atrophy on different molecular levels. Our results thus underpin TREML2's role in neurodegeneration, might point to its involvement in AD and WMH via different biological mechanisms, and highlight TREML2 as a worthwhile target for disentangling the two pathologies.

White Blood Cells and Blood Pressure: A Mendelian Randomization Study.

BACKGROUND: High blood pressure (BP) is a risk factor for cardiovascular morbidity and mortality. While BP is regulated by the function of kidney, vasculature, and sympathetic nervous system, recent experimental data suggest that immune cells may play a role in hypertension. METHODS: We studied the relationship between major white blood cell types and blood pressure in the UK Biobank population and used Mendelian randomization (MR) analyses using the ≈750 000 UK-Biobank/International Consortium of Blood Pressure-Genome-Wide Association Studies to examine which leukocyte populations may be causally linked to BP. RESULTS: A positive association between quintiles of lymphocyte, monocyte, and neutrophil counts, and increased systolic BP, diastolic BP, and pulse pressure was observed (eg, adjusted systolic BP mean±SE for 1st versus 5th quintile respectively: 140.13±0.08 versus 141.62±0.07 mm Hg for lymphocyte, 139.51±0.08 versus 141.84±0.07 mm Hg for monocyte, and 137.96±0.08 versus 142.71±0.07 mm Hg for neutrophil counts; all P<10-50). Using 121 single nucleotide polymorphisms in MR, implemented through the inverse-variance weighted approach, we identified a potential causal relationship of lymphocyte count with systolic BP and diastolic BP (causal estimates: 0.69 [95% CI, 0.19-1.20] and 0.56 [95% CI, 0.23-0.90] of mm Hg per 1 SD genetically elevated lymphocyte count, respectively), which was directionally concordant to the observational findings. These inverse-variance weighted estimates were consistent with other robust MR methods. The exclusion of rs3184504 SNP in the SH2B3 locus attenuated the magnitude of the signal in some of the MR analyses. MR in the reverse direction found evidence of positive effects of BP indices on counts of monocytes, neutrophils, and eosinophils but not lymphocytes or basophils. Subsequent MR testing of lymphocyte count in the context of genetic correlation with renal function or resting and postexercise heart rate demonstrated a positive association of lymphocyte count with urine albumin-to-creatinine ratio. CONCLUSIONS: Observational and genetic analyses demonstrate a concordant, positive and potentially causal relationship of lymphocyte count with systolic BP and diastolic BP.

Meta-analysis uncovers genome-wide significant variants for rapid kidney function decline.

Rapid decline of glomerular filtration rate estimated from creatinine (eGFRcrea) is associated with severe clinical endpoints. In contrast to cross-sectionally assessed eGFRcrea, the genetic basis for rapid eGFRcrea decline is largely unknown. To help define this, we meta-analyzed 42 genome-wide association studies from the Chronic Kidney Diseases Genetics Consortium and United Kingdom Biobank to identify genetic loci for rapid eGFRcrea decline. Two definitions of eGFRcrea decline were used: 3 mL/min/1.73m2/year or more ("Rapid3"; encompassing 34,874 cases, 107,090 controls) and eGFRcrea decline 25% or more and eGFRcrea under 60 mL/min/1.73m2 at follow-up among those with eGFRcrea 60 mL/min/1.73m2 or more at baseline ("CKDi25"; encompassing 19,901 cases, 175,244 controls). Seven independent variants were identified across six loci for Rapid3 and/or CKDi25: consisting of five variants at four loci with genome-wide significance (near UMOD-PDILT (2), PRKAG2, WDR72, OR2S2) and two variants among 265 known eGFRcrea variants (near GATM, LARP4B). All these loci were novel for Rapid3 and/or CKDi25 and our bioinformatic follow-up prioritized variants and genes underneath these loci. The OR2S2 locus is novel for any eGFRcrea trait including interesting candidates. For the five genome-wide significant lead variants, we found supporting effects for annual change in blood urea nitrogen or cystatin-based eGFR, but not for GATM or LARP4B. Individuals at high compared to those at low genetic risk (8-14 vs. 0-5 adverse alleles) had a 1.20-fold increased risk of acute kidney injury (95% confidence interval 1.08-1.33). Thus, our identified loci for rapid kidney function decline may help prioritize therapeutic targets and identify mechanisms and individuals at risk for sustained deterioration of kidney function.

Liver Phenotypes of European Adults Heterozygous or Homozygous for Pi∗Z Variant of AAT (Pi∗MZ vs Pi∗ZZ genotype) and Noncarriers.

BACKGROUND & AIMS: Homozygosity for the Pi∗Z variant of the gene that encodes the alpha-1 antitrypsin peptide (AAT), called the Pi∗ZZ genotype, causes a liver and lung disease called alpha-1 antitrypsin deficiency. Heterozygosity (the Pi∗MZ genotype) is a risk factor for cirrhosis in individuals with liver disease. Up to 4% of Europeans have the Pi∗MZ genotype; we compared features of adults with and without Pi∗MZ genotype among persons without preexisting liver disease. METHODS: We analyzed data from the European Alpha-1 Liver Cohort, from 419 adults with the Pi∗MZ genotype, 309 adults with the Pi∗ZZ genotype, and 284 individuals without the variant (noncarriers). All underwent a comprehensive evaluation; liver stiffness measurements (LSMs) were made by transient elastography. Liver biopsies were analyzed to define histologic and biochemical features associated with the Pi∗Z variant. Levels of serum transaminases were retrieved from 444,642 participants, available in the United Kingdom biobank. RESULTS: In the UK biobank database, levels of serum transaminases were increased in subjects with the Pi∗MZ genotype compared with noncarriers. In the Alpha-1 Liver Cohort, adults with Pi∗MZ had lower levels of gamma-glutamyl transferase in serum and lower LSMs than adults with the Pi∗ZZ variant, but these were higher than in noncarriers. Ten percent of subjects with the Pi∗MZ genotype vs 4% of noncarriers had LSMs of 7.1 kPa or more (adjusted odds ratio, 4.8; 95% confidence interval, 2.0-11.8). Obesity and diabetes were the most important factors associated with LSMs ≥7.1 kPa in subjects with the Pi∗MZ genotype. AAT inclusions were detected in liver biopsies of 63% of subjects with the Pi∗MZ genotype, vs 97% of subjects with the Pi∗ZZ genotype, and increased with liver fibrosis stages. Subjects with the Pi∗MZ genotype did not have increased hepatic levels of AAT, whereas levels of insoluble AAT varied among individuals. CONCLUSIONS: Adults with the Pi∗MZ genotype have lower levels of serum transaminases, fewer AAT inclusions in liver, and lower liver stiffness than adults with the Pi∗ZZ genotype, but higher than adults without the Pi∗Z variant. These findings should help determine risk of subjects with the Pi∗MZ genotype and aid in counseling.

Mendelian randomization provides evidence for a causal effect of higher serum IGF-1 concentration on risk of hip and knee osteoarthritis.

OBJECTIVES: How insulin-like growth factor-1 (IGF-1) is related to OA is not well understood. We determined relationships between IGF-1 and hospital-diagnosed hand, hip and knee OA in UK Biobank, using Mendelian randomization (MR) to determine causality. METHODS: Serum IGF-1 was assessed by chemiluminescent immunoassay. OA was determined using Hospital Episode Statistics. One-sample MR (1SMR) was performed using two-stage least-squares regression, with an unweighted IGF-1 genetic risk score as an instrument. Multivariable MR included BMI as an additional exposure (instrumented by BMI genetic risk score). MR analyses were adjusted for sex, genotyping chip and principal components. We then performed two-sample MR (2SMR) using summary statistics from Cohorts for Heart and Aging Research in Genetic Epidemiology (CHARGE) (IGF-1, N = 30 884) and the recent genome-wide association study meta-analysis (N = 455 221) of UK Biobank and Arthritis Research UK OA Genetics (arcOGEN). RESULTS: A total of 332 092 adults in UK Biobank had complete data. Their mean (s.d.) age was 56.5 (8.0) years and 54% were female. IGF-1 was observationally related to a reduced odds of hand OA [odds ratio per doubling = 0.87 (95% CI 0.82, 0.93)], and an increased odds of hip OA [1.04 (1.01, 1.07)], but was unrelated to knee OA [0.99 (0.96, 1.01)]. Using 1SMR, we found strong evidence for an increased risk of hip [odds ratio per s.d. increase = 1.57 (1.21, 2.01)] and knee [1.30 (1.07, 1.58)] OA with increasing IGF-1 concentration. By contrast, we found no evidence for a causal effect of IGF-1 concentration on hand OA [0.98 (0.57, 1.70)]. Results were consistent when estimated using 2SMR and in multivariable MR analyses accounting for BMI. CONCLUSION: We have found evidence that increased serum IGF-1 is causally related to higher risk of hip and knee OA.

Thyroid function, sex hormones and sexual function: a Mendelian randomization study.

Hypothyroidism and hyperthyroidism are observationally associated with sex hormone concentrations and sexual dysfunction, but causality is unclear. We investigated whether TSH, fT4, hypo- and hyperthyroidism are causally associated with sex hormones and sexual function. We used publicly available summary statistics from genome-wide association studies on TSH and fT4 and hypo- and hyperthyroidism from the ThyroidOmics Consortium (N ≤ 54,288). Outcomes from UK Biobank (women ≤ 194,174/men ≤ 167,020) and ReproGen (women ≤ 252,514) were sex hormones (sex hormone binding globulin [SHBG], testosterone, estradiol, free androgen index [FAI]) and sexual function (ovulatory function in women: duration of menstrual period, age at menarche and menopause, reproductive lifespan, and erectile dysfunction in men). We performed two-sample Mendelian randomization (MR) analyses on summary level, and unweighted genetic risk score (GRS) analysis on individual level data. One SD increase in TSH was associated with a 1.332 nmol/L lower (95% CI: - 0.717,- 1.946; p = 2 × 10-5) SHBG and a 0.103 nmol/l lower (- 0.051,V0.154; p = 9 × 10-5) testosterone in two-sample MR, supported by the GRS approach. Genetic predisposition to hypothyroidism was associated with decreased and genetic predisposition to hyperthyroidism with increased SHBG and testosterone in both approaches. The GRS for fT4 was associated with increased testosterone and estradiol in women only. The GRS for TSH and hypothyroidism were associated with increased and the GRS for hyperthyroidism with decreased FAI in men only. While genetically predicted thyroid function was associated with sex hormones, we found no association with sexual function.

Meta-analysis of epigenome-wide association studies of carotid intima-media thickness.

Common carotid intima-media thickness (cIMT) is an index of subclinical atherosclerosis that is associated with ischemic stroke and coronary artery disease (CAD). We undertook a cross-sectional epigenome-wide association study (EWAS) of measures of cIMT in 6400 individuals. Mendelian randomization analysis was applied to investigate the potential causal role of DNA methylation in the link between atherosclerotic cardiovascular risk factors and cIMT or clinical cardiovascular disease. The CpG site cg05575921 was associated with cIMT (beta = -0.0264, p value = 3.5 × 10-8) in the discovery panel and was replicated in replication panel (beta = -0.07, p value = 0.005). This CpG is located at chr5:81649347 in the intron 3 of the aryl hydrocarbon receptor repressor gene (AHRR). Our results indicate that DNA methylation at cg05575921 might be in the pathway between smoking, cIMT and stroke. Moreover, in a region-based analysis, 34 differentially methylated regions (DMRs) were identified of which a DMR upstream of ALOX12 showed the strongest association with cIMT (p value = 1.4 × 10-13). In conclusion, our study suggests that DNA methylation may play a role in the link between cardiovascular risk factors, cIMT and clinical cardiovascular disease.