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.     

牛分支杆菌与肺结核分支杆菌基因组的比较

通过比较基因组学的方法研究发现,牛分支杆菌与肺结核杆菌基因组的同源性为99．95％,但在牛分枝杆菌基因组中有11个缺失区,大小从1kb到12．7kb,遗传信息的缺失引起牛分枝杆菌的基因组减小;牛分枝杆菌与肺结核分枝杆菌H37Rv间存在着2437个单核苷酸多态性（SNPs）,与肺结核分枝杆菌CDC1551间存在着2423个单核苷酸多态性（SNPs）,牛分支杆菌与肺结核分枝杆菌在编码细胞壁和分泌蛋白上变异程度也是巨大的。研究结果揭示了牛分支杆菌与肺结核分枝杆菌的遗传关系,为研究分支杆菌疫苗和诊断试剂提供理论依据,对牛肺结核病的防治有着非常重要的意义。     

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.     

两株超级广泛耐药结核分枝杆菌全基因组比较分析

【目的】结合现有数据,通过对两株临床超级广泛耐药的结核分枝杆菌全基因组的测序和分析,发现其型别相关的突变位点,解释发生广泛耐药的基因组突变机制。【方法】利用Solexa第二代测序技术对两株广泛耐药结核分枝杆菌(FJ05194和GuangZ0019)进行全基因组测序分析。以H37Rv为参考序列得到两株广泛耐药菌株的单核苷酸多态性(SNPs),构建系统发育树鉴定菌株型别,判断突变位点中型别相关和非型别相关的SNPs。定位SNPs所在的基因组区域,对型别相关的突变基因进行KEGG通路的富集分析,对非型别相关的突变基因和间隔区判断是否与耐药相关。【结果】两株广泛耐药菌株分别属于Lineage2和Lineage4型别,两菌株在碱基替换方面存在差异性,Lineage2型别相关的基因功能富集于ABC转运蛋白和核苷酸切除修复的通路。耐药方面,发现了已知的耐药相关基因的突变(rpoB、katG、rpsl、gyrA、gyrB、embB和ethA等),但卷曲霉素和卡那霉素相关的rrs、tlyA和eis启动子区域未发生突变,不足以解释其耐药性的产生。与最新报道的候选耐药基因比较,发现了卷曲霉素和卡那霉素相关的突变(Rv1393c、Rv0265c和narX等)和外排泵相关的pstB、Rv2333c和Rv2687c突变。【结论】结核分枝杆菌Lineage2型别相关的SNPs中含有影响结核分枝杆菌突变率和耐药性的突变。对于两株超级广泛耐药的结核菌,已知的激活药物或药靶相关的单耐药基因突变集合不能完全解释其广泛耐药性,还涉及新候选结核耐药基因、外排泵和补偿等其他潜在机制的相关基因突变。     

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.     

Whole genome identification of Mycobacterium tuberculosis vaccine candidates by comprehensive data mining and bioinformatic analyses

Background

Mycobacterium tuberculosis, the causative agent of tuberculosis (TB), infects ~8 million annually culminating in ~2 million deaths. Moreover, about one third of the population is latently infected, 10% of which develop disease during lifetime. Current approved prophylactic TB vaccines (BCG and derivatives thereof) are of variable efficiency in adult protection against pulmonary TB (0%–80%), and directed essentially against early phase infection.

Methods

A genome-scale dataset was constructed by analyzing published data of: (1) global gene expression studies under conditions which simulate intra-macrophage stress, dormancy, persistence and/or reactivation; (2) cellular and humoral immunity, and vaccine potential. This information was compiled along with revised annotation/bioinformatic characterization of selected gene products and in silico mapping of T-cell epitopes. Protocols for scoring, ranking and prioritization of the antigens were developed and applied.

Results

Cross-matching of literature and in silico-derived data, in conjunction with the prioritization scheme and biological rationale, allowed for selection of 189 putative vaccine candidates from the entire genome. Within the 189 set, the relative distribution of antigens in 3 functional categories differs significantly from their distribution in the whole genome, with reduction in the Conserved hypothetical category (due to improved annotation) and enrichment in Lipid and in Virulence categories. Other prominent representatives in the 189 set are the PE/PPE proteins; iron sequestration, nitroreductases and proteases, all within the Intermediary metabolism and respiration category; ESX secretion systems, resuscitation promoting factors and lipoproteins, all within the Cell wall category. Application of a ranking scheme based on qualitative and quantitative scores, resulted in a list of 45 best-scoring antigens, of which: 74% belong to the dormancy/reactivation/resuscitation classes; 30% belong to the Cell wall category; 13% are classical vaccine candidates; 9% are categorized Conserved hypotheticals, all potentially very potent T-cell antigens.

Conclusion

The comprehensive literature and in silico-based analyses allowed for the selection of a repertoire of 189 vaccine candidates, out of the whole-genome 3989 ORF products. This repertoire, which was ranked to generate a list of 45 top-hits antigens, is a platform for selection of genes covering all stages of M. tuberculosis infection, to be incorporated in rBCG or subunit-based vaccines.     

Complete annotated genome sequence of Mycobacterium tuberculosis Erdman

We report the completely annotated genome sequence of Mycobacterium tuberculosis Erdman (TMC 107; ATCC 35801), which is a well-known laboratory strain of M. tuberculosis.     

The preponderance of P450s in the Mycobacterium tuberculosis genome

The genome of Mycobacterium tuberculosis (Mtb) encodes 20 different cytochrome P450 enzymes (P450s). P450s are mono-oxygenases, which are historically considered to facilitate prokaryotic usage of unusual carbon sources. However, their preponderance in Mtb strongly indicates crucial physiological functions, as does the fact that polycyclic azoles (known P450 inhibitors) have potent anti-mycobacterial effects. Recent structural and enzyme characterization data reveal novel features for at least two Mtb P450s (CYP121 and CYP51). Genome analysis, knockout studies and structural comparisons signify important roles in cell biology and pathogenesis for various P450s and redox partner enzymes in Mtb. Elucidation of structure, function and metabolic roles will be essential in targeting the P450s as an 'Achilles heel' in this major human pathogen.     

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.     

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.     

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.