DISTMIX: direct imputation of summary statistics for unmeasured SNPs from mixed ethnicity cohorts

Motivation: To increase the signal resolution for large-scale meta-analyses of genome-wide association studies, genotypes at unmeasured single nucleotide polymorphisms (SNPs) are commonly imputed using large multi-ethnic reference panels. However, the ever increasing size and ethnic diversity of both reference panels and cohorts makes genotype imputation computationally challenging for moderately sized computer clusters. Moreover, genotype imputation requires subject-level genetic data, which unlike summary statistics provided by virtually all studies, is not publicly available. While there are much less demanding methods which avoid the genotype imputation step by directly imputing SNP statistics, e.g. Directly Imputing summary STatistics (DIST) proposed by our group, their implicit assumptions make them applicable only to ethnically homogeneous cohorts. Results: To decrease computational and access requirements for the analysis of cosmopolitan cohorts, we propose DISTMIX, which extends DIST capabilities to the analysis of mixed ethnicity cohorts. The method uses a relevant reference panel to directly impute unmeasured SNP statistics based only on statistics at measured SNPs and estimated/user-specified ethnic proportions. Simulations show that the proposed method adequately controls the Type I error rates. The 1000 Genomes panel imputation of summary statistics from the ethnically diverse Psychiatric Genetic Consortium Schizophrenia Phase 2 suggests that, when compared to genotype imputation methods, DISTMIX offers comparable imputation accuracy for only a fraction of computational resources

Profiling allele-specific gene expression in brains from individuals with autism spectrum disorder reveals preferential minor allele usage.

One fundamental but understudied mechanism of gene regulation in disease is allele-specific expression (ASE), the preferential expression of one allele. We leveraged RNA-sequencing data from human brain to assess ASE in autism spectrum disorder (ASD). When ASE is observed in ASD, the allele with lower population frequency (minor allele) is preferentially more highly expressed than the major allele, opposite to the canonical pattern. Importantly, genes showing ASE in ASD are enriched in those downregulated in ASD postmortem brains and in genes harboring de novo mutations in ASD. Two regions, 14q32 and 15q11, containing all known orphan C/D box small nucleolar RNAs (snoRNAs), are particularly enriched in shifts to higher minor allele expression. We demonstrate that this allele shifting enhances snoRNA-targeted splicing changes in ASD-related target genes in idiopathic ASD and 15q11-q13 duplication syndrome. Together, these results implicate allelic imbalance and dysregulation of orphan C/D box snoRNAs in ASD pathogenesis

Introduction to This Special Issue on Open Design at the Intersection of Making and Manufacturing

What is ‘open design’ and who gets to say what it is? In the emerging body of literature on open design, there is a clear alignment to the values and practices of free culture and open source software and hardware. Yet this same literature includes multiple, sometimes even contradictory strands of technology practice and research. These different perspectives can be traced back to free culture advocates from the 1970s to the 1990s who formulated the ideal of the internet as inherently empowering, democratizing, and countercultural. However, more recent approaches include feminist and critical interventions into hacking and making as well as corporate strategies of “open innovation” that bring end-users and consumers into the design process. What remains today seems to fall into two schools of thought. On one hand, we have the celebratory endorsements of ‘openness’ as applied to technology and design. On the other hand, we have a continuous and expanding critique of these very ideals and questions, where that critique identifies persisting forms of racial, gender, age, and class-based exclusions, and questions about the relationship between open design, labor and power remain largely unanswered

Genome-wide association analysis and fine mapping of NT-proBNP level provide novel insight into the role of the MTHFR-CLCN6-NPPA-NPPB gene cluster

High blood concentration of the N-terminal cleavage product of the B-type natriuretic peptide (NT-proBNP) is strongly associated with cardiac dysfunction and is increasingly used for heart failure diagnosis. To identify genetic variants associated with NT-proBNP level, we performed a genome-wide association analysis in 1325 individuals from South Tyrol, Italy, and followed up the most significant results in 1746 individuals from two German population-based studies. A genome-wide significant signal in the MTHFR-CLCN6-NPPA-NPPB gene cluster was replicated, after correction for multiple testing (replication one-sided P-value = 8.4 × 10−10). A conditional regression analysis of 128 single-nucleotide polymorphisms in the region of interest identified novel variants in the CLCN6 gene as independently associated with NT-proBNP. In this locus, four haplotypes were associated with increased NT-proBNP levels (haplotype-specific combined P-values from 8.3 × 10−03 to 9.3 × 10−11). The observed increase in the NT-proBNP level was proportional to the number of haplotype copies present (i.e. dosage effect), with an increase associated with two copies that varied between 20 and 100 pg/ml across populations. The identification of novel variants in the MTHFR-CLCN6-NPPA-NPPB cluster provides new insights into the biological mechanisms of cardiac dysfunction

Genetic association study of QT interval highlights role for calcium signaling pathways in myocardial repolarization.

The QT interval, an electrocardiographic measure reflecting myocardial repolarization, is a heritable trait. QT prolongation is a risk factor for ventricular arrhythmias and sudden cardiac death (SCD) and could indicate the presence of the potentially lethal mendelian long-QT syndrome (LQTS). Using a genome-wide association and replication study in up to 100,000 individuals, we identified 35 common variant loci associated with QT interval that collectively explain ∼8-10% of QT-interval variation and highlight the importance of calcium regulation in myocardial repolarization. Rare variant analysis of 6 new QT interval-associated loci in 298 unrelated probands with LQTS identified coding variants not found in controls but of uncertain causality and therefore requiring validation. Several newly identified loci encode proteins that physically interact with other recognized repolarization proteins. Our integration of common variant association, expression and orthogonal protein-protein interaction screens provides new insights into cardiac electrophysiology and identifies new candidate genes for ventricular arrhythmias, LQTS and SCD

KCND3 potassium channel gene variant confers susceptibility to electrocardiographic early repolarization pattern.

BACKGROUNDThe presence of an early repolarization pattern (ERP) on the surface ECG is associated with risk of ventricular fibrillation and sudden cardiac death. Family studies have shown that ERP is a highly heritable trait, but molecular genetic determinants are unknown.METHODSTo identify genetic susceptibility loci for ERP, we performed a GWAS and meta-analysis in 2,181 cases and 23,641 controls of European ancestry.RESULTSWe identified a genome-wide significant (P < 5 × 10-8) locus in the potassium voltage-gated channel subfamily D member 3 (KCND3) gene that was successfully replicated in additional 1,124 cases and 12,510 controls. A subsequent joint meta-analysis of the discovery and replication cohorts identified rs1545300 as the lead SNP at the KCND3 locus (OR 0.82 per minor T allele, P = 7.7 × 10-12) but did not reveal additional loci. Colocalization analyses indicate causal effects of KCND3 gene expression levels on ERP in both cardiac left ventricle and tibial artery.CONCLUSIONSIn this study, we identified for the first time to our knowledge a genome-wide significant association of a genetic variant with ERP. Our findings of a locus in the KCND3 gene provide insights not only into the genetic determinants but also into the pathophysiological mechanism of ERP, discovering a promising candidate for functional studies.FUNDINGThis project was funded by the German Center for Cardiovascular Research (DZHK Shared Expertise SE081 - STATS). For detailed funding information per study, see the Supplemental Acknowledgments

Childhood behaviour problems show the greatest gap between DNA-based and twin heritability

For most complex traits, DNA-based heritability (‘SNP heritability’) is roughly half that of twin-based heritability. A previous report from the Twins Early Development Study suggested that this heritability gap is much greater for childhood behaviour problems than for other domains. If true, this finding is important because SNP heritability, not twin heritability, is the ceiling for genome-wide association studies. With twice the sample size as the previous report, we estimated SNP heritabilities (N up to 4653 unrelated individuals) and compared them with twin heritabilities from the same sample (N up to 4724 twin pairs) for diverse domains of childhood behaviour problems as rated by parents, teachers, and children themselves at ages 12 and 16. For 37 behaviour problem measures, the average twin heritability was 0.52, whereas the average SNP heritability was just 0.06. In contrast, results for cognitive and anthropometric traits were more typical (average twin and SNP heritabilities were 0.58 and 0.28, respectively). Future research should continue to investigate the reasons why SNP heritabilities for childhood behaviour problems are so low compared with twin estimates, and find ways to maximise SNP heritability for genome-wide association studies

Identification and Functional Characterization of G6PC2 Coding Variants Influencing Glycemic Traits Define an Effector Transcript at the G6PC2-ABCB11 Locus

Genome wide association studies (GWAS) for fasting glucose (FG) and insulin (FI) have identified common variant signals which explain 4.8% and 1.2% of trait variance, respectively. It is hypothesized that low-frequency and rare variants could contribute substantially to unexplained genetic variance. To test this, we analyzed exome-array data from up to 33,231 non-diabetic individuals of European ancestry. We found exome-wide significant (P&lt;5&times;10-7) evidence for two loci not previously highlighted by common variant GWAS: GLP1R (p.Ala316Thr, minor allele frequency (MAF)=1.5%) influencing FG levels, and URB2 (p.Glu594Val, MAF = 0.1%) influencing FI levels. Coding variant associations can highlight potential effector genes at (non-coding) GWAS signals. At the G6PC2/ABCB11 locus, we identified multiple coding variants in G6PC2 (p.Val219Leu, p.His177Tyr, and p.Tyr207Ser) influencing FG levels, conditionally independent of each other and the non-coding GWAS signal. In vitro assays demonstrate that these associated coding alleles result in reduced protein abundance via proteasomal degradation, establishing G6PC2 as an effector gene at this locus. Reconciliation of single-variant associations and functional effects was only possible when haplotype phase was considered. In contrast to earlier reports suggesting that, paradoxically, glucose-raising alleles at this locus are protective against type 2 diabetes (T2D), the p.Val219Leu G6PC2 variant displayed a modest but directionally consistent association with T2D risk. Coding variant associations for glycemic traits in GWAS signals highlight PCSK1, RREB1, and ZHX3 as likely effector transcripts. These coding variant association signals do not have a major impact on the trait variance explained, but they do provide valuable biological insights