Sanger Institute - Publications 2008

Number of papers published in 2008: 142

  • EYS, encoding an ortholog of Drosophila spacemaker, is mutated in autosomal recessive retinitis pigmentosa.

    Abd El-Aziz MM, Barragan I, O'Driscoll CA, Goodstadt L, Prigmore E, Borrego S, Mena M, Pieras JI, El-Ashry MF, Safieh LA, Shah A, Cheetham ME, Carter NP, Chakarova C, Ponting CP, Bhattacharya SS and Antinolo G

    Department of Molecular Genetics, Institute of Ophthalmology, London EC1V 9EL, UK.

    Using a positional cloning approach supported by comparative genomics, we have identified a previously unreported gene, EYS, at the RP25 locus on chromosome 6q12 commonly mutated in autosomal recessive retinitis pigmentosa. Spanning over 2 Mb, this is the largest eye-specific gene identified so far. EYS is independently disrupted in four other mammalian lineages, including that of rodents, but is well conserved from Drosophila to man and is likely to have a role in the modeling of retinal architecture.

    Funded by: Medical Research Council: MC_U137761446; Wellcome Trust: 077008

    Nature genetics 2008;40;11;1285-7

  • Contemporary approaches for modifying the mouse genome.

    Adams DJ and van der Weyden L

    Experimental Cancer Genetics, The Wellcome Trust Sanger Institute, Hinxton, United Kingdom.

    The mouse is a premiere experimental organism that has contributed significantly to our understanding of vertebrate biology. Manipulation of the mouse genome via embryonic stem (ES) cell technology makes it possible to engineer an almost limitless repertoire of mutations to model human disease and assess gene function. In this review we outline recent advances in mouse experimental genetics and provide a "how-to" guide for those people wishing to access this technology. We also discuss new technologies, such as transposon-mediated mutagenesis, and resources of targeting vectors and ES cells, which are likely to dramatically accelerate the pace with which we can assess gene function in vivo, and the progress of forward and reverse genetic screens in mice.

    Funded by: Cancer Research UK; Wellcome Trust

    Physiological genomics 2008;34;3;225-38

  • Genomic-scale prioritization of drug targets: the TDR Targets database.

    Agüero F, Al-Lazikani B, Aslett M, Berriman M, Buckner FS, Campbell RK, Carmona S, Carruthers IM, Chan AW, Chen F, Crowther GJ, Doyle MA, Hertz-Fowler C, Hopkins AL, McAllister G, Nwaka S, Overington JP, Pain A, Paolini GV, Pieper U, Ralph SA, Riechers A, Roos DS, Sali A, Shanmugam D, Suzuki T, Van Voorhis WC and Verlinde CL

    Instituto de Investigaciones Biotecnológicas, Universidad Nacional de General San Martín, San Martín 1650, Buenos Aires, Argentina.

    The increasing availability of genomic data for pathogens that cause tropical diseases has created new opportunities for drug discovery and development. However, if the potential of such data is to be fully exploited, the data must be effectively integrated and be easy to interrogate. Here, we discuss the development of the TDR Targets database (, which encompasses extensive genetic, biochemical and pharmacological data related to tropical disease pathogens, as well as computationally predicted druggability for potential targets and compound desirability information. By allowing the integration and weighting of this information, this database aims to facilitate the identification and prioritization of candidate drug targets for pathogens.

    Funded by: NIGMS NIH HHS: R01 GM054762

    Nature reviews. Drug discovery 2008;7;11;900-7

  • Data growth and its impact on the SCOP database: new developments.

    Andreeva A, Howorth D, Chandonia JM, Brenner SE, Hubbard TJ, Chothia C and Murzin AG

    MRC Centre for Protein Engineering, Hills Road, Cambridge CB2 0QH, UK.

    The Structural Classification of Proteins (SCOP) database is a comprehensive ordering of all proteins of known structure, according to their evolutionary and structural relationships. The SCOP hierarchy comprises the following levels: Species, Protein, Family, Superfamily, Fold and Class. While keeping the original classification scheme intact, we have changed the production of SCOP in order to cope with a rapid growth of new structural data and to facilitate the discovery of new protein relationships. We describe ongoing developments and new features implemented in SCOP. A new update protocol supports batch classification of new protein structures by their detected relationships at Family and Superfamily levels in contrast to our previous sequential handling of new structural data by release date. We introduce pre-SCOP, a preview of the SCOP developmental version that enables earlier access to the information on new relationships. We also discuss the impact of worldwide Structural Genomics initiatives, which are producing new protein structures at an increasing rate, on the rates of discovery and growth of protein families and superfamilies. SCOP can be accessed at

    Funded by: Medical Research Council: G0100305, MC_U105192716; NIGMS NIH HHS: R01 GM073109, R01-GM073109; Wellcome Trust: 077198

    Nucleic acids research 2008;36;Database issue;D419-25

  • Consistently replicating locus linked to migraine on 10q22-q23.

    Anttila V, Nyholt DR, Kallela M, Artto V, Vepsäläinen S, Jakkula E, Wennerström A, Tikka-Kleemola P, Kaunisto MA, Hämäläinen E, Widén E, Terwilliger J, Merikangas K, Montgomery GW, Martin NG, Daly M, Kaprio J, Peltonen L, Färkkilä M, Wessman M and Palotie A

    Biomedicum Helsinki, Research Program in Molecular Medicine, University of Helsinki, 00290 Helsinki, Finland.

    Here, we present the results of two genome-wide scans in two diverse populations in which a consistent use of recently introduced migraine-phenotyping methods detects and replicates a locus on 10q22-q23, with an additional independent replication. No genetic variants have been convincingly established in migraine, and although several loci have been reported, none of them has been consistently replicated. We employed the three known migraine-phenotyping methods (clinical end diagnosis, latent-class analysis, and trait-component analysis) with robust multiple testing correction in a large sample set of 1675 individuals from 210 migraine families from Finland and Australia. Genome-wide multipoint linkage analysis that used the Kong and Cox exponential model in Finns detected a locus on 10q22-q23 with highly significant evidence of linkage (LOD 7.68 at 103 cM in female-specific analysis). The Australian sample showed a LOD score of 3.50 at the same locus (100 cM), as did the independent Finnish replication study (LOD score 2.41, at 102 cM). In addition, four previously reported loci on 8q21, 14q21, 18q12, and Xp21 were also replicated. A shared-segment analysis of 10q22-q23 linked Finnish families identified a 1.6-9.5 cM segment, centered on 101 cM, which shows in-family homology in 95% of affected Finns. This region was further studied with 1323 SNPs. Although no significant association was observed, four regions warranting follow-up studies were identified. These results support the use of symptomology-based phenotyping in migraine and suggest that the 10q22-q23 locus probably contains one or more migraine susceptibility variants.

    Funded by: NCRR NIH HHS: U54 RR020278; NIAAA NIH HHS: AA013320, AA013326, AA014041, AA07728, AA11998, P50 AA011998, R01 AA007535, R01 AA007728, R01 AA010249, R01 AA013320, R01 AA013326, R01 AA014041, R37 AA007728; NINDS NIH HHS: R01 NS037675, R01 NS37675

    American journal of human genetics 2008;82;5;1051-63

  • Loci influencing lipid levels and coronary heart disease risk in 16 European population cohorts.

    Aulchenko YS, Ripatti S, Lindqvist I, Boomsma D, Heid IM, Pramstaller PP, Penninx BW, Janssens AC, Wilson JF, Spector T, Martin NG, Pedersen NL, Kyvik KO, Kaprio J, Hofman A, Freimer NB, Jarvelin MR, Gyllensten U, Campbell H, Rudan I, Johansson A, Marroni F, Hayward C, Vitart V, Jonasson I, Pattaro C, Wright A, Hastie N, Pichler I, Hicks AA, Falchi M, Willemsen G, Hottenga JJ, de Geus EJ, Montgomery GW, Whitfield J, Magnusson P, Saharinen J, Perola M, Silander K, Isaacs A, Sijbrands EJ, Uitterlinden AG, Witteman JC, Oostra BA, Elliott P, Ruokonen A, Sabatti C, Gieger C, Meitinger T, Kronenberg F, Döring A, Wichmann HE, Smit JH, McCarthy MI, van Duijn CM, Peltonen L and ENGAGE Consortium

    [1] Department of Epidemiology and Biostatistics, Erasmus University Medical Center, P.O. Box 2040, 3000 CA Rotterdam, The Netherlands. [2] These authors contributed equally to this work.

    Recent genome-wide association (GWA) studies of lipids have been conducted in samples ascertained for other phenotypes, particularly diabetes. Here we report the first GWA analysis of loci affecting total cholesterol (TC), low-density lipoprotein (LDL) cholesterol, high-density lipoprotein (HDL) cholesterol and triglycerides sampled randomly from 16 population-based cohorts and genotyped using mainly the Illumina HumanHap300-Duo platform. Our study included a total of 17,797-22,562 persons, aged 18-104 years and from geographic regions spanning from the Nordic countries to Southern Europe. We established 22 loci associated with serum lipid levels at a genome-wide significance level (P < 5 x 10(-8)), including 16 loci that were identified by previous GWA studies. The six newly identified loci in our cohort samples are ABCG5 (TC, P = 1.5 x 10(-11); LDL, P = 2.6 x 10(-10)), TMEM57 (TC, P = 5.4 x 10(-10)), CTCF-PRMT8 region (HDL, P = 8.3 x 10(-16)), DNAH11 (LDL, P = 6.1 x 10(-9)), FADS3-FADS2 (TC, P = 1.5 x 10(-10); LDL, P = 4.4 x 10(-13)) and MADD-FOLH1 region (HDL, P = 6 x 10(-11)). For three loci, effect sizes differed significantly by sex. Genetic risk scores based on lipid loci explain up to 4.8% of variation in lipids and were also associated with increased intima media thickness (P = 0.001) and coronary heart disease incidence (P = 0.04). The genetic risk score improves the screening of high-risk groups of dyslipidemia over classical risk factors.

    Funded by: Chief Scientist Office: CZB/4/710; Medical Research Council: MC_U127561128; NHLBI NIH HHS: 5R01HL087679-02, R01 HL087679; Wellcome Trust: 089061

    Nature genetics 2008;41;1;47-55

  • Dynamic nature of the proximal AZFc region of the human Y chromosome: multiple independent deletion and duplication events revealed by microsatellite analysis.

    Balaresque P, Bowden GR, Parkin EJ, Omran GA, Heyer E, Quintana-Murci L, Roewer L, Stoneking M, Nasidze I, Carvalho-Silva DR, Tyler-Smith C, de Knijff P and Jobling MA

    Department of Genetics, University of Leicester, Leicester, United Kingdom.

    The human Y chromosome shows frequent structural variants, some of which are selectively neutral, while others cause impaired fertility due to the loss of spermatogenic genes. The large-scale use of multiple Y-chromosomal microsatellites in forensic and population genetic studies can reveal such variants, through the absence or duplication of specific markers in haplotypes. We describe Y chromosomes in apparently normal males carrying null and duplicated alleles at the microsatellite DYS448, which lies in the proximal part of the azoospermia factor c (AZFc) region, important in spermatogenesis, and made up of "ampliconic" repeats that act as substrates for nonallelic homologous recombination (NAHR). Physical mapping in 26 DYS448 deletion chromosomes reveals that only three cases belong to a previously described class, representing independent occurrences of an approximately 1.5-Mb deletion mediated by recombination between the b1 and b3 repeat units. The remainder belong to five novel classes; none appears to be mediated through homologous recombination, and all remove some genes, but are likely to be compatible with normal fertility. A combination of deletion analysis with binary-marker and microsatellite haplotyping shows that the 26 deletions represent nine independent events. Nine DYS448 duplication chromosomes can be explained by four independent events. Some lineages have risen to high frequency in particular populations, in particular a deletion within haplogroup (hg) C(*)(xC3a,C3c) found in 18 Asian males. The nonrandom phylogenetic distribution of duplication and deletion events suggests possible structural predisposition to such mutations in hgs C and G.

    Funded by: Wellcome Trust: 057559, 077009

    Human mutation 2008;29;10;1171-80

  • Genomic complexity of the Y-STR DYS19: inversions, deletions and founder lineages carrying duplications.

    Balaresque P, Parkin EJ, Roewer L, Carvalho-Silva DR, Mitchell RJ, van Oorschot RA, Henke J, Stoneking M, Nasidze I, Wetton J, de Knijff P, Tyler-Smith C and Jobling MA

    Department of Genetics, University of Leicester, University Road, Leicester, LE1 7RH, UK.

    The Y-STR DYS19 is firmly established in the repertoire of Y-chromosomal markers used in forensic analysis yet is poorly understood at the molecular level, lying in a complex genomic environment and exhibiting null alleles, as well as duplications and occasional triplications in population samples. Here, we analyse three null alleles and 51 duplications and show that DYS19 can also be involved in inversion events, so that even its location within the short arm of the Y chromosome is uncertain. Deletion mapping in the three chromosomes carrying null alleles shows that their deletions are less than approximately 300 kb in size. Haplotypic analysis with binary markers shows that they belong to three different haplogroups and so represent independent events. In contrast, a collection of 51 DYS19 duplication chromosomes belong to only four haplogroups: two are singletons and may represent somatic mutation in lymphoblastoid cell lines, but two, in haplogroups G and C3c, represent founder lineages that have spread widely in Central Europe/West Asia and East Asia, respectively. Consideration of candidate mechanisms underlying both deletions and duplications provides no evidence for the involvement of non-allelic homologous recombination, and they are likely to represent sporadic events with low mutation rates. Understanding the basis and population distribution of these DYS19 alleles will aid in the utilisation and interpretation of profiles that contain them.

    Funded by: Wellcome Trust: 057559, 077009

    International journal of legal medicine 2008;123;1;15-23

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

    Barnes C, Plagnol V, Fitzgerald T, Redon R, Marchini J, Clayton D and Hurles ME

    Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, UK.

    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.

    Funded by: Wellcome Trust: 061860

    Nature genetics 2008;40;10;1245-52

  • Genome-wide association defines more than 30 distinct susceptibility loci for Crohn's disease.

    Barrett JC, Hansoul S, Nicolae DL, Cho JH, Duerr RH, Rioux JD, Brant SR, Silverberg MS, Taylor KD, Barmada MM, Bitton A, Dassopoulos T, Datta LW, Green T, Griffiths AM, Kistner EO, Murtha MT, Regueiro MD, Rotter JI, Schumm LP, Steinhart AH, Targan SR, Xavier RJ, NIDDK IBD Genetics Consortium, Libioulle C, Sandor C, Lathrop M, Belaiche J, Dewit O, Gut I, Heath S, Laukens D, Mni M, Rutgeerts P, Van Gossum A, Zelenika D, Franchimont D, Hugot JP, de Vos M, Vermeire S, Louis E, Belgian-French IBD Consortium, Wellcome Trust Case Control Consortium, Cardon LR, Anderson CA, Drummond H, Nimmo E, Ahmad T, Prescott NJ, Onnie CM, Fisher SA, Marchini J, Ghori J, Bumpstead S, Gwilliam R, Tremelling M, Deloukas P, Mansfield J, Jewell D, Satsangi J, Mathew CG, Parkes M, Georges M and Daly MJ

    Bioinformatics and Statistical Genetics, Wellcome Trust Centre for Human Genetics, University of Oxford, Roosevelt Drive, Oxford OX3 7BN, UK.

    Several risk factors for Crohn's disease have been identified in recent genome-wide association studies. To advance gene discovery further, we combined data from three studies on Crohn's disease (a total of 3,230 cases and 4,829 controls) and carried out replication in 3,664 independent cases with a mixture of population-based and family-based controls. The results strongly confirm 11 previously reported loci and provide genome-wide significant evidence for 21 additional loci, including the regions containing STAT3, JAK2, ICOSLG, CDKAL1 and ITLN1. The expanded molecular understanding of the basis of this disease offers promise for informed therapeutic development.

    Funded by: Chief Scientist Office: CZB/4/540; Medical Research Council: G0000934, G0600329, G0800759; NIAID NIH HHS: AI06277, R01 AI062773, R01 AI062773-02; NIDDK NIH HHS: DK064869, DK62413, DK62420, DK62422, DK62423, DK62429, DK62431, DK62432, P30 DK040561, P30 DK040561-13, P30 DK063491, P30 DK063491-019004, P30 DK063491-029004, P30 DK063491-039004, P30 DK063491-049004, P30 DK063491-05, R01 DK064869, R01 DK064869-04, U01 DK062413, U01 DK062413-06, U01 DK062420, U01 DK062420-01, U01 DK062420-02, U01 DK062420-03, U01 DK062420-04, U01 DK062420-05, U01 DK062420-06, U01 DK062422, U01 DK062422-07, U01 DK062423, U01 DK062423-06, U01 DK062429, U01 DK062429-07, U01 DK062431, U01 DK062431-06, U01 DK062432; Wellcome Trust: 068545/Z/02

    Nature genetics 2008;40;8;955-62

  • Population-specific risk of type 2 diabetes conferred by HNF4A P2 promoter variants: a lesson for replication studies.

    Barroso I, Luan J, Wheeler E, Whittaker P, Wasson J, Zeggini E, Weedon MN, Hunt S, Venkatesh R, Frayling TM, Delgado M, Neuman RJ, Zhao J, Sherva R, Glaser B, Walker M, Hitman G, McCarthy MI, Hattersley AT, Permutt MA, Wareham NJ and Deloukas P

    Wellcome Trust Sanger Institute, Hinxton, Cambridge, UK.

    Objective: Single nucleotide polymorphisms (SNPs) in the P2 promoter region of HNF4A were originally shown to be associated with predisposition for type 2 diabetes in Finnish, Ashkenazi, and, more recently, Scandinavian populations, but they generated conflicting results in additional populations. We aimed to investigate whether data from a large-scale mapping approach would replicate this association in novel Ashkenazi samples and in U.K. populations and whether these data would allow us to refine the association signal.

    Research design and methods: Using a dense linkage disequilibrium map of 20q, we selected SNPs from a 10-Mb interval centered on HNF4A. In a staged approach, we first typed 4,608 SNPs in case-control populations from four U.K. populations and an Ashkenazi population (n = 2,516). In phase 2, a subset of 763 SNPs was genotyped in 2,513 additional samples from the same populations.

    Results: Combined analysis of both phases demonstrated association between HNF4A P2 SNPs (rs1884613 and rs2144908) and type 2 diabetes in the Ashkenazim (n = 991; P < 1.6 x 10(-6)). Importantly, these associations are significant in a subset of Ashkenazi samples (n = 531) not previously tested for association with P2 SNPs (odds ratio [OR] approximately 1.7; P < 0.002), thus providing replication within the Ashkenazim. In the U.K. populations, this association was not significant (n = 4,022; P > 0.5), and the estimate for the OR was much smaller (OR 1.04; [95%CI 0.91-1.19]).

    Conclusions: These data indicate that the risk conferred by HNF4A P2 is significantly different between U.K. and Ashkenazi populations (P < 0.00007), suggesting that the underlying causal variant remains unidentified. Interactions with other genetic or environmental factors may also contribute to this difference in risk between populations.

    Funded by: Medical Research Council: MC_U106179471; NIDDK NIH HHS: R01 DK049583; PHS HHS: R01K049583; Wellcome Trust: 076113, 077016, 079557

    Diabetes 2008;57;11;3161-5

  • Re-evaluation of putative rheumatoid arthritis susceptibility genes in the post-genome wide association study era and hypothesis of a key pathway underlying susceptibility.

    Barton A, Thomson W, Ke X, Eyre S, Hinks A, Bowes J, Gibbons L, Plant D, Wellcome Trust Case Control Consortium, Wilson AG, Marinou I, Morgan A, Emery P, YEAR consortium, Steer S, Hocking L, Reid DM, Wordsworth P, Harrison P and Worthington J

    Arc-Epidemiology Unit, Stopford Building, The University of Manchester, Manchester, UK.

    Rheumatoid arthritis (RA) is an archetypal, common, complex autoimmune disease with both genetic and environmental contributions to disease aetiology. Two novel RA susceptibility loci have been reported from recent genome-wide and candidate gene association studies. We, therefore, investigated the evidence for association of the STAT4 and TRAF1/C5 loci with RA using imputed data from the Wellcome Trust Case Control Consortium (WTCCC). No evidence for association of variants mapping to the TRAF1/C5 gene was detected in the 1860 RA cases and 2930 control samples tested in that study. Variants mapping to the STAT4 gene did show evidence for association (rs7574865, P = 0.04). Given the association of the TRAF1/C5 locus in two previous large case-control series from populations of European descent and the evidence for association of the STAT4 locus in the WTCCC study, single nucleotide polymorphisms mapping to these loci were tested for association with RA in an independent UK series comprising DNA from >3000 cases with disease and >3000 controls and a combined analysis including the WTCCC data was undertaken. We confirm association of the STAT4 and the TRAF1/C5 loci with RA bringing to 5 the number of confirmed susceptibility loci. The effect sizes are less than those reported previously but are likely to be a more accurate reflection of the true effect size given the larger size of the cohort investigated in the current study.

    Funded by: Arthritis Research UK: 17552, 18475; Medical Research Council: G0000934; Wellcome Trust: 090532

    Human molecular genetics 2008;17;15;2274-9

  • Rheumatoid arthritis susceptibility loci at chromosomes 10p15, 12q13 and 22q13.

    Barton A, Thomson W, Ke X, Eyre S, Hinks A, Bowes J, Plant D, Gibbons LJ, Wellcome Trust Case Control Consortium, YEAR Consortium, BIRAC Consortium, Wilson AG, Bax DE, Morgan AW, Emery P, Steer S, Hocking L, Reid DM, Wordsworth P, Harrison P and Worthington J

    Arthritis Research Campaign, Epidemiology Unit, The University of Manchester, Manchester, UK.

    The WTCCC study identified 49 SNPs putatively associated with rheumatoid arthritis at P = 1 x 10(-4) - 1 x 10(-5) (tier 3). Here we show that three of these SNPs, mapping to chromosome 10p15 (rs4750316), 12q13 (rs1678542) and 22q13 (rs3218253), are also associated (trend P = 4 x 10(-5), P = 4 x 10(-4) and P = 4 x 10(-4), respectively) in a validation study of 4,106 individuals with rheumatoid arthritis and an expanded reference group of 11,238 subjects, confirming them as true susceptibility loci in individuals of European ancestry.

    Funded by: Arthritis Research UK: 17552, 18475; Chief Scientist Office: CZB/4/540; Medical Research Council: G0000934, G0600329; Wellcome Trust: 068545/Z/02, 090532

    Nature genetics 2008;40;10;1156-9

  • Phospholipid scramblases and Tubby-like proteins belong to a new superfamily of membrane tethered transcription factors.

    Bateman A, Finn RD, Sims PJ, Wiedmer T, Biegert A and Söding J

    Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, UK.

    Motivation: Phospholipid scramblases (PLSCRs) constitute a family of cytoplasmic membrane-associated proteins that were identified based upon their capacity to mediate a Ca(2+)-dependent bidirectional movement of phospholipids across membrane bilayers, thereby collapsing the normally asymmetric distribution of such lipids in cell membranes. The exact function and mechanism(s) of these proteins nevertheless remains obscure: data from several laboratories now suggest that in addition to their putative role in mediating transbilayer flip/flop of membrane lipids, the PLSCRs may also function to regulate diverse processes including signaling, apoptosis, cell proliferation and transcription. A major impediment to deducing the molecular details underlying the seemingly disparate biology of these proteins is the current absence of any representative molecular structures to provide guidance to the experimental investigation of their function.

    Results: Here, we show that the enigmatic PLSCR family of proteins is directly related to another family of cellular proteins with a known structure. The Arabidopsis protein At5g01750 from the DUF567 family was solved by X-ray crystallography and provides the first structural model for this family. This model identifies that the presumed C-terminal transmembrane helix is buried within the core of the PLSCR structure, suggesting that palmitoylation may represent the principal membrane anchorage for these proteins. The fold of the PLSCR family is also shared by Tubby-like proteins. A search of the PDB with the HHpred server suggests a common evolutionary ancestry. Common functional features also suggest that tubby and PLSCR share a functional origin as membrane tethered transcription factors with capacity to modulate phosphoinositide-based signaling.

    Funded by: NHLBI NIH HHS: HL036946, HL063819, HL076215; Wellcome Trust: 087656, WT077044/Z/05/Z

    Bioinformatics (Oxford, England) 2008;25;2;159-62

  • A novel 154-bp deletion in the human mitochondrial DNA control region in healthy individuals.

    Behar DM, Blue-Smith J, Soria-Hernanz DF, Tzur S, Hadid Y, Bormans C, Moen A, Tyler-Smith C, Quintana-Murci L, Wells RS and Genographic Consortium

    Molecular Medicine Laboratory, Rambam Health Care Campus, Haifa, Israel.

    The biological role of the mitochondrial DNA (mtDNA) control region in mtDNA replication remains unclear. In a worldwide survey of mtDNA variation in the general population, we have identified a novel large control region deletion spanning positions 16154 to 16307 (m.16154_16307del154). The population prevalence of this deletion is low, since it was only observed in 1 out of over 120,000 mtDNA genomes studied. The deletion is present in a nonheteroplasmic state, and was transmitted by a mother to her two sons with no apparent past or present disease conditions. The identification of this large deletion in healthy individuals challenges the current view of the control region as playing a crucial role in the regulation of mtDNA replication, and supports the existence of a more complex system of multiple or epigenetically-determined replication origins.

    Funded by: Wellcome Trust: 077009

    Human mutation 2008;29;12;1387-91

  • The dawn of human matrilineal diversity.

    Behar DM, Villems R, Soodyall H, Blue-Smith J, Pereira L, Metspalu E, Scozzari R, Makkan H, Tzur S, Comas D, Bertranpetit J, Quintana-Murci L, Tyler-Smith C, Wells RS, Rosset S and Genographic Consortium

    Molecular Medicine Laboratory, Rambam Health Care Campus, Haifa 31096, Israel.

    The quest to explain demographic history during the early part of human evolution has been limited because of the scarce paleoanthropological record from the Middle Stone Age. To shed light on the structure of the mitochondrial DNA (mtDNA) phylogeny at the dawn of Homo sapiens, we constructed a matrilineal tree composed of 624 complete mtDNA genomes from sub-Saharan Hg L lineages. We paid particular attention to the Khoi and San (Khoisan) people of South Africa because they are considered to be a unique relic of hunter-gatherer lifestyle and to carry paternal and maternal lineages belonging to the deepest clades known among modern humans. Both the tree phylogeny and coalescence calculations suggest that Khoisan matrilineal ancestry diverged from the rest of the human mtDNA pool 90,000-150,000 years before present (ybp) and that at least five additional, currently extant maternal lineages existed during this period in parallel. Furthermore, we estimate that a minimum of 40 other evolutionarily successful lineages flourished in sub-Saharan Africa during the period of modern human dispersal out of Africa approximately 60,000-70,000 ybp. Only much later, at the beginning of the Late Stone Age, about 40,000 ybp, did introgression of additional lineages occur into the Khoisan mtDNA pool. This process was further accelerated during the recent Bantu expansions. Our results suggest that the early settlement of humans in Africa was already matrilineally structured and involved small, separately evolving isolated populations.

    Funded by: Wellcome Trust

    American journal of human genetics 2008;82;5;1130-40

  • Genomic 'valleys of death'.

    Bentley S

    Nature reviews. Microbiology 2008;6;4;260-1

  • Genome of the actinomycete plant pathogen Clavibacter michiganensis subsp. sepedonicus suggests recent niche adaptation.

    Bentley SD, Corton C, Brown SE, Barron A, Clark L, Doggett J, Harris B, Ormond D, Quail MA, May G, Francis D, Knudson D, Parkhill J and Ishimaru CA

    Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton CB10 1SA, United Kingdom.

    Clavibacter michiganensis subsp. sepedonicus is a plant-pathogenic bacterium and the causative agent of bacterial ring rot, a devastating agricultural disease under strict quarantine control and zero tolerance in the seed potato industry. This organism appears to be largely restricted to an endophytic lifestyle, proliferating within plant tissues and unable to persist in the absence of plant material. Analysis of the genome sequence of C. michiganensis subsp. sepedonicus and comparison with the genome sequences of related plant pathogens revealed a dramatic recent evolutionary history. The genome contains 106 insertion sequence elements, which appear to have been active in extensive rearrangement of the chromosome compared to that of Clavibacter michiganensis subsp. michiganensis. There are 110 pseudogenes with overrepresentation in functions associated with carbohydrate metabolism, transcriptional regulation, and pathogenicity. Genome comparisons also indicated that there is substantial gene content diversity within the species, probably due to differential gene acquisition and loss. These genomic features and evolutionary dating suggest that there was recent adaptation for life in a restricted niche where nutrient diversity and perhaps competition are low, correlated with a reduced ability to exploit previously occupied complex niches outside the plant. Toleration of factors such as multiplication and integration of insertion sequence elements, genome rearrangements, and functional disruption of many genes and operons seems to indicate that there has been general relaxation of selective pressure on a large proportion of the genome.

    Journal of bacteriology 2008;190;6;2150-60

  • Susceptibility loci for intracranial aneurysm in European and Japanese populations.

    Bilguvar K, Yasuno K, Niemelä M, Ruigrok YM, von Und Zu Fraunberg M, van Duijn CM, van den Berg LH, Mane S, Mason CE, Choi M, Gaál E, Bayri Y, Kolb L, Arlier Z, Ravuri S, Ronkainen A, Tajima A, Laakso A, Hata A, Kasuya H, Koivisto T, Rinne J, Ohman J, Breteler MM, Wijmenga C, State MW, Rinkel GJ, Hernesniemi J, Jääskeläinen JE, Palotie A, Inoue I, Lifton RP and Günel M

    Department of Neurosurgery, Neurobiology, Yale Center for Human Genetics and Genomics, Yale University School of Medicine, New Haven, CT 06510, USA.

    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.

    Funded by: Howard Hughes Medical Institute; NINDS NIH HHS: R01 NS057756, U24 NS051869; Wellcome Trust: 089062

    Nature genetics 2008;40;12;1472-7

  • Comparison of Mascot and X!Tandem performance for low and high accuracy mass spectrometry and the development of an adjusted Mascot threshold.

    Brosch M, Swamy S, Hubbard T and Choudhary J

    The Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge CB10 1SA, United Kingdom.

    It is a major challenge to develop effective sequence database search algorithms to translate molecular weight and fragment mass information obtained from tandem mass spectrometry into high quality peptide and protein assignments. We investigated the peptide identification performance of Mascot and X!Tandem for mass tolerance settings common for low and high accuracy mass spectrometry. We demonstrated that sensitivity and specificity of peptide identification can vary substantially for different mass tolerance settings, but this effect was more significant for Mascot. We present an adjusted Mascot threshold, which allows the user to freely select the best trade-off between sensitivity and specificity. The adjusted Mascot threshold was compared with the default Mascot and X!Tandem scoring thresholds and shown to be more sensitive at the same false discovery rates for both low and high accuracy mass spectrometry data.

    Funded by: Wellcome Trust: 077198

    Molecular & cellular proteomics : MCP 2008;7;5;962-70

  • Identification of variation in the platelet transcriptome associated with glycoprotein 6 haplotype.

    Burns P, Gusnanto A, Macaulay IC, Rankin A, Tom B, Langford CF, Dudbridge F, Ouwehand WH, Watkins NA and Bloodomics Consortium

    Department of Haematology, University of Cambridge and National Health Service Blood and Transplant, Cambridge, UK.

    Platelet Glycoprotein VI (GPVI) is the activatory collagen signalling receptor that transmits an outside-in signal via the FcR gamma-chain. In Caucasians two GP6 haplotypes have been identified which encode GPVI isoforms that differ by five amino-acids. The minor haplotype is associated with a modest but statistically significant reduction in GPVI abundance and reduced downstream signalling events. As GPVI is also expressed on megakaryocytes, different GPVI isoforms may imprint on the platelet transcriptome. We investigated the association of GP6 haplotype with transcription by comparing the transcriptomes of platelets from individuals homozygous for the major ('a') and minor ('b') haplotypes to identify differentially expressed (DE) transcripts. Platelet RNA was isolated from apheresis concentrates from 16 'aa' donors and eight 'bb' donors. mRNA was amplified using a template-switching PCR based protocol and fluorescently labelled. Samples were randomly paired both within and between haplotypes and compared on a cDNA microarray. No consistently DE transcripts were identified within the 'aa' haplotype but 52 significantly DE transcripts were observed between haplotypes. Generally the fold differences were low (two to four-fold) but were confirmed by qRT-PCR for selected transcripts (TUBB1, P = 0.0004; VWF, P = 0.0126). The results of this study indicate that there are subtle differences between the platelet transcriptomes of individuals who differ by GP6 haplotype. The identification of DE genes may identify critical pathways and nodes not previously known to be involved in platelet development and function.

    Funded by: Medical Research Council: MC_U105261167; Wellcome Trust

    Platelets 2008;19;4;258-67

  • Large-scale screening for novel low-affinity extracellular protein interactions.

    Bushell KM, Söllner C, Schuster-Boeckler B, Bateman A and Wright GJ

    Cell Surface Signalling Laboratory, Wellcome Trust Sanger Institute, Hinxton, Cambridge CB10 1HH, United Kingdom.

    Extracellular protein-protein interactions are essential for both intercellular communication and cohesion within multicellular organisms. Approximately a fifth of human genes encode membrane-tethered or secreted proteins, but they are largely absent from recent large-scale protein interaction datasets, making current interaction networks biased and incomplete. This discrepancy is due to the unsuitability of popular high-throughput methods to detect extracellular interactions because of the biochemical intractability of membrane proteins and their interactions. For example, cell surface proteins contain insoluble hydrophobic transmembrane regions, and their extracellular interactions are often highly transient, having half-lives of less than a second. To detect transient extracellular interactions on a large scale, we developed AVEXIS (avidity-based extracellular interaction screen), a high-throughput assay that overcomes these technical issues and can detect very transient interactions (half-lives <or= 0.1 sec) with a low false-positive rate. We used it to systematically screen for receptor-ligand pairs within the zebrafish immunoglobulin superfamily and identified novel ligands for both well-known and orphan receptors. Genes encoding receptor-ligand pairs were often clustered phylogenetically and expressed in the same or adjacent tissues, immediately implying their involvement in similar biological processes. Using AVEXIS, we have determined the first systematic low-affinity extracellular protein interaction network, supported by independent biological data. This technique will now allow large-scale extracellular protein interaction mapping in a broad range of experimental contexts.

    Funded by: Wellcome Trust: 087656

    Genome research 2008;18;4;622-30

  • Subclonal phylogenetic structures in cancer revealed by ultra-deep sequencing.

    Campbell PJ, Pleasance ED, Stephens PJ, Dicks E, Rance R, Goodhead I, Follows GA, Green AR, Futreal PA and Stratton MR

    Wellcome Trust Sanger Institute, Hinxton CB10 1SA, United Kingdom.

    During the clonal expansion of cancer from an ancestral cell with an initiating oncogenic mutation to symptomatic neoplasm, the occurrence of somatic mutations (both driver and passenger) can be used to track the on-going evolution of the neoplasm. All subclones within a cancer are phylogenetically related, with the prevalence of each subclone determined by its evolutionary fitness and the timing of its origin relative to other subclones. Recently developed massively parallel sequencing platforms promise the ability to detect rare subclones of genetic variants without a priori knowledge of the mutations involved. We used ultra-deep pyrosequencing to investigate intraclonal diversification at the Ig heavy chain locus in 22 patients with B-cell chronic lymphocytic leukemia. Analysis of a non-polymorphic control locus revealed artifactual insertions and deletions resulting from sequencing errors and base substitutions caused by polymerase misincorporation during PCR amplification. We developed an algorithm to differentiate genuine haplotypes of somatic hypermutations from such artifacts. This proved capable of detecting multiple rare subclones with frequencies as low as 1 in 5000 copies and allowed the characterization of phylogenetic interrelationships among subclones within each patient. This study demonstrates the potential for ultra-deep resequencing to recapitulate the dynamics of clonal evolution in cancer cell populations.

    Funded by: Wellcome Trust: 088340

    Proceedings of the National Academy of Sciences of the United States of America 2008;105;35;13081-6

  • Identification of somatically acquired rearrangements in cancer using genome-wide massively parallel paired-end sequencing.

    Campbell PJ, Stephens PJ, Pleasance ED, O'Meara S, Li H, Santarius T, Stebbings LA, Leroy C, Edkins S, Hardy C, Teague JW, Menzies A, Goodhead I, Turner DJ, Clee CM, Quail MA, Cox A, Brown C, Durbin R, Hurles ME, Edwards PA, Bignell GR, Stratton MR and Futreal PA

    Wellcome Trust Sanger Institute, Hinxton CB10 1SA, UK.

    Human cancers often carry many somatically acquired genomic rearrangements, some of which may be implicated in cancer development. However, conventional strategies for characterizing rearrangements are laborious and low-throughput and have low sensitivity or poor resolution. We used massively parallel sequencing to generate sequence reads from both ends of short DNA fragments derived from the genomes of two individuals with lung cancer. By investigating read pairs that did not align correctly with respect to each other on the reference human genome, we characterized 306 germline structural variants and 103 somatic rearrangements to the base-pair level of resolution. The patterns of germline and somatic rearrangement were markedly different. Many somatic rearrangements were from amplicons, although rearrangements outside these regions, notably including tandem duplications, were also observed. Some somatic rearrangements led to abnormal transcripts, including two from internal tandem duplications and two fusion transcripts created by interchromosomal rearrangements. Germline variants were predominantly mediated by retrotransposition, often involving AluY and LINE elements. The results demonstrate the feasibility of systematic, genome-wide characterization of rearrangements in complex human cancer genomes, raising the prospect of a new harvest of genes associated with cancer using this strategy.

    Funded by: Wellcome Trust: 077012, 088340

    Nature genetics 2008;40;6;722-9

  • Artemis and ACT: viewing, annotating and comparing sequences stored in a relational database.

    Carver T, Berriman M, Tivey A, Patel C, Böhme U, Barrell BG, Parkhill J and Rajandream MA

    Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge CB10 1SA, UK.

    Motivation: Artemis and Artemis Comparison Tool (ACT) have become mainstream tools for viewing and annotating sequence data, particularly for microbial genomes. Since its first release, Artemis has been continuously developed and supported with additional functionality for editing and analysing sequences based on feedback from an active user community of laboratory biologists and professional annotators. Nevertheless, its utility has been somewhat restricted by its limitation to reading and writing from flat files. Therefore, a new version of Artemis has been developed, which reads from and writes to a relational database schema, and allows users to annotate more complex, often large and fragmented, genome sequences.

    Results: Artemis and ACT have now been extended to read and write directly to the Generic Model Organism Database (GMOD, Chado relational database schema. In addition, a Gene Builder tool has been developed to provide structured forms and tables to edit coordinates of gene models and edit functional annotation, based on standard ontologies, controlled vocabularies and free text.

    Availability: Artemis and ACT are freely available (under a GPL licence) for download (for MacOSX, UNIX and Windows) at the Wellcome Trust Sanger Institute web sites:

    Funded by: Wellcome Trust: 082372

    Bioinformatics (Oxford, England) 2008;24;23;2672-6

  • Priorities for nucleotide trace, sequence and annotation data capture at the Ensembl Trace Archive and the EMBL Nucleotide Sequence Database.

    Cochrane G, Akhtar R, Aldebert P, Althorpe N, Baldwin A, Bates K, Bhattacharyya S, Bonfield J, Bower L, Browne P, Castro M, Cox T, Demiralp F, Eberhardt R, Faruque N, Hoad G, Jang M, Kulikova T, Labarga A, Leinonen R, Leonard S, Lin Q, Lopez R, Lorenc D, McWilliam H, Mukherjee G, Nardone F, Plaister S, Robinson S, Sobhany S, Vaughan R, Wu D, Zhu W, Apweiler R, Hubbard T and Birney E

    EMBL-European Bioinformatics Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, UK.

    The Ensembl Trace Archive ( and the EMBL Nucleotide Sequence Database (, known together as the European Nucleotide Archive, continue to see growth in data volume and diversity. Selected major developments of 2007 are presented briefly, along with data submission and retrieval information. In the face of increasing requirements for nucleotide trace, sequence and annotation data archiving, data capture priority decisions have been taken at the European Nucleotide Archive. Priorities are discussed in terms of how reliably information can be captured, the long-term benefits of its capture and the ease with which it can be captured.

    Funded by: Wellcome Trust: 062023, 077198, 085532

    Nucleic acids research 2008;36;Database issue;D5-12

  • Identifying protein domains with the Pfam database.

    Coggill P, Finn RD and Bateman A

    Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, United Kingdom.

    Pfam is a database of protein domain families, with each family represented by multiple sequence alignments and profile hidden Markov models (HMMs). In addition, each family has associated annotation, literature references, and links to other databases. The entries in Pfam are available via the World Wide Web and in flatfile format. This unit contains detailed information on how to access and utilize the information present in the Pfam database, namely the families, multiple alignments, and annotation. Details on running Pfam, both remotely and locally are presented.

    Funded by: Wellcome Trust: 087656

    Current protocols in bioinformatics / editoral board, Andreas D. Baxevanis ... [et al.] 2008;Chapter 2;Unit 2.5

  • nGASP--the nematode genome annotation assessment project.

    Coghlan A, Fiedler TJ, McKay SJ, Flicek P, Harris TW, Blasiar D, nGASP Consortium and Stein LD

    Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SA, UK.

    Background: While the C. elegans genome is extensively annotated, relatively little information is available for other Caenorhabditis species. The nematode genome annotation assessment project (nGASP) was launched to objectively assess the accuracy of protein-coding gene prediction software in C. elegans, and to apply this knowledge to the annotation of the genomes of four additional Caenorhabditis species and other nematodes. Seventeen groups worldwide participated in nGASP, and submitted 47 prediction sets across 10 Mb of the C. elegans genome. Predictions were compared to reference gene sets consisting of confirmed or manually curated gene models from WormBase.

    Results: The most accurate gene-finders were 'combiner' algorithms, which made use of transcript- and protein-alignments and multi-genome alignments, as well as gene predictions from other gene-finders. Gene-finders that used alignments of ESTs, mRNAs and proteins came in second. There was a tie for third place between gene-finders that used multi-genome alignments and ab initio gene-finders. The median gene level sensitivity of combiners was 78% and their specificity was 42%, which is nearly the same accuracy reported for combiners in the human genome. C. elegans genes with exons of unusual hexamer content, as well as those with unusually many exons, short exons, long introns, a weak translation start signal, weak splice sites, or poorly conserved orthologs posed the greatest difficulty for gene-finders.

    Conclusion: This experiment establishes a baseline of gene prediction accuracy in Caenorhabditis genomes, and has guided the choice of gene-finders for the annotation of newly sequenced genomes of Caenorhabditis and other nematode species. We have created new gene sets for C. briggsae, C. remanei, C. brenneri, C. japonica, and Brugia malayi using some of the best-performing gene-finders.

    Funded by: NHGRI NIH HHS: P41 HG02223; Wellcome Trust

    BMC bioinformatics 2008;9;549

  • Mapping multiprotein complexes by affinity purification and mass spectrometry.

    Collins MO and Choudhary JS

    Proteomic Mass Spectrometry, The Wellcome Trust Sanger Institute, Hinxton, Cambridge CB10 1SA, UK.

    The combination of affinity purification and tandem mass spectrometry (MS) has emerged as a powerful approach to delineate biological processes. In particular, the use of epitope tags has allowed this approach to become scaleable and has bypassed difficulties associated with generation of antibodies. Single epitope tags and tandem affinity purification (TAP) tags have been used to systematically map protein complexes generating protein interaction data at a near proteome-wide scale. Recent developments in the design of tags, optimisation of purification conditions, experimental design and data analysis have greatly improved the sensitivity and specificity of this approach. Concomitant developments in MS, including high accuracy and high-throughput instrumentation together with quantitative MS methods, have facilitated large-scale and comprehensive analysis of multiprotein complexes.

    Current opinion in biotechnology 2008;19;4;324-30

  • Phosphoproteomic analysis of the mouse brain cytosol reveals a predominance of protein phosphorylation in regions of intrinsic sequence disorder.

    Collins MO, Yu L, Campuzano I, Grant SG and Choudhary JS

    Proteomic Mass Spectrometry, The Wellcome Trust Sanger Institute, Hinxton, Cambridge CB101SA, United Kingdom.

    We analyzed the mouse forebrain cytosolic phosphoproteome using sequential (protein and peptide) IMAC purifications, enzymatic dephosphorylation, and targeted tandem mass spectrometry analysis strategies. In total, using complementary phosphoenrichment and LC-MS/MS strategies, 512 phosphorylation sites on 540 non-redundant phosphopeptides from 162 cytosolic phosphoproteins were characterized. Analysis of protein domains and amino acid sequence composition of this data set of cytosolic phosphoproteins revealed that it is significantly enriched in intrinsic sequence disorder, and this enrichment is associated with both cellular location and phosphorylation status. The majority of phosphorylation sites found by MS were located outside of structural protein domains (97%) but were mostly located in regions of intrinsic sequence disorder (86%). 368 phosphorylation sites were located in long regions of disorder (over 40 amino acids long), and 94% of proteins contained at least one such long region of disorder. In addition, we found that 58 phosphorylation sites in this data set occur in 14-3-3 binding consensus motifs, linear motifs that are associated with unstructured regions in proteins. These results demonstrate that in this data set protein phosphorylation is significantly depleted in protein domains and significantly enriched in disordered protein sequences and that enrichment of intrinsic sequence disorder may be a common feature of phosphoproteomes. This supports the hypothesis that disordered regions in proteins allow kinases, phosphatases, and phosphorylation-dependent binding proteins to gain access to target sequences to regulate local protein conformation and activity.

    Funded by: Wellcome Trust

    Molecular & cellular proteomics : MCP 2008;7;7;1331-48

  • Strong genetic evidence for a selective influence of GABAA receptors on a component of the bipolar disorder phenotype.

    Craddock N, Jones L, Jones IR, Kirov G, Green EK, Grozeva D, Moskvina V, Nikolov I, Hamshere ML, Vukcevic D, Caesar S, Gordon-Smith K, Fraser C, Russell E, Norton N, Breen G, St Clair D, Collier DA, Young AH, Ferrier IN, Farmer A, McGuffin P, Holmans PA, Wellcome Trust Case Control Consortium (WTCCC), Donnelly P, Owen MJ and O'Donovan MC

    Department of Psychological Medicine, School of Medicine, Cardiff University, Cardiff, UK.

    Despite compelling evidence for a major genetic contribution to risk of bipolar mood disorder, conclusive evidence implicating specific genes or pathophysiological systems has proved elusive. In part this is likely to be related to the unknown validity of current phenotype definitions and consequent aetiological heterogeneity of samples. In the recent Wellcome Trust Case Control Consortium genome-wide association analysis of bipolar disorder (1868 cases, 2938 controls) one of the most strongly associated polymorphisms lay within the gene encoding the GABA(A) receptor beta1 subunit, GABRB1. Aiming to increase biological homogeneity, we sought the diagnostic subset that showed the strongest signal at this polymorphism and used this to test for independent evidence of association with other members of the GABA(A) receptor gene family. The index signal was significantly enriched in the 279 cases meeting Research Diagnostic Criteria for schizoaffective disorder, bipolar type (P=3.8 x 10(-6)). Independently, these cases showed strong evidence that variation in GABA(A) receptor genes influences risk for this phenotype (independent system-wide P=6.6 x 10(-5)) with association signals also at GABRA4, GABRB3, GABRA5 and GABRR3. [corrected] Our findings have the potential to inform understanding of presentation, pathogenesis and nosology of bipolar disorders. Our method of phenotype refinement may be useful in studies of other complex psychiatric and non-psychiatric disorders.

    Funded by: Medical Research Council: G0000934, G0801418; Wellcome Trust: 079643

    Molecular psychiatry 2008;15;2;146-53

  • The complete genome, comparative and functional analysis of Stenotrophomonas maltophilia reveals an organism heavily shielded by drug resistance determinants.

    Crossman LC, Gould VC, Dow JM, Vernikos GS, Okazaki A, Sebaihia M, Saunders D, Arrowsmith C, Carver T, Peters N, Adlem E, Kerhornou A, Lord A, Murphy L, Seeger K, Squares R, Rutter S, Quail MA, Rajandream MA, Harris D, Churcher C, Bentley SD, Parkhill J, Thomson NR and Avison MB

    Pathogen Sequencing Unit, The Wellcome Trust Sanger Institute, Hinxton, Cambridge, CB10 1SA, UK.

    Background: Stenotrophomonas maltophilia is a nosocomial opportunistic pathogen of the Xanthomonadaceae. The organism has been isolated from both clinical and soil environments in addition to the sputum of cystic fibrosis patients and the immunocompromised. Whilst relatively distant phylogenetically, the closest sequenced relatives of S. maltophilia are the plant pathogenic xanthomonads.

    Results: The genome of the bacteremia-associated isolate S. maltophilia K279a is 4,851,126 bp and of high G+C content. The sequence reveals an organism with a remarkable capacity for drug and heavy metal resistance. In addition to a number of genes conferring resistance to antimicrobial drugs of different classes via alternative mechanisms, nine resistance-nodulation-division (RND)-type putative antimicrobial efflux systems are present. Functional genomic analysis confirms a role in drug resistance for several of the novel RND efflux pumps. S. maltophilia possesses potentially mobile regions of DNA and encodes a number of pili and fimbriae likely to be involved in adhesion and biofilm formation that may also contribute to increased antimicrobial drug resistance.

    Conclusion: The panoply of antimicrobial drug resistance genes and mobile genetic elements found suggests that the organism can act as a reservoir of antimicrobial drug resistance determinants in a clinical environment, which is an issue of considerable concern.

    Funded by: Wellcome Trust

    Genome biology 2008;9;4;R74

  • The RNA WikiProject: community annotation of RNA families.

    Daub J, Gardner PP, Tate J, Ramsköld D, Manske M, Scott WG, Weinberg Z, Griffiths-Jones S and Bateman A

    The online encyclopedia Wikipedia has become one of the most important online references in the world and has a substantial and growing scientific content. A search of Google with many RNA-related keywords identifies a Wikipedia article as the top hit. We believe that the RNA community has an important and timely opportunity to maximize the content and quality of RNA information in Wikipedia. To this end, we have formed the RNA WikiProject ( as part of the larger Molecular and Cellular Biology WikiProject. We have created over 600 new Wikipedia articles describing families of noncoding RNAs based on the Rfam database, and invite the community to update, edit, and correct these articles. The Rfam database now redistributes this Wikipedia content as the primary textual annotation of its RNA families. Users can, therefore, for the first time, directly edit the content of one of the major RNA databases. We believe that this Wikipedia/Rfam link acts as a functioning model for incorporating community annotation into molecular biology databases.

    Funded by: Howard Hughes Medical Institute; NIGMS NIH HHS: R01 GM087721; Wellcome Trust: 077044

    RNA (New York, N.Y.) 2008;14;12;2462-4

  • From gene expression to disease risk.

    Dermitzakis ET

    Nature genetics 2008;40;5;492-3

  • Minimum information specification for in situ hybridization and immunohistochemistry experiments (MISFISHIE).

    Deutsch EW, Ball CA, Berman JJ, Bova GS, Brazma A, Bumgarner RE, Campbell D, Causton HC, Christiansen JH, Daian F, Dauga D, Davidson DR, Gimenez G, Goo YA, Grimmond S, Henrich T, Herrmann BG, Johnson MH, Korb M, Mills JC, Oudes AJ, Parkinson HE, Pascal LE, Pollet N, Quackenbush J, Ramialison M, Ringwald M, Salgado D, Sansone SA, Sherlock G, Stoeckert CJ, Swedlow J, Taylor RC, Walashek L, Warford A, Wilkinson DG, Zhou Y, Zon LI, Liu AY and True LD

    Institute for Systems Biology, 1441 N 34th Street, Seattle, Washington 98103, USA.

    One purpose of the biomedical literature is to report results in sufficient detail that the methods of data collection and analysis can be independently replicated and verified. Here we present reporting guidelines for gene expression localization experiments: the minimum information specification for in situ hybridization and immunohistochemistry experiments (MISFISHIE). MISFISHIE is modeled after the Minimum Information About a Microarray Experiment (MIAME) specification for microarray experiments. Both guidelines define what information should be reported without dictating a format for encoding that information. MISFISHIE describes six types of information to be provided for each experiment: experimental design, biomaterials and treatments, reporters, staining, imaging data and image characterizations. This specification has benefited the consortium within which it was developed and is expected to benefit the wider research community. We welcome feedback from the scientific community to help improve our proposal.

    Funded by: Biotechnology and Biological Sciences Research Council: BB/E025080/1; Medical Research Council: MC_U117532048, MC_U127527203; NHLBI NIH HHS: R33 HL073712; NIDDK NIH HHS: DK63328, DK63400, DK63481, DK63483, DK63630, R01 DK079798, R01 DK079798-01A2

    Nature biotechnology 2008;26;3;305-12

  • X-linked protocadherin 19 mutations cause female-limited epilepsy and cognitive impairment.

    Dibbens LM, Tarpey PS, Hynes K, Bayly MA, Scheffer IE, Smith R, Bomar J, Sutton E, Vandeleur L, Shoubridge C, Edkins S, Turner SJ, Stevens C, O'Meara S, Tofts C, Barthorpe S, Buck G, Cole J, Halliday K, Jones D, Lee R, Madison M, Mironenko T, Varian J, West S, Widaa S, Wray P, Teague J, Dicks E, Butler A, Menzies A, Jenkinson A, Shepherd R, Gusella JF, Afawi Z, Mazarib A, Neufeld MY, Kivity S, Lev D, Lerman-Sagie T, Korczyn AD, Derry CP, Sutherland GR, Friend K, Shaw M, Corbett M, Kim HG, Geschwind DH, Thomas P, Haan E, Ryan S, McKee S, Berkovic SF, Futreal PA, Stratton MR, Mulley JC and Gécz J

    Department of Genetic Medicine, Level 9 Rieger Building, Women's and Children's Hospital, 72 King William Road, North Adelaide, South Australia 5006, Australia.

    Epilepsy and mental retardation limited to females (EFMR) is a disorder with an X-linked mode of inheritance and an unusual expression pattern. Disorders arising from mutations on the X chromosome are typically characterized by affected males and unaffected carrier females. In contrast, EFMR spares transmitting males and affects only carrier females. Aided by systematic resequencing of 737 X chromosome genes, we identified different protocadherin 19 (PCDH19) gene mutations in seven families with EFMR. Five mutations resulted in the introduction of a premature termination codon. Study of two of these demonstrated nonsense-mediated decay of PCDH19 mRNA. The two missense mutations were predicted to affect adhesiveness of PCDH19 through impaired calcium binding. PCDH19 is expressed in developing brains of human and mouse and is the first member of the cadherin superfamily to be directly implicated in epilepsy or mental retardation.

    Funded by: NICHD NIH HHS: N01-HD-4-3368, N01-HD-4-3383; NIGMS NIH HHS: GM061354, P01 GM061354; NIMH NIH HHS: R01 MH 64547, R01 MH064547, R01 MH064547-01, R01 MH064547-01S1, R01 MH064547-02, R01 MH064547-02S1, R01 MH064547-03, R01 MH064547-04, R01 MH064547-05; Wellcome Trust

    Nature genetics 2008;40;6;776-81

  • Efficient targeted transcript discovery via array-based normalization of RACE libraries.

    Djebali S, Kapranov P, Foissac S, Lagarde J, Reymond A, Ucla C, Wyss C, Drenkow J, Dumais E, Murray RR, Lin C, Szeto D, Denoeud F, Calvo M, Frankish A, Harrow J, Makrythanasis P, Vidal M, Salehi-Ashtiani K, Antonarakis SE, Gingeras TR and Guigó R

    Grup de Recerca en Informàtica Biomèdica, Institut Municipal d'Investigació Mèdica/Universitat Pompeu Fabra, Dr. Aiguader 88, 08003 Barcelona, Spain.

    Rapid amplification of cDNA ends (RACE) is a widely used approach for transcript identification. Random clone selection from the RACE mixture, however, is an ineffective sampling strategy if the dynamic range of transcript abundances is large. To improve sampling efficiency of human transcripts, we hybridized the products of the RACE reaction onto tiling arrays and used the detected exons to delineate a series of reverse-transcriptase (RT)-PCRs, through which the original RACE transcript population was segregated into simpler transcript populations. We independently cloned the products and sequenced randomly selected clones. This approach, RACEarray, is superior to direct cloning and sequencing of RACE products because it specifically targets new transcripts and often results in overall normalization of transcript abundance. We show theoretically and experimentally that this strategy leads indeed to efficient sampling of new transcripts, and we investigated multiplexing the strategy by pooling RACE reactions from multiple interrogated loci before hybridization.

    Funded by: NCI NIH HHS: N01-CO-12400, N01CO12400; NHGRI NIH HHS: U01 HG003147, U01 HG003147-01, U01 HG003147-02, U01 HG003147-02S1, U01 HG003147-02S2, U01 HG003147-02S3, U01 HG003150, U01 HG003150-01, U01 HG003150-02, U01 HG003150-03, U01 HG003150-03S1, U01 HG003150-03S2, U01HG003147, U01HG003150, U54 HG004555, U54 HG004557, U54 HG004557-01, U54 HG004557-02, U54 HG004557-02S1, U54 HG004557-03; Wellcome Trust: 077198

    Nature methods 2008;5;7;629-35

  • NestedMICA as an ab initio protein motif discovery tool.

    Doğruel M, Down TA and Hubbard TJ

    Wellcome Trust Sanger Institute, Hinxton, Cambridge CB10 1HH, UK.

    Background: Discovering overrepresented patterns in amino acid sequences is an important step in protein functional element identification. We adapted and extended NestedMICA, an ab initio motif finder originally developed for finding transcription binding site motifs, to find short protein signals, and compared its performance with another popular protein motif finder, MEME. NestedMICA, an open source protein motif discovery tool written in Java, is driven by a Monte Carlo technique called Nested Sampling. It uses multi-class sequence background models to represent different "uninteresting" parts of sequences that do not contain motifs of interest. In order to assess NestedMICA as a protein motif finder, we have tested it on synthetic datasets produced by spiking instances of known motifs into a randomly selected set of protein sequences. NestedMICA was also tested using a biologically-authentic test set, where we evaluated its performance with respect to varying sequence length.

    Results: Generally NestedMICA recovered most of the short (3-9 amino acid long) test protein motifs spiked into a test set of sequences at different frequencies. We showed that it can be used to find multiple motifs at the same time, too. In all the assessment experiments we carried out, its overall motif discovery performance was better than that of MEME.

    Conclusion: NestedMICA proved itself to be a robust and sensitive ab initio protein motif finder, even for relatively short motifs that exist in only a small fraction of sequences.

    Availability: NestedMICA is available under the Lesser GPL open-source license from:

    Funded by: Wellcome Trust: 077198

    BMC bioinformatics 2008;9;19

  • A Bayesian deconvolution strategy for immunoprecipitation-based DNA methylome analysis.

    Down TA, Rakyan VK, Turner DJ, Flicek P, Li H, Kulesha E, Gräf S, Johnson N, Herrero J, Tomazou EM, Thorne NP, Bäckdahl L, Herberth M, Howe KL, Jackson DK, Miretti MM, Marioni JC, Birney E, Hubbard TJ, Durbin R, Tavaré S and Beck S

    Wellcome Trust Cancer Research UK Gurdon Institute, and Department of Genetics, University of Cambridge, Tennis Court Road, Cambridge CB2 1QR, UK.

    DNA methylation is an indispensible epigenetic modification required for regulating the expression of mammalian genomes. Immunoprecipitation-based methods for DNA methylome analysis are rapidly shifting the bottleneck in this field from data generation to data analysis, necessitating the development of better analytical tools. In particular, an inability to estimate absolute methylation levels remains a major analytical difficulty associated with immunoprecipitation-based DNA methylation profiling. To address this issue, we developed a cross-platform algorithm-Bayesian tool for methylation analysis (Batman)-for analyzing methylated DNA immunoprecipitation (MeDIP) profiles generated using oligonucleotide arrays (MeDIP-chip) or next-generation sequencing (MeDIP-seq). We developed the latter approach to provide a high-resolution whole-genome DNA methylation profile (DNA methylome) of a mammalian genome. Strong correlation of our data, obtained using mature human spermatozoa, with those obtained using bisulfite sequencing suggest that combining MeDIP-seq or MeDIP-chip with Batman provides a robust, quantitative and cost-effective functional genomic strategy for elucidating the function of DNA methylation.

    Funded by: Cancer Research UK: C14303/A8646; Wellcome Trust: 077198, 083563, 084071

    Nature biotechnology 2008;26;7;779-85

  • The evolution of the DLK1-DIO3 imprinted domain in mammals.

    Edwards CA, Mungall AJ, Matthews L, Ryder E, Gray DJ, Pask AJ, Shaw G, Graves JA, Rogers J, SAVOIR consortium, Dunham I, Renfree MB and Ferguson-Smith AC

    Department of Physiology, Development, and Neuroscience, University of Cambridge, Cambridge, United Kingdom.

    A comprehensive, domain-wide comparative analysis of genomic imprinting between mammals that imprint and those that do not can provide valuable information about how and why imprinting evolved. The imprinting status, DNA methylation, and genomic landscape of the Dlk1-Dio3 cluster were determined in eutherian, metatherian, and prototherian mammals including tammar wallaby and platypus. Imprinting across the whole domain evolved after the divergence of eutherian from marsupial mammals and in eutherians is under strong purifying selection. The marsupial locus at 1.6 megabases, is double that of eutherians due to the accumulation of LINE repeats. Comparative sequence analysis of the domain in seven vertebrates determined evolutionary conserved regions common to particular sub-groups and to all vertebrates. The emergence of Dlk1-Dio3 imprinting in eutherians has occurred on the maternally inherited chromosome and is associated with region-specific resistance to expansion by repetitive elements and the local introduction of noncoding transcripts including microRNAs and C/D small nucleolar RNAs. A recent mammal-specific retrotransposition event led to the formation of a completely new gene only in the eutherian domain, which may have driven imprinting at the cluster.

    Funded by: Medical Research Council: G0400156

    PLoS biology 2008;6;6;e135

  • Evolutionary expansion and anatomical specialization of synapse proteome complexity.

    Emes RD, Pocklington AJ, Anderson CN, Bayes A, Collins MO, Vickers CA, Croning MD, Malik BR, Choudhary JS, Armstrong JD and Grant SG

    Institute for Science and Technology in Medicine, Keele University, Thornburrow Drive, Hartshill, Stoke-on-Trent ST4 7QB, UK.

    Understanding the origins and evolution of synapses may provide insight into species diversity and the organization of the brain. Using comparative proteomics and genomics, we examined the evolution of the postsynaptic density (PSD) and membrane-associated guanylate kinase (MAGUK)-associated signaling complexes (MASCs) that underlie learning and memory. PSD and MASC orthologs found in yeast carry out basic cellular functions to regulate protein synthesis and structural plasticity. We observed marked changes in signaling complexity at the yeast-metazoan and invertebrate-vertebrate boundaries, with an expansion of key synaptic components, notably receptors, adhesion/cytoskeletal proteins and scaffold proteins. A proteomic comparison of Drosophila and mouse MASCs revealed species-specific adaptation with greater signaling complexity in mouse. Although synaptic components were conserved amongst diverse vertebrate species, mapping mRNA and protein expression in the mouse brain showed that vertebrate-specific components preferentially contributed to differences between brain regions. We propose that the evolution of synapse complexity around a core proto-synapse has contributed to invertebrate-vertebrate differences and to brain specialization.

    Funded by: Medical Research Council: G90/112, G90/93; Wellcome Trust: 077155

    Nature neuroscience 2008;11;7;799-806

  • Evaluating the role of LPIN1 variation in insulin resistance, body weight, and human lipodystrophy in U.K. Populations.

    Fawcett KA, Grimsey N, Loos RJ, Wheeler E, Daly A, Soos M, Semple R, Syddall H, Cooper C, Siniossoglou S, O'Rahilly S, Wareham NJ and Barroso I

    Metabolic Disease Group, Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridgeshire, U.K.

    Objective: Loss of lipin 1 activity causes lipodystrophy and insulin resistance in the fld mouse, and LPIN1 expression and common genetic variation were recently suggested to influence adiposity and insulin sensitivity in humans. We aimed to conduct a comprehensive association study to clarify the influence of common LPIN1 variation on adiposity and insulin sensitivity in U.K. populations and to examine the role of LPIN1 mutations in insulin resistance syndromes.

    Research design and method: Twenty-two single nucleotide polymorphisms tagging common LPIN1 variation were genotyped in Medical Research Council (MRC) Ely (n = 1,709) and Hertfordshire (n = 2,901) population-based cohorts. LPIN1 exons, exon/intron boundaries, and 3' untranslated region were sequenced in 158 patients with idiopathic severe insulin resistance (including 23 lipodystrophic patients) and 48 control subjects.

    Results: We found no association between LPIN1 single nucleotide polymorphisms and fasting insulin but report a nominal association between rs13412852 and BMI (P = 0.042) in a meta-analysis of 8,504 samples from in-house and publicly available studies. Three rare nonsynonymous variants (A353T, R552K, and G582R) were detected in severely insulin-resistant patients. However, these did not cosegregate with disease in affected families, and Lipin1 protein expression and phosphorylation in patients with variants were indistinguishable from those in control subjects.

    Conclusions: Our data do not support a major effect of common LPIN1 variation on metabolic traits and suggest that mutations in LPIN1 are not a common cause of lipodystrophy in humans. The nominal associations with BMI and other metabolic traits in U.K. cohorts require replication in larger cohorts.

    Funded by: Medical Research Council: G0000934, G0000934(68341), G0701446, MC_U106188470, MC_U147574221, MC_U147585824, MC_UP_A620_1014, U.1475.00.002.00001.01 (85824), U.1475.00.004.00002.01(74221); Wellcome Trust: 068545, 068545/Z/02, 077016, 078986, 078986/Z/06/Z, 080952, 080952/Z/06/Z

    Diabetes 2008;57;9;2527-33

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

    Ferreira MA, O'Donovan MC, Meng YA, Jones IR, Ruderfer DM, Jones L, Fan J, Kirov G, Perlis RH, Green EK, Smoller JW, Grozeva D, Stone J, Nikolov I, Chambert K, Hamshere ML, Nimgaonkar VL, Moskvina V, Thase ME, Caesar S, Sachs GS, Franklin J, Gordon-Smith K, Ardlie KG, Gabriel SB, Fraser C, Blumenstiel B, Defelice M, Breen G, Gill M, Morris DW, Elkin A, Muir WJ, McGhee KA, Williamson R, MacIntyre DJ, MacLean AW, St CD, Robinson M, Van Beck M, Pereira AC, Kandaswamy R, McQuillin A, Collier DA, Bass NJ, Young AH, Lawrence J, Ferrier IN, Anjorin A, Farmer A, Curtis D, Scolnick EM, McGuffin P, Daly MJ, Corvin AP, Holmans PA, Blackwood DH, Gurling HM, Owen MJ, Purcell SM, Sklar P, Craddock N and Wellcome Trust Case Control Consortium

    Department of Psychiatry, Massachusetts General Hospital, Boston, Massachusetts 02114, USA.

    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.

    Funded by: Chief Scientist Office; Medical Research Council: G0500791, G0701003, G9309834, G9623693N; NCRR NIH HHS: U54 RR020278; NIMH NIH HHS: MH062137, MH063445, MH067288, MH63420, N01MH80001, R01 MH062137, R01 MH063420, R01 MH063445, R01 MH067288; Wellcome Trust: 076113, 077011, 082371

    Nature genetics 2008;40;9;1056-8

  • The minimum information about a genome sequence (MIGS) specification.

    Field D, Garrity G, Gray T, Morrison N, Selengut J, Sterk P, Tatusova T, Thomson N, Allen MJ, Angiuoli SV, Ashburner M, Axelrod N, Baldauf S, Ballard S, Boore J, Cochrane G, Cole J, Dawyndt P, De Vos P, DePamphilis C, Edwards R, Faruque N, Feldman R, Gilbert J, Gilna P, Glöckner FO, Goldstein P, Guralnick R, Haft D, Hancock D, Hermjakob H, Hertz-Fowler C, Hugenholtz P, Joint I, Kagan L, Kane M, Kennedy J, Kowalchuk G, Kottmann R, Kolker E, Kravitz S, Kyrpides N, Leebens-Mack J, Lewis SE, Li K, Lister AL, Lord P, Maltsev N, Markowitz V, Martiny J, Methe B, Mizrachi I, Moxon R, Nelson K, Parkhill J, Proctor L, White O, Sansone SA, Spiers A, Stevens R, Swift P, Taylor C, Tateno Y, Tett A, Turner S, Ussery D, Vaughan B, Ward N, Whetzel T, San Gil I, Wilson G and Wipat A

    Natural Environmental Research Council Centre for Ecology and Hydrology, Oxford OX1 3SR, UK.

    With the quantity of genomic data increasing at an exponential rate, it is imperative that these data be captured electronically, in a standard format. Standardization activities must proceed within the auspices of open-access and international working bodies. To tackle the issues surrounding the development of better descriptions of genomic investigations, we have formed the Genomic Standards Consortium (GSC). Here, we introduce the minimum information about a genome sequence (MIGS) specification with the intent of promoting participation in its development and discussing the resources that will be required to develop improved mechanisms of metadata capture and exchange. As part of its wider goals, the GSC also supports improving the 'transparency' of the information contained in existing genomic databases.

    Funded by: Biotechnology and Biological Sciences Research Council: BB/E025080/1; Intramural NIH HHS: Z99 LM999999; Medical Research Council: G8225539; NHGRI NIH HHS: U54 HG004028

    Nature biotechnology 2008;26;5;541-7

  • Genetic determinants of ulcerative colitis include the ECM1 locus and five loci implicated in Crohn's disease.

    Fisher SA, Tremelling M, Anderson CA, Gwilliam R, Bumpstead S, Prescott NJ, Nimmo ER, Massey D, Berzuini C, Johnson C, Barrett JC, Cummings FR, Drummond H, Lees CW, Onnie CM, Hanson CE, Blaszczyk K, Inouye M, Ewels P, Ravindrarajah R, Keniry A, Hunt S, Carter M, Watkins N, Ouwehand W, Lewis CM, Cardon L, Wellcome Trust Case Control Consortium, Lobo A, Forbes A, Sanderson J, Jewell DP, Mansfield JC, Deloukas P, Mathew CG, Parkes M and Satsangi J

    Department of Medical and Molecular Genetics, King's College London School of Medicine, 8th Floor Guy's Tower, Guy's Hospital, London SE1 9RT, UK.

    We report results of a nonsynonymous SNP scan for ulcerative colitis and identify a previously unknown susceptibility locus at ECM1. We also show that several risk loci are common to ulcerative colitis and Crohn's disease (IL23R, IL12B, HLA, NKX2-3 and MST1), whereas autophagy genes ATG16L1 and IRGM, along with NOD2 (also known as CARD15), are specific for Crohn's disease. These data provide the first detailed illustration of the genetic relationship between these common inflammatory bowel diseases.

    Funded by: Chief Scientist Office: CZB/4/540; Medical Research Council: G0000934, G0400874, G0600329, G0800383, G0800759, G0802320, MC_QA137934, MC_U105260799; Wellcome Trust: 076113, 077011, 089120, 090532

    Nature genetics 2008;40;6;710-2

  • The Catalogue of Somatic Mutations in Cancer (COSMIC).

    Forbes SA, Bhamra G, Bamford S, Dawson E, Kok C, Clements J, Menzies A, Teague JW, Futreal PA and Stratton MR

    Wellcome Trust Genome Campus, Hinxton, United Kingdom.

    COSMIC is currently the most comprehensive global resource for information on somatic mutations in human cancer, combining curation of the scientific literature with tumor resequencing data from the Cancer Genome Project at the Sanger Institute, U.K. Almost 4800 genes and 250000 tumors have been examined, resulting in over 50000 mutations available for investigation. This information can be accessed in a number of ways, the most convenient being the Web-based system which allows detailed data mining, presenting the results in easily interpretable formats. This unit describes the graphical system in detail, elaborating an example walkthrough and the many ways that the resulting information can be thoroughly investigated by combining data, respecializing the query, or viewing the results in different ways. Alternate protocols overview the available precompiled data files available for download.

    Funded by: Wellcome Trust: 077012

    Current protocols in human genetics 2008;Chapter 10;Unit 10.11

  • SrfB, a member of the Serum Response Factor family of transcription factors, regulates starvation response and early development in Dictyostelium.

    Galardi-Castilla M, Pergolizzi B, Bloomfield G, Skelton J, Ivens A, Kay RR, Bozzaro S and Sastre L

    Instituto de Investigaciones Biomédicas CSIC/UAM. Arturo Duperier, 4. 28029 Madrid, Spain.

    The Serum Response Factor (SRF) is an important regulator of cell proliferation and differentiation. Dictyostelium discoideum srfB gene codes for an SRF homologue and is expressed in vegetative cells and during development under the control of three alternative promoters, which show different cell-type specific patterns of expression. The two more proximal promoters directed gene transcription in prestalk AB, stalk and lower-cup cells. The generation of a strain where the srfB gene has been interrupted (srfB(-)) has shown that this gene is required for regulation of actin-cytoskeleton-related functions, such as cytokinesis and macropinocytosis. The mutant failed to develop well in suspension, but could be rescued by cAMP pulsing, suggesting a defect in cAMP signaling. srfB(-) cells showed impaired chemotaxis to cAMP and defective lateral pseudopodium inhibition. Nevertheless, srfB(-) cells aggregated on agar plates and nitrocellulose filters 2 h earlier than wild type cells, and completed development, showing an increased tendency to form slug structures. Analysis of wild type and srfB(-) strains detected significant differences in the regulation of gene expression upon starvation. Genes coding for lysosomal and ribosomal proteins, developmentally-regulated genes, and some genes coding for proteins involved in cytoskeleton regulation were deregulated during the first stages of development.

    Funded by: Wellcome Trust

    Developmental biology 2008;316;2;260-74

  • ES cell pluripotency and germ-layer formation require the SWI/SNF chromatin remodeling component BAF250a.

    Gao X, Tate P, Hu P, Tjian R, Skarnes WC and Wang Z

    Cardiovascular Research Center, Massachusetts General Hospital, Harvard Medical School, Richard Simches Research Center, 185 Cambridge Street, Boston, MA 02114, USA.

    ATP-dependent chromatin remodeling complexes are a notable group of epigenetic modifiers that use the energy of ATP hydrolysis to change the structure of chromatin, thereby altering its accessibility to nuclear factors. BAF250a (ARID1a) is a unique and defining subunit of the BAF chromatin remodeling complex with the potential to facilitate chromosome alterations critical during development. Our studies show that ablation of BAF250a in early mouse embryos results in developmental arrest (about embryonic day 6.5) and absence of the mesodermal layer, indicating its critical role in early germ-layer formation. Moreover, BAF250a deficiency compromises ES cell pluripotency, severely inhibits self-renewal, and promotes differentiation into primitive endoderm-like cells under normal feeder-free culture conditions. Interestingly, this phenotype can be partially rescued by the presence of embryonic fibroblast cells. DNA microarray, immunostaining, and RNA analyses revealed that BAF250a-mediated chromatin remodeling contributes to the proper expression of numerous genes involved in ES cell self-renewal, including Sox2, Utf1, and Oct4. Furthermore, the pluripotency defects in BAF250a mutant ES cells appear to be cell lineage-specific. For example, embryoid body-based analyses demonstrated that BAF250a-ablated stem cells are defective in differentiating into fully functional mesoderm-derived cardiomyocytes and adipocytes but are capable of differentiating into ectoderm-derived neurons. Our results suggest that BAF250a is a key component of the gene regulatory machinery in ES cells controlling self-renewal, differentiation, and cell lineage decisions.

    Funded by: Howard Hughes Medical Institute

    Proceedings of the National Academy of Sciences of the United States of America 2008;105;18;6656-61

  • Rfam: updates to the RNA families database.

    Gardner PP, Daub J, Tate JG, Nawrocki EP, Kolbe DL, Lindgreen S, Wilkinson AC, Finn RD, Griffiths-Jones S, Eddy SR and Bateman A

    Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, CB10 1SA, UK.

    Rfam is a collection of RNA sequence families, represented by multiple sequence alignments and covariance models (CMs). The primary aim of Rfam is to annotate new members of known RNA families on nucleotide sequences, particularly complete genomes, using sensitive BLAST filters in combination with CMs. A minority of families with a very broad taxonomic range (e.g. tRNA and rRNA) provide the majority of the sequence annotations, whilst the majority of Rfam families (e.g. snoRNAs and miRNAs) have a limited taxonomic range and provide a limited number of annotations. Recent improvements to the website, methodologies and data used by Rfam are discussed. Rfam is freely available on the Web at

    Funded by: Howard Hughes Medical Institute; Wellcome Trust: 077044

    Nucleic acids research 2008;37;Database issue;D136-40

  • Apheresis donors and platelet function: inherent platelet responsiveness influences platelet quality.

    Garner SF, Jones CI, Stephens J, Burns P, Walton J, Bernard A, Angenent W, Ouwehand WH, Goodall AH and BLOODOMICS Consortium

    Department of Haematology, University of Cambridge, Cambridge, UK.

    Background: Process-induced platelet (PLT) activation occurs with all production methods, including apheresis. Recent studies have highlighted the range and consistence of interindividual variation in the PLT response, but little is known about the contribution of a donors' inherent PLT responsiveness to the activation state of the apheresis PLTs or the effect of frequent apheresis on donors' PLTs.

    Study design and methods: The relationship between the donors' PLT response on the apheresis PLTs was studied in 47 individuals selected as having PLTs with inherently low, intermediate, or high responsiveness. Whole-blood flow cytometry was used to measure PLT activation (levels of bound fibrinogen) before donation and in the apheresis PLTs. The effects of regular apheresis on the activation status of donors' PLTs were studied by comparing the in vivo activation status of PLTs from apheresis (n = 349) and whole-blood donors (n = 157), before donation. The effect of apheresis per se on PLT activation was measured in 10 apheresis donors before and after donation.

    Results: The level of PLT activation in the apheresis packs was generally higher than in the donor, and the most activated PLTs were from high-responder donors. There was no significant difference in PLT activation before donation between the apheresis and whole-blood donors (p = 0.697), and there was no consistent evidence of activation in the donors immediately after apheresis.

    Conclusion: The most activated apheresis PLTs were obtained from donors with more responsive PLTs. Regular apheresis, however, does not lead to PLT activation in the donors.

    Transfusion 2008;48;4;673-80

  • The missing link: Bordetella petrii is endowed with both the metabolic versatility of environmental bacteria and virulence traits of pathogenic Bordetellae.

    Gross R, Guzman CA, Sebaihia M, dos Santos VA, Pieper DH, Koebnik R, Lechner M, Bartels D, Buhrmester J, Choudhuri JV, Ebensen T, Gaigalat L, Herrmann S, Khachane AN, Larisch C, Link S, Linke B, Meyer F, Mormann S, Nakunst D, Rückert C, Schneiker-Bekel S, Schulze K, Vorhölter FJ, Yevsa T, Engle JT, Goldman WE, Pühler A, Göbel UB, Goesmann A, Blöcker H, Kaiser O and Martinez-Arias R

    Chair of Microbiology, Biocenter, University of Würzburg, Am Hubland, D-97074 Würzburg, Germany.

    Background: Bordetella petrii is the only environmental species hitherto found among the otherwise host-restricted and pathogenic members of the genus Bordetella. Phylogenetically, it connects the pathogenic Bordetellae and environmental bacteria of the genera Achromobacter and Alcaligenes, which are opportunistic pathogens. B. petrii strains have been isolated from very different environmental niches, including river sediment, polluted soil, marine sponges and a grass root. Recently, clinical isolates associated with bone degenerative disease or cystic fibrosis have also been described.

    Results: In this manuscript we present the results of the analysis of the completely annotated genome sequence of the B. petrii strain DSMZ12804. B. petrii has a mosaic genome of 5,287,950 bp harboring numerous mobile genetic elements, including seven large genomic islands. Four of them are highly related to the clc element of Pseudomonas knackmussii B13, which encodes genes involved in the degradation of aromatics. Though being an environmental isolate, the sequenced B. petrii strain also encodes proteins related to virulence factors of the pathogenic Bordetellae, including the filamentous hemagglutinin, which is a major colonization factor of B. pertussis, and the master virulence regulator BvgAS. However, it lacks all known toxins of the pathogenic Bordetellae.

    Conclusion: The genomic analysis suggests that B. petrii represents an evolutionary link between free-living environmental bacteria and the host-restricted obligate pathogenic Bordetellae. Its remarkable metabolic versatility may enable B. petrii to thrive in very different ecological niches.

    BMC genomics 2008;9;449

  • A novel streptococcal integrative conjugative element involved in iron acquisition.

    Heather Z, Holden MT, Steward KF, Parkhill J, Song L, Challis GL, Robinson C, Davis-Poynter N and Waller AS

    Centre for Preventive Medicine, Animal Health Trust, Lanwades Park, Kentford, Newmarket, Suffolk, UK.

    In this study, we determined the function of a novel non-ribosomal peptide synthetase (NRPS) system carried by a streptococcal integrative conjugative element (ICE), ICESe2. The NRPS shares similarity with the yersiniabactin system found in the high-pathogenicity island of Yersinia sp. and is the first of its kind to be identified in streptococci. We named the NRPS product 'equibactin' and genes of this locus eqbA-N. ICESe2, although absolutely conserved in Streptococcus equi, the causative agent of equine strangles, was absent from all strains of the closely related opportunistic pathogen Streptococcus zooepidemicus. Binding of EqbA, a DtxR-like regulator, to the eqbB promoter was increased in the presence of cations. Deletion of eqbA resulted in a small-colony phenotype. Further deletion of the irp2 homologue eqbE, or the genes eqbH, eqbI and eqbJ encoding a putative ABC transporter, or addition of the iron chelator nitrilotriacetate, reversed this phenotype, implicating iron toxicity. Quantification of (55)Fe accumulation and sensitivity to streptonigrin suggested that equibactin is secreted by S. equi and that the eqbH, eqbI and eqbJ genes are required for its associated iron import. In agreement with a structure-based model of equibactin synthesis, supplementation of chemically defined media with salicylate was required for equibactin production.

    Molecular microbiology 2008;70;5;1274-92

  • Sequence data swell for nematodes.

    Hertz-Fowler C and Pain A

    Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SA, UK.

    With more than 80,000 described species that are extremely diverse in terms of ecology and biology, the Nematoda phylum is one of the most common animal phyla. This month's Genome Watch describes genomes of several nematodes, including that of the human filarial parasite Brugia malayi.

    Nature reviews. Microbiology 2008;6;11;800-1

  • Telomeric expression sites are highly conserved in Trypanosoma brucei.

    Hertz-Fowler C, Figueiredo LM, Quail MA, Becker M, Jackson A, Bason N, Brooks K, Churcher C, Fahkro S, Goodhead I, Heath P, Kartvelishvili M, Mungall K, Harris D, Hauser H, Sanders M, Saunders D, Seeger K, Sharp S, Taylor JE, Walker D, White B, Young R, Cross GA, Rudenko G, Barry JD, Louis EJ and Berriman M

    Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, United Kingdom.

    Subtelomeric regions are often under-represented in genome sequences of eukaryotes. One of the best known examples of the use of telomere proximity for adaptive purposes are the bloodstream expression sites (BESs) of the African trypanosome Trypanosoma brucei. To enhance our understanding of BES structure and function in host adaptation and immune evasion, the BES repertoire from the Lister 427 strain of T. brucei were independently tagged and sequenced. BESs are polymorphic in size and structure but reveal a surprisingly conserved architecture in the context of extensive recombination. Very small BESs do exist and many functioning BESs do not contain the full complement of expression site associated genes (ESAGs). The consequences of duplicated or missing ESAGs, including ESAG9, a newly named ESAG12, and additional variant surface glycoprotein genes (VSGs) were evaluated by functional assays after BESs were tagged with a drug-resistance gene. Phylogenetic analysis of constituent ESAG families suggests that BESs are sequence mosaics and that extensive recombination has shaped the evolution of the BES repertoire. This work opens important perspectives in understanding the molecular mechanisms of antigenic variation, a widely used strategy for immune evasion in pathogens, and telomere biology.

    Funded by: NIAID NIH HHS: R01 AI021729, R01AI021729; Wellcome Trust: 095161

    PloS one 2008;3;10;e3527

  • A Myo6 mutation destroys coordination between the myosin heads, revealing new functions of myosin VI in the stereocilia of mammalian inner ear hair cells.

    Hertzano R, Shalit E, Rzadzinska AK, Dror AA, Song L, Ron U, Tan JT, Shitrit AS, Fuchs H, Hasson T, Ben-Tal N, Sweeney HL, de Angelis MH, Steel KP and Avraham KB

    Department of Human Molecular Genetics and Biochemistry, Tel Aviv University, Tel Aviv, Israel.

    Myosin VI, found in organisms from Caenorhabditis elegans to humans, is essential for auditory and vestibular function in mammals, since genetic mutations lead to hearing impairment and vestibular dysfunction in both humans and mice. Here, we show that a missense mutation in this molecular motor in an ENU-generated mouse model, Tailchaser, disrupts myosin VI function. Structural changes in the Tailchaser hair bundles include mislocalization of the kinocilia and branching of stereocilia. Transfection of GFP-labeled myosin VI into epithelial cells and delivery of endocytic vesicles to the early endosome revealed that the mutant phenotype displays disrupted motor function. The actin-activated ATPase rates measured for the D179Y mutation are decreased, and indicate loss of coordination of the myosin VI heads or 'gating' in the dimer form. Proper coordination is required for walking processively along, or anchoring to, actin filaments, and is apparently destroyed by the proximity of the mutation to the nucleotide-binding pocket. This loss of myosin VI function may not allow myosin VI to transport its cargoes appropriately at the base and within the stereocilia, or to anchor the membrane of stereocilia to actin filaments via its cargos, both of which lead to structural changes in the stereocilia of myosin VI-impaired hair cells, and ultimately leading to deafness.

    Funded by: Medical Research Council: G0300212, MC_QA137918; NEI NIH HHS: R01 EY012695, R01-EY12695; NIDCD NIH HHS: R01-DC0099100; Wellcome Trust

    PLoS genetics 2008;4;10;e1000207

  • The genome sequence of the fish pathogen Aliivibrio salmonicida strain LFI1238 shows extensive evidence of gene decay.

    Hjerde E, Lorentzen MS, Holden MT, Seeger K, Paulsen S, Bason N, Churcher C, Harris D, Norbertczak H, Quail MA, Sanders S, Thurston S, Parkhill J, Willassen NP and Thomson NR

    Department of Molecular Biotechnology, Institute of Medical Biology, Faculty of Medicine, University of Tromsø, N-9037 Tromsø, Norway.

    Background: The fish pathogen Aliivibrio salmonicida is the causative agent of cold-water vibriosis in marine aquaculture. The Gram-negative bacterium causes tissue degradation, hemolysis and sepsis in vivo.

    Results: In total, 4 286 protein coding sequences were identified, and the 4.6 Mb genome of A. salmonicida has a six partite architecture with two chromosomes and four plasmids. Sequence analysis revealed a highly fragmented genome structure caused by the insertion of an extensive number of insertion sequence (IS) elements. The IS elements can be related to important evolutionary events such as gene acquisition, gene loss and chromosomal rearrangements. New A. salmonicida functional capabilities that may have been aquired through horizontal DNA transfer include genes involved in iron-acquisition, and protein secretion and play potential roles in pathogenicity. On the other hand, the degeneration of 370 genes and consequent loss of specific functions suggest that A. salmonicida has a reduced metabolic and physiological capacity in comparison to related Vibrionaceae species.

    Conclusion: Most prominent is the loss of several genes involved in the utilisation of the polysaccharide chitin. In particular, the disruption of three extracellular chitinases responsible for enzymatic breakdown of chitin makes A. salmonicida unable to grow on the polymer form of chitin. These, and other losses could restrict the variety of carrier organisms A. salmonicida can attach to, and associate with. Gene acquisition and gene loss may be related to the emergence of A. salmonicida as a fish pathogen.

    BMC genomics 2008;9;616

  • The genome of Burkholderia cenocepacia J2315, an epidemic pathogen of cystic fibrosis patients.

    Holden MT, Seth-Smith HM, Crossman LC, Sebaihia M, Bentley SD, Cerdeño-Tárraga AM, Thomson NR, Bason N, Quail MA, Sharp S, Cherevach I, Churcher C, Goodhead I, Hauser H, Holroyd N, Mungall K, Scott P, Walker D, White B, Rose H, Iversen P, Mil-Homens D, Rocha EP, Fialho AM, Baldwin A, Dowson C, Barrell BG, Govan JR, Vandamme P, Hart CA, Mahenthiralingam E and Parkhill J

    The Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Cambridge CB10 1SA, United Kingdom.

    Bacterial infections of the lungs of cystic fibrosis (CF) patients cause major complications in the treatment of this common genetic disease. Burkholderia cenocepacia infection is particularly problematic since this organism has high levels of antibiotic resistance, making it difficult to eradicate; the resulting chronic infections are associated with severe declines in lung function and increased mortality rates. B. cenocepacia strain J2315 was isolated from a CF patient and is a member of the epidemic ET12 lineage that originated in Canada or the United Kingdom and spread to Europe. The 8.06-Mb genome of this highly transmissible pathogen comprises three circular chromosomes and a plasmid and encodes a broad array of functions typical of this metabolically versatile genus, as well as numerous virulence and drug resistance functions. Although B. cenocepacia strains can be isolated from soil and can be pathogenic to both plants and man, J2315 is representative of a lineage of B. cenocepacia rarely isolated from the environment and which spreads between CF patients. Comparative analysis revealed that ca. 21% of the genome is unique in comparison to other strains of B. cenocepacia, highlighting the genomic plasticity of this species. Pseudogenes in virulence determinants suggest that the pathogenic response of J2315 may have been recently selected to promote persistence in the CF lung. The J2315 genome contains evidence that its unique and highly adapted genetic content has played a significant role in its success as an epidemic CF pathogen.

    Funded by: Wellcome Trust

    Journal of bacteriology 2008;191;1;261-77

  • High-throughput sequencing provides insights into genome variation and evolution in Salmonella Typhi.

    Holt KE, Parkhill J, Mazzoni CJ, Roumagnac P, Weill FX, Goodhead I, Rance R, Baker S, Maskell DJ, Wain J, Dolecek C, Achtman M and Dougan G

    The Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge CB10 1SA, UK.

    Isolates of Salmonella enterica serovar Typhi (Typhi), a human-restricted bacterial pathogen that causes typhoid, show limited genetic variation. We generated whole-genome sequences for 19 Typhi isolates using 454 (Roche) and Solexa (Illumina) technologies. Isolates, including the previously sequenced CT18 and Ty2 isolates, were selected to represent major nodes in the phylogenetic tree. Comparative analysis showed little evidence of purifying selection, antigenic variation or recombination between isolates. Rather, evolution in the Typhi population seems to be characterized by ongoing loss of gene function, consistent with a small effective population size. The lack of evidence for antigenic variation driven by immune selection is in contrast to strong adaptive selection for mutations conferring antibiotic resistance in Typhi. The observed patterns of genetic isolation and drift are consistent with the proposed key role of asymptomatic carriers of Typhi as the main reservoir of this pathogen, highlighting the need for identification and treatment of carriers.

    Funded by: Wellcome Trust: 067321

    Nature genetics 2008;40;8;987-93

  • Variation analysis and gene annotation of eight MHC haplotypes: the MHC Haplotype Project.

    Horton R, Gibson R, Coggill P, Miretti M, Allcock RJ, Almeida J, Forbes S, Gilbert JG, Halls K, Harrow JL, Hart E, Howe K, Jackson DK, Palmer S, Roberts AN, Sims S, Stewart CA, Traherne JA, Trevanion S, Wilming L, Rogers J, de Jong PJ, Elliott JF, Sawcer S, Todd JA, Trowsdale J and Beck S

    Wellcome Trust Sanger Institute, Genome Campus, Hinxton, Cambridge, CB10 1SA, UK.

    The human major histocompatibility complex (MHC) is contained within about 4 Mb on the short arm of chromosome 6 and is recognised as the most variable region in the human genome. The primary aim of the MHC Haplotype Project was to provide a comprehensively annotated reference sequence of a single, human leukocyte antigen-homozygous MHC haplotype and to use it as a basis against which variations could be assessed from seven other similarly homozygous cell lines, representative of the most common MHC haplotypes in the European population. Comparison of the haplotype sequences, including four haplotypes not previously analysed, resulted in the identification of >44,000 variations, both substitutions and indels (insertions and deletions), which have been submitted to the dbSNP database. The gene annotation uncovered haplotype-specific differences and confirmed the presence of more than 300 loci, including over 160 protein-coding genes. Combined analysis of the variation and annotation datasets revealed 122 gene loci with coding substitutions of which 97 were non-synonymous. The haplotype (A3-B7-DR15; PGF cell line) designated as the new MHC reference sequence, has been incorporated into the human genome assembly (NCBI35 and subsequent builds), and constitutes the largest single-haplotype sequence of the human genome to date. The extensive variation and annotation data derived from the analysis of seven further haplotypes have been made publicly available and provide a framework and resource for future association studies of all MHC-associated diseases and transplant medicine.

    Funded by: NHGRI NIH HHS: U54 HG004555; Wellcome Trust: 048880, 062023, 077198

    Immunogenetics 2008;60;1;1-18

  • Array painting reveals a high frequency of balanced translocations in breast cancer cell lines that break in cancer-relevant genes.

    Howarth KD, Blood KA, Ng BL, Beavis JC, Chua Y, Cooke SL, Raby S, Ichimura K, Collins VP, Carter NP and Edwards PA

    Department of Pathology, Hutchison-MRC Research Centre, University of Cambridge, Cambridge, UK.

    Chromosome translocations in the common epithelial cancers are abundant, yet little is known about them. They have been thought to be almost all unbalanced and therefore dismissed as mostly mediating tumour suppressor loss. We present a comprehensive analysis by array painting of the chromosome translocations of breast cancer cell lines HCC1806, HCC1187 and ZR-75-30. In array painting, chromosomes are isolated by flow cytometry, amplified and hybridized to DNA microarrays. A total of 200 breakpoints were identified and all were mapped to 1 Mb resolution on bacterial artificial chromosome (BAC) arrays, then 40 selected breakpoints, including all balanced breakpoints, were further mapped on tiling-path BAC arrays or to around 2 kb resolution using oligonucleotide arrays. Many more of the translocations were balanced at 1 Mb resolution than expected, either reciprocal (eight in total) or balanced for at least one participating chromosome (19 paired breakpoints). Second, many of the breakpoints were at genes that are plausible targets of oncogenic translocation, including balanced breaks at CTCF, EP300/p300 and FOXP4. Two gene fusions were demonstrated, TAX1BP1-AHCY and RIF1-PKD1L1. Our results support the idea that chromosome rearrangements may play an important role in common epithelial cancers such as breast cancer.

    Funded by: Cancer Research UK: A4392; Wellcome Trust: 077008

    Oncogene 2008;27;23;3345-59

  • Newly identified genetic risk variants for celiac disease related to the immune response.

    Hunt KA, Zhernakova A, Turner G, Heap GA, Franke L, Bruinenberg M, Romanos J, Dinesen LC, Ryan AW, Panesar D, Gwilliam R, Takeuchi F, McLaren WM, Holmes GK, Howdle PD, Walters JR, Sanders DS, Playford RJ, Trynka G, Mulder CJ, Mearin ML, Verbeek WH, Trimble V, Stevens FM, O'Morain C, Kennedy NP, Kelleher D, Pennington DJ, Strachan DP, McArdle WL, Mein CA, Wapenaar MC, Deloukas P, McGinnis R, McManus R, Wijmenga C and van Heel DA

    Institute of Cell and Molecular Science, Barts and The London School of Medicine and Dentistry, 4 Newark Street, London E1 2AT, UK.

    Our genome-wide association study of celiac disease previously identified risk variants in the IL2-IL21 region. To identify additional risk variants, we genotyped 1,020 of the most strongly associated non-HLA markers in an additional 1,643 cases and 3,406 controls. Through joint analysis including the genome-wide association study data (767 cases, 1,422 controls), we identified seven previously unknown risk regions (P < 5 x 10(-7)). Six regions harbor genes controlling immune responses, including CCR3, IL12A, IL18RAP, RGS1, SH2B3 (nsSNP rs3184504) and TAGAP. Whole-blood IL18RAP mRNA expression correlated with IL18RAP genotype. Type 1 diabetes and celiac disease share HLA-DQ, IL2-IL21, CCR3 and SH2B3 risk regions. Thus, this extensive genome-wide association follow-up study has identified additional celiac disease risk variants in relevant biological pathways.

    Funded by: Medical Research Council: G0000934; Wellcome Trust: 068094, 068545/Z/02, 084743, GR068094MA

    Nature genetics 2008;40;4;395-402

  • The functional impact of structural variation in humans.

    Hurles ME, Dermitzakis ET and Tyler-Smith C

    The Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge CB10 1SA, UK.

    Structural variation includes many different types of chromosomal rearrangement and encompasses millions of bases in every human genome. Over the past 3 years, the extent and complexity of structural variation has become better appreciated. Diverse approaches have been adopted to explore the functional impact of this class of variation. As disparate indications of the important biological consequences of genome dynamism are accumulating rapidly, we review the evidence that structural variation has an appreciable impact on cellular phenotypes, disease and human evolution.

    Funded by: Wellcome Trust: 077009, 077014, 077046

    Trends in genetics : TIG 2008;24;5;238-45

  • Complete genome sequence and comparative genome analysis of enteropathogenic Escherichia coli O127:H6 strain E2348/69.

    Iguchi A, Thomson NR, Ogura Y, Saunders D, Ooka T, Henderson IR, Harris D, Asadulghani M, Kurokawa K, Dean P, Kenny B, Quail MA, Thurston S, Dougan G, Hayashi T, Parkhill J and Frankel G

    Division of Bioenvironmental Science, Frontier Science Research Center, University of Miyazaki, Miyazaki, Japan.

    Enteropathogenic Escherichia coli (EPEC) was the first pathovar of E. coli to be implicated in human disease; however, no EPEC strain has been fully sequenced until now. Strain E2348/69 (serotype O127:H6 belonging to E. coli phylogroup B2) has been used worldwide as a prototype strain to study EPEC biology, genetics, and virulence. Studies of E2348/69 led to the discovery of the locus of enterocyte effacement-encoded type III secretion system (T3SS) and its cognate effectors, which play a vital role in attaching and effacing lesion formation on gut epithelial cells. In this study, we determined the complete genomic sequence of E2348/69 and performed genomic comparisons with other important E. coli strains. We identified 424 E2348/69-specific genes, most of which are carried on mobile genetic elements, and a number of genetic traits specifically conserved in phylogroup B2 strains irrespective of their pathotypes, including the absence of the ETT2-related T3SS, which is present in E. coli strains belonging to all other phylogroups. The genome analysis revealed the entire gene repertoire related to E2348/69 virulence. Interestingly, E2348/69 contains only 21 intact T3SS effector genes, all of which are carried on prophages and integrative elements, compared to over 50 effector genes in enterohemorrhagic E. coli O157. As E2348/69 is the most-studied pathogenic E. coli strain, this study provides a genomic context for the vast amount of existing experimental data. The unexpected simplicity of the E2348/69 T3SS provides the first opportunity to fully dissect the entire virulence strategy of attaching and effacing pathogens in the genomic context.

    Funded by: Medical Research Council: G0700151; Wellcome Trust

    Journal of bacteriology 2008;191;1;347-54

  • A novel CpG island set identifies tissue-specific methylation at developmental gene loci.

    Illingworth R, Kerr A, Desousa D, Jørgensen H, Ellis P, Stalker J, Jackson D, Clee C, Plumb R, Rogers J, Humphray S, Cox T, Langford C and Bird A

    Wellcome Trust Centre for Cell Biology, University of Edinburgh, Edinburgh, United Kingdom.

    CpG islands (CGIs) are dense clusters of CpG sequences that punctuate the CpG-deficient human genome and associate with many gene promoters. As CGIs also differ from bulk chromosomal DNA by their frequent lack of cytosine methylation, we devised a CGI enrichment method based on nonmethylated CpG affinity chromatography. The resulting library was sequenced to define a novel human blood CGI set that includes many that are not detected by current algorithms. Approximately half of CGIs were associated with annotated gene transcription start sites, the remainder being intra- or intergenic. Using an array representing over 17,000 CGIs, we established that 6%-8% of CGIs are methylated in genomic DNA of human blood, brain, muscle, and spleen. Inter- and intragenic CGIs are preferentially susceptible to methylation. CGIs showing tissue-specific methylation were overrepresented at numerous genetic loci that are essential for development, including HOX and PAX family members. The findings enable a comprehensive analysis of the roles played by CGI methylation in normal and diseased human tissues.

    Funded by: Wellcome Trust

    PLoS biology 2008;6;1;e22

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

    Jakkula E, Rehnström K, Varilo T, Pietiläinen OP, Paunio T, Pedersen NL, deFaire U, Järvelin MR, Saharinen J, Freimer N, Ripatti S, Purcell S, Collins A, Daly MJ, Palotie A and Peltonen L

    Department of Molecular Medicine, National Public Health Institute and Institute for Molecular Medicine Finland (FIMM), Helsinki, Finland.

    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.

    Funded by: NCRR NIH HHS: U54 RR020278, U54RR020278; NHLBI NIH HHS: 1R01HL087679-01, R01 HL087679

    American journal of human genetics 2008;83;6;787-94

  • A systematic library for comprehensive overexpression screens in Saccharomyces cerevisiae.

    Jones GM, Stalker J, Humphray S, West A, Cox T, Rogers J, Dunham I and Prelich G

    Department of Molecular Genetics, Albert Einstein College of Medicine, 1300 Morris Park Avenue, Bronx, New York 10461, USA.

    Modern genetic analysis requires the development of new resources to systematically explore gene function in vivo. Overexpression screens are a powerful method to investigate genetic pathways, but the goal of routine and comprehensive overexpression screens has been hampered by the lack of systematic libraries. Here we describe the construction of a systematic collection of the Saccharomyces cerevisiae genome in a high-copy vector and its validation in two overexpression screens.

    Funded by: NIGMS NIH HHS: GM52486; Wellcome Trust

    Nature methods 2008;5;3;239-41

  • Six new loci associated with blood low-density lipoprotein cholesterol, high-density lipoprotein cholesterol or triglycerides in humans.

    Kathiresan S, Melander O, Guiducci C, Surti A, Burtt NP, Rieder MJ, Cooper GM, Roos C, Voight BF, Havulinna AS, Wahlstrand B, Hedner T, Corella D, Tai ES, Ordovas JM, Berglund G, Vartiainen E, Jousilahti P, Hedblad B, Taskinen MR, Newton-Cheh C, Salomaa V, Peltonen L, Groop L, Altshuler DM and Orho-Melander M

    Cardiology Division, Massachusetts General Hospital, Boston, Massachusetts 02114, USA.

    Blood concentrations of lipoproteins and lipids are heritable risk factors for cardiovascular disease. Using genome-wide association data from three studies (n = 8,816 that included 2,758 individuals from the Diabetes Genetics Initiative specific to the current paper as well as 1,874 individuals from the FUSION study of type 2 diabetes and 4,184 individuals from the SardiNIA study of aging-associated variables reported in a companion paper in this issue) and targeted replication association analyses in up to 18,554 independent participants, we show that common SNPs at 18 loci are reproducibly associated with concentrations of low-density lipoprotein (LDL) cholesterol, high-density lipoprotein (HDL) cholesterol, and/or triglycerides. Six of these loci are new (P < 5 x 10(-8) for each new locus). Of the six newly identified chromosomal regions, two were associated with LDL cholesterol (1p13 near CELSR2, PSRC1 and SORT1 and 19p13 near CILP2 and PBX4), one with HDL cholesterol (1q42 in GALNT2) and five with triglycerides (7q11 near TBL2 and MLXIPL, 8q24 near TRIB1, 1q42 in GALNT2, 19p13 near CILP2 and PBX4 and 1p31 near ANGPTL3). At 1p13, the LDL-associated SNP was also strongly correlated with CELSR2, PSRC1, and SORT1 transcript levels in human liver, and a proxy for this SNP was recently shown to affect risk for coronary artery disease. Understanding the molecular, cellular and clinical consequences of the newly identified loci may inform therapy and clinical care.

    Funded by: Wellcome Trust: 089061

    Nature genetics 2008;40;2;189-97

  • iMapper: a web application for the automated analysis and mapping of insertional mutagenesis sequence data against Ensembl genomes.

    Kong J, Zhu F, Stalker J and Adams DJ

    Experimental Cancer Genetics, The Wellcome Trust Sanger Institute, Hinxton, Cambs CB101HH, UK.

    Summary: Insertional mutagenesis is a powerful method for gene discovery. To identify the location of insertion sites in the genome linker based polymerase chain reaction (PCR) methods (such as splinkerette-PCR) may be employed. We have developed a web application called iMapper (Insertional Mutagenesis Mapping and Analysis Tool) for the efficient analysis of insertion site sequence reads against vertebrate and invertebrate Ensembl genomes. Taking linker based sequences as input, iMapper scans and trims the sequence to remove the linker and sequences derived from the insertional mutagen. The software then identifies and removes contaminating sequences derived from chimeric genomic fragments, vector or the transposon concatamer and then presents the clipped sequence reads to a sequence mapping server which aligns them to an Ensembl genome. Insertion sites can then be navigated in Ensembl in the context of genomic features such as gene structures. iMapper also generates test-based format for nucleic acid or protein sequences (FASTA) and generic file format (GFF) files of the clipped sequence reads and provides a graphical overview of the mapped insertion sites against a karyotype. iMapper is designed for high-throughput applications and can efficiently process thousands of DNA sequence reads.

    Availability: iMapper is web based and can be accessed at

    Funded by: Cancer Research UK: C20510/A6997; Wellcome Trust: 76943

    Bioinformatics (Oxford, England) 2008;24;24;2923-5

  • Normal germ line establishment in mice carrying a deletion of the Ifitm/Fragilis gene family cluster.

    Lange UC, Adams DJ, Lee C, Barton S, Schneider R, Bradley A and Surani MA

    Wellcome Trust/Cancer Research UK Gurdon Institute of Cancer & Developmental Biology, University of Cambridge, Tennis Court Road, Cambridge CB2 1QN, United Kingdom.

    The family of interferon-inducible transmembrane proteins (Ifitm) consists of five highly sequence-related cell surface proteins, which are implicated in diverse cellular processes. Ifitm genes are conserved, widely expressed, and characteristically found in genomic clusters, such as the 67-kb Ifitm family locus on mouse chromosome 7. Recently, Ifitm1 and Ifitm3 have been suggested to mediate migration of early primordial germ cells (PGCs), a process that is little understood. To investigate Ifitm function during germ cell development, we used targeted chromosome engineering to generate mutants which either lack the entire Ifitm locus or carry a disrupted Ifitm3 gene only. Here we show that the mutations have no detectable effects on development of the germ line or on the generation of live young. Hence, contrary to previous reports, Ifitm genes are not essential for PGC migration. The Ifitm family is a striking example of a conserved gene cluster which appears to be functionally redundant during development.

    Funded by: Biotechnology and Biological Sciences Research Council; Cancer Research UK; Wellcome Trust: 065601

    Molecular and cellular biology 2008;28;15;4688-96

  • A single gene network accurately predicts phenotypic effects of gene perturbation in Caenorhabditis elegans.

    Lee I, Lehner B, Crombie C, Wong W, Fraser AG and Marcotte EM

    Center for Systems and Synthetic Biology, Department of Chemistry and Biochemistry, Institute for Cellular and Molecular Biology, University of Texas, 2500 Speedway, MBB 3.210, Austin, Texas 78712, USA.

    The fundamental aim of genetics is to understand how an organism's phenotype is determined by its genotype, and implicit in this is predicting how changes in DNA sequence alter phenotypes. A single network covering all the genes of an organism might guide such predictions down to the level of individual cells and tissues. To validate this approach, we computationally generated a network covering most C. elegans genes and tested its predictive capacity. Connectivity within this network predicts essentiality, identifying this relationship as an evolutionarily conserved biological principle. Critically, the network makes tissue-specific predictions-we accurately identify genes for most systematically assayed loss-of-function phenotypes, which span diverse cellular and developmental processes. Using the network, we identify 16 genes whose inactivation suppresses defects in the retinoblastoma tumor suppressor pathway, and we successfully predict that the dystrophin complex modulates EGF signaling. We conclude that an analogous network for human genes might be similarly predictive and thus facilitate identification of disease genes and rational therapeutic targets.

    Funded by: NIGMS NIH HHS: GM06779-01; Wellcome Trust

    Nature genetics 2008;40;2;181-8

  • Evolution of the Rhodococcus equi vap pathogenicity island seen through comparison of host-associated vapA and vapB virulence plasmids.

    Letek M, Ocampo-Sosa AA, Sanders M, Fogarty U, Buckley T, Leadon DP, González P, Scortti M, Meijer WG, Parkhill J, Bentley S and Vázquez-Boland JA

    Division of Microbial Pathogenesis, Centre for Infectious Diseases, Ashworth Laboratories, King's Buildings, University of Edinburgh, Edinburgh EH9 3JT, United Kingdom.

    The pathogenic actinomycete Rhodococcus equi harbors different types of virulence plasmids associated with specific nonhuman hosts. We determined the complete DNA sequence of a vapB(+) plasmid, typically associated with pig isolates, and compared it with that of the horse-specific vapA(+) plasmid type. pVAPB1593, a circular 79,251-bp element, had the same housekeeping backbone as the vapA(+) plasmid but differed over an approximately 22-kb region. This variable region encompassed the vap pathogenicity island (PAI), was clearly subject to selective pressures different from those affecting the backbone, and showed major genetic rearrangements involving the vap genes. The pVAPB1593 PAI harbored five different vap genes (vapB and vapJ to -M, with vapK present in two copies), which encoded products differing by 24 to 84% in amino acid sequence from the six full-length vapA(+) plasmid-encoded Vap proteins, consistent with a role for the specific vap gene complement in R. equi host tropism. Sequence analyses, including interpolated variable-order motifs for detection of alien DNA and reconstruction of Vap family phylogenetic relationships, suggested that the vap PAI was acquired by an ancestor plasmid via lateral gene transfer, subsequently evolving by vap gene duplication and sequence diversification to give different (host-adapted) plasmids. The R. equi virulence plasmids belong to a new family of actinobacterial circular replicons characterized by an ancient conjugative backbone and a horizontally acquired niche-adaptive plasticity region.

    Journal of bacteriology 2008;190;17;5797-805

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

    Lettre G, Jackson AU, Gieger C, Schumacher FR, Berndt SI, Sanna S, Eyheramendy S, Voight BF, Butler JL, Guiducci C, Illig T, Hackett R, Heid IM, Jacobs KB, Lyssenko V, Uda M, Diabetes Genetics Initiative, FUSION, KORA, Prostate, Lung Colorectal and Ovarian Cancer Screening Trial, Nurses' Health Study, SardiNIA, Boehnke M, Chanock SJ, Groop LC, Hu FB, Isomaa B, Kraft P, Peltonen L, Salomaa V, Schlessinger D, Hunter DJ, Hayes RB, Abecasis GR, Wichmann HE, Mohlke KL and Hirschhorn JN

    Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, Cambridge, Massachusetts 02142, USA.

    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.

    Funded by: Intramural NIH HHS; NCI NIH HHS: 5P01CA087969, 5U01CA098233, CA49449; NHGRI NIH HHS: HG02651; NHLBI NIH HHS: HL084729; NIDDK NIH HHS: 5 R01 DK 075787, DK62370, DK72193, R01 DK072193; Wellcome Trust: 089061

    Nature genetics 2008;40;5;584-91

  • Mapping short DNA sequencing reads and calling variants using mapping quality scores.

    Li H, Ruan J and Durbin R

    The Wellcome Trust Sanger Institute, Hinxton CB10 1SA, United Kingdom.

    New sequencing technologies promise a new era in the use of DNA sequence. However, some of these technologies produce very short reads, typically of a few tens of base pairs, and to use these reads effectively requires new algorithms and software. In particular, there is a major issue in efficiently aligning short reads to a reference genome and handling ambiguity or lack of accuracy in this alignment. Here we introduce the concept of mapping quality, a measure of the confidence that a read actually comes from the position it is aligned to by the mapping algorithm. We describe the software MAQ that can build assemblies by mapping shotgun short reads to a reference genome, using quality scores to derive genotype calls of the consensus sequence of a diploid genome, e.g., from a human sample. MAQ makes full use of mate-pair information and estimates the error probability of each read alignment. Error probabilities are also derived for the final genotype calls, using a Bayesian statistical model that incorporates the mapping qualities, error probabilities from the raw sequence quality scores, sampling of the two haplotypes, and an empirical model for correlated errors at a site. Both read mapping and genotype calling are evaluated on simulated data and real data. MAQ is accurate, efficient, versatile, and user-friendly. It is freely available at

    Funded by: Wellcome Trust

    Genome research 2008;18;11;1851-8

  • Lifelong reduction of LDL-cholesterol related to a common variant in the LDL-receptor gene decreases the risk of coronary artery disease--a Mendelian Randomisation study.

    Linsel-Nitschke P, Götz A, Erdmann J, Braenne I, Braund P, Hengstenberg C, Stark K, Fischer M, Schreiber S, El Mokhtari NE, Schaefer A, Schrezenmeir J, Schrezenmeier J, Rubin D, Hinney A, Reinehr T, Roth C, Ortlepp J, Hanrath P, Hall AS, Mangino M, Lieb W, Lamina C, Heid IM, Doering A, Gieger C, Peters A, Meitinger T, Wichmann HE, König IR, Ziegler A, Kronenberg F, Samani NJ, Schunkert H, Wellcome Trust Case Control Consortium (WTCCC) and Cardiogenics Consortium

    Medizinische Klinik II, Universität zu Lübeck, Lübeck, Germany.

    Background: Rare mutations of the low-density lipoprotein receptor gene (LDLR) cause familial hypercholesterolemia, which increases the risk for coronary artery disease (CAD). Less is known about the implications of common genetic variation in the LDLR gene regarding the variability of cholesterol levels and risk of CAD.

    Methods: Imputed genotype data at the LDLR locus on 1 644 individuals of a population-based sample were explored for association with LDL-C level. Replication of association with LDL-C level was sought for the most significant single nucleotide polymorphism (SNP) within the LDLR gene in three European samples comprising 6 642 adults and 533 children. Association of this SNP with CAD was examined in six case-control studies involving more than 15 000 individuals.

    Findings: Each copy of the minor T allele of SNP rs2228671 within LDLR (frequency 11%) was related to a decrease of LDL-C levels by 0.19 mmol/L (95% confidence interval (CI) [0.13-0.24] mmol/L, p = 1.5x10(-10)). This association with LDL-C was uniformly found in children, men, and women of all samples studied. In parallel, the T allele of rs2228671 was associated with a significantly lower risk of CAD (Odds Ratio per copy of the T allele: 0.82, 95% CI [0.76-0.89], p = 2.1x10(-7)). Adjustment for LDL-C levels by logistic regression or Mendelian Randomisation models abolished the significant association between rs2228671 with CAD completely, indicating a functional link between the genetic variant at the LDLR gene locus, change in LDL-C and risk of CAD.

    Conclusion: A common variant at the LDLR gene locus affects LDL-C levels and, thereby, the risk for CAD.

    Funded by: British Heart Foundation; Medical Research Council

    PloS one 2008;3;8;e2986

  • The conserved plant sterility gene HAP2 functions after attachment of fusogenic membranes in Chlamydomonas and Plasmodium gametes.

    Liu Y, Tewari R, Ning J, Blagborough AM, Garbom S, Pei J, Grishin NV, Steele RE, Sinden RE, Snell WJ and Billker O

    Department of Cell Biology, University of Texas Southwestern Medical Center, Dallas, Texas 75390, USA.

    The cellular and molecular mechanisms that underlie species-specific membrane fusion between male and female gametes remain largely unknown. Here, by use of gene discovery methods in the green alga Chlamydomonas, gene disruption in the rodent malaria parasite Plasmodium berghei, and distinctive features of fertilization in both organisms, we report discovery of a mechanism that accounts for a conserved protein required for gamete fusion. A screen for fusion mutants in Chlamydomonas identified a homolog of HAP2, an Arabidopsis sterility gene. Moreover, HAP2 disruption in Plasmodium blocked fertilization and thereby mosquito transmission of malaria. HAP2 localizes at the fusion site of Chlamydomonas minus gametes, yet Chlamydomonas minus and Plasmodium hap2 male gametes retain the ability, using other, species-limited proteins, to form tight prefusion membrane attachments with their respective gamete partners. Membrane dye experiments show that HAP2 is essential for membrane merger. Thus, in two distantly related eukaryotes, species-limited proteins govern access to a conserved protein essential for membrane fusion.

    Funded by: Biotechnology and Biological Sciences Research Council; Medical Research Council: G0501670; NIGMS NIH HHS: GM056778, R01 GM056778; Wellcome Trust

    Genes & development 2008;22;8;1051-68

  • Common variants near MC4R are associated with fat mass, weight and risk of obesity.

    Loos RJ, Lindgren CM, Li S, Wheeler E, Zhao JH, Prokopenko I, Inouye M, Freathy RM, Attwood AP, Beckmann JS, Berndt SI, Prostate, Lung, Colorectal, and Ovarian (PLCO) Cancer Screening Trial, Jacobs KB, Chanock SJ, Hayes RB, Bergmann S, Bennett AJ, Bingham SA, Bochud M, Brown M, Cauchi S, Connell JM, Cooper C, Smith GD, Day I, Dina C, De S, Dermitzakis ET, Doney AS, Elliott KS, Elliott P, Evans DM, Sadaf Farooqi I, Froguel P, Ghori J, Groves CJ, Gwilliam R, Hadley D, Hall AS, Hattersley AT, Hebebrand J, Heid IM, KORA, Lamina C, Gieger C, Illig T, Meitinger T, Wichmann HE, Herrera B, Hinney A, Hunt SE, Jarvelin MR, Johnson T, Jolley JD, Karpe F, Keniry A, Khaw KT, Luben RN, Mangino M, Marchini J, McArdle WL, McGinnis R, Meyre D, Munroe PB, Morris AD, Ness AR, Neville MJ, Nica AC, Ong KK, O'Rahilly S, Owen KR, Palmer CN, Papadakis K, Potter S, Pouta A, Qi L, Nurses' Health Study, Randall JC, Rayner NW, Ring SM, Sandhu MS, Scherag A, Sims MA, Song K, Soranzo N, Speliotes EK, Diabetes Genetics Initiative, Syddall HE, Teichmann SA, Timpson NJ, Tobias JH, Uda M, SardiNIA Study, Vogel CI, Wallace C, Waterworth DM, Weedon MN, Wellcome Trust Case Control Consortium, Willer CJ, FUSION, Wraight, Yuan X, Zeggini E, Hirschhorn JN, Strachan DP, Ouwehand WH, Caulfield MJ, Samani NJ, Frayling TM, Vollenweider P, Waeber G, Mooser V, Deloukas P, McCarthy MI, Wareham NJ, Barroso I, Jacobs KB, Chanock SJ, Hayes RB, Lamina C, Gieger C, Illig T, Meitinger T, Wichmann HE, Kraft P, Hankinson SE, Hunter DJ, Hu FB, Lyon HN, Voight BF, Ridderstrale M, Groop L, Scheet P, Sanna S, Abecasis GR, Albai G, Nagaraja R, Schlessinger D, Jackson AU, Tuomilehto J, Collins FS, Boehnke M and Mohlke KL

    MRC Epidemiology Unit, Addenbrooke's Hospital, Cambridge CB2 0QQ, UK.

    To identify common variants influencing body mass index (BMI), we analyzed genome-wide association data from 16,876 individuals of European descent. After previously reported variants in FTO, the strongest association signal (rs17782313, P = 2.9 x 10(-6)) mapped 188 kb downstream of MC4R (melanocortin-4 receptor), mutations of which are the leading cause of monogenic severe childhood-onset obesity. We confirmed the BMI association in 60,352 adults (per-allele effect = 0.05 Z-score units; P = 2.8 x 10(-15)) and 5,988 children aged 7-11 (0.13 Z-score units; P = 1.5 x 10(-8)). In case-control analyses (n = 10,583), the odds for severe childhood obesity reached 1.30 (P = 8.0 x 10(-11)). Furthermore, we observed overtransmission of the risk allele to obese offspring in 660 families (P (pedigree disequilibrium test average; PDT-avg) = 2.4 x 10(-4)). The SNP location and patterns of phenotypic associations are consistent with effects mediated through altered MC4R function. Our findings establish that common variants near MC4R influence fat mass, weight and obesity risk at the population level and reinforce the need for large-scale data integration to identify variants influencing continuous biomedical traits.

    Funded by: British Heart Foundation; Cancer Research UK; Department of Health: DHCS/07/07/008; Medical Research Council: G0000934, G0000934(68341), G0400874, G0401527, G0600331, G0600705, G0601261, G0701863, G9521010, G9521010(63660), G9824984, G9828345, MC_QA137934, MC_U105161047, MC_U105630924, MC_U106179472, MC_U106188470, MC_U147585824, MC_UP_A620_1014, U.1475.00.002.00001.01 (85824); NIDDK NIH HHS: F32 DK079466, F32 DK079466-01, K23 DK080145, K23 DK080145-01, P30 DK040561, P30 DK040561-13, R01 DK072193; Wellcome Trust: 068545, 076113, 077016, 079557, 084713, 090532

    Nature genetics 2008;40;6;768-75

  • Sites of strong Rec12/Spo11 binding in the fission yeast genome are associated with meiotic recombination and with centromeres.

    Ludin K, Mata J, Watt S, Lehmann E, Bähler J and Kohli J

    Institute of Cell Biology, University of Bern, 3012 Bern, Switzerland.

    Meiotic recombination arises from Rec12/Spo11-dependent formation of DNA double-strand breaks (DSBs) and their subsequent repair. We identified Rec12-binding peaks across the Schizosaccharomyces pombe genome using chromatin immunoprecipitation after reversible formaldehyde cross-linking combined with whole-genome DNA microarrays. Strong Rec12 binding coincided with previously identified DSBs at the recombination hotspots ura4A, mbs1, and mbs2 and correlated with DSB formation at a new site. In addition, Rec12 binding corresponded to eight novel conversion hotspots and correlated with crossover density in segments of chromosome I. Notably, Rec12 binding inversely correlated with guanine-cytosine (GC) content, contrary to findings in Saccharomyces cerevisiae. Although both replication origins and Rec12-binding sites preferred AT-rich gene-free regions, they seemed to exclude each other. We also uncovered a connection between binding sites of Rec12 and meiotic cohesin Rec8. Rec12-binding peaks lay often within 2.5 kb of a Rec8-binding peak. Rec12 binding showed preference for large intergenic regions and was found to bind preferentially near to genes expressed strongly in meiosis. Surprisingly, Rec12 binding was also detected in centromeric core regions, which raises the intriguing possibility that Rec12 plays additional roles in meiotic chromosome dynamics.

    Funded by: Cancer Research UK: A6517, C9546/A6517; Wellcome Trust: 077118

    Chromosoma 2008;117;5;431-44

  • MYO5B mutations cause microvillus inclusion disease and disrupt epithelial cell polarity.

    Müller T, Hess MW, Schiefermeier N, Pfaller K, Ebner HL, Heinz-Erian P, Ponstingl H, Partsch J, Röllinghoff B, Köhler H, Berger T, Lenhartz H, Schlenck B, Houwen RJ, Taylor CJ, Zoller H, Lechner S, Goulet O, Utermann G, Ruemmele FM, Huber LA and Janecke AR

    Department of Pediatrics II, Innsbruck Medical University, 6020 Innsbruck, Austria.

    Following homozygosity mapping in a single kindred, we identified nonsense and missense mutations in MYO5B, encoding type Vb myosin motor protein, in individuals with microvillus inclusion disease (MVID). MVID is characterized by lack of microvilli on the surface of enterocytes and occurrence of intracellular vacuolar structures containing microvilli. In addition, mislocalization of transferrin receptor in MVID enterocytes suggests that MYO5B deficiency causes defective trafficking of apical and basolateral proteins in MVID.

    Funded by: Austrian Science Fund FWF: P 19486

    Nature genetics 2008;40;10;1163-5

  • The neglected role of antibody in protection against bacteremia caused by nontyphoidal strains of Salmonella in African children.

    MacLennan CA, Gondwe EN, Msefula CL, Kingsley RA, Thomson NR, White SA, Goodall M, Pickard DJ, Graham SM, Dougan G, Hart CA, Molyneux ME and Drayson MT

    Malawi-Liverpool-Wellcome Trust Clinical Research Programme, College of Medicine, University of Malawi, Blantyre, Malawi.

    Nontyphoidal strains of Salmonella (NTS) are a common cause of bacteremia among African children. Cell-mediated immune responses control intracellular infection, but they do not protect against extracellular growth of NTS in the blood. We investigated whether antibody protects against NTS bacteremia in Malawian children, because we found this condition mainly occurs before 2 years of age, with relative sparing of infants younger than 4 months old. Sera from all healthy Malawian children tested aged more than 16 months contained anti-Salmonella antibody and successfully killed NTS. Killing was mediated by complement membrane attack complex and not augmented in the presence of blood leukocytes. Sera from most healthy children less than 16 months old lacked NTS-specific antibody, and sera lacking antibody did not kill NTS despite normal complement function. Addition of Salmonella-specific antibody, but not mannose-binding lectin, enabled NTS killing. All NTS strains tested had long-chain lipopolysaccharide and the rck gene, features that resist direct complement-mediated killing. Disruption of lipopolysaccharide biosynthesis enabled killing of NTS by serum lacking Salmonella-specific antibody. We conclude that Salmonella-specific antibody that overcomes the complement resistance of NTS develops by 2 years of life in Malawian children. This finding and the age-incidence of NTS bacteremia suggest that antibody protects against NTS bacteremia and support the development of vaccines against NTS that induce protective antibody.

    Funded by: Wellcome Trust: 067902

    The Journal of clinical investigation 2008;118;4;1553-62

  • Microbiology in the post-genomic era.

    Medini D, Serruto D, Parkhill J, Relman DA, Donati C, Moxon R, Falkow S and Rappuoli R

    Novartis Vaccines and Diagnostics, 53100 Siena, Italy.

    Genomics has revolutionized every aspect of microbiology. Now, 13 years after the first bacterial genome was sequenced, it is important to pause and consider what has changed in microbiology research as a consequence of genomics. In this article, we review the evolving field of bacterial typing and the genomic technologies that enable comparative analysis of multiple genomes and the metagenomes of complex microbial environments, and address the implications of the genomic era for the future of microbiology.

    Nature reviews. Microbiology 2008;6;6;419-30

  • Fission yeast SWI/SNF and RSC complexes show compositional and functional differences from budding yeast.

    Monahan BJ, Villén J, Marguerat S, Bähler J, Gygi SP and Winston F

    Department of Genetics, Harvard Medical School, Boston, Massachusetts 02115, USA.

    SWI/SNF chromatin-remodeling complexes have crucial roles in transcription and other chromatin-related processes. The analysis of the two members of this class in Saccharomyces cerevisiae, SWI/SNF and RSC, has heavily contributed to our understanding of these complexes. To understand the in vivo functions of SWI/SNF and RSC in an evolutionarily distant organism, we have characterized these complexes in Schizosaccharomyces pombe. Although core components are conserved between the two yeasts, the compositions of S. pombe SWI/SNF and RSC differ from their S. cerevisiae counterparts and in some ways are more similar to metazoan complexes. Furthermore, several of the conserved proteins, including actin-like proteins, are markedly different between the two yeasts with respect to their requirement for viability. Finally, phenotypic and microarray analyses identified widespread requirements for SWI/SNF and RSC on transcription including strong evidence that SWI/SNF directly represses iron-transport genes.

    Funded by: Cancer Research UK: A6517, C9546/A6517; NHGRI NIH HHS: HG3456, R01 HG003456, R01 HG003456-03; NIGMS NIH HHS: GM32967, R37 GM032967, R37 GM032967-25; Wellcome Trust: 077118

    Nature structural & molecular biology 2008;15;8;873-80

  • Inferring selection on amino acid preference in protein domains.

    Moses AM and Durbin R

    Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, UK.

    Models that explicitly account for the effect of selection on new mutations have been proposed to account for "codon bias" or the excess of "preferred" codons that results from selection for translational efficiency and/or accuracy. In principle, such models can be applied to any mutation that results in a preferred allele, but in most cases, the fitness effect of a specific mutation cannot be predicted. Here we show that it is possible to assign preferred and unpreferred states to amino acid changing mutations that occur in protein domains. We propose that mutations that lead to more common amino acids (at a given position in a domain) can be considered "preferred alleles" just as are synonymous mutations leading to codons for more abundant tRNAs. We use genome-scale polymorphism data to show that alleles for preferred amino acids in protein domains occur at higher frequencies in the population, as has been shown for preferred codons. We show that this effect is quantitative, such that there is a correlation between the shift in frequency of preferred alleles and the predicted fitness effect. As expected, we also observe a reduction in the numbers of polymorphisms and substitutions at more important positions in domains, consistent with stronger selection at those positions. We examine the derived allele frequency distribution and polymorphism to divergence ratios of preferred and unpreferred differences and find evidence for both negative and positive selections acting to maintain protein domains in the human population. Finally, we analyze a model for selection on amino acid preferences in protein domains and find that it is consistent with the quantitative effects that we observe.

    Funded by: Wellcome Trust: 077192

    Molecular biology and evolution 2008;26;3;527-36

  • Gene-wide analyses of genome-wide association data sets: evidence for multiple common risk alleles for schizophrenia and bipolar disorder and for overlap in genetic risk.

    Moskvina V, Craddock N, Holmans P, Nikolov I, Pahwa JS, Green E, Wellcome Trust Case Control Consortium, Owen MJ and O'Donovan MC

    Department of Psychological Medicine, School of Medicine, Cardiff University, Cardiff, UK.

    Genome-wide association (GWAS) analyses have identified susceptibility loci for many diseases, but most risk for any complex disorder remains unattributed. There is therefore scope for complementary approaches to these data sets. Gene-wide approaches potentially offer additional insights. They might identify association to genes through multiple signals. Also, by providing support for genes rather than single nucleotide polymorphisms (SNPs), they offer an additional opportunity to compare the results across data sets. We have undertaken gene-wide analysis of two GWAS data sets: schizophrenia and bipolar disorder. We performed two forms of analysis, one based on the smallest P-value per gene, the other on a truncated product of P method. For each data set and at a range of statistical thresholds, we observed significantly more SNPs within genes (P(min) for excess<0.001) showing evidence for association than expected whereas this was not true for extragenic SNPs (P(min) for excess>0.1). At a range of thresholds of significance, we also observed substantially more associated genes than expected (P(min) for excess in schizophrenia=1.8 x 10(-8), in bipolar=2.4 x 10(-6)). Moreover, an excess of genes showed evidence for association across disorders. Among those genes surpassing thresholds highly enriched for true association, we observed evidence for association to genes reported in other GWAS data sets (CACNA1C) or to closely related family members of those genes including CSF2RB, CACNA1B and DGKI. Our analyses show that association signals are enriched in and around genes, large numbers of genes contribute to both disorders and gene-wide analyses offer useful complementary approaches to more standard methods.

    Funded by: Medical Research Council: G0000934, G0800509, G0801418; Wellcome Trust: 076113, 079643

    Molecular psychiatry 2008;14;3;252-60

  • Dynamic instability of the major urinary protein gene family revealed by genomic and phenotypic comparisons between C57 and 129 strain mice.

    Mudge JM, Armstrong SD, McLaren K, Beynon RJ, Hurst JL, Nicholson C, Robertson DH, Wilming LG and Harrow JL

    Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge CB101SA, UK.

    Background: The major urinary proteins (MUPs) of Mus musculus domesticus are deposited in urine in large quantities, where they bind and release pheromones and also provide an individual 'recognition signal' via their phenotypic polymorphism. Whilst important information about MUP functionality has been gained in recent years, the gene cluster is poorly studied in terms of structure, genic polymorphism and evolution.

    Results: We combine targeted sequencing, manual genome annotation and phylogenetic analysis to compare the Mup clusters of C57BL/6J and 129 strains of mice. We describe organizational heterogeneity within both clusters: a central array of cassettes containing Mup genes highly similar at the protein level, flanked by regions containing Mup genes displaying significantly elevated divergence. Observed genomic rearrangements in all regions have likely been mediated by endogenous retroviral elements. Mup loci with coding sequences that differ between the strains are identified--including a gene/pseudogene pair--suggesting that these inbred lineages exhibit variation that exists in wild populations. We have characterized the distinct MUP profiles in the urine of both strains by mass spectrometry. The total MUP phenotype data is reconciled with our genomic sequence data, matching all proteins identified in urine to annotated genes.

    Conclusion: Our observations indicate that the MUP phenotypic polymorphism observed in wild populations results from a combination of Mup gene turnover coupled with currently unidentified mechanisms regulating gene expression patterns. We propose that the structural heterogeneity described within the cluster reflects functional divergence within the Mup gene family.

    Funded by: Biotechnology and Biological Sciences Research Council: S19816; NHGRI NIH HHS: U54 HG004555; Wellcome Trust: 077198

    Genome biology 2008;9;5;R91

  • Organ-specific requirements for Hdac1 in liver and pancreas formation.

    Noël ES, Casal-Sueiro A, Busch-Nentwich E, Verkade H, Dong PD, Stemple DL and Ober EA

    National Institute for Medical Research, Division of Developmental Biology, The Ridgeway, Mill Hill, London, UK.

    Liver, pancreas and lung originate from the presumptive foregut in temporal and spatial proximity. This requires precisely orchestrated transcriptional activation and repression of organ-specific gene expression within the same cell. Here, we show distinct roles for the chromatin remodelling factor and transcriptional repressor Histone deacetylase 1 (Hdac1) in endodermal organogenesis in zebrafish. Loss of Hdac1 causes defects in timely liver specification and in subsequent differentiation. Mosaic analyses reveal a cell-autonomous requirement for hdac1 within the hepatic endoderm. Our studies further reveal specific functions for Hdac1 in pancreas development. Loss of hdac1 causes the formation of ectopic endocrine clusters anteriorly to the main islet, as well as defects in exocrine pancreas specification and differentiation. In addition, we observe defects in extrahepatopancreatic duct formation and morphogenesis. Finally, loss of hdac1 results in an expansion of the foregut endoderm in the domain from which the liver and pancreas originate. Our genetic studies demonstrate that Hdac1 is crucial for regulating distinct steps in endodermal organogenesis. This suggests a model in which Hdac1 may directly or indirectly restrict foregut fates while promoting hepatic and exocrine pancreatic specification and differentiation, as well as pancreatic endocrine islet morphogenesis. These findings establish zebrafish as a tractable system to investigate chromatin remodelling factor functions in controlling gene expression programmes in vertebrate endodermal organogenesis.

    Funded by: Medical Research Council: MC_U117581329

    Developmental biology 2008;322;2;237-50

  • Mutations in mRNA export mediator GLE1 result in a fetal motoneuron disease.

    Nousiainen HO, Kestilä M, Pakkasjärvi N, Honkala H, Kuure S, Tallila J, Vuopala K, Ignatius J, Herva R and Peltonen L

    Department of Molecular Medicine, National Public Health Institute, Helsinki 00290, Finland.

    The most severe forms of motoneuron disease manifest in utero are characterized by marked atrophy of spinal cord motoneurons and fetal immobility. Here, we report that the defective gene underlying lethal motoneuron syndrome LCCS1 is the mRNA export mediator GLE1. Our finding of mutated GLE1 exposes a common pathway connecting the genes implicated in LCCS1, LCCS2 and LCCS3 and elucidates mRNA processing as a critical molecular mechanism in motoneuron development and maturation.

    Funded by: NIEHS NIH HHS: P01 ES11253-03; Wellcome Trust: 089061

    Nature genetics 2008;40;2;155-7

  • Post-genomic challenges for collaborative research in infectious diseases.

    Okeke IN and Wain J

    Department of Biology, Haverford College, Haverford, Pennsylvania 19041, USA.

    Although high-burden pathogens have been prioritized for sequencing, genomic research has yet to yield effective vaccines, diagnostics or therapeutics for the infectious diseases that burden developing countries. International research partnerships are needed more today than ever before, and we propose that increased participation by scientists in endemic areas would overcome current roadblocks and is an essential path towards translational research outcomes.

    Funded by: Wellcome Trust

    Nature reviews. Microbiology 2008;6;11;858-64

  • Genomic adaptation: a fungal perspective.

    Pain A and Hertz-Fowler C

    Nature reviews. Microbiology 2008;6;8;572-3

  • The genome of the simian and human malaria parasite Plasmodium knowlesi.

    Pain A, Böhme U, Berry AE, Mungall K, Finn RD, Jackson AP, Mourier T, Mistry J, Pasini EM, Aslett MA, Balasubrammaniam S, Borgwardt K, Brooks K, Carret C, Carver TJ, Cherevach I, Chillingworth T, Clark TG, Galinski MR, Hall N, Harper D, Harris D, Hauser H, Ivens A, Janssen CS, Keane T, Larke N, Lapp S, Marti M, Moule S, Meyer IM, Ormond D, Peters N, Sanders M, Sanders S, Sargeant TJ, Simmonds M, Smith F, Squares R, Thurston S, Tivey AR, Walker D, White B, Zuiderwijk E, Churcher C, Quail MA, Cowman AF, Turner CM, Rajandream MA, Kocken CH, Thomas AW, Newbold CI, Barrell BG and Berriman M

    Wellcome Trust Sanger Institute, Genome Campus, Hinxton, Cambridgeshire CB10 1SA, UK.

    Plasmodium knowlesi is an intracellular malaria parasite whose natural vertebrate host is Macaca fascicularis (the 'kra' monkey); however, it is now increasingly recognized as a significant cause of human malaria, particularly in southeast Asia. Plasmodium knowlesi was the first malaria parasite species in which antigenic variation was demonstrated, and it has a close phylogenetic relationship to Plasmodium vivax, the second most important species of human malaria parasite (reviewed in ref. 4). Despite their relatedness, there are important phenotypic differences between them, such as host blood cell preference, absence of a dormant liver stage or 'hypnozoite' in P. knowlesi, and length of the asexual cycle (reviewed in ref. 4). Here we present an analysis of the P. knowlesi (H strain, Pk1(A+) clone) nuclear genome sequence. This is the first monkey malaria parasite genome to be described, and it provides an opportunity for comparison with the recently completed P. vivax genome and other sequenced Plasmodium genomes. In contrast to other Plasmodium genomes, putative variant antigen families are dispersed throughout the genome and are associated with intrachromosomal telomere repeats. One of these families, the KIRs, contains sequences that collectively match over one-half of the host CD99 extracellular domain, which may represent an unusual form of molecular mimicry.

    Funded by: Wellcome Trust: 085775

    Nature 2008;455;7214;799-803

  • Complete genome sequence of uropathogenic Proteus mirabilis, a master of both adherence and motility.

    Pearson MM, Sebaihia M, Churcher C, Quail MA, Seshasayee AS, Luscombe NM, Abdellah Z, Arrosmith C, Atkin B, Chillingworth T, Hauser H, Jagels K, Moule S, Mungall K, Norbertczak H, Rabbinowitsch E, Walker D, Whithead S, Thomson NR, Rather PN, Parkhill J and Mobley HL

    Department of Microbiology and Immunology, University of Michigan Medical School, Ann Arbor, MI 48109-0620, USA.

    The gram-negative enteric bacterium Proteus mirabilis is a frequent cause of urinary tract infections in individuals with long-term indwelling catheters or with complicated urinary tracts (e.g., due to spinal cord injury or anatomic abnormality). P. mirabilis bacteriuria may lead to acute pyelonephritis, fever, and bacteremia. Most notoriously, this pathogen uses urease to catalyze the formation of kidney and bladder stones or to encrust or obstruct indwelling urinary catheters. Here we report the complete genome sequence of P. mirabilis HI4320, a representative strain cultured in our laboratory from the urine of a nursing home patient with a long-term (> or =30 days) indwelling urinary catheter. The genome is 4.063 Mb long and has a G+C content of 38.88%. There is a single plasmid consisting of 36,289 nucleotides. Annotation of the genome identified 3,685 coding sequences and seven rRNA loci. Analysis of the sequence confirmed the presence of previously identified virulence determinants, as well as a contiguous 54-kb flagellar regulon and 17 types of fimbriae. Genes encoding a potential type III secretion system were identified on a low-G+C-content genomic island containing 24 intact genes that appear to encode all components necessary to assemble a type III secretion system needle complex. In addition, the P. mirabilis HI4320 genome possesses four tandem copies of the zapE metalloprotease gene, genes encoding six putative autotransporters, an extension of the atf fimbrial operon to six genes, including an mrpJ homolog, and genes encoding at least five iron uptake mechanisms, two potential type IV secretion systems, and 16 two-component regulators.

    Funded by: NIAID NIH HHS: AI059722, F32 AI068324, F32 AI068324-01A2, R01 AI059722, T32 AI007528, T32 AI7528; Wellcome Trust

    Journal of bacteriology 2008;190;11;4027-37

  • Copy number variation and evolution in humans and chimpanzees.

    Perry GH, Yang F, Marques-Bonet T, Murphy C, Fitzgerald T, Lee AS, Hyland C, Stone AC, Hurles ME, Tyler-Smith C, Eichler EE, Carter NP, Lee C and Redon R

    School of Human Evolution & Social Change, Arizona State University, Tempe, Arizona 85287, USA.

    Copy number variants (CNVs) underlie many aspects of human phenotypic diversity and provide the raw material for gene duplication and gene family expansion. However, our understanding of their evolutionary significance remains limited. We performed comparative genomic hybridization on a single human microarray platform to identify CNVs among the genomes of 30 humans and 30 chimpanzees as well as fixed copy number differences between species. We found that human and chimpanzee CNVs occur in orthologous genomic regions far more often than expected by chance and are strongly associated with the presence of highly homologous intrachromosomal segmental duplications. By adapting population genetic analyses for use with copy number data, we identified functional categories of genes that have likely evolved under purifying or positive selection for copy number changes. In particular, duplications and deletions of genes with inflammatory response and cell proliferation functions may have been fixed by positive selection and involved in the adaptive phenotypic differentiation of humans and chimpanzees.

    Funded by: Howard Hughes Medical Institute; NCRR NIH HHS: C06 RR014491, C06 RR016483, RR014491, RR015087, RR016483, U42 RR015087; NHGRI NIH HHS: HG004221, P41 HG004221; Wellcome Trust

    Genome research 2008;18;11;1698-710

  • Molecular characterization of the Salmonella enterica serovar Typhi Vi-typing bacteriophage E1.

    Pickard D, Thomson NR, Baker S, Wain J, Pardo M, Goulding D, Hamlin N, Choudhary J, Threfall J and Dougan G

    Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge CB10 1SA, United Kingdom.

    Some bacteriophages target potentially pathogenic bacteria by exploiting surface-associated virulence factors as receptors. For example, phage have been identified that exhibit specificity for Vi capsule producing Salmonella enterica serovar Typhi. Here we have characterized the Vi-associated E1-typing bacteriophage using a number of molecular approaches. The absolute requirement for Vi capsule expression for infectivity was demonstrated using different Vi-negative S. enterica derivatives. The phage particles were shown to have an icosahedral head and a long noncontractile tail structure. The genome is 45,362 bp in length with defined capsid and tail regions that exhibit significant homology to the S. enterica transducing phage ES18. Mass spectrometry was used to confirm the presence of a number of hypothetical proteins in the Vi phage E1 particle and demonstrate that a number of phage proteins are modified posttranslationally. The genome of the Vi phage E1 is significantly related to other bacteriophages belonging to the same serovar Typhi phage-typing set, and we demonstrate a role for phage DNA modification in determining host specificity.

    Funded by: Wellcome Trust

    Journal of bacteriology 2008;190;7;2580-7

  • BAC TransgeneOmics: a high-throughput method for exploration of protein function in mammals.

    Poser I, Sarov M, Hutchins JR, Hériché JK, Toyoda Y, Pozniakovsky A, Weigl D, Nitzsche A, Hegemann B, Bird AW, Pelletier L, Kittler R, Hua S, Naumann R, Augsburg M, Sykora MM, Hofemeister H, Zhang Y, Nasmyth K, White KP, Dietzel S, Mechtler K, Durbin R, Stewart AF, Peters JM, Buchholz F and Hyman AA

    Max Planck Institute for Molecular Cell Biology and Genetics, Pfotenhauerstrasse 108, D-01307 Dresden, Germany.

    The interpretation of genome sequences requires reliable and standardized methods to assess protein function at high throughput. Here we describe a fast and reliable pipeline to study protein function in mammalian cells based on protein tagging in bacterial artificial chromosomes (BACs). The large size of the BAC transgenes ensures the presence of most, if not all, regulatory elements and results in expression that closely matches that of the endogenous gene. We show that BAC transgenes can be rapidly and reliably generated using 96-well-format recombineering. After stable transfection of these transgenes into human tissue culture cells or mouse embryonic stem cells, the localization, protein-protein and/or protein-DNA interactions of the tagged protein are studied using generic, tag-based assays. The same high-throughput approach will be generally applicable to other model systems.

    Funded by: NHGRI NIH HHS: 1R01HG004428-01; Wellcome Trust: 077192

    Nature methods 2008;5;5;409-15

  • Variants in MTNR1B influence fasting glucose levels.

    Prokopenko I, Langenberg C, Florez JC, Saxena R, Soranzo N, Thorleifsson G, Loos RJ, Manning AK, Jackson AU, Aulchenko Y, Potter SC, Erdos MR, Sanna S, Hottenga JJ, Wheeler E, Kaakinen M, Lyssenko V, Chen WM, Ahmadi K, Beckmann JS, Bergman RN, Bochud M, Bonnycastle LL, Buchanan TA, Cao A, Cervino A, Coin L, Collins FS, Crisponi L, de Geus EJ, Dehghan A, Deloukas P, Doney AS, Elliott P, Freimer N, Gateva V, Herder C, Hofman A, Hughes TE, Hunt S, Illig T, Inouye M, Isomaa B, Johnson T, Kong A, Krestyaninova M, Kuusisto J, Laakso M, Lim N, Lindblad U, Lindgren CM, McCann OT, Mohlke KL, Morris AD, Naitza S, Orrù M, Palmer CN, Pouta A, Randall J, Rathmann W, Saramies J, Scheet P, Scott LJ, Scuteri A, Sharp S, Sijbrands E, Smit JH, Song K, Steinthorsdottir V, Stringham HM, Tuomi T, Tuomilehto J, Uitterlinden AG, Voight BF, Waterworth D, Wichmann HE, Willemsen G, Witteman JC, Yuan X, Zhao JH, Zeggini E, Schlessinger D, Sandhu M, Boomsma DI, Uda M, Spector TD, Penninx BW, Altshuler D, Vollenweider P, Jarvelin MR, Lakatta E, Waeber G, Fox CS, Peltonen L, Groop LC, Mooser V, Cupples LA, Thorsteinsdottir U, Boehnke M, Barroso I, Van Duijn C, Dupuis J, Watanabe RM, Stefansson K, McCarthy MI, Wareham NJ, Meigs JB and Abecasis GR

    [1] Oxford Centre for Diabetes, Endocrinology and Metabolism, University of Oxford, Oxford OX3 7LJ, UK. [2] Wellcome Trust Centre for Human Genetics, University of Oxford, Oxford OX3 7BN, UK. [3] These authors contributed equally to this work.

    To identify previously unknown genetic loci associated with fasting glucose concentrations, we examined the leading association signals in ten genome-wide association scans involving a total of 36,610 individuals of European descent. Variants in the gene encoding melatonin receptor 1B (MTNR1B) were consistently associated with fasting glucose across all ten studies. The strongest signal was observed at rs10830963, where each G allele (frequency 0.30 in HapMap CEU) was associated with an increase of 0.07 (95% CI = 0.06-0.08) mmol/l in fasting glucose levels (P = 3.2 x 10(-50)) and reduced beta-cell function as measured by homeostasis model assessment (HOMA-B, P = 1.1 x 10(-15)). The same allele was associated with an increased risk of type 2 diabetes (odds ratio = 1.09 (1.05-1.12), per G allele P = 3.3 x 10(-7)) in a meta-analysis of 13 case-control studies totaling 18,236 cases and 64,453 controls. Our analyses also confirm previous associations of fasting glucose with variants at the G6PC2 (rs560887, P = 1.1 x 10(-57)) and GCK (rs4607517, P = 1.0 x 10(-25)) loci.

    Funded by: Intramural NIH HHS: Z01 HG000024; Medical Research Council: G0000649, G016121, G0500539, G0601261, G0701863, MC_U106179471, MC_U106188470; NCRR NIH HHS: RR-163736; NHGRI NIH HHS: HG-02651, R01 HG002651, R01 HG002651-05; NHLBI NIH HHS: HC-25195, HL-084729, HL-087679, N01 HC025195, N01HC25195, N02-HL-6-4278, R01 HL087679, R01 HL087679-02, U01 HL084729, U01 HL084729-03; NIDA NIH HHS: DA-021519, U54 DA021519, U54 DA021519-04; NIDDK NIH HHS: DK-062370, DK-065978, DK-072193, DK-078616, DK-080140, DK069922, K23 DK065978, K23 DK065978-05, K24 DK080140, K24 DK080140-01, K24 DK080140-02, R01 DK029867, R01 DK062370, R01 DK062370-05, R01 DK069922, R01 DK069922-02, R01 DK072193, R01 DK072193-04, R01 DK078616, R01 DK078616-01A1, R56 DK062370, U01 DK062370, U01 DK078616; NIMH NIH HHS: MH059160, R01 MH059160-04; Wellcome Trust: 076113, 077011, 077016, 079557, 083948, 089061, 090532, GR069224, GR072960

    Nature genetics 2008;41;1;77-81

  • A large genome center's improvements to the Illumina sequencing system.

    Quail MA, Kozarewa I, Smith F, Scally A, Stephens PJ, Durbin R, Swerdlow H and Turner DJ

    Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridgeshire, CB10 1SA, UK.

    The Wellcome Trust Sanger Institute is one of the world's largest genome centers, and a substantial amount of our sequencing is performed with 'next-generation' massively parallel sequencing technologies: in June 2008 the quantity of purity-filtered sequence data generated by our Genome Analyzer (Illumina) platforms reached 1 terabase, and our average weekly Illumina production output is currently 64 gigabases. Here we describe a set of improvements we have made to the standard Illumina protocols to make the library preparation more reliable in a high-throughput environment, to reduce bias, tighten insert size distribution and reliably obtain high yields of data.

    Funded by: Medical Research Council: G0701805; Wellcome Trust: 079643

    Nature methods 2008;5;12;1005-10

  • An integrated resource for genome-wide identification and analysis of human tissue-specific differentially methylated regions (tDMRs).

    Rakyan VK, Down TA, Thorne NP, Flicek P, Kulesha E, Gräf S, Tomazou EM, Bäckdahl L, Johnson N, Herberth M, Howe KL, Jackson DK, Miretti MM, Fiegler H, Marioni JC, Birney E, Hubbard TJ, Carter NP, Tavaré S and Beck S

    Institute of Cell and Molecular Science, Barts and the London, London E1 2AT, United Kingdom.

    We report a novel resource (methylation profiles of DNA, or mPod) for human genome-wide tissue-specific DNA methylation profiles. mPod consists of three fully integrated parts, genome-wide DNA methylation reference profiles of 13 normal somatic tissues, placenta, sperm, and an immortalized cell line, a visualization tool that has been integrated with the Ensembl genome browser and a new algorithm for the analysis of immunoprecipitation-based DNA methylation profiles. We demonstrate the utility of our resource by identifying the first comprehensive genome-wide set of tissue-specific differentially methylated regions (tDMRs) that may play a role in cellular identity and the regulation of tissue-specific genome function. We also discuss the implications of our findings with respect to the regulatory potential of regions with varied CpG density, gene expression, transcription factor motifs, gene ontology, and correlation with other epigenetic marks such as histone modifications.

    Funded by: Cancer Research UK: C14303/A4646; Wellcome Trust: 077198

    Genome research 2008;18;9;1518-29

  • The Protein Feature Ontology: a tool for the unification of protein feature annotations.

    Reeves GA, Eilbeck K, Magrane M, O'Donovan C, Montecchi-Palazzi L, Harris MA, Orchard S, Jimenez RC, Prlic A, Hubbard TJ, Hermjakob H and Thornton JM

    EMBL-European Bioinformatics Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SD, UK.

    Motivation: The advent of sequencing and structural genomics projects has provided a dramatic boost in the number of uncharacterized protein structures and sequences. Consequently, many computational tools have been developed to help elucidate protein function. However, such services are spread throughout the world, often with standalone web pages. Integration of these methods is needed and so far this has not been possible as there was no common vocabulary available that could be used as a standard language.

    Results: The Protein Feature Ontology has been developed to provide a structured controlled vocabulary for features on a protein sequence or structure and comprises approximately 100 positional terms, now integrated into the Sequence Ontology (SO) and 40 non-positional terms which describe features relating to the whole-protein sequence. In addition, post-translational modifications are described by using a pre-existing ontology, the Protein Modification Ontology (MOD). This ontology is being used to integrate over 150 distinct annotations provided by the BioSapiens Network of Excellence, a consortium comprising 19 partner sites in Europe.

    Availability: The Protein Feature Ontology can be browsed by accessing the ontology lookup service at the European Bioinformatics Institute (

    Funded by: Wellcome Trust: 062023, 077198

    Bioinformatics (Oxford, England) 2008;24;23;2767-72

  • Fission yeast MAP kinase Sty1 is recruited to stress-induced genes.

    Reiter W, Watt S, Dawson K, Lawrence CL, Bähler J, Jones N and Wilkinson CR

    Paterson Institute for Cancer Research, University of Manchester, Wilmslow Road, Manchester, UK.

    The stress-induced expression of many fission yeast genes is dependent upon the Sty1 mitogen-activated protein kinase (MAPK) and Atf1 transcription factor. Atf1 is phosphorylated by Sty1 yet this phosphorylation is not required for stress-induced gene expression, suggesting another mechanism exists whereby Sty1 activates transcription. Here we show that Sty1 associates with Atf1-dependent genes and is recruited to both their promoters and coding regions. This occurs in response to various stress conditions coincident with the kinetics of the activation of Sty1. Association with promoters is not a consequence of increased nuclear accumulation of Sty1 nor does it require the phosphorylation of Atf1. However, recruitment is completely abolished in a mutant lacking Sty1 kinase activity. Both Atf1 and its binding partner Pcr1 are required for association of Sty1 with Atf1-dependent promoters, suggesting that this heterodimer must be intact for optimal recruitment of the MAPK. However, many Atf1-dependent genes are still expressed in a pcr1Delta mutant but with significantly delayed kinetics, thus providing an explanation for the relatively mild stress sensitivity displayed by pcr1Delta. Consistent with this delay, Sty1 and Atf1 cannot be detected at these promoters in this condition, suggesting that their association with chromatin is weak or transient in the absence of Pcr1.

    Funded by: Cancer Research UK: A6517; Wellcome Trust: 077118, 098051

    The Journal of biological chemistry 2008;283;15;9945-56

  • Use and misuse of the gene ontology annotations.

    Rhee SY, Wood V, Dolinski K and Draghici S

    Carnegie Institution for Science, Department of Plant Biology, 260 Panama Street, Stanford, California 94305, USA.

    The Gene Ontology (GO) project is a collaboration among model organism databases to describe gene products from all organisms using a consistent and computable language. GO produces sets of explicitly defined, structured vocabularies that describe biological processes, molecular functions and cellular components of gene products in both a computer- and human-readable manner. Here we describe key aspects of GO, which, when overlooked, can cause erroneous results, and address how these pitfalls can be avoided.

    Nature reviews. Genetics 2008;9;7;509-15

  • Male-pattern baldness susceptibility locus at 20p11.

    Richards JB, Yuan X, Geller F, Waterworth D, Bataille V, Glass D, Song K, Waeber G, Vollenweider P, Aben KK, Kiemeney LA, Walters B, Soranzo N, Thorsteinsdottir U, Kong A, Rafnar T, Deloukas P, Sulem P, Stefansson H, Stefansson K, Spector TD and Mooser V

    Department of Twin Research and Genetic Epidemiology, King's College London, London SE1 7EH, UK.

    We conducted a genome-wide association study for androgenic alopecia in 1,125 men and identified a newly associated locus at chromosome 20p11.22, confirmed in three independent cohorts (n = 1,650; OR = 1.60, P = 1.1 x 10(-14) for rs1160312). The one man in seven who harbors risk alleles at both 20p11.22 and AR (encoding the androgen receptor) has a sevenfold-increased odds of androgenic alopecia (OR = 7.12, P = 3.7 x 10(-15)).

    Funded by: Wellcome Trust: 077011

    Nature genetics 2008;40;11;1282-4

  • Maximum-likelihood estimation of site-specific mutation rates in human mitochondrial DNA from partial phylogenetic classification.

    Rosset S, Wells RS, Soria-Hernanz DF, Tyler-Smith C, Royyuru AK, Behar DM and Genographic Consortium

    Department of Statistics and Operations Research, Tel Aviv University, Tel Aviv, Israel.

    The mitochondrial DNA hypervariable segment I (HVS-I) is widely used in studies of human evolutionary genetics, and therefore accurate estimates of mutation rates among nucleotide sites in this region are essential. We have developed a novel maximum-likelihood methodology for estimating site-specific mutation rates from partial phylogenetic information, such as haplogroup association. The resulting estimation problem is a generalized linear model, with a nonstandard link function. We develop inference and bias correction tools for our estimates and a hypothesis-testing approach for site independence. We demonstrate our methodology using 16,609 HVS-I samples from the Genographic Project. Our results suggest that mutation rates among nucleotide sites in HVS-I are highly variable. The 16,400-16,500 region exhibits significantly lower rates compared to other regions, suggesting potential functional constraints. Several loci identified in the literature as possible termination-associated sequences (TAS) do not yield statistically slower rates than the rest of HVS-I, casting doubt on their functional importance. Our tests do not reject the null hypothesis of independent mutation rates among nucleotide sites, supporting the use of site-independence assumption for analyzing HVS-I. Potential extensions of our methodology include its application to estimation of mutation rates in other genetic regions, like Y chromosome short tandem repeats.

    Funded by: Wellcome Trust

    Genetics 2008;180;3;1511-24

  • TreeFam: 2008 Update.

    Ruan J, Li H, Chen Z, Coghlan A, Coin LJ, Guo Y, Hériché JK, Hu Y, Kristiansen K, Li R, Liu T, Moses A, Qin J, Vang S, Vilella AJ, Ureta-Vidal A, Bolund L, Wang J and Durbin R

    Beijing Institute of Genomics of the Chinese Academy of Sciences, Beijing Genomics Institute, Beijing 101300, China.

    TreeFam ( was developed to provide curated phylogenetic trees for all animal gene families, as well as orthologue and paralogue assignments. Release 4.0 of TreeFam contains curated trees for 1314 families and automatically generated trees for another 14,351 families. We have expanded TreeFam to include 25 fully sequenced animal genomes, as well as four genomes from plant and fungal outgroup species. We have also introduced more accurate approaches for automatically grouping genes into families, for building phylogenetic trees, and for inferring orthologues and paralogues. The user interface for viewing phylogenetic trees and family information has been improved. Furthermore, a new perl API lets users easily extract data from the TreeFam mysql database.

    Funded by: Wellcome Trust

    Nucleic acids research 2008;36;Database issue;D735-40

  • Presence of interstereocilial links in waltzer mutants suggests Cdh23 is not essential for tip link formation.

    Rzadzinska AK and Steel KP

    Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SA, United Kingdom.

    Cadherin23 has been proposed to form the upper part of the tip link, an interstereocilial link believed to control opening of transducer channels of sensory hair cells. However, we detect tip link-like links in mouse mutants with null alleles of Cdh23, suggesting the presence of other components that permit formation of a link between the tip of one stereocilium and the side of the adjacent taller stereocilium.

    Funded by: Medical Research Council: G0300212, MC_QA137918; Wellcome Trust

    Neuroscience 2008;158;2;365-8

  • Genome-wide association analysis of metabolic traits in a birth cohort from a founder population.

    Sabatti C, Service SK, Hartikainen AL, Pouta A, Ripatti S, Brodsky J, Jones CG, Zaitlen NA, Varilo T, Kaakinen M, Sovio U, Ruokonen A, Laitinen J, Jakkula E, Coin L, Hoggart C, Collins A, Turunen H, Gabriel S, Elliot P, McCarthy MI, Daly MJ, Järvelin MR, Freimer NB and Peltonen L

    Department of Human Genetics and Los Angeles, Los Angeles, California 90095, USA.

    Genome-wide association studies (GWAS) of longitudinal birth cohorts enable joint investigation of environmental and genetic influences on complex traits. We report GWAS results for nine quantitative metabolic traits (triglycerides, high-density lipoprotein, low-density lipoprotein, glucose, insulin, C-reactive protein, body mass index, and systolic and diastolic blood pressure) in the Northern Finland Birth Cohort 1966 (NFBC1966), drawn from the most genetically isolated Finnish regions. We replicate most previously reported associations for these traits and identify nine new associations, several of which highlight genes with metabolic functions: high-density lipoprotein with NR1H3 (LXRA), low-density lipoprotein with AR and FADS1-FADS2, glucose with MTNR1B, and insulin with PANK1. Two of these new associations emerged after adjustment of results for body mass index. Gene-environment interaction analyses suggested additional associations, which will require validation in larger samples. The currently identified loci, together with quantified environmental exposures, explain little of the trait variation in NFBC1966. The association observed between low-density lipoprotein and an infrequent variant in AR suggests the potential of such a cohort for identifying associations with both common, low-impact and rarer, high-impact quantitative trait loci.

    Funded by: NCRR NIH HHS: U54 RR020278; NHGRI NIH HHS: T32 HG002536; NHLBI NIH HHS: R01 HL087679; NIGMS NIH HHS: GM053275-14, R01 GM053275; NIMH NIH HHS: MH083268, RL1 MH083268; Wellcome Trust: 089061

    Nature genetics 2008;41;1;35-46

  • Pfam 10 years on: 10,000 families and still growing.

    Sammut SJ, Finn RD and Bateman A

    Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton Hall, Hinxton, Cambridgeshire, CB10 1SA, UK.

    Classifications of proteins into groups of related sequences are in some respects like a periodic table for biology, allowing us to understand the underlying molecular biology of any organism. Pfam is a large collection of protein domains and families. Its scientific goal is to provide a complete and accurate classification of protein families and domains. The next release of the database will contain over 10,000 entries, which leads us to reflect on how far we are from completing this work. Currently Pfam matches 72% of known protein sequences, but for proteins with known structure Pfam matches 95%, which we believe represents the likely upper bound. Based on our analysis a further 28,000 families would be required to achieve this level of coverage for the current sequence database. We also show that as more sequences are added to the sequence databases the fraction of sequences that Pfam matches is reduced, suggesting that continued addition of new families is essential to maintain its relevance.

    Funded by: Wellcome Trust: 087656

    Briefings in bioinformatics 2008;9;3;210-9

  • Mendelian randomisation studies of type 2 diabetes: future prospects.

    Sandhu MS, Debenham SL, Barroso I and Loos RJ

    MRC Epidemiology Unit, Strangeways Research Laboratory, Cambridge, UK.

    Funded by: Medical Research Council: MC_U106179471, MC_U106188470; Wellcome Trust: 077016

    Diabetologia 2008;51;2;211-3

  • LDL-cholesterol concentrations: a genome-wide association study.

    Sandhu MS, Waterworth DM, Debenham SL, Wheeler E, Papadakis K, Zhao JH, Song K, Yuan X, Johnson T, Ashford S, Inouye M, Luben R, Sims M, Hadley D, McArdle W, Barter P, Kesäniemi YA, Mahley RW, McPherson R, Grundy SM, Wellcome Trust Case Control Consortium, Bingham SA, Khaw KT, Loos RJ, Waeber G, Barroso I, Strachan DP, Deloukas P, Vollenweider P, Wareham NJ and Mooser V

    Department of Public Health and Primary Care, Strangeways Research Laboratory, University of Cambridge, Cambridge, UK.

    Background: LDL cholesterol has a causal role in the development of cardiovascular disease. Improved understanding of the biological mechanisms that underlie the metabolism and regulation of LDL cholesterol might help to identify novel therapeutic targets. We therefore did a genome-wide association study of LDL-cholesterol concentrations.

    Methods: We used genome-wide association data from up to 11,685 participants with measures of circulating LDL-cholesterol concentrations across five studies, including data for 293 461 autosomal single nucleotide polymorphisms (SNPs) with a minor allele frequency of 5% or more that passed our quality control criteria. We also used data from a second genome-wide array in up to 4337 participants from three of these five studies, with data for 290,140 SNPs. We did replication studies in two independent populations consisting of up to 4979 participants. Statistical approaches, including meta-analysis and linkage disequilibrium plots, were used to refine association signals; we analysed pooled data from all seven populations to determine the effect of each SNP on variations in circulating LDL-cholesterol concentrations.

    Findings: In our initial scan, we found two SNPs (rs599839 [p=1.7x10(-15)] and rs4970834 [p=3.0x10(-11)]) that showed genome-wide statistical association with LDL cholesterol at chromosomal locus 1p13.3. The second genome screen found a third statistically associated SNP at the same locus (rs646776 [p=4.3x10(-9)]). Meta-analysis of data from all studies showed an association of SNPs rs599839 (combined p=1.2x10(-33)) and rs646776 (p=4.8x10(-20)) with LDL-cholesterol concentrations. SNPs rs599839 and rs646776 both explained around 1% of the variation in circulating LDL-cholesterol concentrations and were associated with about 15% of an SD change in LDL cholesterol per allele, assuming an SD of 1 mmol/L.

    Interpretation: We found evidence for a novel locus for LDL cholesterol on chromosome 1p13.3. These results potentially provide insight into the biological mechanisms that underlie the regulation of LDL cholesterol and might help in the discovery of novel therapeutic targets for cardiovascular disease.

    Funded by: Medical Research Council: G0000934, G0701863, MC_QA137934, MC_U105630924, MC_U106188470; Wellcome Trust: 068545/Z/02

    Lancet (London, England) 2008;371;9611;483-91

  • New case of interstitial deletion 12(q15-q21.2) in a girl with facial dysmorphism and mental retardation.

    Schluth C, Gesny R, Borck G, Redon R, Abadie V, Kleinfinger P, Munnich A, Lyonnet S and Colleaux L

    Département de Génétique, Hôpital Necker Enfants Malades, Paris, France.

    Interstitial deletions of the long arm of chromosome 12 are rare rearrangements with only 15 cases reported in the literature. The phenotype may include facial dysmorphism, developmental delay, ectodermal abnormalities, cardiac and renal malformations depending on breakpoints' position. Here, we describe a third case of 12(q15-q21.2) deletion ascertained through CGH-array analyses and provide a 5-year follow-up. The patient presented with pre- and postnatal growth retardation, congenital heart defect, developmental delay, and facial dysmorphism changing with age, underlining the importance of long-term follow-up. We compared this new case with previous observations of 12q deletions in order to propose phenotype-karyotype correlations.

    American journal of medical genetics. Part A 2008;146A;1;93-6

  • Repeated replication and a prospective meta-analysis of the association between chromosome 9p21.3 and coronary artery disease.

    Schunkert H, Götz A, Braund P, McGinnis R, Tregouet DA, Mangino M, Linsel-Nitschke P, Cambien F, Hengstenberg C, Stark K, Blankenberg S, Tiret L, Ducimetiere P, Keniry A, Ghori MJ, Schreiber S, El Mokhtari NE, Hall AS, Dixon RJ, Goodall AH, Liptau H, Pollard H, Schwarz DF, Hothorn LA, Wichmann HE, König IR, Fischer M, Meisinger C, Ouwehand W, Deloukas P, Thompson JR, Erdmann J, Ziegler A, Samani NJ and Cardiogenics Consortium

    Medizinische Klinik II, Universität zu Lübeck, Lübeck, Germany.

    Background: Recently, genome-wide association studies identified variants on chromosome 9p21.3 as affecting the risk of coronary artery disease (CAD). We investigated the association of this locus with CAD in 7 case-control studies and undertook a meta-analysis.

    Methods and results: A single-nucleotide polymorphism (SNP), rs1333049, representing the 9p21.3 locus, was genotyped in 7 case-control studies involving a total of 4645 patients with myocardial infarction or CAD and 5177 controls. The mode of inheritance was determined. In addition, in 5 of the 7 studies, we genotyped 3 additional SNPs to assess a risk-associated haplotype (ACAC). Finally, a meta-analysis of the present data and previously published samples was conducted. A limited fine mapping of the locus was performed. The risk allele (C) of the lead SNP, rs1333049, was uniformly associated with CAD in each study (P<0.05). In a pooled analysis, the odds ratio per copy of the risk allele was 1.29 (95% confidence interval, 1.22 to 1.37; P=0.0001). Haplotype analysis further suggested that this effect was not homogeneous across the haplotypic background (test for interaction, P=0.0079). An autosomal-additive mode of inheritance best explained the underlying association. The meta-analysis of the rs1333049 SNP in 12,004 cases and 28,949 controls increased the overall level of evidence for association with CAD to P=6.04x10(-10) (odds ratio, 1.24; 95% confidence interval, 1.20 to 1.29). Genotyping of 31 additional SNPs in the region identified several with a highly significant association with CAD, but none had predictive information beyond that of the rs1333049 SNP.

    Conclusions: This broad replication provides unprecedented evidence for association between genetic variants at chromosome 9p21.3 and risk of CAD.

    Funded by: British Heart Foundation; Wellcome Trust: 077011, 082371

    Circulation 2008;117;13;1675-84

  • Protein interactions in human genetic diseases.

    Schuster-Böckler B and Bateman A

    Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, CB10 1SA, UK.

    We present a novel method that combines protein structure information with protein interaction data to identify residues that form part of an interaction interface. Our prediction method can retrieve interaction hotspots with an accuracy of 60% (at a 20% false positive rate). The method was applied to all mutations in the Online Mendelian Inheritance in Man (OMIM) database, predicting 1,428 mutations to be related to an interaction defect. Combining predicted and hand-curated sets, we discuss how mutations affect protein interactions in general.

    Funded by: Wellcome Trust: 087656

    Genome biology 2008;9;1;R9

  • Robust, persistent transgene expression in human embryonic stem cells is achieved with AAVS1-targeted integration.

    Smith JR, Maguire S, Davis LA, Alexander M, Yang F, Chandran S, ffrench-Constant C and Pedersen RA

    Department of Surgery, Cambridge Institute for Medical Research, Wellcome Trust/MRC Building, Addenbrooke's Hospital, Hills Road, Cambridge CB2 OXY, United Kingdom.

    Silencing and variegated transgene expression are poorly understood problems that can interfere with gene function studies in human embryonic stem cells (hESCs). We show that transgene expression (enhanced green fluorescent protein [EGFP]) from random integration sites in hESCs is affected by variegation and silencing, with only half of hESCs expressing the transgene, which is gradually lost after withdrawal of selection and differentiation. We tested the hypothesis that a transgene integrated into the adeno-associated virus type 2 (AAV2) target region on chromosome 19, known as the AAVS1 locus, would maintain transgene expression in hESCs. When we used AAV2 technology to target the AAVS1 locus, 4.16% of hESC clones achieved AAVS1-targeted integration. Targeted clones expressed Oct-4, stage-specific embryonic antigen-3 (SSEA3), and Tra-1-60 and differentiated into all three primary germ layers. EGFP expression from the AAVS1 locus showed significantly reduced variegated expression when in selection, with 90% +/- 4% of cells expressing EGFP compared with 57% +/- 32% for randomly integrated controls, and reduced tendency to undergo silencing, with 86% +/- 7% hESCs expressing EGFP 25 days after withdrawal of selection compared with 39% +/- 31% for randomly integrated clones. In addition, quantitative polymerase chain reaction analysis of hESCs also indicated significantly higher levels of EGFP mRNA in AAVS1-targeted clones as compared with randomly integrated clones. Transgene expression from the AAVS1 locus was shown to be stable during hESC differentiation, with more than 90% of cells expressing EGFP after 15 days of differentiation, as compared with approximately 30% for randomly integrated clones. These results demonstrate the utility of transgene integration at the AAVS1 locus in hESCs and its potential clinical application.

    Funded by: Medical Research Council: G0300300, G0300723, G0600275

    Stem cells (Dayton, Ohio) 2008;26;2;496-504

  • Conservation of the H19 noncoding RNA and H19-IGF2 imprinting mechanism in therians.

    Smits G, Mungall AJ, Griffiths-Jones S, Smith P, Beury D, Matthews L, Rogers J, Pask AJ, Shaw G, VandeBerg JL, McCarrey JR, SAVOIR Consortium, Renfree MB, Reik W and Dunham I

    The Babraham Institute, Laboratory of Developmental Genetics and Imprinting, Cambridge CB22 3AT, UK.

    Comparisons between eutherians and marsupials suggest limited conservation of the molecular mechanisms that control genomic imprinting in mammals. We have studied the evolution of the imprinted IGF2-H19 locus in therians. Although marsupial orthologs of protein-coding exons were easily identified, the use of evolutionarily conserved regions and low-stringency Bl2seq comparisons was required to delineate a candidate H19 noncoding RNA sequence. The therian H19 orthologs show miR-675 and exon structure conservation, suggesting functional selection on both features. Transcription start site sequences and poly(A) signals are also conserved. As in eutherians, marsupial H19 is maternally expressed and paternal methylation upstream of the gene originates in the male germline, encompasses a CTCF insulator, and spreads somatically into the H19 gene. The conservation in all therians of the mechanism controlling imprinting of the IGF2-H19 locus suggests a sequential model of imprinting evolution.

    Funded by: Biotechnology and Biological Sciences Research Council: BBS/E/B/0000C214, BBS/E/B/0000S119; Medical Research Council: G0400154; Wellcome Trust

    Nature genetics 2008;40;8;971-6

  • The novel mouse mutation Oblivion inactivates the PMCA2 pump and causes progressive hearing loss.

    Spiden SL, Bortolozzi M, Di Leva F, de Angelis MH, Fuchs H, Lim D, Ortolano S, Ingham NJ, Brini M, Carafoli E, Mammano F and Steel KP

    Wellcome Trust Sanger Institute, Genome Campus, Hinxton, Cambridge, UK.

    Progressive hearing loss is common in the human population, but we have few clues to the molecular basis. Mouse mutants with progressive hearing loss offer valuable insights, and ENU (N-ethyl-N-nitrosourea) mutagenesis is a useful way of generating models. We have characterised a new ENU-induced mouse mutant, Oblivion (allele symbol Obl), showing semi-dominant inheritance of hearing impairment. Obl/+ mutants showed increasing hearing impairment from post-natal day (P)20 to P90, and loss of auditory function was followed by a corresponding base to apex progression of hair cell degeneration. Obl/Obl mutants were small, showed severe vestibular dysfunction by 2 weeks of age, and were completely deaf from birth; sensory hair cells were completely degenerate in the basal turn of the cochlea, although hair cells appeared normal in the apex. We mapped the mutation to Chromosome 6. Mutation analysis of Atp2b2 showed a missense mutation (2630C-->T) in exon 15, causing a serine to phenylalanine substitution (S877F) in transmembrane domain 6 of the PMCA2 pump, the resident Ca(2+) pump of hair cell stereocilia. Transmembrane domain mutations in these pumps generally are believed to be incompatible with normal targeting of the protein to the plasma membrane. However, analyses of hair cells in cultured utricular maculae of Obl/Obl mice and of the mutant Obl pump in model cells showed that the protein was correctly targeted to the plasma membrane. Biochemical and biophysical characterisation showed that the pump had lost a significant portion of its non-stimulated Ca(2+) exporting ability. These findings can explain the progressive loss of auditory function, and indicate the limits in our ability to predict mechanism from sequence alone.

    Funded by: Medical Research Council: G0300212, MC_QA137918; Telethon: GGP04169; Wellcome Trust

    PLoS genetics 2008;4;10;e1000238

  • Large recurrent microdeletions associated with schizophrenia.

    Stefansson H, Rujescu D, Cichon S, Pietiläinen OP, Ingason A, Steinberg S, Fossdal R, Sigurdsson E, Sigmundsson T, Buizer-Voskamp JE, Hansen T, Jakobsen KD, Muglia P, Francks C, Matthews PM, Gylfason A, Halldorsson BV, Gudbjartsson D, Thorgeirsson TE, Sigurdsson A, Jonasdottir A, Jonasdottir A, Bjornsson A, Mattiasdottir S, Blondal T, Haraldsson M, Magnusdottir BB, Giegling I, Möller HJ, Hartmann A, Shianna KV, Ge D, Need AC, Crombie C, Fraser G, Walker N, Lonnqvist J, Suvisaari J, Tuulio-Henriksson A, Paunio T, Toulopoulou T, Bramon E, Di Forti M, Murray R, Ruggeri M, Vassos E, Tosato S, Walshe M, Li T, Vasilescu C, Mühleisen TW, Wang AG, Ullum H, Djurovic S, Melle I, Olesen J, Kiemeney LA, Franke B, GROUP, Sabatti C, Freimer NB, Gulcher JR, Thorsteinsdottir U, Kong A, Andreassen OA, Ophoff RA, Georgi A, Rietschel M, Werge T, Petursson H, Goldstein DB, Nöthen MM, Peltonen L, Collier DA, St Clair D and Stefansson K

    CNS Division, deCODE genetics, Sturlugata 8, IS-101 Reykjavík, Iceland.

    Reduced fecundity, associated with severe mental disorders, places negative selection pressure on risk alleles and may explain, in part, why common variants have not been found that confer risk of disorders such as autism, schizophrenia and mental retardation. Thus, rare variants may account for a larger fraction of the overall genetic risk than previously assumed. In contrast to rare single nucleotide mutations, rare copy number variations (CNVs) can be detected using genome-wide single nucleotide polymorphism arrays. This has led to the identification of CNVs associated with mental retardation and autism. In a genome-wide search for CNVs associating with schizophrenia, we used a population-based sample to identify de novo CNVs by analysing 9,878 transmissions from parents to offspring. The 66 de novo CNVs identified were tested for association in a sample of 1,433 schizophrenia cases and 33,250 controls. Three deletions at 1q21.1, 15q11.2 and 15q13.3 showing nominal association with schizophrenia in the first sample (phase I) were followed up in a second sample of 3,285 cases and 7,951 controls (phase II). All three deletions significantly associate with schizophrenia and related psychoses in the combined sample. The identification of these rare, recurrent risk variants, having occurred independently in multiple founders and being subject to negative selection, is important in itself. CNV analysis may also point the way to the identification of additional and more prevalent risk variants in genes and pathways involved in schizophrenia.

    Funded by: Department of Health: PDA/02/06/016; Medical Research Council: G0901310; NIMH NIH HHS: R01 MH078075, R01MH71425-01A1; Wellcome Trust: 089061

    Nature 2008;455;7210;232-6

  • Insights from the complete genome sequence of Mycobacterium marinum on the evolution of Mycobacterium tuberculosis.

    Stinear TP, Seemann T, Harrison PF, Jenkin GA, Davies JK, Johnson PD, Abdellah Z, Arrowsmith C, Chillingworth T, Churcher C, Clarke K, Cronin A, Davis P, Goodhead I, Holroyd N, Jagels K, Lord A, Moule S, Mungall K, Norbertczak H, Quail MA, Rabbinowitsch E, Walker D, White B, Whitehead S, Small PL, Brosch R, Ramakrishnan L, Fischbach MA, Parkhill J and Cole ST

    Department of Microbiology, Monash University, Clayton 3800, Australia.

    Mycobacterium marinum, a ubiquitous pathogen of fish and amphibia, is a near relative of Mycobacterium tuberculosis, the etiologic agent of tuberculosis in humans. The genome of the M strain of M. marinum comprises a 6,636,827-bp circular chromosome with 5424 CDS, 10 prophages, and a 23-kb mercury-resistance plasmid. Prominent features are the very large number of genes (57) encoding polyketide synthases (PKSs) and nonribosomal peptide synthases (NRPSs) and the most extensive repertoire yet reported of the mycobacteria-restricted PE and PPE proteins, and related-ESX secretion systems. Some of the NRPS genes comprise a novel family and seem to have been acquired horizontally. M. marinum is used widely as a model organism to study M. tuberculosis pathogenesis, and genome comparisons confirmed the close genetic relationship between these two species, as they share 3000 orthologs with an average amino acid identity of 85%. Comparisons with the more distantly related Mycobacterium avium subspecies paratuberculosis and Mycobacterium smegmatis reveal how an ancestral generalist mycobacterium evolved into M. tuberculosis and M. marinum. M. tuberculosis has undergone genome downsizing and extensive lateral gene transfer to become a specialized pathogen of humans and other primates without retaining an environmental niche. M. marinum has maintained a large genome so as to retain the capacity for environmental survival while becoming a broad host range pathogen that produces disease strikingly similar to M. tuberculosis. The work described herein provides a foundation for using M. marinum to better understand the determinants of pathogenesis of tuberculosis.

    Funded by: NIAID NIH HHS: R01 AI036396; Wellcome Trust

    Genome research 2008;18;5;729-41

  • The emerging landscape of breast cancer susceptibility.

    Stratton MR and Rahman N

    Section of Cancer Genetics, Institute of Cancer Research, 15 Cotswold Road, Sutton, Surrey SM2 5NG, UK.

    The genetic basis of inherited predisposition to breast cancer has been assiduously investigated for the past two decades and has been the subject of several recent discoveries. Three reasonably well-defined classes of breast cancer susceptibility alleles with different levels of risk and prevalence in the population have become apparent: rare high-penetrance alleles, rare moderate-penetrance alleles and common low-penetrance alleles. The contribution of each component to breast cancer predisposition is still to be fully explored, as are the phenotypic characteristics of the cancers associated with them, the ways in which they interact, much of their biology and their clinical utility. These recent advances herald a new chapter in the exploration of susceptibility to breast cancer and are likely to provide insights relevant to other common, heterogeneous diseases.

    Nature genetics 2008;40;1;17-22

  • A DNA transposon-based approach to validate oncogenic mutations in the mouse.

    Su Q, Prosser HM, Campos LS, Ortiz M, Nakamura T, Warren M, Dupuy AJ, Jenkins NA, Copeland NG, Bradley A and Liu P

    Wellcome Trust Sanger Institute, Hinxton, Cambridge CB10 1HH, England.

    Large-scale cancer genome projects will soon be able to sequence many cancer genomes to comprehensively identify genetic changes in human cancer. Genome-wide association studies have also identified putative cancer associated loci. Functional validation of these genetic mutations in vivo is becoming a challenge. We describe here a DNA transposon-based platform that permits us to explore the oncogenic potential of genetic mutations in the mouse. Briefly, promoter-less human cancer gene cDNAs were first cloned into Sleeping Beauty (SB) transposons. DNA transposition in the mouse that carried both the transposons and the SB transposase made it possible for the cDNAs to be expressed from an appropriate endogenous genomic locus and in the relevant cell types for tumor development. Consequently, these mice developed a broad spectrum of tumors at very early postnatal stages. This technology thus complements the large-scale cancer genome projects.

    Funded by: Wellcome Trust

    Proceedings of the National Academy of Sciences of the United States of America 2008;105;50;19904-9

  • Promoting coherent minimum reporting guidelines for biological and biomedical investigations: the MIBBI project.

    Taylor CF, Field D, Sansone SA, Aerts J, Apweiler R, Ashburner M, Ball CA, Binz PA, Bogue M, Booth T, Brazma A, Brinkman RR, Michael Clark A, Deutsch EW, Fiehn O, Fostel J, Ghazal P, Gibson F, Gray T, Grimes G, Hancock JM, Hardy NW, Hermjakob H, Julian RK, Kane M, Kettner C, Kinsinger C, Kolker E, Kuiper M, Le Novère N, Leebens-Mack J, Lewis SE, Lord P, Mallon AM, Marthandan N, Masuya H, McNally R, Mehrle A, Morrison N, Orchard S, Quackenbush J, Reecy JM, Robertson DG, Rocca-Serra P, Rodriguez H, Rosenfelder H, Santoyo-Lopez J, Scheuermann RH, Schober D, Smith B, Snape J, Stoeckert CJ, Tipton K, Sterk P, Untergasser A, Vandesompele J and Wiemann S

    European Bioinformatics Institute, Wellcome Trust Genome Campus, Hinxton, Cambridgeshire CB10 1SD, UK.,

    Funded by: Biotechnology and Biological Sciences Research Council: BB/E006248/1, BB/E025080/1, E025080/1; Medical Research Council: MC_U142684171; NHGRI NIH HHS: U54 HG004028; NIBIB NIH HHS: EB005034-01, R01 EB005034, R01 EB005034-04

    Nature biotechnology 2008;26;8;889-96

  • Whole genome-amplified DNA: insights and imputation.

    Teo YY, Inouye M, Small KS, Fry AE, Potter SC, Dunstan SJ, Seielstad M, Barroso I, Wareham NJ, Rockett KA, Kwiatkowski DP and Deloukas P

    Funded by: Medical Research Council: G0600230, G19/9, MC_U106179471; Wellcome Trust: 077011, 077016

    Nature methods 2008;5;4;279-80

  • Maternal footprints of Southeast Asians in North India.

    Thangaraj K, Chaubey G, Kivisild T, Selvi Rani D, Singh VK, Ismail T, Carvalho-Silva D, Metspalu M, Bhaskar LV, Reddy AG, Chandra S, Pande V, Prathap Naidu B, Adarsh N, Verma A, Jyothi IA, Mallick CB, Shrivastava N, Devasena R, Kumari B, Singh AK, Dwivedi SK, Singh S, Rao G, Gupta P, Sonvane V, Kumari K, Basha A, Bhargavi KR, Lalremruata A, Gupta AK, Kaur G, Reddy KK, Rao AP, Villems R, Tyler-Smith C and Singh L

    Centre for Cellular and Molecular Biology, Hyderabad, India.

    We have analyzed 7,137 samples from 125 different caste, tribal and religious groups of India and 99 samples from three populations of Nepal for the length variation in the COII/tRNA(Lys) region of mtDNA. Samples showing length variation were subjected to detailed phylogenetic analysis based on HVS-I and informative coding region sequence variation. The overall frequencies of the 9-bp deletion and insertion variants in South Asia were 1.9 and 0.6%, respectively. We have also defined a novel deep-rooting haplogroup M43 and identified the rare haplogroup H14 in Indian populations carrying the 9-bp deletion by complete mtDNA sequencing. Moreover, we redefined haplogroup M6 and dissected it into two well-defined subclades. The presence of haplogroups F1 and B5a in Uttar Pradesh suggests minor maternal contribution from Southeast Asia to Northern India. The occurrence of haplogroup F1 in the Nepalese sample implies that Nepal might have served as a bridge for the flow of eastern lineages to India. The presence of R6 in the Nepalese, on the other hand, suggests that the gene flow between India and Nepal has been reciprocal.

    Funded by: Wellcome Trust: 077009

    Human heredity 2008;66;1;1-9

  • Characterization of 6q deletions in mature B cell lymphomas and childhood acute lymphoblastic leukemia.

    Thelander EF, Ichimura K, Corcoran M, Barbany G, Nordgren A, Heyman M, Berglund M, Mungall A, Rosenquist R, Collins VP, Grandér D, Larsson C and Lagercrantz S

    Medical Genetics Unit, Department of Molecular Medicine and Surgery, Karolinska Institutet, Stockholm, Sweden.

    The study was undertaken with the aim to outline deletion patterns involving the long arm of chromosome 6, a common abnormality in lymphoproliferative disorders. Using a chromosome 6 specific tile path array, 60 samples from in total 49 cases with mantle cell lymphoma (MCL), de novo diffuse large B-cell lymphoma (DLBCL), transformed DLBCL as well as preceding follicular lymphoma (FL), and childhood acute lymphoblastic leukemia (ALL), were characterized. Twenty-six of the studied cases, representing all diagnoses, showed a 6q deletion among which 85% involved a 3 Mb region in 6q21. The minimal deleted interval in 6q21 encompasses the FOXO3A, PRDM1 and HACE1 candidate genes. The PRDM1 gene was found homozygously deleted in a case of DLBCL. Moreover, in two DLBCL cases, an overlapping homozygous deletion was identified in 6q23.3 - 24.1, encompassing the TNFAIP3 gene among others. Taken together, we refined the deletion pattern within the long arm of chromosome 6 in four different types of hematological malignances, suggesting the location of tumor suppressor genes involved in the tumor progression.

    Leukemia & lymphoma 2008;49;3;477-87

  • Phenotypical characteristics of idiopathic infantile nystagmus with and without mutations in FRMD7.

    Thomas S, Proudlock FA, Sarvananthan N, Roberts EO, Awan M, McLean R, Surendran M, Kumar AS, Farooq SJ, Degg C, Gale RP, Reinecke RD, Woodruff G, Langmann A, Lindner S, Jain S, Tarpey P, Raymond FL and Gottlob I

    Ophthalmology Group, University of Leicester, Leicester, UK.

    Idiopathic infantile nystagmus (IIN) consists of involuntary oscillations of the eyes. The familial form is most commonly X-linked. We recently found mutations in a novel gene FRMD7 (Xq26.2), which provided an opportunity to investigate a genetically defined and homogeneous group of patients with nystagmus. We compared clinical features and eye movement recordings of 90 subjects with mutation in the gene (FRMD7 group) to 48 subjects without mutations but with clinical IIN (non-FRMD7 group). Fifty-eight female obligate carriers of the mutation were also investigated. The median visual acuity (VA) was 0.2 logMAR (Snellen equivalent 6/9) in both groups and most patients had good stereopsis. The prevalence of strabismus was also similar (FRMD7: 7.8%, non-FRMD7: 10%). The presence of anomalous head posture (AHP) was significantly higher in the non-FRMD7 group (P < 0.0001). The amplitude of nystagmus was more strongly dependent on the direction of gaze in the FRMD7 group being lower at primary position (P < 0.0001), compared to non-FRMD7 group (P = 0.83). Pendular nystagmus waveforms were also more frequent in the FRMD7 group (P = 0.003). Fifty-three percent of the obligate female carriers of an FRMD7 mutation were clinically affected. The VA's in affected females were slightly better compared to affected males (P = 0.014). Subnormal optokinetic responses were found in a subgroup of obligate unaffected carriers, which may be interpreted as a sub-clinical manifestation. FRMD7 is a major cause of X-linked IIN. Most clinical and eye movement characteristics were similar in the FRMD7 group and non-FRMD7 group with most patients having good VA and stereopsis and low incidence of strabismus. Fewer patients in the FRMD7 group had AHPs, their amplitude of nystagmus being lower in primary position. Our findings are helpful in the clinical identification of IIN and genetic counselling of nystagmus patients.

    Brain : a journal of neurology 2008;131;Pt 5;1259-67

  • Vive la différence.

    Thomson NR

    Nature reviews. Microbiology 2008;6;7;502-3

  • Comparative genome analysis of Salmonella Enteritidis PT4 and Salmonella Gallinarum 287/91 provides insights into evolutionary and host adaptation pathways.

    Thomson NR, Clayton DJ, Windhorst D, Vernikos G, Davidson S, Churcher C, Quail MA, Stevens M, Jones MA, Watson M, Barron A, Layton A, Pickard D, Kingsley RA, Bignell A, Clark L, Harris B, Ormond D, Abdellah Z, Brooks K, Cherevach I, Chillingworth T, Woodward J, Norberczak H, Lord A, Arrowsmith C, Jagels K, Moule S, Mungall K, Sanders M, Whitehead S, Chabalgoity JA, Maskell D, Humphrey T, Roberts M, Barrow PA, Dougan G and Parkhill J

    The Pathogen Sequencing Unit, The Wellcome Trust Sanger Institute, Hinxton, Cambridge CB10 1SA, United Kingdom.

    We have determined the complete genome sequences of a host-promiscuous Salmonella enterica serovar Enteritidis PT4 isolate P125109 and a chicken-restricted Salmonella enterica serovar Gallinarum isolate 287/91. Genome comparisons between these and other Salmonella isolates indicate that S. Gallinarum 287/91 is a recently evolved descendent of S. Enteritidis. Significantly, the genome of S. Gallinarum has undergone extensive degradation through deletion and pseudogene formation. Comparison of the pseudogenes in S. Gallinarum with those identified previously in other host-adapted bacteria reveals the loss of many common functional traits and provides insights into possible mechanisms of host and tissue adaptation. We propose that experimental analysis in chickens and mice of S. Enteritidis-harboring mutations in functional homologs of the pseudogenes present in S. Gallinarum could provide an experimentally tractable route toward unraveling the genetic basis of host adaptation in S. enterica.

    Funded by: Wellcome Trust

    Genome research 2008;18;10;1624-37

  • Evolutionary plasticity of genetic interaction networks.

    Tischler J, Lehner B and Fraser AG

    The Wellcome Trust Sanger Institute, Hinxton, Cambridge CB10 1HH, UK.

    Non-additive genetic interactions contribute to many genetic disorders, but they are extremely difficult to predict. Here we show that genetic interactions identified in yeast, unlike gene functions or protein interactions, are not highly conserved in animals. Genetic interactions are therefore unlikely to represent simple redundancy between genes or pathways, and genetic interactions from yeast do not directly predict genetic interactions in higher eukaryotes, including humans.

    Funded by: Wellcome Trust

    Nature genetics 2008;40;4;390-1

  • Long-range, high-throughput haplotype determination via haplotype-fusion PCR and ligation haplotyping.

    Turner DJ, Tyler-Smith C and Hurles ME

    Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Cambridge, UK.

    Ligation Haplotyping is a robust, novel method for experimental determination of haplotypes over long distances, which can be applied to assaying both sequence and structural variation. The simplicity and efficacy of the method for genotyping large chromosomal rearrangements and haplotyping SNPs over long distances make it a valuable and powerful addition to the methodological repertoire, which will be beneficial to studies of population genetics and evolution, disease association and inheritance, and genomic variation. We illustrate the versatility of the method both by genotyping a Yp paracentric inversion, found in approximately 60% of Northwest European males, that strongly influences the germline rate of infertility-causing XY translocations and by haplotyping two autosomal SNPs that lie 16.4 kb apart on chromosome 7, and which influence an individual's susceptibility to systemic lupus erythematosus.

    Funded by: Wellcome Trust

    Nucleic acids research 2008;36;13;e82

  • An evolutionary perspective on Y-chromosomal variation and male infertility.

    Tyler-Smith C

    The Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Cambridgeshire, UK.

    Genetic variation on the Y chromosome is one of the best-documented causes of male infertility, but the genes responsible have still not been identified. This review discusses how an evolutionary perspective may help with interpretation of the data available and suggest novel approaches to identify key genes. Comparison with the chimpanzee Y chromosome indicates that USP9Y is dispensable in apes, but that multiple copies of TSPY1 may have an important role. Comparisons between infertile and control groups in search of genetic susceptibility factors are more complex for the Y chromosome than for the rest of the genome because of population stratification and require unusual levels of confirmation. But the extreme population stratification exhibited by the Y also allows populations particularly suitable for some studies to be identified, such as the partial AZFc deletions common in Northern European populations where further dissection of this complex structural region would be facilitated.

    Funded by: Wellcome Trust

    International journal of andrology 2008;31;4;376-82

  • Loss of Rassf1a cooperates with Apc(Min) to accelerate intestinal tumourigenesis.

    van der Weyden L, Arends MJ, Dovey OM, Harrison HL, Lefebvre G, Conte N, Gergely FV, Bradley A and Adams DJ

    Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Cambridge, UK.

    Promoter methylation of the RAS-association domain family 1, isoform A gene (RASSF1A) is one of the most frequent events found in human tumours. In this study we set out to test the hypothesis that loss of Rassf1a can cooperate with inactivation of the adenomatous polyposis coli (Apc) gene to accelerate intestinal tumourigenesis using the Apc-Min (Apc(Min/+)) mouse model, as mutational or deletional inactivation of APC is a frequent early event in the genesis of intestinal cancer. Further, loss of RASSF1A has also been reported to occur in premalignant adenomas of the bowel. RASSF1A has been implicated in an array of pivotal cellular processes, including regulation of the cell cycle, apoptosis, microtubule stability and most recently in the beta-catenin signalling pathway. By interbreeding isoform specific Rassf1a knockout mice with Apc(+/Min) mice, we showed that loss of Rassf1a results in a significant increase in adenomas of the small intestine and accelerated intestinal tumourigenesis leading to the earlier death of adenocarcinoma-bearing mice and decreased overall survival. Comparative genomic hybridization of adenomas from Rassf1a(-/-); Apc(+/Min) mice revealed no evidence of aneuploidy or gross chromosomal instability (no difference to adenomas from Rassf1a(+/+); Apc(+/Min) mice). Immunohistochemical analysis of adenomas revealed increased nuclear beta-catenin accumulation in adenomas from Rassf1a(-/-); Apc(+/Min) mice, compared to those from Rassf1a(+/+); Apc(+/Min) mice, but no differences in proliferation marker (Ki67) staining patterns. Collectively these data demonstrate cooperation between inactivation of Rassf1a and Apc resulting in accelerated intestinal tumourigenesis, with adenomas showing increased nuclear accumulation of beta-catenin, supporting a mechanistic link via loss of the known interaction of Rassf1 with beta-TrCP that usually mediates degradation of beta-catenin.

    Funded by: Cancer Research UK: A8449; Wellcome Trust: 079643

    Oncogene 2008;27;32;4503-8

  • High-resolution mapping of expression-QTLs yields insight into human gene regulation.

    Veyrieras JB, Kudaravalli S, Kim SY, Dermitzakis ET, Gilad Y, Stephens M and Pritchard JK

    Department of Human Genetics, The University of Chicago, Chicago, IL, USA.

    Recent studies of the HapMap lymphoblastoid cell lines have identified large numbers of quantitative trait loci for gene expression (eQTLs). Reanalyzing these data using a novel Bayesian hierarchical model, we were able to create a surprisingly high-resolution map of the typical locations of sites that affect mRNA levels in cis. Strikingly, we found a strong enrichment of eQTLs in the 250 bp just upstream of the transcription end site (TES), in addition to an enrichment around the transcription start site (TSS). Most eQTLs lie either within genes or close to genes; for example, we estimate that only 5% of eQTLs lie more than 20 kb upstream of the TSS. After controlling for position effects, SNPs in exons are approximately 2-fold more likely than SNPs in introns to be eQTLs. Our results suggest an important role for mRNA stability in determining steady-state mRNA levels, and highlight the potential of eQTL mapping as a high-resolution tool for studying the determinants of gene regulation.

    Funded by: Howard Hughes Medical Institute; NHGRI NIH HHS: HG002772, HG02585-01, R01 HG002585, R01 HG002772; NIGMS NIH HHS: GM077959, R01 GM077959

    PLoS genetics 2008;4;10;e1000214

  • Habitual energy expenditure modifies the association between NOS3 gene polymorphisms and blood pressure.

    Vimaleswaran KS, Franks PW, Barroso I, Brage S, Ekelund U, Wareham NJ and Loos RJ

    Medical Research Council Epidemiology Unit, Institute of Metabolic Science, Cambridge, UK.

    Background: The endothelial nitric-oxide synthase (NOS3) gene encodes the enzyme (eNOS) that synthesizes the molecule nitric oxide, which facilitates endothelium-dependent vasodilation in response to physical activity. Thus, energy expenditure may modify the association between the genetic variation at NOS3 and blood pressure.

    Methods: To test this hypothesis, we genotyped 11 NOS3 polymorphisms, capturing all common variations, in 726 men and women from the Medical Research Council (MRC) Ely Study (age (mean +/- s.d.): 55 +/- 10 years, body mass index: 26.4 +/- 4.1 kg/m(2)). Habitual/non-resting energy expenditure (NREE) was assessed via individually calibrated heart rate monitoring over 4 days.

    Results: The intronic variant, IVS25+15 [G-->A], was significantly associated with blood pressure; GG homozygotes had significantly lower levels of diastolic blood pressure (DBP) (-2.8 mm Hg; P = 0.016) and systolic blood pressure (SBP) (-1.9 mm Hg; P = 0.018) than A-allele carriers. The interaction between NREE and IVS25+15 was also significant for both DBP (P = 0.006) and SBP (P = 0.026), in such a way that the effect of the GG-genotype on blood pressure was stronger in individuals with higher NREE (DBP: -4.9 mm Hg, P = 0.02. SBP: -3.8 mm Hg, P= 0.03 for the third tertile). Similar results were observed when the outcome was dichotomously defined as hypertension.

    Conclusions: In summary, the NOS3 IVS25+15 is directly associated with blood pressure and hypertension in white Europeans. However, the associations are most evident in the individuals with the highest NREE. These results need further replication and have to be ideally tested in a trial before being informative for targeted disease prevention. Eventually, the selection of individuals for lifestyle intervention programs could be guided by knowledge of genotype.

    Funded by: Medical Research Council: MC_U106179471, MC_U106179473, MC_U106188470, U.1061.00.001 (79471), U.1061.00.005(79473); Wellcome Trust: 077016, 087636

    American journal of hypertension 2008;21;3;297-302

  • The Gly482Ser genotype at the PPARGC1A gene and elevated blood pressure: a meta-analysis involving 13,949 individuals.

    Vimaleswaran KS, Luan J, Andersen G, Muller YL, Wheeler E, Brito EC, O'Rahilly S, Pedersen O, Baier LJ, Knowler WC, Barroso I, Wareham NJ, Loos RJ and Franks PW

    Department of Public Health & Clinical Medicine, Umeå University Hospital, Umeå, Sweden.

    The protein encoded by the PPARGC1A gene is expressed at high levels in metabolically active tissues and is involved in the control of oxidative stress via reactive oxygen species detoxification. Several recent reports suggest that the PPARGC1A Gly482Ser (rs8192678) missense polymorphism may relate inversely with blood pressure. We used conventional meta-analysis methods to assess the association between Gly482Ser and systolic (SBP) or diastolic blood pressures (DBP) or hypertension in 13,949 individuals from 17 studies, of which 6,042 were previously unpublished observations. The studies comprised cohorts of white European, Asian, and American Indian adults, and adolescents from South America. Stratified analyses were conducted to control for population stratification. Pooled genotype frequencies were 0.47 (Gly482Gly), 0.42 (Gly482Ser), and 0.11 (Ser482Ser). We found no evidence of association between Gly482Ser and SBP [Gly482Gly: mean = 131.0 mmHg, 95% confidence interval (CI) = 130.5-131.5 mmHg; Gly482Ser mean = 133.1 mmHg, 95% CI = 132.6-133.6 mmHg; Ser482Ser: mean = 133.5 mmHg, 95% CI = 132.5-134.5 mmHg; P = 0.409] or DBP (Gly482Gly: mean = 80.3 mmHg, 95% CI = 80.0-80.6 mmHg; Gly482Ser mean = 81.5 mmHg, 95% CI = 81.2-81.8 mmHg; Ser482Ser: mean = 82.1 mmHg, 95% CI = 81.5-82.7 mmHg; P = 0.651). Contrary to previous reports, we did not observe significant effect modification by sex (SBP, P = 0.966; DBP, P = 0.715). We were also unable to confirm the previously reported association between the Ser482 allele and hypertension [odds ratio: 0.97, 95% CI = 0.87-1.08, P = 0.585]. These results were materially unchanged when analyses were focused on whites only. However, statistical evidence of gene-age interaction was apparent for DBP [Gly482Gly: 73.5 (72.8, 74.2), Gly482Ser: 77.0 (76.2, 77.8), Ser482Ser: 79.1 (77.4, 80.9), P = 4.20 x 10(-12)] and SBP [Gly482Gly: 121.4 (120.4, 122.5), Gly482Ser: 125.9 (124.6, 127.1), Ser482Ser: 129.2 (126.5, 131.9), P = 7.20 x 10(-12)] in individuals <50 yr (n = 2,511); these genetic effects were absent in those older than 50 yr (n = 5,088) (SBP, P = 0.41; DBP, P = 0.51). Our findings suggest that the PPARGC1A Ser482 allele may be associated with higher blood pressure, but this is only apparent in younger adults.

    Funded by: Intramural NIH HHS; Medical Research Council: MC_U106179471, MC_U106188470; Wellcome Trust: 077016

    Journal of applied physiology (Bethesda, Md. : 1985) 2008;105;4;1352-8

  • Single-cell genomics.

    Walker A and Parkhill J

    Nature reviews. Microbiology 2008;6;3;176-7

  • Chromosomal transposition of PiggyBac in mouse embryonic stem cells.

    Wang W, Lin C, Lu D, Ning Z, Cox T, Melvin D, Wang X, Bradley A and Liu P

    The Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridgeshire CB10 1SA, United Kingdom.

    Transposon systems are widely used for generating mutations in various model organisms. PiggyBac (PB) has recently been shown to transpose efficiently in the mouse germ line and other mammalian cell lines. To facilitate PB's application in mammalian genetics, we characterized the properties of the PB transposon in mouse embryonic stem (ES) cells. We first measured the transposition efficiencies of PB transposon in mouse embryonic stem cells. We next constructed a PB/SB hybrid transposon to compare PB and Sleeping Beauty (SB) transposon systems and demonstrated that PB transposition was inhibited by DNA methylation. The excision and reintegration rates of a single PB from two independent genomic loci were measured and its ability to mutate genes with gene trap cassettes was tested. We examined PB's integration site distribution in the mouse genome and found that PB transposition exhibited local hopping. The comprehensive information from this study should facilitate further exploration of the potential of PB and SB DNA transposons in mammalian genetics.

    Funded by: NIGMS NIH HHS: 5R21GM079528, R21 GM079528; Wellcome Trust

    Proceedings of the National Academy of Sciences of the United States of America 2008;105;27;9290-5

  • Genome analysis of the platypus reveals unique signatures of evolution.

    Warren WC, Hillier LW, Marshall Graves JA, Birney E, Ponting CP, Grützner F, Belov K, Miller W, Clarke L, Chinwalla AT, Yang SP, Heger A, Locke DP, Miethke P, Waters PD, Veyrunes F, Fulton L, Fulton B, Graves T, Wallis J, Puente XS, López-Otín C, Ordóñez GR, Eichler EE, Chen L, Cheng Z, Deakin JE, Alsop A, Thompson K, Kirby P, Papenfuss AT, Wakefield MJ, Olender T, Lancet D, Huttley GA, Smit AF, Pask A, Temple-Smith P, Batzer MA, Walker JA, Konkel MK, Harris RS, Whittington CM, Wong ES, Gemmell NJ, Buschiazzo E, Vargas Jentzsch IM, Merkel A, Schmitz J, Zemann A, Churakov G, Kriegs JO, Brosius J, Murchison EP, Sachidanandam R, Smith C, Hannon GJ, Tsend-Ayush E, McMillan D, Attenborough R, Rens W, Ferguson-Smith M, Lefèvre CM, Sharp JA, Nicholas KR, Ray DA, Kube M, Reinhardt R, Pringle TH, Taylor J, Jones RC, Nixon B, Dacheux JL, Niwa H, Sekita Y, Huang X, Stark A, Kheradpour P, Kellis M, Flicek P, Chen Y, Webber C, Hardison R, Nelson J, Hallsworth-Pepin K, Delehaunty K, Markovic C, Minx P, Feng Y, Kremitzki C, Mitreva M, Glasscock J, Wylie T, Wohldmann P, Thiru P, Nhan MN, Pohl CS, Smith SM, Hou S, Nefedov M, de Jong PJ, Renfree MB, Mardis ER and Wilson RK

    Genome Sequencing Center, Washington University School of Medicine, Campus Box 8501, 4444 Forest Park Avenue, St Louis, Missouri 63108, USA.

    We present a draft genome sequence of the platypus, Ornithorhynchus anatinus. This monotreme exhibits a fascinating combination of reptilian and mammalian characters. For example, platypuses have a coat of fur adapted to an aquatic lifestyle; platypus females lactate, yet lay eggs; and males are equipped with venom similar to that of reptiles. Analysis of the first monotreme genome aligned these features with genetic innovations. We find that reptile and platypus venom proteins have been co-opted independently from the same gene families; milk protein genes are conserved despite platypuses laying eggs; and immune gene family expansions are directly related to platypus biology. Expansions of protein, non-protein-coding RNA and microRNA families, as well as repeat elements, are identified. Sequencing of this genome now provides a valuable resource for deep mammalian comparative analyses, as well as for monotreme biology and conservation.

    Funded by: Medical Research Council: MC_U137761446; NCI NIH HHS: P01 CA013106, P01 CA013106-37; NHGRI NIH HHS: HG002238, R01 HG002238, R01 HG002385, R01 HG002939, R01 HG004037, R01 HG004037-02, R01HG02385; NIGMS NIH HHS: R01 GM059290, R01 GM59290; Wellcome Trust: 062023

    Nature 2008;453;7192;175-83

  • urg1: a uracil-regulatable promoter system for fission yeast with short induction and repression times.

    Watt S, Mata J, López-Maury L, Marguerat S, Burns G and Bähler J

    Cancer Research United Kingdom Fission Yeast Functional Genomics Group, Wellcome Trust Sanger Institute, Hinxton, Cambridge, United Kingdom.

    Background: The fission yeast Schizosaccharomyces pombe is a popular genetic model organism with powerful experimental tools. The thiamine-regulatable nmt1 promoter and derivatives, which take >15 hours for full induction, are most commonly used for controlled expression of ectopic genes. Given the short cell cycle of fission yeast, however, a promoter system that can be rapidly regulated, similar to the GAL system for budding yeast, would provide a key advantage for many experiments.

    Methodology/principal findings: We used S. pombe microarrays to identify three neighbouring genes (urg1, urg2, and urg3) whose transcript levels rapidly and strongly increased in response to uracil, a condition which otherwise had little effect on global gene expression. We cloned the promoter of urg1 (uracil-regulatable gene) to create several PCR-based gene targeting modules for replacing native promoters with the urg1 promoter (Purg1) in the normal chromosomal locations of genes of interest. The kanMX6 and natMX6 markers allow selection under urg1 induced and repressed conditions, respectively. Some modules also allow N-terminal tagging of gene products placed under urg1 control. Using pom1 as a proof-of-principle, we observed a maximal increase of Purg1-pom1 transcripts after uracil addition within less than 30 minutes, and a similarly rapid decrease after uracil removal. The induced and repressed transcriptional states remained stable over 24-hour periods. RT-PCR comparisons showed that both induced and repressed Purg1-pom1 transcript levels were lower than corresponding P3nmt1-pom1 levels (wild-type nmt1 promoter) but higher than P81nmt1-pom1 levels (weak nmt1 derivative).

    Conclusions/significance: We exploited the urg1 promoter system to rapidly induce pom1 expression at defined cell-cycle stages, showing that ectopic pom1 expression leads to cell branching in G2-phase but much less so in G1-phase. The high temporal resolution provided by the urg1 promoter should facilitate experimental design and improve the genetic toolbox for the fission yeast community.

    Funded by: Cancer Research UK: A6517, C9546/A6517; Wellcome Trust: 077118

    PloS one 2008;3;1;e1428

  • Genome-wide association analysis identifies 20 loci that influence adult height.

    Weedon MN, Lango H, Lindgren CM, Wallace C, Evans DM, Mangino M, Freathy RM, Perry JR, Stevens S, Hall AS, Samani NJ, Shields B, Prokopenko I, Farrall M, Dominiczak A, Diabetes Genetics Initiative, Wellcome Trust Case Control Consortium, Johnson T, Bergmann S, Beckmann JS, Vollenweider P, Waterworth DM, Mooser V, Palmer CN, Morris AD, Ouwehand WH, Cambridge GEM Consortium, Zhao JH, Li S, Loos RJ, Barroso I, Deloukas P, Sandhu MS, Wheeler E, Soranzo N, Inouye M, Wareham NJ, Caulfield M, Munroe PB, Hattersley AT, McCarthy MI and Frayling TM

    Genetics of Complex Traits, Institute of Biomedical and Clinical Science, Peninsula Medical School, Magdalen Road, Exeter EX1 2LU, UK.

    Adult height is a model polygenic trait, but there has been limited success in identifying the genes underlying its normal variation. To identify genetic variants influencing adult human height, we used genome-wide association data from 13,665 individuals and genotyped 39 variants in an additional 16,482 samples. We identified 20 variants associated with adult height (P < 5 x 10(-7), with 10 reaching P < 1 x 10(-10)). Combined, the 20 SNPs explain approximately 3% of height variation, with a approximately 5 cm difference between the 6.2% of people with 17 or fewer 'tall' alleles compared to the 5.5% with 27 or more 'tall' alleles. The loci we identified implicate genes in Hedgehog signaling (IHH, HHIP, PTCH1), extracellular matrix (EFEMP1, ADAMTSL3, ACAN) and cancer (CDK6, HMGA2, DLEU7) pathways, and provide new insights into human growth and developmental processes. Finally, our results provide insights into the genetic architecture of a classic quantitative trait.

    Funded by: British Heart Foundation: FS/05/061/19501, PG/02/128/14470, PG02/128; Medical Research Council: G0600705, G0701863, G9521010, G9521010(63660), G9521010D, MC_U106188470; Wellcome Trust: 076113, 077011, 077016, 090532

    Nature genetics 2008;40;5;575-83

  • Dynamic repertoire of a eukaryotic transcriptome surveyed at single-nucleotide resolution.

    Wilhelm BT, Marguerat S, Watt S, Schubert F, Wood V, Goodhead I, Penkett CJ, Rogers J and Bähler J

    Cancer Research UK Fission Yeast Functional Genomics Group, Wellcome Trust Sanger Institute, Hinxton, Cambridge CB10 1HH, UK.

    Recent data from several organisms indicate that the transcribed portions of genomes are larger and more complex than expected, and that many functional properties of transcripts are based not on coding sequences but on regulatory sequences in untranslated regions or non-coding RNAs. Alternative start and polyadenylation sites and regulation of intron splicing add additional dimensions to the rich transcriptional output. This transcriptional complexity has been sampled mainly using hybridization-based methods under one or few experimental conditions. Here we applied direct high-throughput sequencing of complementary DNAs (RNA-Seq), supplemented with data from high-density tiling arrays, to globally sample transcripts of the fission yeast Schizosaccharomyces pombe, independently from available gene annotations. We interrogated transcriptomes under multiple conditions, including rapid proliferation, meiotic differentiation and environmental stress, as well as in RNA processing mutants to reveal the dynamic plasticity of the transcriptional landscape as a function of environmental, developmental and genetic factors. High-throughput sequencing proved to be a powerful and quantitative method to sample transcriptomes deeply at maximal resolution. In contrast to hybridization, sequencing showed little, if any, background noise and was sensitive enough to detect widespread transcription in >90% of the genome, including traces of RNAs that were not robustly transcribed or rapidly degraded. The combined sequencing and strand-specific array data provide rich condition-specific information on novel, mostly non-coding transcripts, untranslated regions and gene structures, thus improving the existing genome annotation. Sequence reads spanning exon-exon or exon-intron junctions give unique insight into a surprising variability in splicing efficiency across introns, genes and conditions. Splicing efficiency was largely coordinated with transcript levels, and increased transcription led to increased splicing in test genes. Hundreds of introns showed such regulated splicing during cellular proliferation or differentiation.

    Funded by: Cancer Research UK: A6517, C9546/A6517; Wellcome Trust: 077118

    Nature 2008;453;7199;1239-43

  • Sequence differentiation in regions identified by a genome scan for local adaptation.

    Wood HM, Grahame JW, Humphray S, Rogers J and Butlin RK

    Institute of Integrative and Comparative Biology, University of Leeds, Leeds LS29JT, UK.

    Genome scans using large numbers of randomly selected markers have revealed a small proportion of loci that deviate from neutral expectations and so may mark genomic regions that contribute to local adaptation. Measurements of sequence differentiation and identification of genes in these regions is important but difficult, especially in organisms with limited genetic information available. We have followed up a genome scan in the marine gastropod, Littorina saxatilis, by searching a bacterial artificial chromosome library with differentiated and undifferentiated markers, sequencing four bacterial artificial chromosomes and then analysing sequence variation in population samples for fragments at, and close to the original marker polymorphisms. We show that sequence differentiation follows the patterns expected from the original marker frequencies, that differentiated markers identify independent and highly localized sites and that these sites fall outside coding regions. Two differentiated loci are characterized by insertions of putative transposable elements that appear to have increased in frequency recently and which might influence expression of downstream genes. These results provide strong candidate loci for the study of local adaptation in Littorina. They demonstrate an approach that can be applied to follow up genome scans in other taxa and they show that the genome scan approach can lead rapidly to candidate genes in nonmodel organisms.

    Funded by: Biotechnology and Biological Sciences Research Council: G18139

    Molecular ecology 2008;17;13;3123-35

  • Variation of the oxytocin/neurophysin I (OXT) gene in four human populations.

    Xu Y, Xue Y, Asan, Daly A, Wu L and Tyler-Smith C

    The Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambs., CB10 1SA, UK.

    Oxytocin is a short peptide with multiple functions in human biology and has been implicated in autism. We aimed to determine the normal pattern of variation around the oxytocin gene and resequenced it and its flanking regions in 91 individuals from four HapMap populations and one chimpanzee. We identified 14 single nucleotide polymorphisms (SNPs), all noncoding, including eight that were novel. Population genetic analyses were largely consistent with a neutral evolutionary history, but an Hudson-Kreitman-Aguadé (HKA) test revealed more variation within the human population than expected from the level of chimpanzee-human divergence.

    Funded by: Wellcome Trust: 077009

    Journal of human genetics 2008;53;7;637-643

  • Adaptive evolution of UGT2B17 copy-number variation.

    Xue Y, Sun D, Daly A, Yang F, Zhou X, Zhao M, Huang N, Zerjal T, Lee C, Carter NP, Hurles ME and Tyler-Smith C

    The Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton CB10 1SA, UK.

    The human UGT2B17 gene varies in copy number from zero to two per individual and also differs in mean number between populations from Africa, Europe, and East Asia. We show that such a high degree of geographical variation is unusual and investigate its evolutionary history. This required first reinterpreting the reference sequence in this region of the genome, which is misassembled from the two different alleles separated by an artifactual gap. A corrected assembly identifies the polymorphism as a 117 kb deletion arising by nonallelic homologous recombination between approximately 4.9 kb segmental duplications and allows the deletion breakpoint to be identified. We resequenced approximately 12 kb of DNA spanning the breakpoint in 91 humans from three HapMap and one extended HapMap populations and one chimpanzee. Diversity was unusually high and the time to the most recent common ancestor was estimated at approximately 2.4 or approximately 3.0 million years by two different methods, with evidence of balancing selection in Europe. In contrast, diversity was low in East Asia where a single haplotype predominated, suggesting positive selection for the deletion in this part of the world.

    Funded by: Wellcome Trust

    American journal of human genetics 2008;83;3;337-46

  • Y-chromosomal diversity in Lebanon is structured by recent historical events.

    Zalloua PA, Xue Y, Khalife J, Makhoul N, Debiane L, Platt DE, Royyuru AK, Herrera RJ, Hernanz DF, Blue-Smith J, Wells RS, Comas D, Bertranpetit J, Tyler-Smith C and Genographic Consortium

    The Lebanese American University, Chouran, Beirut 1102 2801, Lebanon.

    Lebanon is an eastern Mediterranean country inhabited by approximately four million people with a wide variety of ethnicities and religions, including Muslim, Christian, and Druze. In the present study, 926 Lebanese men were typed with Y-chromosomal SNP and STR markers, and unusually, male genetic variation within Lebanon was found to be more strongly structured by religious affiliation than by geography. We therefore tested the hypothesis that migrations within historical times could have contributed to this situation. Y-haplogroup J*(xJ2) was more frequent in the putative Muslim source region (the Arabian Peninsula) than in Lebanon, and it was also more frequent in Lebanese Muslims than in Lebanese non-Muslims. Conversely, haplogroup R1b was more frequent in the putative Christian source region (western Europe) than in Lebanon and was also more frequent in Lebanese Christians than in Lebanese non-Christians. The most common R1b STR-haplotype in Lebanese Christians was otherwise highly specific for western Europe and was unlikely to have reached its current frequency in Lebanese Christians without admixture. We therefore suggest that the Islamic expansion from the Arabian Peninsula beginning in the seventh century CE introduced lineages typical of this area into those who subsequently became Lebanese Muslims, whereas the Crusader activity in the 11(th)-13(th) centuries CE introduced western European lineages into Lebanese Christians.

    Funded by: Wellcome Trust

    American journal of human genetics 2008;82;4;873-82

  • Meta-analysis of genome-wide association data and large-scale replication identifies additional susceptibility loci for type 2 diabetes.

    Zeggini E, Scott LJ, Saxena R, Voight BF, Marchini JL, Hu T, de Bakker PI, Abecasis GR, Almgren P, Andersen G, Ardlie K, Boström KB, Bergman RN, Bonnycastle LL, Borch-Johnsen K, Burtt NP, Chen H, Chines PS, Daly MJ, Deodhar P, Ding CJ, Doney AS, Duren WL, Elliott KS, Erdos MR, Frayling TM, Freathy RM, Gianniny L, Grallert H, Grarup N, Groves CJ, Guiducci C, Hansen T, Herder C, Hitman GA, Hughes TE, Isomaa B, Jackson AU, Jørgensen T, Kong A, Kubalanza K, Kuruvilla FG, Kuusisto J, Langenberg C, Lango H, Lauritzen T, Li Y, Lindgren CM, Lyssenko V, Marvelle AF, Meisinger C, Midthjell K, Mohlke KL, Morken MA, Morris AD, Narisu N, Nilsson P, Owen KR, Palmer CN, Payne F, Perry JR, Pettersen E, Platou C, Prokopenko I, Qi L, Qin L, Rayner NW, Rees M, Roix JJ, Sandbaek A, Shields B, Sjögren M, Steinthorsdottir V, Stringham HM, Swift AJ, Thorleifsson G, Thorsteinsdottir U, Timpson NJ, Tuomi T, Tuomilehto J, Walker M, Watanabe RM, Weedon MN, Willer CJ, Wellcome Trust Case Control Consortium, Illig T, Hveem K, Hu FB, Laakso M, Stefansson K, Pedersen O, Wareham NJ, Barroso I, Hattersley AT, Collins FS, Groop L, McCarthy MI, Boehnke M and Altshuler D

    Wellcome Trust Centre for Human Genetics, University of Oxford, Oxford, UK.

    Genome-wide association (GWA) studies have identified multiple loci at which common variants modestly but reproducibly influence risk of type 2 diabetes (T2D). Established associations to common and rare variants explain only a small proportion of the heritability of T2D. As previously published analyses had limited power to identify variants with modest effects, we carried out meta-analysis of three T2D GWA scans comprising 10,128 individuals of European descent and approximately 2.2 million SNPs (directly genotyped and imputed), followed by replication testing in an independent sample with an effective sample size of up to 53,975. We detected at least six previously unknown loci with robust evidence for association, including the JAZF1 (P = 5.0 x 10(-14)), CDC123-CAMK1D (P = 1.2 x 10(-10)), TSPAN8-LGR5 (P = 1.1 x 10(-9)), THADA (P = 1.1 x 10(-9)), ADAMTS9 (P = 1.2 x 10(-8)) and NOTCH2 (P = 4.1 x 10(-8)) gene regions. Our results illustrate the value of large discovery and follow-up samples for gaining further insights into the inherited basis of T2D.

    Funded by: Department of Health: DHCS/07/07/008; Intramural NIH HHS: Z01 HG000024; Medical Research Council: G0000649, G0600705, G0601261, MC_U106179471; NCI NIH HHS: CA87969, P01 CA087969; NHGRI NIH HHS: HG002651, N01HG65403, R01 HG002651, U01 HG004171, U01 HG004399; NHLBI NIH HHS: HL084729, U01 HL084729; NIDA NIH HHS: U54 DA021519; NIDDK NIH HHS: DK062370, DK072193, DK58845, P30 DK040561, P30 DK040561-12, R01 DK029867, R01 DK058845, R01 DK062370, R01 DK072193, R56 DK062370, U01 DK062370; Wellcome Trust: 072960, 076113, 077016, 079557, 090532, GR072960

    Nature genetics 2008;40;5;638-45