Collaborative genome-wide association analysis supports a role for ANK3 and CACNA1C in bipolar disorder

To identify susceptibility loci for bipolar disorder, we tested 1.8 million variants in 4,387 cases and 6,209 controls and identified a region of strong association (rs10994336, P = 9.1 x 10(-9)) in ANK3 (ankyrin G). We also found further support for the previously reported CACNA1C (alpha 1C subunit of the L-type voltage-gated calcium channel; combined P = 7.0 x 10(-8), rs1006737). Our results suggest that ion channelopathies may be involved in the pathogenesis of bipolar disorder.     

Susceptibility loci for intracranial aneurysm in European and Japanese populations

Stroke is the world's third leading cause of death. One cause of stroke, intracranial aneurysm, affects approximately 2% of the population and accounts for 500,000 hemorrhagic strokes annually in mid-life (median age 50), most often resulting in death or severe neurological impairment. The pathogenesis of intracranial aneurysm is unknown, and because catastrophic hemorrhage is commonly the first sign of disease, early identification is essential. We carried out a multistage genome-wide association study (GWAS) of Finnish, Dutch and Japanese cohorts including over 2,100 intracranial aneurysm cases and 8,000 controls. Genome-wide genotyping of the European cohorts and replication studies in the Japanese cohort identified common SNPs on chromosomes 2q, 8q and 9p that show significant association with intracranial aneurysm with odds ratios 1.24-1.36. The loci on 2q and 8q are new, whereas the 9p locus was previously found to be associated with arterial diseases, including intracranial aneurysm. Associated SNPs on 8q likely act via SOX17, which is required for formation and maintenance of endothelial cells, suggesting a role in development and repair of the vasculature; CDKN2A at 9p may have a similar role. These findings have implications for the pathophysiology, diagnosis and therapy of intracranial aneurysm.     

A robust statistical method for case-control association testing with copy number variation

Copy number variation (CNV) is pervasive in the human genome and can play a causal role in genetic diseases. The functional impact of CNV cannot be fully captured through linkage disequilibrium with SNPs. These observations motivate the development of statistical methods for performing direct CNV association studies. We show through simulation that current tests for CNV association are prone to false-positive associations in the presence of differential errors between cases and controls, especially if quantitative CNV measurements are noisy. We present a statistical framework for performing case-control CNV association studies that applies likelihood ratio testing of quantitative CNV measurements in cases and controls. We show that our methods are robust to differential errors and noisy data and can achieve maximal theoretical power. We illustrate the power of these methods for testing for association with binary and quantitative traits, and have made this software available as the R package CNVtools.     

Rare independent mutations in renal salt handling genes contribute to blood pressure variation

The effects of alleles in many genes are believed to contribute to common complex diseases such as hypertension. Whether risk alleles comprise a small number of common variants or many rare independent mutations at trait loci is largely unknown. We screened members of the Framingham Heart Study (FHS) for variation in three genes-SLC12A3 (NCCT), SLC12A1 (NKCC2) and KCNJ1 (ROMK)-causing rare recessive diseases featuring large reductions in blood pressure. Using comparative genomics, genetics and biochemistry, we identified subjects with mutations proven or inferred to be functional. These mutations, all heterozygous and rare, produce clinically significant blood pressure reduction and protect from development of hypertension. Our findings implicate many rare alleles that alter renal salt handling in blood pressure variation in the general population, and identify alleles with health benefit that are nonetheless under purifying selection. These findings have implications for the genetic architecture of hypertension and other common complex traits.     

Multiple regions within 8q24 independently affect risk for prostate cancer

After the recent discovery that common genetic variation in 8q24 influences inherited risk of prostate cancer, we genotyped 2,973 SNPs in up to 7,518 men with and without prostate cancer from five populations. We identified seven risk variants, five of them previously undescribed, spanning 430 kb and each independently predicting risk for prostate cancer (P = 7.9 x 10(-19) for the strongest association, and P < 1.5 x 10(-4) for five of the variants, after controlling for each of the others). The variants define common genotypes that span a more than fivefold range of susceptibility to cancer in some populations. None of the prostate cancer risk variants aligns to a known gene or alters the coding sequence of an encoded protein.     

Positional cloning of Sorcs1, a type 2 diabetes quantitative trait locus

We previously mapped the type 2 diabetes mellitus-2 locus (T2dm2), which affects fasting insulin levels, to distal chromosome 19 in a leptin-deficient obese F2 intercross derived from C57BL/6 (B6) and BTBR T+ tf/J (BTBR) mice. Introgression of a 7-Mb segment of the B6 chromosome 19 into the BTBR background (strain 1339A) replicated the reduced insulin linked to T2dm2. The 1339A mice have markedly impaired insulin secretion in vivo and disrupted islet morphology. We used subcongenic strains derived from 1339A to localize the T2dm2 quantitative trait locus (QTL) to a 242-kb segment comprising the promoter, first exon and most of the first intron of the Sorcs1 gene. This was the only gene in the 1339A strain for which we detected amino acid substitutions and expression level differences between mice carrying B6 and BTBR alleles of this insert, thereby identifying variation within the Sorcs1 gene as underlying the phenotype associated with the T2dm2 locus. SorCS1 binds platelet-derived growth factor, a growth factor crucial for pericyte recruitment to the microvasculature, and may thus have a role in expanding or maintaining the islet vasculature. Our identification of the Sorcs1 gene provides insight into the pathway underlying the pathophysiology of obesity-induced type 2 diabetes mellitus.     

A recurrent mutation in the BMP type I receptor ACVR1 causes inherited and sporadic fibrodysplasia ossificans progressiva

Fibrodysplasia ossificans progressiva (FOP) is a rare autosomal dominant disorder of skeletal malformations and progressive extraskeletal ossification. We mapped FOP to chromosome 2q23-24 by linkage analysis and identified an identical heterozygous mutation (617G --> A; R206H) in the glycine-serine (GS) activation domain of ACVR1, a BMP type I receptor, in all affected individuals examined. Protein modeling predicts destabilization of the GS domain, consistent with constitutive activation of ACVR1 as the underlying cause of the ectopic chondrogenesis, osteogenesis and joint fusions seen in FOP.     

Epigenetic characterization of the early embryo with a chromatin immunoprecipitation protocol applicable to small cell populations

Chromatin immunoprecipitation (ChIP) defines the genomic distribution of proteins and their modifications but is limited by the cell numbers required (ideally >10(7)). Here we describe a protocol that uses carrier chromatin and PCR, 'carrier' ChIP (CChIP), to permit analysis of as few as 100 cells. We assayed histone modifications at key regulator genes (such as Nanog, Pou5f1 (also known as Oct4) and Cdx2) by CChIP in mouse embryonic stem (ES) cells and in inner cell mass (ICM) and trophectoderm of cultured blastocysts. Activating and silencing modifications (H4 acetylation and H3K9 methylation) mark active and silent promoters as predicted, and we find close correlation between values derived from CChIP (1,000 ES cells) and conventional ChIP (5 x 10(7) ES cells). Studies on genes silenced in both ICM and ES cells (Cdx2, Cfc1, Hhex and Nkx2-2, also known as Nkx) show that the intensity of silencing marks is relatively diminished in ES cells, indicating a possible relaxation of some components of silencing on adaptation to culture.     

Automating resequencing-based detection of insertion-deletion polymorphisms

Structural and insertion-deletion (indel) variants have received considerable recent attention, partly because of their phenotypic consequences. Among these variants, the most common are small indels ( approximately 1-30 bp). Identifying and genotyping indels using sequence traces obtained from diploid samples requires extensive manual review, which makes large-scale studies inconvenient. We report a new algorithm, implemented in available software (PolyPhred version 6.0), to help automate detection and genotyping of indels from sequence traces. The algorithm identifies heterozygous individuals, which permits the discovery of low-frequency indels. It finds 80% of all indel polymorphisms with almost no false positives and finds 97% with a false discovery rate of 10%. Additionally, genotyping accuracy exceeds 99%, and it correctly infers indel length in 96% of the cases. Using this approach, we identify indels in the HapMap ENCODE regions, providing the first report of these polymorphisms in this data set.