Identification of common variants associated with human hippocampal and intracranial volumes

Identifying genetic variants influencing human brain structures may reveal new biological mechanisms underlying cognition and neuropsychiatric illness. The volume of the hippocampus is a biomarker of incipient Alzheimer's disease and is reduced in schizophrenia, major depression and mesial temporal lobe epilepsy. Whereas many brain imaging phenotypes are highly heritable, identifying and replicating genetic influences has been difficult, as small effects and the high costs of magnetic resonance imaging (MRI) have led to underpowered studies. Here we report genome-wide association meta-analyses and replication for mean bilateral hippocampal, total brain and intracranial volumes from a large multinational consortium. The intergenic variant rs7294919 was associated with hippocampal volume (12q24.22; N = 21,151; P = 6.70 × 10(-16)) and the expression levels of the positional candidate gene TESC in brain tissue. Additionally, rs10784502, located within HMGA2, was associated with intracranial volume (12q14.3; N = 15,782; P = 1.12 × 10(-12)). We also identified a suggestive association with total brain volume at rs10494373 within DDR2 (1q23.3; N = 6,500; P = 5.81 × 10(-7)).     

Age distribution of human gene families shows significant roles of both large- and small-scale duplications in vertebrate evolution

The zebrafish embryo is transparent and can tolerate absence of blood flow because its oxygen is delivered by diffusion rather than by the cardiovascular system. It is therefore possible to attribute cardiac failure directly to particular genes by ruling out the possibility that it is due to a secondary effect of hypoxia. We focus here on pickwickm171 (pikm171), a recessive lethal mutation discovered in a large-scale genetic screen. There are three other alleles in the pik complementation group with this phenotype (pikm242, pikm740, pikm186; ref. 3) and one allele (pikmVO62H) with additional skeletal paralysis. The pik heart develops normally but is poorly contractile from the first beat. Aside from the edema that inevitably accompanies cardiac dysfunction, development is normal during the first three days. We show by positional cloning that the 'causative' mutation is in an alternatively-spliced exon of the gene (ttn) encoding Titin. Titin is the biggest known protein and spans the half-sarcomere from Z-disc to M-line in heart and skeletal muscle. It has been proposed to provide a scaffold for the assembly of thick and thin filaments and to provide elastic recoil engendered by stretch during diastole. We found that nascent myofibrils form in pik mutants, but normal sarcomeres are absent. Mutant cells transplanted to wildtype hearts remain thin and bulge outwards as individual cell aneurysms without affecting nearby wildtype cardiomyocytes, indicating that the contractile deficiency is cell-autonomous. Absence of Titin function thus results in blockage of sarcomere assembly and causes a functional disorder resembling human dilated cardiomyopathies, one form of which is described in another paper in this issue.     

Identification of ten loci associated with height highlights new biological pathways in human growth

Height is a classic polygenic trait, reflecting the combined influence of multiple as-yet-undiscovered genetic factors. We carried out a meta-analysis of genome-wide association study data of height from 15,821 individuals at 2.2 million SNPs, and followed up the strongest findings in >10,000 subjects. Ten newly identified and two previously reported loci were strongly associated with variation in height (P values from 4 x 10(-7) to 8 x 10(-22)). Together, these 12 loci account for approximately 2% of the population variation in height. Individuals with < or =8 height-increasing alleles and > or =16 height-increasing alleles differ in height by approximately 3.5 cm. The newly identified loci, along with several additional loci with strongly suggestive associations, encompass both strong biological candidates and unexpected genes, and highlight several pathways (let-7 targets, chromatin remodeling proteins and Hedgehog signaling) as important regulators of human stature. These results expand the picture of the biological regulation of human height and of the genetic architecture of this classical complex trait.     

Loss-of-function mutations in the EGF-CFC gene CFC1 are associated with human left-right laterality defects

All vertebrates display a characteristic asymmetry of internal organs with the cardiac apex, stomach and spleen towards the left, and the liver and gall bladder on the right. Left-right (L-R) axis abnormalities or laterality defects are common in humans (1 in 8,500 live births). Several genes (such as Nodal, Ebaf and Pitx2) have been implicated in L-R organ positioning in model organisms. In humans, relatively few genes have been associated with a small percentage of human situs defects. These include ZIC3 (ref. 5), LEFTB (formerly LEFTY2; ref. 6) and ACVR2B (encoding activin receptor IIB; ref. 7). The EGF-CFC genes, mouse Cfc1 (encoding the Cryptic protein; ref. 9) and zebrafish one-eyed pinhead (oep; refs 10, 11) are essential for the establishment of the L-R axis. EGF-CFC proteins act as co-factors for Nodal-related signals, which have also been implicated in L-R axis development. Here we identify loss-of-function mutations in human CFC1 (encoding the CRYPTIC protein) in patients with heterotaxic phenotypes (randomized organ positioning). The mutant proteins have aberrant cellular localization in transfected cells and are functionally defective in a zebrafish oep-mutant rescue assay. Our findings indicate that the essential role of EGF-CFC genes and Nodal signalling in left-right axis formation is conserved from fish to humans. Moreover, our results support a role for environmental and/or genetic modifiers in determining the ultimate phenotype in humans.     

Hereditary hypertension caused by chimaeric gene duplications and ectopic expression of aldosterone synthase.

Patients with glucocorticoid-remediable aldosteronism (GRA) from 12 kindreds possess chimaeric gene duplications arising from unequal crossing-over, fusing regulatory sequences of steroid 11 beta-hydroxylase to coding sequences of aldosterone synthase. These chimaeric genes are specific for GRA and explain the biochemistry, physiology and genetics of this form of hypertension. Sites of crossing over range from intron 2 to intron 4. Most mutations have arisen independently from either sister or non-sister chromatid exchange between these genes, which are only 45 kilobases apart. The possibility of a susceptibility allele for GRA of Irish origin is suggested. These findings indicate the utility of a direct genetic test for this disorder.     

Defective colour vision associated with a missense mutation in the human green visual pigment gene.

All red/green colour vision defects described so far have been associated with gross rearrangements within the red/green opsin gene array (Xq28). We now describe a male with severe deuteranomaly without such a rearrangement. A substitution of a highly conserved cysteine by arginine at position 203 in the green opsins presumably accounted for his colour vision defect. Surprisingly, this mutation was fairly common (2%) in the population but apparently was not always expressed. In analogy with nonexpression of some 5'green-red hybrid genes in persons with normal colour vision, we suggest that failure of manifestation occurs when the mutant gene is located at a distal (3') position among several green opsin genes. This mutation might also predispose to certain X-linked retinal dystrophies.     

Identification of foreign gene sequences by transcript filtering against the human genome.

We have developed a computational subtraction approach to detect microbial causes for putative infectious diseases by filtering a set of human tissue-derived sequences against the human genome. We demonstrate the potential of this method by identifying sequences from known pathogens in established expressed-sequence tag libraries.     

Identification of loci associated with schizophrenia by genome-wide association and follow-up

We carried out a genome-wide association study of schizophrenia (479 cases, 2,937 controls) and tested loci with P < 10(-5) in up to 16,726 additional subjects. Of 12 loci followed up, 3 had strong independent support (P < 5 x 10(-4)), and the overall pattern of replication was unlikely to occur by chance (P = 9 x 10(-8)). Meta-analysis provided strongest evidence for association around ZNF804A (P = 1.61 x 10(-7)) and this strengthened when the affected phenotype included bipolar disorder (P = 9.96 x 10(-9)).     

Bayesian inference of ancient human demography from individual genome sequences

Whole-genome sequences provide a rich source of information about human evolution. Here we describe an effort to estimate key evolutionary parameters based on the whole-genome sequences of six individuals from diverse human populations. We used a Bayesian, coalescent-based approach to obtain information about ancestral population sizes, divergence times and migration rates from inferred genealogies at many neutrally evolving loci across the genome. We introduce new methods for accommodating gene flow between populations and integrating over possible phasings of diploid genotypes. We also describe a custom pipeline for genotype inference to mitigate biases from heterogeneous sequencing technologies and coverage levels. Our analysis indicates that the San population of southern Africa diverged from other human populations approximately 108-157 thousand years ago, that Eurasians diverged from an ancestral African population 38-64 thousand years ago, and that the effective population size of the ancestors of all modern humans was ～9,000.     

Identification of low-frequency variants associated with gout and serum uric acid levels

We tested 16 million SNPs, identified through whole-genome sequencing of 457 Icelanders, for association with gout and serum uric acid levels. Genotypes were imputed into 41,675 chip-genotyped Icelanders and their relatives, for effective sample sizes of 968 individuals with gout and 15,506 individuals for whom serum uric acid measurements were available. We identified a low-frequency missense variant (c.1580C>G) in ALDH16A1 associated with gout (OR = 3.12, P = 1.5 × 10(-16), at-risk allele frequency = 0.019) and serum uric acid levels (effect = 0.36 s.d., P = 4.5 × 10(-21)). We confirmed the association with gout by performing Sanger sequencing on 6,017 Icelanders. The association with gout was stronger in males relative to females. We also found a second variant on chromosome 1 associated with gout (OR = 1.92, P = 0.046, at-risk allele frequency = 0.986) and serum uric acid levels (effect = 0.48 s.d., P = 4.5 × 10(-16)). This variant is close to a common variant previously associated with serum uric acid levels. This work illustrates how whole-genome sequencing data allow the detection of associations between low-frequency variants and complex traits.     

Identification of the gene that, when mutated, causes the human obesity syndrome BBS4

Bardet-Biedl syndrome (BBS, MIM 209900) is a heterogeneous autosomal recessive disorder characterized by obesity, pigmentary retinopathy, polydactyly, renal malformations, mental retardation, and hypogenitalism. The disorder is also associated with diabetes mellitus, hypertension, and congenital heart disease. Six distinct BBS loci map to 11q13 (BBS1), 16q21 (BBS2), 3p13-p12 (BBS3), 15q22.3-q23 (BBS4), 2q31 (BBS5), and 20p12 (BBS6). Although BBS is rare in the general population (<1/100,000), there is considerable interest in identifying the genes causing BBS because components of the phenotype, such as obesity and diabetes, are common. We and others have demonstrated that BBS6 is caused by mutations in the gene MKKS (refs. 12,13), mutation of which also causes McKusick-Kaufman syndrome (hydrometrocolpos, post-axial polydactyly, and congenital heart defects). MKKS has sequence homology to the alpha subunit of a prokaryotic chaperonin in the thermosome Thermoplasma acidophilum. We recently identified a novel gene that causes BBS2. The BBS2 protein has no significant similarity to other chaperonins or known proteins. Here we report the positional cloning and identification of mutations in BBS patients in a novel gene designated BBS4.     

Identification of the familial cylindromatosis tumour-suppressor gene

Familial cylindromatosis is an autosomal dominant genetic predisposition to multiple tumours of the skin appendages. The susceptibility gene (CYLD) has previously been localized to chromosome 16q and has the genetic attributes of a tumour-suppressor gene (recessive oncogene). Here we have identified CYLD by detecting germline mutations in 21 cylindromatosis families and somatic mutations in 1 sporadic and 5 familial cylindromas. All mutations predict truncation or absence of the encoded protein. CYLD encodes three cytoskeletal-associated-protein-glycine-conserved (CAP-GLY) domains, which are found in proteins that coordinate the attachment of organelles to microtubules. CYLD also has sequence homology to the catalytic domain of ubiquitin carboxy-terminal hydrolases (UCH).     

Identification of Stk6/STK15 as a candidate low-penetrance tumor-susceptibility gene in mouse and human

Linkage analysis and haplotype mapping in interspecific mouse crosses (Mus musculus x Mus spretus) identified the gene encoding Aurora2 (Stk6 in mouse and STK15 in human) as a candidate skin tumor susceptibility gene. The Stk6 allele inherited from the susceptible M. musculus parent was overexpressed in normal cells and preferentially amplified in tumor cells from F(1) hybrid mice. We identified a common genetic variant in STK15 (resulting in the amino acid substitution F31I) that is preferentially amplified and associated with the degree of aneuploidy in human colon tumors. The Ile31 variant transforms rat1 cells more potently than the more common Phe31 variant. The E2 ubiquitin-conjugating enzyme UBE2N was a preferential binding partner of the 'weak' STK15 Phe31 variant form in yeast two-hybrid screens and in human cells. This interaction results in colocalization of UBE2N with STK15 at the centrosomes during mitosis. These results are consistent with an important role for the Ile31 variant of STK15 in human cancer susceptibility.     

Autosomal recessive lissencephaly with cerebellar hypoplasia is associated with human RELN mutations

Normal development of the cerebral cortex requires long-range migration of cortical neurons from proliferative regions deep in the brain. Lissencephaly ("smooth brain," from "lissos," meaning smooth, and "encephalos," meaning brain) is a severe developmental disorder in which neuronal migration is impaired, leading to a thickened cerebral cortex whose normally folded contour is simplified and smooth. Two identified lissencephaly genes do not account for all known cases, and additional lissencephaly syndromes have been described. An autosomal recessive form of lissencephaly (LCH) associated with severe abnormalities of the cerebellum, hippocampus and brainstem maps to chromosome 7q22, and is associated with two independent mutations in the human gene encoding reelin (RELN). The mutations disrupt splicing of RELN cDNA, resulting in low or undetectable amounts of reelin protein. LCH parallels the reeler mouse mutant (Reln(rl)), in which Reln mutations cause cerebellar hypoplasia, abnormal cerebral cortical neuronal migration and abnormal axonal connectivity. RELN encodes a large (388 kD) secreted protein that acts on migrating cortical neurons by binding to the very low density lipoprotein receptor (VLDLR), the apolipoprotein E receptor 2 (ApoER2; refs 9-11 ), alpha3beta1 integrin and protocadherins. Although reelin was previously thought to function exclusively in brain, some humans with RELN mutations show abnormal neuromuscular connectivity and congenital lymphoedema, suggesting previously unsuspected functions for reelin in and outside of the brain.     

Biochemical defects in ABCR protein variants associated with human retinopathies

Mutations in the gene encoding ABCR (ABCA4), a photoreceptor-specific ATP-binding cassette (ABC) transporter, are responsible for autosomal recessive Stargardt disease (STGD), an early onset macular degeneration, and some forms of autosomal recessive cone-rod dystrophy and autosomal recessive retinitis pigmentosa. Heterozygosity for ABCA4 mutations may also represent a risk factor for age-related macular degeneration (AMD), although this idea is controversial. An ongoing challenge in the analysis of ABCA4-based retinopathies arises from the observation that most of the ABCA4 sequence variants identified so far are missense mutations that are rare in both patient and control populations. With the current sample size of most sequence variants, one cannot determine statistically whether a particular sequence variant is pathogenic or neutral. A related challenge is to determine the degree to which each pathogenic variant impairs ABCR function, as genotype-phenotype analyses indicate that age of onset and disease severity correlate with different ABCA4 alleles. To address these questions, we performed a functional analysis of human ABCR and its variants. These experiments reveal a wide spectrum of biochemical defects in these variants and provide insight into the transport mechanism of ABCR.     

Identification of the gene causing mucolipidosis type IV

Mucolipidosis type IV (MLIV) is an autosomal recessive, neurodegenerative, lysosomal storage disorder characterized by psychomotor retardation and ophthalmological abnormalities including corneal opacities, retinal degeneration and strabismus. Most patients reach a maximal developmental level of 12?15 months. The disease was classified as a mucolipidosis following observations by electron microscopy indicating the lysosomal storage of lipids together with water-soluble, granulated substances. Over 80% of the MLIV patients diagnosed are Ashkenazi Jews, including severely affected and mildly affected patients. The gene causing MLIV was previously mapped to human chromosome 19p13.2-13.3 in a region of approximately 1 cM (ref. 7). Haplotype analysis in the MLIV gene region of over 70 MLIV Ashkenazi chromosomes indicated the existence of two founder chromosomes among 95% of the Ashkenazi MLIV families: a major haplotype in 72% and a minor haplotype in 23% of the MLIV chromosomes (ref. 7, and G.B., unpublished data). The remaining 5% are distinct haplotypes found only in single patients. The basic metabolic defect causing the lysosomal storage in MLIV has not yet been identified. Thus, positional cloning was an alternative to identify the MLIV gene. We report here the identification of a new gene in this human chromosomal region in which MLIV-specific mutations were identified.     

Genome-wide association study identifies a variant in HDAC9 associated with large vessel ischemic stroke

Genetic factors have been implicated in stroke risk, but few replicated associations have been reported. We conducted a genome-wide association study (GWAS) for ischemic stroke and its subtypes in 3,548 affected individuals and 5,972 controls, all of European ancestry. Replication of potential signals was performed in 5,859 affected individuals and 6,281 controls. We replicated previous associations for cardioembolic stroke near PITX2 and ZFHX3 and for large vessel stroke at a 9p21 locus. We identified a new association for large vessel stroke within HDAC9 (encoding histone deacetylase 9) on chromosome 7p21.1 (including further replication in an additional 735 affected individuals and 28,583 controls) (rs11984041; combined P = 1.87 × 10(-11); odds ratio (OR) = 1.42, 95% confidence interval (CI) = 1.28-1.57). All four loci exhibited evidence for heterogeneity of effect across the stroke subtypes, with some and possibly all affecting risk for only one subtype. This suggests distinct genetic architectures for different stroke subtypes.