Quantifying Reproducibility in Computational Biology: The Case of the Tuberculosis Drugome

How easy is it to reproduce the results found in a typical computational biology paper? Either through experience or intuition the reader will already know that the answer is with difficulty or not at all. In this paper we attempt to quantify this difficulty by reproducing a previously published paper for different classes of users (ranging from users with little expertise to domain experts) and suggest ways in which the situation might be improved. Quantification is achieved by estimating the time required to reproduce each of the steps in the method described in the original paper and make them part of an explicit workflow that reproduces the original results. Reproducing the method took several months of effort, and required using new versions and new software that posed challenges to reconstructing and validating the results. The quantification leads to “reproducibility maps” that reveal that novice researchers would only be able to reproduce a few of the steps in the method, and that only expert researchers with advance knowledge of the domain would be able to reproduce the method in its entirety. The workflow itself is published as an online resource together with supporting software and data. The paper concludes with a brief discussion of the complexities of requiring reproducibility in terms of cost versus benefit, and a desiderata with our observations and guidelines for improving reproducibility. This has implications not only in reproducing the work of others from published papers, but reproducing work from one’s own laboratory.

La Evolucion historica del concepto sobre la tuberculosis pulmonar

Separata de t.2, n.3. de Clinica Hispanica

Desarrollo de un sistema de digitalización y análisis de imágenes de microscopía para el diagnóstico de tuberculosis mediante inteligencia colectiva y gamificación

En este Trabajo de Fin de Grado se diseña, implementa y evalúa un sistema se digitalización de muestras de esputo basado en telefonía móvil e integrable con TuberSpot. Además, se proponen técnicas de procesamiento de imagen para el control de calidad del análisis y se implementa un mecanismo para evaluar la eficiencia de la inteligencia colectiva y la gamificación en este contexto. El sistema de adquisición propuesto utiliza smartphones, adaptadores móvil-microscopio y una aplicación Android. El protocolo de adquisición se ha diseñado conforme a un estudio realizado con personal médico cualificado. El control de calidad se basa en la inserción de bacilos simulados en las imágenes. Para la evaluación de eficiencia de TuberSpot se crea, en colaboración con médicos especialistas, un repositorio de imágenes en las que posición y número de bacilos quedan registrados.

Antibody-mediated inhibition of the growth of larvae from an insect causing cutaneous myiasis in a mammalian host

Many insects feed on blood or tissue from mammalian hosts. One potential strategy for the control of these insects is to vaccinate the host with antigens derived from the insect. The larvae of the fly Lucilia cuprina feed on ovine tissue and tissue fluids causing a cutaneous myiasis associated with considerable host morbidity and mortality. A candidate vaccine antigen, peritrophin 95, was purified from the peritrophic membrane, which lines the gut of these larvae. Serum from sheep vaccinated with peritrophin 95 inhibited growth of first-instar L. cuprina larvae that fed on this serum. Growth inhibition was probably caused by antibody-mediated blockage of the normally semipermeable peritrophic membrane and the subsequent development of an impervious layer of undefined composition on the gut lumen side of the peritrophic membrane that restricted access of nutrients to the larvae. The amino acid sequence of peritrophin 95 was determined by cloning the DNA complementary to its mRNA. The deduced amino acid sequence codes for a secreted protein containing a distinct Cys-rich domain of 317 amino acids followed by a mucin-like domain of 139 amino acids. The Cys-rich domain may be involved in binding chitin. This report describes a novel immunological strategy for the potential control of L. cuprina larvae that may have general application to the control of other insect pests.

Exploring drug-induced alterations in gene expression in Mycobacterium tuberculosis by microarray hybridization

Tuberculosis is a chronic infectious disease that is transmitted by cough-propelled droplets that carry the etiologic bacterium, Mycobacterium tuberculosis. Although currently available drugs kill most isolates of M. tuberculosis, strains resistant to each of these have emerged, and multiply resistant strains are increasingly widespread. The growing problem of drug resistance combined with a global incidence of seven million new cases per year underscore the urgent need for new antituberculosis therapies. The recent publication of the complete sequence of the M. tuberculosis genome has made possible, for the first time, a comprehensive genomic approach to the biology of this organism and to the drug discovery process. We used a DNA microarray containing 97% of the ORFs predicted from this sequence to monitor changes in M. tuberculosis gene expression in response to the antituberculous drug isoniazid. Here we show that isoniazid induced several genes that encode proteins physiologically relevant to the drug’s mode of action, including an operonic cluster of five genes encoding type II fatty acid synthase enzymes and fbpC, which encodes trehalose dimycolyl transferase. Other genes, not apparently within directly affected biosynthetic pathways, also were induced. These genes, efpA, fadE23, fadE24, and ahpC, likely mediate processes that are linked to the toxic consequences of the drug. Insights gained from this approach may define new drug targets and suggest new methods for identifying compounds that inhibit those targets.

Attenuation of virulence in Mycobacterium tuberculosis expressing a constitutively active iron repressor

Iron is an essential nutrient for the survival of most organisms and has played a central role in the virulence of many infectious disease pathogens. Mycobacterial IdeR is an iron-dependent repressor that shows 80% identity in the functional domains with its corynebacterial homologue, DtxR (diphtheria toxin repressor). We have transformed Mycobacterium tuberculosis with a vector expressing an iron-independent, positive dominant, corynebacterial dtxR hyperrepressor, DtxR(E175K). Western blots of whole-cell lysates of M. tuberculosis expressing the dtxR(E175K) gene revealed the stable expression of the mutant protein in mycobacteria. BALB/c mice were infected by tail vein injection with 2 × 105 organisms of wild type or M. tuberculosis transformed with the dtxR mutant. At 16 weeks, there was a 1.2 log reduction in bacterial survivors in both spleen (P = 0.0002) and lungs (P = 0.006) with M. tuberculosis DtxR(E175K). A phenotypic difference in colonial morphology between the two strains also was noted. A computerized search of the M. tuberculosis genome for the palindromic consensus sequence to which DtxR and IdeR bind revealed six putative “iron boxes” within 200 bp of an ORF. Using a gel-shift assay we showed that purified DtxR binds to the operator region of five of these boxes. Attenuation of M. tuberculosis can be achieved by the insertion of a plasmid containing a constitutively active, iron-insensitive repressor, DtxR(E175K), which is a homologue of IdeR. Our results strongly suggest that IdeR controls genes essential for virulence in M. tuberculosis.

Restricted structural gene polymorphism in the Mycobacterium tuberculosis complex indicates evolutionarily recent global dissemination

One-third of humans are infected with Mycobacterium tuberculosis, the causative agent of tuberculosis. Sequence analysis of two megabases in 26 structural genes or loci in strains recovered globally discovered a striking reduction of silent nucleotide substitutions compared with other human bacterial pathogens. The lack of neutral mutations in structural genes indicates that M. tuberculosis is evolutionarily young and has recently spread globally. Species diversity is largely caused by rapidly evolving insertion sequences, which means that mobile element movement is a fundamental process generating genomic variation in this pathogen. Three genetic groups of M. tuberculosis were identified based on two polymorphisms that occur at high frequency in the genes encoding catalase-peroxidase and the A subunit of gyrase. Group 1 organisms are evolutionarily old and allied with M. bovis, the cause of bovine tuberculosis. A subset of several distinct insertion sequence IS6110 subtypes of this genetic group have IS6110 integrated at the identical chromosomal insertion site, located between dnaA and dnaN in the region containing the origin of replication. Remarkably, study of ≈6,000 isolates from patients in Houston and the New York City area discovered that 47 of 48 relatively large case clusters were caused by genotypic group 1 and 2 but not group 3 organisms. The observation that the newly emergent group 3 organisms are associated with sporadic rather than clustered cases suggests that the pathogen is evolving toward a state of reduced transmissability or virulence.

Efficient allelic exchange and transposon mutagenesis in Mycobacterium tuberculosis

A better understanding of Mycobacterium tuberculosis virulence mechanisms is highly dependent on the design of efficient mutagenesis systems. A system enabling the positive selection of insertional mutants having lost the delivery vector was developed. It uses ts-sacB vectors, which combine the counterselective properties of the sacB gene and a mycobacterial thermosensitive origin of replication and can therefore be efficiently counterselected on sucrose at 39°C. This methodology allowed the construction of M. tuberculosis transposition mutant libraries. Greater than 106 mutants were obtained, far exceeding the number theoretically required to obtain at least one insertion in every nonessential gene. This system is also efficient for gene exchange mutagenesis as demonstrated with the purC gene: 100% of the selected clones were allelic exchange mutants. Therefore, a single, simple methodology has enabled us to develop powerful mutagenesis systems, the lack of which was a major obstacle to the genetic characterization of M. tuberculosis.

Conditionally replicating mycobacteriophages: A system for transposon delivery to Mycobacterium tuberculosis

Transposon mutagenesis provides a direct selection for mutants and is an extremely powerful technique to analyze genetic functions in a variety of prokaryotes. Transposon mutagenesis of Mycobacterium tuberculosis has been limited in part because of the inefficiency of the delivery systems. This report describes the development of conditionally replicating shuttle phasmids from the mycobacteriophages D29 and TM4 that enable efficient delivery of transposons into both fast- and slow-growing mycobacteria. These shuttle phasmids consist of an Escherichia coli cosmid vector containing either a mini-Tn10(kan) or Tn5367 inserted into a nonessential region of the phage genome. Thermosensitive mutations were created in the mycobacteriophage genome that allow replication at 30°C but not at 37°C (TM4) or 38.5°C (D29). Infection of mycobacteria at the nonpermissive temperature results in highly efficient transposon delivery to the entire population of mycobacterial cells. Transposition of mini-Tn10(kan) occurred in a site-specific fashion in M. smegmatis whereas Tn5367 transposed apparently randomly in M. phlei, Bacille Calmette–Guérin (BCG), and M. tuberculosis. Sequence analysis of the M. tuberculosis and BCG chromosomal regions adjacent to Tn5367 insertions, in combination with M. tuberculosis genomic sequence and physical map data, indicates that the transpositions have occurred randomly in diverse genes in every quadrant of the genome. Using this system, it has been readily possible to generate libraries containing thousands of independent mutants of M. phlei, BCG, and M. tuberculosis.

Identification of differentially expressed mRNA in prokaryotic organisms by customized amplification libraries (DECAL): The effect of isoniazid on gene expression in Mycobacterium tuberculosis

Understanding the effects of the external environment on bacterial gene expression can provide valuable insights into an array of cellular mechanisms including pathogenesis, drug resistance, and, in the case of Mycobacterium tuberculosis, latency. Because of the absence of poly(A)+ mRNA in prokaryotic organisms, studies of differential gene expression currently must be performed either with large amounts of total RNA or rely on amplification techniques that can alter the proportional representation of individual mRNA sequences. We have developed an approach to study differences in bacterial mRNA expression that enables amplification by the PCR of a complex mixture of cDNA sequences in a reproducible manner that obviates the confounding effects of selected highly expressed sequences, e.g., ribosomal RNA. Differential expression using customized amplification libraries (DECAL) uses a library of amplifiable genomic sequences to convert total cellular RNA into an amplified probe for gene expression screens. DECAL can detect 4-fold differences in the mRNA levels of rare sequences and can be performed on as little as 10 ng of total RNA. DECAL was used to investigate the in vitro effect of the antibiotic isoniazid on M. tuberculosis, and three previously uncharacterized isoniazid-induced genes, iniA, iniB, and iniC, were identified. The iniB gene has homology to cell wall proteins, and iniA contains a phosphopantetheine attachment site motif suggestive of an acyl carrier protein. The iniA gene is also induced by the antibiotic ethambutol, an agent that inhibits cell wall biosynthesis by a mechanism that is distinct from isoniazid. The DECAL method offers a powerful new tool for the study of differential gene expression.

A common mechanism for the biosynthesis of methoxy and cyclopropyl mycolic acids in Mycobacterium tuberculosis

Mycobacterium tuberculosis produces three classes of mycolic acids that differ primarily in the presence and nature of oxygen-containing substituents in the distal portion of the meromycolate branch. The methoxymycolate series has a methoxy group adjacent to a methyl branch, in addition to a cyclopropane in the proximal position. Using the gene for the enzyme that introduces the distal cyclopropane (cma1) as a probe, we have cloned and sequenced a cluster of genes coding for four highly homologous methyl transferases (mma1–4). When introduced into Mycobacterium smegmatis, this gene cluster conferred the ability to synthesize methoxymycolates. By determining the structure of the mycolic acids produced following expression of each of these genes individually and in combination, we have elucidated the biosynthetic steps responsible for the production of the major series of methoxymycolates. The mma4 gene product (MMAS-4) catalyzes an unusual S-adenosyl-l-methionine-dependent transformation of the distal cis-olefin into a secondary alcohol with an adjacent methyl branch. MMAS-3 O-methylates this secondary alcohol to form the corresponding methyl ether, and MMAS-2 introduces a cis-cyclopropane in the proximal position of the methoxy series. The similarity of these reactions and the enzymes that catalyze them suggests that some of the structural diversity of mycolic acids results from different chemical fates of a common cationic intermediate, which in turn results from methyl group addition to an olefinic mycolate precursor.

Molecular basis for the exquisite sensitivity of Mycobacterium tuberculosis to isoniazid

The exceptional sensitivity of Mycobacterium tuberculosis to isonicotinic acid hydrazide (INH) lacks satisfactory definition. M. tuberculosis is a natural mutant in oxyR, a central regulator of peroxide stress response. The ahpC gene, which encodes a critical subunit of alkyl hydroperoxide reductase, is one of the targets usually controlled by oxyR in bacteria. Unlike in mycobacterial species less susceptible to INH, the expression of ahpC was below detection limits at the protein level in INH-sensitive M. tuberculosis and Mycobacterium bovis strains. In contrast, AhpC was detected in several series of isogenic INH-resistant (INHr) derivatives. In a demonstration of the critical role of ahpC in sensitivity to INH, insertional inactivation of ahpC on the chromosome of Mycobacterium smegmatis, a species naturally insensitive to INH, dramatically increased its susceptibility to this compound. These findings suggest that AhpC counteracts the action of INH and that the levels of its expression may govern the intrinsic susceptibility of mycobacteria to this front-line antituberculosis drug.

An epidemic of tuberculosis with a high rate of tuberculin anergy among a population previously unexposed to tuberculosis, the Yanomami Indians of the Brazilian Amazon

A survey of an emerging tuberculosis epidemic among the Yanomami Indians of the Amazonian rain forest provided a unique opportunity to study the impact of tuberculosis on a population isolated from contact with the tubercle bacillus for millennia until the mid-1960s. Within the Yanomami population, an extraordinary high prevalence of active tuberculosis (6.4% of 625 individuals clinically examined) was observed, indicating a high susceptibility to disease, even among bacille Calmette–Guérin-vaccinated individuals. Observational studies on cell-mediated and humoral immune responses of the Yanomami Indians compared with contemporary residents of the region suggest profound differences in immunological responsiveness to Mycobacterium tuberculosis infection. Among the Yanomami, a very high prevalence of tuberculin skin test anergy was found. Of patients with active tuberculosis, 46% had purified protein derivative of tuberculosis reactions <10 mm; similarly 58% of recent bacillus Calmette–Guérin vaccines exhibited skin test reactions <5 mm. The Yanomami also had higher titers of antibodies against M. tuberculosis glycolipid antigens (>70%) than the control subjects comprised of Brazilians of European descent (14%). The antibodies were mostly of the IgM isotype. Among the tuberculosis patients who also produced IgG antibodies, the titers of IgG4 were significantly higher among the Yanomami than in the control population. Although it was not possible to analyze T-cell responses or patterns of lymphokine production in vitro because of the remoteness of the villages from laboratory facilities, the results suggest that the first encounter of the Yanomami Indian population with tuberculosis engenders a diminished cell-mediated immune response and an increased production antibody responses, relative to other populations with extensive previous contact with the pathogen. These findings suggest that tuberculosis may represent a powerful selective pressure on human evolution that over centuries has shaped the nature of human immune responses to infection.

Identification of nitric oxide synthase as a protective locus against tuberculosis

Mutagenesis of the host immune system has helped identify response pathways necessary to combat tuberculosis. Several such pathways may function as activators of a common protective gene: inducible nitric oxide synthase (NOS2). Here we provide direct evidence for this gene controlling primary Mycobacterium tuberculosis infection using mice homozygous for a disrupted NOS2 allele. NOS2−/− mice proved highly susceptible, resembling wild-type littermates immunosuppressed by high-dose glucocorticoids, and allowed Mycobacterium tuberculosis to replicate faster in the lungs than reported for other gene-deficient hosts. Susceptibility appeared to be independent of the only known naturally inherited antimicrobial locus, NRAMP1. Progression of chronic tuberculosis in wild-type mice was accelerated by specifically inhibiting NOS2 via administration of N6-(1-iminoethyl)-l-lysine. Together these findings identify NOS2 as a critical host gene for tuberculostasis.

Modeling the impact of global tuberculosis control strategies

An epidemiological model of tuberculosis has been developed and applied to five regions of the world. Globally, 6.7 million new cases of tuberculosis and 2.4 million deaths from tuberculosis are estimated for 1998. Based on current trends in uptake of the World Health Organization’s strategy of directly observed treatment, short-course, we expect a total of 225 million new cases and 79 million deaths from tuberculosis between 1998 and 2030. Active case-finding by using mass miniature radiography could save 23 million lives over this period. A single contact treatment for tuberculosis could avert 24 million cases and 11 million deaths; combined with active screening, it could reduce mortality by nearly 40%. A new vaccine with 50% efficacy could lower incidence by 36 million cases and mortality by 9 million deaths. Support for major extensions to global tuberculosis control strategies will occur only if the size of the problem and the potential for action are recognized more widely.