Genome-wide meta-analysis of ascertainment and symptom structures of major depression in case-enriched and community cohorts.

BACKGROUND: Diagnostic criteria for major depressive disorder allow for heterogeneous symptom profiles but genetic analysis of major depressive symptoms has the potential to identify clinical and etiological subtypes. There are several challenges to integrating symptom data from genetically informative cohorts, such as sample size differences between clinical and community cohorts and various patterns of missing data. METHODS: We conducted genome-wide association studies of major depressive symptoms in three cohorts that were enriched for participants with a diagnosis of depression (Psychiatric Genomics Consortium, Australian Genetics of Depression Study, Generation Scotland) and three community cohorts who were not recruited on the basis of diagnosis (Avon Longitudinal Study of Parents and Children, Estonian Biobank, and UK Biobank). We fit a series of confirmatory factor models with factors that accounted for how symptom data was sampled and then compared alternative models with different symptom factors. RESULTS: The best fitting model had a distinct factor for Appetite/Weight symptoms and an additional measurement factor that accounted for the skip-structure in community cohorts (use of Depression and Anhedonia as gating symptoms). CONCLUSION: The results show the importance of assessing the directionality of symptoms (such as hypersomnia versus insomnia) and of accounting for study and measurement design when meta-analyzing genetic association data.

Gene set enrichment analysis of pathophysiological pathways highlights oxidative stress in psychosis.

Polygenic risk prediction remains an important aim of genetic association studies. Currently, the predictive power of schizophrenia polygenic risk scores (PRSs) is not large enough to allow highly accurate discrimination between cases and controls and thus is not adequate for clinical integration. Since PRSs are rarely used to reveal biological functions or to validate candidate pathways, to fill this gap, we investigated whether their predictive ability could be improved by building genome-wide (GW-PRSs) and pathway-specific PRSs, using distance- or expression quantitative trait loci (eQTLs)- based mapping between genetic variants and genes. We focused on five pathways (glutamate, oxidative stress, GABA/interneurons, neuroimmune/neuroinflammation and myelin) which belong to a critical hub of schizophrenia pathophysiology, centred on redox dysregulation/oxidative stress. Analyses were first performed in the Lausanne Treatment and Early Intervention in Psychosis Program (TIPP) study (n = 340, cases/controls: 208/132), a sample of first-episode of psychosis patients and matched controls, and then validated in an independent study, the epidemiological and longitudinal intervention program of First-Episode Psychosis in Cantabria (PAFIP) (n = 352, 224/128). Our results highlighted two main findings. First, GW-PRSs for schizophrenia were significantly associated with early psychosis status. Second, oxidative stress was the only significantly associated pathway that showed an enrichment in both the TIPP (p = 0.03) and PAFIP samples (p = 0.002), and exclusively when gene-variant linking was done using eQTLs. The results suggest that the predictive accuracy of polygenic risk scores could be improved with the inclusion of information from functional annotations, and through a focus on specific pathways, emphasizing the need to build and study functionally informed risk scores.

Uncovering the genetic architecture of broad antisocial behavior through a genome-wide association study meta-analysis.

Despite the substantial heritability of antisocial behavior (ASB), specific genetic variants robustly associated with the trait have not been identified. The present study by the Broad Antisocial Behavior Consortium (BroadABC) meta-analyzed data from 28 discovery samples (N = 85,359) and five independent replication samples (N = 8058) with genotypic data and broad measures of ASB. We identified the first significant genetic associations with broad ASB, involving common intronic variants in the forkhead box protein P2 (FOXP2) gene (lead SNP rs12536335, p = 6.32 × 10-10). Furthermore, we observed intronic variation in Foxp2 and one of its targets (Cntnap2) distinguishing a mouse model of pathological aggression (BALB/cJ strain) from controls (BALB/cByJ strain). Polygenic risk score (PRS) analyses in independent samples revealed that the genetic risk for ASB was associated with several antisocial outcomes across the lifespan, including diagnosis of conduct disorder, official criminal convictions, and trajectories of antisocial development. We found substantial genetic correlations of ASB with mental health (depression rg = 0.63, insomnia rg = 0.47), physical health (overweight rg = 0.19, waist-to-hip ratio rg = 0.32), smoking (rg = 0.54), cognitive ability (intelligence rg = -0.40), educational attainment (years of schooling rg = -0.46) and reproductive traits (age at first birth rg = -0.58, father's age at death rg = -0.54). Our findings provide a starting point toward identifying critical biosocial risk mechanisms for the development of ASB.

Major depression subtypes are differentially associated with migraine subtype, prevalence and severity.

OBJECTIVE: Migraine and major depressive disorder show a high rate of comorbidity, but little is known about the associations between the subtypes of major depressive disorder and migraine. In this cross-sectional study we aimed at investigating a) the lifetime associations between the atypical, melancholic, combined and unspecified subtype of major depressive disorder and migraine with and without aura and b) the associations between major depressive disorder and its subtypes and the severity of migraine. METHODS: A total of 446 subjects with migraine (migraine without aura: n = 294; migraine with aura: n = 152) and 2511 controls from the population-based CoLaus/PsyCoLaus study, Switzerland, were included. Associations between major depressive disorder subtypes and migraine characteristics were tested using binary logistic or linear regression. RESULTS: Melancholic, combined and unspecified major depressive disorder were associated with increased frequency of migraine with aura, whereas only melancholic major depressive disorder was associated with increased frequency of migraine without aura. Lifetime and unspecified major depressive disorder were associated with severe migraine intensity among subjects with migraine with aura but not migraine without aura, while combined major depressive disorder was associated with higher migraine frequency independently from migraine subtype. CONCLUSION: This study suggests that melancholic but not atypical major depressive disorder is associated with migraine and migraine subtypes. Future studies exploring pathophysiological mechanisms shared between melancholic depression and migraine are warranted.

A pharmacogenetic risk score for the evaluation of major depression severity under treatment with antidepressants.

The severity of symptoms as well as efficacy of antidepressants in major depressive disorder (MDD) is modified by single nucleotide polymorphisms (SNPs) in different genes, which may contribute in an additive or synergistic fashion. We aimed to investigate depression severity in participants with MDD under treatment with antidepressants in relation to the combinatory effect of selected genetic variants combined using a genetic risk score (GRS). The sample included 150 MDD patients on regular AD therapy from the population-based Swiss PsyCoLaus cohort. We investigated 44 SNPs previously associated with antidepressant response by ranking them with regard to their association to the Center for Epidemiologic Studies Short Depression Scale (CES-D) score using random forest. The three top scoring SNPs (rs12248560, rs878567, rs17710780) were subsequently combined into an unweighted GRS, which was included in linear and logistic regression models using the CES-D score, occurrence of a major depressive episode (MDE) during follow-up and regular antidepressant treatment during the 6 months preceding follow-up assessment as outcomes. The GRS was associated with MDE occurrence (p = .02) and ln CES-D score (p = .001). The HTR1A rs878567 variant was associated with ln CES-D after adjustment for demographic and clinical variables [p = .02, lower scores for minor allele (G) carriers]. Additionally, rs12248560 (CYP2C19) CC homozygotes showed a six-fold higher likelihood of regular AD therapy at follow-up compared to minor allele homozygotes [TT; ultrarapid metabolizers (p = .03)]. Our study suggests that the cumulative consideration of pharmacogenetic risk variants more reliably reflects the impact of the genetic background on depression severity than individual SNPs.

Genetic comorbidity between major depression and cardio-metabolic traits, stratified by age at onset of major depression.

It is imperative to understand the specific and shared etiologies of major depression and cardio-metabolic disease, as both traits are frequently comorbid and each represents a major burden to society. This study examined whether there is a genetic association between major depression and cardio-metabolic traits and if this association is stratified by age at onset for major depression. Polygenic risk scores analysis and linkage disequilibrium score regression was performed to examine whether differences in shared genetic etiology exist between depression case control status (N cases = 40,940, N controls = 67,532), earlier (N = 15,844), and later onset depression (N = 15,800) with body mass index, coronary artery disease, stroke, and type 2 diabetes in 11 data sets from the Psychiatric Genomics Consortium, Generation Scotland, and UK Biobank. All cardio-metabolic polygenic risk scores were associated with depression status. Significant genetic correlations were found between depression and body mass index, coronary artery disease, and type 2 diabetes. Higher polygenic risk for body mass index, coronary artery disease, and type 2 diabetes was associated with both early and later onset depression, while higher polygenic risk for stroke was associated with later onset depression only. Significant genetic correlations were found between body mass index and later onset depression, and between coronary artery disease and both early and late onset depression. The phenotypic associations between major depression and cardio-metabolic traits may partly reflect their overlapping genetic etiology irrespective of the age depression first presents.

Discovery of novel heart rate-associated loci using the Exome Chip.

Resting heart rate is a heritable trait, and an increase in heart rate is associated with increased mortality risk. Genome-wide association study analyses have found loci associated with resting heart rate, at the time of our study these loci explained 0.9% of the variation. This study aims to discover new genetic loci associated with heart rate from Exome Chip meta-analyses.Heart rate was measured from either elecrtrocardiograms or pulse recordings. We meta-analysed heart rate association results from 104 452 European-ancestry individuals from 30 cohorts, genotyped using the Exome Chip. Twenty-four variants were selected for follow-up in an independent dataset (UK Biobank, N = 134 251). Conditional and gene-based testing was undertaken, and variants were investigated with bioinformatics methods.We discovered five novel heart rate loci, and one new independent low-frequency non-synonymous variant in an established heart rate locus (KIAA1755). Lead variants in four of the novel loci are non-synonymous variants in the genes C10orf71, DALDR3, TESK2 and SEC31B. The variant at SEC31B is significantly associated with SEC31B expression in heart and tibial nerve tissue. Further candidate genes were detected from long-range regulatory chromatin interactions in heart tissue (SCD, SLF2 and MAPK8). We observed significant enrichment in DNase I hypersensitive sites in fetal heart and lung. Moreover, enrichment was seen for the first time in human neuronal progenitor cells (derived from embryonic stem cells) and fetal muscle samples by including our novel variants.Our findings advance the knowledge of the genetic architecture of heart rate, and indicate new candidate genes for follow-up functional studies.

Psychosocial Stress Over the Lifespan, Psychological Factors, and Cardiometabolic Risk in the Community.

OBJECTIVE: The complex relationship between psychosocial stress over the lifetime, psychological factors, and cardiometabolic risk is still poorly understood. Accordingly, our aims were (1) to independently assess the associations between childhood adversity, life-event stress in remote (earlier than the last 5 years), and recent adulthood and cardiometabolic risk, and (2) to determine the role of psychological factors including personality, coping, and depression in these associations. METHODS: The sample included 2674 adults, aged 35 to 66 years, randomly selected from urban area. Participants underwent a physical examination including the assessment of obesity markers, blood pressure, and blood lipid and glucose levels. Stress during adulthood was determined using the severity scores of 52 stressful life events. Information on adverse childhood experiences and major depressive disorders was collected using semistructured interviews, whereas personality traits and coping mechanisms were evaluated through questionnaires. RESULTS: Both childhood adversity and stress in remote adulthood were associated with elevated body mass index (ß [95% confidence interval {CI}] = 0.249 [0.029 to 0.468]; 0.020 [0.006 to 0.034]), waist circumference (ß [95% CI] = 0.061 [0.024 to 0.099]; 0.08 [0.04 to 0.11]), and the global cardiometabolic risk score (ß [95% CI] = 0.278 [0.017 to 0.540]; 0.017 [0.001 to 0.033]) after adjustment for sociodemographic, lifestyle, and psychological factors. In addition, childhood adversity was associated with low high density lipoprotein levels (ß [95% CI] = -0.021 [-0.042 to 0.000]), as well as increased fat mass and systolic blood pressure levels (ß [95% CI] = 0.506 [0.165 to 0.846]; 0.952 [0.165 to 1.740]) and stress in remote adulthood with apolipoprotein B levels (ß [95% CI] = 0.607 [0.312 to 0.901]). Psychological factors did not account for these associations and were not effect modifiers. CONCLUSIONS: Our data demonstrate that psychosocial stress during childhood and remote adulthood favor adiposity and abnormal lipid metabolism.

Genome-wide gene-environment interaction in depression: A systematic evaluation of candidate genes: The childhood trauma working-group of PGC-MDD.

Gene by environment (GxE) interaction studies have investigated the influence of a number of candidate genes and variants for major depressive disorder (MDD) on the association between childhood trauma and MDD. Most of these studies are hypothesis driven and investigate only a limited number of SNPs in relevant pathways using differing methodological approaches. Here (1) we identified 27 genes and 268 SNPs previously associated with MDD or with GxE interaction in MDD and (2) analyzed their impact on GxE in MDD using a common approach in 3944 subjects of European ancestry from the Psychiatric Genomics Consortium who had completed the Childhood Trauma Questionnaire. (3) We subsequently used the genome-wide SNP data for a genome-wide case-control GxE model and GxE case-only analyses testing for an enrichment of associated SNPs. No genome-wide significant hits and no consistency among the signals of the different analytic approaches could be observed. This is the largest study for systematic GxE interaction analysis in MDD in subjects of European ancestry to date. Most of the known candidate genes/variants could not be supported. Thus, their impact on GxE interaction in MDD may be questionable. Our results underscore the need for larger samples, more extensive assessment of environmental exposures, and greater efforts to investigate new methodological approaches in GxE models for MDD.

A genetic risk score is differentially associated with migraine with and without aura.

Although a number of migraine-associated single-nucleotide polymorphisms (SNP) with small effect size have been identified, little is known about the additive impact of these variants on migraine risk, frequency and severity. We investigated to what extent a genetic risk score (GRS) based on recently published, novel migraine-associated SNPs is associated with migraine prevalence, subtypes and severity in a large population-based sample. The sample comprised 446 subjects with migraine and 2511 controls from the CoLaus/PsyCoLaus study. Fifty-four SNPs earlier associated with migraine were selected. SNPs with a low impact on migraine prevalence in our sample were excluded using random forest. We combined the remaining 21 SNPs into a GRS and analyzed the association with migraine using logistic regression models. The GRS was significantly associated with migraine (OR = 1.56, p = 0.02) and migraine without aura (MWOA) (OR = 2.01, p = 0.003), but not with migraine with aura (MWA). The GRS was not associated with migraine frequency, intensity or interference with daily activities. We show that a GRS combining multiple genetic risk variants is associated with MWOA but not MWA, suggesting a different genetic susceptibility background underlying the two forms of migraine.

New gene functions in megakaryopoiesis and platelet formation.

Platelets are the second most abundant cell type in blood and are essential for maintaining haemostasis. Their count and volume are tightly controlled within narrow physiological ranges, but there is only limited understanding of the molecular processes controlling both traits. Here we carried out a high-powered meta-analysis of genome-wide association studies (GWAS) in up to 66,867 individuals of European ancestry, followed by extensive biological and functional assessment. We identified 68 genomic loci reliably associated with platelet count and volume mapping to established and putative novel regulators of megakaryopoiesis and platelet formation. These genes show megakaryocyte-specific gene expression patterns and extensive network connectivity. Using gene silencing in Danio rerio and Drosophila melanogaster, we identified 11 of the genes as novel regulators of blood cell formation. Taken together, our findings advance understanding of novel gene functions controlling fate-determining events during megakaryopoiesis and platelet formation, providing a new example of successful translation of GWAS to function.

Identification of novel genetic Loci associated with thyroid peroxidase antibodies and clinical thyroid disease.

Autoimmune thyroid diseases (AITD) are common, affecting 2-5% of the general population. Individuals with positive thyroid peroxidase antibodies (TPOAbs) have an increased risk of autoimmune hypothyroidism (Hashimoto's thyroiditis), as well as autoimmune hyperthyroidism (Graves' disease). As the possible causative genes of TPOAbs and AITD remain largely unknown, we performed GWAS meta-analyses in 18,297 individuals for TPOAb-positivity (1769 TPOAb-positives and 16,528 TPOAb-negatives) and in 12,353 individuals for TPOAb serum levels, with replication in 8,990 individuals. Significant associations (P<5×10(-8)) were detected at TPO-rs11675434, ATXN2-rs653178, and BACH2-rs10944479 for TPOAb-positivity, and at TPO-rs11675434, MAGI3-rs1230666, and KALRN-rs2010099 for TPOAb levels. Individual and combined effects (genetic risk scores) of these variants on (subclinical) hypo- and hyperthyroidism, goiter and thyroid cancer were studied. Individuals with a high genetic risk score had, besides an increased risk of TPOAb-positivity (OR: 2.18, 95% CI 1.68-2.81, P = 8.1×10(-8)), a higher risk of increased thyroid-stimulating hormone levels (OR: 1.51, 95% CI 1.26-1.82, P = 2.9×10(-6)), as well as a decreased risk of goiter (OR: 0.77, 95% CI 0.66-0.89, P = 6.5×10(-4)). The MAGI3 and BACH2 variants were associated with an increased risk of hyperthyroidism, which was replicated in an independent cohort of patients with Graves' disease (OR: 1.37, 95% CI 1.22-1.54, P = 1.2×10(-7) and OR: 1.25, 95% CI 1.12-1.39, P = 6.2×10(-5)). The MAGI3 variant was also associated with an increased risk of hypothyroidism (OR: 1.57, 95% CI 1.18-2.10, P = 1.9×10(-3)). This first GWAS meta-analysis for TPOAbs identified five newly associated loci, three of which were also associated with clinical thyroid disease. With these markers we identified a large subgroup in the general population with a substantially increased risk of TPOAbs. The results provide insight into why individuals with thyroid autoimmunity do or do not eventually develop thyroid disease, and these markers may therefore predict which TPOAb-positives are particularly at risk of developing clinical thyroid dysfunction.

Methods for association analysis and meta-analysis of rare variants in families.

Advances in exome sequencing and the development of exome genotyping arrays are enabling explorations of association between rare coding variants and complex traits. To ensure power for these rare variant analyses, a variety of association tests that group variants by gene or functional unit have been proposed. Here, we extend these tests to family-based studies. We develop family-based burden tests, variable frequency threshold tests and sequence kernel association tests. Through simulations, we compare the performance of different tests. We describe situations where family-based studies provide greater power than studies of unrelated individuals to detect rare variants associated with moderate to large changes in trait values. Broadly speaking, we find that when sample sizes are limited and only a modest fraction of all trait-associated variants can be identified, family samples are more powerful. Finally, we illustrate our approach by analyzing the relationship between coding variants and levels of high-density lipoprotein (HDL) cholesterol in 11,556 individuals from the HUNT and SardiNIA studies, demonstrating association for coding variants in the APOC3, CETP, LIPC, LIPG, and LPL genes and illustrating the value of family samples, meta-analysis, and gene-level tests. Our methods are implemented in freely available C++ code.

A Genome-Wide Association Study in isolated populations reveals new genes associated to common food likings.

Food preferences are the first factor driving food choice and thus nutrition. They involve numerous different senses such as taste and olfaction as well as various other factors such as personal experiences and hedonistic aspects. Although it is clear that several of these have a genetic basis, up to now studies have focused mostly on the effects of polymorphisms of taste receptor genes. Therefore, we have carried out one of the first large scale (4611 individuals) GWAS on food likings assessed for 20 specific food likings belonging to 4 different categories (vegetables, fatty, dairy and bitter). A two-step meta-analysis using three different isolated populations from Italy for the discovery step and two populations from The Netherlands and Central Asia for replication, revealed 15 independent genome-wide significant loci (p < 5 × 10(-8)) for 12 different foods. None of the identified genes coded for either taste or olfactory receptors suggesting that genetics impacts in determining food likings in a much broader way than simple differences in taste perception. These results represent a further step in uncovering the genes that underlie liking of common foods that in the end will greatly help understanding the genetics of human nutrition in general.

A meta-analysis of thyroid-related traits reveals novel loci and gender-specific differences in the regulation of thyroid function.

Thyroid hormone is essential for normal metabolism and development, and overt abnormalities in thyroid function lead to common endocrine disorders affecting approximately 10% of individuals over their life span. In addition, even mild alterations in thyroid function are associated with weight changes, atrial fibrillation, osteoporosis, and psychiatric disorders. To identify novel variants underlying thyroid function, we performed a large meta-analysis of genome-wide association studies for serum levels of the highly heritable thyroid function markers TSH and FT4, in up to 26,420 and 17,520 euthyroid subjects, respectively. Here we report 26 independent associations, including several novel loci for TSH (PDE10A, VEGFA, IGFBP5, NFIA, SOX9, PRDM11, FGF7, INSR, ABO, MIR1179, NRG1, MBIP, ITPK1, SASH1, GLIS3) and FT4 (LHX3, FOXE1, AADAT, NETO1/FBXO15, LPCAT2/CAPNS2). Notably, only limited overlap was detected between TSH and FT4 associated signals, in spite of the feedback regulation of their circulating levels by the hypothalamic-pituitary-thyroid axis. Five of the reported loci (PDE8B, PDE10A, MAF/LOC440389, NETO1/FBXO15, and LPCAT2/CAPNS2) show strong gender-specific differences, which offer clues for the known sexual dimorphism in thyroid function and related pathologies. Importantly, the TSH-associated loci contribute not only to variation within the normal range, but also to TSH values outside the reference range, suggesting that they may be involved in thyroid dysfunction. Overall, our findings explain, respectively, 5.64% and 2.30% of total TSH and FT4 trait variance, and they improve the current knowledge of the regulation of hypothalamic-pituitary-thyroid axis function and the consequences of genetic variation for hypo- or hyperthyroidism.

The association of thyroid peroxidase antibody risk loci with susceptibility to and phenotype of Graves' disease.

BACKGROUND: Despite great progress, the genetic basis of Graves' disease (GD) remains poorly understood. Recently, a population-based genomewide association study (GWAS) identified five novel loci (ATXN2/SH2B3, MAGI3, BACH2, TPO and KALRN) as significantly associated with the presence of thyroid peroxidase autoantibodies (TPOAbs), whereas several other loci showed suggestive association. METHODS: In this study, we investigated 16 single nucleotide polymorphisms (SNPs) associated with TPOAbs for the association with susceptibility to and phenotype of GD in a cohort of 647 patients with GD and 769 controls from a Polish Caucasian population. RESULTS: SNPs within/near HCP5 (rs3094228, P = 1·6 × 10(-12) , OR = 1·88), MAGI3 (rs1230666, P = 1·9 × 10(-5) , OR = 1·51) and ATXN2/SH2B3 (rs653178, P = 0·0015, OR = 1·28) loci were significantly associated with susceptibility to GD. Allele frequencies differed significantly in subgroups of patients with GD stratified by age of GD onset for HCP5 (P = 0·0014, OR = 1·50) and showed a suggestive difference for MAGI3 (P = 0·0035, OR = 1·50) SNPs. Although rs11675434 located near TPO showed no association with GD susceptibility, it was significantly associated with the presence of clinically evident Graves' ophthalmopathy (GO, P = 5·2 × 10(-5) , OR = 1·64), and this effect was independent from smoking status, age of GD onset and gender. CONCLUSIONS: This is the first study showing an association of the ATXN2/SH2B3 locus with susceptibility to GD. Furthermore, we observed a novel significant association within the HLA region at a SNP located near HCP5 and confirmed the association of the MAGI3 locus with GD susceptibility. HCP5 and MAGI3 SNPs were further correlated with age of GD onset. Finally, we identified TPO as a new susceptibility locus for GO.

Genome-wide association and functional follow-up reveals new loci for kidney function.

Chronic kidney disease (CKD) is an important public health problem with a genetic component. We performed genome-wide association studies in up to 130,600 European ancestry participants overall, and stratified for key CKD risk factors. We uncovered 6 new loci in association with estimated glomerular filtration rate (eGFR), the primary clinical measure of CKD, in or near MPPED2, DDX1, SLC47A1, CDK12, CASP9, and INO80. Morpholino knockdown of mpped2 and casp9 in zebrafish embryos revealed podocyte and tubular abnormalities with altered dextran clearance, suggesting a role for these genes in renal function. By providing new insights into genes that regulate renal function, these results could further our understanding of the pathogenesis of CKD.

Common variants in UMOD associate with urinary uromodulin levels: a meta-analysis.

Uromodulin is expressed exclusively in the thick ascending limb and is the most abundant protein excreted in normal urine. Variants in UMOD, which encodes uromodulin, are associated with renal function, and urinary uromodulin levels may be a biomarker for kidney disease. However, the genetic factors regulating uromodulin excretion are unknown. We conducted a meta-analysis of urinary uromodulin levels to identify associated common genetic variants in the general population. We included 10,884 individuals of European descent from three genetic isolates and three urban cohorts. Each study measured uromodulin indexed to creatinine and conducted linear regression analysis of approximately 2.5 million single nucleotide polymorphisms using an additive model. We also tested whether variants in genes expressed in the thick ascending limb associate with uromodulin levels. rs12917707, located near UMOD and previously associated with renal function and CKD, had the strongest association with urinary uromodulin levels (P<0.001). In all cohorts, carriers of a G allele of this variant had higher uromodulin levels than noncarriers did (geometric means 10.24, 14.05, and 17.67 µg/g creatinine for zero, one, or two copies of the G allele). rs12446492 in the adjacent gene PDILT (protein disulfide isomerase-like, testis expressed) also reached genome-wide significance (P<0.001). Regarding genes expressed in the thick ascending limb, variants in KCNJ1, SORL1, and CAB39 associated with urinary uromodulin levels. These data indicate that common variants in the UMOD promoter region may influence urinary uromodulin levels. They also provide insights into uromodulin biology and the association of UMOD variants with renal function.

Integration of genome-wide association studies with biological knowledge identifies six novel genes related to kidney function.

In conducting genome-wide association studies (GWAS), analytical approaches leveraging biological information may further understanding of the pathophysiology of clinical traits. To discover novel associations with estimated glomerular filtration rate (eGFR), a measure of kidney function, we developed a strategy for integrating prior biological knowledge into the existing GWAS data for eGFR from the CKDGen Consortium. Our strategy focuses on single nucleotide polymorphism (SNPs) in genes that are connected by functional evidence, determined by literature mining and gene ontology (GO) hierarchies, to genes near previously validated eGFR associations. It then requires association thresholds consistent with multiple testing, and finally evaluates novel candidates by independent replication. Among the samples of European ancestry, we identified a genome-wide significant SNP in FBXL20 (P = 5.6 × 10(-9)) in meta-analysis of all available data, and additional SNPs at the INHBC, LRP2, PLEKHA1, SLC3A2 and SLC7A6 genes meeting multiple-testing corrected significance for replication and overall P-values of 4.5 × 10(-4)-2.2 × 10(-7). Neither the novel PLEKHA1 nor FBXL20 associations, both further supported by association with eGFR among African Americans and with transcript abundance, would have been implicated by eGFR candidate gene approaches. LRP2, encoding the megalin receptor, was identified through connection with the previously known eGFR gene DAB2 and extends understanding of the megalin system in kidney function. These findings highlight integration of existing genome-wide association data with independent biological knowledge to uncover novel candidate eGFR associations, including candidates lacking known connections to kidney-specific pathways. The strategy may also be applicable to other clinical phenotypes, although more testing will be needed to assess its potential for discovery in general.

Common variants in Mendelian kidney disease genes and their association with renal function.

Many common genetic variants identified by genome-wide association studies for complex traits map to genes previously linked to rare inherited Mendelian disorders. A systematic analysis of common single-nucleotide polymorphisms (SNPs) in genes responsible for Mendelian diseases with kidney phenotypes has not been performed. We thus developed a comprehensive database of genes for Mendelian kidney conditions and evaluated the association between common genetic variants within these genes and kidney function in the general population. Using the Online Mendelian Inheritance in Man database, we identified 731 unique disease entries related to specific renal search terms and confirmed a kidney phenotype in 218 of these entries, corresponding to mutations in 258 genes. We interrogated common SNPs (minor allele frequency >5%) within these genes for association with the estimated GFR in 74,354 European-ancestry participants from the CKDGen Consortium. However, the top four candidate SNPs (rs6433115 at LRP2, rs1050700 at TSC1, rs249942 at PALB2, and rs9827843 at ROBO2) did not achieve significance in a stage 2 meta-analysis performed in 56,246 additional independent individuals, indicating that these common SNPs are not associated with estimated GFR. The effect of less common or rare variants in these genes on kidney function in the general population and disease-specific cohorts requires further research.