Disaccharide analogs as probes for glycosyltransferases in Mycobacterium tuberculosis

Glycosyltransferases (GTs) play a crucial role in mycobacterial cell wall biosynthesis and are necessary for the survival of mycobacteria. Hence, these enzymes are potential new drug targets for the treatment of tuberculosis (TB), especially multiple drug-resistant TB (MDR-TB). Herein, we report the efficient syntheses of Araf(alpha 1-->5)Araf, Galf(beta 1-->5)Galf, and Galf(beta 1-->6)Galf disaccharides possessing a 5-N,N-dimethylaminonaphthalene-1-sulfonamidoethyl (dansyl) unit that were prepared as fluorescent disaccharide acceptors for arabinosyl- and galactosyl-transferases, respectively. Such analogs may offer advantages relative to radiolabeled acceptors or donors for studying the enzymes and for assay development and compound screening. Additionally, analogs possessing a 5-azidonaphthalene-1-sulfonamidoethyl unit were prepared as photoaffinity probes for their potential utility in studying active site labeling of the GTs (arabinosyl and galactosyl) in Mycobacterium tuberculosis (MTB). Beyond their preparation, initial biological testing and kinetic analysis of these disaccharides as acceptors toward glycosyltransferases are also presented.     

Recognition of fold and sugar linkage for glycosyltransferases by multivariate sequence analysis

Glycosyltransferases (GTs) are among the largest groups of enzymes found and are usually classified on the basis of sequence comparisons into many families of varying similarity (CAZy systematics). Only two different Rossman-like folds have been detected (GT-A and GT-B) within the small number of established crystal structures. A third uncharacterized fold has been indicated with transmembrane organization (GT-C). We here use a method based on multivariate data analyses (MVDAs) of property patterns in amino acid sequences and can with high accuracy recognize the correct fold in a large data set of GTs. Likewise, a retaining or inverting enzymatic mechanism for attachment of the donor sugar could be properly revealed in the GT-A and GT-B fold group sequences by such analyses. Sequence alignments could be correlated to important variables in MVDA, and the separating amino acid positions could be mapped over the active sites. These seem to be localized to similar positions in space for the alpha/beta/alpha binding motifs in the GT-B fold group structures. Analogous, active-site sequence positions were found for the GT-A fold group. Multivariate property patterns could also easily group most GTs annotated in the genomes of Escherichia coli and Synechocystis to proper fold or organization group, according to benchmarking comparisons at the MetaServer. We conclude that the sequence property patterns revealed by the multivariate analyses seem more conserved than amino acid types for these GT groups, and these patterns are also conserved in the structures. Such patterns may also potentially define substrate preferences.     

Whole genome searching with shotgun proteomic data: applications for genome annotation

High-throughput genome sequencing continues to accelerate the rate at which complete genomes are available for biological research. Many of these new genome sequences have little or no genome annotation currently available and hence rely upon computational predictions of protein coding genes. Evidence of translation from proteomic techniques could facilitate experimental validation of protein coding genes, but the techniques for whole genome searching with MS/MS data have not been adequately developed to date. Here we describe GENQUEST, a novel method using peptide isoelectric focusing and accurate mass to greatly reduce the peptide search space, making fast, accurate, and sensitive whole human genome searching possible on common desktop computers. In an initial experiment, almost all exonic peptides identified in a protein database search were identified when searching genomic sequence. Many peptides identified exclusively in the genome searches were incorrectly identified or could not be experimentally validated, highlighting the importance of orthogonal validation. Experimentally validated peptides exclusive to the genomic searches can be used to reannotate protein coding genes. GENQUEST represents an experimental tool that can be used by the proteomics community at large for validating computational approaches to genome annotation.     

Optimizing multiplex SNP-based data analysis for genotyping of Mycobacterium tuberculosis isolates

Background

Multiplex ligation-dependent probe amplification (MLPA) is a powerful tool to identify genomic polymorphisms. We have previously developed a single nucleotide polymorphism (SNP) and large sequence polymorphisms (LSP)-based MLPA assay using a read out on a liquid bead array to screen for 47 genetic markers in the Mycobacterium tuberculosis genome. In our assay we obtain information regarding the Mycobacterium tuberculosis lineage and drug resistance simultaneously. Previously we called the presence or absence of a genotypic marker based on a threshold signal level. Here we present a more elaborate data analysis method to standardize and streamline the interpretation of data generated by MLPA. The new data analysis method also identifies intermediate signals in addition to classification of signals as positive and negative. Intermediate calls can be informative with respect to identifying the simultaneous presence of sensitive and resistant alleles or infection with multiple different Mycobacterium tuberculosis strains.

Results

To validate our analysis method 100 DNA isolates of Mycobacterium tuberculosis extracted from cultured patient material collected at the National TB Reference Laboratory of the National Center for Tuberculosis and Lung Diseases in Tbilisi, Republic of Georgia were tested by MLPA. The data generated were interpreted blindly and then compared to results obtained by reference methods. MLPA profiles containing intermediate calls are flagged for expert review whereas the majority of profiles, not containing intermediate calls, were called automatically. No intermediate signals were identified in 74/100 isolates and in the remaining 26 isolates at least one genetic marker produced an intermediate signal.

Conclusion

Based on excellent agreement with the reference methods we conclude that the new data analysis method performed well. The streamlined data processing and standardized data interpretation allows the comparison of the Mycobacterium tuberculosis MLPA results between different experiments. All together this will facilitate the implementation of the MLPA assay in different settings.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-572) contains supplementary material, which is available to authorized users.     

Variable human minisatellite-like regions in the Mycobacterium tuberculosis genome

Mycobacterial interspersed repetitive units (MIRUs) are 40-100 bp DNA elements often found as tandem repeats and dispersed in intergenic regions of the Mycobacterium tuberculosis complex genomes. The M. tuberculosis H37Rv chromosome contains 41 MIRU loci. After polymerase chain reaction (PCR) and sequence analyses of these loci in 31 M. tuberculosis complex strains, 12 of them were found to display variations in tandem repeat copy numbers and, in most cases, sequence variations between repeat units as well. These features are reminiscent of those of certain human variable minisatellites. Of the 12 variable loci, only one was found to vary among genealogically distant BCG substrains, suggesting that these interspersed bacterial minisatellite-like structures evolve slowly in mycobacterial populations.     

Dataset of potential targets for Mycobacterium tuberculosis H37Rv through comparative genome analysis

Mycobacterium tuberculosis is the causative agent of the disease, tuberculosis and H37Rv is the most studied clinical strain. We use comparative genome analysis of Mycobacterium tuberculosis H37Rv and human for the identification of potential targets dataset. We used DEG (Database of Essential Genes) to identify essential genes in the H37Rv strain. The analysis shows that 628 of the 3989 genes in Mycobacterium tuberculosis H37Rv were found to be essential of which 324 genes lack similarity to the human genome. Subsequently hypothetical proteins were removed through manual curation. This further resulted in a dataset of 135 proteins with essential function and no homology to human.     

Multivariate exploratory tools for microarray data analysis

The ultimate success of microarray technology in basic and applied biological sciences depends critically on the development of statistical methods for gene expression data analysis. The most widely used tests for differential expression of genes are essentially univariate. Such tests disregard the multidimensional structure of microarray data. Multivariate methods are needed to utilize the information hidden in gene interactions and hence to provide more powerful and biologically meaningful methods for finding subsets of differentially expressed genes. The objective of this paper is to develop methods of multidimensional search for biologically significant genes, considering expression signals as mutually dependent random variables. To attain these ends, we consider the utility of a pertinent distance between random vectors and its empirical counterpart constructed from gene expression data. The distance furnishes exploratory procedures aimed at finding a target subset of differentially expressed genes. To determine the size of the target subset, we resort to successive elimination of smaller subsets resulting from each step of a random search algorithm based on maximization of the proposed distance. Different stopping rules associated with this procedure are evaluated. The usefulness of the proposed approach is illustrated with an application to the analysis of two sets of gene expression data.     

Scanning the genome of Mycobacterium tuberculosis to identify potential lectins

Lectins are carbohydrate binding proteins with important roles in many biological processes such as adhesion. Here we have identified 11 potential lectins from Mycobacterium tuberculosis H37Rv genome, using a comprehensive bioinformatics analysis, which will provide helpful clues in molecular mapping of pathogenesis and perhaps also serve as potential drug targets.     

The glycosyltransferases of Mycobacterium tuberculosis - roles in the synthesis of arabinogalactan, lipoarabinomannan, and other glycoconjugates

Several human pathogens are to be found within the bacterial genus Mycobacterium, notably Mycobacterium tuberculosis, the causative agent of tuberculosis, one of the most threatening of human infectious diseases, with an annual lethality of about two million people. The characteristic mycobacterial cell envelope is the dominant feature of the biology of M. tuberculosis and other mycobacterial pathogens, based on sugars and lipids of exceptional structure. The cell wall consists of a peptidoglycan-arabinogalactan-mycolic acid complex beyond the plasma membrane. Free-standing lipids, lipoglycans, and proteins intercalate within this complex, complement the mycolic acid monolayer and may also appear in a capsular-like arrangement. The consequences of these structural oddities are an extremely robust and impermeable cell envelope. This review reflects on these entities from the perspective of their synthesis, particularly the structural and functional aspects of the glycosyltransferases (GTs) of M. tuberculosis, the dominating group of enzymes responsible for the terminal stages of their biosynthesis. Besides the many nucleotide-sugar dependent GTs with orthologs in prokaryotes and eukaryotes, M. tuberculosis and related species of the order Actinomycetales, in light of the highly lipophilic environment prevailing within the cell envelope, carry a significant number of GTs of the GT-C class dependent on polyprenyl-phosphate-linked sugars. These are of special emphasis in this review.     

Combining data in phylogenetic analysis

Systematists have access to multiple sources of character information in phylogenetic analysis. For example, it is not unusual to have nucleotide sequences from several different genes, or to have molecular and morphological data. How should diverse data be analyzed in phylogenetic analysis? Several methods have been proposed for the treatment of partitioned data: the total evidence, separate analysis, and conditional combination approaches. Here, we review some of the advantages and disadvantages of the different approaches, with special concentration on which methods help us to discern the evolutionary process and provide the most accurate estimates of phylogeny.     

Characterization of three glycosyltransferases involved in the biosynthesis of the phenolic glycolipid antigens from the Mycobacterium tuberculosis complex

Mycobacterium tuberculosis and Mycobacterium leprae, the two main mycobacterial pathogens in humans, produce highly specific long chain beta-diols, the dimycocerosates of phthiocerol, and structurally related phenolic glycolipid (PGL) antigens, which are important virulence factors. In addition, M. tuberculosis also secretes glycosylated p-hydroxybenzoic acid methyl esters (p-HBAD) that contain the same carbohydrate moiety as the species-specific PGL of M. tuberculosis (PGL-tb). The genes involved in the biosynthesis of these compounds in M. tuberculosis are grouped on a 70-kilobase chromosomal fragment containing three genes encoding putative glycosyltransferases: Rv2957, Rv2958c, and Rv2962c. To determine the functions of these genes, three recombinant M. tuberculosis strains, in which these genes were individually inactivated, were constructed and biochemically characterized. Our results demonstrated that (i) the biosynthesis of PGL-tb and p-HBAD involves common enzymatic steps, (ii) the Rv2957, Rv2958c, and Rv2962c genes are involved in the formation of the glycosyl moiety of the two classes of molecules, and (iii) the product of Rv2962c catalyzes the transfer of a rhamnosyl residue onto p-hydroxybenzoic acid ethyl ester or phenolphthiocerol dimycocerosates, whereas the products of Rv2958c and Rv2957 add a second rhamnosyl unit and a fucosyl residue to form the species-specific triglycosyl appendage of PGL-tb and p-HBAD. The recombinant strains produced provide the tools to study the role of the carbohydrate domain of PGL-tb and p-HBAD in M. tuberculosis pathogenesis.     

Combining hierarchical clustering and self-organizing maps for exploratory analysis of gene expression patterns

Self-organizing maps (SOM) constitute an alternative to classical clustering methods because of its linear run times and superior performance to deal with noisy data. Nevertheless, the clustering obtained with SOM is dependent on the relative sizes of the clusters. Here, we show how the combination of SOM with hierarchical clustering methods constitutes an excellent tool for exploratory analysis of massive data like DNA microarray expression patterns.     

Construction of a virtual Mycobacterium tuberculosis consensus genome and its application to data from a next generation sequencer

Background

Although Mycobacterium tuberculosis isolates are consisted of several different lineages and the epidemiology analyses are usually assessed relative to a particular reference genome, M. tuberculosis H37Rv, which might introduce some biased results. Those analyses are essentially based genome sequence information of M. tuberculosis and could be performed in sillico in theory, with whole genome sequence (WGS) data available in the databases and obtained by next generation sequencers (NGSs). As an approach to establish higher resolution methods for such analyses, whole genome sequences of the M. tuberculosis complexes (MTBCs) strains available on databases were aligned to construct virtual reference genome sequences called the consensus sequence (CS), and evaluated its feasibility in in sillico epidemiological analyses.

Results

The consensus sequence (CS) was successfully constructed and utilized to perform phylogenetic analysis, evaluation of read mapping efficacy, which is crucial for detecting single nucleotide polymorphisms (SNPs), and various MTBC typing methods virtually including spoligotyping, VNTR, Long sequence polymorphism and Beijing typing. SNPs detected based on CS, in comparison with H37Rv, were utilized in concatemer-based phylogenetic analysis to determine their reliability relative to a phylogenetic tree based on whole genome alignment as the gold standard. Statistical comparison of phylogenic trees based on CS with that of H37Rv indicated the former showed always better results that that of later. SNP detection and concatenation with CS was advantageous because the frequency of crucial SNPs distinguishing among strain lineages was higher than those of H37Rv. The number of SNPs detected was lower with the consensus than with the H37Rv sequence, resulting in a significant reduction in computational time. Performance of each virtual typing was satisfactory and accorded with those published when those are available.

Conclusions

These results indicated that virtual CS constructed from genome sequence data is an ideal approach as a reference for MTBC studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1368-9) contains supplementary material, which is available to authorized users.