MycoProtease-DB: Useful resource for Mycobacterium tuberculosis complex and nontuberculous mycobacterial proteases

MycoProtease-DB is an online MS SQL and CGI-PERL driven relational database that domiciles protease information of Mycobacterium tuberculosis (MTB) complex and Nontuberculous Mycobacteria (NTM), whose complete genome sequence is available. Our effort is to provide comprehensive information on proteases of 5 strains of Mycobacterium tuberculosis (H₃₇Rv, H₃₇Ra, CDC1551, F11 and KZN 1435), 3 strains of Mycobacterium bovis (AF2122/97, BCG Pasteur 1173P2 and BCG Tokyo 172) and 4 strains of NTM (Mycobacterium avium 104, Mycobacterium smegmatis MC2 155, Mycobacterium avium paratuberculosis K-10 and Nocardia farcinica IFM 10152) at gene, protein and structural level. MycoProtease-DB currently hosts 1324 proteases, which include 906 proteases from MTB complex with 237distinct proteases & 418 from NTM with 404 distinct proteases. Flexible database design and easy expandability & retrieval of information are the main features of MycoProtease-DB. All the data were validated with various online resources and published literatures for reliable serving as comprehensive resources of various Mycobacterial proteases.

Availability

The Database is publicly available at http://www.bicjbtdrc-mgims.in/MycoProtease-DB/     

AKE - the Accelerated k-mer Exploration web-tool for rapid taxonomic classification and visualization

Background

With the advent of low cost, fast sequencing technologies metagenomic analyses are made possible. The large data volumes gathered by these techniques and the unpredictable diversity captured in them are still, however, a challenge for computational biology.

Results

In this paper we address the problem of rapid taxonomic assignment with small and adaptive data models (< 5 MB) and present the accelerated k-mer explorer (AKE). Acceleration in AKE’s taxonomic assignments is achieved by a special machine learning architecture, which is well suited to model data collections that are intrinsically hierarchical. We report classification accuracy reasonably well for ranks down to order, observed on a study on real world data (Acid Mine Drainage, Cow Rumen).

Conclusion

We show that the execution time of this approach is orders of magnitude shorter than competitive approaches and that accuracy is comparable. The tool is presented to the public as a web application (url: https://ani.cebitec.uni-bielefeld.de/ake/, username: bmc, password: bmcbioinfo).

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0384-0) contains supplementary material, which is available to authorized users.     

Differential Immune Responses and Protective Effects in Avirulent Mycobacterial Strains Vaccinated BALB/c Mice

Screening live mycobacterial vaccine candidates is the important strategy to develop new vaccines against adult tuberculosis (TB). In this study, the immunogenicity and protective efficacy of several avirulent mycobacterial strains including Mycobacterium smegmatis, M. vaccae, M. terrae, M. phlei, M. trivial, and M. tuberculosis H37Ra were compared with M. bovis BCG in BALB/c mice. Our results demonstrated that differential immune responses were induced in different mycobacterial species vaccinated mice. As BCG-vaccinated mice did, M. terrae immunization resulted in Th1-type responses in the lung, as well as splenocytes secreting IFN-γ against a highly conserved mycobacterial antigen Ag85A. M. smegmatis also induced the same splenocytes secreting IFN-γ as BCG and M. terrae did. In addition, M. terrae and M. smegmatis-immunized mice predominantly increased expression of IL-10 and TGF-β in the lung. Most importantly, mice vaccinated with H37Ra and M. vaccae could provide the same protection in the lung against virulent M. tuberculosis challenge as BCG. The result may have important implications in developing adult TB vaccine.     

PhyTB: Phylogenetic tree visualisation and sample positioning for M. tuberculosis

Background

Phylogenetic-based classification of M. tuberculosis and other bacterial genomes is a core analysis for studying evolutionary hypotheses, disease outbreaks and transmission events. Whole genome sequencing is providing new insights into the genomic variation underlying intra- and inter-strain diversity, thereby assisting with the classification and molecular barcoding of the bacteria. One roadblock to strain investigation is the lack of user-interactive solutions to interrogate and visualise variation within a phylogenetic tree setting.

Results

We have developed a web-based tool called PhyTB (http://pathogenseq.lshtm.ac.uk/phytblive/index.php) to assist phylogenetic tree visualisation and identification of M. tuberculosis clade-informative polymorphism. Variant Call Format files can be uploaded to determine a sample position within the tree. A map view summarises the geographical distribution of alleles and strain-types. The utility of the PhyTB is demonstrated on sequence data from 1,601 M. tuberculosis isolates.

Conclusion

PhyTB contextualises M. tuberculosis genomic variation within epidemiological, geographical and phylogenic settings. Further tool utility is possible by incorporating large variants and phenotypic data (e.g. drug-resistance profiles), and an assessment of genotype-phenotype associations. Source code is available to develop similar websites for other organisms (http://sourceforge.net/projects/phylotrack).     

A Study of Unusual Pacemaker Infection by Mycobacterium Tuberculosis in Indian Patients

Background

The expanding clinical indications of cardiac rhythm management have led to an increased use of pacemaker implantation which is associated with increased incidence of pacemaker infections. Staphylococcus aureus and epidermidis account for the vast majority of pacemaker infections. Pacemaker infection due to Mycobacterium tuberculosis (M. tuberculosis) is very rare, only few cases having been reported till date.

Methods

We describe here a study of three patients of pacemaker pocket infection with M. tuberculosis.

Conclusion

The possibility of mycobacterial pacemaker infection should always be kept in mind in patients with delayed pacemaker infection.     

Evaluation of amplified rDNA restriction analysis (ARDRA) for the identification of cultured mycobacteria in a diagnostic laboratory

Background

The development of DNA amplification for the direct detection of M. tuberculosis from clinical samples has been a major goal of clinical microbiology during the last ten years. However, the limited sensitivity of most DNA amplification techniques restricts their use to smear positive samples. On the other hand, the development of automated liquid culture has increased the speed and sensitivity of cultivation of mycobacteria. We have opted to combine automated culture with rapid genotypic identification (ARDRA: amplified rDNA restriction analysis) for the detection resp. identification of all mycobacterial species at once, instead of attempting direct PCR based detection from clinical samples of M. tuberculosis only.

Results

During 1998–2000 a total of approx. 3500 clinical samples was screened for the presence of M. tuberculosis. Of the 151 culture positive samples, 61 were M. tuberculosis culture positive. Of the 30 smear positive samples, 26 were M. tuberculosis positive. All but three of these 151 mycobacterial isolates could be identified with ARDRA within on average 36 hours. The three isolates that could not be identified belonged to rare species not yet included in our ARDRA fingerprint library or were isolates with an aberrant pattern.

Conclusions

In our hands, automated culture in combination with ARDRA provides with accurate, practically applicable, wide range identification of mycobacterial species. The existing identification library covers most species, and can be easily updated when new species are studied or described. The drawback is that ARDRA is culture-dependent, since automated culture of M. tuberculosis takes on average 16.7 days (range 6 to 29 days). However, culture is needed after all to assess the antibiotic susceptibility of the strains.

     

Exhibition of persistent and drug-tolerant L-form habit of <Emphasis Type="Italic">Mycobacterium tuberculosis</Emphasis> during infection in rats

A model for studying mycobacterial L-form formation in vivo was established to demonstrate the ability of M. tuberculosis to behave as a drug-tolerant L-form persister. Rats were infected by intranasal (i.n.) and intraperitoneal (i.p.) routes with 1×108 cells/ml of M. tuberculosis. At weekly intervals during a period of five weeks, samples from lung, spleen, liver, kidney, mesenterial and inguinal lymph nodes, broncho-alveolar and peritoneal lavage liquid were plated simultaneously on Löwenstein-Jensen (LJ) medium or inoculated into specially supplemented for L-forms Dubos broth (drug-free and drug-containing variants). The use of liquid media enabled isolation of mycobacterial L-form cultures during the whole period of experiment including the last two weeks, when tubercle bacilli were not isolated on LJ medium. An unique feature of mycobacterial L-forms was their ability to grow faster than the classical tubercle bacilli. Isolation and growth of L-form cultures in primary drug-containing media demonstrated their drug-tolerant properties. Electron microscopy of liquid media isolates showed that they consisted of morphologically heterogenous populations of membrane-bound and of variable sized L-bodies that completely lack cell walls. The identity of the isolated non-acid fast and morphologically modified L-forms as M. tuberculosis was verified by specific spoligotyping test. The results contribute to special aspects concerning the importance of mycobacterial L-form phenomenon for persistence and latency in tuberculosis, phenotypic drug tolerance, as well as for diagnosis of difficult to identify morphologically changed tubercle bacilli which are often mistaken for contaminants.     

Genome sequence and description of Alistipes senegalensis sp. nov.

Alistipes senegalensis strain JC50^T is the type strain of A. senegalensis sp. nov., a new species within the Alistipes genus. This strain, whose genome is described here, was isolated from the fecal flora of an asymptomatic patient. A. senegalensis is an anaerobic Gram-negative rod-shaped bacterium. Here we describe the features of this organism, together with the complete genome sequence and annotation. The 4,017,609 bp long genome (1 chromosome, but no plasmid) contains 3,113 protein-coding and 50 RNA genes, including 5 rRNA genes.     

Genome sequence and description of Aeromicrobium massiliense sp. nov.

Aeromicrobium massiliense strain JC14^Tsp. nov. is the type strain of Aeromicrobium massiliense sp. nov., a new species within the genus Aeromicrobium. This strain, whose genome is described here, was isolated from the fecal microbiota of an asymptomatic patient. Aeromicrobium massiliense is an aerobic rod-shaped gram-positive bacterium. Here we describe the features of this organism, together with the complete genome sequence and annotation. The 3,322,119 bp long genome contains 3,296 protein-coding and 51 RNA genes.     

Complete genome sequence of the sulfate-reducing firmicute Desulfotomaculum ruminis type strain (DLT)

Desulfotomaculum ruminis Campbell and Postgate 1965 is a member of the large genus Desulfotomaculum which contains 30 species and is contained in the family Peptococcaceae. This species is of interest because it represents one of the few sulfate-reducing bacteria that have been isolated from the rumen. Here we describe the features of D. ruminis together with the complete genome sequence and annotation. The 3,969,014 bp long chromosome with a total of 3,901 protein-coding and 85 RNA genes is the second completed genome sequence of a type strain of the genus Desulfotomaculum to be published, and was sequenced as part of the DOE Joint Genome Institute Community Sequencing Program 2009.Keywords : anaerobic, motile, sporulating, mesophilic, sulfate-reducer, hydrogen sulfide, incomplete oxidizer, mixotrophic, CSP 2009, Peptococcaceae, Clostridiales     

Comparative Analysis of Pdf-Mediated Circadian Behaviors Between Drosophila melanogaster and D. virilis

PAR proteins (partitioning defective) are major regulators of cell polarity and asymmetric cell division. One of the par genes, par-1, encodes a Ser/Thr kinase that is conserved from yeast to mammals. In Caenorhabditis elegans, par-1 governs asymmetric cell division by ensuring the polar distribution of cell fate determinants. However the precise mechanisms by which PAR-1 regulates asymmetric cell division in C. elegans remain to be elucidated. We performed a genomewide RNAi screen and identified six genes that specifically suppress the embryonic lethal phenotype associated with mutations in par-1. One of these suppressors is mpk-1, the C. elegans homolog of the conserved mitogen activated protein (MAP) kinase ERK. Loss of function of mpk-1 restored embryonic viability, asynchronous cell divisions, the asymmetric distribution of cell fate specification markers, and the distribution of PAR-1 protein in par-1 mutant embryos, indicating that this genetic interaction is functionally relevant for embryonic development. Furthermore, disrupting the function of other components of the MAPK signaling pathway resulted in suppression of par-1 embryonic lethality. Our data therefore indicates that MAP kinase signaling antagonizes PAR-1 signaling during early C. elegans embryonic polarization.ASYMMETRIC cell division, a process in which a mother cell divides in two different daughter cells, is a fundamental mechanism to achieve cell diversity during development. We use the early embryo of Caenorhabditis elegans as a model system to study asymmetric cell division. The C. elegans one-cell embryo divides asymmetrically along its anteroposterior axis, generating two cells of different sizes and fates: the larger anterior daughter cell will generate somatic tissues while the smaller posterior daughter cell will generate the germline (Sulston et al. 1983).A group of proteins called PAR proteins (partitioning defective) is required for asymmetric cell division in C. elegans (Kemphues et al. 1988). Depletion of any of the seven par genes (par-1 to -6 and pkc-3) leads to defects in asymmetric cell division and embryonic lethality (Kemphues et al. 1988; Kirby et al. 1990; Tabuse et al. 1998; Hung and Kemphues 1999; Hao et al. 2006). PAR-3 and PAR-6 are conserved proteins that contain PDZ-domains and form a complex with PKC-3 (Etemad-Moghadam et al. 1995; Izumi et al. 1998; Tabuse et al. 1998; Hung and Kemphues 1999). This complex becomes restricted to the anterior cortex of the embryo in response to spatially defined actomyosin contractions occurring in the embryo upon fertilization (Goldstein and Hird 1996; Munro et al. 2004). The posterior cortex of the embryo that becomes devoid of the anterior PAR proteins is occupied by the RING protein PAR-2 and the Ser/Thr kinase PAR-1 (Guo and Kemphues 1995; Boyd et al. 1996; Cuenca et al. 2003). Once polarized, the anterior and posterior PAR proteins mutually exclude each other from their respective cortices (Etemad-Moghadam et al. 1995; Boyd et al. 1996; Cuenca et al. 2003; Hao et al. 2006). Loss of function of the gene par-1, as opposed to loss of most other par genes, results in embryos that exhibit only subtle effects on the polarized cortical domains occupied by the other PAR proteins (Cuenca et al. 2003). However defects in this gene are associated with a more symmetric division in size, an aberrant distribution of cell fate specification markers, altered cell fates of the daughter cells of the embryo, and ultimately embryonic lethality (Kemphues et al. 1988; Guo and Kemphues 1995).PAR-1 controls asymmetric cell division and cell fate specification by regulating the localization of the two highly similar CCCH-type zinc-finger proteins MEX-5 and MEX-6 (referred to as MEX-5/6). MEX-5 and MEX-6 are 70% identical in their amino acid sequence and fulfill partially redundant functions in the embryo (Schubert et al. 2000). In wild-type animals, endogenous MEX-5 and GFP fusions of MEX-6 localize primarily to the anterior of the embryo while both proteins are evenly distributed in par-1 mutant embryos (Schubert et al. 2000; Cuenca et al. 2003). This suggests that in wild-type animals, PAR-1 acts in part by restricting MEX-5 and MEX-6 to the anterior of the embryo. The precise mechanism of this regulation is not known, but an elegant study performed for MEX-5 indicates that differential protein mobility in the anterior and posterior cytoplasm of the one-cell embryo contributes to this asymmetry (Tenlen et al. 2008). While increased mobility in the posterior of the one-cell embryo correlates with a par-1- and par-4-dependent phosphorylation on MEX-5, the kinase directly phosphorylating MEX-5 remains to be identified (Tenlen et al. 2008).Some of the phenotypes associated with loss of par-1 function are dependent on the function of mex-5 and mex-6. First, loss of function of par-1 leads to a decreased stability and aberrant localization of the posterior cell fate specification marker PIE-1, a protein that is usually inherited by the posterior daughter cell in wild-type animals and ensures the correct specification of the germline (Mello et al. 1996; Seydoux et al. 1996). This decreased stability is dependent on mex-5/6 function as PIE-1 levels are restored, albeit with symmetrical distribution, in mex-6(RNAi); mex-5(RNAi); par-1(b274) embryos (Schubert et al. 2000; Cuenca et al. 2003; Derenzo et al. 2003). Second, embryos lacking par-1 function exhibit decreased amounts of P granules in the one-cell embryo, while these markers are present in mex-6(pk440); mex-5(zu199); par-1(RNAi) embryos of comparable age (Cheeks et al. 2004). Third, in par-1(RNAi) one-cell embryos the posterior cortical domain occupied by the polarity protein PAR-2 is extended anteriorly, when compared to wild-type embryos (Cuenca et al. 2003). This anterior extension is rescued in embryos deficient for both par-1 and mex-5/6 (Cuenca et al. 2003). Taken together, these results indicate that par-1 acts in the embryo—at least in part—by regulating the localization and/or activity of the proteins MEX-5 and MEX-6. However, it remains unclear whether other proteins can modulate PAR-1 function to affect MEX-5/6 activity.To gain insight into the mechanisms of par-1 function in the embryo, we sought to identify genes that act together with par-1 during embryonic development. We performed an RNAi-based screen for genetic interactors of the temperature-sensitive allele par-1(zu310), using the embryonic lethal phenotype of this mutant as a readout. This method has proven successful in previous screens to identify genes involved in early embryonic processes (Labbé et al. 2006; O''Rourke et al. 2007). We were able to identify six genes that, upon disruption of their function, suppress the embryonic lethal phenotype of par-1 mutant embryos. One of these genes is mpk-1, the C. elegans homolog of the highly conserved MAP kinase ERK. Closer analysis subsequently showed that reduction of function of mpk-1 not only increases viability of par-1 mutant embryos, but also reverts several polarity phenotypes associated with loss of function of par-1. Our data indicate that mpk-1 antagonizes par-1 activity to regulate polarization and asymmetric cell divisions in the early embryo.     

ALLocator: An Interactive Web Platform for the Analysis of Metabolomic LC-ESI-MS Datasets,Enabling Semi-Automated,User-Revised Compound Annotation and Mass Isotopomer Ratio Analysis

Adduct formation, fragmentation events and matrix effects impose special challenges to the identification and quantitation of metabolites in LC-ESI-MS datasets. An important step in compound identification is the deconvolution of mass signals. During this processing step, peaks representing adducts, fragments, and isotopologues of the same analyte are allocated to a distinct group, in order to separate peaks from coeluting compounds. From these peak groups, neutral masses and pseudo spectra are derived and used for metabolite identification via mass decomposition and database matching. Quantitation of metabolites is hampered by matrix effects and nonlinear responses in LC-ESI-MS measurements. A common approach to correct for these effects is the addition of a U-¹³C-labeled internal standard and the calculation of mass isotopomer ratios for each metabolite. Here we present a new web-platform for the analysis of LC-ESI-MS experiments. ALLocator covers the workflow from raw data processing to metabolite identification and mass isotopomer ratio analysis. The integrated processing pipeline for spectra deconvolution “ALLocatorSD” generates pseudo spectra and automatically identifies peaks emerging from the U-¹³C-labeled internal standard. Information from the latter improves mass decomposition and annotation of neutral losses. ALLocator provides an interactive and dynamic interface to explore and enhance the results in depth. Pseudo spectra of identified metabolites can be stored in user- and method-specific reference lists that can be applied on succeeding datasets. The potential of the software is exemplified in an experiment, in which abundance fold-changes of metabolites of the l-arginine biosynthesis in C. glutamicum type strain ATCC 13032 and l-arginine producing strain ATCC 21831 are compared. Furthermore, the capability for detection and annotation of uncommon large neutral losses is shown by the identification of (γ-)glutamyl dipeptides in the same strains. ALLocator is available online at: https://allocator.cebitec.uni-bielefeld.de. A login is required, but freely available.     

Complete genome sequence of the motile actinomycete Actinoplanes missouriensis 431T (= NBRC 102363T)

Actinoplanes missouriensis Couch 1963 is a well-characterized member of the genus Actinoplanes, which is of morphological interest because its members typically produce sporangia containing motile spores. The sporangiospores are motile by means of flagella and exhibit chemotactic properties. It is of further interest that members of Actinoplanes are prolific sources of novel antibiotics, enzymes, and other bioactive compounds. Here, we describe the features of A. missouriensis 431^T, together with the complete genome sequence and annotation. The 8,773,466 bp genome contains 8,125 protein-coding and 79 RNA genes.     

A Novel Role of the N Terminus of Budding Yeast Histone H3 Variant Cse4 in Ubiquitin-Mediated Proteolysis

Understanding the molecular basis of common traits is a primary challenge of modern genetics. One model holds that rare mutations in many genetic backgrounds may often phenocopy one another, together explaining the prevalence of the resulting trait in the population. For the vast majority of phenotypes, the role of rare variants and the evolutionary forces that underlie them are unknown. In this work, we use a population of Saccharomyces paradoxus yeast as a model system for the study of common trait variation. We observed an unusual, flocculation and invasive-growth phenotype in one-third of S. paradoxus strains, which were otherwise unrelated. In crosses with each strain in turn, these morphologies segregated as a recessive Mendelian phenotype, mapping either to IRA1 or to IRA2, yeast homologs of the hypermutable human neurofibromatosis gene NF1. The causal IRA1 and IRA2 haplotypes were of distinct evolutionary origin and, in addition to their morphological effects, associated with hundreds of stress-resistance and growth traits, both beneficial and disadvantageous, across S. paradoxus. Single-gene molecular genetic analyses confirmed variant IRA1 and IRA2 haplotypes as causal for these growth characteristics, many of which were independent of morphology. Our data make clear that common growth and morphology traits in yeast result from a suite of variants in master regulators, which function as a mutation-driven switch between phenotypic states.