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.

New loci associated with birth weight identify genetic links between intrauterine growth and adult height and metabolism

Birth weight within the normal range is associated with a variety of adult-onset diseases, but the mechanisms behind these associations are poorly understood. Previous genome-wide association studies of birth weight identified a variant in the ADCY5 gene associated both with birth weight and type 2 diabetes and a second variant, near CCNL1, with no obvious link to adult traits. In an expanded genome-wide association meta-analysis and follow-up study of birth weight (of up to 69,308 individuals of European descent from 43 studies), we have now extended the number of loci associated at genome-wide significance to 7, accounting for a similar proportion of variance as maternal smoking. Five of the loci are known to be associated with other phenotypes: ADCY5 and CDKAL1 with type 2 diabetes, ADRB1 with adult blood pressure and HMGA2 and LCORL with adult height. Our findings highlight genetic links between fetal growth and postnatal growth and metabolism.

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.

Hundreds of variants clustered in genomic loci and biological pathways affect human height

Most common human traits and diseases have a polygenic pattern of inheritance: DNA sequence variants at many genetic loci influence the phenotype. Genome-wide association (GWA) studies have identified more than 600 variants associated with human traits, but these typically explain small fractions of phenotypic variation, raising questions about the use of further studies. Here, using 183,727 individuals, we show that hundreds of genetic variants, in at least 180 loci, influence adult height, a highly heritable and classic polygenic trait. The large number of loci reveals patterns with important implications for genetic studies of common human diseases and traits. First, the 180 loci are not random, but instead are enriched for genes that are connected in biological pathways (P = 0.016) and that underlie skeletal growth defects (P < 0.001). Second, the likely causal gene is often located near the most strongly associated variant: in 13 of 21 loci containing a known skeletal growth gene, that gene was closest to the associated variant. Third, at least 19 loci have multiple independently associated variants, suggesting that allelic heterogeneity is a frequent feature of polygenic traits, that comprehensive explorations of already-discovered loci should discover additional variants and that an appreciable fraction of associated loci may have been identified. Fourth, associated variants are enriched for likely functional effects on genes, being over-represented among variants that alter amino-acid structure of proteins and expression levels of nearby genes. Our data explain approximately 10% of the phenotypic variation in height, and we estimate that unidentified common variants of similar effect sizes would increase this figure to approximately 16% of phenotypic variation (approximately 20% of heritable variation). Although additional approaches are needed to dissect the genetic architecture of polygenic human traits fully, our findings indicate that GWA studies can identify large numbers of loci that implicate biologically relevant genes and pathways.

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.