Ensembl 2007.

The Ensembl (http://www.ensembl.org/) project provides a comprehensive and integrated source of annotation of chordate genome sequences. Over the past year the number of genomes available from Ensembl has increased from 15 to 33, with the addition of sites for the mammalian genomes of elephant, rabbit, armadillo, tenrec, platypus, pig, cat, bush baby, common shrew, microbat and european hedgehog; the fish genomes of stickleback and medaka and the second example of the genomes of the sea squirt (Ciona savignyi) and the mosquito (Aedes aegypti). Some of the major features added during the year include the first complete gene sets for genomes with low-sequence coverage, the introduction of new strain variation data and the introduction of new orthology/paralog annotations based on gene trees.

Ensembl 2006.

The Ensembl (http://www.ensembl.org/) project provides a comprehensive and integrated source of annotation of large genome sequences. Over the last year the number of genomes available from the Ensembl site has increased from 4 to 19, with the addition of the mammalian genomes of Rhesus macaque and Opossum, the chordate genome of Ciona intestinalis and the import and integration of the yeast genome. The year has also seen extensive improvements to both data analysis and presentation, with the introduction of a redesigned website, the addition of RNA gene and regulatory annotation and substantial improvements to the integration of human genome variation data.

Ensembl 2005.

The Ensembl (http://www.ensembl.org/) project provides a comprehensive and integrated source of annotation of large genome sequences. Over the last year the number of genomes available from the Ensembl site has increased by 7 to 16, with the addition of the six vertebrate genomes of chimpanzee, dog, cow, chicken, tetraodon and frog and the insect genome of honeybee. The majority have been annotated automatically using the Ensembl gene build system, showing its flexibility to reliably annotate a wide variety of genomes. With the increased number of vertebrate genomes, the comparative analysis provided to users has been greatly improved, with new website interfaces allowing annotation of different genomes to be directly compared. The Ensembl software system is being increasingly widely reused in different projects showing the benefits of a completely open approach to software development and distribution.

Improved techniques for the identification of pseudogenes.

MOTIVATION: Pseudogenes are the remnants of genomic sequences of genes which are no longer functional. They are frequent in most eukaryotic genomes, and an important resource for comparative genomics. However, pseudogenes are often mis-annotated as functional genes in sequence databases. Current methods for identifying pseudogenes include methods which rely on the presence of stop codons and frameshifts, as well as methods based on the ratio of non-silent to silent nucleotide substitution rates (dN/dS). A recent survey concluded that 50% of human pseudogenes have no detectable truncation in their pseudo-coding regions, indicating that the former methods lack sensitivity. The latter methods have been used to find sets of genes enriched for pseudogenes, but are not specific enough to accurately separate pseudogenes from expressed genes. RESULTS: We introduce a program called pseudogene inference from loss of constraint (PSILC) which incorporates novel methods for separating pseudogenes from functional genes. The methods calculate the log-odds score that evolution along the final branch of the gene tree to the query gene has been according to the following constraints: A neutral nucleotide model compared to a Pfam domain encoding model (PSILC(nuc/dom)); A protein coding model compared to a Pfam domain encoding model (PSILC(prot/dom)). Using the manual annotation of human chromosome 6, we show that both these methods result in a more accurate classification of pseudogenes than dN/dS when a Pfam domain alignment is available. AVAILABILITY: PSILC is available from http://www.sanger.ac.uk/Software/PSILC

Thiazolidinedione therapy in the prevention/delay of type 2 diabetes in patients with impaired glucose tolerance and insulin resistance.

AIM: The second-generation thiazolidinediones (TZDs), rosiglitazone and pioglitazone, significantly decrease fasting plasma glucose and glycosylated haemoglobin (HbA(1c)) levels in patients with diabetes. Recent studies suggest that early treatment with TZDs may prevent the progression from insulin resistance (IR) to type 2 diabetes mellitus (T2DM). This prospective analysis examined the effect of early TZD treatment in the prevention or delay of T2DM in a multiethnic population with impaired glucose tolerance (IGT) and IR. METHODS: The analysis included 172 patients (aged 29-86 years) with IGT and IR (normal or borderline HbA(1c), C-peptide levels > 2 mg/ml, fasting blood sugar 100-125 mg/dl, and 2-h postprandial blood glucose levels 140-200 mg/dl). Patients in the active treatment group (n = 101) had received troglitazone for an average of 10 months before being randomly switched to rosiglitazone (4 mg/day) or pioglitazone (30 mg/day). Patients were switched when troglitazone was withdrawn from the US market because of liver toxicity concerns. Patients with IGT and IR who received no antidiabetic medication served as a control group (n = 71). HbA(1c) and C-peptide levels were measured at baseline (2 years) and study end point (3 years). Kaplan-Meier testing, using time to outcome as the main outcome variable, determined risk reduction in the TZD group relative to the control group. RESULTS: Mean HbA(1c) and C-peptide levels decreased for patients receiving either TZD at the 2-year assessment, and reductions were maintained at study end point. After 2 years, none of the patients receiving TZD therapy progressed to T2DM; three patients progressed to T2DM by study end point. In the control group, 11 patients became diabetic after 2 years and 19 patients became diabetic by the end of the study. The incidence (risk reduction) of diabetes after 3 years was 88.9% lower in the TZD group compared with the control group (p < 0.001). CONCLUSIONS: The TZDs, rosiglitazone and pioglitazone, were effective in reducing HbA(1c) and C-peptide levels in patients with IGT/IR. Progression of IR/IGT to T2DM appears to be significantly delayed or prevented with early TZD treatment.

Ensembl 2004.

The Ensembl (http://www.ensembl.org/) database project provides a bioinformatics framework to organize biology around the sequences of large genomes. It is a comprehensive and integrated source of annotation of large genome sequences, available via interactive website, web services or flat files. As well as being one of the leading sources of genome annotation, Ensembl is an open source software engineering project to develop a portable system able to handle very large genomes and associated requirements. The facilities of the system range from sequence analysis to data storage and visualization and installations exist around the world both in companies and at academic sites. With a total of nine genome sequences available from Ensembl and more genomes to follow, recent developments have focused mainly on closer integration between genomes and external data.

Ensembl 2002: accommodating comparative genomics.

The Ensembl (http://www.ensembl.org/) database project provides a bioinformatics framework to organise biology around the sequences of large genomes. It is a comprehensive source of stable automatic annotation of human, mouse and other genome sequences, available as either an interactive web site or as flat files. Ensembl also integrates manually annotated gene structures from external sources where available. As well as being one of the leading sources of genome annotation, Ensembl is an open source software engineering project to develop a portable system able to handle very large genomes and associated requirements. These range from sequence analysis to data storage and visualisation and installations exist around the world in both companies and at academic sites. With both human and mouse genome sequences available and more vertebrate sequences to follow, many of the recent developments in Ensembl have focusing on developing automatic comparative genome analysis and visualisation.

PKR protection against intranasal vesicular stomatitis virus infection is mouse strain dependent.

The interferon-induced antiviral state is mediated by interferon-stimulated genes that are upregulated in concert after stimulation by type I interferons. Because so many viruses encode strategies to inactivate the interferon-inducible double-stranded RNA (dsRNA)-dependent protein kinase PKR, this protein is likely to be a major player in antiviral defense. Here we demonstrate the increased susceptibility of PKR-/- animals to vesicular stomatitis virus (VSV) by the intranasal route, but also demonstrate that the protective effects of PKR are mouse strain dependent. We have found the difference between wild-type-BALB/c and 129SvEv animals to be on the order of 5 logs, with high levels of virus present in the lungs of BALB/c but not 129SvEv animals. To evaluate the sensitivity of PKR-/- mice to VSV clearly, the PKR mutation was bred onto the resistant 129SvEv background. The increased sensitivity of PKR-/- mice, compared to PKR+/+ strain-matched controls, is on the order of 10-fold as measured by median lethal dose (LD50). PKR-/- 129 mice support VSV replication in the lung unlike controls. While this result clearly demonstrates an important role for PKR in protection against VSV infection of the lung, it also underlines the importance of other host factors in containing a viral infection.

The Ensembl genome database project.

The Ensembl (http://www.ensembl.org/) database project provides a bioinformatics framework to organise biology around the sequences of large genomes. It is a comprehensive source of stable automatic annotation of the human genome sequence, with confirmed gene predictions that have been integrated with external data sources, and is available as either an interactive web site or as flat files. It is also an open source software engineering project to develop a portable system able to handle very large genomes and associated requirements from sequence analysis to data storage and visualisation. The Ensembl site is one of the leading sources of human genome sequence annotation and provided much of the analysis for publication by the international human genome project of the draft genome. The Ensembl system is being installed around the world in both companies and academic sites on machines ranging from supercomputers to laptops.

A computational scan for U12-dependent introns in the human genome sequence.

U12-dependent introns are found in small numbers in most eukaryotic genomes, but their scarcity makes accurate characterisation of their properties challenging. A computational search for U12-dependent introns was performed using the draft version of the human genome sequence. Human expressed sequences confirmed 404 U12-dependent introns within the human genome, a 6-fold increase over the total number of non-redundant U12-dependent introns previously identified in all genomes. Although most of these introns had AT-AC or GT-AG terminal dinucleotides, small numbers of introns with a surprising diversity of termini were found, suggesting that many of the non-canonical introns found in the human genome may be variants of U12-dependent introns and, thus, spliced by the minor spliceosome. Comparisons with U2-dependent introns revealed that the U12-dependent intron set lacks the 'short intron' peak characteristic of U2-dependent introns. Analysis of this U12-dependent intron set confirmed reports of a biased distribution of U12-dependent introns in the genome and allowed the identification of several alternative splicing events as well as a surprising number of apparent splicing errors. This new larger reference set of U12-dependent introns will serve as a resource for future studies of both the properties and evolution of the U12 spliceosome.

The SH integral membrane protein of the paramyxovirus simian virus 5 is required to block apoptosis in MDBK cells.

In some cell types the paramyxovirus simian virus 5 (SV5) causes little cytopathic effect (CPE) and infection continues productively for long periods of time; e.g., SV5 can be produced from MDBK cells for up to 40 days with little CPE. SV5 differs from most paramyxoviruses in that it encodes a small (44-amino-acid) hydrophobic integral membrane protein (SH). When MDBK cells were infected with a recombinant SV5 containing a deletion of the SH gene (rSV5DeltaSH), the MDBK cells exhibited an increase in CPE compared to cells infected with wild-type SV5 (recovered from cDNA; rSV5). The increased CPE correlated with an increase in apoptosis in rSV5DeltaSH-infected cells over mock-infected and rSV5-infected cells when assayed for annexin V binding, DNA content (propidium iodide staining), and DNA fragmentation (terminal deoxynucleotidyltransferase-mediated dUTP-biotin nick end labeling assay). In rSV5DeltaSH-infected MDBK cells an increase in caspase-2 and caspase-3 activities was observed. By using peptide inhibitors of individual caspases it was found that caspase-2 and caspase-3 were activated separately in rSV5DeltaSH-infected cells. Expression of caspase-2 and -3 in rSV5DeltaSH-infected MDBK cells appeared not to require STAT1 protein, as STAT1 protein could not be detected in SV5-infected MDBK cells. When mutant mice homologous for a targeted disruption of STAT1 were used as a model animal system and infected with the viruses it was found that rSV5DeltaSH caused less mortality than wild-type rSV5, consistent with the notion of clearance of apoptotic cells in a host species.

Initial sequencing and analysis of the human genome.

The human genome holds an extraordinary trove of information about human development, physiology, medicine and evolution. Here we report the results of an international collaboration to produce and make freely available a draft sequence of the human genome. We also present an initial analysis of the data, describing some of the insights that can be gleaned from the sequence.

The physical maps for sequencing human chromosomes 1, 6, 9, 10, 13, 20 and X.

We constructed maps for eight chromosomes (1, 6, 9, 10, 13, 20, X and (previously) 22), representing one-third of the genome, by building landmark maps, isolating bacterial clones and assembling contigs. By this approach, we could establish the long-range organization of the maps early in the project, and all contig extension, gap closure and problem-solving was simplified by containment within local regions. The maps currently represent more than 94% of the euchromatic (gene-containing) regions of these chromosomes in 176 contigs, and contain 96% of the chromosome-specific markers in the human gene map. By measuring the remaining gaps, we can assess chromosome length and coverage in sequenced clones.

The InterPro database, an integrated documentation resource for protein families, domains and functional sites.

Signature databases are vital tools for identifying distant relationships in novel sequences and hence for inferring protein function. InterPro is an integrated documentation resource for protein families, domains and functional sites, which amalgamates the efforts of the PROSITE, PRINTS, Pfam and ProDom database projects. Each InterPro entry includes a functional description, annotation, literature references and links back to the relevant member database(s). Release 2.0 of InterPro (October 2000) contains over 3000 entries, representing families, domains, repeats and sites of post-translational modification encoded by a total of 6804 different regular expressions, profiles, fingerprints and Hidden Markov Models. Each InterPro entry lists all the matches against SWISS-PROT and TrEMBL (more than 1,000,000 hits from 462,500 proteins in SWISS-PROT and TrEMBL). The database is accessible for text- and sequence-based searches at http://www.ebi.ac.uk/interpro/. Questions can be emailed to interhelp@ebi.ac.uk.

WormBase: network access to the genome and biology of Caenorhabditis elegans.

WormBase (http://www.wormbase.org) is a web-based resource for the Caenorhabditis elegans genome and its biology. It builds upon the existing ACeDB database of the C.elegans genome by providing data curation services, a significantly expanded range of subject areas and a user-friendly front end.

Epidemiology of Reye's syndrome, United States, 1991-1994: comparison of CDC surveillance and hospital admission data.

This investigation describes the epidemiology of Reye's syndrome (RS) during 1991-1994 and compares two different sources of information in the United States. Estimates of the incidence of RS from the Centers for Disease Control and Prevention (CDC) are compared with hospital inpatient data from approximately one third of the hospitals from HCIA, Inc. During 1991-1994, 48 RS cases were reported to the CDC and 93 RS hospitalizations based on HCIA data. When the HCIA data are extrapolated to the US population, there were an estimated 284 hospitalizations. Cases reported from both data sources were similar in distribution by onset, age, and sex. CDC data probably underestimate the incidence of RS due to incomplete reporting and HCIA data may overestimate it because not all cases were known to meet the CDC case definition. The true annual incidence of RS during the study years was probably between 0.2 and 1.1 cases per million population <18 years of age.

Alfresco--a workbench for comparative genomic sequence analysis.

Comparative analysis of genomic sequences provides a powerful tool for identifying regions of potential biologic function; by comparing corresponding regions of genomes from suitable species, protein coding or regulatory regions can be identified by their homology. This requires the use of several specific types of computational analysis tools. Many programs exist for these types of analysis; not many exist for overall view/control of the results, which is necessary for large-scale genomic sequence analysis. Using Java, we have developed a new visualization tool that allows effective comparative genome sequence analysis. The program handles a pair of sequences from putatively homologous regions in different species. Results from various different existing external analysis programs, such as database searching, gene prediction, repeat masking, and alignment programs, are visualized and used to find corresponding functional sequence domains in the two sequences. The user interacts with the program through a graphic display of the genome regions, in which an independently scrollable and zoomable symbolic representation of the sequences is shown. As an example, the analysis of two unannotated orthologous genomic sequences from human and mouse containing parts of the UTY locus is presented.

Using GeneWise in the Drosophila annotation experiment.

The GeneWise method for combining gene prediction and homology searches was applied to the 2.9-Mb region from Drosophila melanogaster. The results from the Genome Annotation Assessment Project (GASP) showed that GeneWise provided reasonably accurate gene predictions. Further investigation indicates that many of the incorrect gene predictions from GeneWise were due to transposons with valid protein-coding genes and the remaining cases are pseudogenes or possible annotation oversights.

InterPro--an integrated documentation resource for protein families, domains and functional sites.

MOTIVATION: InterPro is a new integrated documentation resource for protein families, domains and functional sites, developed initially as a means of rationalising the complementary efforts of the PROSITE, PRINTS, Pfam and ProDom database projects. RESULTS: Merged annotations from PRINTS, PROSITE and Pfam form the InterPro core. Each combined InterPro entry includes functional descriptions and literature references, and links are made back to the relevant parent database(s), allowing users to see at a glance whether a particular family or domain has associated patterns, profiles, fingerprints, etc. Merged and individual entries (i.e. those that have no counterpart in the companion resources) are assigned unique accession numbers. Release 1.2 of InterPro (June 2000) contains over 3000 entries, representing families, domains, repeats and sites of post-translational modification (PTMs) encoded by 6581 different regular expressions, profiles, fingerprints and Hidden Markov Models (HMMs). Each InterPro entry lists all the matches against SWISS-PROT and TrEMBL (more than 1000000 hits from 264333 different proteins out of 384572 in SWISS-PROT and TrEMBL).

The Pfam protein families database.

Pfam is a large collection of protein multiple sequence alignments and profile hidden Markov models. Pfam is available on the WWW in the UK at http://www.sanger.ac.uk/Software/Pfam/, in Sweden at http://www.cgr.ki.se/Pfam/ and in the US at http://pfam.wustl.edu/. The latest version (4.3) of Pfam contains 1815 families. These Pfam families match 63% of proteins in SWISS-PROT 37 and TrEMBL 9. For complete genomes Pfam currently matches up to half of the proteins. Genomic DNA can be directly searched against the Pfam library using the Wise2 package.