High-level expression of pseudogenes in Mycobacterium leprae

Pathogenicity and virulence of Mycobacterium leprae

Leprosy is caused by Mycobacterium leprae (M. leprae) and M. lepromatosis, an obligate intracellular organism, and over 200,000 new cases occur every year. M. leprae parasitizes histiocytes (skin macrophages) and Schwann cells in the peripheral nerves. Although leprosy can be treated by multidrug therapy, some patients relapse or have a prolonged clinical course and/or experience leprosy reaction. These varying outcomes depend on host factors such as immune responses against bacterial components that determine a range of symptoms. To understand these host responses, knowledge of the mechanisms by which M. leprae parasitizes host cells is important. This article describes the characteristics of leprosy through bacteriology, genetics, epidemiology, immunology, animal models, routes of infection, and clinical findings. It also discusses recent diagnostic methods, treatment, and measures according to the World Health Organization (WHO), including prevention. Recently, the antibacterial activities of anti-hyperlipidaemia agents against other pathogens, such as M. tuberculosis and Staphylococcus aureus have been investigated. Our laboratory has been focused on the metabolism of lipids which constitute the cell wall of M. leprae. Our findings may be useful for the development of future treatments.

Rv1915 and Rv1916 from Mycobacterium tuberculosis H37Rv form in vitro protein-protein complex

BACKGROUND

Mycobacterium tuberculosis (Mtb) isocitrate lyase (ICL) is an established drug target that facilitates Mtb persistence. Unlike other mycobacterial strains, where ICL2 is a single gene product, H37Rv has a split event, resulting in two tandemly coded icls - rv1915 and rv1916. Our recent report on functionality of individual Rv1915 and Rv1916, led to postulate the cooperative role of these proteins in pathogen's survival under nutrient-limiting conditions. This study investigates the possibility of Rv1915 and Rv1916 interacting and forming a complex.

METHODS

Pull down assay, activity assay, mass spectrometry and site directed mutagenesis was employed to investigate and validate Rv1915-Rv1916 complex formation.

RESULTS

Rv1915 and Rv1916 form a stable complex in vitro, with enhanced ICL/MICL activities as opposed to individual proteins. Further, activities monitored in the presence of acetyl-CoA show significant increase for Rv1916 and the complex but not of Rv0467 and Rv1915Δ90CT. Both full length and truncated Rv1915Δ90CT can form complex, implying the absence of its C-terminal disordered region in complex formation. Further, in silico analysis and site-directed mutagenesis studies reveal Y64 and Y65 to be crucial residues for Rv1915-Rv1916 complex formation.

CONCLUSIONS

This study uncovers the association between Rv1915 and Rv1916 and supports the role of acetyl-CoA in escalating the ICL/MICL activities of Rv1916 and Rv1915Δ90CT-Rv1916 complex.

GENERAL SIGNIFICANCE

Partitioning of ICL2 into Rv1915 and Rv1916 that associates to form a complex in Mtb H37Rv, suggests its importance in signaling and regulation of metabolic pathway particularly in carbon assimilation.

Mycobacterium leprae promotes triacylglycerol de novo synthesis through induction of GPAT3 expression in human premonocytic THP-1 cells

Mycobacterium leprae (M. leprae) is the etiological agent of leprosy, and the skin lesions of lepromatous leprosy are filled with numerous foamy or xanthomatous histiocytes that are parasitized by M. leprae. Lipids are an important nutrient for the intracellular survival of M. leprae. In this study, we attempted to determine the intracellular lipid composition and underlying mechanisms for changes in host cell lipid metabolism induced by M. leprae infection. Using high-performance thin-layer chromatography (HPTLC), we demonstrated specific induction of triacylglycerol (TAG) production in human macrophage THP-1 cells following M. leprae infection. We then used [14C] stearic acid tracing to show incorporation of this newly synthesized host cell TAG into M. leprae. In parallel with TAG accumulation, expression of host glycerol-3-phosphate acyltransferase 3 (GPAT3), a key enzyme in de novo TAG synthesis, was significantly increased in M. leprae-infected cells. CRISPR/Cas9 genome editing of GPAT3 in THP-1 cells (GPAT3 KO) dramatically reduced accumulation of TAG following M. leprae infection, intracellular mycobacterial load, and bacteria viability. These results together suggest that M. leprae induces host GPAT3 expression to facilitate TAG accumulation within macrophages to maintain a suitable environment that is crucial for intracellular survival of these bacilli.

Pathogenesis and Host Immune Response in Leprosy

Leprosy is an ancient insidious disease caused by Mycobacterium leprae, where the skin and peripheral nerves undergo chronic granulomatous infections, leading to sensory and motor impairment with characteristic deformities. Susceptibility to leprosy and its disease state are determined by the manifestation of innate immune resistance mediated by cells of monocyte lineage. Due to insufficient innate resistance, granulomatous infection is established, influencing the specific cellular immunity. The clinical presentation of leprosy ranges between two stable polar forms (tuberculoid to lepromatous) and three unstable borderline forms. The tuberculoid form involves Th1 response, characterized by a well demarcated granuloma, infiltrated by CD4⁺ T lymphocytes, containing epitheloid and multinucleated giant cells. In the lepromatous leprosy, there is no characteristic granuloma but only unstructured accumulation of ineffective macrophages containing engulfed pathogens. Th1 response, characterised by IFN-γ and IL-2 production, activates macrophages in order to kill intracellular pathogens. Conversely, a Th2 response, characterized by the production of IL-4, IL-5 and IL-10, helps in antibody production and consequently downregulates the cell-mediated immunity induced by the Th1 response. M. lepare has a long generation time and its inability to grow in culture under laboratory conditions makes its study challenging. The nine-banded armadillo still remains the best clinical and immunological model to study host-pathogen interaction in leprosy. In this chapter, we present cellular morphology and the genomic uniqueness of M. leprae, and how the pathogen shows tropism for Schwann cells, macrophages and dendritic cells.

Complete genome sequence and annotation of the laboratory reference strain Shigella flexneri serotype 5a M90T and genome-wide transcriptional start site determination

Background

Shigella is a Gram-negative facultative intracellular bacterium that causes bacillary dysentery in humans. Shigella invades cells of the colonic mucosa owing to its virulence plasmid-encoded Type 3 Secretion System (T3SS), and multiplies in the target cell cytosol. Although the laboratory reference strain S. flexneri serotype 5a M90T has been extensively used to understand the molecular mechanisms of pathogenesis, its complete genome sequence is not available, thereby greatly limiting studies employing high-throughput sequencing and systems biology approaches.

Results

We have sequenced, assembled, annotated and manually curated the full genome of S. flexneri 5a M90T. This yielded two complete circular contigs, the chromosome and the virulence plasmid (pWR100). To obtain the genome sequence, we have employed long-read PacBio DNA sequencing followed by polishing with Illumina RNA-seq data. This provides a new hybrid strategy to prepare gapless, highly accurate genome sequences, which also cover AT-rich tracks or repetitive sequences that are transcribed. Furthermore, we have performed genome-wide analysis of transcriptional start sites (TSS) and determined the length of 5′ untranslated regions (5′-UTRs) at typical culture conditions for the inoculum of in vitro infection experiments. We identified 6723 primary TSS (pTSS) and 7328 secondary TSS (sTSS). The S. flexneri 5a M90T annotated genome sequence and the transcriptional start sites are integrated into RegulonDB (http://regulondb.ccg.unam.mx) and RSAT (http://embnet.ccg.unam.mx/rsat/) databases to use their analysis tools in the S. flexneri 5a M90T genome.

Conclusions

We provide the first complete genome for S. flexneri serotype 5a, specifically the laboratory reference strain M90T. Our work opens the possibility of employing S. flexneri M90T in high-quality systems biology studies such as transcriptomic and differential expression analyses or in genome evolution studies. Moreover, the catalogue of TSS that we report here can be used in molecular pathogenesis studies as a resource to know which genes are transcribed before infection of host cells. The genome sequence, together with the analysis of transcriptional start sites, is also a valuable tool for precise genetic manipulation of S. flexneri 5a M90T. Further, we present a new hybrid strategy to prepare gapless, highly accurate genome sequences. Unlike currently used hybrid strategies combining long- and short-read DNA sequencing technologies to maximize accuracy, our workflow using long-read DNA sequencing and short-read RNA sequencing provides the added value of using non-redundant technologies, which yield distinct, exploitable datasets.

Mycobacterium leprae's evolution and environmental adaptation

Leprosy is an ancient disease caused by the acid-fast bacillus Mycobacterium leprae, also known as Hansen's bacillus. M. leprae is an obligate intracellular microorganism with a marked Schwann cell tropism and is the only human pathogen capable of invading the superficial peripheral nerves. The transmission mechanism of M. leprae is not fully understood; however, the nasal mucosa is accepted as main route of M. leprae entry to the human host. The complete sequencing and the comparative genome analysis show that M. leprae underwent a genome reductive evolution process, as result of lifestyle change and adaptation to different environments; some of lost genes are homologous to those of host cells. Thus, M. leprae reduced its genome size to 3.3 Mbp, contributing to obtain the lowest GC content (approximately 58%) among mycobacteria. The M. leprae genome contains 1614 open reading frames coding for functional proteins, and 1310 pseudogenes corresponding to 41% of the genome, approximately. Comparative analyses to different microorganisms showed that M. leprae possesses the highest content of pseudogenes among pathogenic and non-pathogenic bacteria and archaea. The pathogen adaptation into host cells, as the Schwann cells, brought about the reduction of the genome and induced multiple gene inactivation. The present review highlights the characteristics of genome's reductive evolution that M. leprae experiences in the genetic aspects compared with other pathogens. The possible mechanisms of pseudogenes formation are discussed.

Novel coordination of lipopolysaccharide modifications in Vibrio cholerae promotes CAMP resistance

In the environment and during infection, the human intestinal pathogen Vibrio cholerae must overcome noxious compounds that damage the bacterial outer membrane. The El Tor and classical biotypes of O1 V. cholerae show striking differences in their resistance to membrane disrupting cationic antimicrobial peptides (CAMPs), such as polymyxins. The classical biotype is susceptible to CAMPs, but current pandemic El Tor biotype isolates gain CAMP resistance by altering the net charge of their cell surface through glycine modification of lipid A. Here we report a second lipid A modification mechanism that only functions in the V. cholerae El Tor biotype. We identify a functional EptA ortholog responsible for the transfer of the amino-residue phosphoethanolamine (pEtN) to the lipid A of V. cholerae El Tor that is not functional in the classical biotype. We previously reported that mildly acidic growth conditions (pH 5.8) downregulate expression of genes encoding the glycine modification machinery. In this report, growth at pH 5.8 increases expression of eptA with concomitant pEtN modification suggesting coordinated regulation of these LPS modification systems. Similarly, efficient pEtN lipid A substitution is seen in the absence of lipid A glycinylation. We further demonstrate EptA orthologs from non-cholerae Vibrio species are functional.

Decoding the similarities and differences among mycobacterial species

Mycobacteriaceae comprises pathogenic species such as Mycobacterium tuberculosis, M. leprae and M. abscessus, as well as non-pathogenic species, for example, M. smegmatis and M. thermoresistibile. Genome comparison and annotation studies provide insights into genome evolutionary relatedness, identify unique and pathogenicity-related genes in each species, and explore new targets that could be used for developing new diagnostics and therapeutics. Here, we present a comparative analysis of ten-mycobacterial genomes with the objective of identifying similarities and differences between pathogenic and non-pathogenic species. We identified 1080 core orthologous clusters that were enriched in proteins involved in amino acid and purine/pyrimidine biosynthetic pathways, DNA-related processes (replication, transcription, recombination and repair), RNA-methylation and modification, and cell-wall polysaccharide biosynthetic pathways. For their pathogenicity and survival in the host cell, pathogenic species have gained specific sets of genes involved in repair and protection of their genomic DNA. M. leprae is of special interest owing to its smallest genome (1600 genes and ~1300 psuedogenes), yet poor genome annotation. More than 75% of the pseudogenes were found to have a functional ortholog in the other mycobacterial genomes and belong to protein families such as transferases, oxidoreductases and hydrolases.

Members of the Mycobacteriaceae family, which are known to adapt to different environmental niches, comprise bacterial species with varied genome sizes. They are unique in their cell-wall composition, which is remarkably thick and lipid-rich as compared to other bacteria. We performed a comparative analysis at the proteome level for ten mycobacterial species that differ in their pathogenicity, genome size and environmental niches. A total of 1080 orthologous clusters with representation from all ten species were obtained, and these were further examined for their domain annotations, domain architecture similarities and enriched GO terms. These core orthologous clusters are enriched in various biosynthetic pathways. The proteins that are specific to each of the ten species were also investigated for their GO functions. The M. leprae genome has a large number of pseudogenes and we searched for their functional orthologs in other mycobacterial species in order to understand the functions that are lost from the M. leprae genome. The proteins present exclusively in M. leprae genome were studied in more detail, in order to predict putative drug targets and diagnostic markers. These findings, which have implications in understanding evolution of mycobacterial genomes, identify species-specific proteins that have potential for use in developing new diagnostic tools and therapeutics.

Strand-specific RNA-seq analysis of the Lactobacillus delbrueckii subsp. bulgaricus transcriptome

Lactobacillus delbrueckii subsp. bulgaricus 2038 (Lb. bulgaricus 2038) is an industrial bacterium that is used as a starter for dairy products. We proposed several hypotheses concerning its industrial features previously. Here, we utilized RNA-seq to explore the transcriptome of Lb. bulgaricus 2038 from four different growth phases under whey conditions. The most abundantly expressed genes in the four stages were mainly involved in translation (for the logarithmic stage), glycolysis (for control/lag stages), lactic acid production (all the four stages), and 10-formyl tetrahydrofolate production (for the stationary stage). The high expression of genes like d-lactate dehydrogenase was thought as a result of energy production, and consistent expression of EPS synthesis genes, the restriction-modification (RM) system and the CRISPR/Cas system were validated for explaining the advantage of this strain in yoghurt production. Several postulations, like NADPH production through GapN bypass, converting aspartate into carbon-skeleton intermediates, and formate production through degrading GTP, were proved not working under these culture conditions. The high expression of helicase genes and co-expressed amino acids/oligopeptides transporting proteins indicated that the helicase might mediate the strain obtaining nitrogen source from the environment. The transport system of Lb. bulgaricus 2038 was found to be regulated by antisense RNA, hinting the potential application of non-coding RNA in regulating lactic acid bacteria (LAB) gene expression. Our study has primarily uncovered Lb. bulgaricus 2038 transcriptome, which could gain a better understanding of the regulation system in Lb. bulgaricus and promote its industrial application.

To Be or Not to Be a Pseudogene: A Molecular Epidemiological Approach to the mclx Genes and Its Impact in Tuberculosis

Tuberculosis presents a myriad of symptoms, progression routes and propagation patterns not yet fully understood. Whereas for a long time research has focused solely on the patient immunity and overall susceptibility, it is nowadays widely accepted that the genetic diversity of its causative agent, Mycobacterium tuberculosis, plays a key role in this dynamic. This study focuses on a particular family of genes, the mclxs (Mycobacteriumcyclase/LuxR-like genes), which codify for a particular and nearly mycobacterial-exclusive combination of protein domains. mclxs genes were found to be pseudogenized by frameshift-causing insertion(s)/deletion(s) in a considerable number of M. tuberculosis complex strains and clinical isolates. To discern the functional implications of the pseudogenization, we have analysed the pattern of frameshift-causing mutations in a group of M. tuberculosis isolates while taking into account their microbial-, patient- and disease-related traits. Our logistic regression-based analyses have revealed disparate effects associated with the transcriptional inactivation of two mclx genes. In fact, mclx2 (Rv1358) pseudogenization appears to be primarily driven by the microbial phylogenetic background, being mainly related to the Euro-American (EAm) lineage; on the other hand, mclx3 (Rv2488c) presents a higher tendency for pseudogenization among isolates from patients born on the Western Pacific area, and from isolates causing extra-pulmonary infections. These results contribute to the overall knowledge on the biology of M. tuberculosis infection, whereas at the same time launch the necessary basis for the functional assessment of these so far overlooked genes.

Taking the pseudo out of pseudogenes

Pseudogenes are defined as fragments of once-functional genes that have been silenced by one or more nonsense, frameshift or missense mutations. Despite continuing increases in the speed of sequencing and annotating bacterial genomes, the identification and categorisation of pseudogenes remains problematic. Even when identified, pseudogenes are considered to be rare and tend to be ignored. On the contrary, pseudogenes are surprisingly prevalent and can persist for long evolutionary time periods, representing a record of once-functional genetic characteristics. Most importantly, pseudogenes provide an insight into prokaryotic evolutionary history as a record of phenotypic traits that have been lost. Focusing on the intracellular and symbiotic bacteria in which pseudogenes predominate, this review discusses the importance of identifying pseudogenes to fully understand the abilities of bacteria, and to understand prokaryotes within their evolutionary context.

Salmonella Typhi shdA: pseudogene or allelic variant?

ShdA from Salmonella Typhimurium (ShdASTm) is a large outer membrane protein that specifically recognizes and binds to fibronectin. ShdASTm is involved in the colonization of the cecum and the Peyer's patches of terminal ileum in mice. On the other hand, shdA gene from Salmonella Typhi (shdASTy) has been considered a pseudogene (i.e. a nonfunctional sequence of genomic DNA) due to the presence of deletions and mutations that gave rise to premature stop codons. In this work we show that, despite the deletions and mutations, shdASTy is fully functional. S. Typhi ΔshdA mutants presented an impaired adherence and invasion of HEp-2 pre-treated with TGF-β1, an inducer of fibronectin production. Moreover, shdA from S. Typhi and S. Typhimurium seem to be equivalent since shdASTm restored the adherence and invasion of S. Typhi ΔshdA mutant to wild type levels. In addition, anti-FLAG mAbs interfered with the adherence and invasion of the S. Typhi shdA-3xFLAG strain. Finally, shdASTy encodes a detectable protein when heterologously expressed in Escherichia coli DH5α. The data presented here show that shdASTy is not a pseudogene, but a different functional allele compared with shdASTm.

Design of tiling arrays and their application to bacterial transcriptome analysis

Whole-genome sequencing and annotation have clarified total gene number and structure in a variety of organisms. Microarrays have facilitated transcriptome analysis through the use of probes that target a large number of genes based on genomic information. However, microarrays are limited in that they can only examine known or predicted genes; non-annotated genes and noncoding regions cannot be accounted for.Recent advances in technology have led to the design of tiling arrays, which contain a vastly increased number of spotted probes, and at higher density. Tiling arrays cover the entire genome of a prokaryotic species in an unbiased fashion by designing a large number of probes. Upon hybridization of total RNA, all the transcribed regions of the genome, irrespective of gene annotation, can be detected. As opposed to next-generation sequencing, tiling arrays are cost-effective, easy to analyze, and have been used for experiments as diverse as transcriptome analysis, ChIP-chip, and DNA sequence variation detection. In this chapter, the methods for bacterial tiling array slide design, RNA sample preparation, hybridization, and data analysis are described.

Pseudogenes are not pseudo any more

Recent significant progress toward understanding the function of pseudogenes in protozoa (Trypanosoma brucei), metazoa (mouse) and plants, make it pertinent to provide a brief overview on what has been learned about this fascinating subject. We discuss the regulatory mechanisms of pseudogenes at the post-transcriptional level and advance new ideas toward understanding the evolution of these, sometimes called "garbage genes" or "junk DNA," seeking to stimulate the interest of scientists and additional research on the subject. We hope this point-of-view can be helpful to scientists working or seeking to work on these and related issues.

Chimpanzees used for medical research shed light on the pathoetiology of leprosy

Leprosy is a chronic infectious disorder caused by Mycobacterium leprae, which mainly affects skin and peripheral nerves. It is classified as either paucibacillary or multibacillary based upon clinical manifestations and slit-skin smear results. It is speculated that leprosy develops after a long latency period following M. leprae infection. However, the actual time of infection and the duration of latency have never been proven in human patients. To date, four cases of spontaneous leprosy have been reported in chimpanzees who were caught in West Africa in infancy and used for medical research in the USA and Japan. One of these chimpanzees was extensively studied in Japan, and single-nucleotide polymorphism analysis for the M. leprae genome was conducted. This analysis revealed that the chimpanzee was infected with M. leprae during infancy in West Africa and the pathognomonic signs of leprosy appeared after at least 30 years of incubation. Analysis of leprosy in chimpanzees can contribute not only to medical research but also to the understanding of the pathoetiology of leprosy.

Current status of leprosy: epidemiology, basic science and clinical perspectives

Leprosy has affected humans for millennia and remains an important health problem worldwide, as evidenced by nearly 250 000 new cases detected every year. It is a chronic infectious disorder, caused by Mycobacterium leprae, that primarily affects the skin and peripheral nerves. Recent advances in basic science have improved our knowledge of the disease. Variation in the cellular immune response is the basis of a range of clinical manifestations. The introduction of multidrug therapy has significantly contributed to a decrease in the prevalence of the disease. However, leprosy control activities, including monitoring and prevention programs, must be maintained.

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.

Pseudogenes: pseudo-functional or key regulators in health and disease?

Pseudogenes have long been labeled as "junk" DNA, failed copies of genes that arise during the evolution of genomes. However, recent results are challenging this moniker; indeed, some pseudogenes appear to harbor the potential to regulate their protein-coding cousins. Far from being silent relics, many pseudogenes are transcribed into RNA, some exhibiting a tissue-specific pattern of activation. Pseudogene transcripts can be processed into short interfering RNAs that regulate coding genes through the RNAi pathway. In another remarkable discovery, it has been shown that pseudogenes are capable of regulating tumor suppressors and oncogenes by acting as microRNA decoys. The finding that pseudogenes are often deregulated during cancer progression warrants further investigation into the true extent of pseudogene function. In this review, we describe the ways in which pseudogenes exert their effect on coding genes and explore the role of pseudogenes in the increasingly complex web of noncoding RNA that contributes to normal cellular regulation.

Mycobacterium leprae-host-cell interactions and genetic determinants in leprosy: an overview

Leprosy, also known as Hansen's disease, is a chronic infectious disease caused by Mycobacterium leprae in which susceptibility to the mycobacteria and its clinical manifestations are attributed to the host immune response. Even though leprosy prevalence has decreased dramatically, the high number of new cases indicates active transmission. Owing to its singular features, M. leprae infection is an attractive model for investigating the regulation of human immune responses to pathogen-induced disease. Leprosy is one of the most common causes of nontraumatic peripheral neuropathy worldwide. The proportion of patients with disabilities is affected by the type of leprosy and delay in diagnosis. This article briefly reviews the clinical features as well as the immunopathological mechanisms related to the establishment of the different polar forms of leprosy, the mechanisms related to M. leprae-host cell interactions and prophylaxis and diagnosis of this complex disease. Host genetic factors are summarized and the impact of the development of interventions that prevent, reverse or limit leprosy-related nerve impairments are discussed.

Detection of Mycobacterium leprae DNA from archaeological skeletal remains in Japan using whole genome amplification and polymerase chain reaction

Background

Identification of pathogen DNA from archaeological human remains is a powerful tool in demonstrating that the infectious disease existed in the past. However, it is very difficult to detect trace amounts of DNA remnants attached to the human skeleton, especially from those buried in a humid atmosphere with a relatively high environmental temperature such as in Asia.

Methodology/Principal Findings

Here we demonstrate Mycobacterium leprae DNA from archaeological skeletal remains in Japan by polymerase chain reaction, DNA sequencing and single nucleotide polymorphism (SNP) analysis. In addition, we have established a highly sensitive method of detecting DNA using a combination of whole genome amplification and polymerase chain reaction, or WGA-PCR, which provides superior sensitivity and specificity in detecting DNA from trace amounts of skeletal materials.

Conclusion/Significance

We have detected M. leprae DNA in archaeological skeletal remains for the first time in the Far East. Its SNP genotype corresponded to type 1; the first detected case worldwide of ancient M. leprae DNA. We also developed a highly sensitive method to detect ancient DNA by utilizing whole genome amplification.

Analysis of Mycobacterium leprae gene expression using DNA microarray

Mycobacterium leprae, the causative agent of leprosy, does not grow under in vitro condition, making molecular analysis of this bacterium difficult. For this reason, bacteriological information regarding M. leprae gene function is limited compared with other mycobacterium species. In this study, we performed DNA microarray analysis to clarify the RNA expression profile of the Thai53 strain of M. leprae grown in footpads of hypertensive nude rats (SHR/NCrj-rnu). Of 1605 M. leprae genes, 315 showed signal intensity twofold higher than the median. These genes include Acyl-CoA metabolic enzymes and drug metabolic enzymes, which might be related to the virulence of M. leprae. In addition, consecutive RNA expression profile and in silico analyses enabled identification of possible operons within the M. leprae genome. The present results will shed light on M. leprae gene function and further our understanding of the pathogenesis of leprosy.

The Frankia alni symbiotic transcriptome

The actinobacteria Frankia spp. are able to induce the formation of nodules on the roots of a large spectrum of actinorhizal plants, where they convert dinitrogen to ammonia in exchange for plant photosynthates. In the present study, transcriptional analyses were performed on nitrogen-replete free-living Frankia alni cells and on Alnus glutinosa nodule bacteria, using whole-genome microarrays. Distribution of nodule-induced genes on the genome was found to be mostly over regions with high synteny between three Frankia spp. genomes, while nodule-repressed genes, which were mostly hypothetical and not conserved, were spread around the genome. Genes known to be related to nitrogen fixation were highly induced, nif (nitrogenase), hup2 (hydrogenase uptake), suf (sulfur-iron cluster), and shc (hopanoids synthesis). The expression of genes involved in ammonium assimilation and transport was strongly modified, suggesting that bacteria ammonium assimilation was limited. Genes involved in particular in transcriptional regulation, signaling processes, protein drug export, protein secretion, lipopolysaccharide, and peptidoglycan biosynthesis that may play a role in symbiosis were also identified. We also showed that this Frankia symbiotic transcriptome was highly similar among phylogenetically distant plant families Betulaceae and Myricaceae. Finally, comparison with rhizobia transcriptome suggested that F. alni is metabolically more active in symbiosis than rhizobia.

Validating divergent ORF annotation of the Mycobacterium leprae genome through a full translation data set and peptide identification by tandem mass spectrometry

Mycobacterium leprae has undergone extensive degenerative evolution, with a large number of pseudogenes. It is also the organism with the greatest divergence between gene annotations from independent institutes. Therefore, M. leprae is a good model to verify the currently predicted coding sequence regions between different annotations, to identify new ones and to investigate the expression of pseudogenes. We submitted a total extract of the bacteria isolated from Armadillo to Gel-LC-MS/MS using a linear quadrupole ion trap-Orbitrap mass spectrometer. Spectra were analyzed using the Leproma (1614 genes and 1133 pseudogenes) and TIGR (5446 genes) databases and a database containing the full genome translation. We identified a total of 1046 proteins, including five proteins encoded by previously predicted pseudogenes, which upon closer inspection appeared to be proper genes. Only 11 of the additional annotations by TIGR were verified. We also identified six tryptic peptides from five proteins from regions not considered to be coding sequences, in addition to peptides from two unannotated gene candidates that overlap with other genes. Our data show that the Leproma annotation of M. leprae is quite accurate, and there were no peptide observations corresponding to true pseudogenes, except for a new gene candidate, overlapping with an essential enolase on the complementary strand.

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.

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.

[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.

[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.

Expression of adipose differentiation-related protein (ADRP) and perilipin in macrophages infected with Mycobacterium leprae

Mycobacterium leprae survives and replicates within a lipid droplet stored in the enlarged phagosome of histiocytes, a typical feature of lepromatous leprosy that is thought to be an important nutrient source for the bacillus. However, the underlying mechanisms by which lipids accumulate within phagosomes remain unclear. Recently, it was revealed that the lipid droplet-associated proteins, including ADRP and perilipin, play essential roles in lipid accumulation in adipocytes or macrophages. Therefore, we attempted to examine the role of these proteins in leprosy pathogenesis. ADRP and perilipin localized to the phagosomal membrane, which contains M. leprae in skin biopsy specimens of lepromatous leprosy. ADRP expression was transiently increased after phagocytosis in THP-1 cells. However, high levels of ADRP expression persisted only when live M. leprae, but not dead bacilli or latex beads, was added. Furthermore, although peptidoglycan, a Toll-like receptor 2 ligand, suppressed the expression levels of ADRP and perilipin, M. leprae infection inhibited this suppression. These results suggest that live M. leprae has the ability to actively induce and support ADRP/perilipin expression to facilitate the accumulation of lipids within the phagosome and to further maintain a suitable environment for the intracellular survival within the macrophage.