Similarities in eukaryotic genomes

1. The degree of overlap between the human genome and that of other eukaryotes is considered. Biochemical and molecular studies have shown that all eukaryotic organisms evolved from a common progenator that lived several billion years ago. 2. From a geneological point of view, all eukaryotes are related and their genes are all descended from common ancestors. 3. However, most of the DNA in eukaryotic genomes is not transcribed and has been free to drift in nucleotide sequence. Therefore, the question of overlap can only be applied meaningfully to the few per cent of the genome that is expressed. 4. During the last billion years many genes have duplicated and diverged and new genes have been formed by accretion of domains copied from other genes (exon shuffling). 5. The rate of genetic divergence has been such that only a few portions coding for pieces of highly conserved proteins are still shared by all eukaryotes including those that diverged over 600 million years ago. 6. On the other hand, a fairly large number of shared genes can be recognized among species that separated within the last few hundred million years. 7. Human genes have a high degree of identity with homologs in closely related organisms such as other mammals and a decreasing level of identity with their homologs in more distantly related species.     

DoriC: a database of oriC regions in bacterial genomes

Replication origins (oriCs) of bacterial genomes currently available in GenBank have been predicted by using a systematic method comprising the Z-curve analysis for nucleotide distribution asymmetry, DnaA box distribution, genes adjacent to candidate oriCs and phylogenetic relationships. These oriCs are organized into a MySQL database, DoriC, which provides extensive information and graphical views of the oriC regions. In addition, users can Blast a query sequence or even a whole genome against DoriC to find a homologous one. DoriC will be updated timely and the latest version is DoriC 1.8, in which oriCs of 425 genomes (468 chromosomes) are identified. AVAILABILITY: DoriC can be accessed from http://tubic.tju.edu.cn/doric/. SUPPLEMENTARY INFORMATION: Supplementary data are available at http://tubic.tju.edu.cn/doric/supplementary.htm.     

Differentiation of regions with atypical oligonucleotide composition in bacterial genomes

Background  

Complete sequencing of bacterial genomes has become a common technique of present day microbiology. Thereafter, data mining in the complete sequence is an essential step. New in silico methods are needed that rapidly identify the major features of genome organization and facilitate the prediction of the functional class of ORFs. We tested the usefulness of local oligonucleotide usage (OU) patterns to recognize and differentiate types of atypical oligonucleotide composition in DNA sequences of bacterial genomes.     

Noise-reduction filtering for accurate detection of replication termini in bacterial genomes

Bacterial chromosomes are highly polarized in their nucleotide composition through mutational selection related to replication. Using compositional skews such as the GC skew, replication origin and terminus can be predicted in silico by observing the shift points. However, the genome sequence is affected by myriad functional requirements and selection on numerous subgenomic features, and elimination of this "noise" should lead to better predictions. Here, we present a noise-reduction approach that uses low-pass filtering through Fast Fourier transform coupled with cumulative skew graphs. It increases the prediction accuracy of the replication termini compared with previously documented methods based on genomic base composition.     

Effects of mutations of the initiation nucleotides on hepatitis C virus RNA replication in the cell

Replication of nearly all RNA viruses depends on a virus-encoded RNA-dependent RNA polymerase (RdRp). Our earlier work found that purified recombinant hepatitis C virus (HCV) RdRp (NS5B) was able to initiate RNA synthesis de novo by using purine (A and G) but not pyrimidine (C and U) nucleotides (G. Luo et al., J. Virol. 74:851-863, 2000). For most human RNA viruses, the initiation nucleotides of both positive- and negative-strand RNAs were found to be either an adenylate (A) or guanylate (G). To determine the nucleotide used for initiation and control of HCV RNA replication, a genetic mutagenesis analysis of the nucleotides at the very 5' and 3' ends of HCV RNAs was performed by using a cell-based HCV replicon replication system. Either a G or an A at the 5' end of HCV genomic RNA was able to efficiently induce cell colony formation, whereas a nucleotide C at the 5' end dramatically reduced the efficiency of cell colony formation. Likewise, the 3'-end nucleotide U-to-C mutation did not significantly affect the efficiency of cell colony formation. In contrast, a U-to-G mutation at the 3' end caused a remarkable decrease in cell colony formation, and a U-to-A mutation resulted in a complete abolition of cell colony formation. Sequence analysis of the HCV replicon RNAs recovered from G418-resistant Huh7 cells revealed several interesting findings. First, the 5'-end nucleotide G of the replicon RNA was changed to an A upon multiple rounds of replication. Second, the nucleotide A at the 5' end was stably maintained among all replicon RNAs isolated from Huh7 cells transfected with an RNA with a 5'-end A. Third, initiation of HCV RNA replication with a CTP resulted in a >10-fold reduction in the levels of HCV RNAs, suggesting that initiation of RNA replication with CTP was very inefficient. Fourth, the 3'-end nucleotide U-to-C and -G mutations were all reverted back to a wild-type nucleotide U. In addition, extra U and UU residues were identified at the 3' ends of revertants recovered from Huh7 cells transfected with an RNA with a nucleotide G at the 3' end. We also determined the 5'-end nucleotide of positive-strand RNA of some clinical HCV isolates. Either G or A was identified at the 5' end of HCV RNA genome depending on the specific HCV isolate. Collectively, these findings demonstrate that replication of positive-strand HCV RNA was preferentially initiated with purine nucleotides (ATP and GTP), whereas the negative-strand HCV RNA replication is invariably initiated with an ATP.     

Crystal structures of two intermediates in the assembly of the papillomavirus replication initiation complex

Initiation of DNA replication of the papillomavirus genome is a multi-step process involving the sequential loading of viral E1 protein subunits onto the origin of replication. Here we have captured structural snapshots of two sequential steps in the assembly process. Initially, an E1 dimer binds to adjacent major grooves on one face of the double helix; a second dimer then binds to another face of the helix. Each E1 monomer has two DNA-binding modules: a DNA-binding loop, which binds to one DNA strand and a DNA-binding helix, which binds to the opposite strand. The nature of DNA binding suggests a mechanism for the transition between double- and single-stranded DNA binding that is implicit in the progression to a functional helicase.     

Ordered assembly of nucleoprotein structures at the bacteriophage lambda replication origin during the initiation of DNA replication

Replication of the chromosome of bacteriophage lambda depends on the cooperative action of two phage-coded proteins and seven replication and heat shock proteins from its Escherichia coli host. As previously described, the first stage in this process is the binding of multiple copies of the lambda O initiator to the lambda replication origin (ori lambda) to form the nucleosomelike O-some. The O-some serves to localize subsequent protein-protein and protein-DNA interactions involved in the initiation of lambda DNA replication to ori lambda. To study these interactions, we have developed a sensitive immunoblotting protocol that permits the protein constituents of complex nucleoprotein structures to be identified. Using this approach, we have defined a series of sequential protein assembly and protein disassembly events that occur at ori lambda during the initiation of lambda DNA replication. A second-stage ori lambda.O (lambda O protein).P (lambda P protein).DnaB nucleoprotein structure is formed when O, P, and E. coli DnaB helicase are incubated with ori lambda DNA. In a third-stage reaction the E. coli DnaJ heat shock protein specifically binds to the second-stage structure to form an ori lambda.O.P.DnaB.DnaJ complex. Each of the nucleoprotein structures formed in the first three stages was isolated and shown to be a physiological intermediate in the initiation of lambda DNA replication. The E. coli DnaK heat shock protein can bind to any of these early stage nucleoprotein structures, and in a fourth-stage reaction a complete ori lambda.O.P.DnaB.DnaJ.DnaK initiation complex is assembled. Addition of ATP to the reaction enables the DnaK and DnaJ heat shock proteins to mediate a partial disassembly of the fourth-stage complex. These protein disassembly reactions activate the intrinsic helicase activity of DnaB and result in localized unwinding of the ori lambda template. The protein disassembly reactions are described in the accompanying articles.     

Identification of variable regions in the genomes of tubercle bacilli using bacterial artificial chromosome arrays

Whole-genome comparisons of the tubercle bacilli were undertaken using ordered bacterial artificial chromosome (BAC) libraries of Mycobacterium tuberculosis and the vaccine strain, Mycobacterium bovis BCG-Pasteur, together with the complete genome sequence of M. tuberculosis H37Rv. Restriction-digested BAC arrays of M. tuberculosis H37Rv were used in hybridization experiments with radiolabelled M. bovis BCG genomic DNA to reveal the presence of 10 deletions (RD1-RD10) relative to M. tuberculosis. Seven of these regions, RD4-RD10, were also found to be deleted from M. bovis, with the three M. bovis BCG-specific deletions being identical to the RD1-RD3 loci described previously. The distribution of RD4-RD10 in Mycobacterium africanum resembles that of M. tuberculosis more closely than that of M. bovis, whereas an intermediate arrangement was found in Mycobacterium microti, suggesting that the corresponding genes may affect host range and virulence of the various tubercle bacilli. Among the known products encoded by these loci are a copy of the proposed mycobacterial invasin Mce, three phospholipases, several PE, PPE and ESAT-6 proteins, epoxide hydrolase and an insertion sequence. In a complementary approach, direct comparison of BACs uncovered a third class of deletions consisting of two M. tuberculosis H37Rv loci, RvD1 and RvD2, deleted from the genome relative to M. bovis BCG and M. bovis. These deletions affect a further seven genes, including a fourth phospholipase, plcD. In summary, the insertions and deletions described here have important implications for our understanding of the evolution of the tubercle complex.     

Recognizing the pseudogenes in bacterial genomes

Pseudogenes are now known to be a regular feature of bacterial genomes and are found in particularly high numbers within the genomes of recently emerged bacterial pathogens. As most pseudogenes are recognized by sequence alignments, we use newly available genomic sequences to identify the pseudogenes in 11 genomes from 4 bacterial genera, each of which contains at least 1 human pathogen. The numbers of pseudogenes range from 27 in Staphylococcus aureus MW2 to 337 in Yersinia pestis CO92 (e.g. 1–8% of the annotated genes in the genome). Most pseudogenes are formed by small frameshifting indels, but because stop codons are A + T-rich, the two low-G + C Gram-positive taxa (Streptococcus and Staphylococcus) have relatively high fractions of pseudogenes generated by nonsense mutations when compared with more G + C-rich genomes. Over half of the pseudogenes are produced from genes whose original functions were annotated as ‘hypothetical’ or ‘unknown’; however, several broadly distributed genes involved in nucleotide processing, repair or replication have become pseudogenes in one of the sequenced Vibrio vulnificus genomes. Although many of our comparisons involved closely related strains with broadly overlapping gene inventories, each genome contains a largely unique set of pseudogenes, suggesting that pseudogenes are formed and eliminated relatively rapidly from most bacterial genomes.     

The most deviated codon position in AT-rich bacterial genomes: a function related analysis

We have performed systematic study on more than 120 archaeal and bacterial genomes. Based on the index proposed in the current paper, clear patterns are observed showing the relation between the base compositional deviation at three codon positions and the genomic GC content. For AT-rich genomes, the Most Deviated Codon Position (MDCP) is the 1st codon position, while for GC-rich genomes, MDCP appears at the 2nd or 3rd codon position alternatively. According to MDCP, the CDSs of a genome can be classified into two types: typical and atypical. In AT-rich genomes the typical represent the majority and account for about 3/4 of all the CDSs. Based on the functional classification of COG database, the two types of CDSs are examined. An apparent bias of distribution is observed that the CDSs with the function of 'information processing' are more likely to present in typical type.     

Open chromatin structures regulate the efficiencies of pre-RC formation and replication initiation in Epstein-Barr virus

Whether or not metazoan replication initiates at random or specific but flexible sites is an unsolved question. The lack of sequence specificity in origin recognition complex (ORC) DNA binding complicates genome-scale chromatin immunoprecipitation (ChIP)-based studies. Epstein-Barr virus (EBV) persists as chromatinized minichromosomes that are replicated by the host replication machinery. We used EBV to investigate the link between zones of pre-replication complex (pre-RC) assembly, replication initiation, and micrococcal nuclease (MNase) sensitivity at different cell cycle stages in a genome-wide fashion. The dyad symmetry element (DS) of EBV's latent origin, a well-established and very efficient pre-RC assembly region, served as an internal control. We identified 64 pre-RC zones that correlate spatially with 57 short nascent strand (SNS) zones. MNase experiments revealed that pre-RC and SNS zones were linked to regions of increased MNase sensitivity, which is a marker of origin strength. Interestingly, although spatially correlated, pre-RC and SNS zones were characterized by different features. We propose that pre-RCs are formed at flexible but distinct sites, from which only a few are activated per single genome and cell cycle.     

Polypurine.polypyrimidine sequences in complete bacterial genomes: preference for polypurines in protein-coding regions

The genomes of Methanococcus jannaschii, Mycoplasma genitalium, Haemophilus influenzae, Archaeoglobus fulgidus, Helicobacter pylori, Treponema pallidum, Borrelia burgdorferri, Rickettsia prowazekeii, Mycobacterium tuberculosis, Methanobacterium thermoautotrophicum, Synechocystis sp. PCC6803, Bacillus subtilis, Chlamydia trachomatis, Pyrococcus horikoshii, Aquifex aeolicus, Mycoplasma pneumoniae and Escherichia coli have been analysed for the presence of polypurine.polypyrimidine tracts, in order to understand their distribution in these genomes. We observed a variation in abundance of such sequences in these bacteria, with the archaeal genomes forming a high-abundance group and the canonical eubacteria forming a low-abundance group. The genomes of M. tuberculosis and A. aeolicus are unique among the organisms analysed here in the abnormal underrepresentation and overrepresentation of polypurine.polypyrimidine, respectively. We also observe a strand bias, i.e., a preferential occurrence of polypurines in coding strands. It varies widely among the bacteria, from the very high bias in M. jannaschii to the slightly inverse bias in the parasitic genomes of T. pallidum and C. trachomatis. The extent of strand bias, however, cannot be explained on the basis of the GC-content of the genome, use of all-purine codons or an excess in the amino acids that are encoded by such codons. The probable causes and effects of this phenomenon are discussed.     

Evolution of bacterial genomes

This review examines evolution of bacterial genomes with an emphasis on RNA based life, the transition to functional DNA and small evolving genomes (possibly plasmids) that led to larger, functional bacterial genomes.     

Presence of isochore structures in reptile genomes suggested by the relationship between GC contents of intron regions and those of coding regions

Vertebrate genomes are mosaics of isochores. On the assumption that marked differences exist in the isochore structure between warm-blooded and cold-blooded animals, variations among vertebrates were previously attributed to adaptation to homeothermy. However, based on the data of coding regions from representatives of extant vertebrates, including a turtle, a crocodile (Archosauromorpha) and a few kinds of snakes (Lepidosauromorpha), it was recently hypothesized that the common ancestors of mammals, birds and extant reptiles already had the "warm-blooded" isochore structure. To test this hypothesis, the nucleotide sequences of alpha-globin genes including non-coding regions (introns) from two snakes, N. kaouthia and E. climacophora, were determined (accession number: AB104824, AB104825). The correlation between the GC contents in the introns and exons of alpha-globin genes from snakes and those from other vertebrates supports the above hypothesis. Similar analysis using data for exons and introns of other genes obtained from the GenBank (Release 131) also support the above hypothesis.     

Cell-cycle-specific initiation of replication

The following characteristics are relevant when replication of chromosomes and plasmids is discussed in relation to the cell cycle: the timing or replication, the selection of molecules for replication, and the coordination of multiple initiation events within a single cell cycle. Several fundamentally different methods have been used to study these processes: Meselson—Stahl density-shift experiments, experiments with the so-called‘baby machine', sorting of cells according to size, and flow cytometry. The evidence for precise timing and co-ordination of chromosome replication in Escherichia coli is overwhelming. Similarly, the high-copy-number plasmid ColE1 and the low-copy-number plasmids R1/R100 without any doubt replicate randomly throughout the cell cycle. Data about the low-copy-number plasmids F and P1 are conflicting. This calls for new types of experiments and for a better understanding of how these plasmids control their replication and partitioning.     

Computer analysis of regulatory signals in complete bacterial genomes. Translation initiation of ribosomal protein operons]

Signals of translation initiation of operons of Haemophilus influenzae ribosomal proteins were predicted. This process is regulated by the formation of secondary RNA structures to which one of the proteins encoded in a particular operon binds. In some cases, these structures imitate the region of protein binding to rRNA. Predictions are made by comparing with homologous operons of Escherichia coli and analogous regions of rRNA and by estimating the energy of secondary structure formation. It is shown that this regulatory mechanism occurs: in operons L11, S10, S15, spc, and alpha of H.influenzae and, probably, in operon S15 of Helicobacter pylori, Bacillus subtilis, and Mycoplasma genitalium.     

Identification of replication origins in prokaryotic genomes

The availability of hundreds of complete bacterial genomes has created new challenges and simultaneously opportunities for bioinformatics. In the area of statistical analysis of genomic sequences, the studies of nucleotide compositional bias and gene bias between strands and replichores paved way to the development of tools for prediction of bacterial replication origins. Only a few (about 20) origin regions for eubacteria and archaea have been proven experimentally. One reason for that may be that this is now considered as an essentially bioinformatics problem, where predictions are sufficiently reliable not to run labor-intensive experiments, unless specifically needed. Here we describe the main existing approaches to the identification of replication origin (oriC) and termination (terC) loci in prokaryotic chromosomes and characterize a number of computational tools based on various skew types and other types of evidence. We also classify the eubacterial and archaeal chromosomes by predictability of their replication origins using skew plots. Finally, we discuss possible combined approaches to the identification of the oriC sites that may be used to improve the prediction tools, in particular, the analysis of DnaA binding sites using the comparative genomic methods.