Deletions in the genomes of fifteen inbred mouse lines and their possible implications for fat accumulation

Deletions in the genomes of fifteen inbred mouse lines and their possible implications for fat accumulation 1 first-author †E-mail:armin.schmitt@agrar.hu-berlin.de †E-mail:armin.schmitt@agrar.hu-berlin.de armin.schmitt@agrar.hu-berlin.de 2 1 Institute for Animal Sciences, Humboldt-Universität zu Berlin, Invalidenstraße 42, 10115 Berlin, Germany Institute for Animal Sciences, Humboldt-Universität zu Berlin, Invalidenstraße 42, 10115 Berlin, Germany Institute for Medical Biometry and Epidemiology, Phillips-Universität Marburg, Bunsenstraße 3, 35037 Marburg, Germany Institute for Medical Biometry and Epidemiology, Phillips-Universität Marburg, Bunsenstraße 3, 35037 Marburg, Germany Copy number variants (CNVs) are pieces of genomic DNA of 1000 base pairs or longer which occur in a given genome at a different frequency than in a reference genome. Their importance as a source for phenotypic variability has been recognized only in the last couple of years. Chromosomal deletions can be seen as a special case of CNVs where stretches of DNA are missing in certain lines when compared to the reference genome of the mouse line C57BL/6, for example. Based upon more than 8 million single nucleotide polymorphisms (SNPs) in the fifteen inbred mouse lines which were determined in a whole genome chip based resequencing project by Perlegen Sciences, we detected 20 166 such long chromosomal deletions. They cover altogether between 4.4 million and 8.8 million base pairs, depending on the mouse line. Thus, their extent is comparable to that of SNPs. The chromosomal deletions were found by searching for clusters of missing values in the genotyping data by applying bioinformatics and biostatistical methods. In contrast to isolated missing values, clusters are likely the consequence of missing DNA probe rather than of a failed hybridization or deficient oligos. We analyzed these deletion sites in various ways. Twenty-two percent of these deletion sites overlap with exons; they could therefore affect a gene’s functioning. The corresponding genes seem to exist in alternative forms, a phenomenon that reminds of the alternative forms of mRNA generated during gene splicing. We furthermore detected statistically significant association between hundreds of deletion sites and fat weight at the age of eight weeks. Copy number variants (CNVs) Chromosomal deletions Single nucleotide polymorphisms (SNPs) Resequencing Cluster analysis Association between genotype and phenotype INTRODUCTION 1 Copy number variants (CNVs) gained much attraction in the last couple of years. In two seminal papers 76 and 255 CNVs, respectively, were identified in a genome wide study in human (Sebat et al., 2004 2004 ; Iafrate et al., 2004 2004 ). On a chromosome wide scale, CNVs were identified using the technique of comparative genome hybridization (CGH) about a decade ago (Pinkel et al., 1998 1998 ). CNVs were found to influence disease susceptibility, e.g., for AIDS if the number of copies of the chemokine receptor gene CCL4L1 is too low (Gonzalez et al., 2005 2005 ) or to cause disease directly. One example for the latter would be the well known Prader-Willi-syndrome, where a partially deleted chromosome 15 causes a range of severe physical and mental disablements in human. Deletions were also found to be the key driving force in the divergence between humans and chimpanzees (Britten, 2002 2002 ). They should therefore be seen as one of the main factors driving evolution in general. Deletions in the genome can be considered as a special case of CNVs, where one individual has no copy of a certain piece of DNA and other individuals have one or several copies. Pieces of DNA can above all get lost through unequal crossing over during replication (Graur and Li, 1991 1991 ). Other mechanisms leading to deletions are presented by Stankiewicz and Lupski ( 2002 2002 ). The deletions and other CNVs found in various studies are collected in the Database of Genomic Variants (Zhang et al., 2006 2006 ). While uncalled alleles, so called null genotypes, are usually excluded from further analysis in medical studies, it was shown that they can be exploited to derive deletions from single nucleotide polymorphism (SNP) data (McCarroll et al., 2006 2006 ; Conrad et al., 2006 2006 ). The basic idea is that uncalled alleles should cluster in an ordered genotype data set if they are caused by missing DNA probe, but they should be randomly scattered in the data set if they are caused by occasionally failed oligos or other sources of noise (Fig. 1 1 ). Here, we apply a clustering algorithm to distinguish between the two types of uncalled alleles. Obvious questions that we address are: Can putatively functional DNA get lost? Can deletions affect the phenotpye? MATERIALS AND METHODS 1 SNP data 2 We downloaded more than 8 million SNPs which were detected in a resequencing project at Perlegen Sciences (Frazer et al., 2007 2007 ) ( http://mouse.perlegen.com/mouse http://mouse.perlegen.com/mouse ). In this project, all of the unique mouse genome was probed with oligomers derived from the line C57BL/6. Repetitive sequence cannot be examined this way, so 1 500 megabases resulted. Unmapped SNPs and SNPs in amplificates whose primer pairs were marked as failed by Perlegen Scinences were excluded. The fifteen mouse lines whose genomes were resequenced are standard mouse lines (Fig. 2 2 ). Algorithm to identify deletion sites 2 We considered uncalled alleles as caused by missing DNA probe due to deletions if the following three criteria were met: 1. The same mouse lines have three or more uncalled alleles at consecutive SNP positions. 2. The uncalled alleles span 1 kb ore more. 3. The distance between neighbouring SNPs does not exceed 1 kb. The last requirement had to be met in order not to bridge too much of non-genotyped DNA. We mapped the genotypes onto a two-letter alphabet. One symbol stands for the four nucleotides A, C, G, and T. The other stands for uncalled alleles. We were then looking for stretches of at least three SNPs where the identical mouse lines have uncalled alleles. All such stretches were extracted and subsequently criteria 2 and 3 were applied to extract the deletion sites. In order to assess if deletion sites could also be observed in the case where the uncalled alleles were randomly scattered across the mouse lines, we generated a dataset of scrambled genotypes, i.e., we fixed the number of uncalled alleles that were found at a specific chromosomal position, but assigned them randomly to the mouse lines. Then the search algorithm was applied to this scrambled dataset. Association between deletion sites and quantitative traits 2 We downloaded the measurements for average fat content at the age of eight weeks from the Mouse Phenome Database ( http://www.jax.org/phenome http://www.jax.org/phenome ) (Bogue et al., 2007 2007 ). The symbol for this trait was ‘tissuemass_fat8’ and the filename containing all the measurements was ‘summstats.zip’. We extracted the corresponding measurements for male animals. Phenotypic measurements for all genotyped lines were available except for the line PWD/PhJ. For each deletion site two groups of mouse lines were formed, one group with the deletion and the other without the deletion. We tested all ~20 000 deletion sites for association using the t -test for clustered data (Donner et al., 1981 1981 ), where each line was considered as one cluster. To this end, the R-function t .test.cluster in library Hmisc by Frank Harrell was adopted for our needs. The P -values were corrected using the method of Bonferroni. Analysis of the chromosomal location of deletion sites 2 In order to examine if exons were affected by deletions, we compared their positions with the positions of all exons in the genome database of Ensembl release 46 ( http://www.ensembl.org http://www.ensembl.org ) using bioperl modules. The position of the deletion sites are defined by the chromosomal position of their first and their last SNP. We confirmed these positions by sequence comparison of the SNP flanking sequence of 100 bp upstream and 100 bp downstream of the SNP against the mouse genome (National Center for Biotechnology Information (NCBI) build 36, Ensembl database release 46) using BLAST (Altschul et al., 1997 1997 ). When the first and the last SNP flanking sequences could be mapped to the same chromosomal locations as indicated in the original data file with an E -value of 10 −50 or better, we accepted a deletion site for further analysis. RESULTS 1 Applying the algorithm sketched in the section MATERIALS AND METHODS we identified 20 885 deletion sites. The position of 20 166 deletion sites with a total length of 31.8 Mb could be confirmed by mapping the sequence of the first and the last SNPs of each deletion site onto the current assembly of the mouse genome (NCBI build 36, Ensembl database release 46). In our scrambled data set we only found 52 deletion sites with a total length of 65 kb. This is only a minute proportion of the results obtained from real data. There is, thus, good evidence that the deletion sites that we found are no artefacts. On average, 5.4 Mb were lost in a mouse line, with a minimum of 4.4 Mb in the line A/J and a maximum of 9.3 Mb in the line PWD/PhJ. Next, we were interested if exclusively so called junk DNA gets lost or if also presumably functional DNA is affected. To this end, we compared the positions of the deletions sites, as defined by the position of the first and the last SNPs in it, with the positions of the exons in the genome. Surprisingly, 2 594 (12.9%) of the deletion sites overlap with at least one exon. Altogether, exons of 2 104 out of 28 545 confirmed and predicted mouse genes (7.3%) and a total of 4 121 murine exons (1.75% of all 235 412 exons, including predicted ones), are affected by deletions in at least one mouse line. We tested all 20 166 deletion sites with confirmed location for association with the fat weight at the age of eight weeks. Two thousand two hundred and sixty-six deletion sites had a significant correlation with this trait (significance level α =0.001 after Bonferroni correction). As an example, we show the fat distribution for the two groups of lines, one with the deletion, the other without the deletion, at the chromosomal location between 93 339 270 bp and 93 340 305 bp on chromosome 13 (Fig. 3 3 ). The first SNP of this deletion site is ss52681234, the last ss52681232. This deletion site overlaps the gene Msh3 , which was previously associated with DNA mismatch repair (Edelmann et al., 2000 2000 ). DISCUSSION AND CONCLUSION 1 We have shown that uncalled alleles obtained in a genome-wide SNP experiment lend themselves to the identification of genomic deletions. A straightforward algorithm which is largely based upon the coherence between uncalled alleles in mouse lines, is capable of discerning clusters of uncalled alleles from uncalled alleles which are scattered across the data. As a first validation, we scrambled the genotype data and obtained only a tiny fraction of the deletion sites that we obtained in the real data. The extent of the deletions that we extracted in 15 inbred mouse lines, on average 5.4 Mb, is comparable to that of the SNPs, which was 8 million. Of course, we are aware that a rigorous validation of the deletion sites that we identified would include sequencing of mouse lines with and without the predicted deletion. A further improvement of this analysis would consist in the merging of neighbouring deletion sites. Since we wanted to be on the safe side, we did not bridge regions that were not genotyped. This leads to a strong fragmentation of the putatively longer deletion sites. Another validation would be the comparison of the deletion sites that we detected with those that Li et al.( 2004 2004 ) found in 14 inbred lines using whole-genome BAC arrays. We found furthermore that a significant proportion of the deletion sites overlaps with exons, which should lead to functional consequences. As an example, we were looking for association between the deletion sites and the fat content of mice at the age of eight weeks. About 11% of the deletion sites showed a significant association with this trait. There is, therefore, evidence that chromosomal deletions play an important role for quantitative traits in general. Next to SNPs and epigenetic changes, their careful analyses should gain more attention. References SF Altschul , TL Madden , AA Schaffer , J Zhang , Z Zhang , W Miller , DJ Lipman . Gapped BLAST and PSI-BLAST: a new generation of protein database search programs . Nucleic Acids Res , 1997 . 25 ( 17 ): 3389 - 3402 . DOI: 10.1093/nar/25.17.3389 http://doi.org/10.1093/nar/25.17.3389 . Gapped BLAST and PSI-BLAST: a new generation of protein database search programs MA Bogue , SC Grubb , TP Maddatu , CJ Bult . Mouse Phenome Database (MPD) . Nucleic Acids Res , 2007 . 35 ( Database issue ): 643 - 649 . DOI: 10.1093/nar/gkl1049 http://doi.org/10.1093/nar/gkl1049 . Mouse Phenome Database (MPD) RJ Britten . Divergence between samples of chimpanzee and human DNA sequences is 5%, counting indels . PNAS , 2002 . 99 ( 21 ): 13633 - 13635 . DOI: 10.1073/pnas.172510699 http://doi.org/10.1073/pnas.172510699 . Divergence between samples of chimpanzee and human DNA sequences is 5%, counting indels FD Conrad , TD Andrews , NP Carter , ME Hurles , JK Pritchard . A high-resolution survey of deletion polymorphism in the human genome . Nat Genet , 2006 . 38 ( 1 ): 75 - 81 . DOI: 10.1038/ng1697 http://doi.org/10.1038/ng1697 . A high-resolution survey of deletion polymorphism in the human genome A Donner , N Birkett , C Buck . Randomization by cluster . Am J Epidemiol , 1981 . 114 ( 6 ): 906 - 914 . . Randomization by cluster W Edelmann , A Umar , K Yang , J Heyer , M Kucherlapati , M Lia , B Kneitz , E Avdievich , K Fan , E Wong , 等 . The DNA mismatch repair genes Msh3 and Msh6 cooperate in intestinal tumor suppression . Cancer Res , 2000 . 60 ( 4 ): 803 - 807 . . The DNA mismatch repair genes Msh3 and Msh6 cooperate in intestinal tumor suppression KA Frazer , E Eskin , HM Kang , MA Bogue , DA Hinds , EJ Beilharz , RV Gupta , J Montgomery , MM Morenzoni , GB Nilsen , 等 . A sequence-based variation map of 8.27 million SNPs in inbred mouse strains . Nature , 2007 . 448 ( 7157 ): 1050 - 1053 . DOI: 10.1038/nature06067 http://doi.org/10.1038/nature06067 . A sequence-based variation map of 8.27 million SNPs in inbred mouse strains C Gonzalez , H Kulkarni , H Bolivar , A Mangano , R Sanchez , G Catano , RJ Nibbs , BI Freedman , MP Quinones , MJ Bamshad , 等 . The influence of CCL3L1 gene-containing segmental duplications on HIV-1/AIDS susceptibility . Science , 2005 . 307 ( 5714 ): 1434 - 1440 . DOI: 10.1126/science.1101160 http://doi.org/10.1126/science.1101160 . The influence of CCL3L1 gene-containing segmental duplications on HIV-1/AIDS susceptibility D Graur , WH Li . Fundamentals of Molecular Evolution , 1991 . : Sunderland, MA, USA Sinauer Associates , 15 . Fundamentals of Molecular Evolution AJ Iafrate , L Feuk , MN Rivera , ML Listewnik , PD Donahoe , Y Qi , SW Scherer , C Lee . Detection of large-scale variation in the human genome . Nat Genet , 2004 . 36 ( 9 ): 949 - 951 . DOI: 10.1038/ng1416 http://doi.org/10.1038/ng1416 . Detection of large-scale variation in the human genome J Li , T Jiang , JH Mao , A Balmain , L Peterson , C Harris , PH Rao , P Havlak , R Gibbs , WW Cai . Genomic segmental polymorphisms in inbred mouse strains . Nat Genet , 2004 . 36 ( 9 ): 952 - 954 . DOI: 10.1038/ng1417 http://doi.org/10.1038/ng1417 . Genomic segmental polymorphisms in inbred mouse strains SA McCarroll , TN Hadnott , GH Perry , PC Sabeti , MC Zody , JC Barrett , S Dalaire , SB Gariel , C Lee , MJ Daly , 等 . Common deletion polymorphims in the human genome . Nat Genet , 2006 . 38 ( 1 ): 86 - 92 . DOI: 10.1038/ng1696 http://doi.org/10.1038/ng1696 . Common deletion polymorphims in the human genome D Pinkel , R Segraves , D Sudar , S Clark , I Poole , D Kowbel , C Collins , WL Kuo , C Chen , Y Zhai , 等 . High resolution analysis of DNA copy number variation using comparative genomic hybridization to microarrays . Nat Genet , 1998 . 20 ( 2 ): 207 - 211 . DOI: 10.1038/2524 http://doi.org/10.1038/2524 . High resolution analysis of DNA copy number variation using comparative genomic hybridization to microarrays J Sebat , B Lakshmi , J Troge , J Alexander , J Young , P Lundin , S Mänér , H Massa , M Walker , M Chi , 等 . Large-scale copy number polymorphism in the human genome . Science , 2004 . 305 ( 5683 ): 525 - 528 . DOI: 10.1126/science.1098918 http://doi.org/10.1126/science.1098918 . Large-scale copy number polymorphism in the human genome P Stankiewicz , JR Lupski . Genome architecture, rearrangements and genomic disorders . TRENDS in Genetics , 2002 . 18 ( 2 ): 74 - 82 . DOI: 10.1016/S0168-9525(02)02592-1 http://doi.org/10.1016/S0168-9525(02)02592-1 . Genome architecture, rearrangements and genomic disorders J Zhang , L Feuk , GE Duggan , R Khaja , SW Scherer . Development of bioinformatics resources for display and analysis of copy number and other structural variants in the human genome . Cytogenet Genome Res , 2006 . 115 ( 3-4 ): 205 - 214 . DOI: 10.1159/000095916 http://doi.org/10.1159/000095916 . Development of bioinformatics resources for display and analysis of copy number and other structural variants in the human genome We thank Ivo Große (University of Halle, Germany), Richard Mott (Wellcome Trust Centre for Human Genetics, Oxford, UK), Geoff Nilsen and Brian Karlak (both at Perlegen Sciences, Mountain View, USA), Bert Overduin and Xose Fernandez (both at European Bioinformatics Institute, Hinxton, UK) for stimulating discussions and valuable assistance. 10.1631/jzus.2007.B0777 Project supported by the German Ministry of Education and Research (BMBF) through the National Genome Research Network (NGFN) (Nos. 01GS0486 and 01GR0460) and the Deutsche Forschungsgemeinschaft (DFG) for a Travel Grant to Armin O. Schmitt 2007-09-16 2007-09-18 2022-01-05 Copyright © 2007, Journal of Zhejiang University Science 2007 Journal of Zhejiang University-SCIENCE B (Biomedicine and Biotechnology) 11 08