Genotype, haplotype and copy-number variation in worldwide human populations.

Genome-wide patterns of variation across individuals provide a powerful source of data for uncovering the history of migration, range expansion, and adaptation of the human species. However, high-resolution surveys of variation in genotype, haplotype and copy number have generally focused on a small number of population groups. Here we report the analysis of high-quality genotypes at 525,910 single-nucleotide polymorphisms (SNPs) and 396 copy-number-variable loci in a worldwide sample of 29 populations. Analysis of SNP genotypes yields strongly supported fine-scale inferences about population structure. Increasing linkage disequilibrium is observed with increasing geographic distance from Africa, as expected under a serial founder effect for the out-of-Africa spread of human populations. New approaches for haplotype analysis produce inferences about population structure that complement results based on unphased SNPs. Despite a difference from SNPs in the frequency spectrum of the copy-number variants (CNVs) detected--including a comparatively large number of CNVs in previously unexamined populations from Oceania and the Americas--the global distribution of CNVs largely accords with population structure analyses for SNP data sets of similar size. Our results produce new inferences about inter-population variation, support the utility of CNVs in human population-genetic research, and serve as a genomic resource for human-genetic studies in diverse worldwide populations.

Abundant quantitative trait loci exist for DNA methylation and gene expression in human brain.

A fundamental challenge in the post-genome era is to understand and annotate the consequences of genetic variation, particularly within the context of human tissues. We present a set of integrated experiments that investigate the effects of common genetic variability on DNA methylation and mRNA expression in four human brain regions each from 150 individuals (600 samples total). We find an abundance of genetic cis regulation of mRNA expression and show for the first time abundant quantitative trait loci for DNA CpG methylation across the genome. We show peak enrichment for cis expression QTLs to be approximately 68,000 bp away from individual transcription start sites; however, the peak enrichment for cis CpG methylation QTLs is located much closer, only 45 bp from the CpG site in question. We observe that the largest magnitude quantitative trait loci occur across distinct brain tissues. Our analyses reveal that CpG methylation quantitative trait loci are more likely to occur for CpG sites outside of islands. Lastly, we show that while we can observe individual QTLs that appear to affect both the level of a transcript and a physically close CpG methylation site, these are quite rare. We believe these data, which we have made publicly available, will provide a critical step toward understanding the biological effects of genetic variation.

A two-stage genome-wide association study of sporadic amyotrophic lateral sclerosis.

The cause of sporadic amyotrophic lateral sclerosis (ALS) is largely unknown, but genetic factors are thought to play a significant role in determining susceptibility to motor neuron degeneration. To identify genetic variants altering risk of ALS, we undertook a two-stage genome-wide association study (GWAS): we followed our initial GWAS of 545 066 SNPs in 553 individuals with ALS and 2338 controls by testing the 7600 most associated SNPs from the first stage in three independent cohorts consisting of 2160 cases and 3008 controls. None of the SNPs selected for replication exceeded the Bonferroni threshold for significance. The two most significantly associated SNPs, rs2708909 and rs2708851 [odds ratio (OR) = 1.17 and 1.18, and P-values = 6.98 x 10(-7) and 1.16 x 10(-6)], were located on chromosome 7p13.3 within a 175 kb linkage disequilibrium block containing the SUNC1, HUS1 and C7orf57 genes. These associations did not achieve genome-wide significance in the original cohort and failed to replicate in an additional independent cohort of 989 US cases and 327 controls (OR = 1.18 and 1.19, P-values = 0.08 and 0.06, respectively). Thus, we chose to cautiously interpret our data as hypothesis-generating requiring additional confirmation, especially as all previously reported loci for ALS have failed to replicate successfully. Indeed, the three loci (FGGY, ITPR2 and DPP6) identified in previous GWAS of sporadic ALS were not significantly associated with disease in our study. Our findings suggest that ALS is more genetically and clinically heterogeneous than previously recognized. Genotype data from our study have been made available online to facilitate such future endeavors.

A genome-wide association study identifies protein quantitative trait loci (pQTLs).

There is considerable evidence that human genetic variation influences gene expression. Genome-wide studies have revealed that mRNA levels are associated with genetic variation in or close to the gene coding for those mRNA transcripts - cis effects, and elsewhere in the genome - trans effects. The role of genetic variation in determining protein levels has not been systematically assessed. Using a genome-wide association approach we show that common genetic variation influences levels of clinically relevant proteins in human serum and plasma. We evaluated the role of 496,032 polymorphisms on levels of 42 proteins measured in 1200 fasting individuals from the population based InCHIANTI study. Proteins included insulin, several interleukins, adipokines, chemokines, and liver function markers that are implicated in many common diseases including metabolic, inflammatory, and infectious conditions. We identified eight Cis effects, including variants in or near the IL6R (p = 1.8x10(-57)), CCL4L1 (p = 3.9x10(-21)), IL18 (p = 6.8x10(-13)), LPA (p = 4.4x10(-10)), GGT1 (p = 1.5x10(-7)), SHBG (p = 3.1x10(-7)), CRP (p = 6.4x10(-6)) and IL1RN (p = 7.3x10(-6)) genes, all associated with their respective protein products with effect sizes ranging from 0.19 to 0.69 standard deviations per allele. Mechanisms implicated include altered rates of cleavage of bound to unbound soluble receptor (IL6R), altered secretion rates of different sized proteins (LPA), variation in gene copy number (CCL4L1) and altered transcription (GGT1). We identified one novel trans effect that was an association between ABO blood group and tumour necrosis factor alpha (TNF-alpha) levels (p = 6.8x10(-40)), but this finding was not present when TNF-alpha was measured using a different assay , or in a second study, suggesting an assay-specific association. Our results show that protein levels share some of the features of the genetics of gene expression. These include the presence of strong genetic effects in cis locations. The identification of protein quantitative trait loci (pQTLs) may be a powerful complementary method of improving our understanding of disease pathways.

A duplication at chromosome 11q12.2-11q12.3 is associated with spinocerebellar ataxia type 20.

Spinocerebellar ataxia type 20 (SCA20) has been linked to chromosome 11q12, but the underlying genetic defect has yet to be identified. We applied single-nucleotide polymorphism genotyping to detect structural alterations in the genomic DNA of patients with SCA20. We found a 260 kb duplication within the previously linked SCA20 region, which was confirmed by quantitative polymerase chain reaction and fiber fluorescence in situ hybridization, the latter also showing its direct orientation. The duplication spans 10 known and 2 unknown genes, and is present in all affected individuals in the single reported SCA20 pedigree. While the mechanism whereby this duplication may be pathogenic remains to be established, we speculate that the critical gene within the duplicated segment may be DAGLA, the product of which is normally present at the base of Purkinje cell dendritic spines and contributes to the modulation of parallel fiber-Purkinje cell synapses.

Sequencing analysis of the ITPR1 gene in a pure autosomal dominant spinocerebellar ataxia series.

Spinocerebellar ataxia type 15 and 16 (SCA15/16) are autosomal dominant cerebellar ataxias that are slowly progressive with a predominantly pure ataxia phenotype (ADCA III). The locus for SCA15 was first mapped to 3p24.2-3pter and subsequently full or partial deletions in the inositol 1,4,5-triphosphate receptor type 1 (ITPR1) gene were identified in several ADCA III families that segregated with the disease. A single missense coding variant has been described, but the pathogenicity of this change has not been proven. We sequenced the entire coding region and flanking regions of ITPR1 in unrelated ADCA III families (n = 38) that were negative for large deletions on whole genome arrays, and for which SCAs 1, 2, 3, 6, 7, 8, 11, 12, 14, 17 and the Friedreich's ataxia expansion were excluded in all probands. Mutation at SCA5, 10, and 27 was also excluded in some families. A number of coding and noncoding polymorphisms were identified but no ITPR1 mutations were found. The results indicate that point mutations in ITPR1 are at best a rare cause of ADCA III.

Deletion at ITPR1 underlies ataxia in mice and spinocerebellar ataxia 15 in humans.

We observed a severe autosomal recessive movement disorder in mice used within our laboratory. We pursued a series of experiments to define the genetic lesion underlying this disorder and to identify a cognate disease in humans with mutation at the same locus. Through linkage and sequence analysis we show here that this disorder is caused by a homozygous in-frame 18-bp deletion in Itpr1 (Itpr1(Delta18/Delta18)), encoding inositol 1,4,5-triphosphate receptor 1. A previously reported spontaneous Itpr1 mutation in mice causes a phenotype identical to that observed here. In both models in-frame deletion within Itpr1 leads to a decrease in the normally high level of Itpr1 expression in cerebellar Purkinje cells. Spinocerebellar ataxia 15 (SCA15), a human autosomal dominant disorder, maps to the genomic region containing ITPR1; however, to date no causal mutations had been identified. Because ataxia is a prominent feature in Itpr1 mutant mice, we performed a series of experiments to test the hypothesis that mutation at ITPR1 may be the cause of SCA15. We show here that heterozygous deletion of the 5' part of the ITPR1 gene, encompassing exons 1-10, 1-40, and 1-44 in three studied families, underlies SCA15 in humans.

Twin-rainbow metrology. I. Measurement of the thickness of a thin liquid film draining under gravity.

We describe twin-rainbow metrology, a new optical technique used to measure the thickness of thin films in a cylindrical geometry. We also present an application of the technique: measurement of the thickness of a Newtonian fluid draining under gravity. We compare these measurements with fluid mechanics models.