A comparison of whole-genome shotgun-derived mouse chromosome 16 and the human genome

Establishment of a sepsis model following implantation ofKlebsiella pneumoniae-infected fibrin clot into the peritoneal cavity of mice

Successful establishment of sepsis by entrapping a dose of 150 colony forming units of Klebsiella pneumoniae in a fibrin clot following implantation into the peritoneal cavity of mice is reported. The dose in the fibrin clot gave 50% mortality in mice, spread over a period of one week. All the infected mice showed positive blood culture up to 6 d post-infection; histopathology revealed inflammatory changes in both liver and spleen. Introduction of K. pneumoniae into experimental mice without entrapment in fibrin clot caused no mortality and blood culture remained positive only up to 2 d; histopathology of liver and spleen throughout the period of study showed relatively mild inflammatory changes, which almost cleared during 14 d post-infection. The use of the fibrin-clot model may thus be considered to be useful in studying both the initial and the persisting stage of infection in the peritoneum, whence a slow release of bacteria into the blood takes place which finally leads to sepsis and septicemia.

Automated Gene Ontology annotation for anonymous sequence data

Gene Ontology (GO) is the most widely accepted attempt to construct a unified and structured vocabulary for the description of genes and their products in any organism. Annotation by GO terms is performed in most of the current genome projects, which besides generality has the advantage of being very convenient for computer based classification methods. However, direct use of GO in small sequencing projects is not easy, especially for species not commonly represented in public databases. We present a software package (GOblet), which performs annotation based on GO terms for anonymous cDNA or protein sequences. It uses the species independent GO structure and vocabulary together with a series of protein databases collected from various sites, to perform a detailed GO annotation by sequence similarity searches. The sensitivity and the reference protein sets can be selected by the user. GOblet runs automatically and is available as a public service on our web server. The paper also addresses the reliability of automated GO annotations by using a reference set of more than 6000 human proteins. The GOblet server is accessible at http://goblet.molgen.mpg.de.

Structural dynamics of eukaryotic chromosome evolution

Large-scale genome sequencing is providing a comprehensive view of the complex evolutionary forces that have shaped the structure of eukaryotic chromosomes. Comparative sequence analyses reveal patterns of apparently random rearrangement interspersed with regions of extraordinarily rapid, localized genome evolution. Numerous subtle rearrangements near centromeres, telomeres, duplications, and interspersed repeats suggest hotspots for eukaryotic chromosome evolution. This localized chromosomal instability may play a role in rapidly evolving lineage-specific gene families and in fostering large-scale changes in gene order. Computational algorithms that take into account these dynamic forces along with traditional models of chromosomal rearrangement show promise for reconstructing the natural history of eukaryotic chromosomes.

Characterization of the genomic organization of the region bordering the centromere of chromosome V of Podospora anserina by direct sequencing

A Podospora anserina BAC library of 4800 clones has been constructed in the vector pBHYG allowing direct selection in fungi. Screening of the BAC collection for centromeric sequences of chromosome V allowed the recovery of clones localized on either sides of the centromere, but no BAC clone was found to contain the centromere. Seven BAC clones containing 322,195 and 156,244bp from either sides of the centromeric region were sequenced and annotated. One 5S rRNA gene, 5 tRNA genes, and 163 putative coding sequences (CDS) were identified. Among these, only six CDS seem specific to P. anserina. The gene density in the centromeric region is approximately one gene every 2.8kb. Extrapolation of this gene density to the whole genome of P. anserina suggests that the genome contains about 11,000 genes. Synteny analyses between P. anserina and Neurospora crassa show that co-linearity extends at the most to a few genes, suggesting rapid genome rearrangements between these two species.

The Feline Genome Project

The compilation of a dense gene map and eventually a whole genome sequence (WGS) of the domestic cat holds considerable value for human genome annotation, for veterinary medicine, and for insight into the evolution of genome organization among mammals. Human association and veterinary studies of the cat, its domestic breeds, and its charismatic wild relatives of the family Felidae have rendered the species a powerful model for human hereditary diseases, for infectious disease agents, for adaptive evolutionary divergence, for conservation genetics, and for forensic applications. Here we review the advantages, rationale, and present strategy of a feline genome project, and we describe the disease models, comparative genomics, and biological applications posed by the full resolution of the cat's genome.

Examination of sequence homology between human chromosome 20 and the mouse genome: intense conservation of many genomic elements

The conservation of genomic organization of mammalian species has been of interest for its usefulness in characterizing the genetics of traits and diseases and as one tool for examining evolution. The recent rough draft sequencing of the mouse and human genomes provides the opportunity for more detailed analyses. The current study examines the extent of homology between human chromosome 20 and the mouse genome by comparing putative coding and non-coding sequence to provide insight into organizational and sequence similarities between the species. The relative position of each of 460 putative coding orthologues was the same in both species, except for a single genomic segment rearrangement. The similarity extended to exon/intron structure, the size of introns, as well as strong evidence for the conservation of position of ancient LINE-1, LINE-2 and LTR repetitive sequence and the subtelomeric region of the long arm of human chromosome 20 and that of mouse chromosome 2. There was also evidence for conservation of a limited amount of non-coding single-copy sequence. Together these data provide additional insight into the extent of conservation of mammalian genomic organization and sequence.

Assessing the Drosophila melanogaster and Anopheles gambiae genome annotations using genome-wide sequence comparisons

We performed genome-wide sequence comparisons at the protein coding level between the genome sequences of Drosophila melanogaster and Anopheles gambiae. Such comparisons detect evolutionarily conserved regions (ecores) that can be used for a qualitative and quantitative evaluation of the available annotations of both genomes. They also provide novel candidate features for annotation. The percentage of ecores mapping outside annotations in the A. gambiae genome is about fourfold higher than in D. melanogaster. The A. gambiae genome assembly also contains a high proportion of duplicated ecores, possibly resulting from artefactual sequence duplications in the genome assembly. The occurrence of 4063 ecores in the D. melanogaster genome outside annotations suggests that some genes are not yet or only partially annotated. The present work illustrates the power of comparative genomics approaches towards an exhaustive and accurate establishment of gene models and gene catalogues in insect genomes.

SLAM web server for comparative gene finding and alignment

SLAM is a program that simultaneously aligns and annotates pairs of homologous sequences. The SLAM web server integrates SLAM with repeat masking tools and the AVID alignment program to allow for rapid alignment and gene prediction in user submitted sequences. Along with annotations and alignments for the submitted sequences, users obtain a list of predicted conserved non-coding sequences (and their associated alignments). The web site also links to whole genome annotations of the human, mouse and rat genomes produced with the SLAM program. The server can be accessed at http://bio.math.berkeley.edu/slam.

Altered astrocyte calcium homeostasis and proliferation in theTs65Dn mouse, a model of Down syndrome

Genes from the Down syndrome (DS) critical region of human chromosome 21, which contribute to the pathology of DS, are also found on mouse chromosome 16. Several animal models of DS with triplication of genes from the DS critical region have been generated, including mouse trisomy 16 (Ts16) and a partial trisomic mouse, Ts65Dn. Using computer-assisted imaging of fura-2 fluorescence, we found an elevation of intracellular cytoplasmic calcium in cortical astrocytes from neonatal Ts65Dn mouse brain, similar to that observed previously in embryonic Ts16 astrocytes. Furthermore, astrocytes from both Ts65Dn and Ts16 cortex fail to respond to the anti-proliferative actions of glutamate. These results suggest that defective regulation of cell proliferation and cellular calcium can result from triplication of DS critical region genes.

Human and mouse genomic sequences reveal extensive breakpoint reuse in mammalian evolution

The human and mouse genomic sequences provide evidence for a larger number of rearrangements than previously thought and reveal extensive reuse of breakpoints from the same short fragile regions. Breakpoint clustering in regions implicated in cancer and infertility have been reported in previous studies; we report here on breakpoint clustering in chromosome evolution. This clustering reveals limitations of the widely accepted random breakage theory that has remained unchallenged since the mid-1980s. The genome rearrangement analysis of the human and mouse genomes implies the existence of a large number of very short "hidden" synteny blocks that were invisible in the comparative mapping data and ignored in the random breakage model. These blocks are defined by closely located breakpoints and are often hard to detect. Our results suggest a model of chromosome evolution that postulates that mammalian genomes are mosaics of fragile regions with high propensity for rearrangements and solid regions with low propensity for rearrangements.

Inhibition of tumor growth with a vaccine based on xenogeneic homologous fibroblast growth factor receptor-1 in mice

Angiogenesis is important for the growth of solid tumors. The breaking of the immune tolerance against the molecule associated with angiogenesis should be a useful approach for cancer therapy. However, the immunity to self-molecules is difficult to elicit by a vaccine based on autologous or syngeneic molecules due to immune tolerance. Basic fibroblast growth factor (bFGF) is a specific and potent angiogenic factor implicated in tumor growth. The biological activity of bFGF is mediated through interaction with its high-affinity receptor, fibroblast growth factor receptor-1 (FGFR-1). In this study, we selected Xenopus FGFR-1 as a model antigen by the breaking of immune tolerance to explore the feasibility of cancer therapy in murine tumor models. We show here that vaccination with Xenopus FGFR-1 (pxFR1) is effective at antitumor immunity in three murine models. FGFR-1-specific autoantibodies in sera of pxFR1-immunized mice could be found in Western blotting analysis. The purified immunoglobulins were effective at the inhibition of endothelial cell proliferation in vitro and at the antitumor activity in vivo. The antitumor activity and production of FGFR-1-specific autoantibodies could be abrogated by depletion of CD4+ T lymphocytes. Histological examination revealed that the autoantibody was deposited on the endothelial cells within tumor tissues from pxFR1-immunized mice, and intratumoral angiogenesis was significantly suppressed. Furthermore, the inhibition of angiogenesis could also be found in alginate-encapsulate tumor cell assay. These observations may provide a new vaccine strategy for cancer therapy through the induction of autoimmunity against FGFR-1 associated with angiogenesis in a cross-reaction.

Human-mouse gene identification by comparative evidence integration and evolutionary analysis

The identification of genes in the human genome remains a challenge, as the actual predictions appear to disagree tremendously and vary dramatically on the basis of the specific gene-finding methodology used. Because the pattern of conservation in coding regions is expected to be different from intronic or intergenic regions, a comparative computational analysis can lead, in principle, to an improved computational identification of genes in the human genome by using a reference, such as mouse genome. However, this comparative methodology critically depends on three important factors: (1) the selection of the most appropriate reference genome. In particular, it is not clear whether the mouse is at the correct evolutionary distance from the human to provide sufficiently distinctive conservation levels in different genomic regions, (2) the selection of comparative features that provide the most benefit to gene recognition, and (3) the selection of evidence integration architecture that effectively interprets the comparative features. We address the first question by a novel evolutionary analysis that allows us to explicitly correlate the performance of the gene recognition system with the evolutionary distance (time) between the two genomes. Our simulation results indicate that there is a wide range of reference genomes at different evolutionary time points that appear to deliver reasonable comparative prediction of human genes. In particular, the evolutionary time between human and mouse generally falls in the region of good performance; however, better accuracy might be achieved with a reference genome further than mouse. To address the second question, we propose several natural comparative measures of conservation for identifying exons and exon boundaries. Finally, we experiment with Bayesian networks for the integration of comparative and compositional evidence.

G protein-coupled receptor genes in the FANTOM2 database

G protein-coupled receptors (GPCRs) comprise the largest family of receptor proteins in mammals and play important roles in many physiological and pathological processes. Gene expression of GPCRs is temporally and spatially regulated, and many splicing variants are also described. In many instances, different expression profiles of GPCR gene are accountable for the changes of its biological function. Therefore, it is intriguing to assess the complexity of the transcriptome of GPCRs in various mammalian organs. In this study, we took advantage of the FANTOM2 (Functional Annotation Meeting of Mouse cDNA 2) project, which aimed to collect full-length cDNAs inclusively from mouse tissues, and found 410 candidate GPCR cDNAs. Clustering of these clones into transcriptional units (TUs) reduced this number to 213. Out of these, 165 genes were represented within the known 308 GPCRs in the Mouse Genome Informatics (MGI) resource. The remaining 48 genes were new to mouse, and 14 of them had no clear mammalian ortholog. To dissect the detailed characteristics of each transcript, tissue distribution pattern and alternative splicing were also ascertained. We found many splicing variants of GPCRs that may have a relevance to disease occurrence. In addition, the difficulty in cloning tissue-specific and infrequently transcribed GPCRs is discussed further.

Mouse proteome analysis

A general overview of the protein sequence set for the mouse transcriptome produced during the FANTOM2 sequencing project is presented here. We applied different algorithms to characterize protein sequences derived from a nonredundant representative protein set (RPS) and a variant protein set (VPS) of the mouse transcriptome. The functional characterization and assignment of Gene Ontology terms was done by analysis of the proteome using InterPro. The Superfamily database analyses gave a detailed structural classification according to SCOP and provide additional evidence for the functional characterization of the proteome data. The MDS database analysis revealed new domains which are not presented in existing protein domain databases. Thus the transcriptome gives us a unique source of data for the detection of new functional groups. The data obtained for the RPS and VPS sets facilitated the comparison of different patterns of protein expression. A comparison of other existing mouse and human protein sequence sets (e.g., the International Protein Index) demonstrates the common patterns in mammalian proteomes. The analysis of the membrane organization within the transcriptome of multiple eukaryotes provides valuable statistics about the distribution of secretory and transmembrane proteins

The SCAN domain defines a large family of zinc finger transcription factors

The SCAN domain is a highly conserved dimerization motif that is vertebrate-specific and found near the N-terminus of C(2)H(2) zinc finger proteins (SCAN-ZFP). Although the function of most SCAN-ZFPs is unknown, some have been implicated in the transcriptional regulation of growth factors, genes involved in lipid metabolism, as well as other genes involved in cell survival and differentiation. Here we utilize a bioinformatics approach to define the structures and gene locations of the 71 members of the human SCAN domain family, as well as to assess the conserved syntenic segments in the mouse genome and identify potential orthologs. The genes encoding SCAN domains are clustered, often in tandem arrays, in both the human and mouse genomes and are capable of generating isoforms that may affect the function of family members. Twenty-three members of the mouse SCAN family appear to be orthologous with human family members, and human-specific cluster expansions were observed. Remarkably, the SCAN domains in lower vertebrates are not associated with C(2)H(2) zinc finger genes, but are contained in large retrovirus-like polyproteins. Collectively, these studies define a large family of vertebrate-specific transcriptional regulators that may have rapidly expanded during recent evolution.

Conservation of human alternative splice events in mouse

Human and mouse genomes share similar long-range sequence organization, and have most of their genes being homologous. As alternative splicing is a frequent and important aspect of gene regulation, it is of interest to assess the level of conservation of alternative splicing. We examined mouse transcript data sets (EST and mRNA) for the presence of transcripts that both make spliced-alignment with the draft mouse genome sequence and demonstrate conservation of human transcript-confirmed alternative and constitutive splice junctions. This revealed 15% of alternative and 67% of constitutive splice junctions as conserved; however, these numbers are patently dependent on the extent of transcript coverage. Transcript coverage of conserved splice patterns is found to correlate well between human and mouse. A model, which extrapolates from observed levels of conservation at increasing levels of transcript support, estimates overall conservation of 61% of alternative and 74% of constitutive splice junctions, albeit with broad confidence intervals. Observed numbers of conserved alternative splicing events agreed with those expected on the basis of the model. Thus, it is apparent that many, and probably most, alternative splicing events are conserved between human and mouse. This, combined with the preservation of alternative frame stop codons in conserved frame breaking events, indicates a high level of commonality in patterns of gene expression between these two species.

Human chromosome 7: DNA sequence and biology

DNA sequence and annotation of the entire human chromosome 7, encompassing nearly 158 million nucleotides of DNA and 1917 gene structures, are presented. To generate a higher order description, additional structural features such as imprinted genes, fragile sites, and segmental duplications were integrated at the level of the DNA sequence with medical genetic data, including 440 chromosome rearrangement breakpoints associated with disease. This approach enabled the discovery of candidate genes for developmental diseases including autism.

Quantitative estimates of sequence divergence for comparative analyses of mammalian genomes

Comparative sequence analyses on a collection of carefully chosen mammalian genomes could facilitate identification of functional elements within the human genome and allow quantification of evolutionary constraint at the single nucleotide level. High-resolution quantification would be informative for determining the distribution of important positions within functional elements and for evaluating the relative importance of nucleotide sites that carry single nucleotide polymorphisms (SNPs). Because the level of resolution in comparative sequence analyses is a direct function of sequence diversity, we propose that the information content of a candidate mammalian genome be defined as the sequence divergence it would add relative to already-sequenced genomes. We show that reliable estimates of genomic sequence divergence can be obtained from small genomic regions. On the basis of a multiple sequence alignment of approximately 1.4 megabases each from eight mammals, we generate such estimates for five unsequenced mammals. Estimates of the neutral divergence in these data suggest that a small number of diverse mammalian genomes in addition to human, mouse, and rat would allow single nucleotide resolution in comparative sequence analyses.

Structure and chromosomal localization of the murine coxsackievirus and adenovirus receptor gene

We analyzed BAC genomic clones encoding the murine coxsackievirus and adenovirus receptor (mCAR). The mCAR gene is situated on the distal portion of murine chromosome 16, and is composed of at least eight exons, with intron-exon boundaries similar to those reported for the human CAR gene. We previously described two cDNAs encoding mCAR isoforms: the extracellular and transmembrane portions of both are encoded by exons 1-6; the cytoplasmic domain of mCAR 1 is encoded by exon 7, whereas mCAR 2 results from an RNA splice linking the proximal portion of exon 7 to an alternative exon 8. RT-PCR analysis of the mCAR RNA 5'-terminus suggests that transcription may begin 141-161 nucleotides upstream of the ATG translational start site.

Comparative genomics: genome-wide analysis in metazoan eukaryotes

The increasing number of complete and nearly complete metazoan genome sequences provides a significant amount of material for large-scale comparative genomic analysis. Finding new effective methods to analyse such enormous datasets has been the object of intense research. Three main areas in comparative genomics have recently shown important developments: whole-genome alignment, gene prediction and regulatory-region prediction. Each of these areas improves the methods of deciphering long genomic sequences and uncovering what lies hidden in them.

Genome-wide single-nucleotide polymorphism analysis defines haplotype patterns in mouse

The nature and organization of polymorphisms, or differences, between genomes of individuals are of great interest, because these variations can be associated with or even underlie phenotypic traits, including disease susceptibility. To gain insight into the genetic and evolutionary factors influencing such biological variation, we have examined the arrangement (haplotype) of single-nucleotide polymorphisms across the genomes of eight inbred strains of mice. These analyses define blocks of high or low diversity, often extending across tens of megabases that are delineated by abrupt transitions. These observations provide a striking contrast to the haplotype structure of the human genome.

Molecular evolution meets the genomics revolution

Changes in technology in the past decade have had such an impact on the way that molecular evolution research is done that it is difficult now to imagine working in a world without genomics or the Internet. In 1992, GenBank was less than a hundredth of its current size and was updated every three months on a huge spool of tape. Homology searches took 30 minutes and rarely found a hit. Now it is difficult to find sequences with only a few homologs to use as examples for teaching bioinformatics. For molecular evolution researchers, the genomics revolution has showered us with raw data and the information revolution has given us the wherewithal to analyze it. In broad terms, the most significant outcome from these changes has been our newfound ability to examine the evolution of genomes as a whole, enabling us to infer genome-wide evolutionary patterns and to identify subsets of genes whose evolution has been in some way atypical.

A review and comparison of the murine alpha1-antitrypsin and alpha1-antichymotrypsin multigene clusters with the human clade A serpins

The major human plasma protease inhibitors, alpha(1)-antitrypsin and alpha(1)-antichymotrypsin, are each encoded by a single gene, whereas in the mouse they are represented by clusters of 5 and 14 genes, respectively. Although there is a high degree of overall sequence similarity within these groupings, the reactive-center loop (RCL) domain, which determines target protease specificity, is markedly divergent. The literature dealing with members of these mouse serine protease inhibitor (serpin) clusters has been complicated by inconsistent nomenclature. Furthermore, some investigators, unaware of the complexity of the family, have failed to distinguish between closely related genes when measuring expression levels or functional activity. We have reviewed the literature dealing with the mouse equivalents of human alpha(1)-antitrypsin and alpha(1)-antichymotrypsin and made use of the recently completed mouse genome sequence to propose a systematic nomenclature. We have also examined the extended mouse clade "a" serpin cluster at chromosome 12F1 and compared it with the syntenic region at human chromosome 14q32. In summarizing the literature and suggesting a standardized nomenclature, we aim to provide a logical structure on which future research may be based.

Analysis of primate genomic variation reveals a repeat-driven expansion of the human genome

We performed a detailed analysis of both single-nucleotide and large insertion/deletion events based on large-scale comparison of 10.6 Mb of genomic sequence from lemur, baboon, and chimpanzee to human. Using a human genomic reference, optimal global alignments were constructed from large (>50-kb) genomic sequence clones. These alignments were examined for the pattern, frequency, and nature of mutational events. Whereas rates of single-nucleotide substitution remain relatively constant (1-2 x 10(-9) substitutions/site/year), rates of retrotransposition vary radically among different primate lineages. These differences have lead to a 15%-20% expansion of human genome size over the last 50 million years of primate evolution, 90% of it due to new retroposon insertions. Orthologous comparisons with the chimpanzee suggest that the human genome continues to significantly expand due to shifts in retrotransposition activity. Assuming that the primate genome sequence we have sampled is representative, we estimate that human euchromatin has expanded 30 Mb and 550 Mb compared to the primate genomes of chimpanzee and lemur, respectively.

Massive parallelism, randomness and genomic advances

In reviewing the past decade, it is clear that genomics was, and still is, driven by innovative technologies, perhaps more so than any other scientific area in recent memory. From the outset, computing, mathematics and new automated laboratory techniques have been key components in allowing the field to move forward rapidly. We highlight some key innovations that have come together to nurture the explosive growth that makes a new era of genomics a reality. We also document how these new approaches have fueled further innovations and discoveries.

Novel quantitative trait loci controlling development of experimental autoimmune encephalomyelitis and proportion of lymphocyte subpopulations

The B10.RIII mouse strain (H-2(r)) develops chronic experimental autoimmune encephalomyelitis (EAE) upon immunization with the myelin basic protein 89-101 peptide. EAE was induced and studied in a backcross between B10.RIII and the EAE-resistant RIIIS/J strain (H-2(r)), and a complete genome scan with microsatellite markers was performed. Five loci were significantly linked to different traits and clinical subtypes of EAE on chromosomes 1, 5, 11, 15, and 16, three of the loci having sex specificity. The quantitative trait locus on chromosome 15 partly overlapped with the Eae2 locus, previously identified in crosses between the B10.RIII and RIIIS/J mouse strains. The loci on chromosomes 11 and 16 overlapped with Eae loci identified in other mouse crosses. By analyzing the backcross animals for lymphocyte phenotypes, the proportion of B and T cells in addition to the levels of CD4(+)CD8(-) and CD4(-)CD8(+) T cells and the CD4(+)/CD8(+) ratio in spleen were linked to different loci on chromosomes 1, 2, 3, 5, 6, 11, and 15. On chromosome 16, we found significant linkage to spleen cell proliferation. Several linkages overlapped with the quantitative trait loci for disease phenotypes. The identification of subphenotypes that are linked to the same loci as disease traits could be most useful in the search for candidate genes and biological pathways involved in the pathological process.

A novel class of ligand-gated ion channel is activated by Zn2+

In mammals, the superfamily of "Cys loop," ligand-gated ion channels (LGICs), is assembled from a pool of more than 40 homologous subunits. These subunits have been classified into four families representing channels that are gated by acetylcholine, serotonin, gamma-aminobutyric acid, or glycine. By searching anonymous genomic sequence data for exons that encode characteristic motifs of the channel subunits, we have identified a novel LGIC that defines a fifth family member. Putative exons were used to isolate a full-length cDNA that encodes a protein of 411 amino acid residues. This protein (ZAC) contains all of the motifs that are characteristic of Cys loop channel subunits but cannot be assigned to any of the four established families on the basis of sequence similarity. Genes for ZAC are present in human and dog but appear to have been lost from mouse and rat genomes. Transcripts of ZAC subunits were detected in human placenta, trachea, spinal cord, stomach, and fetal brain. Transfection of human embryonic kidney cells with ZAC subunit cDNA caused expression of spontaneous current. By screening with a broad range of potential agonists and antagonists, we determined that tubocurarine inhibits the spontaneous current whereas Zn(2+) activates the expressed receptors. The absence of Zn(2+)-activated channels in rats and mice may explain why this fifth member of the LGIC superfamily has evaded detection until now.

Animal models of mental retardation: from gene to cognitive function

About 2-3% of all children are affected by mental retardation, and genetic conditions rank among the leading causes of mental retardation. Alterations in the information encoded by genes that regulate critical steps of brain development can disrupt the normal course of development, and have profound consequences on mental processes. Genetically modified mouse models have helped to elucidate the contribution of specific gene alterations and gene-environment interactions to the phenotype of several forms of mental retardation. Mouse models of several neurodevelopmental pathologies, such as Down and Rett syndromes and X-linked forms of mental retardation, have been developed. Because behavior is the ultimate output of brain, behavioral phenotyping of these models provides functional information that may not be detectable using molecular, cellular or histological evaluations. In particular, the study of ontogeny of behavior is recommended in mouse models of disorders having a developmental onset. Identifying the role of specific genes in neuropathologies provides a framework in which to understand key stages of human brain development, and provides a target for potential therapeutic intervention.

BayGenomics: a resource of insertional mutations in mouse embryonic stem cells

The BayGenomics gene-trap resource (http://baygenomics.ucsf.edu) provides researchers with access to thousands of mouse embryonic stem (ES) cell lines harboring characterized insertional mutations in both known and novel genes. Each cell line contains an insertional mutation in a specific gene. The identity of the gene that has been interrupted can be determined from a DNA sequence tag. Approximately 75% of our cell lines contain insertional mutations in known mouse genes or genes that share strong sequence similarities with genes that have been identified in other organisms. These cell lines readily transmit the mutation to the germline of mice and many mutant lines of mice have already been generated from this resource. BayGenomics provides facile access to our entire database, including sequence tags for each mutant ES cell line, through the World Wide Web. Investigators can browse our resource, search for specific entries, download any portion of our database and BLAST sequences of interest against our entire set of cell line sequence tags. They can then obtain the mutant ES cell line for the purpose of generating knockout mice.

Whole-genome sequence assembly for mammalian genomes: Arachne 2

We previously described the whole-genome assembly program Arachne, presenting assemblies of simulated data for small to mid-sized genomes. Here we describe algorithmic adaptations to the program, allowing for assembly of mammalian-size genomes, and also improving the assembly of smaller genomes. Three principal changes were simultaneously made and applied to the assembly of the mouse genome, during a six-month period of development: (1) Supercontigs (scaffolds) were iteratively broken and rejoined using several criteria, yielding a 64-fold increase in length (N50), and apparent elimination of all global misjoins; (2) gaps between contigs in supercontigs were filled (partially or completely) by insertion of reads, as suggested by pairing within the supercontig, increasing the N50 contig length by 50%; (3) memory usage was reduced fourfold. The outcome of this mouse assembly and its analysis are described in (Mouse Genome Sequencing Consortium 2002).

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.

Knockouts model the 100 best-selling drugs--will they model the next 100?

The biopharmaceutical industry is currently faced with a tremendous number of potential drug targets identified through the sequencing of the human genome. The challenge ahead is to delineate those targets with the greatest value for therapeutic intervention. Here, we critically evaluate mouse-knockout technology for target discovery and validation. A retrospective evaluation of the knockout phenotypes for the targets of the 100 best-selling drugs indicates that these phenotypes correlate well with known drug efficacy, illuminating a productive path forward for discovering future drug targets. Prospective mining of the druggable genome is being catalysed by large-scale mouse knockout programs combined with phenotypic screens focused on identifying targets that modulate mammalian physiology in a therapeutically relevant manner.

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.

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.

Initial sequencing and comparative analysis of the mouse genome

The sequence of the mouse genome is a key informational tool for understanding the contents of the human genome and a key experimental tool for biomedical research. Here, we report the results of an international collaboration to produce a high-quality draft sequence of the mouse genome. We also present an initial comparative analysis of the mouse and human genomes, describing some of the insights that can be gleaned from the two sequences. We discuss topics including the analysis of the evolutionary forces shaping the size, structure and sequence of the genomes; the conservation of large-scale synteny across most of the genomes; the much lower extent of sequence orthology covering less than half of the genomes; the proportions of the genomes under selection; the number of protein-coding genes; the expansion of gene families related to reproduction and immunity; the evolution of proteins; and the identification of intraspecies polymorphism.

Analysis of the mouse transcriptome based on functional annotation of 60,770 full-length cDNAs

Only a small proportion of the mouse genome is transcribed into mature messenger RNA transcripts. There is an international collaborative effort to identify all full-length mRNA transcripts from the mouse, and to ensure that each is represented in a physical collection of clones. Here we report the manual annotation of 60,770 full-length mouse complementary DNA sequences. These are clustered into 33,409 'transcriptional units', contributing 90.1% of a newly established mouse transcriptome database. Of these transcriptional units, 4,258 are new protein-coding and 11,665 are new non-coding messages, indicating that non-coding RNA is a major component of the transcriptome. 41% of all transcriptional units showed evidence of alternative splicing. In protein-coding transcripts, 79% of splice variations altered the protein product. Whole-transcriptome analyses resulted in the identification of 2,431 sense-antisense pairs. The present work, completely supported by physical clones, provides the most comprehensive survey of a mammalian transcriptome so far, and is a valuable resource for functional genomics.

Numerous potentially functional but non-genic conserved sequences on human chromosome 21

The use of comparative genomics to infer genome function relies on the understanding of how different components of the genome change over evolutionary time. The aim of such comparative analysis is to identify conserved, functionally transcribed sequences such as protein-coding genes and non-coding RNA genes, and other functional sequences such as regulatory regions, as well as other genomic features. Here, we have compared the entire human chromosome 21 with syntenic regions of the mouse genome, and have identified a large number of conserved blocks of unknown function. Although previous studies have made similar observations, it is unknown whether these conserved sequences are genes or not. Here we present an extensive experimental and computational analysis of human chromosome 21 in an effort to assign function to sequences conserved between human chromosome 21 (ref. 8) and the syntenic mouse regions. Our data support the presence of a large number of potentially functional non-genic sequences, probably regulatory and structural. The integration of the properties of the conserved components of human chromosome 21 to the rapidly accumulating functional data for this chromosome will improve considerably our understanding of the role of sequence conservation in mammalian genomes.

In silico discovery of gene-coding variants in murine quantitative trait loci using strain-specific genome sequence databases

BACKGROUND

The identification of genes underlying complex traits has been aided by quantitative trait locus (QTL) mapping approaches, which in turn have benefited from advances in mammalian genome research. Most recently, whole-genome draft sequences and assemblies have been generated for mouse strains that have been used for a large fraction of QTL mapping studies. Here we show how such strain-specific mouse genome sequence databases can be used as part of a high-throughput pipeline for the in silico discovery of gene-coding variations within murine QTLs. As a test of this approach we focused on two QTLs on mouse chromosomes 1 and 13 that are involved in physical dependence on alcohol.

RESULTS

Interstrain alignment of sequences derived from the relevant mouse strain genome sequence databases for 199 QTL-localized genes spanning 210,020 base-pairs of coding sequence identified 21 genes with different coding sequences for the progenitor strains. Several of these genes, including four that exhibit strong phenotypic links to chronic alcohol withdrawal, are promising candidates to underlie these QTLs.

CONCLUSIONS

This approach has wide general utility, and should be applicable to any of the several hundred mouse QTLs, encompassing over 60 different complex traits, that have been identified using strains for which relatively complete genome sequences are available.

Characterization of transcription factors by mass spectrometry and the role of SELDI-MS

Over the last decade, much progress has been made in the field of biological mass spectrometry, with numerous advances in technology, resolution, and affinity capture. The field of genomics has also been transformed by the sequencing and characterization of entire genomes. Some of the next challenges lie in understanding the relationship between the genome and the proteome, the protein complement of the genome, and in characterizing the regulatory processes involved in progressing from gene to functional protein. In this new age of proteomics, development of mass spectrometry methods to characterize transcription factors promises to add greatly to our understanding of regulatory networks that govern expression. However, at this time, regulatory networks of transcription factors are mostly uncharted territory. In this review, we summarize the latest advances in characterization of transcription factors by mass spectrometry including affinity capture, identification of complexes of DNA-binding proteins, structural characterization, determination of protein-DNA and protein-protein interactions, assessment of modification sites and metal binding, studies of functional activity, and the latest chip technologies that use SELDI-MS that allow the rapid capture and identification of transcription factors.

The genome sequence of the malaria mosquito Anopheles gambiae

Anopheles gambiae is the principal vector of malaria, a disease that afflicts more than 500 million people and causes more than 1 million deaths each year. Tenfold shotgun sequence coverage was obtained from the PEST strain of A. gambiae and assembled into scaffolds that span 278 million base pairs. A total of 91% of the genome was organized in 303 scaffolds; the largest scaffold was 23.1 million base pairs. There was substantial genetic variation within this strain, and the apparent existence of two haplotypes of approximately equal frequency ("dual haplotypes") in a substantial fraction of the genome likely reflects the outbred nature of the PEST strain. The sequence produced a conservative inference of more than 400,000 single-nucleotide polymorphisms that showed a markedly bimodal density distribution. Analysis of the genome sequence revealed strong evidence for about 14,000 protein-encoding transcripts. Prominent expansions in specific families of proteins likely involved in cell adhesion and immunity were noted. An expressed sequence tag analysis of genes regulated by blood feeding provided insights into the physiological adaptations of a hematophagous insect.

Genome sequence and comparative analysis of the model rodent malaria parasite Plasmodium yoelii yoelii

Species of malaria parasite that infect rodents have long been used as models for malaria disease research. Here we report the whole-genome shotgun sequence of one species, Plasmodium yoelii yoelii, and comparative studies with the genome of the human malaria parasite Plasmodium falciparum clone 3D7. A synteny map of 2,212 P. y. yoelii contiguous DNA sequences (contigs) aligned to 14 P. falciparum chromosomes reveals marked conservation of gene synteny within the body of each chromosome. Of about 5,300 P. falciparum genes, more than 3,300 P. y. yoelii orthologues of predominantly metabolic function were identified. Over 800 copies of a variant antigen gene located in subtelomeric regions were found. This is the first genome sequence of a model eukaryotic parasite, and it provides insight into the use of such systems in the modelling of Plasmodium biology and disease.