Highly Reduced Genome of the New Species Mycobacterium uberis, the Causative Agent of Nodular Thelitis and Tuberculoid Scrotitis in Livestock and a Close Relative of the Leprosy Bacilli.

Nodular thelitis is a chronic enzootic infection affecting dairy cows and goats. The causative agent was recently shown to be related to the leprosy-causing bacilli Mycobacterium leprae and Mycobacterium lepromatosis In this study, the genome of this pathogen was sequenced and analyzed. Phylogenomic analyses confirmed that the pathogen present in nodular thelitis and tuberculoid scrotitis is a distinct species related to the leprosy bacilli and Mycobacterium haemophilum Because the pathogen was originally isolated from a bovine udder, it was named "Mycobacterium uberis" The genome of "M. uberis" is only 3.12 Mb in length, which represents the smallest mycobacterial genome identified so far but which is close to that of leprosy bacilli in size. The genome contains 1,759 protein-coding genes and 1,081 pseudogenes, indicative of extensive reductive evolution and likely the reason that M. uberis cannot be grown axenically. The pseudogenization and genome reduction in M. uberis seem to have been to some extent independent from the results determined for the genomes of the leprosy bacilli.IMPORTANCEM. uberis is an emerging skin pathogen in dairy animals. Its genome underwent massive reduction and gene decay, leading to a minimal set of genes required for an obligatory intracellular lifestyle, which highly resembles the evolution of the leprosy agents M. leprae and M. lepromatosis The genomic similarity between M. uberis and the leprosy bacilli can help in identifying key virulence factors of these closely related species or in identifying genes responsible for the distinct differences between thelitis or scrotitis and leprosy with respect to clinical manifestations. Specific DNA markers can now be developed for quick detection of this pathogen.

Analysis of a novel multiplex polymerase chain reaction assay as a sensitive tool for the diagnosis of indeterminate and tuberculoid forms of leprosy.

OBJECTIVE/BACKGROUND: Clinical diagnosis of indeterminate and tuberculoid leprosy is often difficult due to limited and confounding signs and symptoms. In the current study, we evaluated the utility of new multiplex polymerase chain reaction (PCR) using Mycobacterium leprae-specific DNA sequences in the pseudogene regions of ML1545, ML2180, and ML2179 for PCR-based diagnosis of indeterminate leprosy (IND) and leprosy cases across the immunological spectrum. The sensitivity was compared with that of RLEP PCR. METHODS: DNA was extracted from paraffin-embedded skin biopsy specimens of 220 leprosy cases, which were divided into IND (41), tuberculoid form (3), borderline tuberculoid (42), midborderline (3), borderline lepromatous (n=59), and lepromatous leprosy (72) cases. PCR positivity of both multiplex and RLEP PCR were compared in all the samples. A decision tree was constructed using the classification and regression trees algorithm to predict the probability of PCR positivity with the new multiplex PCR scheme in various clinical groups of leprosy. Sensitivity of each pseudogene target was determined using real-time PCR assays, and specificity was confirmed by PCR amplification of DNA extracted from three other mycobacterial species and skin biopsies of 44 non-leprosy cases. RESULTS: A multiplex PCR positivity of 75.61% was noted in IND cases when compared to that of 58.54% using RLEP PCR (P < 0.05). Enhanced multiplex PCR positivity was noted across various clinical groups in comparison to RLEP PCR. The decision tree classifier has predicted statistically significant probability for multiplex PCR positivity among RLEP-PCR negative group and clinical groups with a low bacillary load. CONCLUSION: This new multiplex PCR scheme can support the diagnosis of indeterminate and tuberculoid forms of leprosy with limited clinical manifestations and can be implemented in basic clinical/diagnostic setting that possess conventional PCR facilities.

Whole-genome expression analysis of Mycobacterium leprae and its clinical application.

The whole-genome sequence analysis of Mycobacterium leprae, which was completed in 2001, revealed the characteristics of this microbe's genomic structure. Half of the M. leprae genome consists of a limited number of protein-coding genes and the rest comprises non-coding regions and pseudogenes. We performed membrane array and tiling array analyses to analyze the gene-expression profile of the M. leprae genome and found that pseudogenes and non-coding regions were expressed similarly to coding regions at the RNA level. The RNA expressions were confirmed by real-time PCR analysis. Expression of these RNAs in clinical samples showed varying patterns among patients, thus indicating that the analysis of RNA expression patterns, including non-coding regions and pseudogenes, may be useful for understanding the pathological state, prognosis, and assessment of therapeutic progress in leprosy.

Comparative sequence analysis of Mycobacterium leprae and the new leprosy-causing Mycobacterium lepromatosis.

Mycobacterium lepromatosis is a newly discovered leprosy-causing organism. Preliminary phylogenetic analysis of its 16S rRNA gene and a few other gene segments revealed significant divergence from Mycobacterium leprae, a well-known cause of leprosy, that justifies the status of M. lepromatosis as a new species. In this study we analyzed the sequences of 20 genes and pseudogenes (22,814 nucleotides). Overall, the level of matching of these sequences with M. leprae sequences was 90.9%, which substantiated the species-level difference; the levels of matching for the 16S rRNA genes and 14 protein-encoding genes were 98.0% and 93.1%, respectively, but the level of matching for five pseudogenes was only 79.1%. Five conserved protein-encoding genes were selected to construct phylogenetic trees and to calculate the numbers of synonymous substitutions (dS values) and nonsynonymous substitutions (dN values) in the two species. Robust phylogenetic trees constructed using concatenated alignment of these genes placed M. lepromatosis and M. leprae in a tight cluster with long terminal branches, implying that the divergence occurred long ago. The dS and dN values were also much higher than those for other closest pairs of mycobacteria. The dS values were 14 to 28% of the dS values for M. leprae and Mycobacterium tuberculosis, a more divergent pair of species. These results thus indicate that M. lepromatosis and M. leprae diverged approximately 10 million years ago. The M. lepromatosis pseudogenes analyzed that were also pseudogenes in M. leprae showed nearly neutral evolution, and their relative ages were similar to those of M. leprae pseudogenes, suggesting that they were pseudogenes before divergence. Taken together, the results described above indicate that M. lepromatosis and M. leprae diverged from a common ancestor after the massive gene inactivation event described previously for M. leprae.

Whole-genome tiling array analysis of Mycobacterium leprae RNA reveals high expression of pseudogenes and noncoding regions.

Whole-genome sequence analysis of Mycobacterium leprae has revealed a limited number of protein-coding genes, with half of the genome composed of pseudogenes and noncoding regions. We previously showed that some M. leprae pseudogenes are transcribed at high levels and that their expression levels change following infection. In order to clarify the RNA expression profile of the M. leprae genome, a tiling array in which overlapping 60-mer probes cover the entire 3.3-Mbp genome was designed. The array was hybridized with M. leprae RNA from the SHR/NCrj-rnu nude rat, and the results were compared to results from an open reading frame array and confirmed by reverse transcription-PCR. RNA expression was detected from genes, pseudogenes, and noncoding regions. The signal intensities obtained from noncoding regions were higher than those from pseudogenes. Expressed noncoding regions include the M. leprae unique repetitive sequence RLEP and other sequences without any homology to known functional noncoding RNAs. Although the biological functions of RNA transcribed from M. leprae pseudogenes and noncoding regions are not known, RNA expression analysis will provide insights into the bacteriological significance of the species. In addition, our study suggests that M. leprae will be a useful model organism for the study of the molecular mechanism underlying the creation of pseudogenes and the role of microRNAs derived from noncoding regions.

[Detection of RNA expression on whole genome analysis of Mycobacterium leprae by tiling array].

Completion of Mycobacterium leprae genome sequence revealed that there are many pseudogenes and non-coding regions, but rather small numbers of protein-coding genes. Although it was thought that pseudogenes and non-coding regions were silent and junk, our previous studies indicated that RNA expression was detected from these regions. To elucidate comprehensive RNA expression pattern on M. leprae whole genome, tiling array was designed and total RNA of M. leprae Thai-53 strain was analyzed. As a result, highly expressed regions were detected among not only the gene regions but also pseudogenes and non-coding regions. Since some of the RNA expression levels were modulated by MDT, evaluation of RNA expression pattern might be a good indicator for the treatment of leprosy.

[Comprehensive analysis of RNA expression of Mycobacterium leprae and clinical and biological significance].

Completion of Mycobacterium leprae genome sequence revealed that there are many pseudogenes and non-coding regions, but rather small numbers of protein-coding genes. This result indicates that M. leprae is a very unique organism, and this future is important to understand the biological nature and/or pathogenicity of M. leprae, which remain unclear. We attempted to find the biological nature of M. leprae by detecting the gene and pseudogene regions transcribed at high level. We detected the genomic regions including pseudogenes and demonstrated that six out of twelve high expression regions were pseudogenes. In addition, its transcription level was changed when M. leprae infects macrophage. RNA was detected from genes, pseudogenes and non-coding regions. The expression levels of these regions were different among patients and a part of them is disappeared just after treatment. These results suggested that RNA derived from pseudogene and non-coding region have some function concerning the infection and/or intracellular parasitism and that the analysis of pseudogene and non-coding region expression pattern of M. leprae is available as a criterion for therapeutic effect and disease type of leprosy, and a prognostic marker.

Where are the pseudogenes in bacterial genomes?

Most bacterial genomes have very few pseudogenes; notable exceptions include the genomes of the intracellular parasites Rickettsia prowazekii and Mycobacterium leprae. As DNA can be introduced into microbial genomes in many ways, the compact nature of these genomes suggests that the rate of DNA influx is balanced by the rate of DNA deletion. We propose that the influx of dangerous genetic elements such as transposons and bacteriophages selects for the maintenance of relatively high deletion rates in most bacteria; the sheltered lifestyle of intracellular parasites removes this threat, leading to reduced deletion rates and larger pseudogene loads.

On the catalase-peroxidase gene, katG, of Mycobacterium leprae and the implications for treatment of leprosy with isoniazid.

The toxicity of the potent tuberculocidal agent, isoniazid, is mediated by the heme-containing enzyme, catalase-peroxidase, encoded by the katG gene. Although isoniazid has been used for the treatment of leprosy, it is shown here that the katG gene of Mycobacterium leprae is a pseudogene, which has probably been inactivated by multiple mutations. Inactive genes were detected by the polymerase chain reaction in several isolates of M. leprae, of different geographical origins, and attempts to complement an isoniazid-resistant strain of Mycobacterium smegmatis with the katG pseudogene were unsuccessful. Isoniazid is thus likely to be of no therapeutic benefit to leprosy patients.