Lineage-Specific Genome Architecture Links Enhancers and Non-coding Disease Variants to Target Gene Promoters.

Long-range interactions between regulatory elements and gene promoters play key roles in transcriptional regulation. The vast majority of interactions are uncharted, constituting a major missing link in understanding genome control. Here, we use promoter capture Hi-C to identify interacting regions of 31,253 promoters in 17 human primary hematopoietic cell types. We show that promoter interactions are highly cell type specific and enriched for links between active promoters and epigenetically marked enhancers. Promoter interactomes reflect lineage relationships of the hematopoietic tree, consistent with dynamic remodeling of nuclear architecture during differentiation. Interacting regions are enriched in genetic variants linked with altered expression of genes they contact, highlighting their functional role. We exploit this rich resource to connect non-coding disease variants to putative target promoters, prioritizing thousands of disease-candidate genes and implicating disease pathways. Our results demonstrate the power of primary cell promoter interactomes to reveal insights into genomic regulatory mechanisms underlying common diseases.

Familial pseudohyperkalemia induces significantly higher levels of extracellular potassium in early storage of red cell concentrates without affecting other standard measures of quality: A case control and allele frequency study.

BACKGROUND: Familial pseudohyperkalemia (FP) is characterized by an increased rate of potassium leakage in refrigerated red cells and is associated with the minor allele of the single nucleotide polymorphism rs148211042 (R723Q) in the ABCB6 gene. The study aims were to obtain the minor allele frequencies of ABCB6 variants and to measure supernatant potassium accumulation, and other red cell storage parameters, in red cell concentrates (RCC) from carriers of variant rs148211042 under standard blood bank conditions. STUDY DESIGN: Whole blood units were collected from 6 FP individuals and 11 controls and processed into RCC in additive solution. RCC were sampled and tested over cold storage for full blood count, extracellular potassium, glucose, lactate, microvesicle release, deformability, hemolysis, pH, adenosine triphosphate, and 2,3-diphosphoglycerate. RESULTS: Screening of genotyped cohorts identified that variant rs148211042 is present in 1 in 394 British citizens of European ancestry. FP RCC had significantly higher supernatant potassium at all time points from day 3 onwards (p < .001) and higher mean cell volume (p = .032) than controls. The initial rate of potassium release was higher in FP RCC; supernatant potassium reached 46.0 (23.8-57.6) mmol/L (mean [range]) by day 5, increasing to 68.9 (58.8-73.7) mmol/L by day 35. Other quality parameters were not significantly different between FP RCC and controls. CONCLUSION: These data suggest that if a blood donor has FP, reducing the RCC shelf-life to 5 days may be insufficient to reduce the risk of hyperkalemia in clinical scenarios such as neonatal large volume transfusion.

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⁻8) including splice variants in LPA that lowered plasma lipoprotein(a) levels (P = 1.5×10⁻¹¹7). 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⁻4), 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.

Identification of small exonic CNV from whole-exome sequence data and application to autism spectrum disorder.

Copy number variation (CNV) is an important determinant of human diversity and plays important roles in susceptibility to disease. Most studies of CNV carried out to date have made use of chromosome microarray and have had a lower size limit for detection of about 30 kilobases (kb). With the emergence of whole-exome sequencing studies, we asked whether such data could be used to reliably call rare exonic CNV in the size range of 1-30 kilobases (kb), making use of the eXome Hidden Markov Model (XHMM) program. By using both transmission information and validation by molecular methods, we confirmed that small CNV encompassing as few as three exons can be reliably called from whole-exome data. We applied this approach to an autism case-control sample (n = 811, mean per-target read depth = 161) and observed a significant increase in the burden of rare (MAF ≤1%) 1-30 kb CNV, 1-30 kb deletions, and 1-10 kb deletions in ASD. CNV in the 1-30 kb range frequently hit just a single gene, and we were therefore able to carry out enrichment and pathway analyses, where we observed enrichment for disruption of genes in cytoskeletal and autophagy pathways in ASD. In summary, our results showed that XHMM provided an effective means to assess small exonic CNV from whole-exome data, indicated that rare 1-30 kb exonic deletions could contribute to risk in up to 7% of individuals with ASD, and implicated a candidate pathway in developmental delay syndromes.

Intracranial aneurysm risk locus 5q23.2 is associated with elevated systolic blood pressure.

Although genome-wide association studies (GWAS) have identified hundreds of complex trait loci, the pathomechanisms of most remain elusive. Studying the genetics of risk factors predisposing to disease is an attractive approach to identify targets for functional studies. Intracranial aneurysms (IA) are rupture-prone pouches at cerebral artery branching sites. IA is a complex disease for which GWAS have identified five loci with strong association and a further 14 loci with suggestive association. To decipher potential underlying disease mechanisms, we tested whether there are IA loci that convey their effect through elevating blood pressure (BP), a strong risk factor of IA. We performed a meta-analysis of four population-based Finnish cohorts (n(FIN)  =  11 266) not selected for IA, to assess the association of previously identified IA candidate loci (n  =  19) with BP. We defined systolic BP (SBP), diastolic BP, mean arterial pressure, and pulse pressure as quantitative outcome variables. The most significant result was further tested for association in the ICBP-GWAS cohort of 200 000 individuals. We found that the suggestive IA locus at 5q23.2 in PRDM6 was significantly associated with SBP in individuals of European descent (p(FIN)  =  3.01E-05, p(ICBP-GWAS)  =  0.0007, p(ALL)  =  8.13E-07). The risk allele of IA was associated with higher SBP. PRDM6 encodes a protein predominantly expressed in vascular smooth muscle cells. Our study connects a complex disease (IA) locus with a common risk factor for the disease (SBP). We hypothesize that common variants in PRDM6 can contribute to altered vascular wall structure, hence increasing SBP and predisposing to IA. True positive associations often fail to reach genome-wide significance in GWAS. Our findings show that analysis of traditional risk factors as intermediate phenotypes is an effective tool for deciphering hidden heritability. Further, we demonstrate that common disease loci identified in a population isolate may bear wider significance.

Phenotype mining in CNV carriers from a population cohort.

Phenotype mining is a novel approach for elucidating the genetic basis of complex phenotypic variation. It involves a search of rich phenotype databases for measures correlated with genetic variation, as identified in genome-wide genotyping or sequencing studies. An initial implementation of phenotype mining in a prospective unselected population cohort, the Northern Finland 1966 Birth Cohort (NFBC1966), identifies neurodevelopment-related traits-intellectual deficits, poor school performance and hearing abnormalities-which are more frequent among individuals with large (>500 kb) deletions than among other cohort members. Observation of extensive shared single nucleotide polymorphism haplotypes around deletions suggests an opportunity to expand phenotype mining from cohort samples to the populations from which they derive.

Platelet surface receptor glycoprotein VI-dimer is overexpressed in stroke: The Glycoprotein VI in Stroke (GYPSIE) study results.

OBJECTIVES: Platelet activation underpins thrombus formation in ischemic stroke. The active, dimeric form of platelet receptor glycoprotein (GP) VI plays key roles by binding platelet ligands collagen and fibrin, leading to platelet activation. We investigated whether patients presenting with stroke expressed more GPVI on their platelet surface and had more active circulating platelets as measured by platelet P-selectin exposure. METHODS: 129 ischemic or hemorrhagic stroke patients were recruited within 8h of symptom onset. Whole blood was analyzed for platelet-surface expression of total GPVI, GPVI-dimer, and P-selectin by flow cytometry at admission and day-90 post-stroke. Results were compared against a healthy control population (n = 301). RESULTS: The platelets of stroke patients expressed significantly higher total GPVI and GPVI-dimer (P<0.0001) as well as demonstrating higher resting P-selectin exposure (P<0.0001), a measure of platelet activity, compared to the control group, suggesting increased circulating platelet activation. GPVI-dimer expression was strongly correlated circulating platelet activation [r2 = 0.88, P<0.0001] in stroke patients. Furthermore, higher platelet surface GPVI expression was associated with increased stroke severity at admission. At day-90 post-stroke, GPVI-dimer expression and was further raised compared to the level at admission (P<0.0001) despite anti-thrombotic therapy. All ischemic stroke subtypes and hemorrhagic strokes expressed significantly higher GPVI-dimer compared to controls (P<0.0001). CONCLUSIONS: Stroke patients express more GPVI-dimer on their platelet surface at presentation, lasting at least until day-90 post-stroke. Small molecule GPVI-dimer inhibitors are currently in development and the results of this study validate that GPVI-dimer as an anti-thrombotic target in ischemic stroke.

Transcriptional, epigenetic and metabolic signatures in cardiometabolic syndrome defined by extreme phenotypes.

BACKGROUND: This work is aimed at improving the understanding of cardiometabolic syndrome pathophysiology and its relationship with thrombosis by generating a multi-omic disease signature. METHODS/RESULTS: We combined classic plasma biochemistry and plasma biomarkers with the transcriptional and epigenetic characterisation of cell types involved in thrombosis, obtained from two extreme phenotype groups (morbidly obese and lipodystrophy) and lean individuals to identify the molecular mechanisms at play, highlighting patterns of abnormal activation in innate immune phagocytic cells. Our analyses showed that extreme phenotype groups could be distinguished from lean individuals, and from each other, across all data layers. The characterisation of the same obese group, 6 months after bariatric surgery, revealed the loss of the abnormal activation of innate immune cells previously observed. However, rather than reverting to the gene expression landscape of lean individuals, this occurred via the establishment of novel gene expression landscapes. NETosis and its control mechanisms emerge amongst the pathways that show an improvement after surgical intervention. CONCLUSIONS: We showed that the morbidly obese and lipodystrophy groups, despite some differences, shared a common cardiometabolic syndrome signature. We also showed that this could be used to discriminate, amongst the normal population, those individuals with a higher likelihood of presenting with the disease, even when not displaying the classic features.

Linkage and linkage disequilibrium scan for autism loci in an extended pedigree from Finland.

Population isolates, such as Finland, have proved beneficial in mapping rare causative genetic variants due to a limited number of founders resulting in reduced genetic heterogeneity and extensive linkage disequilibrium (LD). We have here used this special opportunity to identify rare alleles in autism by genealogically tracing 20 autism families into one extended pedigree with verified genealogical links reaching back to the 17th century. In this unique pedigree, we performed a dense microsatellite marker genome-wide scan of linkage and LD and followed initial findings with extensive fine-mapping. We identified a putative autism susceptibility locus at 19p13.3 and obtained further evidence for previously identified loci at 1q23 and 15q11-q13. Most promising candidate genes were TLE2 and TLE6 clustered at 19p13 and ATP1A2 at 1q23.

The genome-wide patterns of variation expose significant substructure in a founder population.

Although high-density SNP genotyping platforms generate a momentum for detailed genome-wide association (GWA) studies, an offshoot is a new insight into population genetics. Here, we present an example in one of the best-known founder populations by scrutinizing ten distinct Finnish early- and late-settlement subpopulations. By determining genetic distances, homozygosity, and patterns of linkage disequilibrium, we demonstrate that population substructure, and even individual ancestry, is detectable at a very high resolution and supports the concept of multiple historical bottlenecks resulting from consecutive founder effects. Given that genetic studies are currently aiming at identifying smaller and smaller genetic effects, recognizing and controlling for population substructure even at this fine level becomes imperative to avoid confounding and spurious associations. This study provides an example of the power of GWA data sets to demonstrate stratification caused by population history even within a seemingly homogeneous population, like the Finns. Further, the results provide interesting lessons concerning the impact of population history on the genome landscape of humans, as well as approaches to identify rare variants enriched in these subpopulations.

Allelic variants in HTR3C show association with autism.

Autism spectrum disorders (ASDs) are severe neurodevelopmental disorders with a strong genetic component. Only a few predisposing genes have been identified so far. We have previously performed a genome-wide linkage screen for ASDs in Finnish families where the most significant linkage peak was identified at 3q25-27. Here, 11 positional and functionally relevant candidate genes at 3q25-27 were tested for association with autistic disorder. Genotypes of 125 single nucleotide polymorphisms (SNPs) were determined in 97 families with at least one individual affected with autistic disorder. The most significant association was observed using two non-synonymous SNPs in HTR3C, rs6766410 and rs6807362, both resulting in P = 0.0012 in family-based association analysis. In addition, the haplotype C-C corresponding to amino acids N163-A405 was overtransmitted to affected individuals (P = 0.006). Sequencing revealed no other variants in the coding region or splice sites of HTR3C. Based on the association analysis results in a previously identified linkage region, we propose that HTR3C represents a novel candidate locus for ASDs and should be tested in other populations.

Identification of common genetic risk variants for autism spectrum disorder.

Autism spectrum disorder (ASD) is a highly heritable and heterogeneous group of neurodevelopmental phenotypes diagnosed in more than 1% of children. Common genetic variants contribute substantially to ASD susceptibility, but to date no individual variants have been robustly associated with ASD. With a marked sample-size increase from a unique Danish population resource, we report a genome-wide association meta-analysis of 18,381 individuals with ASD and 27,969 controls that identified five genome-wide-significant loci. Leveraging GWAS results from three phenotypes with significantly overlapping genetic architectures (schizophrenia, major depression, and educational attainment), we identified seven additional loci shared with other traits at equally strict significance levels. Dissecting the polygenic architecture, we found both quantitative and qualitative polygenic heterogeneity across ASD subtypes. These results highlight biological insights, particularly relating to neuronal function and corticogenesis, and establish that GWAS performed at scale will be much more productive in the near term in ASD.

No association between common variants in glyoxalase 1 and autism spectrum disorders.

The autism spectrum disorders (ASDs) are complex diseases with a strong genetic component. Numerous candidate gene studies have tested association between various functional and positional candidate genes and autism, but no common variation predisposing for autism has been identified to date. It has been previously proposed, that glyoxalase 1 (GLO1) might be involved in the pathogenesis of autism as GLO1 protein polarity was significantly changed in the brains of autism patients compared to controls. GLO1 harbors a functional polymorphism that affects the polarity and the enzymatic activity of the protein. In the same study, this polymorphism showed a suggestive association to autism. To investigate whether common variants in GLO1 predispose to autism in the Finnish population, we have genotyped six polymorphisms in GLO1 in families with more than 230 individuals affected with ASDs and carried out both linkage and association analyses. We did not observe significant linkage or association between any SNP and ASDs. Therefore, we suggest that common variants in GLO1 are not significant susceptibility factors for ASDs in the Finnish population.

Increased DNA methylation variability in type 1 diabetes across three immune effector cell types.

The incidence of type 1 diabetes (T1D) has substantially increased over the past decade, suggesting a role for non-genetic factors such as epigenetic mechanisms in disease development. Here we present an epigenome-wide association study across 406,365 CpGs in 52 monozygotic twin pairs discordant for T1D in three immune effector cell types. We observe a substantial enrichment of differentially variable CpG positions (DVPs) in T1D twins when compared with their healthy co-twins and when compared with healthy, unrelated individuals. These T1D-associated DVPs are found to be temporally stable and enriched at gene regulatory elements. Integration with cell type-specific gene regulatory circuits highlight pathways involved in immune cell metabolism and the cell cycle, including mTOR signalling. Evidence from cord blood of newborns who progress to overt T1D suggests that the DVPs likely emerge after birth. Our findings, based on 772 methylomes, implicate epigenetic changes that could contribute to disease pathogenesis in T1D.

Family-based association study of DYX1C1 variants in autism.

DYX1C1: was recently identified as a candidate gene for developmental dyslexia, which is characterized by an unexpected difficulty in learning to read and write despite adequate intelligence, motivation, and education. It will be important to clarify, whether the phenotype caused by DYX1C1 extends to other language-related or comorbid disorders. Impaired language development is one of the essential features in autism. Therefore, we analyzed the allelic distribution of the DYX1C1 gene by family-based association method in 100 Finnish autism families. No evidence for association was observed with any intragenic marker or with haplotypes constructed from alleles of several adjacent markers. No evidence for deviated allelic diversity was either observed: the frequency of expected dyslexia risk haplotype was comparable to its frequency in Finnish controls. Thus it seems unlikely that DYX1C1 gene would be involved in the genetic etiology of autism in Finnish patients.

Analysis of four neuroligin genes as candidates for autism.

Neuroligins are cell-adhesion molecules located at the postsynaptic side of the synapse. Neuroligins interact with beta-neurexins and this interaction is involved in the formation of functional synapses. Mutations in two X-linked neuroligin genes, NLGN3 and NLGN4, have recently been implicated in pathogenesis of autism. The neuroligin gene family consists of five members (NLGN1 at 3q26, NLGN2 at 17p13, NLGN3 at Xq13, NLGN4 at Xp22, and NLGN4Y at Yq11), of which NLGN1 and NLGN3 are located within the best loci observed in our previous genome-wide scan for autism in the Finnish sample. Here, we report a detailed molecular genetic analysis of NLGN1, NLGN3, NLGN4, and NLNG4Y in the Finnish autism sample. Mutation analysis of 30 probands selected from families producing linkage evidence for Xq13 and/or 3q26 loci revealed several polymorphisms, but none of these seemed to be functional. Family-based association analysis in 100 families with autism spectrum disorders yielded only modest associations at NLGN1 (rs1488545, P=0.002), NLGN3 (DXS7132, P=0.014), and NLGN4 (DXS996, P=0.031). We conclude that neuroligin mutations most probably represent rare causes of autism and that it is unlikely that the allelic variants in these genes would be major risk factors for autism.

A Pooled Genome-Wide Association Study of Asperger Syndrome.

Asperger Syndrome (AS) is a neurodevelopmental condition characterized by impairments in social interaction and communication, alongside the presence of unusually repetitive, restricted interests and stereotyped behaviour. Individuals with AS have no delay in cognitive and language development. It is a subset of Autism Spectrum Conditions (ASC), which are highly heritable and has a population prevalence of approximately 1%. Few studies have investigated the genetic basis of AS. To address this gap in the literature, we performed a genome-wide pooled DNA association study to identify candidate loci in 612 individuals (294 cases and 318 controls) of Caucasian ancestry, using the Affymetrix GeneChip Human Mapping version 6.0 array. We identified 11 SNPs that had a p-value below 1x10-5. These SNPs were independently genotyped in the same sample. Three of the SNPs (rs1268055, rs7785891 and rs2782448) were nominally significant, though none remained significant after Bonferroni correction. Two of our top three SNPs (rs7785891 and rs2782448) lie in loci previously implicated in ASC. However, investigation of the three SNPs in the ASC genome-wide association dataset from the Psychiatric Genomics Consortium indicated that these three SNPs were not significantly associated with ASC. The effect sizes of the variants were modest, indicating that our study was not sufficiently powered to identify causal variants with precision.

A genome wide association study of mathematical ability reveals an association at chromosome 3q29, a locus associated with autism and learning difficulties: a preliminary study.

Mathematical ability is heritable, but few studies have directly investigated its molecular genetic basis. Here we aimed to identify specific genetic contributions to variation in mathematical ability. We carried out a genome wide association scan using pooled DNA in two groups of U.K. samples, based on end of secondary/high school national academic exam achievement: high (n = 419) versus low (n = 183) mathematical ability while controlling for their verbal ability. Significant differences in allele frequencies between these groups were searched for in 906,600 SNPs using the Affymetrix GeneChip Human Mapping version 6.0 array. After meeting a threshold of p<1.5×10(-5), 12 SNPs from the pooled association analysis were individually genotyped in 542 of the participants and analyzed to validate the initial associations (lowest p-value 1.14 ×10(-6)). In this analysis, one of the SNPs (rs789859) showed significant association after Bonferroni correction, and four (rs10873824, rs4144887, rs12130910 rs2809115) were nominally significant (lowest p-value 3.278 × 10(-4)). Three of the SNPs of interest are located within, or near to, known genes (FAM43A, SFT2D1, C14orf64). The SNP that showed the strongest association, rs789859, is located in a region on chromosome 3q29 that has been previously linked to learning difficulties and autism. rs789859 lies 1.3 kbp downstream of LSG1, and 700 bp upstream of FAM43A, mapping within the potential promoter/regulatory region of the latter. To our knowledge, this is only the second study to investigate the association of genetic variants with mathematical ability, and it highlights a number of interesting markers for future study.

A framework for the interpretation of de novo mutation in human disease.

Spontaneously arising (de novo) mutations have an important role in medical genetics. For diseases with extensive locus heterogeneity, such as autism spectrum disorders (ASDs), the signal from de novo mutations is distributed across many genes, making it difficult to distinguish disease-relevant mutations from background variation. Here we provide a statistical framework for the analysis of excesses in de novo mutation per gene and gene set by calibrating a model of de novo mutation. We applied this framework to de novo mutations collected from 1,078 ASD family trios, and, whereas we affirmed a significant role for loss-of-function mutations, we found no excess of de novo loss-of-function mutations in cases with IQ above 100, suggesting that the role of de novo mutations in ASDs might reside in fundamental neurodevelopmental processes. We also used our model to identify ∼1,000 genes that are significantly lacking in functional coding variation in non-ASD samples and are enriched for de novo loss-of-function mutations identified in ASD cases.

A genome-wide analysis of populations from European Russia reveals a new pole of genetic diversity in northern Europe.

Several studies examined the fine-scale structure of human genetic variation in Europe. However, the European sets analyzed represent mainly northern, western, central, and southern Europe. Here, we report an analysis of approximately 166,000 single nucleotide polymorphisms in populations from eastern (northeastern) Europe: four Russian populations from European Russia, and three populations from the northernmost Finno-Ugric ethnicities (Veps and two contrast groups of Komi people). These were compared with several reference European samples, including Finns, Estonians, Latvians, Poles, Czechs, Germans, and Italians. The results obtained demonstrated genetic heterogeneity of populations living in the region studied. Russians from the central part of European Russia (Tver, Murom, and Kursk) exhibited similarities with populations from central-eastern Europe, and were distant from Russian sample from the northern Russia (Mezen district, Archangelsk region). Komi samples, especially Izhemski Komi, were significantly different from all other populations studied. These can be considered as a second pole of genetic diversity in northern Europe (in addition to the pole, occupied by Finns), as they had a distinct ancestry component. Russians from Mezen and the Finnic-speaking Veps were positioned between the two poles, but differed from each other in the proportions of Komi and Finnic ancestries. In general, our data provides a more complete genetic map of Europe accounting for the diversity in its most eastern (northeastern) populations.