Genomewide pattern of synonymous nucleotide substitution in two complete genomes of Mycobacterium tuberculosis

The paradoxes of Mycobacterium tuberculosis molecular evolution and consequences for the inference of tuberculosis emergence date

The date of Mycobacterium tuberculosis complex emergence has been the subject of long debates. New studies joining archaeological efforts with sequencing methods raise high hopes for solving whether this emergence is closer to 70,000 or to 6000 years before present. Inferring the date of emergence of this pathogen based on sequence data requires a molecular clock. Several clocks inferred from different types of loci and/or different samples, using both sound reasoning and reliable data, are actually very different, which we refer to as the paradoxes of M. tuberculosis molecular evolution. After having presented these paradoxes, we will remind the limits of the molecular clocks used in the different studies such as the assumption of homogeneous substitution rate. We will then review recent results that shed new light on the characteristics of M. tuberculosis molecular evolution: traces of diverse selection pressures, the impact of host dynamics, etc. We provide some ideas on what to do next to get nearer to a reliable dating of Mycobacterium tuberculosis complex emergence. Among them, the collection of additional remains from ancient tuberculosis seems still essential.

Macrophage activation highlight an important role for NER proteins in the survival, latency and multiplication of Mycobacterium tuberculosis

Nucleotide excision repair (NER) is one of the most extensively studied DNA repair processes in both prokaryotes and eukaryotes. The NER pathway is a highly conserved, ATP-dependent multi-step process involving several proteins/enzymes that function in a concerted manner to recognize and excise a wide spectrum of helix-distorting DNA lesions and bulky adducts by nuclease cleavage on either side of the damaged bases. As such, the NER pathway of Mycobacterium tuberculosis (Mtb) is essential for its survival within the hostile environment of macrophages and disease progression. This review focuses on present published knowledge about the crucial roles of Mtb NER proteins in the survival and multiplication of the pathogen within the macrophages and as potential targets for drug discovery.

Genome-wide estimation of recombination, mutation and positive selection enlightens diversification drivers of Mycobacterium bovis

Genome sequencing has reinvigorated the infectious disease research field, shedding light on disease epidemiology, pathogenesis, host-pathogen interactions and also evolutionary processes exerted upon pathogens. Mycobacterium tuberculosis complex (MTBC), enclosing M. bovis as one of its animal-adapted members causing tuberculosis (TB) in terrestrial mammals, is a paradigmatic model of bacterial evolution. As other MTBC members, M. bovis is postulated as a strictly clonal, slowly evolving pathogen, with apparently no signs of recombination or horizontal gene transfer. In this work, we applied comparative genomics to a whole genome sequence (WGS) dataset composed by 70 M. bovis from different lineages (European and African) to gain insights into the evolutionary forces that shape genetic diversification in M. bovis. Three distinct approaches were used to estimate signs of recombination. Globally, a small number of recombinant events was identified and confirmed by two independent methods with solid support. Still, recombination reveals a weaker effect on M. bovis diversity compared with mutation (overall r/m = 0.037). The differential r/m average values obtained across the clonal complexes of M. bovis in our dataset are consistent with the general notion that the extent of recombination may vary widely among lineages assigned to the same taxonomical species. Based on this work, recombination in M. bovis cannot be excluded and should thus be a topic of further effort in future comparative genomics studies for which WGS of large datasets from different epidemiological scenarios across the world is crucial. A smaller M. bovis dataset (n = 42) from a multi-host TB endemic scenario was then subjected to additional analyses, with the identification of more than 1,800 sites wherein at least one strain showed a single nucleotide polymorphism (SNP). The majority (87.1%) was located in coding regions, with the global ratio of non-synonymous upon synonymous alterations (dN/dS) exceeding 1.5, suggesting that positive selection is an important evolutionary force exerted upon M. bovis. A higher percentage of SNPs was detected in genes enriched into "lipid metabolism", "cell wall and cell processes" and "intermediary metabolism and respiration" functional categories, revealing their underlying importance in M. bovis biology and evolution. A closer look on genes prone to horizontal gene transfer in the MTBC ancestor and included in the 3R (DNA repair, replication and recombination) system revealed a global average negative value for Taijima's D neutrality test, suggesting that past selective sweeps and population expansion after a recent bottleneck remain as major evolutionary drivers of the obligatory pathogen M. bovis in its struggle with the host.

Mycobacterium tuberculosis complex and human coadaptation: a two-way street complicating host susceptibility to TB

For centuries, the Mycobacterium tuberculosis complex (MTBC) has infected numerous populations, both human and non-human, causing symptomatic tuberculosis (TB) in some hosts. Research investigating the MTBC and how it has evolved with its host over time is sparse and has not resulted in many significant findings. There are even fewer studies investigating adaptation of the human host susceptibility to TB and these have largely focused on genome-wide association and candidate gene association studies. However, results emanating from these association studies are rarely replicated and appear to be population specific. It is, therefore, necessary to relook at the approach taken to investigate the relationship between the MTBC and the human host. Understanding that the evolution of the pathogen is coupled to the evolution of the host might be the missing link needed to effectively investigate their relationship. We hypothesize that this knowledge will bolster future efforts in combating the disease.

Comparative analysis of corrected tiger genome provides clues to its neuronal evolution

The availability of completed and draft genome assemblies of tiger, leopard, and other felids provides an opportunity to gain comparative insights on their unique evolutionary adaptations. However, genome-wide comparative analyses are susceptible to errors in genome sequences and thus require accurate genome assemblies for reliable evolutionary insights. In this study, while analyzing the tiger genome, we found almost one million erroneous substitutions in the coding and non-coding region of the genome affecting 4,472 genes, hence, biasing the current understanding of tiger evolution. Moreover, these errors produced several misleading observations in previous studies. Thus, to gain insights into the tiger evolution, we corrected the erroneous bases in the genome assembly and gene set of tiger using ‘SeqBug’ approach developed in this study. We sequenced the first Bengal tiger genome and transcriptome from India to validate these corrections. A comprehensive evolutionary analysis was performed using 10,920 orthologs from nine mammalian species including the corrected gene sets of tiger and leopard and using five different methods at three hierarchical levels, i.e. felids, Panthera, and tiger. The unique genetic changes in tiger revealed that the genes showing signatures of adaptation in tiger were enriched in development and neuronal functioning. Specifically, the genes belonging to the Notch signalling pathway, which is among the most conserved pathways involved in embryonic and neuronal development, were found to have significantly diverged in tiger in comparison to the other mammals. Our findings suggest the role of adaptive evolution in neuronal functions and development processes, which correlates well with the presence of exceptional traits such as sensory perception, strong neuro-muscular coordination, and hypercarnivorous behaviour in tiger.

Mycobacterial Evolution Intersects With Host Tolerance

Over the past 200 years, tuberculosis (TB) has caused more deaths than any other infectious disease, likely infecting more people than it has at any other time in human history. Mycobacterium tuberculosis (Mtb), the etiologic agent of TB, is an obligate human pathogen that has evolved through the millennia to become an archetypal human-adapted pathogen. This review focuses on the evolutionary framework by which Mtb emerged as a specialized human pathogen and applies this perspective to the emergence of specific lineages that drive global TB burden. We consider how evolutionary pressures, including transmission dynamics, host tolerance, and human population patterns, may have shaped the evolution of diverse mycobacterial genomes.

Turkish and Japanese Mycobacterium tuberculosis sublineages share a remote common ancestor

Two geographically distant M. tuberculosis sublineages, Tur from Turkey and T3-Osaka from Japan, exhibit partially identical genotypic signatures (identical 12-loci MIRU-VNTR profiles, distinct spoligotyping patterns). We investigated T3-Osaka and Tur sublineages characteristics and potential genetic relatedness, first using MIRU-VNTR locus analysis on 21 and 25 samples of each sublineage respectively, and second comparing Whole Genome Sequences of 8 new samples to public data from 45 samples uncovering human tuberculosis diversity. We then tried to date their Most Recent Common Ancestor (MRCA) using three calibrations of SNP accumulation rate (long-term=0.03SNP/genome/year, derived from a tuberculosis ancestor of around 70,000years old; intermediate=0.2SNP/genome/year derived from a Peruvian mummy; short-term=0.5SNP/genome/year). To disentangle between these scenarios, we confronted the corresponding divergence times with major human history events and knowledge on human genetic divergence. We identified relatively high intrasublineage diversity for both T3-Osaka and Tur. We definitively proved their monophyly; the corresponding super-sublineage (referred to as "T3-Osa-Tur") shares a common ancestor with T3-Ethiopia and Ural sublineages but is only remotely related to other Euro-American sublineages such as X, LAM, Haarlem and S. The evolutionary scenario based on long-term evolution rate being valid until T3-Osa-Tur MRCA was not supported by Japanese fossil data. The evolutionary scenario relying on short-term evolution rate since T3-Osa-Tur MRCA was contradicted by human history and potential traces of past epidemics. T3-Osaka and Tur sublineages were found likely to have diverged between 800y and 2000years ago, potentially at the time of Mongol Empire. Altogether, this study definitively proves a strong genetic link between Turkish and Japanese tuberculosis. It provides a first hypothesis for calibrating TB Euro-American lineage molecular clock; additional studies are needed to reliably date events corresponding to intermediate depths in tuberculosis phylogeny.

Co-evolution of Mycobacterium tuberculosis and Homo sapiens

The causative agent of human tuberculosis (TB), Mycobacterium tuberculosis, is an obligate pathogen that evolved to exclusively persist in human populations. For M. tuberculosis to transmit from person to person, it has to cause pulmonary disease. Therefore, M. tuberculosis virulence has likely been a significant determinant of the association between M. tuberculosis and humans. Indeed, the evolutionary success of some M. tuberculosis genotypes seems at least partially attributable to their increased virulence. The latter possibly evolved as a consequence of human demographic expansions. If co-evolution occurred, humans would have counteracted to minimize the deleterious effects of M. tuberculosis virulence. The fact that human resistance to infection has a strong genetic basis is a likely consequence of such a counter-response. The genetic architecture underlying human resistance to M. tuberculosis remains largely elusive. However, interactions between human genetic polymorphisms and M. tuberculosis genotypes have been reported. Such interactions are consistent with local adaptation and allow for a better understanding of protective immunity in TB. Future 'genome-to-genome' studies, in which locally associated human and M. tuberculosis genotypes are interrogated in conjunction, will help identify new protective antigens for the development of better TB vaccines.

Genotyping of ancient Mycobacterium tuberculosis strains reveals historic genetic diversity

The evolutionary history of the Mycobacterium tuberculosis complex (MTBC) has previously been studied by analysis of sequence diversity in extant strains, but not addressed by direct examination of strain genotypes in archaeological remains. Here, we use ancient DNA sequencing to type 11 single nucleotide polymorphisms and two large sequence polymorphisms in the MTBC strains present in 10 archaeological samples from skeletons from Britain and Europe dating to the second–nineteenth centuries AD. The results enable us to assign the strains to groupings and lineages recognized in the extant MTBC. We show that at least during the eighteenth–nineteenth centuries AD, strains of M. tuberculosis belonging to different genetic groups were present in Britain at the same time, possibly even at a single location, and we present evidence for a mixed infection in at least one individual. Our study shows that ancient DNA typing applied to multiple samples can provide sufficiently detailed information to contribute to both archaeological and evolutionary knowledge of the history of tuberculosis.

Evolution of smooth tubercle Bacilli PE and PE_PGRS genes: evidence for a prominent role of recombination and imprint of positive selection

BACKGROUND

PE and PE_PGRS are two mycobateria-restricted multigene families encoding membrane associated and secreted proteins that have expanded mainly in the pathogenic species, notably the Mycobacterium tuberculosis complex (MTBC). Several lines of evidence attribute to PE and PE_PGRS genes critical roles in mycobacterial pathogenicity. To get more insight into the nature of these genes, we sought to address their evolutionary trajectories in the group of smooth tubercle bacilli (STB), the putative ancestor of the clonal MTBC.

METHODOLOGY/PRINCIPAL FINDINGS

By focussing on six polymorphic STB PE/PE_PGRS genes, we demonstrate significant incongruence among single gene genealogies and detect strong signals of recombination using various approaches. Coalescent-based estimation of population recombination and mutation rates (ρ and θ, respectively) indicates that the two mechanisms are of roughly equal importance in generating diversity (ρ/θ = 1.457), a finding in a marked contrast to house keeping genes (HKG) whose evolution is chiefly brought about by mutation (ρ/θ = 0.012). In comparison to HKG, we found 15 times higher mean rate of nonsynonymous substitutions, with strong evidence of positive selection acting on PE_PGRS62 (dN/dS = 1.42), a gene that has previously been shown to be essential for mycobacterial survival in macrophages and granulomas. Imprint of positive selection operating on specific amino acid residues or along branches of PE_PGRS62 phylogenetic tree was further demonstrated using maximum likelihood- and covarion-based approaches, respectively. Strikingly, PE_PGR62 proved highly conserved in present-day MTBC strains.

CONCLUSIONS/SIGNIFICANCE

Overall the data indicate that, in STB, PE/PE_PGRS genes have undergone a strong diversification process that is speeded up by recombination, with evidence of positive selection. The finding that positive selection involved an essential PE_PGRS gene whose sequence appears to be driven to fixation in present-day MTBC strains lends further support to the critical role of PE/PE_PGRS genes in the evolution of mycobacterial pathogenicity.

Genetic diversity in Mycobacterium tuberculosis

Recent years have witnessed an increased appreciation of the extent and relevance of strain-to-strain variation in Mycobacterium tuberculosis. This paradigm shift can largely be attributed to an improved understanding of the global population structure of this organism, and to the realisation that the various members of the M. tuberculosis complex (MTBC) harbour more genetic diversity than previously realised. Moreover, many studies using experimental models of infection have demonstrated that MTBC diversity translates into significant differences in immunogenecity and virulence . However, linking these experimental phenotypes to relevant clinical phenotypes has been difficult, and to date, largely unsuccessful. Nevertheless, emerging high-throughput technologies, in particular next-generation sequencing , offer new opportunities, and have already lead to important new insights. Given the complexity of the host-pathogen interaction in tuberculosis, systems approaches will be key to define the role of MTBC diversity in the fight against one of humankind's most important pathogens.

Comparative analysis of Mycobacterium tuberculosis pe and ppe genes reveals high sequence variation and an apparent absence of selective constraints

Mycobacterium tuberculosis complex (MTBC) genomes contain 2 large gene families termed pe and ppe. The function of pe/ppe proteins remains enigmatic but studies suggest that they are secreted or cell surface associated and are involved in bacterial virulence. Previous studies have also shown that some pe/ppe genes are polymorphic, a finding that suggests involvement in antigenic variation. Using comparative sequence analysis of 18 publicly available MTBC whole genome sequences, we have performed alignments of 33 pe (excluding pe_pgrs) and 66 ppe genes in order to detect the frequency and nature of genetic variation. This work has been supplemented by whole gene sequencing of 14 pe/ppe (including 5 pe_pgrs) genes in a cohort of 40 diverse and well defined clinical isolates covering all the main lineages of the M. tuberculosis phylogenetic tree. We show that nsSNP's in pe (excluding pgrs) and ppe genes are 3.0 and 3.3 times higher than in non-pe/ppe genes respectively and that numerous other mutation types are also present at a high frequency. It has previously been shown that non-pe/ppe M. tuberculosis genes display a remarkably low level of purifying selection. Here, we also show that compared to these genes those of the pe/ppe families show a further reduction of selection pressure that suggests neutral evolution. This is inconsistent with the positive selection pressure of "classical" antigenic variation. Finally, by analyzing such a large number of genes we were able to detect large differences in mutation type and frequency between both individual genes and gene sub-families. The high variation rates and absence of selective constraints provides valuable insights into potential pe/ppe function. Since pe/ppe proteins are highly antigenic and have been studied as potential vaccine components these results should also prove informative for aspects of M. tuberculosis vaccine design.

Bacterial genomics in infectious disease and the clinical pathology laboratory

CONTEXT

Throughout history, technologic advancements have fueled the engine of innovation, which, in turn, has driven discovery. Accordingly, recent advancements in DNA sequencing technology are revolutionizing bacterial genomics.

OBJECTIVE

To review important developments from the literature. The current state of bacterial genomics, with an emphasis on human pathogens and the clinical pathology laboratory, will be discussed.

DATA SOURCES

A comprehensive review was performed of the relevant literature indexed in PubMed (National Library of Medicine) and referenced medical texts.

CONCLUSIONS

Many important discoveries bearing on infectious disease research and pathology laboratory practice have been achieved through whole-genome sequencing strategies. Bacterial genomics has improved our understanding of molecular pathogenesis, host-pathogen interactions, and antibiotic-resistance mechanisms. Bacterial genomics has also facilitated the study of population structures, epidemics and outbreaks, and newly identified pathogens. Many opportunities now exist for clinical pathologists to contribute to bacterial genomics, including in the design of new diagnostic tests, therapeutic agents, and vaccines.

After the bottleneck: Genome-wide diversification of the Mycobacterium tuberculosis complex by mutation, recombination, and natural selection

Many of the most virulent bacterial pathogens show low genetic diversity and sexual isolation. Accordingly, Mycobacterium tuberculosis, the deadliest human pathogen, is thought to be clonal and evolve by genetic drift. Yet, its genome shows few of the concomitant signs of genome degradation. We analyzed 24 genomes and found an excess of genetic diversity in regions encoding key adaptive functions including the type VII secretion system and the ancient horizontally transferred virulence-related regions. Four different approaches showed evident signs of recombination in M. tuberculosis. Recombination tracts add a high density of polymorphisms, and many are thus predicted to arise from outside the clade. Some of these tracts match Mycobacterium canettii sequences. Recombination introduced an excess of non-synonymous diversity in general and even more in genes expected to be under positive or diversifying selection, e.g., cell wall component genes. Mutations leading to non-synonymous SNPs are effectively purged in MTBC, which shows dominance of purifying selection. MTBC mutation bias toward AT nucleotides is not compensated by biased gene conversion, suggesting the action of natural selection also on synonymous changes. Together, all of these observations point to a strong imprint of recombination and selection in the genome affecting both non-synonymous and synonymous positions. Hence, contrary to some other pathogens and previous proposals concerning M. tuberculosis, this lineage may have come out of its ancestral bottleneck as a very successful pathogen that is rapidly diversifying by the action of mutation, recombination, and natural selection.

Independent large scale duplications in multiple M. tuberculosis lineages overlapping the same genomic region

Mycobacterium tuberculosis, the causative agent of most human tuberculosis, infects one third of the world's population and kills an estimated 1.7 million people a year. With the world-wide emergence of drug resistance, and the finding of more functional genetic diversity than previously expected, there is a renewed interest in understanding the forces driving genome evolution of this important pathogen. Genetic diversity in M. tuberculosis is dominated by single nucleotide polymorphisms and small scale gene deletion, with little or no evidence for large scale genome rearrangements seen in other bacteria. Recently, a single report described a large scale genome duplication that was suggested to be specific to the Beijing lineage. We report here multiple independent large-scale duplications of the same genomic region of M. tuberculosis detected through whole-genome sequencing. The duplications occur in strains belonging to both M. tuberculosis lineage 2 and 4, and are thus not limited to Beijing strains. The duplications occur in both drug-resistant and drug susceptible strains. The duplicated regions also have substantially different boundaries in different strains, indicating different originating duplication events. We further identify a smaller segmental duplication of a different genomic region of a lab strain of H37Rv. The presence of multiple independent duplications of the same genomic region suggests either instability in this region, a selective advantage conferred by the duplication, or both. The identified duplications suggest that large-scale gene duplication may be more common in M. tuberculosis than previously considered.

DNA fingerprinting of Mycobacterium tuberculosis: from phage typing to whole-genome sequencing

Current typing methods for Mycobacterium tuberculosis complex evolved from simple phenotypic approaches like phage typing and drug susceptibility profiling to DNA-based strain typing methods, such as IS6110-restriction fragment length polymorphisms (RFLP) and variable number of tandem repeats (VNTR) typing. Examples of the usefulness of molecular typing are source case finding and epidemiological linkage of tuberculosis (TB) cases, international transmission of MDR/XDR-TB, the discrimination between endogenous reactivation and exogenous re-infection as a cause of relapses after curative treatment of tuberculosis, the evidence of multiple M. tuberculosis infections, and the disclosure of laboratory cross-contaminations. Simultaneously, phylogenetic analyses were developed based on single nucleotide polymorphisms (SNPs), genomic deletions usually referred to as regions of difference (RDs) and spoligotyping which served both strain typing and phylogenetic analysis. National and international initiatives that rely on the application of these typing methods have brought significant insight into the molecular epidemiology of tuberculosis. However, current DNA fingerprinting methods have important limitations. They can often not distinguish between genetically closely related strains and the turn-over of these markers is variable. Moreover, the suitability of most DNA typing methods for phylogenetic reconstruction is limited as they show a high propensity of convergent evolution or misinfer genetic distances. In order to fully explore the possibilities of genotyping in the molecular epidemiology of tuberculosis and to study the phylogeny of the causative bacteria reliably, the application of whole-genome sequencing (WGS) analysis for all M. tuberculosis isolates is the optimal, although currently still a costly solution. In the last years WGS for typing of pathogens has been explored and yielded important additional information on strain diversity in comparison to the classical DNA typing methods. With the ongoing cost reduction of DNA sequencing it is possible that WGS will become the sole diagnostic tool in the secondary laboratory diagnosis of tuberculosis for identification, drug susceptibility testing and genetic characterization.

Single nucleotide polymorphism in the genes of mce1 and mce4 operons of Mycobacterium tuberculosis: analysis of clinical isolates and standard reference strains

Background

The presence of four mammalian cell entry (mce) operons in Mycobacterium tuberculosis suggests the essentiality of the functions of the genes in these operons. The differential expression of the four mce operons in different phases of in vitro growth and in infected animals reported earlier from our laboratory further justifies the apparent redundancy for these genes in the genome.

Here we investigate the extent of polymorphism in eight genes in the mce1 and mce4 operons of M. tuberculosis from four standard reference strains (H37Rv, H37Ra, LVS (Low Virulent Strain) and BCG) and 112 clinical isolates varying in their drug susceptibility profile, analysed by direct sequencing and Sequenom MassARRAY platform.

Results

We discovered 20 single nucleotide polymorphisms (SNPs) in the two operons. The comparative analysis of the genes of mce1 and mce4 operons revealed that yrbE1A [Rv0167] was most polymorphic in mce1 operon while yrbE4A [Rv3501c] and lprN [Rv3495c] had the highest number of SNPs in the mce4 operon. Of 20 SNPs, 12 were found to be nonsynonymous and were further analysed for their pathological relevance to M. tuberculosis using web servers PolyPhen and PMut, which predicted five deleterious nonsynonymous SNPs. A mutation from proline to serine at position 359 of the native Mce1A protein was most deleterious as predicted by both PolyPhen and PMut servers. Energy minimization of the structure of native Mce1A protein and mutated protein was performed using InsightII. The mutated Mce1A protein showed structural changes that could account for the effects of this mutation.

Conclusions

Our results show that SNPs in the coding sequences of mce1 and mce4 operons in clinical isolates can be significantly high. Moreover, mce4 operon is significantly more polymorphic than mce1 operon (p < 0.001). However, the frequency of nonsynonymous substitutions is higher in mce1 operon and synonymous substitutions are more in mce4 operon. In silico modeling predict that nonsynonymous SNP at mce1A [Rv0169], a virulence gene could play a pivotal role in causing functional changes in M. tuberculosis that may reflect upon the biology of the bacteria.

DNA repair systems and the pathogenesis of Mycobacterium tuberculosis: varying activities at different stages of infection

Mycobacteria, including most of all MTB (Mycobacterium tuberculosis), cause pathogenic infections in humans and, during the infectious process, are exposed to a range of environmental insults, including the host's immune response. From the moment MTB is exhaled by infected individuals, through an active and latent phase in the body of the new host, until the time they reach the reactivation stage, MTB is exposed to many types of DNA-damaging agents. Like all cellular organisms, MTB has efficient DNA repair systems, and these are believed to play essential roles in mycobacterial pathogenesis. As different stages of infection have great variation in the conditions in which mycobacteria reside, it is possible that different repair systems are essential for progression to specific phases of infection. MTB possesses homologues of DNA repair systems that are found widely in other species of bacteria, such as nucleotide excision repair, base excision repair and repair by homologous recombination. MTB also possesses a system for non-homologous end-joining of DNA breaks, which appears to be widespread in prokaryotes, although its presence is sporadic within different species within a genus. However, MTB does not possess homologues of the typical mismatch repair system that is found in most bacteria. Recent studies have demonstrated that DNA repair genes are expressed differentially at each stage of infection. In the present review, we focus on different DNA repair systems from mycobacteria and identify questions that remain in our understanding of how these systems have an impact upon the infection processes of these important pathogens.

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.

The evolutionary biology of poxviruses

The poxviruses (family Poxviridae) are a family of double-stranded viruses including several species that infect humans and their domestic animals, most notably Variola virus (VARV), the causative agent of smallpox. The evolutionary biology of these viruses poses numerous questions, for which we have only partial answers at present. Here we review evidence regarding the origin of poxviruses, the frequency of host transfer in poxvirus history, horizontal transfer of host genes to poxviruses, and the population processes accounting for patterns of nucleotide sequence polymorphism.

Genomic diversity and evolution of Mycobacterium ulcerans revealed by next-generation sequencing

Mycobacterium ulcerans is the causative agent of Buruli ulcer, the third most common mycobacterial disease after tuberculosis and leprosy. It is an emerging infectious disease that afflicts mainly children and youths in West Africa. Little is known about the evolution and transmission mode of M. ulcerans, partially due to the lack of known genetic polymorphisms among isolates, limiting the application of genetic epidemiology. To systematically profile single nucleotide polymorphisms (SNPs), we sequenced the genomes of three M. ulcerans strains using 454 and Solexa technologies. Comparison with the reference genome of the Ghanaian classical lineage isolate Agy99 revealed 26,564 SNPs in a Japanese strain representing the ancestral lineage. Only 173 SNPs were found when comparing Agy99 with two other Ghanaian isolates, which belong to the two other types previously distinguished in Ghana by variable number tandem repeat typing. We further analyzed a collection of Ghanaian strains using the SNPs discovered. With 68 SNP loci, we were able to differentiate 54 strains into 13 distinct SNP haplotypes. The average SNP nucleotide diversity was low (average 0.06–0.09 across 68 SNP loci), and 96% of the SNP locus pairs were in complete linkage disequilibrium. We estimated that the divergence of the M. ulcerans Ghanaian clade from the Japanese strain occurred 394 to 529 thousand years ago. The Ghanaian subtypes diverged about 1000 to 3000 years ago, or even much more recently, because we found evidence that they evolved significantly faster than average. Our results offer significant insight into the evolution of M. ulcerans and provide a comprehensive report on genetic diversity within a highly clonal M. ulcerans population from a Buruli ulcer endemic region, which can facilitate further epidemiological studies of this pathogen through the development of high-resolution tools.

Mycobacterium ulcerans is the causative agent of Buruli ulcer (BU), a necrotizing skin disease and the third most common mycobacterial disease after tuberculosis and leprosy. It is an emerging infectious disease that afflicts mainly children and youths in West Africa. The disease is also found in tropical and subtropical regions of Asia, the Western Pacific, and Latin America. Limited knowledge of this neglected tropical disease is partially due to the lack of known genetic polymorphisms among isolates, which hinder the study of transmission, epidemiology, and evolution of M. ulcerans. Our aim is to systematically profile genetic diversity among M. ulcerans isolates by sequencing and comparing the genomes of selected strains. We identified single nucleotide polymorphisms (SNPs) within a highly clonal M. ulcerans population from a Buruli ulcer endemic region. Based on the SNPs discovered, we developed SNP typing assays and were able to differentiate a collection of M. ulcerans isolates from this Buruli ulcer endemic region into 13 SNP haplotypes. Our results lay the ground for developing a highly discriminatory and cost-effective tool to study M. ulcerans evolution and epidemiology at a population level.

Diet, tuberculosis, and the paleopathological record

Osseous manifestation of infectious disease is of paramount importance to paleopathologists seeking to interpret ancient health, but the relationships among infectious agent exposure, development of disease, and skeletal involvement are complex. The outcome of an exposure strongly depends on multiple factors, including ecology, diet, nutrition, immune function, and the genetics of pathogen and host. Mycobacterial diseases are often studied in ancient remains but also are especially influenced by these factors; individual and population differences in severity and course are apparent following onset of active disease. The osteological record for these diseases represents the complex interplay of host and pathogen characteristics influencing within- and among-individual skeletal lesion prevalence and distribution. However, many of these characteristics may be assessed independently through the archaeological record. Here, we explore the contributions of dietary protein and iron to immune function, particularly the course and outcome of infection with Mycobacterium tuberculosis. We emphasize how nutrition may influence the dissemination of bacilli to the skeleton and subsequent formation of diagnostic lesions. We then generate models and hypotheses informed by this interplay and apply them to four prehistoric New World areas. Finally, discrepancies between our expectations and the observed record are explored as a basis for new hypotheses.

Genomic comparison of PE and PPE genes in the Mycobacterium avium complex

The Mycobacterium avium complex (MAC) comprises genomically similar but phenotypically divergent bacteria that inhabit diverse environments and that cause disease in different hosts. In this study, a whole-genome approach was used to examine the polymorphic PE (Pro-Glu) and PPE (Pro-Pro-Glu) gene families, implicated in immunostimulation and virulence. The four major groups of MAC organisms were examined, including the newly sequenced type strains of M. intracellulare and M. avium subsp. avium, plus M. avium subsp. paratuberculosis and M. avium subsp. hominissuis, for the purpose of finding genetic differences that could be exploited to design diagnostic tests specific to these groups and that could help explain their divergence in pathogenesis and host specificity. Unique and missing PPE genes were found in all MAC members except M. avium subsp. avium. Only M. intracellulare had a unique PE gene. Apart from this, most PE and PPE sequences were conserved, with average nucleotide sequence identities of 99.1 and 98.1%, respectively, among the M. avium subspecies, but only 82.9 and 79.7% identities with the PE and PPE sequences of M. intracellulare, respectively. A detailed analysis of the amino acid sequences was performed between M. avium subsp. paratuberculosis and M. avium subsp. hominissuis. Most differences were detected in the PPE proteins, with amino acid substitutions and frame shifts leading to unique amino acid sequences. In conclusion, several unique PPE proteins were identified in MAC organisms next to numerous polymorphisms in both the PE and PPE gene families. These substantial differences could help explain the divergence in phenotypes within the MAC and could lead to diagnostic tests with better discriminatory abilities.

High functional diversity in Mycobacterium tuberculosis driven by genetic drift and human demography

Mycobacterium tuberculosis infects one third of the human world population and kills someone every 15 seconds. For more than a century, scientists and clinicians have been distinguishing between the human- and animal-adapted members of the M. tuberculosis complex (MTBC). However, all human-adapted strains of MTBC have traditionally been considered to be essentially identical. We surveyed sequence diversity within a global collection of strains belonging to MTBC using seven megabase pairs of DNA sequence data. We show that the members of MTBC affecting humans are more genetically diverse than generally assumed, and that this diversity can be linked to human demographic and migratory events. We further demonstrate that these organisms are under extremely reduced purifying selection and that, as a result of increased genetic drift, much of this genetic diversity is likely to have functional consequences. Our findings suggest that the current increases in human population, urbanization, and global travel, combined with the population genetic characteristics of M. tuberculosis described here, could contribute to the emergence and spread of drug-resistant tuberculosis.

Detection and molecular characterization of 9,000-year-old Mycobacterium tuberculosis from a Neolithic settlement in the Eastern Mediterranean

Background

Mycobacterium tuberculosis is the principal etiologic agent of human tuberculosis. It has no environmental reservoir and is believed to have co-evolved with its host over millennia. This is supported by skeletal evidence of the disease in early humans, and inferred from M. tuberculosis genomic analysis. Direct examination of ancient human remains for M. tuberculosis biomarkers should aid our understanding of the nature of prehistoric tuberculosis and the host/pathogen relationship.

Methodology/Principal Findings

We used conventional PCR to examine bone samples with typical tuberculosis lesions from a woman and infant, who were buried together in the now submerged site of Atlit-Yam in the Eastern Mediterranean, dating from 9250-8160 years ago. Rigorous precautions were taken to prevent contamination, and independent centers were used to confirm authenticity of findings. DNA from five M tuberculosis genetic loci was detected and had characteristics consistent with extant genetic lineages. High performance liquid chromatography was used as an independent method of verification and it directly detected mycolic acid lipid biomarkers, specific for the M. tuberculosis complex.

Conclusions/Significance

Human tuberculosis was confirmed by morphological and molecular methods in a population living in one of the first villages with evidence of agriculture and animal domestication. The widespread use of animals was not a source of infection but may have supported a denser human population that facilitated transmission of the tubercle bacillus. The similarity of the M. tuberculosis genetic signature with those of today gives support to the theory of a long-term co-existence of host and pathogen.

Origin, spread and demography of the Mycobacterium tuberculosis complex

The evolutionary timing and spread of the Mycobacterium tuberculosis complex (MTBC), one of the most successful groups of bacterial pathogens, remains largely unknown. Here, using mycobacterial tandem repeat sequences as genetic markers, we show that the MTBC consists of two independent clades, one composed exclusively of M. tuberculosis lineages from humans and the other composed of both animal and human isolates. The latter also likely derived from a human pathogenic lineage, supporting the hypothesis of an original human host. Using Bayesian statistics and experimental data on the variability of the mycobacterial markers in infected patients, we estimated the age of the MTBC at 40,000 years, coinciding with the expansion of "modern" human populations out of Africa. Furthermore, coalescence analysis revealed a strong and recent demographic expansion in almost all M. tuberculosis lineages, which coincides with the human population explosion over the last two centuries. These findings thus unveil the dynamic dimension of the association between human host and pathogen populations.

Sequence polymorphisms in a surface PPE protein distinguish types I, II, and III of Mycobacterium avium subsp. paratuberculosis

In the last 2 decades, a variety of different molecular typing methods have been developed to differentiate strains of Mycobacterium avium subsp. paratuberculosis. The most successful techniques are based on insertion sequences, repetitive loci, comparative genomics, or single nucleotide polymorphisms. In the present study, we chose to examine whether a single M. avium subsp. paratuberculosis gene could serve as a means of differentiation of a variety of isolates. The MAP1506 gene locus encodes a member of the polymorphic PPE protein family that has putative roles relevant to M. avium subsp. paratuberculosis pathogenicity. The MAP1506 locus was sequenced from a collection of 58 M. avium subsp. paratuberculosis isolates from different sources, hosts, and typing profiles. Following sequence alignment and analysis, it was found that bovine (type II) strains of M. avium subsp. paratuberculosis consistently differed from ovine (type I) and intermediate (type III) strains in seven and eight nucleotides, respectively. Polymorphic regions of the MAP1506 locus were selected for analysis by denaturing gradient gel electrophoresis, allowing visual discrimination of the three subtypes of M. avium subsp. paratuberculosis isolates. This is the first report describing the use of PCR and denaturing gradient gel electrophoresis on a single gene as a method to distinguish types I, II, and III of M. avium subsp. paratuberculosis.

Genomics and the evolution, pathogenesis, and diagnosis of tuberculosis

Tuberculosis kills nearly 2 million people annually, and current approaches to tuberculosis control are expensive, have limited efficacy, and are vulnerable to being overcome by extensively drug-resistant strains of Mycobacterium tuberculosis. Determination of the genome sequence of M. tuberculosis has revolutionized tuberculosis research, contributed to major advances in the understanding of the evolution and pathogenesis of M. tuberculosis, and facilitated development of new diagnostic tests with increased specificity for tuberculosis. In this review, we describe some of the major progress in tuberculosis research that has resulted from knowledge of the genome sequence and note some of the problems that remain unsolved.

The role of IS6110 in the evolution of Mycobacterium tuberculosis

Members of the Mycobacterium tuberculosis complex contain the transposable element IS6110 which, due to its high numerical and positional polymorphism, has become a widely used marker in epidemiological studies. Here, we review the evidence that IS6110 is not simply a passive or 'junk' DNA sequence, but that, through its transposable activity, it is able to generate genotypic variation that translates into strain-specific phenotypic variation. We also speculate on the role that this variation has played in the evolution of M. tuberculosis and conclude that the presence of a moderate IS6110 copy number within the genome may provide the pathogen with a selective advantage that has aided its virulence.

Global phylogeography of Mycobacterium tuberculosis and implications for tuberculosis product development

New tools for controlling tuberculosis are urgently needed. Despite our emerging understanding of the biogeography of Mycobacterium tuberculosis, the implications for development of new diagnostics, drugs, and vaccines is unknown. M tuberculosis has a clonal genetic population structure that is geographically constrained. Evidence suggests strain-specific differences in virulence and immunogenicity in light of this global phylogeography. We propose a strain selection framework, based on robust phylogenetic markers, which will allow for systematic and comprehensive evaluation of new tools for tuberculosis control.

Evaluating candidate agents of selective pressure for cystic fibrosis

Cystic fibrosis is the most common lethal single-gene mutation in people of European descent, with a carrier frequency upwards of 2%. Based upon molecular research, resistances in the heterozygote to cholera and typhoid fever have been proposed to explain the persistence of the mutation. Using a population genetic model parameterized with historical demographic and epidemiological data, we show that neither cholera nor typhoid fever provided enough historical selective pressure to produce the modern incidence of cystic fibrosis. However, we demonstrate that the European tuberculosis pandemic beginning in the seventeenth century would have provided sufficient historical, geographically appropriate selective pressure under conservative assumptions. Tuberculosis has been underappreciated as a possible selective agent in producing cystic fibrosis but has clinical, molecular and now historical, geographical and epidemiological support. Implications for the future trajectory of cystic fibrosis are discussed. Our result supports the importance of novel investigations into the role of arylsulphatase B deficiency in cystic fibrosis and tuberculosis.

Homologous recombination and the pattern of nucleotide substitution in Ehrlichia ruminantium

Patterns of nucleotide substitution at orthologous loci were examined between three genomes of Ehrlichia ruminantium, the causative agent of heartwater disease of ruminants. The most recent common ancestor of two genomes (Erwe and Erwo) belonging to the Welgevonden strain was estimated to have occurred 26,500-57,000 years ago, while the most recent common ancestor of these two genomes and the Erga genome (Gardel strain) was estimated to have occurred 2.1-4.7 million years ago. The search for genes showing extremely high values of the number of synonymous substitutions per site was used to identify genes involved in past homologous recombination. The most striking case involved the map1 gene, encoding major antigenic protein-1; evidence for homologous recombination is consistent with previous phylogenetic analysis of map1 alleles. At this and certain other loci, homologous recombination may have contributed to the evolution of host-pathogen interactions. In addition, comparison of the patterns of synonymous and nonsynonymous substitution provided evidence for positive selection favoring a high level of amino acid change between the Welgevonden and Gardel strains at a locus of unknown function (designated Erum4340 in the Erwo genome).

Population-based study of deletions in five different genomic regions of Mycobacterium tuberculosis and possible clinical relevance of the deletions

Regions of difference (RDs) have been described in clinical isolates of Mycobacterium tuberculosis, but the potential epidemiological and clinical relevance of the genotypes of these RDs remains to be investigated. We screened a population-based sample of 648 isolates for the deletion of five RDs, designated RD105, RD181, RD142, RD150, and RD239, using microarray-based hybridization, PCR, and DNA sequencing and assessed the associations between the RD deletions and the clinical characteristics of the patients using chi-square analysis and multivariate logistic regression model. Of the 648 isolates, 18 (2.8%) had the RD239 deletion and 39 (6.0%) had the RD105 deletion. The deletions of RD142, RD150, and RD181 subdivided the isolates with the RD105 deletion into four groups comprising a group with concurrent deletions of RD105, RD181, and RD142 (n = 13); a group with concurrent deletions of RD105, RD181, and RD150 (n = 5); a group with concurrent deletions of RD105 and RD181 (n = 13); and a group with a deletion of RD105 only (n = 8). Extrathoracic tuberculosis is statistically significantly associated with infection with the isolates with concurrent deletions of RD105, RD181, and RD142 (adjusted odds ratio [OR] = 3.05; 95% confidence interval [CI] = 1.58, 5.90) and the isolates with concurrent deletions of RD105, RD181, and RD150 (adjusted OR = 11.09; 95% CI = 4.27, 28.80), after controlling for the previously identified risk factors for extrathoracic tuberculosis (human immunodeficiency virus serostatus, race, gender, and the genotype of the plcD gene). These two combinations of RD deletions have the potential for predicting the clinical presentation of M. tuberculosis infection in the human host.

Bottlenecks and broomsticks: the molecular evolution of Mycobacterium bovis

Mycobacterium bovis is the cause of tuberculosis in cattle and is a member of the Mycobacterium tuberculosis complex. In contrast to many other pathogenic bacterial species, there is little evidence for the transfer and recombination of genes between cells. The clonality of this group of organisms indicates that the population structure is dominated by reductions in diversity, caused either by population bottlenecks or selective sweeps as entire chromosomes become fixed in the population. We describe how these forces have shaped not only the phylogeny of this group but also, at a very local level, the population structure of Mycobacterium bovis in the British Isles. We also discuss the practical implications of applying this knowledge to understanding the spread of infection and the development of improved vaccines and diagnostic tests.

Genomic diversity in Mycobacterium avium: single nucleotide polymorphisms between the S and C strains of M. avium subsp. paratuberculosis and with M. a. avium

Following identification of large genomic polymorphisms in a previous study, a polymerase chain reaction and sequencing strategy was used to identify single nucleotide polymorphisms (SNPs) in 25 genes in the sheep (S) and cattle (C) strains of Mycobacterium avium subsp. paratuberculosis (M. a. paratuberculosis) and between M. a. paratuberculosis and M. a. avium. From 12,117 bp of sequence representing 26 loci across 25 genes, 11 SNPs were identified between the S and C strains in eight genes: hsp65, sodA, dnaA, dnaN, recF, gyrB, inhA, and pks8. An in silico comparison of these M. a. paratuberculosis sequences and the M. a. avium 104 genome revealed 86 SNPs, which corresponded well with similar studies of SNPs in the M. avium complex.

Novel genetic polymorphisms that further delineate the phylogeny of the Mycobacterium tuberculosis complex

In a previous report, we described a PCR protocol for the differentiation of the various species of the Mycobacterium tuberculosis complex (MTC) on the basis of genomic deletions (R. C. Huard, L. C. de Oliveira Lazzarini, W. R. Butler, D. van Soolingen, and J. L. Ho, J. Clin. Microbiol. 41:1637-1650, 2003). That report also provided a broad cross-comparison of several previously identified, phylogenetically relevant, long-sequence and single-nucleotide polymorphisms (LSPs and SNPs, respectively). In the present companion report, we expand upon the previous work (i) by continuing the evaluation of known MTC phylogenetic markers in a larger collection of tubercle bacilli (n = 125), (ii) by evaluating additional recently reported MTC species-specific and interspecific polymorphisms, and (iii) by describing the identification and distribution of a number of novel LSPs and SNPs. Notably, new genomic deletions were found in various Mycobacterium tuberculosis strains, new species-specific SNPs were identified for "Mycobacterium canettii," Mycobacterium microti, and Mycobacterium pinnipedii, and, for the first time, intraspecific single-nucleotide DNA differences were discovered for the dassie bacillus, the oryx bacillus, and the two Mycobacterium africanum subtype I variants. Surprisingly, coincident polymorphisms linked one M. africanum subtype I genotype with the dassie bacillus and M. microti with M. pinnipedii, thereby suggesting closer evolutionary ties within each pair of species than had been previously thought. Overall, the presented data add to the genetic definitions of several MTC organisms as well as fine-tune current models for the evolutionary history of the MTC.

Probing genomic diversity and evolution of Escherichia coli O157 by single nucleotide polymorphisms

Infections by Shiga toxin-producing Escherichia coli O157:H7 (STEC O157) are the predominant cause of bloody diarrhea and hemolytic uremic syndrome in the United States. In silico comparison of the two complete STEC O157 genomes (Sakai and EDL933) revealed a strikingly high level of sequence identity in orthologous protein-coding genes, limiting the use of nucleotide sequences to study the evolution and epidemiology of this bacterial pathogen. To systematically examine single nucleotide polymorphisms (SNPs) at a genome scale, we designed comparative genome sequencing microarrays and analyzed 1199 chromosomal genes (a total of 1,167,948 bp) and 92,721 bp of the large virulence plasmid (pO157) of eleven outbreak-associated STEC O157 strains. We discovered 906 SNPs in 523 chromosomal genes and observed a high level of DNA polymorphisms among the pO157 plasmids. Based on a uniform rate of synonymous substitution for Escherichia coli and Salmonella enterica (4.7x10(-9) per site per year), we estimate that the most recent common ancestor of the contemporary beta-glucuronidase-negative, non-sorbitolfermenting STEC O157 strains existed ca. 40 thousand years ago. The phylogeny of the STEC O157 strains based on the informative synonymous SNPs was compared to the maximum parsimony trees inferred from pulsed-field gel electrophoresis and multilocus variable numbers of tandem repeats analysis. The topological discrepancies indicate that, in contrast to the synonymous mutations, parts of STEC O157 genomes have evolved through different mechanisms with highly variable divergence rates. The SNP loci reported here will provide useful genetic markers for developing high-throughput methods for fine-resolution genotyping of STEC O157. Functional characterization of nucleotide polymorphisms should shed new insights on the evolution, epidemiology, and pathogenesis of STEC O157 and related pathogens.

Global phylogeny of Mycobacterium tuberculosis based on single nucleotide polymorphism (SNP) analysis: insights into tuberculosis evolution, phylogenetic accuracy of other DNA fingerprinting systems, and recommendations for a minimal standard SNP set

We analyzed a global collection of Mycobacterium tuberculosis strains using 212 single nucleotide polymorphism (SNP) markers. SNP nucleotide diversity was high (average across all SNPs, 0.19), and 96% of the SNP locus pairs were in complete linkage disequilibrium. Cluster analyses identified six deeply branching, phylogenetically distinct SNP cluster groups (SCGs) and five subgroups. The SCGs were strongly associated with the geographical origin of the M. tuberculosis samples and the birthplace of the human hosts. The most ancestral cluster (SCG-1) predominated in patients from the Indian subcontinent, while SCG-1 and another ancestral cluster (SCG-2) predominated in patients from East Asia, suggesting that M. tuberculosis first arose in the Indian subcontinent and spread worldwide through East Asia. Restricted SCG diversity and the prevalence of less ancestral SCGs in indigenous populations in Uganda and Mexico suggested a more recent introduction of M. tuberculosis into these regions. The East African Indian and Beijing spoligotypes were concordant with SCG-1 and SCG-2, respectively; X and Central Asian spoligotypes were also associated with one SCG or subgroup combination. Other clades had less consistent associations with SCGs. Mycobacterial interspersed repetitive unit (MIRU) analysis provided less robust phylogenetic information, and only 6 of the 12 MIRU microsatellite loci were highly differentiated between SCGs as measured by GST. Finally, an algorithm was devised to identify two minimal sets of either 45 or 6 SNPs that could be used in future investigations to enable global collaborations for studies on evolution, strain differentiation, and biological differences of M. tuberculosis.

Ancient origin and gene mosaicism of the progenitor of Mycobacterium tuberculosis

The highly successful human pathogen Mycobacterium tuberculosis has an extremely low level of genetic variation, which suggests that the entire population resulted from clonal expansion following an evolutionary bottleneck around 35,000 y ago. Here, we show that this population constitutes just the visible tip of a much broader progenitor species, whose extant representatives are human isolates of tubercle bacilli from East Africa. In these isolates, we detected incongruence among gene phylogenies as well as mosaic gene sequences, whose individual elements are retrieved in classical M. tuberculosis. Therefore, despite its apparent homogeneity, the M. tuberculosis genome appears to be a composite assembly resulting from horizontal gene transfer events predating clonal expansion. The amount of synonymous nucleotide variation in housekeeping genes suggests that tubercle bacilli were contemporaneous with early hominids in East Africa, and have thus been coevolving with their human host much longer than previously thought. These results open novel perspectives for unraveling the molecular bases of M. tuberculosis evolutionary success.

Mycobacterium tuberculosis, the agent of tuberculosis, is a highly successful human pathogen and kills nearly 3 million persons each year. This pathogen and its close relatives sum up in a single and compact clonal group dating back only a few tens of thousands of years. Using genetic data, the researchers have discovered that human tubercle bacilli from East Africa represent extant bacteria of a much broader progenitor species from which the M. tuberculosis clonal group evolved. They estimate that this progenitor species is as old as 3 million years. This suggests that our remote hominid ancestors may well have already suffered from tuberculosis. In addition, the researchers show that tubercle bacilli are able to exchange parts of their genome with other strains, a process that is known to play a crucial role in adaptation of pathogens to their hosts. Thus, the M. tuberculosis genome appears to be a composite assembly, resulting from ancient horizontal DNA exchanges before its clonal expansion. These findings open novel perspectives for unraveling the origin and the molecular bases of M. tuberculosis evolutionary success, and lead to reconsideration of the impact of tuberculosis on human natural selection.

Nucleotide substitution and recombination at orthologous loci in Staphylococcus aureus

The pattern of nucleotide substitution was examined at 2,129 orthologous loci among five genomes of Staphylococcus aureus, which included two sister pairs of closely related genomes (MW2/MSSA476 and Mu50/N315) and the more distantly related MRSA252. A total of 108 loci were unusual in lacking any synonymous differences among the five genomes; most of these were short genes encoding proteins highly conserved at the amino acid sequence level (including many ribosomal proteins) or unknown predicted genes. In contrast, 45 genes were identified that showed anomalously high divergence at synonymous sites. The latter genes were evidently introduced by homologous recombination from distantly related genomes, and in many cases, the pattern of nucleotide substitution made it possible to reconstruct the most probable recombination event involved. These recombination events introduced genes encoding proteins that differed in amino acid sequence and thus potentially in function. Several of the proteins are known or likely to be involved in pathogenesis (e.g., staphylocoagulase, exotoxin, Ser-Asp fibrinogen-binding bone sialoprotein-binding protein, fibrinogen and keratin-10 binding surface-anchored protein, fibrinogen-binding protein ClfA, and enterotoxin P). Therefore, the results support the hypothesis that exchange of homologous genes among S. aureus genomes can play a role in the evolution of pathogenesis in this species.

Patterns of sequence divergence in 5' intergenic spacers and linked coding regions in 10 species of pathogenic bacteria reveal distinct recombinational histories

We compared the pattern of nucleotide difference in 8034 genes and in their 5' intergenic spacers between conspecific pairs of genomes from 10 species of pathogenic bacteria. Certain genes or spacers showed much greater sequence divergence between the genotypes compared to others; such divergent regions plausibly originated by recombinational events by which a gene and/or spacers was donated from a divergent genome. Different patterns of divergence in genes and spacers identified different recombinational patterns. For example, in Chlamydophila pneumoniae, there were examples of both unusually divergent spacers and unusually divergent genes, but there were no cases in which a gene and its spacer were both unusually divergent. This pattern suggests that, in C. pneumoniae, recombination events have broken up the linkage between genes and 5' spacers. By contrast, in Streptococcus agalactiae, there were a number of cases in which both spacer and gene were unusually divergent, indicating that a number of large-scale recombination events that included both genes and 5' spacers have occurred; there was evidence of at least two large-scale recombination events in the genomic region including the pur genes in S. agalactiae.

Size of the protein-coding genome and rate of molecular evolution

In diploid populations of size N, there will be 2 Nmu mutations per nucleotide (nt) site (or per locus) per generation (mu stands for mutation rate). If either the population or the coding genome double in size, one expects 4 Nmu mutations. What is important is not the population size per se but the number of genes (coding sites), the two being often interconverted. Here we compared the total physical length of protein-coding genomes (n) with the corresponding absolute rates of synonymous substitution (K(S)), an empirical neutral reference. In the classical occupancy problem and in the coupons collector (CC) problem, n was expressed as the mean rate of change (K(CC)). Despite inherently very low power of the approaches involving averaging of rates, the mode of molecular evolution of the total size phenotype of the coding genome could be evidenced through differences between the genomic estimates of K(CC) [K(CC)=1/(ln n + 0.57721) n] and rate of molecular evolution, K(S). We found that (1) the estimates of n and K(S) are reciprocally correlated across taxa (r=0.812; p<< 0.001); (2) the gamete-cell division hypothesis (Chang et al. Proc Natl Acad Sci USA 91:827-831, 1994) can be confirmed independently in terms of K(CC)/K(S) ratios; (3) the time scale of molecular evolution changes with change in mutation rate, as previously shown by Takahata (Proc Natl Acad Sci USA 87:2419-2423, 1990), Takahata et al. (Genetics 130:925-938, 1992), and Vekemans and Slatkin (Genetics 137:1157-1165, 1994); (4) the generation time and population size (Lynch and Conery, Science 302:1401-1404, 2003) effects left their "signatures" at the level of the size phenotype of the protein-coding genome.

Effects of genetic variability of Mycobacterium tuberculosis strains on the presentation of disease

The nature of the variability in the clinical and epidemiological consequences of Mycobacterium tuberculosis infection remains poorly understood. Environmental and host factors that contribute to the outcome of infection and disease presentation are well recognised, but the role of bacterial factors has been more elusive. The rapid increase in the understanding of the molecular basis of M tuberculosis over the past decades has revived research into its pathogenesis. DNA fingerprinting techniques have been used to distinguish between strains of M tuberculosis, and efforts to characterise the strains present within populations have led to increased understanding of their global distribution. This research has shown that in certain areas a small number of strains are causing a disproportionate number of cases of the disease. The sequencing of the complete genome of M tuberculosis has accelerated the development of molecular techniques to differentiate strains according to their genetic polymorphisms. Investigation into the reasons why some strains are predominant by genetic strain-typing techniques may clarify which bacterial factors contribute to disease. This knowledge has the potential to influence control and prevention strategies for tuberculosis in the future. However, there are still limitations in these techniques and their results. This review discusses molecular epidemiology and genetic studies, and their contribution to the understanding of the links between genotypic and phenotypic characteristics of M tuberculosis strains.

Common evolutionary origin for the unstable virulence plasmid pMUM found in geographically diverse strains of Mycobacterium ulcerans

The 174-kb virulence plasmid pMUM001 in Mycobacterium ulcerans epidemic strain Agy99 harbors three very large and homologous genes that encode giant polyketide synthases (PKS) responsible for the synthesis of the lipid toxin mycolactone. Deeper investigation of M. ulcerans Agy99 resulted in identification of two types of spontaneous deletion variants of pMUM001 within a population of cells that also contained the intact plasmid. These variants arose from recombination between two 8-kb sections of the same plasmid sequence, resulting in the loss of a 65-kb region bearing two of the three mycolactone PKS genes. Investigation of nine diverse M. ulcerans strains by using PCR and Southern hybridization for eight pMUM001 gene sequences confirmed the presence of pMUM001-like elements (collectively called pMUM) in all M. ulcerans strains. Physical mapping of these plasmids revealed that like M. ulcerans Agy99, three strains had undergone major deletions in their mycolactone PKS loci. Online liquid chromatography-sequential mass spectrometry analysis of lipid extracts confirmed that strains with PKS deletions were unable to produce mycolactone or any related cometabolites. Interstrain comparisons of the plasmid gene sequences revealed greater than 98% nucleotide identity, and the phylogeny inferred from these sequences closely mimicked the phylogeny from a previous multilocus sequence typing study in which chromosomally encoded loci were used, a result that is consistent with the hypothesis that M. ulcerans diverged from the closely related organism Mycobacterium marinum by acquiring pMUM. Our results suggest that pMUM is a defining characteristic of M. ulcerans but that in the absence of purifying selection, deletion of plasmid sequences and a corresponding loss of mycolactone production readily arise.