A phylogenomic study quantifies competing mechanisms for pseudogenization in prokaryotes-The Mycobacterium leprae case.

BACKGROUND: Pseudogenes are non-functional sequences in the genome with homologous sequences that are functional (i.e. genes). They are abundant in eukaryotes where they have been extensively investigated, while in prokaryotes they are significantly scarcer and less well studied. Here we conduct a comprehensive analysis of the evolution of orthologs of Mycobacterium leprae pseudogenes in prokaryotes. The leprosy pathogen M. leprae is of particular interest since it contains an unusually large number of pseudogenes, comprising approximately 40% of its entire genome. The analysis is conducted in both broad and narrow phylogenetic ranges. RESULTS: We have developed an informatics-based approach to characterize the evolution of pseudogenes. This approach combines tools from phylogenomics, genomics, and transcriptomics. The results we obtain are used to assess the contributions of two mechanisms for pseudogene formation: failed horizontal gene transfer events and disruption of native genes. CONCLUSIONS: We conclude that, although it was reported that in most bacteria the former is most likely responsible for the majority of pseudogenization events, in mycobacteria, and in particular in M. leprae with its exceptionally high pseudogene numbers, the latter predominates. We believe that our study sheds new light on the evolution of pseudogenes in bacteria, by utilizing new methodologies that are applied to the unusually abundant M. leprae pseudogenes and their orthologs.

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.

Identification of novel open reading frames in the intergenic regions of Mycobacterium leprae genome and detection of transcript by qRT-PCR.

Mycobacterium leprae is an unculturable obligate pathogen and causative agent for debilitating human disease leprosy. Due to reductive genome evolution M leprae genome harbours large number of pseudogenes and small number of genes (∼1600 genes and ∼1300 pseudogenes). How M leprae remained a successful human parasite with small set of genes remains poorly understood and provided us the impetus to investigate the intergenic regions of M leprae genome for the presence of possible open reading frames (ORFs). In this work, we have manually scanned all the intergenic regions of M leprae genome and identified 106 potential ORFs. Among these, 12 are large ORFs: encoding hypothetical proteins (HP) of more than 100 amino acids. We have also found 67 ORFs encoding 50-100 amino acids proteins and another 27 ORFs for 30-50 amino acids peptides. We have validated the presence of transcripts for large HPs by quantitative reverse transcriptase PCR (qRT-PCR). Our results suggest that some of the M leprae large HPs are indeed expressed at low level in leprosy patients. The present results will shed light on the intergenic ORFs of M leprae and further our understanding of the pathogenesis of leprosy.

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.

Mycobacterium leprae: genes, pseudogenes and genetic diversity.

Leprosy, which has afflicted human populations for millenia, results from infection with Mycobacterium leprae, an unculturable pathogen with an exceptionally long generation time. Considerable insight into the biology and drug resistance of the leprosy bacillus has been obtained from genomics. M. leprae has undergone reductive evolution and pseudogenes now occupy half of its genome. Comparative genomics of four different strains revealed remarkable conservation of the genome (99.995% identity) yet uncovered 215 polymorphic sites, mainly single nucleotide polymorphisms, and a handful of new pseudogenes. Mapping these polymorphisms in a large panel of strains defined 16 single nucleotide polymorphism-subtypes that showed strong geographical associations and helped retrace the evolution of M. leprae.

The pseudogenes of Mycobacterium leprae reveal the functional relevance of gene order within operons.

Almost 50 years following the discovery of the prokaryotic operon, the functional relevance of gene order within operons remains unclear. In this work, we take advantage of the eroded genome of Mycobacterium leprae to add evidence supporting the notion that functionally less important genes have a tendency to be located at the end of its operons. M. leprae's genome includes 1133 pseudogenes and 1614 protein-coding genes and can be compared with the close genome of M. tuberculosis. Assuming M. leprae's pseudogenes to represent dispensable genes, we have studied the position of these pseudogenes in the operons of M. leprae and of their orthologs in M. tuberculosis. We observed that both tend to be located in the 3' (downstream) half of the operon (P-values of 0.03 and 0.18, respectively). Analysis of pseudogenes in all available prokaryotic genomes confirms this trend (P-value of 7.1 × 10(-7)). In a complementary analysis, we found a significant tendency for essential genes to be located at the 5' (upstream) half of the operon (P-value of 0.006). Our work provides an indication that, in prokarya, functionally less important genes have a tendency to be located at the end of operons, while more relevant genes tend to be located toward operon starts.

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.

Implications of high level pseudogene transcription in Mycobacterium leprae.

BACKGROUND: The Mycobacterium leprae genome has less than 50% coding capacity and 1,133 pseudogenes. Preliminary evidence suggests that some pseudogenes are expressed. Therefore, defining pseudogene transcriptional and translational potentials of this genome should increase our understanding of their impact on M. leprae physiology. RESULTS: Gene expression analysis identified transcripts from 49% of all M. leprae genes including 57% of all ORFs and 43% of all pseudogenes in the genome. Transcribed pseudogenes were randomly distributed throughout the chromosome. Factors resulting in pseudogene transcription included: 1) co-orientation of transcribed pseudogenes with transcribed ORFs within or exclusive of operon-like structures; 2) the paucity of intrinsic stem-loop transcriptional terminators between transcribed ORFs and downstream pseudogenes; and 3) predicted pseudogene promoters. Mechanisms for translational "silencing" of pseudogene transcripts included the lack of both translational start codons and strong Shine-Dalgarno (SD) sequences. Transcribed pseudogenes also contained multiple "in-frame" stop codons and high Ka/Ks ratios, compared to that of homologs in M. tuberculosis and ORFs in M. leprae. A pseudogene transcript containing an active promoter, strong SD site, a start codon, but containing two in frame stop codons yielded a protein product when expressed in E. coli. CONCLUSION: Approximately half of M. leprae's transcriptome consists of inactive gene products consuming energy and resources without potential benefit to M. leprae. Presently it is unclear what additional detrimental affect(s) this large number of inactive mRNAs has on the functional capability of this organism. Translation of these pseudogenes may play an important role in overall energy consumption and resultant pathophysiological characteristics of M. leprae. However, this study also demonstrated that multiple translational "silencing" mechanisms are present, reducing additional energy and resource expenditure required for protein production from the vast majority of these transcripts.

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.

Detection of RNA expression from pseudogenes and non-coding genomic regions of Mycobacterium leprae.

We have previously reported that some pseudogenes are expressed in Mycobacterium leprae (M. leprae), the causative agent of leprosy, and that their expression levels alter upon infection of macrophages. We attempted to further examine the expression of pseudogene and non-coding genomic region in M. leprae, in this study. 19 Pseudogenes, 17 non-coding genomic regions, and 21 coding genes expression in M. leprae maintained in the footpads of the hypertensive nude rat (SHR/NCrj-rnu) were examined by reverse transcriptase polymerase chain reaction (RT-PCR). The expression of some of these pseudogenes, non-coding genomic regions and coding genes were also examined in M. leprae from skin smear specimens obtained from patients with lepromatous leprosy by RT-PCR. Transcripts from pseudogenes, non-coding genomic regions and coding genes examined in this study were clearly observed in M. leprae. The expression patterns of some of these transcripts vary greatly among different leprosy patients. These results indicate that some of pseudogenes and non-coding genomic regions are transcribed in M. leprae and analysis of RNA expression patterns including pseudogene and non-coding genomic region in M. leprae may be useful in understanding the pathological states of infected patients.

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.

Long simple sequence repeats in host-adapted pathogens localize near genes encoding antigens, housekeeping genes, and pseudogenes.

Simple sequence repeats (SSRs) in DNA sequences are tandem iterations of a single nucleotide or a short oligonucleotide. SSRs are subject to slipped-strand mutations and a common source of phase variation in bacteria and antigenic variation in pathogens. Significantly long SSRs are generally rare in prokaryotic genomes, and long SSRs composed of iterations of mono-, di-, tri-, and tetranucleotides are mostly restricted to host-adapted pathogens. We present new results concerning associations between long SSRs and genes related to different cellular functions in genomes of host-adapted pathogens. We found that in the majority of the analyzed genomes, at least some of the genes associated with SSRs encode potential antigens, which is expected if the primary function of SSRs is their contribution to antigenic variation. However, we also found a number of long SSRs associated with housekeeping genes, including rRNA and tRNA genes, genes encoding ribosomal proteins, amino acyl-tRNA synthetases, chaperones, and important metabolic enzymes. Many of these genes are probably essential and it is unlikely that they are phase-variable. Few statistically significant associations between SSRs and gene functional classifications were detected, suggesting that most long SSRs are not related to a particular cellular function or process. Long SSRs in Mycobacterium leprae are mostly associated with pseudogenes and may be contributing to gene loss following the adaptation to an obligate pathogenic lifestyle. We speculate that LSSRs may have played a similar role in genome reduction of other host-adapted pathogens.

Promiscuous DNA in the nuclear genomes of hemiascomycetous yeasts.

Abstract Transfer of fragments of mtDNA to the nuclear genome is a general phenomenon that gives rise to NUMTs (NUclear sequences of MiTochondrial origin). We present here the first comparative analysis of the NUMT content of entirely sequenced species belonging to a monophyletic group, the hemiascomycetous yeasts (Candida glabrata, Kluyveromyces lactis, Kluyveromyces thermotolerans, Debaryomyces hansenii and Yarrowia lipolytica, along with the updated NUMT content of Saccharomyces cerevisiae). This study revealed a huge diversity in NUMT number and organization across the six species. Debaryomyces hansenii harbors the highest number of NUMTs (145), half of which are distributed in numerous large mosaics of up to eight NUMTs arising from multiple noncontiguous mtDNA fragments inserted at the same chromosomal locus. Most NUMTs, in all species, are found within intergenic regions including seven NUMTs in pseudogenes. However, five NUMTs overlap a gene, suggesting a positive impact of NUMTs on protein evolution. Contrary to the other species, K. lactis and K. thermotolerans harbor only a few diverged NUMTs, suggesting that mitochondrial transfer to the nuclear genome has decreased or ceased in these phylogenetic branches. The dynamics of NUMT acquisition and loss are illustrated here by their species-specific distribution.

Visualization of pseudogenes in intracellular bacteria reveals the different tracks to gene destruction.

BACKGROUND: Pseudogenes reveal ancestral gene functions. Some obligate intracellular bacteria, such as Mycobacterium leprae and Rickettsia spp., carry substantial fractions of pseudogenes. Until recently, horizontal gene transfers were considered to be rare events in obligate host-associated bacteria. RESULTS: We present a visualization tool that displays the relationships and positions of degraded and partially overlapping gene sequences in multiple genomes. With this tool we explore the origin and deterioration patterns of the Rickettsia pseudogenes and find that variably present genes and pseudogenes tend to have been acquired more recently, are more divergent in sequence, and exhibit a different functional profile compared with genes conserved across all species. Overall, the origin of only one-quarter of the variable genes and pseudogenes can be traced back to the common ancestor of Rickettsia and the outgroup genera Orientia and Wolbachia. These sequences contain only a few disruptive mutations and show a broad functional distribution profile, much like the core genes. The remaining genes and pseudogenes are extensively degraded or solely present in a single species. Their functional profile was heavily biased toward the mobile gene pool and genes for components of the cell wall and the lipopolysaccharide. CONCLUSION: Reductive evolution of the vertically inherited genomic core accounts for 25% of the predicted genes in the variable segments of the Rickettsia genomes, whereas 75% stems from the flux of the mobile gene pool along with genes for cell surface structures. Thus, most of the variably present genes and pseudogenes in Rickettsia have arisen from recent acquisitions.

Molecular basis of the defective heat stress response in Mycobacterium leprae.

Mycobacterium leprae, a major human pathogen, grows poorly at 37 degrees C. The basis for its inability to survive at elevated temperatures was investigated. We determined that M. leprae lacks a protective heat shock response as a result of the lack of transcriptional induction of the alternative sigma factor genes sigE and sigB and the major heat shock operons, HSP70 and HSP60, even though heat shock promoters and regulatory circuits for these genes appear to be intact. M. leprae sigE was found to be capable of complementing the defective heat shock response of mycobacterial sigE knockout mutants only in the presence of a functional mycobacterial sigH, which orchestrates the mycobacterial heat shock response. Since the sigH of M. leprae is a pseudogene, these data support the conclusion that a key aspect of the defective heat shock response in M. leprae is the absence of a functional sigH. In addition, 68% of the genes induced during heat shock in M. tuberculosis were shown to be either absent from the M. leprae genome or were present as pseudogenes. Among these is the hsp/acr2 gene, whose product is essential for M. tuberculosis survival during heat shock. Taken together, these results suggest that the reduced ability of M. leprae to survive at elevated temperatures results from the lack of a functional transcriptional response to heat shock and the absence of a full repertoire of heat stress response genes, including sigH.

High-level expression of pseudogenes in Mycobacterium leprae.

Recent studies have revealed that some RNAs are transcribed from noncoding DNA regions, including pseudogenes, and are functional as riboregulators. We have attempted to assess the gene expression profile throughout the Mycobacterium leprae genome using an array technique. Twelve highly expressed gene regions were identified that show an alteration in expression levels upon infection. Six of these were pseudogenes. Although M. leprae has an exceptional number and proportion of pseudogenes among species, our results suggest that some of the M. leprae pseudogenes are not just 'decayed' genes, but may have a functional role.

Inferring the pattern of spontaneous mutation from the pattern of substitution in unitary pseudogenes of Mycobacterium leprae and a comparison of mutation patterns among distantly related organisms.

The pattern of spontaneous mutation can be inferred from the pattern of substitution in pseudogenes, which are known to be under very weak or no selective constraint. We modified an existing method (Gojobori T, et al., J Mol Evol 18:360, 1982) to infer the pattern of mutation in bacteria by using 569 pseudogenes from Mycobacterium leprae. In Gojobori et al.'s method, the pattern is inferred by using comparisons involving a pseudogene, a conspecific functional paralog, and an outgroup functional ortholog. Because pseudogenes in M. leprae are unitary, we replaced the missing paralogs by functional orthologs from M. tuberculosis. Functional orthologs from Streptomyces coelicolor served as outgroups. We compiled a database consisting of 69,378 inferred mutations. Transitional mutations were found to constitute more than 56% of all mutations. The transitional bias was mainly due to C-->T and G-->A, which were also the most frequent mutations on the leading strand and the only ones that were significantly more frequent than the random expectation. The least frequent mutations on the leading strand were A-->T and T-->A, each with a relative frequency of less than 3%. The mutation pattern was found to differ between the leading and the lagging strands. This asymmetry is thought to be the cause for the typical chirochoric structure of bacterial genomes. The physical distance of the pseudogene from the origin of replication (ori) was found to have almost no effect on the pattern of mutation. A surprising similarity was found between the mutation pattern in M. leprae and previously inferred patterns for such distant taxa as human and Drosophila. The mutation pattern on the leading strand of M. leprae was also found to share some common features with the pattern inferred for the heavy strand of the human mitochondrial genome. These findings indicate that taxon-specific factors may only play secondary roles in determining patterns of mutation.

On the origin of leprosy.

Leprosy, a chronic human disease with potentially debilitating neurological consequences, results from infection with Mycobacterium leprae. This unculturable pathogen has undergone extensive reductive evolution, with half of its genome now occupied by pseudogenes. Using comparative genomics, we demonstrated that all extant cases of leprosy are attributable to a single clone whose dissemination worldwide can be retraced from analysis of very rare single-nucleotide polymorphisms. The disease seems to have originated in Eastern Africa or the Near East and spread with successive human migrations. Europeans or North Africans introduced leprosy into West Africa and the Americas within the past 500 years.