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 genera

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

The genetics of blood pressure regulation and its target organs from association studies in 342,415 individuals

To dissect the genetic architecture of blood pressure and assess effects on target-organ damage, we analyzed 128,272 SNPs from targeted and genome-wide arrays in 201,529 individuals of European ancestry and genotypes from an additional 140,886 individuals were used for validation. We identified 66 blood pressure loci, of which 17 were novel and 15 harbored multiple distinct association signals. The 66 index SNPs were enriched for cis-regulatory elements, particularly in vascular endothelial cells, consistent with a primary role in blood pressure control through modulation of vascular tone across multiple tissues. The 66 index SNPs combined in a risk score showed comparable effects in 64,421 individuals of non-European descent. The 66-SNP blood pressure risk score was significantly associated with target-organ damage in multiple tissues, with minor effects in the kidney. Our findings expand current knowledge of blood pressure pathways and highlight tissues beyond the classic renal system in blood pressure regulation

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, here we conduct genome-wide association meta-analyses of traits related to waist and hip circumferences in up to 224,459 individuals. We identify 49 loci (33 new) associated with waist-to-hip ratio adjusted for body mass index (BMI), and an additional 19 loci newly associated with related waist and hip circumference measures (P < 5 × 10(-8)). In total, 20 of the 49 waist-to-hip ratio adjusted for BMI loci show significant sexual dimorphism, 19 of which display 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

Domain-wide regulation of gene expression in the human genome

Transcription factor complexes bind to regulatory sequences of genes, providing a system of individual expression regulation. Targets of distinct transcription factors usually map throughout the genome, without clustering. Nevertheless, highly and weakly expressed genes do cluster in separate chromosomal domains with an average size of 80–90 genes. We therefore asked whether, besides transcription factors, an additional level of gene expression regulation exists that acts on chromosomal domains. Here we show that identical green fluorescent protein (GFP) reporter constructs integrated at 90 different chromosomal positions obtain expression levels that correspond to the activity of the domains of integration. These domains are up to 80 genes long and can exert an eightfold effect on the expression levels of integrated genes. 3D-FISH shows that active domains of integration have a more open chromatin structure than integration domains with weak activity. These results reveal a novel domain-wide regulatory mechanism that, together with transcription factors, exerts a dual control over gene transcription

Enhancer of zeste homologue 2 plays an important role in neuroblastoma cell survival independent of its histone methyltransferase activity

Neuroblastoma is predominantly characterised by chromosomal rearrangements. Next to V-Myc Avian Myelocytomatosis Viral Oncogene Neuroblastoma Derived Homolog (MYCN) amplification, chromosome 7 and 17q gains are frequently observed. We identified a neuroblastoma patient with a regional 7q36 gain, encompassing the enhancer of zeste homologue 2 (EZH2) gene. EZH2 is the histone methyltransferase of lysine 27 of histone H3 (H3K27me3) that forms the catalytic subunit of the polycomb repressive complex 2. H3K27me3 is commonly associated with the silencing of genes involved in cellular processes such as cell cycle regulation, cellular differentiation and cancer. High EZH2 expression correlated with poor prognosis and overall survival independent of MYCN amplification status. Unexpectedly, treatment of 3 EZH2-high expressing neuroblastoma cell lines (IMR32, CHP134 and NMB), with EZH2-specific inhibitors (GSK126 and EPZ6438) resulted in only a slight G1 arrest, despite maximum histone methyltransferase activity inhibition. Furthermore, colony formation in cell lines treated with the inhibitors was reduced only at concentrations much higher than necessary for complete inhibition of EZH2 histone methyltransferase activity. Knockdown of the complete protein with three independent shRNAs resulted in a strong apoptotic response and decreased cyclin D1 levels. This apoptotic response could be rescued by overexpressing EZH2ΔSET, a truncated form of wild-type EZH2 lacking the SET transactivation domain necessary for histone methyltransferase activity. Our findings suggest that high EZH2 expression, at least in neuroblastoma, has a survival function independent of its methyltransferase activity. This important finding highlights the need for studies on EZH2 beyond its methyltransferase function and the requirement for compounds that will target EZH2 as a complete protei

Identification and Functional Verification of Archaeal-Type Phosphoenolpyruvate Carboxylase, a Missing Link in Archaeal Central Carbohydrate Metabolism

Despite the fact that phosphoenolpyruvate carboxylase (PEPC) activity has been measured and in some cases even purified from some Archaea, the gene responsible for this activity has not been elucidated. Using sensitive sequence comparison methods, we detected a highly conserved, uncharacterized archaeal gene family that is distantly related to the catalytic core of the canonical PEPC. To verify the predicted function of this archaeal gene family, we cloned a representative from the hyperthermophilic acidophile Sulfolobus solfataricus and functionally produced the corresponding enzyme as a fusion with the Escherichia coli maltose-binding protein. The purified fusion protein indeed displayed highly thermostable PEPC activity. The structural and biochemical properties of the characterized archaeal-type PEPC (atPEPC) from S. solfataricus are in good agreement with previously reported biochemical analyses of other archaeal PEPC enzymes. The newly identified atPEPC, with its distinct properties, constitutes yet another example of the versatility of the enzymes of the central carbon metabolic pathways in the archaeal domain

The Three-Dimensional Structure of Human Interphase Chromosomes Is Related to the Transcriptome Map▿

The three-dimensional (3D) organization of the chromosomal fiber in the human interphase nucleus is an important but poorly understood aspect of gene regulation. Here we quantitatively analyze and compare the 3D structures of two types of genomic domains as defined by the human transcriptome map. While ridges are gene dense and show high expression levels, antiridges, on the other hand, are gene poor and carry genes that are expressed at low levels. We show that ridges are in general less condensed, more irregularly shaped, and located more closely to the nuclear center than antiridges. Six human cell lines that display different gene expression patterns and karyotypes share these structural parameters of chromatin. This shows that the chromatin structures of these two types of genomic domains are largely independent of tissue-specific variations in gene expression and differentiation state. Moreover, we show that there is remarkably little intermingling of chromatin from different parts of the same chromosome in a chromosome territory, neither from adjacent nor from distant parts. This suggests that the chromosomal fiber has a compact structure that sterically suppresses intermingling. Together, our results reveal novel general aspects of 3D chromosome architecture that are related to genome structure and function

Whole-Genome Sequencing of the World’s Oldest People

<div><p>Supercentenarians (110 years or older) are the world’s oldest people. Seventy four are alive worldwide, with twenty two in the United States. We performed whole-genome sequencing on 17 supercentenarians to explore the genetic basis underlying extreme human longevity. We found no significant evidence of enrichment for a single rare protein-altering variant or for a gene harboring different rare protein altering variants in supercentenarian compared to control genomes. We followed up on the gene most enriched for rare protein-altering variants in our cohort of supercentenarians, TSHZ3, by sequencing it in a second cohort of 99 long-lived individuals but did not find a significant enrichment. The genome of one supercentenarian had a pathogenic mutation in DSC2, known to predispose to arrhythmogenic right ventricular cardiomyopathy, which is recommended to be reported to this individual as an incidental finding according to a recent position statement by the American College of Medical Genetics and Genomics. Even with this pathogenic mutation, the proband lived to over 110 years. The entire list of rare protein-altering variants and DNA sequence of all 17 supercentenarian genomes is available as a resource to assist the discovery of the genetic basis of extreme longevity in future studies.</p></div

Pipeline to test supercentenarians for enrichment of rare protein-altering variants or genes harboring them.

<p>All female Caucasian supercentenarian genomes were annotated for protein-altering variants. (A) To test for enrichment of a single variant, we filtered against dbSNP131 and compared each remaining rare protein-altering variant against 1000G EUR. No single variant was significantly enriched. (B) To test for enrichment of a gene with rare protein-altering variants, we collapsed all variants in to their respective genes and filtered against 1000G EUR (MAF<0.015). We tested for enrichment against 34 control genomes from PGP using the RVT1 burden test or a gene-based Fisher’s Exact (for recessive model). No gene was significantly enriched for rare protein-altering variants in supercentenarians. We then Sanger validated TSHZ3 as the best candidate from our burden-test for follow-up.</p