Paweł Mackiewicz

Chromalveolate plastids: direct descent or multiple endosymbioses?

We are grateful to Katja Bargum for valuable comments on the manuscript. A.B. is supported by funds from Wroclaw University grant BS/1018/2008.

The relationships between the isoelectric point and: length of proteins, taxonomy and ecology of organisms

BackgroundThe distribution of isoelectric point (pI) of proteins in a proteome is universal for all organisms. It is bimodal dividing the proteome into two sets of acidic and basic proteins. Different species however have different abundance of acidic and basic proteins that may be correlated with taxonomy, subcellular localization, ecological niche of organisms and proteome size.ResultsWe have analysed 1784 proteomes encoded by chromosomes of Archaea, Bacteria, Eukaryota, and also mitochondria, plastids, prokaryotic plasmids, phages and viruses. We have found significant correlation in more than 95% of proteomes between the protein length and pI in proteomes – positive for acidic proteins and negative for the basic ones. Plastids, viruses and plasmids encode more basic proteomes while chromosomes of Archaea, Bacteria, Eukaryota, mitochondria and phages more acidic ones. Mitochondrial proteomes of Viridiplantae, Protista and Fungi are more basic than Metazoa. It results from the presence of basic proteins in the former proteomes and their absence from the latter ones and is related with reduction of metazoan genomes. Significant correlation was found between the pI bias of proteomes encoded by prokaryotic chromosomes and proteomes encoded by plasmids but there is no correlation between eukaryotic nuclear-coded proteomes and proteomes encoded by organelles. Detailed analyses of prokaryotic proteomes showed significant relationships between pI distribution and habitat, relation to the host cell and salinity of the environment, but no significant correlation with oxygen and temperature requirements. The salinity is positively correlated with acidicity of proteomes. Host-associated organisms and especially intracellular species have more basic proteomes than free-living ones. The higher rate of mutations accumulation in the intracellular parasites and endosymbionts is responsible for the basicity of their tiny proteomes that explains the observed positive correlation between the decrease of genome size and the increase of basicity of proteomes. The results indicate that even conserved proteins subjected to strong selectional constraints follow the global trend in the pI distribution.ConclusionThe distribution of pI of proteins in proteomes shows clear relationships with length of proteins, subcellular localization, taxonomy and ecology of organisms. The distribution is also strongly affected by mutational pressure especially in intracellular organisms.

Flip-flop around the origin and terminus of replication in prokaryotic genomes

A response to Evidence for symmetric chromosomal inversions around the replication origin in bacteria by JA Eisen, JF Heidelberg, O White, SL Salzberg. Genome Biology 2000, 1:research0011.1-0011.9.

Early steps in plastid evolution: current ideas and controversies

Some nuclear‐encoded proteins are imported into higher plant plastids via the endomembrane (EM) system. Compared with multi‐protein Toc and Tic translocons required for most plastid protein import, the relatively uncomplicated nature of EM trafficking led to suggestions that it was the original transport mechanism for nuclear‐encoded endosymbiont proteins, and critical for the early stages of plastid evolution. Its apparent simplicity disappears, however, when EM transport is considered in light of selective constraints likely encountered during the conversion of stable endosymbionts into fully integrated organelles. From this perspective it is more parsimonious to presume the early evolution of post‐translational protein import via simpler, ancestral forms of modern Toc and Tic plastid translocons, with EM trafficking arising later to accommodate glycosylation and/or protein targeting to multiple cellular locations. This hypothesis is supported by both empirical and comparative data, and is consistent with the relative paucity of EM‐based transport to modern primary plastids.

A fragment of chloroplast DNA was transferred horizontally, probably from non-eudicots, to mitochondrial genome of Phaseolus

The mitochondrial genomes of some Phaseolus species contain a fragment of chloroplast trnA gene intron, named pvs-trnA for its location within the Phaseolus vulgaris sterility sequence (pvs). The purpose of this study was to determine the type of transfer (intracellular or horizontal) that gave rise to pvs-trnA. Using a PCR approach we could not find the respective portion of the trnA gene as a part of pvs outside the Phaseolus genus. However, a BLAST search revealed longer fragments of trnA present in the mitochondrial genomes of some Citrus species, Helianthus annuus and Zea mays. Basing on the identity or near-identity between these mitochondrial sequences and their chloroplast counterparts we concluded that they had relocated from chloroplasts to mitochondria via recent, independent, intracellular DNA transfers. In contrast, pvs-trnA displayed a relatively higher sequence divergence when compared with its chloroplast counterpart from Phaseolus vulgaris. Alignment of pvs-trnA with corresponding trnAfragments from 35 plant species as well as phylogenetic analysis revealed that pvs-trnA grouped with non-eudicot sequences and was well separated from all Fabalessequences. In conclusion, we propose that pvs-trnA arose via horizontal transfer of a trnA intron fragment from chloroplast of a non-eudicot plant to Phaseolus mitochondria. This is the first example of horizontal transfer of a chloroplast sequence to the mitochondrial genome in higher plants.

Isolation and characterisation of KP34—a novel φKMV-like bacteriophage for Klebsiella pneumoniae

Bacteriophage KP34 is a novel virus belonging to the subfamily Autographivirinae lytic for extended-spectrum β-lactamase-producing Klebsiella pneumoniae strains. Its biological features, morphology, susceptibility to chemical and physical agents, burst size, host specificity and activity spectrum were determined. As a potential antibacterial agent used in therapy, KP34 molecular features including genome sequence and protein composition were examined. Phylogenetic analyses and clustering of KP34 phage genome sequences revealed its clear relationships with “phiKMV-like viruses”. Simultaneously, whole-genome analyses permitted clustering and classification of all phages, with completely sequenced genomes, belonging to the Podoviridae.

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.

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

BackgroundAny 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.ResultsWe 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.ConclusionsParameters 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.

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. Copyright

Evolution Rates of Genes on Leading and Lagging DNA Strands

Abstract. 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.