Naturally occurring horizontal gene transfer and homologous recombination in Mycobacterium

Unveiling metabolic characteristics of an uncultured Gammaproteobacterium responsible for in situ PAH biodegradation in petroleum polluted soil

Exploring the metabolic characteristics of indigenous PAH degraders is critical to understanding the PAH bioremediation mechanism in the natural environment. While stable-isotopic probing (SIP) is a viable method to identify functional microorganisms in complex environments, the metabolic characteristics of uncultured degraders are still elusive. Here, we investigated the naphthalene (NAP) biodegradation of petroleum polluted soils by combining SIP, amplicon sequencing and metagenome binning. Based on the SIP and amplicon sequencing results, an uncultured Gammaproteobacterium sp. was identified as the key NAP degrader. Additionally, the assembled genome of this uncultured degrader was successfully obtained from the ¹³ C-DNA metagenomes by matching its 16S rRNA gene with the SIP identified OTU sequence. Meanwhile, a number of NAP degrading genes encoding naphthalene/PAH dioxygenases were identified in this genome, further confirming the direct involvement of this indigenous degrader in the NAP degradation. The degrader contained genes related to the metabolisms of several carbon sources, energy substances and vitamins, illuminating potential reasons for why microorganisms cannot be cultivated and finally realize their cultivation. Our findings provide novel information on the mechanisms of in situ PAH biodegradation and add to our current knowledge on the cultivation of non-culturable microorganisms by combining both SIP and metagenome binning.

Direct PCR on Tissue Samples To Detect Mycobacterium tuberculosis Complex: an Alternative to the Bacteriological Culture

Bovine tuberculosis (bTB) is an ongoing issue in several countries within the European Union. Microbiological culture is the official confirmation technique for the presence of Mycobacterium tuberculosis complex (MTBC) members in bovine tissues, but several methodological issues, such as moderate sensitivity and long incubation times, require the development of more sensitive and rapid techniques. This study evaluates the analytical and diagnostic performance, comparative to culture, of a real-time PCR targeting the MTBC-specific IS6110 transposon using a panel of bovine tissue samples sourced from the Spanish bTB eradication campaign. Robustness and repeatability were evaluated in an interlaboratory trial between European Union National Reference Laboratories. The limit of detection with 95% confidence was established at 65 fg/reaction of purified genomic equivalents. Diagnostic sensitivity (Se) and specificity (Sp) were, respectively, 96.45% and 93.66%, and the overall agreement (κ) was 0.88. Cross-reactivity was detected against two mycobacterial isolates identified as Mycobacterium marinum and "Mycobacterium avium subsp. hominissuis," and whole-genome sequencing (WGS) analysis of the latter isolate revealed an IS6110-like sequence with 83% identity. An identical IS-like element was found in other Mycobacterium avium complex species in the NCBI nucleotide and WGS databases. Despite this finding, this methodology is considered a valuable alternative to culture, and the strategy of use should be defined depending on the control or eradication programs.

New Mycobacteroides abscessus subsp. massiliense strains with recombinant hsp65 gene laterally transferred from Mycobacteroides abscessus subsp. abscessus: Potential for misidentification of M. abscessus strains with the hsp65-based method

It has been reported that lateral gene transfer (LGT) events among Mycobacteroides abscessus strains are prevalent. The hsp65 gene, a chronometer gene for bacterial phylogenetic analysis, is resistant to LGT events, particularly among mycobacterial strains, rendering the hsp65-targeting method the most widely used method for mycobacterial detection. To determine the prevalence of M. abscessus strains that are subject to hsp65 LGT, we applied rpoB typing to 100 clinically isolated Korean strains of M. abscessus that had been identified by hsp65 sequence analysis. The analysis indicated the presence of 2 rough strains, showing a discrepancy between the 2 typing methods. MLST analysis based on the partial sequencing of seven housekeeping genes, erm(41) PCR and further hsp65 PCR-restriction enzyme and polymorphism analysis (PRA) were conducted to identify the two strains. The MLST results showed that the two strains belong to M. abscessus subsp. massiliense and not to M. abscessus subsp. abscessus, as indicated by the rpoB-based analysis, suggesting that their hsp65 genes are subject to LGT from M. abscessus subsp. abscessus. Further analysis of these strains using the hsp65 PRA method indicated that these strains possess a PRA pattern identical to that of M. abscessus subsp. abscessus and distinct from that of M. abscessus subsp. massiliense. In conclusion, we identified two M. abscessus subsp. massiliense rough strains from Korean patients with hsp65 genes that might be laterally transferred from M. abscessus subsp. abscessus. To the best of our knowledge, this is the first demonstration of possible LGT events associated with the hsp65 gene in mycobacteria. Our results also suggest that there is the potential for misidentification when the hsp65-based protocol is used for mycobacterial identification.

RecA requires two molecules of Mg2+ ions for its optimal strand exchange activity in vitro

Mg²⁺ ion stimulates the DNA strand exchange reaction catalyzed by RecA, a key step in homologous recombination. To elucidate the molecular mechanisms underlying the role of Mg²⁺ and the strand exchange reaction itself, we investigated the interaction of RecA with Mg²⁺ and sought to determine which step of the reaction is affected. Thermal stability, intrinsic fluorescence, and native mass spectrometric analyses of RecA revealed that RecA binds at least two Mg²⁺ ions with K_D ≈ 2 mM and 5 mM. Deletion of the C-terminal acidic tail of RecA made its thermal stability and fluorescence characteristics insensitive to Mg²⁺ and similar to those of full-length RecA in the presence of saturating Mg²⁺. These observations, together with the results of a molecular dynamics simulation, support the idea that the acidic tail hampers the strand exchange reaction by interacting with other parts of RecA, and that binding of Mg²⁺ to the tail prevents these interactions and releases RecA from inhibition. We observed that binding of the first Mg²⁺ stimulated joint molecule formation, whereas binding of the second stimulated progression of the reaction. Thus, RecA is actively involved in the strand exchange step as well as bringing the two DNAs close to each other.

Life at Home and on the Roam: Genomic Adaptions Reflect the Dual Lifestyle of an Intracellular, Facultative Symbiont

"Candidatus Synechococcus feldmannii" is a facultative intracellular symbiont of the Atlanto-Mediterranean sponge Petrosia ficiformis. Genomic information of sponge-associated cyanobacteria derives thus far from the obligate and extracellular symbiont "Candidatus Synechococcus spongiarum." Here we utilized a differential methylation-based approach for bacterial DNA enrichment combined with metagenomics to obtain the first draft genomes of "Ca. Synechococcus feldmannii." By comparative genomics, we revealed that some genomic features (e.g., iron transport mediated by siderophores, eukaryotic-like proteins, and defense mechanisms, like CRISPR-Cas [clustered regularly interspaced short palindromic repeats-associated proteins]) are unique to both symbiont types and absent or rare in the genomes of taxonomically related free-living cyanobacteria. These genomic features likely enable life under the conditions found inside the sponge host. Interestingly, there are many genomic features that are shared by "Ca. Synechococcus feldmannii" and free-living cyanobacteria, while they are absent in the obligate symbiont "Ca. Synechococcus spongiarum." These include genes related to cell surface structures, genetic regulation, and responses to environmental stress, as well as the composition of photosynthetic genes and DNA metabolism. We speculate that the presence of these genes confers on "Ca. Synechococcus feldmannii" its facultative nature (i.e., the ability to respond to a less stable environment when free-living). Our comparative analysis revealed that distinct genomic features depend on the nature of the symbiotic interaction: facultative and intracellular versus obligate and extracellular. IMPORTANCE Given the evolutionary position of sponges as one of the earliest phyla to depart from the metazoan stem lineage, studies on their distinct and exceptionally diverse microbial communities should yield a better understanding of the origin of animal-bacterium interactions. While genomes of several extracellular sponge symbionts have been published, the intracellular symbionts have, so far, been elusive. Here we compare the genomes of two unicellular cyanobacterial sponge symbionts that share an ancestor but followed different evolutionary paths-one became intracellular and the other extracellular. Counterintuitively, the intracellular cyanobacteria are facultative, while the extracellular ones are obligate. By sequencing the genomes of the intracellular cyanobacteria and comparing them to the genomes of the extracellular symbionts and related free-living cyanobacteria, we show how three different cyanobacterial lifestyles are reflected by adaptive genomic features.

A statistical method to identify recombination in bacterial genomes based on SNP incompatibility

BACKGROUND

Phylogeny estimation for bacteria is likely to reflect their true evolutionary histories only if they are highly clonal. However, recombination events could occur during evolution for some species. The reconstruction of phylogenetic trees from an alignment without considering recombination could be misleading, since the relationships among strains in some parts of the genome might be different than in others. Using a single, global tree can create the appearance of homoplasy in recombined regions. Hence, the identification of recombination breakpoints is essential to better understand the evolutionary relationships of isolates among a bacterial population.

RESULTS

Previously, we have developed a method (called ACR) to detect potential breakpoints in an alignment by evaluating compatibility of polymorphic sites in a sliding window. To assess the statistical significance of candidate breakpoints, we propose an extension of the algorithm (ptACR) that applies a permutation test to generate a null distribution for comparing the average local compatibility. The performance of ptACR is evaluated on both simulated and empirical datasets. ptACR is shown to have similar sensitivity (true positive rate) but a lower false positive rate and higher F1 score compared to basic ACR. When used to analyze a collection of clinical isolates of Staphylococcus aureus, ptACR finds clear evidence of recombination events in this bacterial pathogen, and is able to identify statistically significant boundaries of chromosomal regions with distinct phylogenies.

CONCLUSIONS

ptACR is an accurate and efficient method for identifying genomic regions affected by recombination in bacterial genomes.

Population Structure and Local Adaptation of MAC Lung Disease Agent Mycobacterium avium subsp. hominissuis

Mycobacterium avium subsp. hominissuis (MAH) is one of the most common nontuberculous mycobacterial species responsible for chronic lung disease in humans. Despite increasing worldwide incidence, little is known about the genetic mechanisms behind the population evolution of MAH. To elucidate the local adaptation mechanisms of MAH, we assessed genetic population structure, the mutual homologous recombination, and gene content for 36 global MAH isolates, including 12 Japanese isolates sequenced in the present study. We identified five major MAH lineages and found that extensive mutual homologous recombination occurs among them. Two lineages (MahEastAsia1 and MahEastAsia2) were predominant in the Japanese isolates. We identified alleles unique to these two East Asian lineages in the loci responsible for trehalose biosynthesis (treS and mak) and in one mammalian cell entry operon, which presumably originated from as yet undiscovered mycobacterial lineages. Several genes and alleles unique to East Asian strains were located in the fragments introduced via recombination between East Asian lineages, suggesting implication of recombination in local adaptation. These patterns of MAH genomes are consistent with the signature of distribution conjugative transfer, a mode of sexual reproduction reported for other mycobacterial species.

Blending genomes: distributive conjugal transfer in mycobacteria, a sexier form of HGT

This review discusses a novel form of horizontal gene transfer (HGT) found in mycobacteria called Distributive Conjugal Transfer (DCT). While satisfying the criteria for conjugation, DCT occurs by a mechanism so distinct from oriT-mediated conjugation that it could be considered a fourth category of HGT. DCT involves the transfer of chromosomal DNA between mycobacteria and, most significantly, generates transconjugants with mosaic genomes of the parental strains. Multiple segments of donor chromosomal DNA can be co-transferred regardless of their location or the genetic selection and, as a result, the transconjugant genome contains many donor-derived segments; hence the name DCT. This distinguishing feature of DCT separates it from the other known mechanisms of HGT, which generally result in the introduction of a single, defined segment of DNA into the recipient chromosome (Fig. ). Moreover, these mosaic progeny are generated from a single conjugal event, which provides enormous capacity for rapid adaptation and evolution, again distinguishing it from the three classical modes of HGT. Unsurprisingly, the unusual mosaic products of DCT are generated by a conjugal mechanism that is also unusual. Here, we will describe the unique features of DCT and contrast those to other mechanisms of HGT, both from a mechanistic and an evolutionary perspective. Our focus will be on transfer of chromosomal DNA, as opposed to plasmid mobilization, because DCT mediates transfer of chromosomal DNA and is a chromosomally encoded process.

Comparative genomics of the wheat fungal pathogen Pyrenophora tritici-repentis reveals chromosomal variations and genome plasticity

Background

Pyrenophora tritici-repentis (Ptr) is a necrotrophic fungal pathogen that causes the major wheat disease, tan spot. We set out to provide essential genomics-based resources in order to better understand the pathogenicity mechanisms of this important pathogen.

Results

Here, we present eight new Ptr isolate genomes, assembled and annotated; representing races 1, 2 and 5, and a new race. We report a high quality Ptr reference genome, sequenced by PacBio technology with Illumina paired-end data support and optical mapping. An estimated 98% of the genome coverage was mapped to 10 chromosomal groups, using a two-enzyme hybrid approach. The final reference genome was 40.9 Mb and contained a total of 13,797 annotated genes, supported by transcriptomic and proteogenomics data sets.

Conclusions

Whole genome comparative analysis revealed major chromosomal segmental rearrangements and fusions, highlighting intraspecific genome plasticity in this species. Furthermore, the Ptr race classification was not supported at the whole genome level, as phylogenetic analysis did not cluster the ToxA producing isolates. This expansion of available Ptr genomics resources will directly facilitate research aimed at controlling tan spot disease.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4680-3) contains supplementary material, which is available to authorized users.

Phylogenetic analysis of Mycobacterium massiliense strains having recombinant rpoB gene laterally transferred from Mycobacterium abscessus

Recent multi locus sequence typing (MLST) and genome based studies indicate that lateral gene transfer (LGT) events in the rpoB gene are prevalent between Mycobacterium abscessus complex strains. To check the prevalence of the M. massiliense strains subject to rpoB LGT (Rec-mas), we applied rpoB typing (711 bp) to 106 Korean strains of M. massiliense infection that had already been identified by hsp65 sequence analysis (603 bp). The analysis indicated 6 smooth strains in M. massiliense Type I (10.0%, 6/60) genotypes but no strains in M. massiliense Type II genotypes (0%, 0/46), showing a discrepancy between the 2 typing methods. Further MLST analysis based on the partial sequencing of seven housekeeping genes, argH, cya, glpK, gnd, murC, pta and purH, as well as erm(41) PCR proved that these 6 Rec-mas strains consisted of two distinct genotypes belonging to M. massiliense and not M. abscessus. The complete rpoB sequencing analysis showed that these 6 Rec-mas strains have an identical hybrid rpoB gene, of which a 478 bp partial rpoB fragment may be laterally transferred from M. abscessus. Notably, five of the 6 Rec-mas strains showed complete identical sequences in a total of nine genes, including the seven MLST genes, hsp65, and rpoB, suggesting their clonal propagation in South Korea. In conclusion, we identified 6 M. massiliense smooth strains of 2 phylogenetically distinct genotypes with a specific hybrid rpoB gene laterally transferred from M. abscessus from Korean patients. Their clinical relevance and bacteriological traits remain to be elucidated.

Assessment of Metabolic Changes in Mycobacterium smegmatis Wild-Type and alr Mutant Strains: Evidence of a New Pathway of d-Alanine Biosynthesis

In mycobacteria, d-alanine is an essential precursor for peptidoglycan biosynthesis. The only confirmed enzymatic pathway to form d-alanine is through the racemization of l-alanine by alanine racemase (Alr, EC 5.1.1.1). Nevertheless, the essentiality of Alr in Mycobacterium tuberculosis and Mycobacterium smegmatis for cell survivability in the absence of d-alanine has been a point of controversy with contradictory results reported in the literature. To address this issue, we examined the effects of alr inactivation on the cellular metabolism of M. smegmatis. The M. smegmatis alr insertion mutant TAM23 exhibited essentially identical growth to wild-type mc²155 in the absence of d-alanine. NMR metabolomics revealed drastically distinct phenotypes between mc²155 and TAM23. A metabolic switch was observed for TAM23 as a function of supplemented d-alanine. In the absence of d-alanine, the metabolic response directed carbon through an unidentified transaminase to provide the essential d-alanine required for survival. The process is reversed when d-alanine is available, in which the d-alanine is directed to peptidoglycan biosynthesis. Our results provide further support for the hypothesis that Alr is not an essential function of M. smegmatis and that specific Alr inhibitors will have no bactericidal action.

Comparative functional pan-genome analyses to build connections between genomic dynamics and phenotypic evolution in polycyclic aromatic hydrocarbon metabolism in the genus Mycobacterium

Background

The bacterial genus Mycobacterium is of great interest in the medical and biotechnological fields. Despite a flood of genome sequencing and functional genomics data, significant gaps in knowledge between genome and phenome seriously hinder efforts toward the treatment of mycobacterial diseases and practical biotechnological applications. In this study, we propose the use of systematic, comparative functional pan-genomic analysis to build connections between genomic dynamics and phenotypic evolution in polycyclic aromatic hydrocarbon (PAH) metabolism in the genus Mycobacterium.

Results

Phylogenetic, phenotypic, and genomic information for 27 completely genome-sequenced mycobacteria was systematically integrated to reconstruct a mycobacterial phenotype network (MPN) with a pan-genomic concept at a network level. In the MPN, mycobacterial phenotypes show typical scale-free relationships. PAH degradation is an isolated phenotype with the lowest connection degree, consistent with phylogenetic and environmental isolation of PAH degraders. A series of functional pan-genomic analyses provide conserved and unique types of genomic evidence for strong epistatic and pleiotropic impacts on evolutionary trajectories of the PAH-degrading phenotype. Under strong natural selection, the detailed gene gain/loss patterns from horizontal gene transfer (HGT)/deletion events hypothesize a plausible evolutionary path, an epistasis-based birth and pleiotropy-dependent death, for PAH metabolism in the genus Mycobacterium. This study generated a practical mycobacterial compendium of phenotypic and genomic changes, focusing on the PAH-degrading phenotype, with a pan-genomic perspective of the evolutionary events and the environmental challenges.

Conclusions

Our findings suggest that when selection acts on PAH metabolism, only a small fraction of possible trajectories is likely to be observed, owing mainly to a combination of the ambiguous phenotypic effects of PAHs and the corresponding pleiotropy- and epistasis-dependent evolutionary adaptation. Evolutionary constraints on the selection of trajectories, like those seen in PAH-degrading phenotypes, are likely to apply to the evolution of other phenotypes in the genus Mycobacterium.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0302-8) contains supplementary material, which is available to authorized users.

Comparative genomic analysis of Mycobacterium iranicum UM_TJL against representative mycobacterial species suggests its environmental origin

Mycobacterium iranicum is a newly reported mycobacterial species. We present the first comparative study of M. iranicum UM_TJL and other mycobacteria. We found M. iranicum to have a close genetic association with environmental mycobacteria infrequently associated with human infections. Nonetheless, UM_TJL is also equipped with many virulence genes (some of which appear to be the consequence of transduction-related gene transfer) that have been identified in established human pathogens. Taken all together, our data suggest that M. iranicum is an environmental bacterium adapted for pathogenicity in the human host. This comparative study provides important clues and forms the basis for future functional studies on this mycobacterium.

Genomic signatures of distributive conjugal transfer among mycobacteria

Distributive conjugal transfer (DCT) is a newly described mechanism of lateral gene transfer (LGT) that results in a mosaic transconjugant structure, similar to the products of meiosis. We have tested popular LGT detection methods on whole-genome sequence data from experimental DCT transconjugants and used the best performing methods to compare genomic signatures of DCT with those of LGT through natural transformation, conjugative plasmids, and mobile genetic elements (MGE). We found that DCT results in transfer of larger chromosomal segments, that these segments are distributed more broadly around the chromosome, and that a greater proportion of the chromosome is affected by DCT than by other mechanisms of LGT. We used the best performing methods to characterize LGT in Mycobacterium canettii, the mycobacterial species most closely related to Mycobacterium tuberculosis. Patterns of LGT among M. canettii were highly distinctive. Gene flow appeared unidirectional, from lineages with minimal evidence of LGT to isolates with a substantial proportion (6–13%) of sites identified as recombinant. Among M. canettii isolates with evidence of LGT, recombinant fragments were larger and more evenly distributed relative to bacteria that undergo LGT through natural transformation, conjugative plasmids, and MGE. Spatial bias in M. canettii was also unusual in that patterns of recombinant fragment sharing mirrored overall phylogenetic structure. Based on the proportion of recombinant sites, the size of recombinant fragments, their spatial distribution and lack of association with MGE, as well as unidirectionality of DNA transfer, we conclude that DCT is the predominant mechanism of LGT among M. canettii.

Distributive Conjugal Transfer: New Insights into Horizontal Gene Transfer and Genetic Exchange in Mycobacteria

The last decade has seen an explosion in the application of genomic tools across all biological disciplines. This is also true for mycobacteria, where whole genome sequences are now available for pathogens and non-pathogens alike. Genomes within the Mycobacterium tuberculosis Complex (MTBC) bear the hallmarks of horizontal gene transfer (HGT). Conjugation is the form of HGT with the highest potential capacity and evolutionary influence. Donor and recipient strains of Mycobacterium smegmatis actively conjugate upon co-culturing in biofilms and on solid media. Whole genome sequencing of the transconjugant progeny demonstrated the incredible scale and range of genomic variation that conjugation generates. Transconjugant genomes are complex mosaics of the parental strains. Some transconjugant genomes are up to one-quarter donor-derived, distributed over 30 segments. Transferred segments range from ~50 bp to ~225,000 bp in length, and are exchanged with their recipient orthologs all around the genome. This unpredictable genome-wide infusion of DNA sequences is called Distributive Conjugal Transfer (DCT), to distinguish it from traditional oriT-based conjugation. The mosaicism generated in a single transfer event resembles that seen from meiotic recombination in sexually reproducing organisms, and contrasts with traditional models of HGT. This similarity allowed the application of a GWAS-like approach to map the donor genes that confer a donor mating identity phenotype. The mating identity genes map to the esx1 locus, expanding the central role of ESX-1 function in conjugation. The potential for DCT to instantaneously blend genomes will affect how we view mycobacterial evolution, and provide new tools for the facile manipulation of mycobacterial genomes.

Evidence for diversifying selection in a set of Mycobacterium tuberculosis genes in response to antibiotic- and nonantibiotic-related pressure

Tuberculosis (TB) is a global health problem estimated to kill 1.4 million people per year. Recent advances in the genomics of the causative agents of TB, bacteria known as the Mycobacterium tuberculosis complex (MTBC), have allowed a better comprehension of its population structure and provided the foundation for molecular evolution analyses. These studies are crucial for a better understanding of TB, including the variation of vaccine efficacy and disease outcome, together with the emergence of drug resistance. Starting from the analysis of 73 publicly available genomes from all the main MTBC lineages, we have screened for evidences of positive selection, a set of 576 genes previously associated with drug resistance or encoding membrane proteins. As expected, because antibiotics constitute strong selective pressure, some of the codons identified correspond to the position of confirmed drug-resistance-associated substitutions in the genes embB, rpoB, and katG. Furthermore, we identified diversifying selection in specific codons of the genes Rv0176 and Rv1872c coding for MCE1-associated transmembrane protein and a putative l-lactate dehydrogenase, respectively. Amino acid sequence analyses showed that in Rv0176, sites undergoing diversifying selection were in a predicted antigen region that varies between "modern" lineages and "ancient" MTBC/BCG strains. In Rv1872c, some of the sites under selection are predicted to impact protein function and thus might result from metabolic adaptation. These results illustrate that diversifying selection in MTBC is happening as a consequence of both antibiotic treatment and other evolutionary pressures.

The apoptogenic toxin AIP56 is a metalloprotease A-B toxin that cleaves NF-κb P65

AIP56 (apoptosis-inducing protein of 56 kDa) is a major virulence factor of Photobacterium damselae piscicida (Phdp), a Gram-negative pathogen that causes septicemic infections, which are among the most threatening diseases in mariculture. The toxin triggers apoptosis of host macrophages and neutrophils through a process that, in vivo, culminates with secondary necrosis of the apoptotic cells contributing to the necrotic lesions observed in the diseased animals. Here, we show that AIP56 is a NF-κB p65-cleaving zinc-metalloprotease whose catalytic activity is required for the apoptogenic effect. Most of the bacterial effectors known to target NF-κB are type III secreted effectors. In contrast, we demonstrate that AIP56 is an A-B toxin capable of acting at distance, without requiring contact of the bacteria with the target cell. We also show that the N-terminal domain cleaves NF-κB at the Cys³⁹-Glu⁴⁰ peptide bond and that the C-terminal domain is involved in binding and internalization into the cytosol.

The apoptosis inducing protein of 56 kDa (AIP56) is a key virulence factor secreted by Photobacterium damselae piscicida (Phdp), a Gram-negative bacterium that causes septicaemic infections in economically important marine fish species. It is known that AIP56 induces massive destruction of the phagocytic cells of the infected host, allowing the extracellular multiplication of the bacteria and contributing to the genesis of the pathology. Here we show that AIP56 acts by cleaving NF-κB p65. The NF-κB family of transcription factors is evolutionarily conserved and plays a central role in the host responses to microbial pathogen invasion, regulating the expression of inflammatory and anti-apoptotic genes. Pathogenic bacteria have evolved complex strategies to interfere with NF-κB signalling, usually by injecting protein effectors directly into the cell's cytosol through bacterial secretion machineries that require contact with host cells. In contrast, AIP56 acts at distance and has an intrinsic ability to reach the cytosol due to the presence of a C-terminal domain that functions as “delivery module.”

Difference in virulence of Mycobacterium avium isolates sharing indistinguishable DNA fingerprint determined in murine model of lung infection

BACKGROUND

Opportunistic Mycobacterium avium typically causes disease in immunocompromised patients and in some groups of apparently healthy individuals. The high virulence of some bacterial lineages increases the disease risk. High-resolution molecular genotyping studies of M. avium clinical isolates demonstrated that some genotype patterns were more prevalent than others, suggesting that close genetic relatedness of these successful isolates sharing a similar genotype could determine similar biological properties associated with high virulence.

METHODS AND FINDINGS

In this study, we aimed to compare the virulence and pathogenic properties of two epidemiologically unrelated M. avium isolates sharing an indistinguishable DNA fingerprint in a well-characterized model of pulmonary infection in mice, resistant or susceptible to mycobacteria. The mice, C57BL/6 wild- type or IFN-gamma gene disrupted (GKO), respectively, were intratracheally infected with two isolates, H27 (human blood isolate) and P104 (pig lymph node isolate), and the lungs were examined for bacterial loads, histopathology and cytokine gene expression. The obtained data demonstrated significant differences in the virulence properties of these strains. Although the H27 strain grew significantly faster than P104 in the early stage of infection, this bacterium induced protective immunity that started to reduce bacterial numbers in the wild-type mice, whereas the P104 strain established a chronic infection. In the GKO mice, both strains were capable of causing a chronic infection, associated with higher bacterial burdens and severe lung pathology, in a similar manner.

CONCLUSIONS/SIGNIFICANCE

The results demonstrated that the studied isolates differed in the pathogenic properties although were indistinguishable by actually widely used genotyping techniques demonstrating that the genotype similarity does not predict similarity in virulence of M. avium isolates.

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.

Evidence of extensive homologous recombination in the core genome of rickettsia

The important role of homologous recombination has been extensively demonstrated to be fundamental for genetic variation in bacterial genomes. In contrast to extracellular or facultative intracellular bacteria, obligate intracellular bacteria are considered to be less prone to recombination, especially for their core genomes. In Rickettsia, only antigen-related genes were identified to have experienced homologous recombination. In this study, we employed evolutionary genomic approaches to investigate the impact of recombination on the core genome of Rickettsia. Phylogenetic network and phylogenetic compatibility matrix analyses are clearly consistent with the hypothesis that recombination has occurred frequently during Rickettsia evolution. 28% of Rickettsia core genes (194 out of 690) are found to present the evidence of recombination under four independent statistical methods. Further functional classification shows that these recombination events occur across all functional categories, with a significant overrepresentation in the cell wall/membrane/envelope biogenesis, which may provide a molecular basis for the parasite adaptation to host immunity. This evolutionary genomic analysis provides insight into the substantial role of recombination in the evolution of the intracellular pathogenic bacteria Rickettsia.

Frequent homologous recombination events in Mycobacterium tuberculosis PE/PPE multigene families: potential role in antigenic variability

The PE and PPE (PE/PPE) multigene families of Mycobacterium tuberculosis are particularly GC-rich and share extensive homologous repetitive sequences. We hypothesized that they may undergo homologous recombination events, a mechanism rarely described in the natural evolution of mycobacteria. To test our hypothesis, we developed a specific oligonucleotide-based microarray targeting nearly all of the PE/PPE genes, aimed at detecting signals for homologous recombination. Such a microarray has never before been reported due to the multiplicity and highly repetitive and homologous nature of these sequences. Application of the microarray to a collection of M. tuberculosis clinical isolates (n = 33) representing prevalent spoligotype strain families in Tunisia allowed successful detection of six deleted genomic regions involving a total of two PE and seven PPE genes. Some of these deleted genes are known to be immunodominant or involved in virulence. The four precisely determined deletions were flanked by 400- to 500-bp stretches of nearly identical sequences lying mainly at the conserved N-terminal region of the PE/PPE genes. These highly homologous sequences thus serve as substrates to mediate both intergenic and intragenic homologous recombination events, indicating an important function in generating strain variation. Importantly, all recombination events yielded a new in-frame fusion chimeric gene. Hence, homologous recombination within and between PE/PPE genes likely increased their antigenic variability, which may have profound implications in pathogenicity and/or host adaptation. The finding of high prevalence (approximately 45% and approximately 58%) for at least two of the genomic deletions suggests that they likely confer advantageous biological attributes.

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.

The mycobacterial glycopeptidolipids: structure, function, and their role in pathogenesis

Glycopeptidolipids (GPLs) are a class of glycolipids produced by several nontuberculosis-causing members of the Mycobacterium genus including pathogenic and nonpathogenic species. GPLs are expressed in different forms with production of highly antigenic, typeable serovar-specific GPLs in members of the Mycobacterium avium complex (MAC). M. avium and M. intracellulare, which comprise this complex, are slow-growing mycobacteria noted for producing disseminated infections in AIDS patients and pulmonary infections in non-AIDS patients. Previous studies have defined the gene cluster responsible for GPL biosynthesis and more recent work has characterized the function of the individual genes. Current research has also focused on the GPL's role in colony morphology, sliding motility, biofilm formation, immune modulation and virulence. These topics, along with new information on the enzymes involved in GPL biosynthesis, are the subject of this review.

Mycobacterium avium subsp. paratuberculosis and M. avium subsp. avium are independently evolved pathogenic clones of a much broader group of M. avium organisms

Mycobacterium avium comprises organisms that share the same species designation despite considerable genomic and phenotypic variability. To determine the degree and nature of variability between subspecies and strains of M. avium, we used multilocus sequencing analysis, studying 56 genetically diverse strains of M. avium that included all described subspecies. In total, 8,064 bp of sequence from 10 gene loci were studied, with 205 (2.5%) representing variable positions. The majority (149/205) of these variations were found among M. avium subsp. hominissuis organisms. Recombination was also evident in this subspecies. In contrast, there was comparatively little variability and no evidence of recombination within the pathogenic subspecies, M. avium subsp. paratuberculosis, M. avium subsp. avium, and M. avium subsp. silvaticum. Phylogenetic analysis showed that M. avium subsp. avium and M. avium subsp. silvaticum strains clustered together on one branch, while a distinct branch defined M. avium subsp. paratuberculosis organisms. Despite the independent origin of these pathogenic subspecies, an analysis of their rates of nonsynonymous (dN) to synonymous (dS) substitutions showed increased dN/dS ratios for both: 0.67 for M. avium subsp. paratuberculosis and 0.50 for M. avium subsp. avium/M. avium subsp. silvaticum, while the value was 0.08 for M. avium subsp. hominissuis organisms. In conclusion, M. avium subsp. hominissuis represents a diverse group of organisms from which two pathogenic clones (M. avium subsp. paratuberculosis and M. avium subsp. avium/M. avium subsp. silvaticum) have evolved independently.

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.

Terminal reassortment drives the quantum evolution of type III effectors in bacterial pathogens

Many bacterial pathogens employ a type III secretion system to deliver type III secreted effectors (T3SEs) into host cells, where they interact directly with host substrates to modulate defense pathways and promote disease. This interaction creates intense selective pressures on these secreted effectors, necessitating rapid evolution to overcome host surveillance systems and defenses. Using computational and evolutionary approaches, we have identified numerous mosaic and truncated T3SEs among animal and plant pathogens. We propose that these secreted virulence genes have evolved through a shuffling process we have called “terminal reassortment.” In terminal reassortment, existing T3SE termini are mobilized within the genome, creating random genetic fusions that result in chimeric genes. Up to 32% of T3SE families in species with relatively large and well-characterized T3SE repertoires show evidence of terminal reassortment, as compared to only 7% of non-T3SE families. Terminal reassortment may permit the near instantaneous evolution of new T3SEs and appears responsible for major modifications to effector activity and function. Because this process plays a more significant role in the evolution of T3SEs than non-effectors, it provides insight into the evolutionary origins of T3SEs and may also help explain the rapid emergence of new infectious agents.

Many pathogenic bacteria rely on specialized virulence proteins to cause disease. These proteins, known as type III secreted effectors (T3SEs), are directly injected into the host's cells and facilitate the disease process by interacting with host proteins and interfering with the defense response. Although most T3SEs lack any sequence similarity, several T3SEs share a common terminus, suggesting that part of these proteins was derived from the same sequence. The authors propose an evolutionary mechanism, called “terminal reassortment,” in which the termini of T3SEs reassort with other genetic information to create new chimeric proteins. This study shows that this process has given rise to T3SEs with new virulence functions and that it may influence bacterial host specificity. Chimeric T3SEs are present in eight different genera and in some cases are present in as many as 32% of known T3SE families. This is significantly more than what is found in other protein families, suggesting that terminal reassortment plays a disproportionately important role in the evolution of T3SE. Terminal reassortment may lead to the very rapid evolution of new T3SEs, thereby contributing to the emergence of new infectious diseases.

Comparative analysis and insights into the evolution of gene clusters for glycopeptide antibiotic biosynthesis

The bal, cep, dbv, sta and tcp gene clusters specify the biosynthesis of the glycopeptide antibiotics balhimycin, chloroeremomycin, A40926, A47934 and teicoplanin, respectively. These structurally related compounds share a similar mechanism of action in their inhibition of bacterial cell wall formation. Comparative sequence analysis was performed on the five gene clusters. Extensive conserved synteny was observed between the bal and cep clusters, which direct the synthesis of very similar compounds but originate from two different species of the genus Amycolatopsis. All other cluster pairs show a limited degree of conserved synteny, involving biosynthetically functional gene cassettes: these include those involved in the synthesis of the carbon backbone of two non-proteinogenic amino acids; in the linkage of amino acids 1--3 and 4--7 in the heptapeptide; and in the formation of the aromatic cross-links. Furthermore, these segments of conserved synteny are often preceded by conserved intergenic regions. Phylogenetic analysis of protein families shows several instances in which relatedness in the chemical structure of the glycopeptides is not reflected in the extent of the relationship of the corresponding polypeptides. Coherent branchings are observed for all polypeptides encoded by the syntenous gene cassettes. These results suggest that the acquisition of distinct, functional genetic elements has played a significant role in the evolution of glycopeptide gene clusters, giving them a mosaic structure. In addition, the synthesis of the structurally similar compounds A40926 and teicoplanin appears as the result of convergent evolution.

Phylogeny of Mycobacterium avium strains inferred from glycopeptidolipid biosynthesis pathway genes

The Mycobacterium avium complex (MAC) encompasses two species, M. avium and Mycobacterium intracellulare, which are opportunistic pathogens of humans and animals. The standard method of MAC strain differentiation is serotyping based on a variation in the antigenic glycopeptidolipid (GPL) composition. To elucidate the relationships among M. avium serotypes a phylogenetic analysis of 13 reference and clinical M. avium strains from 8 serotypes was performed using as markers two genomic regions (890 bp of the gtfB gene and 2150 bp spanning the rtfA-mtfC genes) which are associated with the strains' serological properties. Strains belonging to three other known M. avium serotypes were not included in the phylogeny inference due to apparent lack of the marker sequences in their genomes, as revealed by PCR and Southern blot analysis. These studies suggest that serotypes prevalent in AIDS patients have multiple origins. In trees inferred from both markers, serotype 1 strains, known to have the simplest and shortest GPLs among all other serotypes, were polyphyletic. Likewise, comparisons of the inferred phylogenies with the molecular typing results imply that the existing tools used in epidemiological studies may be poor estimators of M. avium strain relatedness. Additionally, trees inferred from each marker had significantly incongruent topologies due to a well supported alternative placement of strain 2151, suggesting a complex evolutionary history of this genomic region.