Fatty acid biosynthesis in Mycobacterium tuberculosis: lateral gene transfer, adaptive evolution, and gene duplication

Tuberculosis and leprosy in perspective

Two of humankind's most socially and psychologically devastating diseases, tuberculosis and leprosy, have been the subject of intensive paleopathological research due to their antiquity, a presumed association with human settlement and subsistence patterns, and their propensity to leave characteristic lesions on skeletal and mummified remains. Despite a long history of medical research and the development of effective chemotherapy, these diseases remain global health threats even in the 21st century, and as such, their causative agents Mycobacterium tuberculosis and M. leprae, respectively, have recently been the subject of molecular genetics research. The new genome-level data for several mycobacterial species have informed extensive phylogenetic analyses that call into question previously accepted theories concerning the origins and antiquity of these diseases. Of special note is the fact that all new models are in broad agreement that human TB predated that in other animals, including cattle and other domesticates, and that this disease originated at least 35,000 years ago and probably closer to 2.6 million years ago. In this work, we review current phylogenetic and biogeographic models derived from molecular biology and explore their implications for the global development of TB and leprosy, past and present. In so doing, we also briefly review the skeletal evidence for TB and leprosy, explore the current status of these pathogens, critically consider current methods for identifying ancient mycobacterial DNA, and evaluate coevolutionary models.

Recurring cluster and operon assembly for Phenylacetate degradation genes

Background

A large number of theories have been advanced to explain why genes involved in the same biochemical processes are often co-located in genomes. Most of these theories have been dismissed because empirical data do not match the expectations of the models. In this work we test the hypothesis that cluster formation is most likely due to a selective pressure to gradually co-localise protein products and that operon formation is not an inevitable conclusion of the process.

Results

We have selected an exemplar well-characterised biochemical pathway, the phenylacetate degradation pathway, and we show that its complex history is only compatible with a model where a selective advantage accrues from moving genes closer together. This selective pressure is likely to be reasonably weak and only twice in our dataset of 102 genomes do we see independent formation of a complete cluster containing all the catabolic genes in the pathway. Additionally, de novo clustering of genes clearly occurs repeatedly, even though recombination should result in the random dispersal of such genes in their respective genomes. Interspecies gene transfer has frequently replaced in situ copies of genes resulting in clusters that have similar content but very different evolutionary histories.

Conclusion

Our model for cluster formation in prokaryotes, therefore, consists of a two-stage selection process. The first stage is selection to move genes closer together, either because of macromolecular crowding, chromatin relaxation or transcriptional regulation pressure. This proximity opportunity sets up a separate selection for co-transcription.

Horizontally acquired genomic islands in the tubercle bacilli

Most mycobacteria are environmental species, causing disease only occasionally when they encounter a susceptible human or animal host. A few species, such as Mycobacterium tuberculosis and Mycobacterium avium, have acquired the ability to parasitize host macrophages during the course of evolution and have become major pathogens. Recent genetic studies in these two species have suggested that early episodes of horizontal transfer of genomic islands from surrounding environmental species might have contributed to the evolution towards this virulence phenotype, possibly by helping bacilli to persist in protozoa and, subsequently, in mammalian phagocytes. A better understanding of the function of the proteins encoded by these genomic islands in mycobacterial metabolism might help to define novel targets for the development of future antimicrobials.

Insights from the complete genome sequence of Mycobacterium marinum on the evolution of Mycobacterium tuberculosis

Mycobacterium marinum, a ubiquitous pathogen of fish and amphibia, is a near relative of Mycobacterium tuberculosis, the etiologic agent of tuberculosis in humans. The genome of the M strain of M. marinum comprises a 6,636,827-bp circular chromosome with 5424 CDS, 10 prophages, and a 23-kb mercury-resistance plasmid. Prominent features are the very large number of genes (57) encoding polyketide synthases (PKSs) and nonribosomal peptide synthases (NRPSs) and the most extensive repertoire yet reported of the mycobacteria-restricted PE and PPE proteins, and related-ESX secretion systems. Some of the NRPS genes comprise a novel family and seem to have been acquired horizontally. M. marinum is used widely as a model organism to study M. tuberculosis pathogenesis, and genome comparisons confirmed the close genetic relationship between these two species, as they share 3000 orthologs with an average amino acid identity of 85%. Comparisons with the more distantly related Mycobacterium avium subspecies paratuberculosis and Mycobacterium smegmatis reveal how an ancestral generalist mycobacterium evolved into M. tuberculosis and M. marinum. M. tuberculosis has undergone genome downsizing and extensive lateral gene transfer to become a specialized pathogen of humans and other primates without retaining an environmental niche. M. marinum has maintained a large genome so as to retain the capacity for environmental survival while becoming a broad host range pathogen that produces disease strikingly similar to M. tuberculosis. The work described herein provides a foundation for using M. marinum to better understand the determinants of pathogenesis of tuberculosis.

Characterization of a Novel Cell Wall-anchored Protein with Carboxylesterase Activity Required for Virulence in Mycobacterium tuberculosis

Pooled mutant competition assays have shown that the Mycobacterium tuberculosis MT2282 gene (Rv2224c, annotated as encoding a proteinase) is required for bacterial survival in mice. To understand the mechanism of this requirement, we conducted a genetic and biochemical study of the MT2282 gene and its product. MT2282 encodes a member of the microbial esterase/lipase family with active site consensus sequences of G-X-S-X-G, and we have concluded that the MT2282 protein is, in fact, a cell wall-associated carboxylesterase rather than a proteinase, as initially annotated. The MT2282 gene product preferentially hydrolyzes ester bonds of substrates with intermediate carbon chain length. Purified MT2282 is a monomer with enzymatic catalysis properties that fit in the Michaelis-Menten kinetic model. Esterase activity was inhibited by paraoxon and dichlorvos. Replacement of Ser215, Asp450, and His477 by Ala in the consensus motifs completely abolishes esterase activity, suggesting that Ser215-Asp450-His477 forms a catalytic triad with Ser215 as an active site residue. To evaluate the role of the MT2282 in pathogenesis, the gene was deleted from the M. tuberculosis genome. BALB/c mouse aerosol infections showed reduced colony-forming unit loads in lungs and spleens and less lung pathology for the DeltaMT2282 mutant. High dose intravenous infection of mice with the mutant resulted in a significantly delayed time to death compared with the wild type or complemented mutant. These results indicate that MT2282 encodes a cell wall-associated carboxylesterase, which is required for full virulence of M. tuberculosis. We propose that MT2282 (Rv2224c) and its adjacent paralogous gene MT2281 (Rv2223c) be named caeA and caeB respectively, for carboxylesterase A and B.

Comparative genomics of metabolic pathways in Mycobacterium species: gene duplication, gene decay and lateral gene transfer

The genus Mycobacterium comprises significant pathogenic species that infect both humans and animals. One species within this genus, Mycobacterium tuberculosis, is the primary killer of humans resulting from bacterial infections. Five mycobacterial genomes belonging to four different species (M. tuberculosis, Mycobacterium bovis, Mycobacterium leprae and Mycobacterium avium ssp. paratuberculosis) have been sequenced to date and another 14 mycobacterial genomes are at various stages of completion. A comparative analysis of the gene products of key metabolic pathways revealed that the major differences among these species are in the gene products constituting the cell wall and the gene families encoding the acidic glycine-rich (PE/PPE/PGRS) proteins. Mycobacterium leprae has evolved by retaining a minimal gene set for most of the gene families, whereas M. avium ssp. paratuberculosis has acquired some of the virulence factors by lateral gene transfer.

Identification of substrates of the Mycobacterium tuberculosis proteasome

The putative proteasome-associated proteins Mpa (Mycobaterium proteasomal ATPase) and PafA (proteasome accessory factor A) of the human pathogen Mycobacterium tuberculosis (Mtb) are essential for virulence and resistance to nitric oxide. However, a direct link between the proteasome protease and Mpa or PafA has never been demonstrated. Furthermore, protein degradation by bacterial proteasomes in vitro has not been accomplished, possibly due to the failure to find natural degradation substrates or other necessary proteasome co-factors. In this work, we identify the first bacterial proteasome substrates, malonyl Co-A acyl carrier protein transacylase and ketopantoate hydroxymethyltransferase, enzymes that are required for the biosynthesis of fatty acids and polyketides that are essential for the pathogenesis of Mtb. Maintenance of the physiological levels of these enzymes required Mpa and PafA in addition to proteasome protease activity. Mpa levels were also regulated in a proteasome-dependent manner. Finally, we found that a conserved tyrosine of Mpa was essential for function. Thus, these results suggest that Mpa, PafA, and the Mtb proteasome degrade bacterial proteins that are important for virulence in mice.

Lateral gene transfer in Mycobacterium avium subspecies paratuberculosis

Lateral gene transfer is an integral part of genome evolution in most bacteria. Bacteria can readily change the contents of their genomes to increase adaptability to ever-changing surroundings and to generate evolutionary novelty. Here, we report instances of lateral gene transfer in Mycobacterium avium subsp. paratuberculosis, a pathogenic bacteria that causes Johne's disease in cattle. A set of 275 genes are identified that are likely to have been recently acquired by lateral gene transfer. The analysis indicated that 53 of the 275 genes were acquired after the divergence of M. avium subsp. paratuberculosis from M. avium subsp. avium, whereas the remaining 222 genes were possibly acquired by a common ancestor of M. avium subsp. paratuberculosis and M. avium subsp. avium after its divergence from the ancestor of M. tuberculosis complex. Many of the acquired genes were from proteobacteria or soil dwelling actinobacteria. Prominent among the predicted laterally transferred genes is the gene rsbR, a possible regulator of sigma factor, and the genes designated MAP3614 and MAP3757, which are similar to genes in eukaryotes. The results of this study suggest that like most other bacteria, lateral gene transfers seem to be a common feature in M. avium subsp. paratuberculosis and that the proteobacteria contribute most of these genetic exchanges.

Is there an emerging endosymbiotic relationship between mycobacteria and the human host based on horizontal transfer of genetic sequences?

While not negating the seriousness of tuberculosis and the need to prevent and combat the disease effectively, the large percentage of infected, apparently healthy individuals who harbour latent infections warrants consideration whether an endosymbiotic relationship is being established between mycobacteria and man. By means of a gene decay process eliminating their most metabolically important pathogenic genes associated with an increasing need for host gene products during prolonged intracellular survival, mycobacteria appears to be undergoing a process of establishing a less dangerous relationship with its host. To have tolerated this relationship over time, humans must have benefited. This is suggested to have occurred via changes in DNA higher order structure altering combinatorially regulated gene expression allowing increased cerebrodiversity. It can be expected that, beyond a certain threshold, negative effects ensued, leading to neuropathology and increased susceptibility for certain psychiatric disorders. These processes have probably been happening since the earliest contact with mycobacteria, but recently may have become modified by the emergence of epidemic tuberculosis and waves of increased oxidative stress following the circumstances associated with the Industrial Revolution and the more recent AIDS pandemic. The organism seems to have uniquely exploited the normal stress reaction of the host. Genomic stresses include changes associated with glucocorticoid effects as well as upregulated reactive oxygen species and stress/(heat shock) protein production, the latter two of which result in host cell cycle delay. Subsequently replication dependent chromosomal fragile sites appear in the host genome and together with upregulated chaperonins and mobile element activation, the scene is set for sequence exchange between the organism and host. If proven, these events raise the possibility of modifying chromatin epigenetically to retain the proposed advantages while silencing pathogenicity factors.

Horizontal Transfer of a Virulence Operon to the Ancestor of Mycobacterium tuberculosis

The contribution of interspecies horizontal gene transfer (HGT) to the evolution and virulence of Mycobacterium tuberculosis, the agent of tuberculosis in humans, has been barely investigated. Here we have studied the evolutionary history of the M. tuberculosis Rv0986-8 virulence operon recently identified, through functional genomics approaches, as playing an important role in parasitism of host phagocytic cells. We showed that among actinobacteria, this operon is specific to the M. tuberculosis complex and to ancestral Mycobacterium prototuberculosis species. These data, together with phylogenetic reconstruction and other in silico analyses, provided strong evidence that this operon has been acquired horizontally by the ancestor of M. tuberculosis, before the recent evolutionary bottleneck that preceded the clonal-like evolution of the M. tuberculosis complex. Genomic signature profiling further suggested that the transfer was plasmid mediated and that the operon originated from a gamma-proteobacterium donor species. Our study points out for the first time the contribution of HGT to the emergence of M. tuberculosis and close relatives as major pathogens. In addition, our data underline the importance of deciphering gene transfer networks in M. tuberculosis in order to better understand the evolutionary mechanisms involved in mycobacterial virulence.

Mycobacterium tuberculosis complex genetic diversity: mining the fourth international spoligotyping database (SpolDB4) for classification, population genetics and epidemiology

Background

The Direct Repeat locus of the Mycobacterium tuberculosis complex (MTC) is a member of the CRISPR (Clustered regularly interspaced short palindromic repeats) sequences family. Spoligotyping is the widely used PCR-based reverse-hybridization blotting technique that assays the genetic diversity of this locus and is useful both for clinical laboratory, molecular epidemiology, evolutionary and population genetics. It is easy, robust, cheap, and produces highly diverse portable numerical results, as the result of the combination of (1) Unique Events Polymorphism (UEP) (2) Insertion-Sequence-mediated genetic recombination. Genetic convergence, although rare, was also previously demonstrated. Three previous international spoligotype databases had partly revealed the global and local geographical structures of MTC bacilli populations, however, there was a need for the release of a new, more representative and extended, international spoligotyping database.

Results

The fourth international spoligotyping database, SpolDB4, describes 1939 shared-types (STs) representative of a total of 39,295 strains from 122 countries, which are tentatively classified into 62 clades/lineages using a mixed expert-based and bioinformatical approach. The SpolDB4 update adds 26 new potentially phylogeographically-specific MTC genotype families. It provides a clearer picture of the current MTC genomes diversity as well as on the relationships between the genetic attributes investigated (spoligotypes) and the infra-species classification and evolutionary history of the species. Indeed, an independent Naïve-Bayes mixture-model analysis has validated main of the previous supervised SpolDB3 classification results, confirming the usefulness of both supervised and unsupervised models as an approach to understand MTC population structure. Updated results on the epidemiological status of spoligotypes, as well as genetic prevalence maps on six main lineages are also shown. Our results suggests the existence of fine geographical genetic clines within MTC populations, that could mirror the passed and present Homo sapiens sapiens demographical and mycobacterial co-evolutionary history whose structure could be further reconstructed and modelled, thereby providing a large-scale conceptual framework of the global TB Epidemiologic Network.

Conclusion

Our results broaden the knowledge of the global phylogeography of the MTC complex. SpolDB4 should be a very useful tool to better define the identity of a given MTC clinical isolate, and to better analyze the links between its current spreading and previous evolutionary history. The building and mining of extended MTC polymorphic genetic databases is in progress.

Testing of a whole genome PCR scanning approach to identify genomic variability in four different species of lactic acid bacteria

Genomes can be markedly heterogeneous in conspecific bacterial strains. Genome sequences can be used to analyze genome plasticity via a PCR(2) (plasticity of chromosome revealed by PCR) approach. Small-sized chromosomes can indeed be fully amplified by long-range PCR with a set of primers designed using a reference strain and then applied to several other strains. Analysis of the resulting patterns can reveal genome plasticity. GenoFrag, a software package for the design of primers optimized for PCR(2) [N. Ben Zakour, M. Gautier, R. Andonov, D. Lavenier, M.F. Cochet, P. Veber, A. Sorokin, Y. Le Loir, GenoFrag: Software to design primers optimized for whole genome scanning by long-range PCR amplification, Nucleic Acids Res. 32 (2004) 17-24] was developed for the analysis of bacterial genome plasticity by whole genome amplification in approximately 10-kb-long fragments. By applying GenoFrag, we provide herewith evidence that genome plasticity can be analyzed in lactic acid bacteria using a PCR(2) approach. The genome sequences of Lactococcus lactis IL1403, Lactobacillus plantarum WCFS1, Lactobacillus bulgaricus ATCC11842 and Bifidobacterium longum NCC2705 were used to design four sets of primers. Each set was evaluated in silico to check that it ensured optimum coverage of the bacterial chromosome. To validate the primers generated by GenoFrag, a subset of primers was successfully used in LR-PCR experiments on genomic DNA from four L. bulgaricus strains.

Genome phylogenies indicate a meaningful alpha-proteobacterial phylogeny and support a grouping of the mitochondria with the Rickettsiales

Placement of the mitochondrial branch on the tree of life has been problematic. Sparse sampling, the uncertainty of how lateral gene transfer might overwrite phylogenetic signals, and the uncertainty of phylogenetic inference have all contributed to the issue. Here we address this issue using a supertree approach and completed genomic sequences. We first determine that a sensible alpha-proteobacterial phylogenetic tree exists and that it can confidently be inferred using orthologous genes. We show that congruence across these orthologous gene trees is significantly better than might be expected by random chance. There is some evidence of horizontal gene transfer within the alpha-proteobacteria, but it appears to be restricted to a minority of genes ( approximately 23%) most of whom ( approximately 74%) can be categorized as operational. This means that placement of the mitochondrion should not be excessively hampered by interspecies gene transfer. We then show that there is a consistently strong signal for placement of the mitochondrion on this tree and that this placement is relatively insensitive to methodological approach or data set. A concatenated alignment was created consisting of 15 mitochondrion-encoded proteins that are unlikely to have undergone any lateral gene transfer in the timeline under consideration. This alignment infers that the sister group of the mitochondria, for the taxa that have been sampled, is the order Rickettsiales.

Evidence of positive Darwinian selection in putative meningococcal vaccine antigens

Meningococcal meningitidis is a life-threatening disease. In Europe and the United States the majority of cases are caused by virulent meningococcal strains belonging to serogroup B. Presently there is no effective vaccine against serogroup B strains, as traditional vaccine antigens such as polysaccharide capsules are unusable as they lead to autoimmunity. The year 2000 saw the publication of the complete genome of Neisseria meningitidis MC58, a virulent serogroup B bacterium. Working in conjunction with the sequencing project, researchers endeavored to locate highly conserved membrane-associated proteins that elicit an immune response. It is hoped that these proteins will provide a basis for novel vaccines against serogroup B strains. A number of potential vaccine antigens have been located and are presently in phase I clinical trials. Recently many reports pertaining to the evidence of positive Darwinian selection in membrane proteins of pathogens have been reported. This study utilized in silico methods to test for evidence of historical positive Darwinian selection in seven such vaccine candidates. We found that two of these proteins show signatures of adaptive evolution, while the remaining proteins show evidence of strong purifying selection. This has significant implications for the design of a vaccine against serogroup B strains, as it has been shown that vaccines that target epitopes that are under strong purifying selection are better than those that target variable epitopes.

Evidence of positive Darwinian selection in Omp85, a highly conserved bacterial outer membrane protein essential for cell viability

Omp85 is a highly conserved outer membrane protein found in all gram-negative bacteria. It is essential for bacterial cell viability and plays an integral function in the positioning and folding of other outer membrane proteins into the bacterial outer membrane. We have employed a maximum likelihood and a maximum parsimony approach to detect evidence of positive Darwinian selection in Omp85 homologues from 10 delta-proteobacteria and have identified 14 amino acid sites that show evidence of being under the influence of adaptive evolution. Interestingly all sites bar one are concentrated within surface loops of the protein that most likely interact with host immune response or the surrounding environment. Alternatively amino acids within membrane-spanning regions of the protein are found to be under purifying selection most likely as a result of structural constraints.

The comparative method rules! Codon volatility cannot detect positive Darwinian selection using a single genome sequence

All established methods for detecting positive selection at the molecular level rely on comparisons between nucleotide sequences. An exceptional method that purports to detect selection on the basis of a single genomic sequence has recently been proposed. This method uses a measure called "codon volatility," defined for each codon as the ratio between the number of nonsynonymous codons that differ from the codon under study at a single nucleotide position and the number of sense codons that differ from the codon under study at a single nucleotide position. Here, we examine various properties of codon volatility and its derivatives and use simulation of evolutionary processes to determine whether they can be used to detect selective pressures. Codons for only four amino acids (glycine, leucine, arginine, and serine) show any variation in codon volatility. Thus, codon volatility is mainly a proxy for amino acid usage, rather than for codon usage, with 65% of all synonymous changes and 27% of all nonsynonymous changes being undetectable by this measure. Genes identified by the volatility method as being subject to positive selection tend to have idiosyncratic amino acid compositions (e.g., they are glycine rich or arginine poor). An additional property of codon volatility is the near zero variance of its mean expectation, which translates into overestimated statistical significance estimates, especially in the absence of corrections for multiple comparisons. A comparison with measures of selection inferred through comparative methodology reveals no relationship between the results of the two methods. Finally, we show that codon volatility can increase in the absence of positive Darwinian selection; that is, increased codon volatility is not indicative of positive selection.

Naturally occurring horizontal gene transfer and homologous recombination in Mycobacterium

Acquisition of genetic information through horizontal gene transfer (HGT) is an important evolutionary process by which micro-organisms gain novel phenotypic characteristics. In pathogenic bacteria, for example, it facilitates maintenance and enhancement of virulence and spread of drug resistance. In the genus Mycobacterium, to which several primary human pathogens belong, HGT has not been clearly demonstrated. The few existing reports suggesting this process are based on circumstantial evidence of similarity of sequences found in distantly related species. Here, direct evidence of HGT between strains of Mycobacterium avium representing two different serotypes is presented. Conflicting evolutionary histories of genes encoding elements of the glycopeptidolipid (GPL) biosynthesis pathway led to an analysis of the GPL cluster genomic sequences from four Mycobacterium avium strains. The sequence of M. avium strain 2151 appeared to be a mosaic consisting of three regions having alternating identities to either M. avium strains 724 or 104. Maximum-likelihood estimation of two breakpoints allowed a approximately 4100 bp region horizontally transferred into the strain 2151 genome to be pinpointed with confidence. The maintenance of sequence continuity at both breakpoints and the lack of insertional elements at these sites strongly suggest that the integration of foreign DNA occurred by homologous recombination. To our knowledge, this is the first report to demonstrate naturally occurring homologous recombination in Mycobacterium. This previously undiscovered mechanism of genetic exchange may have major implications for the understanding of Mycobacterium pathogenesis.

Overexpression of isocitrate lyase is an important strategy in the survival of Pseudomonas fluorescens exposed to aluminum

Isocitrate lyase, ICL (EC 4.1.3.1), an enzyme that cleaves isocitrate into succinate, and glyoxylate appears to play a pivotal role in the detoxification of aluminum (Al) in Pseudomonas fluorescens. Here, we present evidence that the 4-fold increase in ICL activity observed in Al-stressed cells is due to the overexpression of this enzyme. Blue-Native-PAGE, Western blotting, and spectrophotometric experiments revealed that ICL is optimally expressed at 35 h of growth in Al-stressed cells. However, following the immobilization of Al, at 60 h of growth, the level of the enzyme decreases markedly. This enzyme that exists as a homotetramer with a molecular mass of approximately 133 kDa appears to be transcriptionally regulated. The overexpression of ICL may be a specific response to Al-stress as P. fluorescens grown in the presence of such metals as Ga3+, Pb2+, and Ca2+ does not undergo any significant increase in ICL activity. Thus, these findings support the notion that the overexpression of ICL plays a pivotal role in the survival and in the increased oxalogenesis observed in Al-stressed P. fluorescens.