Characterization of a cross-reactive, immunodominant and HLA-promiscuous epitope of Mycobacterium tuberculosis-specific major antigenic protein PPE68.

PPE68 (Rv3873), a major antignic protein encoded by Mycobacteriun tuberculosis-specific genomic region of difference (RD)1, is a strong stimulator of peripheral blood mononuclear cells (PBMCs) obtained from tuberculosis patients and Mycobacterium bovis bacillus Calmette Guerin (BCG)-vaccianted healthy subjects in T helper (Th)1 cell assays, i.e. antigen-induced proliferation and interferon-gamma (IFN-γ) secretion. To confirm the antigen-specific recognition of PPE68 by T cells in IFN-γ assays, antigen-induced human T-cell lines were established from PBMCs of M. Bovis BCG-vaccinated and HLA-heterogeneous healthy subjects and tested with peptide pools of RD1 proteins. The results showed that PPE68 was recognized by antigen-specific T-cell lines from HLA-heteregeneous subjects. To further identify the immunodominant and HLA-promiscuous Th1-1 cell epitopes present in PPE68, 24 synthetic peptides covering the sequence of PPE68 were indivdually analyzed for HLA-DR binding prediction analysis and tested with PBMCs from M. bovis BCG-vaccinated and HLA-heterogeuous healthy subjects in IFN-γ assays. The results identified the peptide P9, i.e. aa 121-VLTATNFFGINTIPIALTEMDYFIR-145, as an immunodominant and HLA-DR promiscuous peptide of PPE68. Furthermore, by using deletion peptides, the immunodominant and HLA-DR promiscuous core sequence was mapped to aa 127-FFGINTIPIA-136. Interestingly, the core sequence is present in several PPE proteins of M. tuberculosis, and conserved in all sequenced strains/species of M. tuberculosis and M. tuberculosis complex, and several other pathogenic mycobacterial species, including M. leprae and M. avium-intracellulalae complex. These results suggest that the peptide aa 121-145 may be exploited as a peptide-based vaccine candidate against tuberculosis and other mycobacterial diseases.

Comparative genomics for mycobacterial peptidoglycan remodelling enzymes reveals extensive genetic multiplicity.

BACKGROUND: Mycobacteria comprise diverse species including non-pathogenic, environmental organisms, animal disease agents and human pathogens, notably Mycobacterium tuberculosis. Considering that the mycobacterial cell wall constitutes a significant barrier to drug penetration, the aim of this study was to conduct a comparative genomics analysis of the repertoire of enzymes involved in peptidoglycan (PG) remodelling to determine the potential of exploiting this area of bacterial metabolism for the discovery of new drug targets. RESULTS: We conducted an in silico analysis of 19 mycobacterial species/clinical strains for the presence of genes encoding resuscitation promoting factors (Rpfs), penicillin binding proteins, endopeptidases, L,D-transpeptidases and N-acetylmuramoyl-L-alanine amidases. Our analysis reveals extensive genetic multiplicity, allowing for classification of mycobacterial species into three main categories, primarily based on their rpf gene complement. These include the M. tuberculosis Complex (MTBC), other pathogenic mycobacteria and environmental species. The complement of these genes within the MTBC and other mycobacterial pathogens is highly conserved. In contrast, environmental strains display significant genetic expansion in most of these gene families. Mycobacterium leprae retains more than one functional gene from each enzyme family, underscoring the importance of genetic multiplicity for PG remodelling. Notably, the highest degree of conservation is observed for N-acetylmuramoyl-L-alanine amidases suggesting that these enzymes are essential for growth and survival. CONCLUSION: PG remodelling enzymes in a range of mycobacterial species are associated with extensive genetic multiplicity, suggesting functional diversification within these families of enzymes to allow organisms to adapt.

The antibiotic resistance arrow of time: efflux pump induction is a general first step in the evolution of mycobacterial drug resistance.

We hypothesize that low-level efflux pump expression is the first step in the development of high-level drug resistance in mycobacteria. We performed 28-day azithromycin dose-effect and dose-scheduling studies in our hollow-fiber model of disseminated Mycobacterium avium-M. intracellulare complex. Both microbial kill and resistance emergence were most closely linked to the within-macrophage area under the concentration-time curve (AUC)/MIC ratio. Quantitative PCR revealed that subtherapeutic azithromycin exposures over 3 days led to a 56-fold increase in expression of MAV_3306, which encodes a putative ABC transporter, and MAV_1406, which encodes a putative major facilitator superfamily pump, in M. avium. By day 7, a subpopulation of M. avium with low-level resistance was encountered and exhibited the classic inverted U curve versus AUC/MIC ratios. The resistance was abolished by an efflux pump inhibitor. While the maximal microbial kill started to decrease after day 7, a population with high-level azithromycin resistance appeared at day 28. This resistance could not be reversed by efflux pump inhibitors. Orthologs of pumps encoded by MAV_3306 and MAV_1406 were identified in Mycobacterium tuberculosis, Mycobacterium leprae, Mycobacterium marinum, Mycobacterium abscessus, and Mycobacterium ulcerans. All had highly conserved protein secondary structures. We propose that induction of several efflux pumps is the first step in a general pathway to drug resistance that eventually leads to high-level chromosomal-mutation-related resistance in mycobacteria as ordered events in an "antibiotic resistance arrow of time."

Homology modeling and docking analyses of M. leprae Mur ligases reveals the common binding residues for structure based drug designing to eradicate leprosy.

Multi drug resistance capacity for Mycobacterium leprae (MDR-Mle) demands the profound need for developing new anti-leprosy drugs. Since most of the drugs target a single enzyme, mutation in the active site renders the antibiotic ineffective. However, structural and mechanistic information on essential bacterial enzymes in a pathway could lead to the development of antibiotics that targets multiple enzymes. Peptidoglycan is an important component of the cell wall of M. leprae. The biosynthesis of bacterial peptidoglycan represents important targets for the development of new antibacterial drugs. Biosynthesis of peptidoglycan is a multi-step process that involves four key Mur ligase enzymes: MurC (EC:6.3.2.8), MurD (EC:6.3.2.9), MurE (EC:6.3.2.13) and MurF (EC:6.3.2.10). Hence in our work, we modeled the three-dimensional structure of the above Mur ligases using homology modeling method and analyzed its common binding features. The residues playing an important role in the catalytic activity of each of the Mur enzymes were predicted by docking these Mur ligases with their substrates and ATP. The conserved sequence motifs significant for ATP binding were predicted as the probable residues for structure based drug designing. Overall, the study was successful in listing significant and common binding residues of Mur enzymes in peptidoglycan pathway for multi targeted therapy.

An enlarged, adaptable active site in CYP164 family P450 enzymes, the sole P450 in Mycobacterium leprae.

CYP164 family P450 enzymes are found in only a subset of mycobacteria and include CYP164A1, which is the sole P450 found in Mycobacterium leprae, the causative agent of leprosy. This has previously led to interest in this enzyme as a potential drug target. Here we describe the first crystal structure of a CYP164 enzyme, CYP164A2 from Mycobacterium smegmatis. CYP164A2 has a distinctive, enlarged hydrophobic active site that extends above the porphyrin ring toward the access channels. Unusually, we find that CYP164A2 can simultaneously bind two econazole molecules in different regions of the enlarged active site and is accompanied by the rearrangement and ordering of the BC loop. The primary location is through a classic interaction of the azole group with the porphyrin iron. The second econazole molecule is bound to a unique site and is linked to a tetracoordinated metal ion complexed to one of the heme carboxylates and to the side chains of His 105 and His 364. All of these features are preserved in the closely homologous M. leprae CYP164A1. The computational docking of azole compounds to a homology model of CYP164A1 suggests that these compounds will form effective inhibitors and is supported by the correlation of parallel docking with experimental binding studies of CYP164A2. The binding of econazole to CYP164A2 occurs primarily through the high-spin "open" conformation of the enzyme (K(d) [dissociation constant] of 0.1 µM), with binding to the low-spin "closed" form being significantly hindered (K(d) of 338 µM). These studies support previous suggestions that azole derivatives may provide an effective strategy to improve the treatment of leprosy.

Identification, characterization, and azole-binding properties of Mycobacterium smegmatis CYP164A2, a homolog of ML2088, the sole cytochrome P450 gene of Mycobacterium leprae.

The genome sequence of Mycobacterium leprae revealed a single open reading frame, ML2088 (CYP164A1), encoding a putative full-length cytochrome P450 monooxygenase and 12 pseudogenes. We have identified a homolog of ML2088 in Mycobacterium smegmatis and report here the cloning, expression, purification, and azole-binding characteristics of this cytochrome P450 (CYP164A2). CYP164A2 is 1,245 bp long and encodes a protein of 414 amino acids and molecular mass of 45 kDa. CYP164A2 has 60% identity with Mycobacterium leprae CYP161A1 and 66 to 69% identity with eight other mycobacterial CYP164A1 homologs, with three identified highly conserved motifs. Recombinant CYP164A2 has the typical spectral characteristics of a cytochrome P450 monooxygenase, predominantly in the ferric low-spin state. Unusually, the spin state was readily modulated by increasing ionic strength at pH 7.5, with 50% high-spin occupancy achieved with 0.14 M NaCl. CYP164A2 bound clotrimazole, econazole, and miconazole strongly (K(d), 1.2 to 2.5 muM); however, strong binding with itraconazole, ketoconazole, and voriconazole was only observed in the presence of 0.5 M NaCl. Fluconazole did not bind to CYP164A2 at pH 7.5 and no discernible type II binding spectrum was observed.

Genome scale portrait of cAMP-receptor protein (CRP) regulons in mycobacteria points to their role in pathogenesis.

cAMP Receptor Protein (CRP)/Fumarate Nitrate Reductase Regulator (FNR) family proteins are ubiquitous regulators of cell stress in eubacteria. These proteins are commonly associated with maintenance of intracellular oxygen levels, redox-state, oxidative and nitrosative stresses, and extreme temperature conditions by regulating expression of target genes that contain regulatory cognate DNA elements. We describe the use of informatics enabled comparative genomics to identify novel genes under the control of CRP regulator in Mycobacterium tuberculosis (M.tb). An inventory of CRP regulated genes and their operon context in important mycobacterial species such as M. leprae, M. avium subsp. paratuberculosis and M. smegmatis and several common genes within this genus including the important cellular functions, mainly, cell-wall biogenesis, cAMP signaling and metabolism associated with such regulons were identified. Our results provide a possible theoretical framework for better understanding of the stress response in mycobacteria. The conservation of the CRP regulated genes in pathogenic mycobacteria, as opposed to non-pathogenic ones, highlights the importance of CRP-regulated genes in pathogenesis.

The reductase that catalyzes mycolic motif synthesis is required for efficient attachment of mycolic acids to arabinogalactan.

Mycolic acids are essential components of the cell walls of bacteria belonging to the suborder Corynebacterineae, including the important human pathogens Mycobacterium tuberculosis and Mycobacterium leprae. Mycolic acid biosynthesis is complex and the target of several frontline antimycobacterial drugs. The condensation of two fatty acids to form a 2-alkyl-3-keto mycolate precursor and the subsequent reduction of this precursor represent two key and highly conserved steps in this pathway. Although the enzyme catalyzing the condensation step has recently been identified, little is known about the putative reductase. Using an extensive bioinformatic comparison of the genomes of M. tuberculosis and Corynebacterium glutamicum, we identified NCgl2385, the orthologue of Rv2509 in M. tuberculosis, as a potential reductase candidate. Deletion of the gene in C. glutamicum resulted in a slow growing strain that was deficient in arabinogalactan-linked mycolates and synthesized abnormal forms of the mycolate-containing glycolipids trehalose dicorynomycolate and trehalose monocorynomycolate. Analysis of the native and acetylated trehalose glycolipids by MALDI-TOF mass spectrometry indicated that these novel glycolipids contained an unreduced beta-keto ester. This was confirmed by analysis of sodium borodeuteride-reduced mycolic acids by gas chromatography mass spectrometry. Reintroduction of the NCgl2385 gene into the mutant restored the transfer of mature mycolic acids to both the trehalose glycolipids and cell wall arabinogalactan. These data indicate that NCgl2385, which we have designated CmrA, is essential for the production of mature trehalose mycolates and subsequent covalent attachment of mycolic acids onto the cell wall, thus representing a focus for future structural and pathogenicity studies.

Regulatory genes controlling fatty acid catabolism and peroxisomal functions in the filamentous fungus Aspergillus nidulans.

The catabolism of fatty acids is important in the lifestyle of many fungi, including plant and animal pathogens. This has been investigated in Aspergillus nidulans, which can grow on acetate and fatty acids as sources of carbon, resulting in the production of acetyl coenzyme A (CoA). Acetyl-CoA is metabolized via the glyoxalate bypass, located in peroxisomes, enabling gluconeogenesis. Acetate induction of enzymes specific for acetate utilization as well as glyoxalate bypass enzymes is via the Zn2-Cys6 binuclear cluster activator FacB. However, enzymes of the glyoxalate bypass as well as fatty acid beta-oxidation and peroxisomal proteins are also inducible by fatty acids. We have isolated mutants that cannot grow on fatty acids. Two of the corresponding genes, farA and farB, encode two highly conserved families of related Zn2-Cys6 binuclear proteins present in filamentous ascomycetes, including plant pathogens. A single ortholog is found in the yeasts Candida albicans, Debaryomyces hansenii, and Yarrowia lipolytica, but not in the Ashbya, Kluyveromyces, Saccharomyces lineage. Northern blot analysis has shown that deletion of the farA gene eliminates induction of a number of genes by both short- and long-chain fatty acids, while deletion of the farB gene eliminates short-chain induction. An identical core 6-bp in vitro binding site for each protein has been identified in genes encoding glyoxalate bypass, beta-oxidation, and peroxisomal functions. This sequence is overrepresented in the 5' region of genes predicted to be fatty acid induced in other filamentous ascomycetes, C. albicans, D. hansenii, and Y. lipolytica, but not in the corresponding genes in Saccharomyces cerevisiae.

The RNA polymerase III-dependent family of genes in hemiascomycetes: comparative RNomics, decoding strategies, transcription and evolutionary implications.

We present the first comprehensive analysis of RNA polymerase III (Pol III) transcribed genes in ten yeast genomes. This set includes all tRNA genes (tDNA) and genes coding for SNR6 (U6), SNR52, SCR1 and RPR1 RNA in the nine hemiascomycetes Saccharomyces cerevisiae, Saccharomyces castellii, Candida glabrata, Kluyveromyces waltii, Kluyveromyces lactis, Eremothecium gossypii, Debaryomyces hansenii, Candida albicans, Yarrowia lipolytica and the archiascomycete Schizosaccharomyces pombe. We systematically analysed sequence specificities of tRNA genes, polymorphism, variability of introns, gene redundancy and gene clustering. Analysis of decoding strategies showed that yeasts close to S.cerevisiae use bacterial decoding rules to read the Leu CUN and Arg CGN codons, in contrast to all other known Eukaryotes. In D.hansenii and C.albicans, we identified a novel tDNA-Leu (AAG), reading the Leu CUU/CUC/CUA codons with an unusual G at position 32. A systematic 'p-distance tree' using the 60 variable positions of the tRNA molecule revealed that most tDNAs cluster into amino acid-specific sub-trees, suggesting that, within hemiascomycetes, orthologous tDNAs are more closely related than paralogs. We finally determined the bipartite A- and B-box sequences recognized by TFIIIC. These minimal sequences are nearly conserved throughout hemiascomycetes and were satisfactorily retrieved at appropriate locations in other Pol III genes.

A surprisingly large RNase P RNA in Candida glabrata.

We have found an extremely large ribonuclease P (RNase P) RNA (RPR1) in the human pathogen Candida glabrata and verified that this molecule is expressed and present in the active enzyme complex of this hemiascomycete yeast. A structural alignment of the C. glabrata sequence with 36 other hemiascomycete RNase P RNAs (abbreviated as P RNAs) allows us to characterize the types of insertions. In addition, 15 P RNA sequences were newly characterized by searching in the recently sequenced genomes Candida albicans, C. glabrata, Debaryomyces hansenii, Eremothecium gossypii, Kluyveromyces lactis, Kluyveromyces waltii, Naumovia castellii, Saccharomyces kudriavzevii, Saccharomyces mikatae, and Yarrowia lipolytica; and by PCR amplification for other Candida species (Candida guilliermondii, Candida krusei, Candida parapsilosis, Candida stellatoidea, and Candida tropicalis). The phylogenetic comparative analysis identifies a hemiascomycete secondary structure consensus that presents a conserved core in all species with variable insertions or deletions. The most significant variability is found in C. glabrata P RNA in which three insertions exceeding in total 700 nt are present in the Specificity domain. This P RNA is more than twice the length of any other homologous P RNAs known in the three domains of life and is eight times the size of the smallest. RNase P RNA, therefore, represents one of the most diversified noncoding RNAs in terms of size variation and structural diversity.

The acyl-AMP ligase FadD32 and AccD4-containing acyl-CoA carboxylase are required for the synthesis of mycolic acids and essential for mycobacterial growth: identification of the carboxylation product and determination of the acyl-CoA carboxylase components.

Mycolic acids are major and specific long-chain fatty acids of the cell envelope of several important human pathogens such as Mycobacterium tuberculosis, M. leprae, and Corynebacterium diphtheriae. Their biosynthesis is essential for mycobacterial growth and represents an attractive target for developing new antituberculous drugs. We have previously shown that the pks13 gene encodes condensase, the enzyme that performs the final condensation step of mycolic acid biosynthesis and is flanked by two genes, fadD32 and accD4. To determine the functions of the gene products we generated two mutants of C. glutamicum with an insertion/deletion within either fadD32 or accD4. The two mutant strains were deficient in mycolic acid production and exhibited the colony morphology that typifies the mycolate-less mutants of corynebacteria. Application of multiple analytical approaches to the analysis of the mutants demonstrated the accumulation of a tetradecylmalonic acid in the DeltafadD32::km mutant and its absence from the DeltaaccD4::km strain. The parental corynebacterial phenotype was restored upon the transfer of the wild-type fadD32 and accD4 genes in the mutants. These data demonstrated that both FadD32 and AccD4-containing acyl-CoA carboxylase are required for the production of mycolic acids. They also prove that the proteins catalyze, respectively, the activation of one fatty acid substrate and the carboxylation of the other substrate, solving the long-debated question of the mechanism involved in the condensation reaction. We used comparative genomics and applied a combination of molecular biology and proteomic technologies to the analysis of proteins that co-immunoprecipitated with AccD4. This resulted in the identification of AccA3 and AccD5 as subunits of the acyl-CoA carboxylase. Finally, we used conditionally replicative plasmids to show that both the fadD32 and accD4 genes are essential for the survival of M. smegmatis. Thus, in addition to Pks13, FadD32 and AccD4 are promising targets for the development of new antimicrobial drugs against pathogenic species of mycobacteria and related microorganisms.

zPicture: dynamic alignment and visualization tool for analyzing conservation profiles.

Comparative sequence analysis has evolved as an essential technique for identifying functional coding and noncoding elements conserved throughout evolution. Here, we introduce zPicture, an interactive Web-based sequence alignment and visualization tool for dynamically generating conservation profiles and identifying evolutionarily conserved regions (ECRs). zPicture is highly flexible, because critical parameters can be modified interactively, allowing users to differentially predict ECRs in comparisons of sequences of different phylogenetic distances and evolutionary rates. We demonstrate the application of this module to identify a known regulatory element in the HOXD locus, in which functional ECRs are difficult to discern against the highly conserved genomic background. zPicture also facilitates transcription factor binding-site analysis via the rVista tool portal. We present an example of the HBB complex when zPicture/rVista combination specifically pinpoints to two ECRs containing GATA-1, NF-E2, and TAL1/E47 binding sites that were identified previously as transcriptional enhancers. In addition, zPicture is linked to the UCSC Genome Browser, allowing users to automatically extract sequences and gene annotations for any recorded locus. Finally, we describe how this tool can be efficiently applied to the analysis of nonvertebrate genomes, including those of microbial organisms.

Isolation of an acetyl-CoA synthetase gene (ZbACS2) from Zygosaccharomyces bailii.

A gene homologous to Saccharomyces cerevisiae ACS genes, coding for acetyl-CoA synthetase, has been cloned from the yeast Zygosaccharomyces bailii ISA 1307, by using reverse genetic approaches. A probe obtained by PCR amplification from Z. bailii DNA, using primers derived from two conserved regions of yeast ACS proteins, RIGAIHSVVF (ScAcs1p; 210-219) and RVDDVVNVSG (ScAcs1p; 574-583), was used for screening a Z. bailii genomic library. Nine clones with partially overlapping inserts were isolated. The sequenced DNA fragment contains a complete ORF of 2027 bp (ZbACS2) and the deduced polypeptide shares significant homologies with the products of ACS2 genes from S. cerevisiae and Kluyveromyces lactis (81% and 82% identity and 84% and 89% similarity, respectively). Phylogenetic analysis shows that the sequence of Zbacs2 is more closely related to the sequences from Acs2 than to those from Acs1 proteins. Moreover, this analysis revealed that the gene duplication producing Acs1 and Acs2 proteins has occurred in the common ancestor of S. cerevisiae, K. lactis, Candida albicans, C. glabrata and Debaryomyces hansenii lineages. Additionally, the cloned gene allowed growth of S. cerevisiae Scacs2 null mutant, in medium containing glucose as the only carbon and energy source, indicating that it encodes a functional acetyl-CoA synthetase. Also, S. cerevisiae cells expressing ZbACS2 have a shorter lag time, in medium containing glucose (2%, w/v) plus acetic acid (0.1-0.35%, v/v). No differences in cell response to acetic acid stress were detected both by specific growth and death rates. The mode of regulation of ZbACS2 appears to be different from ScACS2 and KlACS2, being subject to repression by a glucose pulse in acetic acid-grown cells.

Establishment of a functional genomics platform for Leifsonia xyli subsp. xyli.

Leifsonia xyli subsp. xyli, the causal agent of ratoon stunting disease in sugarcane, is a xylem-limited, nutritionally fastidious, slow growing, gram-positive coryneform bacterium. Because of the difficulties in growing this bacterium in pure culture, little is known about the molecular mechanisms of pathogenesis. Currently, the genome sequence of L. xyli subsp. xyli is being completed by the Agronomical and Environmental Genomes group from the Organization for Nucleotide Sequencing and Analysis in Brazil. To complement this work, we produced 712 Lxx::Tn4431 transposon mutants and sequenced flanking regions from 383 of these, using a rapid polymerase chain reaction-based approach. Tn4431 insertions appeared to be widespread throughout the L. xyli subsp. xyli genome; however, there were regions that had significantly higher concentrations of insertions. The Tn4431 mutant library was screened for individuals unable to colonize sugarcane, and one noncolonizing mutant was found. The mutant contained a transposon insertion disrupting two open reading frames (ORF), one of which had homology to an integral membrane protein from Mycobacterium leprae. Sequencing of the surrounding regions revealed two operons, pro and cyd, both of which are believed to play roles in disease. Complementation studies were carried out using the noncolonizing Lxx::Tn4431 mutant. The noncolonizing mutant was transformed with a cosmid containing 40 kbp of wild-type sequence, which included the two ORF disrupted in the mutant, and several transformants were subsequently able to colonize sugarcane. However, analysis of each of these transformants, before and after colonization, suggests that they have all undergone various recombinant events, obscuring the roles of these ORF in L. xyli subsp. xyli pathogenesis.

Identification of the first eubacterial endonuclease coded by an intein allele in the pps1 gene of mycobacteria.

A survey of a vast range of mycobacterial strains led us to discover a new Pps1 intein allele in Mycobacterium gastri which differs from those of Mycobacterium tuberculosis and Mycobacterium leprae in both its sequence and insertion site. While little is known about Pps1, except that it belongs to the YC24 family of ABC transporters, we show that, unlike the other inteins described so far from Eubacteria, the MgaPps1 intein possesses a specific endonuclease activity. The intein is the first eubacterial intein to be characterised as an endonuclease. Like other intein endonucleases, its minimal sequence for recognition and cleavage is quite large, with 22 bp spanning the Pps1-c site. The fact that an active endonuclease is found among the mycobacterial inteins supports the concept of a cyclical model of invasion by horizontal transfer of these genes, followed by degeneration and loss until a new invasion event, thus explaining their long-term persistence in closely related eubacterial species.

Erp, an extracellular protein family specific to mycobacteria.

Erp (exported repeated protein) was originally characterized as a virulence factor in Mycobacterium tuberculosis and was thought to be present only in Mycobacterium leprae and members of the TB complex. Here it is shown that Erp is a ubiquitous extracellular protein found in all of the mycobacterial species tested. Erp proteins have a modular organization and contain three domains: a highly conserved amino-terminal domain which includes a signal sequence, a central variable region containing repeats based on the motif PGLTS, and a conserved carboxy-terminal domain rich in proline and alanine. The number and fidelity of PGLTS repeats of the central region differ considerably between mycobacterial species. This region is, however, identical in all of the clinical M. tuberculosis strains tested. In addition, it is shown here that a Mycobacterium smegmatis erp::aph mutant displays altered colony morphology which is complemented by all the Erp orthologues tested. The genome sequence flanking the erp gene includes cell-wall-related ORFs and displays extensive conservation between saprophytic and pathogenic mycobacteria.

Unique expression of a highly conserved mycobacterial gene in IS901(+) Mycobacterium avium.

Expression of a gene encoding a novel protein antigen of 40 kDa (p40) was detected in IS901(+) strains of Mycobacterium avium, but not in any other species or subspecies of Mycobacterium tested, including IS901(-) M. avium and the other members of the M. avium complex. Although Southern hybridization revealed that the p40 gene is widely distributed within the genus, expression of the antigen could not be detected on Western blots of mycobacterial cell lysates. Nucleotide sequence analysis of the cloned p40 gene, and a database search, revealed high levels of sequence identity with a homologous gene in IS901(-) M. avium, M. avium subsp. paratuberculosis, Mycobacterium bovis, Mycobacterium leprae, Mycobacterium smegmatis and Mycobacterium tuberculosis. Further analysis of upstream sequences identified a putative promoter region. The p40 gene is the first example of a gene that is widely distributed within the genus Mycobacterium but expressed only in association with the presence of a genomic insertion element, in this case IS901, in strains of M. avium isolated from birds and domestic livestock.

Identification and characterization of a gene encoding a 35-kDa protein from Mycobacterium avium subspecies paratuberculosis.

Mycobacterium avium subspecies paratuberculosis is the causative agent of Johne's disease, a chronic enteritis in ruminants. A gene homologous to that of 35-kDa antigen of Mycobacterium leprae was cloned and sequenced from Mycobacterium paratuberculosis. The database searches revealed 82.79% and 95.67% similarities of its nucleotide sequence, with those of immunodominant 35-kDa protein of M. leprae and M. avium, respectively.

Gene expression, amino acid conservation, and hydrophobicity are the main factors shaping codon preferences in Mycobacterium tuberculosis and Mycobacterium leprae.

Mycobacterium tuberculosis and Mycobacterium leprae are the ethiological agents of tuberculosis and leprosy, respectively. After performing extensive comparisons between genes from these two GC-rich bacterial species, we were able to construct a set of 275 homologous genes. Since these two bacterial species also have a very low growth rate, translational selection could not be so determinant in their codon preferences as it is in other fast-growing bacteria. Indeed, principal-components analysis of codon usage from this set of homologous genes revealed that the codon choices in M. tuberculosis and M. leprae are correlated not only with compositional constraints and translational selection, but also with the degree of amino acid conservation and the hydrophobicity of the encoded proteins. Finally, significant correlations were found between GC3 and synonymous distances as well as between synonymous and nonsynonymous distances.