Initial sequencing and comparative analysis of the mouse genome

Mice, microbes and models of infection

We urgently need animal models to study infectious disease. Mice are susceptible to a similar range of microbial infections as humans. Marked differences between inbred strains of mice in their response to pathogen infection can be exploited to analyse the genetic basis of infections. In addition, the genetic tools that are available in the laboratory mouse, and new techniques to monitor the expression of bacterial genes in vivo, make it the principal experimental animal model for studying mechanisms of infection and immunity.

Congenic mice: cutting tools for complex immune disorders

Autoimmune diseases are, in general, under complex genetic control and subject to strong interactions between genetics and the environment. Greater knowledge of the underlying genetics will provide immunologists with a framework for study of the immune dysregulation that occurs in such diseases. Ascertaining the number of genes that are involved and their characterization have, however, proven to be difficult. Improved methods of genetic analysis and the availability of a draft sequence of the complete mouse genome have markedly improved the outlook for such research, and they have emphasized the advantages of mice as a model system. In this review, we provide an overview of the genetic analysis of autoimmune diseases and of the crucial role of congenic and consomic mouse strains in such research.

Phylogenetic shadowing of primate sequences to find functional regions of the human genome

Nonhuman primates represent the most relevant model organisms to understand the biology of Homo sapiens. The recent divergence and associated overall sequence conservation between individual members of this taxon have nonetheless largely precluded the use of primates in comparative sequence studies. We used sequence comparisons of an extensive set of Old World and New World monkeys and hominoids to identify functional regions in the human genome. Analysis of these data enabled the discovery of primate-specific gene regulatory elements and the demarcation of the exons of multiple genes. Much of the information content of the comprehensive primate sequence comparisons could be captured with a small subset of phylogenetically close primates. These results demonstrate the utility of intraprimate sequence comparisons to discover common mammalian as well as primate-specific functional elements in the human genome, which are unattainable through the evaluation of more evolutionarily distant species.

Smooth muscle uses another promoter to express primarily a form of human Cav1.2 L-type calcium channel different from the principal heart form

Several different first exons and amino termini have been reported for the cardiac Ca channel known as alpha(1C) or Ca(V)1.2. The aim of this study was to investigate whether the expression of this channel is regulated by different promoters in smooth muscle cells and in heart in humans. Ribonuclease protection assay (RPA) indicates that the longer first exon 1a is found in certain human smooth muscle-containing tissues, notably bladder and fetal aorta, but that it is not expressed to any significant degree in lung or intestine. On the other hand, all four smooth muscle-containing tissues examined strongly express transcripts containing exon 1b, first reported cloned from human fibroblast cells. In addition, primary cultures of human colonic myocytes and coronary artery smooth muscle cells express predominantly transcripts containing exon 1b. The promoter immediately upstream of exon 1b was cloned, and it displays functional promoter activity when luciferase-expressing constructs were transfected into three different cultured smooth muscle cells: primary human coronary artery smooth muscles cells, primary human colonocytes, and the fetal rat aorta-derived A7r5 cell line. These results indicate that expression in smooth muscle is primarily driven by a promoter different from that which drives expression in cardiac myocytes.

Global insights into protein complexes through integrated analysis of the reliable interactome and knockout lethality

We performed an integrated computational analysis of data derived from a comprehensive set of protein-protein interactions (interactome) and a phenotype dataset on lethality in Saccharomyces cerevisiae. For the analysis, we selected reliable interactome data using our previous 'interaction generality,' a computational approach to assess reliability of interactions. Those efforts gave clear evidence that proteins with lethal phenotypes in knockout studies (lethal proteins) may interact with each other to form functional protein complexes to perform their cellular roles. However, our analysis indicates that interactions between lethal proteins are rather restricted to the same cellular pathway or function, and it is quite unlikely that they interact with other lethal proteins functioning in different cellular roles. Furthermore, our results allowed us predictions on the functions of thus far uncharacterized lethal proteins with an estimated 93% accuracy. Thus, the analysis described in here can provide global insights into the biological features of the protein complexes.

Comparison of mouse and human genomes followed by experimental verification yields an estimated 1,019 additional genes

A primary motivation for sequencing the mouse genome was to accelerate the discovery of mammalian genes by using sequence conservation between mouse and human to identify coding exons. Achieving this goal proved challenging because of the large proportion of the mouse and human genomes that is apparently conserved but apparently does not code for protein. We developed a two-stage procedure that exploits the mouse and human genome sequences to produce a set of genes with a much higher rate of experimental verification than previously reported prediction methods. RT-PCR amplification and direct sequencing applied to an initial sample of mouse predictions that do not overlap previously known genes verified the regions flanking one intron in 139 predictions, with verification rates reaching 76%. On average, the confirmed predictions show more restricted expression patterns than the mouse orthologs of known human genes, and two-thirds lack homologs in fish genomes, demonstrating the sensitivity of this dual-genome approach to hard-to-find genes. We verified 112 previously unknown homologs of known proteins, including two homeobox proteins relevant to developmental biology, an aquaporin, and a homolog of dystrophin. We estimate that transcription and splicing can be verified for >1,000 gene predictions identified by this method that do not overlap known genes. This is likely to constitute a significant fraction of the previously unknown, multiexon mammalian genes.

Tag, you're hit: retroviral insertions identify genes involved in cancer

Retroviral integrations have been used for many years to identify genes involved in cancer. The recently published mouse genome sequence has allowed large-scale identification of potential human cancer genes and their classification into distinct signaling pathways.

PRISM, a generic large scale proteomic investigation strategy for mammals

We have developed a systematic analytical approach, termed PRISM (Proteomic Investigation Strategy for Mammals), that permits routine, large scale protein expression profiling of mammalian cells and tissues. PRISM combines subcellular fractionation, multidimensional liquid chromatography-tandem mass spectrometry-based protein shotgun sequencing, and two newly developed computer algorithms, STATQUEST and GOClust, as a means to rapidly identify, annotate, and categorize thousands of expressed mammalian proteins. The application of PRISM to adult mouse lung and liver resulted in the high confidence identification of over 2,100 unique proteins including more than 100 integral membrane proteins, 400 nuclear proteins, and 500 uncharacterized proteins, the largest proteome study carried out to date on this important model organism. Automated clustering of the identified proteins into Gene Ontology annotation groups allowed for streamlined analysis of the large data set, revealing interesting and physiologically relevant patterns of tissue and organelle specificity. PRISM therefore offers an effective platform for in-depth investigation of complex mammalian proteomes.

Genomic approaches to identification of tumour-susceptibility genes using mouse models

Individual susceptibility to cancer in humans is determined by complex interactions between germline genetic variation and levels of exposure to environmental carcinogens or tumour promoters. Only a small fraction of cancer susceptibility is inherited in a Mendelian manner (high-penetrance familial cancer), and most tumours result from the combined effects of many gene-gene and gene-environment interactions. The sequencing of the mouse genome provides new approaches to one of the most challenging tasks of cancer genetics today.

Refinement of a chimpanzee pericentric inversion breakpoint to a segmental duplication cluster

BACKGROUND

Pericentric inversions are the most common euchromatic chromosomal differences among humans and the great apes. The human and chimpanzee karyotype differs by nine such events, in addition to several constitutive heterochromatic increases and one chromosomal fusion event. Reproductive isolation and subsequent speciation are thought to be the potential result of pericentric inversions, as reproductive boundaries form as a result of hybrid sterility.

RESULTS

Here we employed a comparative fluorescence in situ hybridization approach, using probes selected from a combination of physical mapping, genomic sequence, and segmental duplication analyses to narrow the breakpoint interval of a pericentric inversion in chimpanzee involving the orthologous human 15q11-q13 region. We have refined the inversion breakpoint of this chimpanzee-specific rearrangement to a 600 kilobase (kb) interval of the human genome consisting of entirely duplicated material. Detailed analysis of the underlying sequence indicated that this region comprises multiple segmental duplications, including a previously characterized duplication of the alpha7 neuronal nicotinic acetylcholine receptor subunit gene (CHRNA7) in 15q13.3 and several Golgin-linked-to-PML, or LCR15, duplications.

CONCLUSIONS

We conclude that, on the basis of experimental data excluding the CHRNA7 duplicon as the site of inversion, and sequence analysis of regional duplications, the most likely rearrangement site is within a GLP/LCR15 duplicon. This study further exemplifies the genomic plasticity due to the presence of segmental duplications and highlights their importance for a complete understanding of genome evolution.

Of mice and men: phylogenetic footprinting aids the discovery of regulatory elements

Phylogenetic footprinting is an approach to finding functionally important sequences in the genome that relies on detecting their high degrees of conservation across different species. A new study shows how much it improves the prediction of gene-regulatory elements in the human genome.

Time scale of eutherian evolution estimated without assuming a constant rate of molecular evolution

Controversies over the molecular clock hypothesis were reviewed. Since it is evident that the molecular clock does not hold in an exact sense, accounting for evolution of the rate of molecular evolution is a prerequisite when estimating divergence times with molecular sequences. Recently proposed statistical methods that account for this rate variation are overviewed and one of these procedures is applied to the mitochondrial protein sequences and to the nuclear gene sequences from many mammalian species in order to estimate the time scale of eutherian evolution. This Bayesian method not only takes account of the variation of molecular evolutionary rate among lineages and among genes, but it also incorporates fossil evidence via constraints on node times. With denser taxonomic sampling and a more realistic model of molecular evolution, this Bayesian approach is expected to increase the accuracy of divergence time estimates.

The mouse genome sequence--the end of the tail, or just the beginning?

The recent flurry of papers on the mouse genome includes the description of the full genome assembly, analysis of the mouse transcriptome, the origin of interstrain variation, initial analysis of conserved non-coding regions and high-throughput expression analysis of a subset of genes. Each illustrates how the availability of the genome sequence will change the way mouse biologists do business in future.

The phusion assembler

The Phusion assembler has assembled the mouse genome from the whole-genome shotgun (WGS) dataset collected by the Mouse Genome Sequencing Consortium, at ~7.5x sequence coverage, producing a high-quality draft assembly 2.6 gigabases in size, of which 90% of these bases are in 479 scaffolds. For the mouse genome, which is a large and repeat-rich genome, the input dataset was designed to include a high proportion of paired end sequences of various size selected inserts, from 2-200 kbp lengths, into various host vector templates. Phusion uses sequence data, called reads, and information about reads that share common templates, called read pairs, to drive the assembly of this large genome to highly accurate results. The preassembly stage, which clusters the reads into sensible groups, is a key element of the entire assembler, because it permits a simple approach to parallelization of the assembly stage, as each cluster can be treated independent of the others. In addition to the application of Phusion to the mouse genome, we will also present results from the WGS assembly of Caenorhabditis briggsae sequenced to about 11x coverage. The C. briggsae assembly was accessioned through EMBL, http://www.ebi.ac.uk/services/index.html, using the series CAAC01000001-CAAC01000578, however, the Phusion mouse assembly described here was not accessioned. The mouse data was generated by the Mouse Genome Sequencing Consortium. The C. briggsae sequence was generated at The Wellcome Trust Sanger Institute and the Genome Sequencing Center, Washington University School of Medicine.

Seeing chordate evolution through the Ciona genome sequence

A draft sequence of the compact genome of the sea squirt Ciona intestinalis, a non-vertebrate chordate that diverged very early from other chordates, including vertebrates, illuminates how chordates originated and how vertebrate developmental innovations evolved.

Covariation in frequencies of substitution, deletion, transposition, and recombination during eutherian evolution

Six measures of evolutionary change in the human genome were studied, three derived from the aligned human and mouse genomes in conjunction with the Mouse Genome Sequencing Consortium, consisting of (1) nucleotide substitution per fourfold degenerate site in coding regions, (2) nucleotide substitution per site in relics of transposable elements active only before the human-mouse speciation, and (3) the nonaligning fraction of human DNA that is nonrepetitive or in ancestral repeats; and three derived from human genome data alone, consisting of (4) SNP density, (5) frequency of insertion of transposable elements, and (6) rate of recombination. Features 1 and 2 are measures of nucleotide substitutions at two classes of "neutral" sites, whereas 4 is a measure of recent mutations. Feature 3 is a measure dominated by deletions in mouse, whereas 5 represents insertions in human. It was found that all six vary significantly in megabase-sized regions genome-wide, and many vary together. This indicates that some regions of a genome change slowly by all processes that alter DNA, and others change faster. Regional variation in all processes is correlated with, but not completely accounted for, by GC content in human and the difference between GC content in human and mouse.

Recent segmental and gene duplications in the mouse genome

BACKGROUND

The high quality of the mouse genome draft sequence and its associated annotations are an invaluable biological resource. Identifying recent duplications in the mouse genome, especially in regions containing genes, may highlight important events in recent murine evolution. In addition, detecting recent sequence duplications can reveal potentially problematic regions of the genome assembly. We use BLAST-based computational heuristics to identify large (>/= 5 kb) and recent (>/= 90% sequence identity) segmental duplications in the mouse genome sequence. Here we present a database of recently duplicated regions of the mouse genome found in the mouse genome sequencing consortium (MGSC) February 2002 and February 2003 assemblies.

RESULTS

We determined that 33.6 Mb of 2,695 Mb (1.2%) of sequence from the February 2003 mouse genome sequence assembly is involved in recent segmental duplications, which is less than that observed in the human genome (around 3.5-5%). From this dataset, 8.9 Mb (26%) of the duplication content consisted of 'unmapped' chromosome sequence. Moreover, we suspect that an additional 18.5 Mb of sequence is involved in duplication artifacts arising from sequence misassignment errors in this genome assembly. By searching for genes that are located within these regions, we identified 675 genes that mapped to duplicated regions of the mouse genome. Sixteen of these genes appear to have been duplicated independently in the human genome. From our dataset we further characterized a 42 kb recent segmental duplication of Mater, a maternal-effect gene essential for embryogenesis in mice.

CONCLUSION

Our results provide an initial analysis of the recently duplicated sequence and gene content of the mouse genome. Many of these duplicated loci, as well as regions identified to be involved in potential sequence misassignment errors, will require further mapping and sequencing to achieve accuracy. A Genome Browser database was set up to display the identified duplication content presented in this work. This data will also be relevant to the growing number of investigators who use the draft genome sequence for experimental design and analysis.

Comparative gene prediction in human and mouse

The completion of the sequencing of the mouse genome promises to help predict human genes with greater accuracy. While current ab initio gene prediction programs are remarkably sensitive (i.e., they predict at least a fragment of most genes), their specificity is often low, predicting a large number of false-positive genes in the human genome. Sequence conservation at the protein level with the mouse genome can help eliminate some of those false positives. Here we describe SGP2, a gene prediction program that combines ab initio gene prediction with TBLASTX searches between two genome sequences to provide both sensitive and specific gene predictions. The accuracy of SGP2 when used to predict genes by comparing the human and mouse genomes is assessed on a number of data sets, including single-gene data sets, the highly curated human chromosome 22 predictions, and entire genome predictions from ENSEMBL. Results indicate that SGP2 outperforms purely ab initio gene prediction methods. Results also indicate that SGP2 works about as well with 3x shotgun data as it does with fully assembled genomes. SGP2 provides a high enough specificity that its predictions can be experimentally verified at a reasonable cost. SGP2 was used to generate a complete set of gene predictions on both the human and mouse by comparing the genomes of these two species. Our results suggest that another few thousand human and mouse genes currently not in ENSEMBL are worth verifying experimentally.

Hemoglobinopathies

The outlook for patients with sickle cell disease has improved steadily during the last two decades. In spite of these improvements, curative therapies are currently available only to a small minority of patients. The main theme of this chapter is to describe new therapeutic options that are at different stages of development that might result in further improvements in the outlook for patients with these disorders.

Dr. Joseph DeSimone and his colleagues had previously made the important observation that the hypomethylating agent 5-azacytidine can reverse the switch from adult to fetal hemoglobin in adult baboons. Although similar activity was demonstrated in patients with sickle cell disease and β-thalassemia, concern about the toxicity of 5-azacytidine prevented its widespread use in these disorders. In Section I, Dr. DeSimone discusses the role of DNA methylation in globin gene regulation and describe recent clinical experience with decitabine (an analogue of 5-azacytidine) in patients with sickle cell disease. These encouraging studies demonstrate significant fetal hemoglobin inducing activity of decitabine in patients who fail to respond to hydroxyurea.

In Section II, Dr. George Atweh continues the same theme by describing recent progress in the study of butyrate, another inducer of fetal hemoglobin, in patients with sickle cell disease and β-thalassemia. The main focus of his section is on the use of a combination of butyrate and hydroxyurea to achieve higher levels of fetal hemoglobin that might be necessary for complete amelioration of the clinical manifestations of these disorders. Dr. Atweh also describes novel laboratory studies that shed new light on the mechanisms of fetal hemoglobin induction by butyrate.

In Section III, Dr. Ronald Nagel discusses the different available transgenic sickle mice as experimental models for human sickle cell disease. These experimental models have already had a significant impact on our understanding of the pathophysiology of sickle cell disease. Dr. Nagel describes more recent studies in which transgenic sickle mice provide the first proof of principle that globin gene transfer into hematopoietic stem cells inhibits in vivo sickling and ameliorates the severity of the disease.

Although stroke in adult patients with sickle cell disease is not as common as in children, adult hematologists, like their pediatric colleagues, need to make management decisions in adult patients with a stroke or a history of stroke. Dr. Robert Adams has led several large clinical studies that investigated the role of transfusions in the prevention of stroke in children with sickle cell disease. Much less is known, however, about the prevention of first or subsequent strokes in adult patients with sickle cell disease. In Section IV, Dr. Adams provides some general guidelines for the management of adult patients with stroke while carefully distinguishing between recommendations that are evidence-based and those that are anecdotal in nature.

Genome rearrangements in mammalian evolution: lessons from human and mouse genomes

Although analysis of genome rearrangements was pioneered by Dobzhansky and Sturtevant 65 years ago, we still know very little about the rearrangement events that produced the existing varieties of genomic architectures. The genomic sequences of human and mouse provide evidence for a larger number of rearrangements than previously thought and shed some light on previously unknown features of mammalian evolution. In particular, they reveal that a large number of microrearrangements is required to explain the differences in draft human and mouse sequences. Here we describe a new algorithm for constructing synteny blocks, study arrangements of synteny blocks in human and mouse, derive a most parsimonious human-mouse rearrangement scenario, and provide evidence that intrachromosomal rearrangements are more frequent than interchromosomal rearrangements. Our analysis is based on the human-mouse breakpoint graph, which reveals related breakpoints and allows one to find a most parsimonious scenario. Because these graphs provide important insights into rearrangement scenarios, we introduce a new visualization tool that allows one to view breakpoint graphs superimposed with genomic dot-plots.

Strategies and tools for whole-genome alignments

The availability of the assembled mouse genome makes possible, for the first time, an alignment and comparison of two large vertebrate genomes. We investigated different strategies of alignment for the subsequent analysis of conservation of genomes that are effective for assemblies of different quality. These strategies were applied to the comparison of the working draft of the human genome with the Mouse Genome Sequencing Consortium assembly, as well as other intermediate mouse assemblies. Our methods are fast and the resulting alignments exhibit a high degree of sensitivity, covering more than 90% of known coding exons in the human genome. We obtained such coverage while preserving specificity. With a view towards the end user, we developed a suite of tools and Web sites for automatically aligning and subsequently browsing and working with whole-genome comparisons. We describe the use of these tools to identify conserved non-coding regions between the human and mouse genomes, some of which have not been identified by other methods.

11th Intelligent Systems for Molecular Biology 2003 (ISMB 2003)

This report profiles the keynote talks given at ISMB03 in Brisbane, Australia by Ron Shamir, David Haussler, John Mattick, Yoshihide Hayashizaki, Sydney Brenner, the Overton Prize winner, Jim Kent, and the ISCB Senior Accomplishment Awardee, David Sankov.

AVID: A global alignment program

In this paper we describe a new global alignment method called AVID. The method is designed to be fast, memory efficient, and practical for sequence alignments of large genomic regions up to megabases long. We present numerous applications of the method, ranging from the comparison of assemblies to alignment of large syntenic genomic regions and whole genome human/mouse alignments. We have also performed a quantitative comparison of AVID with other popular alignment tools. To this end, we have established a format for the representation of alignments and methods for their comparison. These formats and methods should be useful for future studies. The tools we have developed for the alignment comparisons, as well as the AVID program, are publicly available. See Web Site References section for AVID Web address and Web addresses for other programs discussed in this paper.

Cross-species sequence comparisons: a review of methods and available resources

With the availability of whole-genome sequences for an increasing number of species, we are now faced with the challenge of decoding the information contained within these DNA sequences. Comparative analysis of DNA sequences from multiple species at varying evolutionary distances is a powerful approach for identifying coding and functional noncoding sequences, as well as sequences that are unique for a given organism. In this review, we outline the strategy for choosing DNA sequences from different species for comparative analyses and describe the methods used and the resources publicly available for these studies.

Leveraging the mouse genome for gene prediction in human: from whole-genome shotgun reads to a global synteny map

The availability of draft sequences for both the mouse and human genomes makes it possible, for the first time, to annotate whole mammalian genomes using comparative methods. TWINSCAN is a gene-prediction system that combines the methods of single-genome predictors like GENSCAN with information derived from genome comparison, thereby improving accuracy. Because TWINSCAN uses genomic sequence only, it is less biased toward highly and/or ubiquitously expressed genes than GENEWISE, GENOMESCAN, and other methods based on evidence derived from transcripts. We show that TWINSCAN improves gene prediction in human using intermediate products from various stages of the sequencing and analysis of the mouse genome, from low-redundancy, whole-genome shotgun reads to the draft assembly and the synteny map. TWINSCAN improves on the prior state of the art even when alignments from only 1X coverage of the mouse genome are available. Gene prediction accuracy improves steadily from 1X through 3X, more slowly from 3X to 4X, and relatively little thereafter. The assembly and the synteny map greatly speed the computations, however. Our human annotation using the mouse assembly is conservative, predicting only 25,622 genes, and appears to be one of the best de novo annotations of the human genome to date.

Distinguishing regulatory DNA from neutral sites

We explore several computational approaches to analyzing interspecies genomic sequence alignments, aiming to distinguish regulatory regions from neutrally evolving DNA. Human-mouse genomic alignments were collected for three sets of human regions: (1) experimentally defined gene regulatory regions, (2) well-characterized exons (coding sequences, as a positive control), and (3) interspersed repeats thought to have inserted before the human-mouse split (a good model for neutrally evolving DNA). Models that potentially could distinguish functional noncoding sequences from neutral DNA were evaluated on these three data sets, as well as bulk genome alignments. Our analyses show that discrimination based on frequencies of individual nucleotide pairs or gaps (i.e., of possible alignment columns) is only partially successful. In contrast, scoring procedures that include the alignment context, based on frequencies of short runs of alignment columns, dramatically improve separation between regulatory and neutral features. Such scoring functions should aid in the identification of putative regulatory regions throughout the human genome.

Impact of Pharmacogenomics on Neglected Diseases of the Developing World

Pharmacogenomics promises to have an important impact on the major health problems of the developing world, especially on neglected infectious diseases such as malaria, tuberculosis, and HIV/AIDS. Its capacity to identify new targets for drug development, together with its potential application in identifying populations who will respond favorably to a particular drug, gives it a unique place as a technology to bridge the genomics divide between rich and poor nations. To realize its true potential, however, significant scientific, legal, ethical, political, and economic challenges need to be overcome. For this to occur, an innovative global approach based on strong collaboration between industry, academia, non-governmental, and international organizations will be required. Simultaneously, more equitable and active participation from developing country researchers themselves is critical in overcoming these challenges.

Human-mouse alignments with BLASTZ

The Mouse Genome Analysis Consortium aligned the human and mouse genome sequences for a variety of purposes, using alignment programs that suited the various needs. For investigating issues regarding genome evolution, a particularly sensitive method was needed to permit alignment of a large proportion of the neutrally evolving regions. We selected a program called BLASTZ, an independent implementation of the Gapped BLAST algorithm specifically designed for aligning two long genomic sequences. BLASTZ was subsequently modified, both to attain efficiency adequate for aligning entire mammalian genomes and to increase its sensitivity. This work describes BLASTZ, its modifications, the hardware environment on which we run it, and several empirical studies to validate its results.