GWAS for the composite traits of hematuria and albuminuria.

Our GWAS of hematuria in the UK Biobank identified 6 loci, some of which overlap with loci for albuminuria suggesting pleiotropy. Since clinical syndromes are often defined by combinations of traits, generating a combined phenotype can improve power to detect loci influencing multiple characteristics. Thus the composite trait of hematuria and albuminuria was chosen to enrich for glomerular pathologies. Cases had both hematuria defined by ICD codes and albuminuria defined as uACR > 3 mg/mmol. Controls had neither an ICD code for hematuria nor an uACR > 3 mg/mmol. 2429 cases and 343,509 controls from the UK Biobank were included. eGFR was lower in cases compared to controls, with the exception of the comparison in females using CKD-EPI after age adjustment. Variants at 4 loci met genome-wide significance with the following nearest genes: COL4A4, TRIM27, ETV1 and CUBN. TRIM27 is part of the extended MHC locus. All loci with the exception of ETV1 were replicated in the Geisinger MyCode cohort. The previous GWAS of hematuria reported COL4A3-COL4A4 variants and HLA-B*0801 within MHC, which is in linkage disequilibrium with the TRIM27 variant (D' = 0.59). TRIM27 is highly expressed in the tubules. Additional loci included a coding sequence variant in CUBN (p.Ala2914Val, MAF = 0.014 (A), p = 3.29E-8, OR = 2.09, 95% CI = 1.61-2.72). Overall, GWAS for the composite trait of hematuria and albuminuria identified 4 loci, 2 of which were not previously identified in a GWAS of hematuria.

Genetically predicted waist-to-hip circumference ratio and coronary artery disease: A sex-specific Mendelian randomization study.

Coronary artery disease (CAD) affects millions of individuals worldwide and results in a substantial burden to healthcare systems. Although it is established that CAD affects females differently than males, differences between the sexes are not routinely accounted for. Body mass index is a known risk factor for CAD. However, more accurate metrics of body fat, including waist-to-hip circumference ratio (WHR), could be more meaningful clinically. WHR exhibits sex differences due to sex hormones, differing effects at genetic risk loci, and other factors. It is unclear if WHR is a causal factor for CAD in one or both sexes, but this information will be crucial for improving heart health. Causal inference, however, can be challenging. Large-scale cohorts with genetic data allow for Mendelian randomization, which, given certain assumptions, tests whether there is a causal relationship between an exposure and the outcome using genetic variants. We conducted sex-specific, one-sample MR analyses using two-stage least-squares regression in the UK Biobank with genetic variants robustly associated with WHR. We found evidence of a causal relationship between WHR and CAD risk in females (OR [95% CI] = 1.16 [1.06-1.26]; p value = 7.5E-4), whereas in males, we did not find evidence of a causal relationship (OR [95% CI] = 1.40 [0.98-2.01]; p value = 0.063). Results were supported by two additional MR approaches (using a genetic risk score and two-sample MR using the inverse variance weighted approach). We encourage future work assessing sex-specific effects using causal inference techniques to better understand factors contributing to complex disease risk.

Regional genetic correlations highlight relationships between neurodegenerative disease loci and the immune system.

Neurodegenerative diseases, including Alzheimer's and Parkinson's disease, are devastating complex diseases resulting in physical and psychological burdens on patients and their families. There have been important efforts to understand their genetic basis leading to the identification of disease risk-associated loci involved in several molecular mechanisms, including immune-related pathways. Regional, in contrast to genome-wide, genetic correlations between pairs of immune and neurodegenerative traits have not been comprehensively explored, but could uncover additional immune-mediated risk-associated loci. Here, we systematically assess the role of the immune system in five neurodegenerative diseases by estimating regional genetic correlations between these diseases and immune-cell-derived single-cell expression quantitative trait loci (sc-eQTLs). We also investigate correlations between diseases and protein levels. We observe significant (FDR < 0.01) correlations between sc-eQTLs and neurodegenerative diseases across 151 unique genes, spanning both the innate and adaptive immune systems, across most diseases tested. With Parkinson's, for instance, RAB7L1 in CD4+ naïve T cells is positively correlated and KANSL1-AS1 is negatively correlated across all adaptive immune cell types. Follow-up colocalization highlight candidate causal risk genes. The outcomes of this study will improve our understanding of the immune component of neurodegeneration, which can warrant repurposing of existing immunotherapies to slow disease progression.

Oligogenic Architecture of Rare Noncoding Variants Distinguishes 4 Congenital Heart Disease Phenotypes.

BACKGROUND: Congenital heart disease (CHD) is highly heritable, but the power to identify inherited risk has been limited to analyses of common variants in small cohorts. METHODS: We performed reimputation of 4 CHD cohorts (n=55 342) to the TOPMed reference panel (freeze 5), permitting meta-analysis of 14 784 017 variants including 6 035 962 rare variants of high imputation quality as validated by whole genome sequencing. RESULTS: Meta-analysis identified 16 novel loci, including 12 rare variants, which displayed moderate or large effect sizes (median odds ratio, 3.02) for 4 separate CHD categories. Analyses of chromatin structure link 13 of the genome-wide significant loci to key genes in cardiac development; rs373447426 (minor allele frequency, 0.003 [odds ratio, 3.37 for Conotruncal heart disease]; P=1.49×10-8) is predicted to disrupt chromatin structure for 2 nearby genes BDH1 and DLG1 involved in Conotruncal development. A lead variant rs189203952 (minor allele frequency, 0.01 [odds ratio, 2.4 for left ventricular outflow tract obstruction]; P=1.46×10-8) is predicted to disrupt the binding sites of 4 transcription factors known to participate in cardiac development in the promoter of SPAG9. A tissue-specific model of chromatin conformation suggests that common variant rs78256848 (minor allele frequency, 0.11 [odds ratio, 1.4 for Conotruncal heart disease]; P=2.6×10-8) physically interacts with NCAM1 (PFDR=1.86×10-27), a neural adhesion molecule acting in cardiac development. Importantly, while each individual malformation displayed substantial heritability (observed h2 ranging from 0.26 for complex malformations to 0.37 for left ventricular outflow tract obstructive disease) the risk for different CHD malformations appeared to be separate, without genetic correlation measured by linkage disequilibrium score regression or regional colocalization. CONCLUSIONS: We describe a set of rare noncoding variants conferring significant risk for individual heart malformations which are linked to genes governing cardiac development. These results illustrate that the oligogenic basis of CHD and significant heritability may be linked to rare variants outside protein-coding regions conferring substantial risk for individual categories of cardiac malformation.

Local genetic correlations exist among neurodegenerative and neuropsychiatric diseases.

Genetic correlation ([Formula: see text]) between traits can offer valuable insight into underlying shared biological mechanisms. Neurodegenerative diseases overlap neuropathologically and often manifest comorbid neuropsychiatric symptoms. However, global [Formula: see text] analyses show minimal [Formula: see text] among neurodegenerative and neuropsychiatric diseases. Importantly, local [Formula: see text] s can exist in the absence of global relationships. To investigate this possibility, we applied LAVA, a tool for local [Formula: see text] analysis, to genome-wide association studies of 3 neurodegenerative diseases (Alzheimer's disease, Lewy body dementia and Parkinson's disease) and 3 neuropsychiatric disorders (bipolar disorder, major depressive disorder and schizophrenia). We identified several local [Formula: see text] s missed in global analyses, including between (i) all 3 neurodegenerative diseases and schizophrenia and (ii) Alzheimer's and Parkinson's disease. For those local [Formula: see text] s identified in genomic regions containing disease-implicated genes, such as SNCA, CLU and APOE, incorporation of expression quantitative trait loci identified genes that may drive genetic overlaps between diseases. Collectively, we demonstrate that complex genetic relationships exist among neurodegenerative and neuropsychiatric diseases, highlighting putative pleiotropic genomic regions and genes. These findings imply sharing of pathogenic processes and the potential existence of common therapeutic targets.

The contribution of Neanderthal introgression and natural selection to neurodegenerative diseases.

Humans are thought to be more susceptible to neurodegeneration than equivalently-aged primates. It is not known whether this vulnerability is specific to anatomically-modern humans or shared with other hominids. The contribution of introgressed Neanderthal DNA to neurodegenerative disorders remains uncertain. It is also unclear how common variants associated with neurodegenerative disease risk are maintained by natural selection in the population despite their deleterious effects. In this study, we aimed to quantify the genome-wide contribution of Neanderthal introgression and positive selection to the heritability of complex neurodegenerative disorders to address these questions. We used stratified-linkage disequilibrium score regression to investigate the relationship between five SNP-based signatures of natural selection, reflecting different timepoints of evolution, and genome-wide associated variants of the three most prevalent neurodegenerative disorders: Alzheimer's disease, amyotrophic lateral sclerosis and Parkinson's disease. We found no evidence for enrichment of positively-selected SNPs in the heritability of Alzheimer's disease, amyotrophic lateral sclerosis and Parkinson's disease, suggesting that common deleterious disease variants are unlikely to be maintained by positive selection. There was no enrichment of Neanderthal introgression in the SNP-heritability of these disorders, suggesting that Neanderthal admixture is unlikely to have contributed to disease risk. These findings provide insight into the origins of neurodegenerative disorders within the evolution of Homo sapiens and addresses a long-standing debate, showing that Neanderthal admixture is unlikely to have contributed to common genetic risk of neurodegeneration in anatomically-modern humans.

A fast linkage method for population GWAS cohorts with related individuals.

Linkage analysis, a class of methods for detecting co-segregation of genomic segments and traits in families, was used to map disease-causing genes for decades before genotyping arrays and dense SNP genotyping enabled genome-wide association studies in population samples. Population samples often contain related individuals, but the segregation of alleles within families is rarely used because traditional linkage methods are computationally inefficient for larger datasets. Here, we describe Population Linkage, a novel application of Haseman-Elston regression as a method of moments estimator of variance components and their standard errors. We achieve additional computational efficiency by using modern methods for detection of IBD segments and variance component estimation, efficient preprocessing of input data, and minimizing redundant numerical calculations. We also refined variance component models to account for the biases in population-scale methods for IBD segment detection. We ran Population Linkage on four blood lipid traits in over 70,000 individuals from the HUNT and SardiNIA studies, successfully detecting 25 known genetic signals. One notable linkage signal that appeared in both was for low-density lipoprotein (LDL) cholesterol levels in the region near the gene APOE (LOD = 29.3, variance explained = 4.1%). This is the region where the missense variants rs7412 and rs429358, which together make up the ε2, ε3, and ε4 alleles each account for 2.4% and 0.8% of variation in circulating LDL cholesterol. Our results show the potential for linkage analysis and other large-scale applications of method of moments variance components estimation.

Evaluation of polygenic risk scores to differentiate between type 1 and type 2 diabetes.

Polygenic risk scores (PRS) quantify the genetic liability to disease and are calculated using an individual's genotype profile and disease-specific genome-wide association study (GWAS) summary statistics. Type 1 (T1D) and type 2 (T2D) diabetes both are determined in part by genetic loci. Correctly differentiating between types of diabetes is crucial for accurate diagnosis and treatment. PRS have the potential to address possible misclassification of T1D and T2D. Here we evaluated PRS models for T1D and T2D in European genetic ancestry participants from the UK Biobank (UKB) and then in the Michigan Genomics Initiative (MGI). Specifically, we investigated the utility of T1D and T2D PRS to discriminate between T1D, T2D, and controls in unrelated UKB individuals of European ancestry. We derived PRS models using external non-UKB GWAS. The T1D PRS model with the best discrimination between T1D cases and controls (area under the receiver operator curve [AUC] = 0.805) also yielded the best discrimination of T1D from T2D cases in the UKB (AUC = 0.792) and separation in MGI (AUC = 0.686). In contrast, the best T2D model did not discriminate between T1D and T2D cases (AUC = 0.527). Our analysis suggests that a T1D PRS model based on independent single nucleotide polymorphisms may help differentiate between T1D, T2D, and controls in individuals of European genetic ancestry.

APOE alleles are associated with sex-specific structural differences in brain regions affected in Alzheimer's disease and related dementia.

Alzheimer's disease is marked by intracellular tau aggregates in the medial temporal lobe (MTL) and extracellular amyloid aggregates in the default network (DN). Here, we examined codependent structural variations between the MTL's most vulnerable structure, the hippocampus (HC), and the DN at subregion resolution in individuals with Alzheimer's disease and related dementia (ADRD). By leveraging the power of the approximately 40,000 participants of the UK Biobank cohort, we assessed impacts from the protective APOE ɛ2 and the deleterious APOE ɛ4 Alzheimer's disease alleles on these structural relationships. We demonstrate ɛ2 and ɛ4 genotype effects on the inter-individual expression of HC-DN co-variation structural patterns at the population level. Across these HC-DN signatures, recurrent deviations in the CA1, CA2/3, molecular layer, fornix's fimbria, and their cortical partners related to ADRD risk. Analyses of the rich phenotypic profiles in the UK Biobank cohort further revealed male-specific HC-DN associations with air pollution and female-specific associations with cardiovascular traits. We also showed that APOE ɛ2/2 interacts preferentially with HC-DN co-variation patterns in estimating social lifestyle in males and physical activity in females. Our structural, genetic, and phenotypic analyses in this large epidemiological cohort reinvigorate the often-neglected interplay between APOE ɛ2 dosage and sex and link APOE alleles to inter-individual brain structural differences indicative of ADRD familial risk.

Lipid-lowering drug targets and Parkinson's disease: A sex-specific Mendelian randomization study.

Parkinson's disease (PD) affects millions of individuals worldwide, and it is the second most common late-onset neurodegenerative disorder. There is no cure and current treatments only alleviate symptoms. Modifiable risk factors have been explored as possible options for decreasing risk or developing drug targets to treat PD, including low-density lipoprotein cholesterol (LDL-C). There is evidence of sex differences for cholesterol levels as well as for PD risk. Genetic datasets of increasing size are permitting association analyses with increased power, including sex-stratified analyses. These association results empower Mendelian randomization (MR) studies, which, given certain assumptions, test whether there is a causal relationship between the risk factor and the outcome using genetic instruments. Sex-specific causal inference approaches could highlight sex-specific effects that may otherwise be masked by sex-agnostic approaches. We conducted a sex-specific two-sample cis-MR analysis based on genetic variants in LDL-C target encoding genes to assess the impact of lipid-lowering drug targets on PD risk. To complement the cis-MR analysis, we also conducted a sex-specific standard MR analysis (using genome-wide independent variants). We did not find evidence of a causal relationship between LDL-C levels and PD risk in females [OR (95% CI) = 1.01 (0.60, 1.69), IVW random-effects] or males [OR (95% CI) = 0.93 (0.55, 1.56)]. The sex-specific standard MR analysis also supported this conclusion. We encourage future work assessing sex-specific effects using causal inference techniques to better understand factors that may contribute to complex disease risk differently between the sexes.

Extent to which array genotyping and imputation with large reference panels approximate deep whole-genome sequencing.

Understanding the genetic basis of human diseases and traits is dependent on the identification and accurate genotyping of genetic variants. Deep whole-genome sequencing (WGS), the gold standard technology for SNP and indel identification and genotyping, remains very expensive for most large studies. Here, we quantify the extent to which array genotyping followed by genotype imputation can approximate WGS in studies of individuals of African, Hispanic/Latino, and European ancestry in the US and of Finnish ancestry in Finland (a population isolate). For each study, we performed genotype imputation by using the genetic variants present on the Illumina Core, OmniExpress, MEGA, and Omni 2.5M arrays with the 1000G, HRC, and TOPMed imputation reference panels. Using the Omni 2.5M array and the TOPMed panel, ≥90% of bi-allelic single-nucleotide variants (SNVs) are well imputed (r2 > 0.8) down to minor-allele frequencies (MAFs) of 0.14% in African, 0.11% in Hispanic/Latino, 0.35% in European, and 0.85% in Finnish ancestries. There was little difference in TOPMed-based imputation quality among the arrays with >700k variants. Individual-level imputation quality varied widely between and within the three US studies. Imputation quality also varied across genomic regions, producing regions where even common (MAF > 5%) variants were consistently not well imputed across ancestries. The extent to which array genotyping and imputation can approximate WGS therefore depends on reference panel, genotype array, sample ancestry, and genomic location. Imputation quality by variant or genomic region can be queried with our new tool, RsqBrowser, now deployed on the Michigan Imputation Server.

Genome-wide meta-analysis of iron status biomarkers and the effect of iron on all-cause mortality in HUNT.

Iron is essential for many biological processes, but iron levels must be tightly regulated to avoid harmful effects of both iron deficiency and overload. Here, we perform genome-wide association studies on four iron-related biomarkers (serum iron, serum ferritin, transferrin saturation, total iron-binding capacity) in the Trøndelag Health Study (HUNT), the Michigan Genomics Initiative (MGI), and the SardiNIA study, followed by their meta-analysis with publicly available summary statistics, analyzing up to 257,953 individuals. We identify 123 genetic loci associated with iron traits. Among 19 novel protein-altering variants, we observe a rare missense variant (rs367731784) in HUNT, which suggests a role for DNAJC13 in transferrin recycling. We further validate recently published results using genetic risk scores for each biomarker in HUNT (6% variance in serum iron explained) and present linear and non-linear Mendelian randomization analyses of the traits on all-cause mortality. We find evidence of a harmful effect of increased serum iron and transferrin saturation in linear analyses that estimate population-averaged effects. However, there was weak evidence of a protective effect of increasing serum iron at the very low end of its distribution. Our findings contribute to our understanding of the genes affecting iron status and its consequences on human health.

Rare coding variants in ten genes confer substantial risk for schizophrenia.

Rare coding variation has historically provided the most direct connections between gene function and disease pathogenesis. By meta-analysing the whole exomes of 24,248 schizophrenia cases and 97,322 controls, we implicate ultra-rare coding variants (URVs) in 10 genes as conferring substantial risk for schizophrenia (odds ratios of 3-50, P < 2.14 × 10-6) and 32 genes at a false discovery rate of <5%. These genes have the greatest expression in central nervous system neurons and have diverse molecular functions that include the formation, structure and function of the synapse. The associations of the NMDA (N-methyl-D-aspartate) receptor subunit GRIN2A and AMPA (α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid) receptor subunit GRIA3 provide support for dysfunction of the glutamatergic system as a mechanistic hypothesis in the pathogenesis of schizophrenia. We observe an overlap of rare variant risk among schizophrenia, autism spectrum disorders1, epilepsy and severe neurodevelopmental disorders2, although different mutation types are implicated in some shared genes. Most genes described here, however, are not implicated in neurodevelopment. We demonstrate that genes prioritized from common variant analyses of schizophrenia are enriched in rare variant risk3, suggesting that common and rare genetic risk factors converge at least partially on the same underlying pathogenic biological processes. Even after excluding significantly associated genes, schizophrenia cases still carry a substantial excess of URVs, which indicates that more risk genes await discovery using this approach.

GWAS of Hematuria.

BACKGROUND AND OBJECTIVES: Glomerular hematuria has varied causes but can have a genetic basis, including Alport syndrome and IgA nephropathy. DESIGN, SETTING, PARTICIPANTS, & MEASUREMENTS: We used summary statistics to identify genetic variants associated with hematuria in White British UK Biobank participants. Individuals with glomerular hematuria were enriched by excluding participants with genitourinary conditions. A strongly associated locus on chromosome 2 (COL4A4-COL4A3) was identified. The region was reimputed using the Trans-Omics for Precision Medicine Program followed by sequential rounds of regional conditional analysis, conditioning on previous genetic signals. Similarly, we applied conditional analysis to identify independent variants in the MHC region on chromosome 6 using imputed HLA haplotypes. RESULTS: In total, 16,866 hematuria cases and 391,420 controls were included. Cases had higher urinary albumin-creatinine compared with controls (women: 13.01 mg/g [8.05-21.33] versus 12.12 mg/g [7.61-19.29]; P<0.001; men: 8.85 mg/g [5.66-16.19] versus 7.52 mg/g [5.04-12.39]; P<0.001) and lower eGFR (women: 88±14 versus 90±13 ml/min per 1.72 m2; P<0.001; men: 87±15 versus 90±13 ml/min per 1.72 m2; P<0.001), supporting enrichment of glomerular hematuria. Variants at six loci (PDPN, COL4A4-COL4A3, HLA-B, SORL1, PLLP, and TGFB1) met genome-wide significance (P<5E-8). At chromosome 2, COL4A4 p.Ser969X (rs35138315; minor allele frequency=0.00035; P<7.95E-35; odds ratio, 87.3; 95% confidence interval, 47.9 to 159.0) had the most significant association, and two variants in the locus remained associated with hematuria after conditioning for this variant: COL4A3 p.Gly695Arg (rs200287952; minor allele frequency=0.00021; P<2.16E-7; odds ratio, 45.5; 95% confidence interval, 11.8 to 168.0) and a common COL4A4 intron 25 variant (not previously reported; rs58261427; minor allele frequency=0.214; P<2.00E-9; odds ratio, 1.09; 95% confidence interval, 1.06 to 1.12). Of the HLA haplotypes, HLA-B (*0801; minor allele frequency=0.14; P<4.41E-24; odds ratio, 0.84; 95% confidence interval, 0.82 to 0.88) displayed the most statistically significant association. For remaining loci, we identified three novel associations, which were replicated in the deCODE dataset for dipstick hematuria (nearest genes: PDPN, SORL1, and PLLP). CONCLUSIONS: Our study identifies six loci associated with hematuria, including independent variants in COL4A4-COL4A3 and HLA-B. Additionally, three novel loci are reported, including an association with an intronic variant in PDPN expressed in the podocyte. PODCAST: This article contains a podcast at https://www.asn-online.org/media/podcast/CJASN/2022_04_26_CJN13711021.mp3.

Genome-wide association study of cardiac troponin I in the general population.

Circulating cardiac troponin proteins are associated with structural heart disease and predict incident cardiovascular disease in the general population. However, the genetic contribution to cardiac troponin I (cTnI) concentrations and its causal effect on cardiovascular phenotypes are unclear. We combine data from two large population-based studies, the Trøndelag Health Study and the Generation Scotland Scottish Family Health Study, and perform a genome-wide association study of high-sensitivity cTnI concentrations with 48 115 individuals. We further use two-sample Mendelian randomization to investigate the causal effects of circulating cTnI on acute myocardial infarction (AMI) and heart failure (HF). We identified 12 genetic loci (8 novel) associated with cTnI concentrations. Associated protein-altering variants highlighted putative functional genes: CAND2, HABP2, ANO5, APOH, FHOD3, TNFAIP2, KLKB1 and LMAN1. Phenome-wide association tests in 1688 phecodes and 83 continuous traits in UK Biobank showed associations between a genetic risk score for cTnI and cardiac arrhythmias, metabolic and anthropometric measures. Using two-sample Mendelian randomization, we confirmed the non-causal role of cTnI in AMI (5948 cases, 355 246 controls). We found indications for a causal role of cTnI in HF (47 309 cases and 930 014 controls), but this was not supported by secondary analyses using left ventricular mass as outcome (18 257 individuals). Our findings clarify the biology underlying the heritable contribution to circulating cTnI and support cTnI as a non-causal biomarker for AMI in the general population. Using genetically informed methods for causal inference helps inform the role and value of measuring cTnI in the general population.

Human-lineage-specific genomic elements are associated with neurodegenerative disease and APOE transcript usage.

Knowledge of genomic features specific to the human lineage may provide insights into brain-related diseases. We leverage high-depth whole genome sequencing data to generate a combined annotation identifying regions simultaneously depleted for genetic variation (constrained regions) and poorly conserved across primates. We propose that these constrained, non-conserved regions (CNCRs) have been subject to human-specific purifying selection and are enriched for brain-specific elements. We find that CNCRs are depleted from protein-coding genes but enriched within lncRNAs. We demonstrate that per-SNP heritability of a range of brain-relevant phenotypes are enriched within CNCRs. We find that genes implicated in neurological diseases have high CNCR density, including APOE, highlighting an unannotated intron-3 retention event. Using human brain RNA-sequencing data, we show the intron-3-retaining transcript to be more abundant in Alzheimer's disease with more severe tau and amyloid pathological burden. Thus, we demonstrate potential association of human-lineage-specific sequences in brain development and neurological disease.