The COG database: new developments in phylogenetic classification of proteins from complete genomes

The database of Clusters of Orthologous Groups of proteins (COGs), which represents an attempt on a phylogenetic classification of the proteins encoded in complete genomes, currently consists of 2791 COGs including 45 350 proteins from 30 genomes of bacteria, archaea and the yeast Saccharomyces cerevisiae (http://www.ncbi.nlm.nih.gov/COG). In addition, a supplement to the COGs is available, in which proteins encoded in the genomes of two multicellular eukaryotes, the nematode Caenorhabditis elegans and the fruit fly Drosophila melanogaster, and shared with bacteria and/or archaea were included. The new features added to the COG database include information pages with structural and functional details on each COG and literature references, improvements of the COGNITOR program that is used to fit new proteins into the COGs, and classification of genomes and COGs constructed by using principal component analysis.

Genomic evidence for a complete sexual cycle in Candida albicans

Candida albicans is a diploid fungus that has become a medically important opportunistic pathogen in immunocompromised individuals. We have sequenced the C. albicans genome to 10.4-fold coverage and performed a comparative genomic analysis between C. albicans and Saccharomyces cerevisiae with the objective of assessing whether Candida possesses a genetic repertoire that could support a complete sexual cycle. Analyzing over 500 genes important for sexual differentiation in S. cerevisiae, we find many homologues of genes that are implicated in the initiation of meiosis, chromosome recombination, and the formation of synaptonemal complexes. However, others are striking in their absence. C. albicans seems to have homologues of all of the elements of a functional pheromone response pathway involved in mating in S. cerevisiae but lacks many homologues of S. cerevisiae genes for meiosis. Other meiotic gene homologues in organisms ranging from filamentous fungi to Drosophila melanogaster and Caenorhabditis elegans were also found in the C. albicans genome, suggesting potential alternative mechanisms of genetic exchange.

Complete genomic sequence of Pasteurella multocida,Pm70

We present here the complete genome sequence of a common avian clone of Pasteurella multocida, Pm70. The genome of Pm70 is a single circular chromosome 2,257,487 base pairs in length and contains 2,014 predicted coding regions, 6 ribosomal RNA operons, and 57 tRNAs. Genome-scale evolutionary analyses based on pairwise comparisons of 1,197 orthologous sequences between P. multocida, Haemophilus influenzae, and Escherichia coli suggest that P. multocida and H. influenzae diverged ≈270 million years ago and the γ subdivision of the proteobacteria radiated about 680 million years ago. Two previously undescribed open reading frames, accounting for ≈1% of the genome, encode large proteins with homology to the virulence-associated filamentous hemagglutinin of Bordetella pertussis. Consistent with the critical role of iron in the survival of many microbial pathogens, in silico and whole-genome microarray analyses identified more than 50 Pm70 genes with a potential role in iron acquisition and metabolism. Overall, the complete genomic sequence and preliminary functional analyses provide a foundation for future research into the mechanisms of pathogenesis and host specificity of this important multispecies pathogen.

Complete genome sequence of Caulobacter crescentus

The complete genome sequence of Caulobacter crescentus was determined to be 4,016,942 base pairs in a single circular chromosome encoding 3,767 genes. This organism, which grows in a dilute aquatic environment, coordinates the cell division cycle and multiple cell differentiation events. With the annotated genome sequence, a full description of the genetic network that controls bacterial differentiation, cell growth, and cell cycle progression is within reach. Two-component signal transduction proteins are known to play a significant role in cell cycle progression. Genome analysis revealed that the C. crescentus genome encodes a significantly higher number of these signaling proteins (105) than any bacterial genome sequenced thus far. Another regulatory mechanism involved in cell cycle progression is DNA methylation. The occurrence of the recognition sequence for an essential DNA methylating enzyme that is required for cell cycle regulation is severely limited and shows a bias to intergenic regions. The genome contains multiple clusters of genes encoding proteins essential for survival in a nutrient poor habitat. Included are those involved in chemotaxis, outer membrane channel function, degradation of aromatic ring compounds, and the breakdown of plant-derived carbon sources, in addition to many extracytoplasmic function sigma factors, providing the organism with the ability to respond to a wide range of environmental fluctuations. C. crescentus is, to our knowledge, the first free-living α-class proteobacterium to be sequenced and will serve as a foundation for exploring the biology of this group of bacteria, which includes the obligate endosymbiont and human pathogen Rickettsia prowazekii, the plant pathogen Agrobacterium tumefaciens, and the bovine and human pathogen Brucella abortus.

The human XPG gene: gene architecture, alternative splicing and single nucleotide polymorphisms

Defects in the XPG DNA repair endonuclease gene can result in the cancer-prone disorders xeroderma pigmentosum (XP) or the XP–Cockayne syndrome complex. While the XPG cDNA sequence was known, determination of the genomic sequence was required to understand its different functions. In cells from normal donors, we found that the genomic sequence of the human XPG gene spans 30 kb, contains 15 exons that range from 61 to 1074 bp and 14 introns that range from 250 to 5763 bp. Analysis of the splice donor and acceptor sites using an information theory-based approach revealed three splice sites with low information content, which are components of the minor (U12) spliceosome. We identified six alternatively spliced XPG mRNA isoforms in cells from normal donors and from XPG patients: partial deletion of exon 8, partial retention of intron 8, two with alternative exons (in introns 1 and 6) and two that retained complete introns (introns 3 and 9). The amount of alternatively spliced XPG mRNA isoforms varied in different tissues. Most alternative splice donor and acceptor sites had a relatively high information content, but one has the U12 spliceosome sequence. A single nucleotide polymorphism has allele frequencies of 0.74 for 3507G and 0.26 for 3507C in 91 donors. The human XPG gene contains multiple splice sites with low information content in association with multiple alternatively spliced isoforms of XPG mRNA.

Complete genome sequence of an M1 strain of Streptococcus pyogenes

The 1,852,442-bp sequence of an M1 strain of Streptococcus pyogenes, a Gram-positive pathogen, has been determined and contains 1,752 predicted protein-encoding genes. Approximately one-third of these genes have no identifiable function, with the remainder falling into previously characterized categories of known microbial function. Consistent with the observation that S. pyogenes is responsible for a wider variety of human disease than any other bacterial species, more than 40 putative virulence-associated genes have been identified. Additional genes have been identified that encode proteins likely associated with microbial “molecular mimicry” of host characteristics and involved in rheumatic fever or acute glomerulonephritis. The complete or partial sequence of four different bacteriophage genomes is also present, with each containing genes for one or more previously undiscovered superantigen-like proteins. These prophage-associated genes encode at least six potential virulence factors, emphasizing the importance of bacteriophages in horizontal gene transfer and a possible mechanism for generating new strains with increased pathogenic potential.

Does recombination improve selection on codon usage? Lessons from nematode and fly complete genomes

Understanding the factors responsible for variations in mutation patterns and selection efficacy along chromosomes is a prerequisite for deciphering genome sequences. Population genetics models predict a positive correlation between the efficacy of selection at a given locus and the local rate of recombination because of Hill–Robertson effects. Codon usage is considered one of the most striking examples that support this prediction at the molecular level. In a wide range of species including Caenorhabditis elegans and Drosophila melanogaster, codon usage is essentially shaped by selection acting for translational efficiency. Codon usage bias correlates positively with recombination rate in Drosophila, apparently supporting the hypothesis that selection on codon usage is improved by recombination. Here we present an exhaustive analysis of codon usage in C. elegans and D. melanogaster complete genomes. We show that in both genomes there is a positive correlation between recombination rate and the frequency of optimal codons. However, we demonstrate that in both species, this effect is due to a mutational bias toward G and C bases in regions of high recombination rate, possibly as a direct consequence of the recombination process. The correlation between codon usage bias and recombination rate in these species appears to be essentially determined by recombination-dependent mutational patterns, rather than selective effects. This result highlights that it is necessary to take into account the mutagenic effect of recombination to understand the evolutionary role and impact of recombination.

The complete genome of the crenarchaeon Sulfolobus solfataricus P2

The genome of the crenarchaeon Sulfolobus solfataricus P2 contains 2,992,245 bp on a single chromosome and encodes 2,977 proteins and many RNAs. One-third of the encoded proteins have no detectable homologs in other sequenced genomes. Moreover, 40% appear to be archaeal-specific, and only 12% and 2.3% are shared exclusively with bacteria and eukarya, respectively. The genome shows a high level of plasticity with 200 diverse insertion sequence elements, many putative nonautonomous mobile elements, and evidence of integrase-mediated insertion events. There are also long clusters of regularly spaced tandem repeats. Different transfer systems are used for the uptake of inorganic and organic solutes, and a wealth of intracellular and extracellular proteases, sugar, and sulfur metabolizing enzymes are encoded, as well as enzymes of the central metabolic pathways and motility proteins. The major metabolic electron carrier is not NADH as in bacteria and eukarya but probably ferredoxin. The essential components required for DNA replication, DNA repair and recombination, the cell cycle, transcriptional initiation and translation, but not DNA folding, show a strong eukaryal character with many archaeal-specific features. The results illustrate major differences between crenarchaea and euryarchaea, especially for their DNA replication mechanism and cell cycle processes and their translational apparatus.

The cranium of Parapithecus grangeri, an Egyptian Oligocene anthropoidean primate

A nearly complete skull of Parapithecus grangeri from the early Oligocene of Egypt is described. The specimen is relatively undistorted and is undoubtedly the most complete higher primate skull yet found in the African Oligocene, which also makes it the most complete Oligocene primate cranium worldwide. Belonging in superfamily Parapithecoidea, a group regarded by some as the sister group to all other Anthropoidea, this skull reveals important information about the radiation of stem anthropoideans. This cranium is about 15% larger than size estimates based on a fragmentary cranium of its contemporary and close relative Apidium phiomense. It is about the same size as that of the gray gentle lemur, Hapalemur griseus, or of platyrrhines such as the owl monkey, Aotus trivirgatus, or the titi monkey, Callicebus torquatus. Comparatively small orbits and size differences in jaws and teeth show it was both diurnal and dimorphic. This is the only specimen of the species that shows (from sockets) that there were four small upper incisors. Several mandibular specimens of the species establish that there were no permanent lower incisors and that the symphysis was fused. Like other early anthropoideans this species possessed a lower encephalization quotient and less-developed orbital frontality than later anthropoideans. There is full postorbital closure and fusion of the metopic suture, and the ectotympanic forms a rim to the auditory aperture. A probable frontal/alisphenoid contact is a potentially derived resemblance to Catarrhini. A proposed separate genus for the species P. grangeri is not sustained.

Complete Cytolysis and Neonatal Lethality in Keratin 5 Knockout Mice Reveal Its Fundamental Role in Skin Integrity and in Epidermolysis Bullosa Simplex

In human patients, a wide range of mutations in keratin (K) 5 or K14 lead to the blistering skin disorder epidermolysis bullosa simplex. Given that K14 deficiency does not lead to the ablation of a basal cell cytoskeleton because of a compensatory role of K15, we have investigated the requirement for the keratin cytoskeleton in basal cells by inactivating the K5 gene in mice. We report that the K5−/− mice die shortly after birth, lack keratin filaments in the basal epidermis, and are more severely affected than K14−/− mice. In contrast to the K14−/− mice, we detected a strong induction of the wound-healing keratin K6 in the suprabasal epidermis of cytolyzed areas of postnatal K5−/− mice. In addition, K5 and K14 mice differed with respect to tongue lesions. Moreover, we show that in the absence of K5 and other type II keratins, residual K14 and K15 aggregated along hemidesmosomes, demonstrating that individual keratins without a partner are stable in vivo. Our data indicate that K5 may be the natural partner of K15 and K17. We suggest that K5 null mutations may be lethal in human epidermolysis bullosa simplex patients.

Bacteriophage T4 gene 41 helicase and gene 59 helicase-loading protein: A versatile couple with roles in replication and recombination

Bacteriophage T4 uses two modes of replication initiation: origin-dependent replication early in infection and recombination-dependent replication at later times. The same relatively simple complex of T4 replication proteins is responsible for both modes of DNA synthesis. Thus the mechanism for loading the T4 41 helicase must be versatile enough to allow it to be loaded on R loops created by transcription at several origins, on D loops created by recombination, and on stalled replication forks. T4 59 helicase-loading protein is a small, basic, almost completely α-helical protein whose N-terminal domain has structural similarity to high mobility group family proteins. In this paper we review recent evidence that 59 protein recognizes specific structures rather than specific sequences. It binds and loads the helicase on replication forks and on three- and four-stranded (Holliday junction) recombination structures, without sequence specificity. We summarize our experiments showing that purified T4 enzymes catalyze complete unidirectional replication of a plasmid containing the T4 ori(uvsY) origin, with a preformed R loop at the position of the R loop identified at this origin in vivo. This replication depends on the 41 helicase and is strongly stimulated by 59 protein. Moreover, the helicase-loading protein helps to coordinate leading and lagging strand synthesis by blocking replication on the ori(uvsY) R loop plasmid until the helicase is loaded. The T4 enzymes also can replicate plasmids with R loops that do not have a T4 origin sequence, but only if the R loops are within an easily unwound DNA sequence.

A minimal gene set for cellular life derived by comparison of complete bacterial genomes.

The recently sequenced genome of the parasitic bacterium Mycoplasma genitalium contains only 468 identified protein-coding genes that have been dubbed a minimal gene complement [Fraser, C.M., Gocayne, J.D., White, O., Adams, M.D., Clayton, R.A., et al. (1995) Science 270, 397-403]. Although the M. genitalium gene complement is indeed the smallest among known cellular life forms, there is no evidence that it is the minimal self-sufficient gene set. To derive such a set, we compared the 468 predicted M. genitalium protein sequences with the 1703 protein sequences encoded by the other completely sequenced small bacterial genome, that of Haemophilus influenzae. M. genitalium and H. influenzae belong to two ancient bacterial lineages, i.e., Gram-positive and Gram-negative bacteria, respectively. Therefore, the genes that are conserved in these two bacteria are almost certainly essential for cellular function. It is this category of genes that is most likely to approximate the minimal gene set. We found that 240 M. genitalium genes have orthologs among the genes of H. influenzae. This collection of genes falls short of comprising the minimal set as some enzymes responsible for intermediate steps in essential pathways are missing. The apparent reason for this is the phenomenon that we call nonorthologous gene displacement when the same function is fulfilled by nonorthologous proteins in two organisms. We identified 22 nonorthologous displacements and supplemented the set of orthologs with the respective M. genitalium genes. After examining the resulting list of 262 genes for possible functional redundancy and for the presence of apparently parasite-specific genes, 6 genes were removed. We suggest that the remaining 256 genes are close to the minimal gene set that is necessary and sufficient to sustain the existence of a modern-type cell. Most of the proteins encoded by the genes from the minimal set have eukaryotic or archaeal homologs but seven key proteins of DNA replication do not. We speculate that the last common ancestor of the three primary kingdoms had an RNA genome. Possibilities are explored to further reduce the minimal set to model a primitive cell that might have existed at a very early stage of life evolution.

Primordial emergence of the recombination activating gene 1 (RAG1): sequence of the complete shark gene indicates homology to microbial integrases.

The rearrangement of antibody and T-cell receptor gene segments is indispensable to the vertebrate immune response. All extant jawed vertebrates can rearrange these gene segments. This ability is conferred by the recombination activating genes I and II (RAG I and RAG II). To elucidate their origin and function, the cDNA encoding RAG I from a member of the most ancient class of extant gnathostomes, the Carcharhine sharks, was characterized. Homology domains identified within shark RAG I prompted sequence comparison analyses that suggested similarity of the RAG I and II genes, respectively, to the integrase family genes and integration host factor genes of the bacterial site-specific recombination system. Thus, the apparent explosive evolution (or "big bang") of the ancestral immune system may have been initiated by a transfer of microbial site-specific recombinases.

Complete replication of an animal virus and maintenance of expression vectors derived from it in Saccharomyces cerevisiae.

Here we describe the first instances to our knowledge of animal virus genome replication, and of de novo synthesis of infectious virions by a nonendogenous virus, in the yeast Saccharomyces cerevisiae, whose versatile genetics offers significant advantages for studying viral replication and virus-host interactions. Flock house virus (FHV) is the most extensively studied member of the Nodaviridae family of (+) strand RNA animal viruses. Transfection of yeast with FHV genomic RNA induced viral RNA replication, transcription, and assembly of infectious virions. Genome replication and virus synthesis were robust: all replicating FHV RNA species were readily detected in yeast by Northern blot analysis and yields of virions per cell were similar to those from Drosophila cells. We also describe in vivo expression and maintenance of a selectable yeast marker gene from an engineered FHV RNA derivative dependent on FHV-directed RNA replication. Use of these approaches with FHV and their possible extension to other viruses should facilitate identification and characterization of host factors required for genomic replication, gene expression, and virion assembly.

Complete structure of the human alpha-albumin gene, a new member of the serum albumin multigene family.

The nucleotide sequence of the human alpha-albumin gene, including 887 bp of the 5'-flanking region and 1311 bp of the 3-flanking region (24,454 in total), was determined from three overlapping lambda phage clones. The sequence spans 22,256 bp from the cap site to the polyadenylylation site, revealing a gene structure of 15 exons separated by 14 introns. The methionine initiation codon ATG is within exon 1; the termination codon TGA is within exon 14. Exon 15 is entirely untranslated and contains the polyadenylylation signal AATAAA. The deduced polypeptide chain is composed of a 21-amino-acid leader peptide, followed by 578 amino acids of the mature protein. There are seven repetitive DNA elements (Alu and Kpn) in the introns and 3-flanking region. The sizes of the 15 alpha-albumin exons match closely those of the albumin, alpha-fetoprotein, and vitamin D-binding protein genes. The exons are symmetrically placed within the three domains of the individual proteins, and they share a characteristic codon splitting pattern that is conserved among members of the gene family. The results provide strong evidence that alpha-albumin belongs to, and most likely completes with, the serum albumin gene family. Based on structural similarity, alpha-albumin appears to be most closely related to alpha-fetoprotein. The complete structure of this family of four tandemly linked genes provides a well-characterized approximately 200 kb locus in the 4q subcentromeric region of the human genome.