Fur homolog regulates Porphyromonas gingivalis virulence under low-iron/heme conditions through a complex regulatory network.

Porphyromonas gingivalis is a key pathogen responsible for initiation and progression of chronic periodontitis. Little is known about the regulatory mechanisms of iron and heme uptake that allow P. gingivalis to express virulence factors and survive in the hostile environment of the oral cavity, so we initiated characterization of a P. gingivalis Fur homolog (PgFur). Many Fur paralogs found in microbial genomes, including Bacteroidetes, confirm that Fur proteins have a tendency to be subjected to a sub- or even neofunctionalization process. PgFur revealed extremely high sequence divergence, which could be associated with its functional dissimilarity in comparison with other Fur homologs. A fur mutant strain constructed by insertional inactivation exhibited retarded growth during the early growth phase and a significantly lower tendency to form a homotypic biofilm on abiotic surfaces. The mutant also showed significantly weaker adherence and invasion to epithelial cells and macrophages. Transcripts of many differentially regulated genes identified in the fur mutant strain were annotated as hypothetical proteins, suggesting that PgFur can play a novel role in the regulation of gene expression. Inactivation of the fur gene resulted in decreased hmuY gene expression, increased expression of other hmu components and changes in the expression of genes encoding hemagglutinins and proteases (mainly gingipains), HtrA, some extracytoplasmic sigma factors and two-component systems. Our data suggest that PgFur can influence in vivo growth and virulence, at least in part by affecting iron/heme acquisition, allowing efficient infection through a complex regulatory network.

Comparative genomic studies suggest that the cyanobacterial endosymbionts of the amoeba Paulinella chromatophora possess an import apparatus for nuclear-encoded proteins.

Plastids evolved from free-living cyanobacteria through a process of primary endosymbiosis. The most widely accepted hypothesis derives three ancient lineages of primary plastids, i.e. those of glaucophytes, red algae and green plants, from a single cyanobacterial endosymbiosis. This hypothesis was originally predicated on the assumption that transformations of endosymbionts into organelles must be exceptionally rare because of the difficulty in establishing efficient protein trafficking between a host cell and incipient organelle. It turns out, however, that highly integrated endosymbiotic associations are more common than once thought. Among them is the amoeba Paulinella chromatophora, which harbours independently acquired cyanobacterial endosymbionts functioning as plastids. Sequencing of the Paulinella endosymbiont genome revealed an absence of essential genes for protein trafficking, suggesting their residence in the host nucleus and import of protein products back into the endosymbiont. To investigate this hypothesis, we searched the Paulinella endosymbiont genome for homologues of higher plant translocon proteins that form the import apparatus in two-membrane envelopes of primary plastids. We found homologues of Toc12, Tic21 and Tic32, but genes for other key translocon proteins (e.g. Omp85/Toc75 and Tic20) were missing. We propose that these missing genes were transferred to the Paulinella nucleus and their products are imported and integrated into the endosymbiont envelope membranes, thereby creating an effective protein import apparatus. We further suggest that other bacterial/cyanobacterial endosymbionts found in protists, plants and animals could have evolved efficient protein import systems independently and, therefore, reached the status of true cellular organelles.

Were class C iron-containing superoxide dismutases of trypanosomatid parasites initially imported into a complex plastid? A hypothesis based on analyses of their N-terminal targeting signals.

Trypanosomatid parasites possess 2 distinct iron-containing superoxide dismutases (Fe-SODs) designated SODA and SODC, both of which are targeted to their mitochondria. In contrast to SODAs that carry typical mitochondrial transit peptides, SODCs have highly unusual mitochondrial targeting signals. Our analyses clearly show that these pre-sequences are bipartite possessing a signal peptide-like domain followed by a transit peptide-like domain. Consequently, they resemble N-terminal extensions of proteins targeted to multi-membrane plastids, suggesting that trypanosomatids once contained a eukaryotic alga-derived plastid. Further support for this hypothesis comes from striking similarities in length, hydropathy profile, and amino acid composition of SODC pre-sequences to those of Euglena and dinoflagellate plastid proteins. To account for these data, we propose that the Trypanosomatidae initially possessed a gene encoding a mitochondrial Fe-SOD with a classical mitochondrial transit peptide. Before or after plastid acquisition, a gene duplication event gave rise to SODA and SODC. In a subsequent evolutionary step a signal peptide was linked to SODC, enabling its import into the plastid. When the trypanosomatid plastid subsequently was lost, natural selection favoured adaptation of the SODC N-terminal signal as a mitochondrial transit peptide and re-targeting to the mitochondrion.

The differential killing of genes by inversions in prokaryotic genomes.

We have elaborated a method which has allowed us to estimate the direction of translocation of orthologs which have changed, during the phylogeny, their positions on chromosome in respect to the leading or lagging role of DNA strands. We have shown that the relative number of translocations which have switched positions of genes from the leading to the lagging DNA strand is lower than the number of translocations which have transferred genes from the lagging strand to the leading strand of prokaryotic genomes. This paradox could be explained by assuming that the stronger mutation pressure and selection after inversion preferentially eliminate genes transferred from the leading to the lagging DNA strand.

Evolution rates of genes on leading and lagging DNA strands.

One of the main causes of bacterial chromosome asymmetry is replication-associated mutational pressure. Different rates of nucleotide substitution accumulation on leading and lagging strands implicate qualitative and quantitative differences in the accumulation of mutations in protein coding sequences lying on different DNA strands. We show that the divergence rate of orthologs situated on leading strands is lower than the divergence rate of those situated on lagging strands. The ratio of the mutation accumulation rate for sequences lying on lagging strands to that of sequences lying on leading strands is rather stable and time-independent. The divergence rate of sequences which changed their positions, with respect to the direction of replication fork movement, is not stable-sequences which have recently changed their positions are the most prone to mutation accumulation. This effect may influence estimations of evolutionary distances between species and the topology of phylogenetic trees.

Identification of a putative chromosomal replication origin from Helicobacter pylori and its interaction with the initiator protein DnaA.

The key elements of the initiation of Helicobacter pylori chromosome replication, DnaA protein and putative oriC region, have been characterized. The gene arrangement in the H.pylori dnaA region differs from that found in many other eubacterial dnaA regions (rnpA-rmpH-dnaA-dnaN-recF-gyrB). Helicobacter pylori dnaA is flanked by two open reading frames with unknown function, while dnaN-gyrB and rnpA-rmpH loci are separated from the dnaA gene by 600 and 90 kb, respectively. We show that the dnaA gene encoding initiator protein DnaA is expressed in H.pylori cells. The H.pylori DnaA protein, like other DnaA proteins, can be divided into four domains. Here we demonstrate that the C-terminal domain of H.pylori DnaA protein is responsible for DNA binding. Using in silico and in vitro studies, the putative oriC region containing five DnaA boxes has been located upstream of the dnaA gene. DNase I and gel retardation analyses show that the C-terminal domain of H.pylori DnaA protein specifically binds each of five DnaA boxes.

High correlation between the turnover of nucleotides under mutational pressure and the DNA composition.

BACKGROUND: Any DNA sequence is a result of compromise between the selection and mutation pressures exerted on it during evolution. It is difficult to estimate the relative influence of each of these pressures on the rate of accumulation of substitutions. However, it is important to discriminate between the effect of mutations, and the effect of selection, when studying the phylogenic relations between taxa. RESULTS: We have tested in computer simulations, and analytically, the available substitution matrices for many genomes, and we have found that DNA strands in equilibrium under mutational pressure have unique feature: the fraction of each type of nucleotide is linearly dependent on the time needed for substitution of half of nucleotides of a given type, with a correlation coefficient close to 1. Substitution matrices found for sequences under selection pressure do not have this property. A substitution matrix for the leading strand of the Borrelia burgdorferi genome, having reached equilibrium in computer simulation, gives a DNA sequence with nucleotide composition and asymmetry corresponding precisely to the third positions in codons of protein coding genes located on the leading strand. CONCLUSIONS: Parameters of mutational pressure allow us to count DNA composition in equilibrium with this mutational pressure. Comparing any real DNA sequence with the sequence in equilibrium it is possible to estimate the distance between these sequences, which could be used as a measure of the selection pressure. Furthermore, the parameters of the mutational pressure enable direct estimation of the relative mutation rates in any DNA sequence in the studied genome.

DNA asymmetry and the replicational mutational pressure.

The mode of replication and organisation of bacterial genomes impose asymmetry on their nucleotide composition. The asymmetry is seen in coding and non-coding sequences and is reflected in the amino acid composition of proteins. The mechanisms generating asymmetry include: unequal mutation rates connected with replication and transcription, selection forces positioning genes and signal sequences nonrandomly in the genome, and protein coding constraints on coding sequences. There are different methods of visualising and measuring the asymmetry. Some of them can assess the contribution of individual mechanisms to the observed asymmetry and those have been described in greater detail. Asymmetric mutational and selection pressures differentiate the rates of evolution of genes on leading and lagging strands. The genes relocated to the opposite strand have to adapt to a different mutational pressure or are eliminated. Translocations from leading to lagging strands are more often selected against than from lagging to leading strands. Comparison of intergenic sequences that have lost the coding function to the original genes enables finding the frequencies of the twelve substitution rates in sequences free from selection. In the absence of selection, the half-time of substitution of a given type of nucleotide is linearly correlated with the fraction of that nucleotide in the sequence.

Is there replication-associated mutational pressure in the Saccharomyces cerevisiae genome?

Compositional bias of yeast chromosomes was analysed using detrended DNA walks. Unlike eubacterial chromosomes, the yeast chromosomes did not show the specific asymmetry correlated with origin and terminus of replication. It is probably a result of a relative excess of autonomously replicating sequences (ARS) and of random choice of these sequences in each replication cycle. Nevertheless, the last ARS from both ends of chromosomes are responsible for unidirectional replication of subtelomeric sequences with pre-established leading/lagging roles of DNA strands. In these sequences a specific asymmetry is observed, resembling the asymmetry introduced by replication-associated mutational pressure into eubacterial chromosomes.

Total number of coding open reading frames in the yeast genome.

At the end of 1996 we approximated the total number of protein coding ORFs in the Saccharomyces cerevisiae genome, based on their properties, as 4700-4800. The number is much smaller than the 5800 which is widely accepted. According to our calculations, there remain about 200-300 orphans-ORFs without known function or homology to already discovered genes, which is only about 5% of the total number of genes. Our results would be questionable if the analysed set of known genes was not a statistically representative sample of the whole set of protein coding genes in the S. cerevisiae genome. Therefore, we repeated our estimation using recently updated databases. In the course of the last 18 months, previously unknown functions of about 500 genes have been found. We have used these to check our method, former results and conclusions. Our previous estimation of the total number of coding ORFs was confirmed.

Origin and properties of non-coding ORFs in the yeast genome.

In a recent paper we have estimated the total number of protein coding open reading frames (ORFs) in the Saccharomyces cerevisiae genome, based on their properties, at about 4800. This number is much smaller than the 5800-6000 which is widely accepted. In this paper we analyse differences between the set of ORFs with known phenotypes annotated in the Munich Information Centre for Protein Sequences (MIPS) database and ORFs for which the probability of coding, counted by us, is very low. We have found that many of the latter ORFs have properties of antisense sequences of coding ORFs, which suggests that they could have been generated by duplication of coding sequences. Since coding sequences generate ORFs inside themselves, with especially high frequency in the antisense sequences, we have looked for homology between known proteins and hypothetical polypeptides generated by ORFs under consideration in all the six phases. For many ORFs we have found paralogues and orthologues in phases different than the phase which had been assumed in the MIPS database as coding.

The role of the genetic code in generating new coding sequences inside existing genes.

The genetic code has a very interesting property--it generates an open reading frame (ORF) inside a coding sequence, in a specific phase of the antisense strand with much higher probability than in the random DNA sequences. Furthermore, these antisense ORFs (A-ORFs) possess the same features as real genes--the asymmetry in the nucleotide composition at the first and second positions in codons. About two thirds of the 2997 overlapping ORFs in the yeast genome possess this feature. Thus, the question arises: has this feature of the genetic code been exploited in the evolution of genes? We have searched the FASTA data bases for homologies with the antisense translation products of a specific class of genes and we have found some sequences with relatively high homology. Many of them have scores which could be randomly found in the searched data bases with a probability lower than 10(-6). We conclude that some genes could arise by positioning a copy of the original gene under a promoter in the opposite direction in such a way that both, the original gene and its copy initially use the same nucleotides in the third, degenerated positions in codons.

Asymmetry of coding versus noncoding strand in coding sequences of different genomes.

We have used the asymmetry between the coding and noncoding strands in different codon positions of coding sequences of DNA as a parameter to evaluate the coding probability for open reading frames (ORFs). The method enables an approximation of the total number of coding ORFs in the set of analyzed sequences as well as an estimation of the coding probability for the ORFs. The asymmetry observed in the nucleotide composition of codons in coding sequences has been used successfully for analysis of the genomes completed at the time of this analysis.