The Journey to smORFland

The genome sequences completed so far contain more than 20 000 genes with unknown function and no similarity to genes in other genomes. The origin and evolution of the orphan genes is an enigma. Here, we discuss the suggestion that some orphan genes may represent pseudogenes or short fragments of genes that were functional in the genome of a common ancestor. These may be the remains of unsuccessful duplication or horizontal gene transfer events, in which the acquired sequences have entered the fragmentation process and thereby lost their similarity to genes in other species. This scenario is supported by a recent case study of orphan genes in several closely related species of Rickettsia, where full-length ancestral genes were reconstructed from sets of short, overlapping orphan genes. One of these was found to display similarity to genes encoding proteins with ankyrin-repeat domains.     

Genome degradation is an ongoing process in Rickettsia.

To study reductive evolutionary processes in bacterial genomes, we examine sequences in the Rickettsia genomes which are unconstrained by selection and evolve as pseudogenes, one of which is the metK gene, which codes for AdoMet synthetase. Here, we sequenced the metK gene and three surrounding genes in eight different species of the genus Rickettsia. The metK gene was found to contain a high incidence of deletions in six lineages, while the three genes in its surroundings were functionally conserved in all eight lineages. A more drastic example of gene degradation was identified in the metK downstream region, which contained an open reading frame in Rickettsia felis. Remnants of this open reading frame could be reconstructed in five additional species by eliminating sites of frameshift mutations and termination codons. A detailed examination of the two reconstructed genes revealed that deletions strongly predominate over insertions and that there is a strong transition bias for point mutations which is coupled to an excess of GC-to-AT substitutions. Since the molecular evolution of these inactive genes should reflect the rates and patterns of neutral mutations, our results strongly suggest that there is a high spontaneous rate of deletions as well as a strong mutation bias toward AT pairs in the Rickettsia genomes. This may explain the low genomic G + C content (29%), the small genome size (1.1 Mb), and the high noncoding content (24%), as well as the presence of several pseudogenes in the Rickettsia prowazekii genome.     

Cloning of a novel pilin-like gene from Bordetella pertussis: homology to the fim2 gene

A search for pilin genes in a Bordetella pertussis (Bp) genomic library has led to the identification of several clones which hybridize to synthetic oligonucleotides with sequences derived from amino acid sequences of Bp fimbrial subunits. One of these clones (corresponding to a gene we have named fimX) contains an open reading frame encoding a protein with a molecular weight of about 20 kD and a sequence similar but not identical to the fimbrial subunit fim2 and to other fimbrial protein sequences. In this communication we present the cloning and nucleotide sequence of the fimX gene and its homology to the fim2 gene. A genomic analysis on the positional relationship between the two genes is also presented.     

Proliferation and deterioration of Rickettsia palindromic elements

It has been suggested that Rickettsia Palindromic Elements (RPEs) have evolved as selfish DNA that mediate protein sequence evolution by being targeted to genes that code for RNA and proteins. Here, we have examined the phylogenetic depth of two RPEs that are located close to the genes encoding elongation factors Tu (tuf) and G (fus) in Rickettsia. An exceptional organization of the elongation factor genes was found in all 11 species examined, with complete or partial RPEs identified downstream of the tuf gene (RPE-tuf) in six species and of the fus gene (RPE-fus) in 10 species. A phylogenetic reconstruction shows that both RPE-tuf and RPE-fus have evolved in a manner that is consistent with the expected species divergence. The analysis provides evidence for independent loss of RPE-tuf in several species, possibly mediated by short repetitive sequences flanking the site of excision. The remaining RPE-tuf sequences evolve as neutral sequences in different stages of deterioration. Likewise, highly fragmented remnants of the RPE-fus sequence were identified in two species. This suggests that genome-specific differences in the content of RPEs are the result of recent loss rather than recent proliferation.     

Sequence analysis of the gene encoding a spotted fever group-specific intracytoplasmic protein PS120 of Rickettsia japonica

The 3,438-nucleotide (nt) sequence containing a 3,054-nt open reading frame of the gene (rps120) encoding an antigenic, intracytoplasmic, spotted fever group-specific and heat-stable 120-kilodalton protein (PS120) of Rickettsia japonica was determined. The nt and deduced 1,018 amino-acid (aa) sequences were compared to those of R. conorii since only those of this species had been determined among SFG rickettsiae. The homologies of these sequences between R. japonica and R. conorii were considerably high at 97 and 95%, respectively. These high homologies were comparable to those of beta-peptides encoded by the ompB genes among SFG rickettsiae. It was also found that the genome of R. prowazekii contained a nt sequence with 68% homology to that of the rps120 gene of R. japonica.     

Complete DNA sequence of yeast chromosome II.

In the framework of the EU genome-sequencing programmes, the complete DNA sequence of the yeast Saccharomyces cerevisiae chromosome II (807 188 bp) has been determined. At present, this is the largest eukaryotic chromosome entirely sequenced. A total of 410 open reading frames (ORFs) were identified, covering 72% of the sequence. Similarity searches revealed that 124 ORFs (30%) correspond to genes of known function, 51 ORFs (12.5%) appear to be homologues of genes whose functions are known, 52 others (12.5%) have homologues the functions of which are not well defined and another 33 of the novel putative genes (8%) exhibit a degree of similarity which is insufficient to confidently assign function. Of the genes on chromosome II, 37-45% are thus of unpredicted function. Among the novel putative genes, we found several that are related to genes that perform differentiated functions in multicellular organisms of are involved in malignancy. In addition to a compact arrangement of potential protein coding sequences, the analysis of this chromosome confirmed general chromosome patterns but also revealed particular novel features of chromosomal organization. Alternating regional variations in average base composition correlate with variations in local gene density along chromosome II, as observed in chromosomes XI and III. We propose that functional ARS elements are preferably located in the AT-rich regions that have a spacing of approximately 110 kb. Similarly, the 13 tRNA genes and the three Ty elements of chromosome II are found in AT-rich regions. In chromosome II, the distribution of coding sequences between the two strands is biased, with a ratio of 1.3:1. An interesting aspect regarding the evolution of the eukaryotic genome is the finding that chromosome II has a high degree of internal genetic redundancy, amounting to 16% of the coding capacity.     

Analysis of fragment homology among DNA sequences from cytochrome P450 families 4 and 6

Cytochrome P450s comprise a diverse superfamily of proteins that often share as little as 12% amino acid identity. Accordingly, the identification of novel gene families, subfamilies and alleles has been based primarily on ‘benchmark’ levels of global amino acid identity and, more recently, phylogenetic relatedness has been used to resolve ambiguous relationships. However, PCR-based cloning strategies have resulted in a large increase in the number of short DNA sequences, particularly among insects. Many of these fragments remain unnamed and even their gene family membership remains unknown due to the uncertainty as to whether these fragments accurately reflect the levels of sequence identity or patterns of evolutionary divergence exhibited by groups of full-length P450 sequences. As a result, the nature of P450 diversity among insects remains obscure. In this paper, P450 sequences belonging to the two major gene families, that is, families CYP4 and CYP6, are analyzed by comparing segments to full-length sequences. A parameter called the Segmental Divergence Index is used to characterize segments of P450s with respect to the degree to which they mirror the divergence of full-length sequences. The reliability of these fragments in phylogenetic analyses is also tested. The results of these analyses support the use of some commonly amplified P450 gene fragments for naming new P450 genes and for studies of P450 diversity and evolution. This revised version was published online in July 2006 with corrections to the Cover Date.     

Structure and expression of the Chinese hamster thymidine kinase gene.

My colleagues and I have cloned a nearly full-length Chinese hamster thymidine kinase (TK) cDNA in a lambda gt10 vector and characterized this cDNA by nucleotide sequencing. The hamster TK protein is encoded in this cDNA by a 702-base-pair open reading frame which specifies a 25,625-dalton protein closely homologous to the previously described human and chicken TK proteins. Using cDNA nucleotide sequence data in conjunction with sequence data derived from selected subclones of the hamster TK gene recombinant phage lambda HaTK.5, we have resolved the structure of the TK gene, finding the 1,219 base pairs of the cDNA sequence to be distributed through 11.2 kilobases of genomic DNA in at least seven exon segments. In addition, we have constructed a variety of Chinese hamster TK minigenes and exonuclease III-S1 derivatives of these genes which have permitted us to define the limits of the Chinese hamster TK gene promoter and demonstrate that efficient TK transformation of Ltk- cells by TK minigenes depends on the presence of both TK intervening sequences and sequences 3' to the site of mRNA polyadenylation.     

斑马鱼DrSpin-1基因的克隆和特征分析

Spin蛋白家族是具有Spin/Ssty保守结构域并在配子发生过程中发挥关键作用的一类分子。研究利用简并引物PCR,从斑马鱼成熟卵母细胞SMART cDNA文库中筛选到260 bp的DrSpin-1和DrSpin-2部分序列,经序列同源性比对,斑马鱼DrSpin-1的部分氨基酸序列与银鲫CagSpin一致性高达81%。利用RACEPCR从该cDNA文库中获得斑马鱼DrSpin-1的全长cDNA序列。序列分析表明,DrSpin-1全长cDNA为1082 bp,开放阅读框771 bp,编码257个氨基酸,具有三个Spin/Ssty保守域,8个可能的磷酸化位点,初步确定斑马鱼DrSpin-1是Spin基因家族成员。斑马鱼DrSpin-1蛋白与已报道的鱼类Spin蛋白多重序列比对表明,DrSpin-1蛋白与银鲫CagSpin蛋白同源性最高。可以推测克隆得到的斑马鱼DrSpin-1与已知功能的银鲫CagSpin具有相近的表达谱和生物学功能,可能在配子发生和受精过程中发挥重要作用。     

Identification of new protein kinase-related genes in three herpesviruses, herpes simplex virus, varicella-zoster virus, and Epstein-Barr virus.

By using amino acid sequence patterns (motifs) diagnostic of conserved regions within the catalytic domains of protein kinases, homologous open reading frames of three herpesviruses were identified as protein kinase-related genes. The three sequences, herpes simplex virus gene UL13, varicella-zoster virus gene 47, and Epstein-Barr virus gene BGLF4, resemble serine/threonine kinases rather than tyrosine kinases.     

A hyper-recombination mutation in S. cerevisiae identifies a novel eukaryotic topoisomerase

A hyper-recombination mutation was isolated that causes an increase in recombination between short repeated delta sequences surrounding the SUP4-omicron gene in S. cerevisiae. The wild-type copy of this gene was cloned by complementation of one of its pleiotropic phenotypes, slow growth. DNA sequence of the clone revealed a 656 amino acid open reading frame capable of encoding a protein homologous to the bacterial type I topoisomerase. No homology was detected with previously identified eukaryotic topoisomerases. Construction of double mutants with either of the two known yeast topoisomerase genes revealed synergistic effects on growth suggesting overlapping functions. Expression of bacterial topoisomerase I in yeast can fully complement the slow growth defect of a null mutation. We have named this locus TOP3 and suggest that it defines a novel eukaryotic topoisomerase gene.