Effects of genetically predicted posttraumatic stress disorder on autoimmune phenotypes

Observational studies suggest that posttraumatic stress disorder (PTSD) increases risk for various autoimmune diseases. Insights into shared biology and causal relationships between these diseases may inform intervention approaches to PTSD and co-morbid autoimmune conditions. We investigated the shared genetic contributions and causal relationships between PTSD, 18 autoimmune diseases, and 3 immune/inflammatory biomarkers. Univariate MiXeR was used to contrast the genetic architectures of phenotypes. Genetic correlations were estimated using linkage disequilibrium score regression. Bi-directional, two-sample Mendelian randomization (MR) was performed using independent, genome-wide significant single nucleotide polymorphisms; inverse variance weighted and weighted median MR estimates were evaluated. Sensitivity analyses for uncorrelated (MR PRESSO) and correlated horizontal pleiotropy (CAUSE) were also performed. PTSD was considerably more polygenic (10,863 influential variants) than autoimmune diseases (median 255 influential variants). However, PTSD evidenced significant genetic correlation with nine autoimmune diseases and three inflammatory biomarkers. PTSD had putative causal effects on autoimmune thyroid disease (p = 0.00009) and C-reactive protein (CRP) (p = 4.3 × 10-7). Inferences were not substantially altered by sensitivity analyses. Additionally, the PTSD-autoimmune thyroid disease association remained significant in multivariable MR analysis adjusted for genetically predicted inflammatory biomarkers as potential mechanistic pathway variables. No autoimmune disease had a significant causal effect on PTSD (all p values > 0.05). Although causal effect models were supported for associations of PTSD with CRP, shared pleiotropy was adequate to explain a putative causal effect of CRP on PTSD (p = 0.18). In summary, our results suggest a significant genetic overlap between PTSD, autoimmune diseases, and biomarkers of inflammation. PTSD has a putative causal effect on autoimmune thyroid disease, consistent with existing epidemiologic evidence. A previously reported causal effect of CRP on PTSD is potentially confounded by shared genetics. Together, results highlight the nuanced links between PTSD, autoimmune disorders, and associated inflammatory signatures, and suggest the importance of targeting related pathways to protect against disease and disability.

Contributions of common genetic variants to constitutional delay of puberty and idiopathic hypogonadotropic hypogonadism

CONTEXT

Constitutional delay of puberty (CDP) is highly heritable, but the genetic basis for CDP is largely unknown. Idiopathic hypogonadotropic hypogonadism (IHH) can be caused by rare genetic variants, but in about half of cases, no rare-variant cause is found.

OBJECTIVE

To determine whether common genetic variants that influence pubertal timing contribute to CDP and IHH.

DESIGN

Case-control study.

PARTICIPANTS

80 individuals with CDP; 301 with normosmic IHH, and 348 with Kallmann syndrome; control genotyping data from unrelated studies.

MAIN OUTCOME MEASURES

Polygenic scores (PGS) based on genome-wide association studies for timing of male pubertal hallmarks and age at menarche (AAM).

RESULTS

The CDP cohort had higher PGS for male pubertal hallmarks and for AAM compared to controls (for male hallmarks, Cohen's d = 0.85, p = 1 × 10-16; for AAM, d = 0.67, p = 1 × 10-10). The normosmic IHH cohort also had higher PGS for male hallmarks compared to controls, but the difference was smaller (male hallmarks d = 0.20, p = 0.003; AAM d = 0.10, p = 0.055). No differences were seen for the KS cohort compared to controls (male hallmarks d = 0.04, p = 0.45; AAM d = -0.03, p = 0.86).

CONCLUSIONS

Common genetic variants that influence pubertal timing in the general population contribute strongly to the genetics of CDP, weakly to normosmic IHH, and potentially not at all to KS. These findings demonstrate that the common-variant genetics of CDP and normosmic IHH are largely but not entirely distinct.

Functional EPAS1/HIF2A missense variant is associated with hematocrit in Andean highlanders

Hypoxia-inducible factor pathway genes are linked to adaptation in both human and nonhuman highland species. EPAS1, a notable target of hypoxia adaptation, is associated with relatively lower hemoglobin concentration in Tibetans. We provide evidence for an association between an adaptive EPAS1 variant (rs570553380) and the same phenotype of relatively low hematocrit in Andean highlanders. This Andean-specific missense variant is present at a modest frequency in Andeans and absent in other human populations and vertebrate species except the coelacanth. CRISPR-base-edited human cells with this variant exhibit shifts in hypoxia-regulated gene expression, while metabolomic analyses reveal both genotype and phenotype associations and validation in a lowland population. Although this genocopy of relatively lower hematocrit in Andean highlanders parallels well-replicated findings in Tibetans, it likely involves distinct pathway responses based on a protein-coding versus noncoding variants, respectively. These findings illuminate how unique variants at EPAS1 contribute to the same phenotype in Tibetans and a subset of Andean highlanders despite distinct evolutionary trajectories.

Pediatric Features of Genetic Predisposition to Polycystic Ovary Syndrome

CONTEXT

Polycystic ovary syndrome (PCOS) has historically been conceptualized as a disorder of the reproductive system in women. However, offspring of women with PCOS begin to show metabolic features of PCOS in childhood, suggestive of childhood manifestations.

OBJECTIVE

To identify childhood manifestations of genetic risk for PCOS.

METHODS

We calculated a PCOS polygenic risk score (PRS) for 12 350 girls and boys in 4 pediatric cohorts-ALSPAC (UK), COPSAC (Denmark), Project Viva (USA), and The HOLBÆK Study (Denmark). We tested for association of the PRS with PCOS-related phenotypes throughout childhood and with age at pubarche and age at peak height velocity and meta-analyzed effects across cohorts using fixed-effect models.

RESULTS

Higher PRS for PCOS was associated with higher body mass index in midchildhood (0.05 kg/m2 increase per 1 SD of PRS, 95% CI 0.03, 0.07, P = 3 × 10-5) and higher risk of obesity in early childhood (OR 1.34, 95% CI 1.13, 1.59, P = .0009); both persisted through late adolescence (P all ≤.03). Higher PCOS PRS was associated with earlier age at pubarche (0.85-month decrease per 1 SD of PRS, 95% CI -1.44, -0.26, P = .005) and younger age at peak height velocity (0.64-month decrease per 1 SD of PRS, 95% CI -0.94, -0.33, P = 4 × 10-5).

CONCLUSION

Genetic risk factors for PCOS are associated with alterations in metabolic, growth, and developmental traits in childhood. Thus, PCOS may not simply be a condition that affects women of reproductive age but, rather, a possible manifestation of an underlying condition that affects both sexes starting in early life.

Autoimmune alleles at the major histocompatibility locus modify melanoma susceptibility

Autoimmunity and cancer represent two different aspects of immune dysfunction. Autoimmunity is characterized by breakdowns in immune self-tolerance, while impaired immune surveillance can allow for tumorigenesis. The class I major histocompatibility complex (MHC-I), which displays derivatives of the cellular peptidome for immune surveillance by CD8+ T cells, serves as a common genetic link between these conditions. As melanoma-specific CD8+ T cells have been shown to target melanocyte-specific peptide antigens more often than melanoma-specific antigens, we investigated whether vitiligo- and psoriasis-predisposing MHC-I alleles conferred a melanoma-protective effect. In individuals with cutaneous melanoma from both The Cancer Genome Atlas (n = 451) and an independent validation set (n = 586), MHC-I autoimmune-allele carrier status was significantly associated with a later age of melanoma diagnosis. Furthermore, MHC-I autoimmune-allele carriers were significantly associated with decreased risk of developing melanoma in the Million Veteran Program (OR = 0.962, p = 0.024). Existing melanoma polygenic risk scores (PRSs) did not predict autoimmune-allele carrier status, suggesting these alleles provide orthogonal risk-relevant information. Mechanisms of autoimmune protection were neither associated with improved melanoma-driver mutation association nor improved gene-level conserved antigen presentation relative to common alleles. However, autoimmune alleles showed higher affinity relative to common alleles for particular windows of melanocyte-conserved antigens and loss of heterozygosity of autoimmune alleles caused the greatest reduction in presentation for several conserved antigens across individuals with loss of HLA alleles. Overall, this study presents evidence that MHC-I autoimmune-risk alleles modulate melanoma risk unaccounted for by current PRSs.

Genetics of skeletal proportions in two different populations

Human height can be divided into sitting height and leg length, reflecting growth of different parts of the skeleton whose relative proportions are captured by the ratio of sitting to total height (as sitting height ratio, SHR). Height is a highly heritable trait, and its genetic basis has been well-studied. However, the genetic determinants of skeletal proportion are much less well-characterized. Expanding substantially on past work, we performed a genome-wide association study (GWAS) of SHR in ∼450,000 individuals with European ancestry and ∼100,000 individuals with East Asian ancestry from the UK and China Kadoorie Biobanks. We identified 565 loci independently associated with SHR, including all genomic regions implicated in prior GWAS in these ancestries. While SHR loci largely overlap height-associated loci (P < 0.001), the fine-mapped SHR signals were often distinct from height. We additionally used fine-mapped signals to identify 36 credible sets with heterogeneous effects across ancestries. Lastly, we used SHR, sitting height, and leg length to identify genetic variation acting on specific body regions rather than on overall human height.

Germline modifiers of the tumor immune microenvironment implicate drivers of cancer risk and immunotherapy response

With the continued promise of immunotherapy for treating cancer, understanding how host genetics contributes to the tumor immune microenvironment (TIME) is essential to tailoring cancer screening and treatment strategies. Here, we study 1084 eQTLs affecting the TIME found through analysis of The Cancer Genome Atlas and literature curation. These TIME eQTLs are enriched in areas of active transcription, and associate with gene expression in specific immune cell subsets, such as macrophages and dendritic cells. Polygenic score models built with TIME eQTLs reproducibly stratify cancer risk, survival and immune checkpoint blockade (ICB) response across independent cohorts. To assess whether an eQTL-informed approach could reveal potential cancer immunotherapy targets, we inhibit CTSS, a gene implicated by cancer risk and ICB response-associated polygenic models; CTSS inhibition results in slowed tumor growth and extended survival in vivo. These results validate the potential of integrating germline variation and TIME characteristics for uncovering potential targets for immunotherapy.

Differential Methylation of Telomere-Related Genes Is Associated with Kidney Disease in Individuals with Type 1 Diabetes

Diabetic kidney disease (DKD) represents a major global health problem. Accelerated ageing is a key feature of DKD and, therefore, characteristics of accelerated ageing may provide useful biomarkers or therapeutic targets. Harnessing multi-omics, features affecting telomere biology and any associated methylome dysregulation in DKD were explored. Genotype data for nuclear genome polymorphisms in telomere-related genes were extracted from genome-wide case-control association data (n = 823 DKD/903 controls; n = 247 end-stage kidney disease (ESKD)/1479 controls). Telomere length was established using quantitative polymerase chain reaction. Quantitative methylation values for 1091 CpG sites in telomere-related genes were extracted from epigenome-wide case-control association data (n = 150 DKD/100 controls). Telomere length was significantly shorter in older age groups (p = 7.6 × 10^-6). Telomere length was also significantly reduced (p = 6.6 × 10^-5) in DKD versus control individuals, with significance remaining after covariate adjustment (p = 0.028). DKD and ESKD were nominally associated with telomere-related genetic variation, with Mendelian randomisation highlighting no significant association between genetically predicted telomere length and kidney disease. A total of 496 CpG sites in 212 genes reached epigenome-wide significance (p ≤ 10^-8) for DKD association, and 412 CpG sites in 193 genes for ESKD. Functional prediction revealed differentially methylated genes were enriched for Wnt signalling involvement. Harnessing previously published RNA-sequencing datasets, potential targets where epigenetic dysregulation may result in altered gene expression were revealed, useful as potential diagnostic and therapeutic targets for intervention.

Associations of Accelerometer-Measured Physical Activity and Sedentary Time With All-Cause Mortality by Genetic Predisposition for Longevity

The goal of this study was to examine associations between accelerometer-measured physical activity (PA) and sedentary time (ST) with mortality by a genetic risk score (GRS) for longevity. Among 5,446 women, (mean [SD]: age, 78.2 [6.6] years), 1,022 deaths were observed during 33,350 person-years of follow-up. Using multivariable Cox proportional hazards models, higher light PA and moderate to vigorous PA were associated with lower mortality across all GRS for longevity categories (low/medium/high; all ptrend < .001). Higher ST was associated with higher mortality (ptrend across all GRS categories < .001). Interaction tests for PA and ST with the GRS were not statistically significant. Findings support the importance of higher PA and lower ST for reducing mortality risk in older women, regardless of genetic predisposition for longevity.

Epigenome-wide meta-analysis identifies DNA methylation biomarkers associated with diabetic kidney disease

Type 1 diabetes affects over nine million individuals globally, with approximately 40% developing diabetic kidney disease. Emerging evidence suggests that epigenetic alterations, such as DNA methylation, are involved in diabetic kidney disease. Here we assess differences in blood-derived genome-wide DNA methylation associated with diabetic kidney disease in 1304 carefully characterised individuals with type 1 diabetes and known renal status from two cohorts in the United Kingdom-Republic of Ireland and Finland. In the meta-analysis, we identify 32 differentially methylated CpGs in diabetic kidney disease in type 1 diabetes, 18 of which are located within genes differentially expressed in kidneys or correlated with pathological traits in diabetic kidney disease. We show that methylation at 21 of the 32 CpGs predict the development of kidney failure, extending the knowledge and potentially identifying individuals at greater risk for diabetic kidney disease in type 1 diabetes.

Rare copy number variation in posttraumatic stress disorder

Posttraumatic stress disorder (PTSD) is a heritable (h2 = 24-71%) psychiatric illness. Copy number variation (CNV) is a form of rare genetic variation that has been implicated in the etiology of psychiatric disorders, but no large-scale investigation of CNV in PTSD has been performed. We present an association study of CNV burden and PTSD symptoms in a sample of 114,383 participants (13,036 cases and 101,347 controls) of European ancestry. CNVs were called using two calling algorithms and intersected to a consensus set. Quality control was performed to remove strong outlier samples. CNVs were examined for association with PTSD within each cohort using linear or logistic regression analysis adjusted for population structure and CNV quality metrics, then inverse variance weighted meta-analyzed across cohorts. We examined the genome-wide total span of CNVs, enrichment of CNVs within specified gene-sets, and CNVs overlapping individual genes and implicated neurodevelopmental regions. The total distance covered by deletions crossing over known neurodevelopmental CNV regions was significant (beta = 0.029, SE = 0.005, P = 6.3 × 10-8). The genome-wide neurodevelopmental CNV burden identified explains 0.034% of the variation in PTSD symptoms. The 15q11.2 BP1-BP2 microdeletion region was significantly associated with PTSD (beta = 0.0206, SE = 0.0056, P = 0.0002). No individual significant genes interrupted by CNV were identified. 22 gene pathways related to the function of the nervous system and brain were significant in pathway analysis (FDR q < 0.05), but these associations were not significant once NDD regions were removed. A larger sample size, better detection methods, and annotated resources of CNV are needed to explore this relationship further.

Genome-wide meta-analysis and omics integration identifies novel genes associated with diabetic kidney disease

AIMS/HYPOTHESIS

Diabetic kidney disease (DKD) is the leading cause of kidney failure and has a substantial genetic component. Our aim was to identify novel genetic factors and genes contributing to DKD by performing meta-analysis of previous genome-wide association studies (GWAS) on DKD and by integrating the results with renal transcriptomics datasets.

METHODS

We performed GWAS meta-analyses using ten phenotypic definitions of DKD, including nearly 27,000 individuals with diabetes. Meta-analysis results were integrated with estimated quantitative trait locus data from human glomerular (N=119) and tubular (N=121) samples to perform transcriptome-wide association study. We also performed gene aggregate tests to jointly test all available common genetic markers within a gene, and combined the results with various kidney omics datasets.

RESULTS

The meta-analysis identified a novel intronic variant (rs72831309) in the TENM2 gene associated with a lower risk of the combined chronic kidney disease (eGFR<60 ml/min per 1.73 m2) and DKD (microalbuminuria or worse) phenotype (p=9.8×10-9; although not withstanding correction for multiple testing, p>9.3×10-9). Gene-level analysis identified ten genes associated with DKD (COL20A1, DCLK1, EIF4E, PTPRN-RESP18, GPR158, INIP-SNX30, LSM14A and MFF; p<2.7×10-6). Integration of GWAS with human glomerular and tubular expression data demonstrated higher tubular AKIRIN2 gene expression in individuals with vs without DKD (p=1.1×10-6). The lead SNPs within six loci significantly altered DNA methylation of a nearby CpG site in kidneys (p<1.5×10-11). Expression of lead genes in kidney tubules or glomeruli correlated with relevant pathological phenotypes (e.g. TENM2 expression correlated positively with eGFR [p=1.6×10-8] and negatively with tubulointerstitial fibrosis [p=2.0×10-9], tubular DCLK1 expression correlated positively with fibrosis [p=7.4×10-16], and SNX30 expression correlated positively with eGFR [p=5.8×10-14] and negatively with fibrosis [p<2.0×10-16]).

CONCLUSIONS/INTERPRETATION

Altogether, the results point to novel genes contributing to the pathogenesis of DKD.

DATA AVAILABILITY

The GWAS meta-analysis results can be accessed via the type 1 and type 2 diabetes (T1D and T2D, respectively) and Common Metabolic Diseases (CMD) Knowledge Portals, and downloaded on their respective download pages ( https://t1d.hugeamp.org/downloads.html ; https://t2d.hugeamp.org/downloads.html ; https://hugeamp.org/downloads.html ).

Abstract 3825: Germline modifiers of the tumor immune microenvironment reveal drivers of immunotherapy response

With the continued promise of immunotherapy as an avenue for treating cancer, understanding how host genetics contributes to the tumor immune microenvironment (TIME) is essential to tailoring cancer risk screening and treatment strategies. Using genotypes from over 8,000 European individuals in The Cancer Genome Atlas and 137 heritable tumor immune phenotype components (IP components), we identified and investigated 532 TIME-SNPs. Focusing on 77 variants that were relevant to cancer risk, survival, or treatment response, we explored their potential to reveal novel targets for immunotherapy. Many variants overlapped regions with histone marks indicating active transcription, and influenced gene activities in specific immune cell subsets, such as macrophages and dendritic cells. TIME-SNPs implicated genes such as LAIR1, TREX1, CTSS, CTSW and LILRB2 were differentially expressed between responders and non-responders to immune-checkpoint blockade (ICB) in preclinical studies. Of these, LILRB2 and LAIR1 have already been identified as putative targets for immunotherapy. Here we found that inhibition of CTSS led to better tumor control and survival in murine models, alone or in combination with anti-PD-1. Collectively we show that through an integrative approach, it is possible to link host genetics to TIME characteristics, informing novel biomarkers for cancer risk and target identification in immunotherapy.

Citation Format: Meghana Pagadala, Victoria Wu, Eva Pérez-Guijarro, Hyo Kim, Andrea Castro, James Talwar, Cristian Gonzalez-Colin, Steven Cao, Benjamin J. Schmiedel, Timothy Sears, Shervin Goudarzi, Divya Kirani, Rany M. Salem, Gerald P. Morris, Olivier Harismendy, Sandip P. Patel, Jill P. Mesirov, Maurizio Zanetti, Chi-Ping Day, Chun C. Fan, Wesley K. Thompson, Glenn Merlino, J. Silvio Gutkind, Pandurangan Vijayanand, Hannah Carter. Germline modifiers of the tumor immune microenvironment reveal drivers of immunotherapy response [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 3825.

Markers of kidney function, genetic variation related to cognitive function, and cognitive performance in the UK Biobank

BACKGROUND

Chronic kidney disease has been linked to worse cognition. However, this association may be dependent on the marker of kidney function used, and studies assessing modification by genetics are lacking. This study examined associations between multiple measures of kidney function and assessed effect modification by a polygenic score for general cognitive function.

METHODS

In this cross-sectional study of up to 341,208 European ancestry participants from the UK Biobank study, we examined associations between albuminuria and estimated glomerular filtration rate based on creatinine (eGFRcre) or cystatin C (eGFRcys) with cognitive performance on tests of verbal-numeric reasoning, reaction time and visual memory. Adjustment for confounding factors was performed using multivariate regression and propensity-score matching. Interaction between kidney function markers and a polygenic risk score for general cognitive function was also assessed.

RESULTS

Albuminuria was associated with worse performance on tasks of verbal-numeric reasoning (β(points) = -0.09, p < 0.001), reaction time (β(milliseconds) = 7.06, p < 0.001) and visual memory (β(log errors) = 0.013, p = 0.01). A polygenic score for cognitive function modified the association between albuminuria and verbal-numeric reasoning with significantly lower scores in those with albuminuria and a lower polygenic score (p = 0.009). Compared to participants with eGFRcre ≥ 60 ml/min, those with eGFRcre < 60 ml/min had lower verbal-numeric reasoning scores and slower mean reaction times (verbal numeric reasoning β = -0.11, p < 0.001 and reaction time β = 6.08, p < 0.001 for eGFRcre < 60 vs eGFRcre ≥ 60). Associations were stronger using cystatin C-based eGFR than creatinine-based eGFR (verbal numeric reasoning β = -0.21, p < 0.001 and reaction time β = 11.21, p < 0.001 for eGFRcys < 60 vs eGFRcys ≥ 60).

CONCLUSIONS

Increased urine albumin is associated with worse cognition, but this may depend on genetic risk. Cystatin C-based eGFR may better predict cognitive performance than creatinine-based estimates.

Enhancing Discovery of Genetic Variants for Posttraumatic Stress Disorder Through Integration of Quantitative Phenotypes and Trauma Exposure Information

BACKGROUND

Posttraumatic stress disorder (PTSD) is heritable and a potential consequence of exposure to traumatic stress. Evidence suggests that a quantitative approach to PTSD phenotype measurement and incorporation of lifetime trauma exposure (LTE) information could enhance the discovery power of PTSD genome-wide association studies (GWASs).

METHODS

A GWAS on PTSD symptoms was performed in 51 cohorts followed by a fixed-effects meta-analysis (N = 182,199 European ancestry participants). A GWAS of LTE burden was performed in the UK Biobank cohort (N = 132,988). Genetic correlations were evaluated with linkage disequilibrium score regression. Multivariate analysis was performed using Multi-Trait Analysis of GWAS. Functional mapping and annotation of leading loci was performed with FUMA. Replication was evaluated using the Million Veteran Program GWAS of PTSD total symptoms.

RESULTS

GWASs of PTSD symptoms and LTE burden identified 5 and 6 independent genome-wide significant loci, respectively. There was a 72% genetic correlation between PTSD and LTE. PTSD and LTE showed largely similar patterns of genetic correlation with other traits, albeit with some distinctions. Adjusting PTSD for LTE reduced PTSD heritability by 31%. Multivariate analysis of PTSD and LTE increased the effective sample size of the PTSD GWAS by 20% and identified 4 additional loci. Four of these 9 PTSD loci were independently replicated in the Million Veteran Program.

CONCLUSIONS

Through using a quantitative trait measure of PTSD, we identified novel risk loci not previously identified using prior case-control analyses. PTSD and LTE have a high genetic overlap that can be leveraged to increase discovery power through multivariate methods.

Causal association of cognitive reserve on Alzheimer's disease with putative sex difference

INTRODUCTION

Sex-dependent risk factors may underlie sex differences in Alzheimer's disease (AD).

METHODS

Using sex-stratified genome-wide association studies (GWAS) of AD, we evaluated associations of 12 traits with AD through polygenic risk scores (PRS) and Mendelian randomization (MR), and explored joint genetic architecture among significant traits by genomic structural equation modeling and network analysis.

RESULTS

AD was associated with lower PRS for premorbid cognitive performance, intelligence, and educational attainment. MR showed a causal role for the cognition-related traits in AD, particularly among females. Their joint genetic components encompassed RNA processing, neuron projection development, and cell cycle pathways that overlap with cellular senescence. Cholesterol and C-reactive protein showed pleiotropy but no causality with AD.

DISCUSSION

Lower cognitive reserve is causally related to AD. The stronger causal link between cognitive performance and AD in females, despite similar PRS between sexes, suggest these differences may result from gene-environmental interactions accumulated over the lifespan.

Biomarkers of kidney function and cognitive ability: A Mendelian randomization study

BACKGROUND

Estimated glomerular filtration rate (eGFR), albuminuria and serum uric acid (SUA) are markers of kidney function that have been associated with cognitive ability. However, whether these associations are causal is unclear.

METHODS

We performed one-sample Mendelian randomization (MR) to estimate the effects of kidney function markers on cognitive performance using data from the UK Biobank. Polygenic scores for SUA, urine albumin to creatinine ratio (ACR), estimated glomerular filtration rate based on serum creatinine (eGFRcre) and serum cystatin C (eGFRcys) were used as instrumental variables, and cognitive function outcomes included tests of verbal-numeric reasoning, reaction time, visual memory, and numeric memory.

RESULTS

We found no evidence of a causal effect of genetically determined SUA, eGFRcre or eGFRcys on cognitive function outcomes. There was no association between a polygenic score for ACR and verbal-numeric reasoning or numeric memory. However, there was suggestive evidence of a relationship between genetically increased ACR and slower reaction time and worse visual memory. ACR was no longer significantly associated with visual memory in analyses using an unweighted polygenic score and in analyses stratified by sex and age category. Pleiotropy adjusted estimates were directionally consistent with those of the principal analysis but overlapped with the null.

CONCLUSIONS

This MR study does not support causal effects of SUA, eGFRcre or eGFRcys on cognitive performance. Genetically increased ACR was associated with slower processing speed and visual memory, but results need confirmation in independent samples.

Similar Genetic Architecture of Alzheimer's Disease and Differential APOE Effect Between Sexes

Sex differences have been observed in the clinical manifestations of Alzheimer’s disease (AD) and elucidating their genetic basis is an active research topic. Based on autosomal genotype data of 7,216 men and 10,680 women, including 8,136 AD cases and 9,760 controls, we explored sex-related genetic heterogeneity in AD by investigating SNP heritability, genetic correlation, as well as SNP- and gene-based genome-wide analyses. We found similar SNP heritability (men: 19.5%; women: 21.5%) and high genetic correlation (R_g = 0.96) between the sexes. The heritability of APOE ε4-related risks for AD, after accounting for effects of all SNPs excluding chromosome 19, was nominally, but not significantly, higher in women (10.6%) than men (9.7%). In age-stratified analyses, ε3/ε4 was associated with a higher risk of AD among women than men aged 65–75 years, but not in the full sample. Apart from APOE, no new significant locus was identified in sex-stratified gene-based analyses. Our result of the high genetic correlation indicates overall similar genetic architecture of AD in both sexes at the genome-wide averaged level. Our study suggests that clinically observed sex differences may arise from sex-specific variants with small effects or more complicated mechanisms involving epigenetic alterations, sex chromosomes, or gene-environment interactions.

Assessment of differentially methylated loci in individuals with end-stage kidney disease attributed to diabetic kidney disease: an exploratory study

BACKGROUND

A subset of individuals with type 1 diabetes mellitus (T1DM) are predisposed to developing diabetic kidney disease (DKD), the most common cause globally of end-stage kidney disease (ESKD). Emerging evidence suggests epigenetic changes in DNA methylation may have a causal role in both T1DM and DKD. The aim of this exploratory investigation was to assess differences in blood-derived DNA methylation patterns between individuals with T1DM-ESKD and individuals with long-duration T1DM but no evidence of kidney disease upon repeated testing to identify potential blood-based biomarkers. Blood-derived DNA from individuals (107 cases, 253 controls and 14 experimental controls) were bisulphite treated before DNA methylation patterns from both groups were generated and analysed using Illumina's Infinium MethylationEPIC BeadChip arrays (n = 862,927 sites). Differentially methylated CpG sites (dmCpGs) were identified (false discovery rate adjusted p ≤ × 10^-8 and fold change ± 2) by comparing methylation levels between ESKD cases and T1DM controls at single site resolution. Gene annotation and functionality was investigated to enrich and rank methylated regions associated with ESKD in T1DM.

RESULTS

Top-ranked genes within which several dmCpGs were located and supported by functional data with methylation look-ups in other cohorts include: AFF3, ARID5B, CUX1, ELMO1, FKBP5, HDAC4, ITGAL, LY9, PIM1, RUNX3, SEPTIN9 and UPF3A. Top-ranked enrichment pathways included pathways in cancer, TGF-β signalling and Th17 cell differentiation.

CONCLUSIONS

Epigenetic alterations provide a dynamic link between an individual's genetic background and their environmental exposures. This robust evaluation of DNA methylation in carefully phenotyped individuals has identified biomarkers associated with ESKD, revealing several genes and implicated key pathways associated with ESKD in individuals with T1DM.

Markers of Kidney Function and Longitudinal Cognitive Ability Among Older Community-Dwelling Adults: The Rancho Bernardo Study

BACKGROUND

Reduced kidney function has been associated with cognitive decline. Most studies have examined a single marker of kidney function and have limited duration of follow-up.

OBJECTIVE

This study evaluated associations between markers of kidney function (urine albumin, estimated glomerular filtration rate [eGFR], and hyperuricemia) with cognitive performance over time.

METHODS

This is a longitudinal study of 1,634 community-dwelling adults (mean age = 71.7 years), with kidney function markers and cognitive ability measured at baseline (1992-1996) and at up to five additional time points with a maximum of 23.4 years (mean = 8.1 years) of follow-up. Associations between kidney function and cognitive performance were assessed using linear mixed effects models. Testing for interaction by sex was conducted.

RESULTS

Albuminuria (urine albumin-to-creatinine ratio [ACR]≥30 mg/g) was associated with steeper annual declines in global cognitive function (MMSE, β= -0.12, p = 0.003), executive function (Trails B, β= 4.50, p < 0.0001) and episodic memory (Buschke total recall, β= -0.62, p = 0.02) scores in men. Results were similar when cognitive test scores were regressed on latent trajectory classes of ACR. In men, hyperuricemia (serum uric acid [SUA]≥6.8 mg/dl for men and SUA≥6.0 mg/dl for women) was associated with lower baseline MMSE (β= -0.70, p = 0.009) scores but not with MMSE change over time. No such associations were detected in women. There were no significant associations between eGFR and cognitive performance for either sex.

CONCLUSION

In older men, albuminuria is an independent predictor of subsequent cognitive decline. More investigations are needed to explain the observed sex differences and the potential relationship between hyperuricemia and poorer global cognition.

Genome-Wide Association Study of Diabetic Kidney Disease Highlights Biology Involved in Glomerular Basement Membrane Collagen

BACKGROUND

Although diabetic kidney disease demonstrates both familial clustering and single nucleotide polymorphism heritability, the specific genetic factors influencing risk remain largely unknown.

METHODS

To identify genetic variants predisposing to diabetic kidney disease, we performed genome-wide association study (GWAS) analyses. Through collaboration with the Diabetes Nephropathy Collaborative Research Initiative, we assembled a large collection of type 1 diabetes cohorts with harmonized diabetic kidney disease phenotypes. We used a spectrum of ten diabetic kidney disease definitions based on albuminuria and renal function.

RESULTS

Our GWAS meta-analysis included association results for up to 19,406 individuals of European descent with type 1 diabetes. We identified 16 genome-wide significant risk loci. The variant with the strongest association (rs55703767) is a common missense mutation in the collagen type IV alpha 3 chain (COL4A3) gene, which encodes a major structural component of the glomerular basement membrane (GBM). Mutations in COL4A3 are implicated in heritable nephropathies, including the progressive inherited nephropathy Alport syndrome. The rs55703767 minor allele (Asp326Tyr) is protective against several definitions of diabetic kidney disease, including albuminuria and ESKD, and demonstrated a significant association with GBM width; protective allele carriers had thinner GBM before any signs of kidney disease, and its effect was dependent on glycemia. Three other loci are in or near genes with known or suggestive involvement in this condition (BMP7) or renal biology (COLEC11 and DDR1).

CONCLUSIONS

The 16 diabetic kidney disease-associated loci may provide novel insights into the pathogenesis of this condition and help identify potential biologic targets for prevention and treatment.

Genetic Evidence That Carbohydrate-Stimulated Insulin Secretion Leads to Obesity

BACKGROUND

A fundamental precept of the carbohydrate-insulin model of obesity is that insulin secretion drives weight gain. However, fasting hyperinsulinemia can also be driven by obesity-induced insulin resistance. We used genetic variation to isolate and estimate the potentially causal effect of insulin secretion on body weight.

METHODS

Genetic instruments of variation of insulin secretion [assessed as insulin concentration 30 min after oral glucose (insulin-30)] were used to estimate the causal relationship between increased insulin secretion and body mass index (BMI), using bidirectional Mendelian randomization analysis of genome-wide association studies. Data sources included summary results from the largest published metaanalyses of predominantly European ancestry for insulin secretion (n = 26037) and BMI (n = 322154), as well as individual-level data from the UK Biobank (n = 138541). Data from the Cardiology and Metabolic Patient Cohort study at Massachusetts General Hospital (n = 1675) were used to validate genetic associations with insulin secretion and to test the observational association of insulin secretion and BMI.

RESULTS

Higher genetically determined insulin-30 was strongly associated with higher BMI (β = 0.098, P = 2.2 × 10^-21), consistent with a causal role in obesity. Similar positive associations were noted in sensitivity analyses using other genetic variants as instrumental variables. By contrast, higher genetically determined BMI was not associated with insulin-30.

CONCLUSIONS

Mendelian randomization analyses provide evidence for a causal relationship of glucose-stimulated insulin secretion on body weight, consistent with the carbohydrate-insulin model of obesity.

Interrogation of human hematopoiesis at single-cell and single-variant resolution

Widespread linkage disequilibrium and incomplete annotation of cell-to-cell state variation represent substantial challenges to elucidating mechanisms of trait-associated genetic variation. Here, we perform genetic fine-mapping for blood cell traits in the UK Biobank to identify putative causal variants. These variants are enriched in genes encoding for proteins in trait-relevant biological pathways and in accessible chromatin of hematopoietic progenitors. For regulatory variants, we explore patterns of developmental enhancer activity, predict molecular mechanisms, and identify likely target genes. In several instances, we localize multiple independent variants to the same regulatory element or gene. We further observe that variants with pleiotropic effects preferentially act in common progenitor populations to direct the production of distinct lineages. Finally, we leverage fine-mapped variants in conjunction with continuous epigenomic annotations to identify trait-cell type enrichments within closely related populations and in single cells. Our study provides a comprehensive framework for single-variant and single-cell analyses of genetic associations.

Fine mapping of blood cell traits in UK Biobank identifies putative causal variants and enrichment of fine-mapped variants in accessible chromatin of hematopoietic progenitor cells. The study provides an analytical framework for single-variant and single-cell analyses of genetic associations.

Evolutionary Pressure against MHC Class II Binding Cancer Mutations

The anti-cancer immune response against mutated peptides of potential immunological relevance (neoantigens) is primarily attributed to MHC-I-restricted cytotoxic CD8+ T cell responses. MHC-II-restricted CD4+ T cells also drive anti-tumor responses, but their relation to neoantigen selection and tumor evolution has not been systematically studied. Modeling the potential of an individual's MHC-II genotype to present 1,018 driver mutations in 5,942 tumors, we demonstrate that the MHC-II genotype constrains the mutational landscape during tumorigenesis in a manner complementary to MHC-I. Mutations poorly bound to MHC-II are positively selected during tumorigenesis, even more than mutations poorly bound to MHC-I. This emphasizes the importance of CD4+ T cells in anti-tumor immunity. In addition, we observed less inter-patient variation in mutation presentation for MHC-II than for MHC-I. These differences were reflected by age at diagnosis, which was correlated with presentation by MHC-I only. Collectively, our results emphasize the central role of MHC-II presentation in tumor evolution.

Prospective associations of C-reactive protein (CRP) levels and CRP genetic risk scores with risk of total knee and hip replacement for osteoarthritis in a diverse cohort

OBJECTIVE

To examine associations of high-sensitivity C-reactive protein (CRP) levels and polygenic CRP genetic risk scores (GRS) with risk of end-stage hip or knee osteoarthritis (OA), defined as incident total hip (THR) or knee replacement (TKR) for OA.

DESIGN

This study included a cohort of postmenopausal white, African American, and Hispanic women from the Women's Health Initiative. Women were followed from baseline to date of THR or TKR, death, or December 31, 2014. Medicare claims data identified THR and TKR. Hs-CRP and genotyping data were collected at baseline. Three CRP GRS were constructed: 1) a 4-SNP GRS comprised of genetic variants representing variation in the CRP gene among European populations; 2) a multilocus 18-SNP GRS of genetic variants significantly associated with CRP levels in a meta-analysis of genome-wide association studies; and 3) a 5-SNP GRS of genetic variants significantly associated with CRP levels among African American women.

RESULTS

In analyses conducted separately among each race and ethnic group, there were no significant associations of ln hs-CRP with risk of THR or TKR, after adjusting for age, body mass index, lifestyle characteristics, chronic diseases, hormone therapy use, and non-steroidal anti-inflammatory drug use. CRP GRS were not associated with risk of THR or TKR in any ethnic group.

CONCLUSIONS

Serum levels of ln hs-CRP and genetically-predicted CRP levels were not associated with risk of THR or TKR for OA among a diverse cohort of women.

A Genome-Wide Association Study of Diabetic Kidney Disease in Subjects With Type 2 Diabetes

Identification of sequence variants robustly associated with predisposition to diabetic kidney disease (DKD) has the potential to provide insights into the pathophysiological mechanisms responsible. We conducted a genome-wide association study (GWAS) of DKD in type 2 diabetes (T2D) using eight complementary dichotomous and quantitative DKD phenotypes: the principal dichotomous analysis involved 5,717 T2D subjects, 3,345 with DKD. Promising association signals were evaluated in up to 26,827 subjects with T2D (12,710 with DKD). A combined T1D+T2D GWAS was performed using complementary data available for subjects with T1D, which, with replication samples, involved up to 40,340 subjects with diabetes (18,582 with DKD). Analysis of specific DKD phenotypes identified a novel signal near GABRR1 (rs9942471, P = 4.5 × 10^-8) associated with microalbuminuria in European T2D case subjects. However, no replication of this signal was observed in Asian subjects with T2D or in the equivalent T1D analysis. There was only limited support, in this substantially enlarged analysis, for association at previously reported DKD signals, except for those at UMOD and PRKAG2, both associated with estimated glomerular filtration rate. We conclude that, despite challenges in addressing phenotypic heterogeneity, access to increased sample sizes will continue to provide more robust inference regarding risk variant discovery for DKD.

The Genetic Landscape of Renal Complications in Type 1 Diabetes

Diabetes is the leading cause of ESRD. Despite evidence for a substantial heritability of diabetic kidney disease, efforts to identify genetic susceptibility variants have had limited success. We extended previous efforts in three dimensions, examining a more comprehensive set of genetic variants in larger numbers of subjects with type 1 diabetes characterized for a wider range of cross-sectional diabetic kidney disease phenotypes. In 2843 subjects, we estimated that the heritability of diabetic kidney disease was 35% (P=6.4×10^-3). Genome-wide association analysis and replication in 12,540 individuals identified no single variants reaching stringent levels of significance and, despite excellent power, provided little independent confirmation of previously published associated variants. Whole-exome sequencing in 997 subjects failed to identify any large-effect coding alleles of lower frequency influencing the risk of diabetic kidney disease. However, sets of alleles increasing body mass index (P=2.2×10^-5) and the risk of type 2 diabetes (P=6.1×10^-4) associated with the risk of diabetic kidney disease. We also found genome-wide genetic correlation between diabetic kidney disease and failure at smoking cessation (P=1.1×10^-4). Pathway analysis implicated ascorbate and aldarate metabolism (P=9.0×10^-6), and pentose and glucuronate interconversions (P=3.0×10^-6) in pathogenesis of diabetic kidney disease. These data provide further evidence for the role of genetic factors influencing diabetic kidney disease in those with type 1 diabetes and highlight some key pathways that may be responsible. Altogether these results reveal important biology behind the major cause of kidney disease.

Determinants of Power in Gene-Based Burden Testing for Monogenic Disorders

Whole-exome sequencing has enabled new approaches for discovering genes associated with monogenic disorders. One such approach is gene-based burden testing, in which the aggregate frequency of "qualifying variants" is compared between case and control subjects for each gene. Despite substantial successes of this approach, the genetic causes for many monogenic disorders remain unknown or only partially known. It is possible that particular genetic architectures lower rates of discovery, but the influence of these factors on power has not been rigorously evaluated. Here, we leverage large-scale exome-sequencing data to create an empirically based simulation framework to evaluate the impact of key parameters (background variation rates, locus heterogeneity, mode of inheritance, penetrance) on power in gene-based burden tests in the context of monogenic disorders. Our results demonstrate that across genes, there is a wide range in sample sizes needed to achieve power due to differences in the background rate of rare variants in each gene. Increasing locus heterogeneity results in rapid increases in sample sizes needed to achieve adequate power, particularly when individual genes contribute to less than 5% of cases under a dominant model. Interestingly, incomplete penetrance as low as 10% had little effect on power due to the low prevalence of monogenic disorders. Our results suggest that moderate incomplete penetrance is not an obstacle in this gene-based burden testing approach but that dominant disorders with high locus heterogeneity will require large sample sizes. Our simulations also provide guidance on sample size needs and inform study design under various genetic architectures.

Recurring exon deletions in the HP (haptoglobin) gene contribute to lower blood cholesterol levels

Two exons of the human haptoglobin (HP) gene exhibit copy number variation that affects HP multimerization and underlies one of the first protein polymorphisms identified in humans. The evolutionary origins and medical significance of this polymorphism have been uncertain. Here we show that this variation has likely arisen from the recurring reversion of an ancient hominin-specific duplication of these exons. Though this polymorphism has been largely invisible to genome-wide genetic studies to date, we describe a way to analyze it by imputation from SNP haplotypes and find among 22,288 individuals that these HP exonic deletions associate with reduced LDL and total cholesterol levels. We show that these deletions, and a SNP that affects HP expression, are the likely drivers of the strong but complex association of cholesterol levels to SNPs near HP. Recurring exonic deletions in the haptoglobin gene likely enhance human health by lowering cholesterol levels in the blood.

Genetic Evidence for a Causal Role of Obesity in Diabetic Kidney Disease

Obesity has been posited as an independent risk factor for diabetic kidney disease (DKD), but establishing causality from observational data is problematic. We aimed to test whether obesity is causally related to DKD using Mendelian randomization, which exploits the random assortment of genes during meiosis. In 6,049 subjects with type 1 diabetes, we used a weighted genetic risk score (GRS) comprised of 32 validated BMI loci as an instrument to test the relationship of BMI with macroalbuminuria, end-stage renal disease (ESRD), or DKD defined as presence of macroalbuminuria or ESRD. We compared these results with cross-sectional and longitudinal observational associations. Longitudinal analysis demonstrated a U-shaped relationship of BMI with development of macroalbuminuria, ESRD, or DKD over time. Cross-sectional observational analysis showed no association with overall DKD, higher odds of macroalbuminuria (for every 1 kg/m(2) higher BMI, odds ratio [OR] 1.05, 95% CI 1.03-1.07, P < 0.001), and lower odds of ESRD (OR 0.95, 95% CI 0.93-0.97, P < 0.001). Mendelian randomization analysis showed a 1 kg/m(2) higher BMI conferring an increased risk in macroalbuminuria (OR 1.28, 95% CI 1.11-1.45, P = 0.001), ESRD (OR 1.43, 95% CI 1.20-1.72, P < 0.001), and DKD (OR 1.33, 95% CI 1.17-1.51, P < 0.001). Our results provide genetic evidence for a causal link between obesity and DKD in type 1 diabetes. As obesity prevalence rises, this finding predicts an increase in DKD prevalence unless intervention should occur.

Targeted Application of Human Genetic Variation Can Improve Red Blood Cell Production from Stem Cells

Multipotent and pluripotent stem cells are potential sources for cell and tissue replacement therapies. For example, stem cell-derived red blood cells (RBCs) are a potential alternative to donated blood, but yield and quality remain a challenge. Here, we show that application of insight from human population genetic studies can enhance RBC production from stem cells. The SH2B3 gene encodes a negative regulator of cytokine signaling and naturally occurring loss-of-function variants in this gene increase RBC counts in vivo. Targeted suppression of SH2B3 in primary human hematopoietic stem and progenitor cells enhanced the maturation and overall yield of in-vitro-derived RBCs. Moreover, inactivation of SH2B3 by CRISPR/Cas9 genome editing in human pluripotent stem cells allowed enhanced erythroid cell expansion with preserved differentiation. Our findings therefore highlight the potential for combining human genome variation studies with genome editing approaches to improve cell and tissue production for regenerative medicine.

Genetic risk factors affecting mitochondrial function are associated with kidney disease in people with Type 1 diabetes

AIM

To evaluate the association with diabetic kidney disease of single nucleotide polymorphisms (SNPs) that may contribute to mitochondrial dysfunction.

METHODS

The mitochondrial genome and 1039 nuclear genes that are integral to mitochondrial function were investigated using a case (n = 823 individuals with diabetic kidney disease) vs. control (n = 903 individuals with diabetes and no renal disease) approach. All people included in the analysis were of white European origin and were diagnosed with Type 1 diabetes before the age of 31 years. Replication was conducted in 5093 people with similar phenotypes to those of the discovery collection. Association analyses were performed using the plink genetic analysis toolset, with adjustment for relevant covariates.

RESULTS

A total of 25 SNPs were evaluated in the mitochondrial genome, but none were significantly associated with diabetic kidney disease or end-stage renal disease. A total of 38 SNPs in nuclear genes influencing mitochondrial function were nominally associated with diabetic kidney disease and 16 SNPS were associated with end-stage renal disease, secondary to diabetic kidney disease, with meta-analyses confirming the same direction of effect. Three independent signals (seven SNPs) were common to the replication data for both phenotypes with Type 1 diabetes and persistent proteinuria or end-stage renal disease.

CONCLUSIONS

Our results suggest that SNPs in nuclear genes that influence mitochondrial function are significantly associated with diabetic kidney disease in a white European population.

Directional dominance on stature and cognition in diverse human populations

Homozygosity has long been associated with rare, often devastating, Mendelian disorders, and Darwin was one of the first to recognize that inbreeding reduces evolutionary fitness. However, the effect of the more distant parental relatedness that is common in modern human populations is less well understood. Genomic data now allow us to investigate the effects of homozygosity on traits of public health importance by observing contiguous homozygous segments (runs of homozygosity), which are inferred to be homozygous along their complete length. Given the low levels of genome-wide homozygosity prevalent in most human populations, information is required on very large numbers of people to provide sufficient power. Here we use runs of homozygosity to study 16 health-related quantitative traits in 354,224 individuals from 102 cohorts, and find statistically significant associations between summed runs of homozygosity and four complex traits: height, forced expiratory lung volume in one second, general cognitive ability and educational attainment (P < 1 × 10(-300), 2.1 × 10(-6), 2.5 × 10(-10) and 1.8 × 10(-10), respectively). In each case, increased homozygosity was associated with decreased trait value, equivalent to the offspring of first cousins being 1.2 cm shorter and having 10 months' less education. Similar effect sizes were found across four continental groups and populations with different degrees of genome-wide homozygosity, providing evidence that homozygosity, rather than confounding, directly contributes to phenotypic variance. Contrary to earlier reports in substantially smaller samples, no evidence was seen of an influence of genome-wide homozygosity on blood pressure and low density lipoprotein cholesterol, or ten other cardio-metabolic traits. Since directional dominance is predicted for traits under directional evolutionary selection, this study provides evidence that increased stature and cognitive function have been positively selected in human evolution, whereas many important risk factors for late-onset complex diseases may not have been.

New genetic loci link adipose and insulin biology to body fat distribution

Body fat distribution is a heritable trait and a well-established predictor of adverse metabolic outcomes, independent of overall adiposity. To increase our understanding of the genetic basis of body fat distribution and its molecular links to cardiometabolic traits, we conducted genome-wide association meta-analyses of waist and hip circumference-related traits in up to 224,459 individuals. We identified 49 loci (33 new) associated with waist-to-hip ratio adjusted for body mass index (WHRadjBMI) and an additional 19 loci newly associated with related waist and hip circumference measures (P<5×10⁻⁸). Twenty of the 49 WHRadjBMI loci showed significant sexual dimorphism, 19 of which displayed a stronger effect in women. The identified loci were enriched for genes expressed in adipose tissue and for putative regulatory elements in adipocytes. Pathway analyses implicated adipogenesis, angiogenesis, transcriptional regulation, and insulin resistance as processes affecting fat distribution, providing insight into potential pathophysiological mechanisms.

Efficient Bayesian mixed-model analysis increases association power in large cohorts

Linear mixed models are a powerful statistical tool for identifying genetic associations and avoiding confounding. However, existing methods are computationally intractable in large cohorts and may not optimize power. All existing methods require time cost O(MN(2)) (where N is the number of samples and M is the number of SNPs) and implicitly assume an infinitesimal genetic architecture in which effect sizes are normally distributed, which can limit power. Here we present a far more efficient mixed-model association method, BOLT-LMM, which requires only a small number of O(MN) time iterations and increases power by modeling more realistic, non-infinitesimal genetic architectures via a Bayesian mixture prior on marker effect sizes. We applied BOLT-LMM to 9 quantitative traits in 23,294 samples from the Women's Genome Health Study (WGHS) and observed significant increases in power, consistent with simulations. Theory and simulations show that the boost in power increases with cohort size, making BOLT-LMM appealing for genome-wide association studies in large cohorts.

Novel genetic susceptibility loci for diabetic end-stage renal disease identified through robust naive Bayes classification

AIMS/HYPOTHESIS

Diabetic nephropathy is a major diabetic complication, and diabetes is the leading cause of end-stage renal disease (ESRD). Family studies suggest a hereditary component for diabetic nephropathy. However, only a few genes have been associated with diabetic nephropathy or ESRD in diabetic patients. Our aim was to detect novel genetic variants associated with diabetic nephropathy and ESRD.

METHODS

We exploited a novel algorithm, 'Bag of Naive Bayes', whose marker selection strategy is complementary to that of conventional genome-wide association models based on univariate association tests. The analysis was performed on a genome-wide association study of 3,464 patients with type 1 diabetes from the Finnish Diabetic Nephropathy (FinnDiane) Study and subsequently replicated with 4,263 type 1 diabetes patients from the Steno Diabetes Centre, the All Ireland-Warren 3-Genetics of Kidneys in Diabetes UK collection (UK-Republic of Ireland) and the Genetics of Kidneys in Diabetes US Study (GoKinD US).

RESULTS

Five genetic loci (WNT4/ZBTB40-rs12137135, RGMA/MCTP2-rs17709344, MAPRE1P2-rs1670754, SEMA6D/SLC24A5-rs12917114 and SIK1-rs2838302) were associated with ESRD in the FinnDiane study. An association between ESRD and rs17709344, tagging the previously identified rs12437854 and located between the RGMA and MCTP2 genes, was replicated in independent case-control cohorts. rs12917114 near SEMA6D was associated with ESRD in the replication cohorts under the genotypic model (p < 0.05), and rs12137135 upstream of WNT4 was associated with ESRD in Steno.

CONCLUSIONS/INTERPRETATION

This study supports the previously identified findings on the RGMA/MCTP2 region and suggests novel susceptibility loci for ESRD. This highlights the importance of applying complementary statistical methods to detect novel genetic variants in diabetic nephropathy and, in general, in complex diseases.

Distribution and medical impact of loss-of-function variants in the Finnish founder population

Exome sequencing studies in complex diseases are challenged by the allelic heterogeneity, large number and modest effect sizes of associated variants on disease risk and the presence of large numbers of neutral variants, even in phenotypically relevant genes. Isolated populations with recent bottlenecks offer advantages for studying rare variants in complex diseases as they have deleterious variants that are present at higher frequencies as well as a substantial reduction in rare neutral variation. To explore the potential of the Finnish founder population for studying low-frequency (0.5-5%) variants in complex diseases, we compared exome sequence data on 3,000 Finns to the same number of non-Finnish Europeans and discovered that, despite having fewer variable sites overall, the average Finn has more low-frequency loss-of-function variants and complete gene knockouts. We then used several well-characterized Finnish population cohorts to study the phenotypic effects of 83 enriched loss-of-function variants across 60 phenotypes in 36,262 Finns. Using a deep set of quantitative traits collected on these cohorts, we show 5 associations (p<5×10⁻⁸) including splice variants in LPA that lowered plasma lipoprotein(a) levels (P = 1.5×10⁻¹¹⁷). Through accessing the national medical records of these participants, we evaluate the LPA finding via Mendelian randomization and confirm that these splice variants confer protection from cardiovascular disease (OR = 0.84, P = 3×10⁻⁴), demonstrating for the first time the correlation between very low levels of LPA in humans with potential therapeutic implications for cardiovascular diseases. More generally, this study articulates substantial advantages for studying the role of rare variation in complex phenotypes in founder populations like the Finns and by combining a unique population genetic history with data from large population cohorts and centralized research access to National Health Registers.

An excess of risk-increasing low-frequency variants can be a signal of polygenic inheritance in complex diseases

In most complex diseases, much of the heritability remains unaccounted for by common variants. It has been postulated that lower-frequency variants contribute to the remaining heritability. Here, we describe a method to test for polygenic inheritance from lower-frequency variants by using GWAS summary association statistics. We explored scenarios with many causal low-frequency variants and showed that there is more power to detect risk variants than to detect protective variants, resulting in an increase in the ratio of detected risk to protective variants (R/P ratio). Such an excess can also occur if risk variants are present and kept at lower frequencies because of negative selection. The R/P ratio can be falsely elevated because of reasons unrelated to polygenic inheritance, such as uneven sample sizes or asymmetric population stratification, so precautions to correct for these confounders are essential. We tested our method on published GWAS results and observed a strong signal in some diseases (schizophrenia and type 2 diabetes) but not others. We also explored the shared genetic component in overlapping phenotypes related to inflammatory bowel disease (Crohn disease [CD] and ulcerative colitis [UC]) and diabetic nephropathy (macroalbuminuria and end-stage renal disease [ESRD]). Although the signal was still present when both CD and UC were jointly analyzed, the signal was lost when macroalbuminuria and ESRD were jointly analyzed, suggesting that these phenotypes should best be studied separately. Thus, our method may also help guide the design of future genetic studies of various traits and diseases.

Chromosome 2q31.1 associates with ESRD in women with type 1 diabetes

Sex and genetic variation influence the risk of developing diabetic nephropathy and ESRD in patients with type 1 diabetes. We performed a genome-wide association study in a cohort of 3652 patients from the Finnish Diabetic Nephropathy (FinnDiane) Study with type 1 diabetes to determine whether sex-specific genetic risk factors for ESRD exist. A common variant, rs4972593 on chromosome 2q31.1, was associated with ESRD in women (P<5×10(-8)) but not in men (P=0.77). This association was replicated in the meta-analysis of three independent type 1 diabetes cohorts (P=0.02) and remained significant for women (P<5×10(-8); odds ratio, 1.81 [95% confidence interval, 1.47 to 2.24]) upon combined meta-analysis of the discovery and replication cohorts. rs4972593 is located between the genes that code for the Sp3 transcription factor, which interacts directly with estrogen receptor α and regulates the expression of genes linked to glomerular function and the pathogenesis of nephropathy, and the CDCA7 transcription factor, which regulates cell proliferation. Further examination revealed potential transcription factor-binding sites within rs4972593 and predicted eight estrogen-responsive elements within 5 kb of this locus. Moreover, we found sex-specific differences in the glomerular expression levels of SP3 (P=0.004). Overall, these results suggest that rs4972593 is a sex-specific genetic variant associated with ESRD in patients with type 1 diabetes and may underlie the sex-specific protection against ESRD.

MicroRNA-22 and promoter motif polymorphisms at the Chga locus in genetic hypertension: functional and therapeutic implications for gene expression and the pathogenesis of hypertension

Hypertension is a common hereditary syndrome with unclear pathogenesis. Chromogranin A (Chga), which catalyzes formation and cargo storage of regulated secretory granules in neuroendocrine cells, contributes to blood pressure homeostasis centrally and peripherally. Elevated Chga occurs in spontaneously hypertensive rat (SHR) adrenal glands and plasma, but central expression is unexplored. In this report, we measured SHR and Wistar-Kyoto rat (control) Chga expression in central and peripheral nervous systems, and found Chga protein to be decreased in the SHR brainstem, yet increased in the adrenal and the plasma. By re-sequencing, we systematically identified five promoter, two coding and one 3'-untranslated region (3'-UTR) polymorphism at the SHR (versus WKY or BN) Chga locus. Using HXB/BXH recombinant inbred (RI) strain linkage and correlations, we demonstrated genetic determination of Chga expression in SHR, including a cis-quantitative trait loci (QTLs) (i.e. at the Chga locus), and such expression influenced biochemical determinants of blood pressure, including a cascade of catecholamine biosynthetic enzymes, catecholamines themselves and steroids. Luciferase reporter assays demonstrated that the 3'-UTR polymorphism (which disrupts a microRNA miR-22 motif) and promoter polymorphisms altered gene expression consistent with the decline in SHR central Chga expression. Coding region polymorphisms did not account for changes in Chga expression or function. Thus, we hypothesized that the 3'-UTR and promoter mutations lead to dysregulation (diminution) of Chga in brainstem cardiovascular control nuclei, ultimately contributing to the pathogenesis of hypertension in SHR. Accordingly, we demonstrated that in vivo administration of miR-22 antagomir to SHR causes substantial (∼18 mmHg) reductions in blood pressure, opening a novel therapeutic avenue for hypertension.

Association testing of previously reported variants in a large case-control meta-analysis of diabetic nephropathy

We formed the GEnetics of Nephropathy-an International Effort (GENIE) consortium to examine previously reported genetic associations with diabetic nephropathy (DN) in type 1 diabetes. GENIE consists of 6,366 similarly ascertained participants of European ancestry with type 1 diabetes, with and without DN, from the All Ireland-Warren 3-Genetics of Kidneys in Diabetes U.K. and Republic of Ireland (U.K.-R.O.I.) collection and the Finnish Diabetic Nephropathy Study (FinnDiane), combined with reanalyzed data from the Genetics of Kidneys in Diabetes U.S. Study (U.S. GoKinD). We found little evidence for the association of the EPO promoter polymorphism, rs161740, with the combined phenotype of proliferative retinopathy and end-stage renal disease in U.K.-R.O.I. (odds ratio [OR] 1.14, P = 0.19) or FinnDiane (OR 1.06, P = 0.60). However, a fixed-effects meta-analysis that included the previously reported cohorts retained a genome-wide significant association with that phenotype (OR 1.31, P = 2 × 10(-9)). An expanded investigation of the ELMO1 locus and genetic regions reported to be associated with DN in the U.S. GoKinD yielded only nominal statistical significance for these loci. Finally, top candidates identified in a recent meta-analysis failed to reach genome-wide significance. In conclusion, we were unable to replicate most of the previously reported genetic associations for DN, and significance for the EPO promoter association was attenuated.

A genome-wide association meta-analysis identifies new childhood obesity loci

Multiple genetic variants have been associated with adult obesity and a few with severe obesity in childhood; however, less progress has been made in establishing genetic influences on common early-onset obesity. We performed a North American, Australian and European collaborative meta-analysis of 14 studies consisting of 5,530 cases (≥95th percentile of body mass index (BMI)) and 8,318 controls (<50th percentile of BMI) of European ancestry. Taking forward the eight newly discovered signals yielding association with P < 5 × 10(-6) in nine independent data sets (2,818 cases and 4,083 controls), we observed two loci that yielded genome-wide significant combined P values near OLFM4 at 13q14 (rs9568856; P = 1.82 × 10(-9); odds ratio (OR) = 1.22) and within HOXB5 at 17q21 (rs9299; P = 3.54 × 10(-9); OR = 1.14). Both loci continued to show association when two extreme childhood obesity cohorts were included (2,214 cases and 2,674 controls). These two loci also yielded directionally consistent associations in a previous meta-analysis of adult BMI(1).

The VA Hypertension Primary Care Longitudinal Cohort: Electronic medical records in the post-genomic era

The Veterans Affairs Hypertension Primary Care Longitudinal Cohort (VAHC) was initiated in 2003 as a pilot study designed to link the VA electronic medical record system with individual genetic data. Between June 2003 and December 2004, 1,527 hypertensive participants were recruited. Protected health information (PHI) was extracted from the regional VA data warehouse. Differences between the clinic and mail recruits suggested that clinic recruitment resulted in an over-sampling of African Americans. A review of medical records in a random sample of study participants confirmed that the data warehouse accurately captured most selected diagnoses. Genomic DNA was acquired non-invasively from buccal cells in mouthwash; ~ 96.5 per cent of samples contained DNA suitable for genotyping, with an average DNA yield of 5.02 ± 0.12 micrograms, enough for several thousand genotypes. The coupling of detailed medical databases with genetic information has the potential to facilitate the genetic study of hypertension and other complex diseases.

Methylenetetrahydrofolate reductase (MTHFR) polymorphism A1298C (Glu429Ala) predicts decline in renal function over time in the African-American Study of Kidney Disease and Hypertension (AASK) Trial and Veterans Affairs Hypertension Cohort (VAHC)

BACKGROUND

Hyperhomocysteinemia is associated with increased venous thrombosis and cardiovascular disease (CVD). Mutations in the human methylenetetrahydrofolate reductase (MTHFR) gene have been associated with increased homocysteine levels and risks of CVD in various populations including those with kidney disease. Here, we evaluated the influence of MTHFR variants on progressive loss of kidney function.

METHODS

We analyzed 821 subjects with hypertensive nephrosclerosis from the longitudinal National Institute of Diabetes and Digestive and Kidney Diseases African-American Study of Kidney Disease and Hypertension (AASK) Trial to determine whether decline in glomerular filtration rate (GFR) over ∼4.2 years was predicted by common genetic variation within MTHFR at non-synonymous positions C677T (Ala222Val) and A1298C (Glu429Ala) or by MTHFR haplotypes. The effect on GFR decline was then supported by a study of 1333 subjects from the San Diego Veterans Affairs Hypertension Cohort (VAHC), followed over ∼4.5 years. Linear effect models were utilized to determine both genotype [single-nucleotide polymorphism (SNP)] and genotype (SNP)-by-time interactions.

RESULTS

In AASK, the polymorphism at A1298C predicted the rate of GFR decline: A1298/A1298 major allele homozygosity resulted in a less pronounced decline of GFR, with a significant SNP-by-time interaction. An independent follow-up study in the San Diego VAHC subjects supports that A1298/A1298 homozygotes have the greatest estimated GFR throughout the study. Haplotype analysis with C677T yielded concurring results.

CONCLUSION

We conclude that the MTHFR-coding polymorphism at A1298C is associated with renal decline in African-Americans with hypertensive nephrosclerosis and is supported by a veteran cohort with a primary care diagnosis of hypertension. Further investigation is needed to confirm such findings and to determine what molecular mechanism may contribute to this association.

Direct vasoactive effects of the chromogranin A (CHGA) peptide catestatin in humans in vivo

Catestatin is a bioactive peptide of chromogranin A (CHGA) that is co-released with catecholamines from secretory vesicles. Catestatin may function as a vasodilator and is diminished in hypertension. To evaluate this potential vasodilator in vivo without systemic counterregulation, we infused catestatin to target concentrations of approximately 50, approximately 500, approximately 5000 nM into dorsal hand veins of 18 normotensive men and women, after pharmacologic venoconstriction with phenylephrine. Pancreastatin, another CHGA peptide, was infused as a negative control. After preconstriction to approximately 69%, increasing concentrations of catestatin resulted in dose-dependent vasodilation (P = 0.019), in female subjects (to approximately 44%) predominantly. The EC(50) (approximately 30 nM) for vasodilation induced by catestatin was the same order of magnitude to circulating endogenous catestatin (4.4 nM). No vasodilation occurred during the control infusion with pancreastatin. Plasma CHGA, catestatin, and CHGA-to-catestatin processing were then determined in 622 healthy subjects without hypertension. Female subjects had higher plasma catestatin levels than males (P = 0.001), yet lower CHGA precursor concentrations (P = 0.006), reflecting increased processing of CHGA-to-catestatin (P < 0.001). Our results demonstrate that catestatin dilates human blood vessels in vivo, especially in females. Catestatin may contribute to sex differences in endogenous vascular tone, thereby influencing the complex predisposition to hypertension.

Common charge-shift mutation Glu65Lys in K+ channel β₁-Subunit KCNMB1: pleiotropic consequences for glomerular filtration rate and progressive renal disease

BACKGROUND

Glomerular filtration rate (GFR) is a heritable trait, and hyperfiltration (GFR increment in remnant nephrons) may accelerate renal functional decline in chronic kidney disease (CKD). Mesangial and vascular smooth myocytes control GFR by contraction, dependent on voltage-gated Ca(2+) influx, which is controlled by the regulatory β₁-subunit (KCNMB1) of large-conductance heteromeric K+ ('BK') channels. KCNMB1 gain-of-function variant Glu65Lys results in generalized vasorelaxation and thus protection against systemic hypertension. Here we asked whether the Glu65Lys variant influences GFR, in the basal state or during progressive renal decline.

METHODS

We explored Glu65Lys effects on GFR in three populations spanning two ethnicities and two diseases (hypertension and nephrosclerosis). GFR was either estimated (eGFR from serum creatinine) or directly measured (iothalamate clearance).

RESULTS

The 65Lys variant was relatively common, occurring on ∼5-10% of chromosomes in different biogeographic ancestry groups, and 65Lys carriers exhibited higher eGFR in two primary care populations: extreme BP values in Kaiser clinics (p = 0.029, accounting for ∼0.2% of trait variance), or treated hypertensives in VA clinics (p = 0.017, accounting for ∼0.9% of trait variance). In blacks with progressive renal disease (NIDDK AASK), 65Lys carriers displayed a steeper slope in GFR chronic decline (p = 0.030, accounting for ∼0.4% of trait variance), and Glu65Lys genotype also predicted time of onset of renal failure (log rank p = 0.019).

CONCLUSIONS

Common KCNMB1 gain-of-function variant Glu65Lys influences GFR, and 65Lys carriers exhibit not only elevated baseline GFR, but also more rapid GFR decline (and consequent development of renal failure) in CKD. The results suggest that profiling patients at Glu65Lys can assist in gauging renal prognosis as well as selection of rational therapy in hypertension with progressive renal disease.

Longitudinal genome-wide association of cardiovascular disease risk factors in the Bogalusa heart study

Cardiovascular disease (CVD) is the leading cause of death worldwide. Recent genome-wide association (GWA) studies have pinpointed many loci associated with CVD risk factors in adults. It is unclear, however, if these loci predict trait levels at all ages, if they are associated with how a trait develops over time, or if they could be used to screen individuals who are pre-symptomatic to provide the opportunity for preventive measures before disease onset. We completed a genome-wide association study on participants in the longitudinal Bogalusa Heart Study (BHS) and have characterized the association between genetic factors and the development of CVD risk factors from childhood to adulthood. We report 7 genome-wide significant associations involving CVD risk factors, two of which have been previously reported. Top regions were tested for replication in the Young Finns Study (YF) and two associations strongly replicated: rs247616 in CETP with HDL levels (combined P = 9.7×10⁻²⁴), and rs445925 at APOE with LDL levels (combined P = 8.7×10⁻¹⁹). We show that SNPs previously identified in adult cross-sectional studies tend to show age-independent effects in the BHS with effect sizes consistent with previous reports. Previously identified variants were associated with adult trait levels above and beyond those seen in childhood; however, variants with time-dependent effects were also promising predictors. This is the first GWA study to evaluate the role of common genetic variants in the development of CVD risk factors in children as they advance through adulthood and highlights the utility of using longitudinal studies to identify genetic predictors of adult traits in children.

We have studied the association between genetic factors on a whole genome level and cardiovascular disease (CVD) risk factors in a population of individuals studied from childhood through adulthood. The longitudinal study design has enabled the investigation of genetic variation influencing trait values over time. We have identified DNA variants that are associated with CVD trait values consistently over time, and a second set of variants that are associated with CVD trait values in a time-dependent manner. We also show that variants previously identified in adult populations have consistent effects within our population and that these effects are usually similar across childhood through adulthood. The discovery of time-dependent variants that influence CVD trait values over time can potentially be used to screen young individuals who are pre-symptomatic and provide the opportunity for preventive measures decades before disease onset.

Curve-based multivariate distance matrix regression analysis: application to genetic association analyses involving repeated measures

Most, if not all, human phenotypes exhibit a temporal, dosage-dependent, or age effect. Despite this fact, it is rare that data are collected over time or in sequence in relevant studies of the determinants of these phenotypes. The costs and organizational sophistication necessary to collect repeated measurements or longitudinal data for a given phenotype are clearly impediments to this, but greater efforts in this area are needed if insights into human phenotypic expression are to be obtained. Appropriate data analysis methods for genetic association studies involving repeated or longitudinal measures are also needed. We consider the use of longitudinal profiles obtained from fitted functions on repeated data collections from a set of individuals whose similarities are contrasted between sets of individuals with different genotypes to test hypotheses about genetic influences on time-dependent phenotype expression. The proposed approach can accommodate uncertainty of the fitted functions, as well as weighting factors across the time points, and is easily extended to a wide variety of complex analysis settings. We showcase the proposed approach with data from a clinical study investigating human blood vessel response to tyramine. We also compare the proposed approach with standard analytic procedures and investigate its robustness and power via simulation studies. The proposed approach is found to be quite flexible and performs either as well or better than traditional statistical methods.