Characterization of mycobacteria and mycobacteriophages isolated from compost at the São Paulo Zoo Park Foundation in Brazil and creation of the new mycobacteriophage Cluster U
- James Daltro Lima-Junior1,
- Cristina Viana-Niero2,
- Daniel V. Conde Oliveira1,
- Gabriel Esquitini Machado1,
- Michelle Cristiane da Silva Rabello3,
- Joaquim Martins-Junior4,
- Layla Farage Martins4,
- Luciano Antonio Digiampietri5,
- Aline Maria da Silva4,
- João Carlos Setubal4, 6,
- Daniel A. Russell7,
- Deborah Jacobs-Sera7,
- Welkin H. Pope7,
- Graham F. Hatfull7 and
- Sylvia Cardoso Leão1Email author
© The Author(s). 2016
Received: 25 September 2015
Accepted: 8 June 2016
Published: 17 June 2016
A large collection of sequenced mycobacteriophages capable of infecting a single host strain of Mycobacterium smegmatis shows considerable genomic diversity with dozens of distinctive types (clusters) and extensive variation within those sharing evident nucleotide sequence similarity. Here we profiled the mycobacterial components of a large composting system at the São Paulo zoo.
We isolated and sequenced eight mycobacteriophages using Mycobacterium smegmatis mc2155 as a host. None of these eight phages infected any of mycobacterial strains isolated from the same materials. The phage isolates span considerable genomic diversity, including two phages (Barriga, Nhonho) related to Subcluster A1 phages, two Cluster B phages (Pops, Subcluster B1; Godines, Subcluster B2), three Subcluster F1 phages (Florinda, Girafales, and Quico), and Madruga, a relative of phage Patience with which it constitutes the new Cluster U. Interestingly, the two Subcluster A1 phages and the three Subcluster F1 phages have genomic relationships indicating relatively recent evolution within a geographically isolated niche in the composting system.
We predict that composting systems such as those used to obtain these mycobacteriophages will be a rich source for the isolation of additional phages that will expand our view of bacteriophage diversity and evolution.
Bacteriophages (phages) are the most abundant life forms on the planet, and can influence the entire bacterial ecosystem . They are capable of delivering, through horizontal gene transfer, mechanisms of adaptation and resistance, toxins and photosynthesis genes to their hosts. They also contribute to the cycling of organic and inorganic nutrients . According to the International Committee on Taxonomy of Viruses, double-stranded deoxyribonucleic acid (dsDNA) tailed phages constitute one Order (Caudovirales), three Families (Myoviridae, Podoviridae and Siphoviridae) and 39 Genera (http://www.ictvonline.org/virusTaxonomy.asp?taxnode_id=20120366) [3, 4].
Mycobacteriophages infect members of the genus Mycobacterium and were first isolated in the 40’s and 50’s [5, 6]. They were initially used as tools in the study of bacterial genetics and for diagnosis and typing of Mycobacterium tuberculosis [7–11]. More than 800 mycobacteriophage genomes have been sequenced (http://phagesdb.org) a collection that has grown rapidly due in part to the impact of integrated research-education programs [12, 13].
All mycobacteriophages isolated to date are dsDNA tailed phages and morphologically are Siphoviridae or Myoviridae . Mycobacteriophages have been isolated from clinical samples of diseased and healthy people [14, 15], animals , laboratory cultures [17, 18], but primarily from the environment, especially from soil to compost samples [19–21]. Most known mycobacteriophages were isolated using Mycobacterium smegmatis mc2155 as a host, comprising the largest collection of sequenced viruses infecting a common host . However, there is a geographical disequilibrium in mycobacteriophage isolation, and only about 7 % were isolated outside United States of America (USA) (http://phagesdb.org). The analysis of sequenced genomes shows exchange of genome segments among the phage population, revealing mosaicism as the main characteristic of such genomes . Nonetheless, genome comparison allowed the grouping of phages showing significant genomic similarities in clusters and subclusters .
Composting is a controlled process of plant and other living material decomposition by autochthone organisms that results in the production of compost to be used as fertilizer. Macro and microorganisms, as common bacteria, fungi, protozoan, algae, larvae and arthropods can participate in the degradation of organic matter. However, composting processes are typically carried out mainly by complex microbial communities . Such intense microbial activity promotes orchestrated changes in temperature and potential of hydrogen (pH) . The high diversity of organisms present at different times during the composting process offers an important source of environmental microorganisms and phages for study.
The São Paulo Zoo Park Foundation (FPZSP) has an area of 825,000 square meters (m2) and maintains an exhibition fauna of more than 3,200 animals, as well as the free-living fauna living in the neighboring Atlantic rain forest. The park has operated since 2007 as an environmental management program, including an effluent treatment plant, a water treatment plant and an organic compost production unit, which receives four tons/day of organic residues from Zoo activities and solid waste products from the effluent and water treatment stations. The humus-rich compost generated in the unit is used for the production of food for the Zoo animals in the São Paulo Zoo Farm and for fertilization of the park. The Sao Paulo Zoo composting process has been shown to hold a remarkable bacterial diversity, which can fully account for the biomass degradation . The composting operation was started in 2004 and processes 120 ton/month of organic waste for compost production. Its implementation allowed the recovery and preservation of environmental conditions of the park, reduction of fertilization products costs and better control of the food offered to the Zoo animals. Indirectly, the park obtains financial benefits.
The objectives of this work were to evaluate if compost from São Paulo Zoo Park is a suitable source for isolation of mycobacteria and mycobacteriophages, to understand the dynamics of mycobacterial and phage isolation during the compost processing and to compare isolated phages with mycobacteriophages from other countries. We report the genome sequences of eight newly isolated phages providing insights into viral diversity and the first evidence for mycobacteriophage evolution within a geographically isolated niche.
Results and discussion
Isolation and identification of mycobacteria from composting material
Cell composition, temperature (mean ± SD) and pH of composting chambers
chamber 1 vegetables, fruits, manure, beddings and food residues, plant debris and grass clippings, australian goose carcass
chamber 2 vegetables, fruits, manure, beddings and food residues, plant debris and grass clippings, animal carcass and waste water treatment sludge
chamber 3 mixture of partially composted material from chambers 1 and 2
Air T (o C)
T (o C)
T (o C)
T (o C)
57,8 ± (16,9)
55,6 ± (14,7)
56,0 ± (15,7)
63,4 ± (17,5)
53,0 ± (14,5)
65,0 ± (18,2)
51,0 ± (13,6)
55,2 ± (15,7)
53 ± (15,5)
62,4 ± (21,5)
58,2 ± (16,5)
Results of identification by PRA-hsp65, and by rpoB and hsp65 sequencing of 38 mycobacterial isolates recovered from composting chambers
M. fortuitum 1a
M. fortuitum T (JF346874)
M. fortuitum T (AF547833)
M. fortuitum 2a
M. porcinum T (AY262737)
M. houstonense T (DQ987725)
M. peregrinum 1
M. septicum (HM807423)
M. septicum T (AF547873)
M. peregrinum 3
M. alvei T (AY859697)
M. peregrinum T (AF547865)
M. septicum 1b
M. septicum (HM807423)
M. septicum T (AY496142)
M. kumamotonense 1a
M. kumamotonense T (JN571258)
M. kumamotonense T (JF491323)
M. terrae 1
M. senuense (JN571250)
M. algericum T (GU564405)
M. insubricum T (EU022519)
M. insubricum T (JF491319)
The isolate showing the PRA-hsp65 pattern of M. terrae 1 was identified as Mycobacterium senuense by rpoB sequencing and as Mycobacterium algericum by hsp65 sequencing. These species, and M. kumamotonense, whose identification was confirmed by rpoB and hsp65 sequencing, belong to the M. terrae complex [33, 34], which was the second group most prevalent in compost samples, corresponding to 21 % (8/38) of all isolates.
The isolate showing a novel PRA-hsp65 pattern not present in the PRASITE database was identified as Mycobacterium insubricum by rpoB and hsp65 sequencing. This recently described species was isolated from sputum samples but was not considered clinically significant  and from fish . This is the first description of its isolation from environmental samples and PRA-hsp65 pattern.
Isolation and purification of mycobacteriophages
Mycobacteriophages were isolated from chambers 1 and 2 using M. smegmatis mc2155 as host. Crude extracts from the first collection (Mar 11, 2011) yielded 627 lytic plaques with material from chamber 1 and 484 with material from chamber 2. In the second collection (Mar 25, 2011) the number of lytic plaques was reduced to 14 in chamber 1 and 74 in chamber 2. Additional crude extracts collected at intervals of 15 days, were used for phage isolation, but after the fourth collection no lytic plaques were detected, showing that mycobacteriophages able to infect M. smegmatis were more prevalent in the initial phases of the process. The same experiment was repeated using two fresh samples from a different chamber, collected at 15-day intervals in the initial phase of the composting process, and 150 and 290 lytic plaques were obtained, respectively. Eight mycobacteriophages randomly selected were purified and stored.
We concluded that compost is a good source of mycobacteriophages capable of infecting M. smegmatis mc2155. To investigate if mycobacterial isolates from compost could also be used as hosts for isolation of mycobacteriophages from compost we used eight isolates representing the different observed PRA-hsp65 patterns (M. fortuitum 1, M. fortuitum 2, M. peregrinum 1, M. peregrinum 3, M. septicum 1, M. kumamotonense 1, M. terrae 1 and M. insubricum) as hosts. Each isolate was infected with the stored material from chambers 1 and 2 and fresh material from two additional chambers. No lytic plaques were obtained in these experiments. Control infections using M. smegmatis mc2155 as host in the same experiments were successful in generating lytic plaques (data not shown). Moreover, the eight purified phages previously isolated from chambers 1 and 2 were not able to infect the mycobacterial isolates from compost. The inability of phages to infect co-existing strains has been also witnessed in other environments, such as seawater cyanophages .
Mycobacteria isolated in this work may not include the natural hosts of mycobacteriophages from compost. Less than 1 % of environmental prokaryotes can grow in controlled conditions  and some mycobacteria can be difficult to recover from the environment. Moreover, few known phages infect the mycobacterial species identified in this work. Phage F-phi WJ-1 is the only described mycobacteriophage isolated from M. fortuitum . Additionally, the protocol used for phage isolation, obtained at http://phagesdb.org/workflow/ included an initial step of filtration of the composting material in 0.22 micrometer (μm) membranes in order to eliminate bacteria and other organisms. This implies that the isolated phages were present as free particles in the material. Occasional phages inside mycobacteria, either free or integrated in the genome, were most probably eliminated in the filtration step. Interestingly, the metagenomic analysis of two composting cells from the FPZSP, distinct from the ones analyzed here, revealed a low amount of viral sequences, corresponding to 0.05 and 0.25 % of the obtained sequences in each chamber .
Sequence of ends (5’–3’)
The eight phages sequenced here were compared to metagenomic reads obtained from a sample that included material from five different days during the composting process of a single chamber at FPZSP (unpublished work). Some hits were observed, particularly with phage Florinda, in sequences of 583 bp, 460 bp and 352 bp, all with at least 99 % identity. The evidence presented suggests that mycobacteriophage sequences can be at least partially recovered by metagenomics.
By comparison with sequences available in the Mycobacteriophages Database at http://phagesdb.org/ and based on the classification in clusters and subclusters proposed by Hatfull et al. , two phages (Nhonho and Barriga) were grouped in Subcluster A1, two phages (Pops and Godines) in Cluster B, one (Pops) from subcluster B1 and the other (Godines) in B2 and three (Florinda, Girafales and Quico) were grouped in Subcluster F1. One phage (Madruga) has high DNA sequence similarity with singleton phage Patience , and these two phages were grouped to form the new Cluster U. Functional assignments from database matches, HHPred results, and syntenic placement can be made for 25.6 % of the genes in these genomes. Across the eight genomes, ten new orphams have been identified, four in Subcluster A1 phages Barriga and Nhonho, and six in Subcluster F1 phages, Florinda, Giraffes, and Quico. The genomic features of each of these groups will be discussed in further detail.
Subcluster A1 mycobacteriophages Nhonho and Barriga
Although Nhonho and Barriga have similar genomic architectures, there are several informative differences between them (Fig. 2). For example, in the Barriga right arm there are eight open reading frames of unknown function (72 – 79) that are absent from Nhonho, but present in a small subset of other Subcluster A1 phages. We also note that Barriga 84 is replaced by gene 73 in Nhonho, an orphan with no close relatives in the mycobacteriophage database. Finally, we note that Barriga gp52, which is a homologue of the Bxb1 Recombination Directionality Factor (RDF) gp47  lacks an intein that is present in the homologue (gene 52) in Nhonho.
Mycobacteriophages Pops and Godines from Cluster B
The Subcluster B2 phage Godines is most closely related to phages Ares and Arbiter (Fig. 4). Godines shares the common Subcluster B2 features such as the lack of Lysin B gene [47–49]. Another characteristic of this group is the presence of six genes between the gene for terminase and the left genome end (Fig. 4), whose products are predicted to be involved in the synthesis of queuosine from guanosine-5’-triphosphate (GTP).
Informative differences between Godines and other B2 genomes include the insertion of a small gene (56) of unknown function between the putative primase/helicase and DNA polymerase genes, which is has only one other relative present in the mycobacteriophages, in the Subcluster B5 phage Baee (Fig. 4). Also, a gene present in both Ares and Arbiter (77 and 74 respectively) is absent from Godines, and a Godines gene (89) is present in Ares (gene 90) but is absent from Arbiter; presumably these genes are not required for lytic growth.
Pops and Godines share only low DNA sequence similarly to each other and share fewer than 25 % of their genes through amino acid sequence comparisons. They are each more closely related to other phages than to each and thus there is no evidence to suggest they have evolved in evolutionary isolation in the geographic area from where they were isolated.
Cluster F mycobacteriophages Florinda, Girafales and Quico
Other notable features of these phages include the 7–8 genes – including the terminase, portal, and capsid subunit genes – at the extreme left end of the genomes (Fig. 5). The Cluster F genomes have two distinct types of these genes that have unrelated DNA sequences but encode similar functions: Florinda has one distinct type (related to genes in phages such as Daenerys and Llij), whereas Girafales and Quico have the alternative type (related to phages Boomer and Ramsey). Secondly, there is a small group of genes in Girafales (58, 59, 62) that substitute for unrelated genes in Florinda and Quico (Fig. 5). We also note that a small gene (115, 55 codons) present at the right end of the Florinda genome is absent from both Quico and Girafales, although found in otherwise unrelated genomes in Subclusters A1, D1 and D2. Its function is unknown (Fig. 5).
Mycobacteriophage Madruga and the novel Cluster U
First, Madruga shares the low GC% content (50.4 %) to that of Patience, and the genomes are of similar lengths (69,377 bp and 70,506 bp in Madruga and Patience respectively). The length difference is primarily accounted for by three genes that are present in Patience but absent from Madruga: Patience genes 5, 7 and 104 (Fig. 7). Patience 5 is a putative histidine-asparagine-histidine motif (HNH) endonuclease and there are distantly related copies in Cluster C phages Myrna and Dandelion (39 % and 33 % aa identity respectively); the protein is expressed both early and late in lytic growth . Patience 7 is related to the Holliday Junction resolvase endonuclease VII which is expressed late in lytic growth  and although this function is likely required for DNA packaging, we note that Patience and Madruga both encode another Holliday Junction resolvase (RusA-like; genes 91 and 93 respectively) elsewhere in the genome that is expressed late in lytic growth and could provide functional redundancy (Fig. 7). Finally, Patience gene 104 is absent from Madruga and was shown to be expressed late in Patience lytic growth but is presumably not required; its function is unknown.
Finally, we note that putative tails proteins in both Madruga and Patience (gp34 and gp36 respectively) contain a C-terminal 150-residue domain that is also present at the C-termini of putative tail proteins encoded by other phages, including Pops (and other Subcluster B1 phages), and phages in Subcluster D1, D2, H1, H2, and Cluster R (Fig. 7). HHPred analysis strongly indicates that this domain has a carbohydrate binding function, perhaps with Endo-1,4-beta-xylanase activity. Presumably this is associated with host surface recognition or capsule degradation to facilitate the early stages of infection.
In conclusion, this work provides new insights into the genetic diversity of mycobacteriophages, and presents supporting evidence for local evolution of groups of phages (Florinda, Quico and Girafales; Barriga and Nhonho) in relative isolation from other mycobacteriophages. Although the database of sequenced mycobacteriophage genomes is quite large, this is the first clear evidence for this phenomenon, and may reflect in part the richness and growth of the microbial communities within the zoo composting systems from which they were isolated. Interestingly, none of the phages isolated here appear to infect any of the mycobacterial strains isolated from the same locations, and thus the natural host for these phages is unclear. We predict that composting systems such as those used to obtain these mycobacteriophages will be a rich source for the isolation of additional phages that will expand our view of bacteriophage diversity and evolution.
Composting material from 8 cubic meters (m3) chambers was collected for this work. Two chambers (1 and 2) were followed and four collections were done from March 11 to April 28, 2011. After this period the partially composted material from both cells was reunited in a third chamber, which was followed until June 10, with three additional collections. This is the regular procedure of compost production carried out at the FPZSP. At each collection, the temperature of five points located 1.5 meters (m) bellow the cell surface was checked and the mean temperature of the five points was calculated. The pH of the collected material was measured using a pH measure strip (EM Science, Germany). The composition, mean temperature and pH of each chamber and the air temperature at each collection are shown in Table 1. At each time point, material from five different points from each chamber was collected and pooled as previously described . The material was aliquoted in two 50 mililitre (mL) tubes, one used for mycobacteria and mycobacteriophage isolation and the other for storage at −20 °C.
M. smegmatis mc2155 was cultivated on solid Middlebrook 7H10 medium (Becton Dickinson and Company, Franklin Lakes, NJ, USA) supplemented with 10 % OADC (oleic acid, albumin, dextrose and catalase - Becton Dickinson and Co). Isolated colonies were then transferred to liquid Middlebrook 7H9 medium (Becton Dickinson and Co) supplemented with 10 % OADC, Tween 80® 0.05 % and CaCl2 1 mM (Merck, Darmstadt, HE, Germany) and cultivated at 37 °C for 24 h. The culture was diluted 1:20 in the same medium without Tween 80 until 0.5 < D.O.600nanometers(nm) < 0.7 (log phase) was reached.
Isolation of mycobacteria from compost
Isolation of mycobacteria was carried out as described . The collected material [10 gram (g)] was diluted in 40 mL of sterile water and maintained in a shaker for 1 h (h) at room temperature (RT). The solution was centrifuged at 600 force of gravity (g) for 5 min (min) at 4 ° C for precipitation of large particles and the supernatant was transferred to sterile tubes and centrifuged at 8000 g for 15 min at 4 ° C. Sodium dodecyl sulfate (SDS) 3 %/NaOH 4 % were added to the pellet, which was incubated for 15 min at RT. The sample was centrifuged again under the same conditions and the pellet was treated with a solution of cetrimide (Henrifarma Produtos Químicos e Farmacêuticos, São Paulo, Brazil) 2 % for 5 min. After centrifugation under the same conditions, the pellet was washed and resuspended in 0.5 mL of deionized water. An aliquot of 0.1 mL was distributed in duplicates on Löwenstein-Jensen (Probac do Brasil, São Paulo, Brazil) slants and on 7H10-OADC plates supplemented with PANTA (antibiotic solution containing final concentrations of 40 U/mL polymyxin, 4 μg/mL amphotericin B, 16 μg/mL nalidixic acid, 4 μg/mL trimethoprim and 4 μg/mL azlocillin) . The cultures were incubated at 30 ° C and 37 ° C and examined every 2 days in the first 10 days and once a week for a total of 90 days.
Identification of mycobacteria
Isolated colonies were visualized by microscopy after Ziehl-Neelsen staining. Acid-fast bacilli were identified by PRA-hsp65. DNA was extracted by boiling a loop full of acid-fast bacilli taken from solid plates in 300 microlitres (μL) of TET [Tris–HCl 10 milimolar (mM), EDTA 1 mM, pH 8.0, Tween 20® 0.05 %] for 10 min. A 441 bp fragment from the hsp65 gene was amplified using primers Tb11 (5’-ACCAACGATGGTGTGTCCAT-3') and Tb12 (5’-CTTGTCGAACCGCATACCCT-3') as previously described . Amplified products were separately digested with BstE II at 60 ° C and Hae III at 37 ° C for 90 min. Digestion products were separated by electrophoresis in 3 % agarose gels stained with ethidium bromide and the resulting bands were analyzed by comparison with restriction profiles at the PRASITE (http://app.chuv.ch/prasite/index.html).
The identification of 11 selected isolates was confirmed by DNA sequencing. Partial sequences of rpoB and hsp65 genes were obtained. The primers used for amplification and sequencing of rpoB were MycoF (5'-GGCAAGGTCACCCCGAAGGG-3') and MycoR (5'-AGCGGCTGCTGGGTGATCATC-3')  and for hsp65 were hsp667-forward (5'-GGCCAAGACAATTGCGTACG-3') and hsp667-reverse (5'-GGAGCTGACCAGCAGGATG-3') . The amplicons were purified using QIAquick PCR purification Kit (Qiagen, Germany). Dideoxy sequencing was performed using BigDye® Terminator v 3.1 Cycle Sequencing Kit (Applied Biosystems, USA) and run in ABI PRISM 3100 DNA Analyzer (Applied Biosystems). The sequences were analysed by comparison with sequences deposited in the GenBank using the Basic Local Alignment Search Tool (BLAST: http://blast.ncbi.nlm.nih.gov/Blast.cgi).
Approximately 30 g of collected material were diluted 1:1 weight/volume percent (w/v) in phage buffer (Tris 10 mM pH 7.5, MgSO4 10 mM, NaCl 68 mM, CaCl2 0.1 mM). After homogenization and incubation at RT for 30 min, the samples were centrifuged at 3,500 g for 15 min (Sorvall, Newtown, CT, USA) and the liquid phase was filtered through a 0.22 μm pore membrane (Techno Plastic Products AG-TPP, Trasadingen, SH, Switzerland). A volume of 10 μL of the filtered material was incubated with 0.5 mL of log-phase cultures of M. smegmatis mc2155 or mycobacterial isolates from compost at RT for 30 min and added to 2.5 mL of Middlebrook top agar (MBTA)/CaCl2 (7H9 liquid medium with agar 0.7 % and CaCl2 0.1 mM) and 2 mL of 7H9-OADC. The final volume was homogenized and poured over solid 7H10-OADC/CaCl2 plates, which were incubated for 16 h at 37 ° C. Negative controls (no addition of compost filtered material) and positive controls (phage D29) were included in each experiment. After the incubation period, putative lytic plaques showing different sizes and morphologies were collected by removing the agar regions using micropipet tips and transferred to microcentrifuge tubes containing 100 μL of phage buffer. Several rounds of purification were performed by infection of M. smegmatis mc2155 until all lytic plaques in each sample showed the same size and morphology. Each purified phage was then amplified and titered. The amplified phage samples were harvested, filtered through a 0.22 μm membrane and ultracentrifuged at 65,000 g at 4 °C for 2 h (Beckman-Coulter Optima XL 100 K, Pasadena, CA, US). The pellet was suspended in 1 mL of ammonium acetate 0.1 molar (M) (Synth, Brazil), pH 7.5 or in phage buffer.
Whole genome sequencing (WGS), assembly and annotation of phages genomes
DNA from purified phages was submitted to pyrosequencing using standard Roche 454 GS FLX Titanium protocols (Roche Applied Science), at the Center for Advanced Technologies in Genomics (CATG), Instituto de Quimica, Universidade de Sao Paulo. Barcoded (tagged) shotgun libraries for each DNA sample were constructed using GS Titanium Rapid Library Prep Kit, pooled and submitted to two sequencing runs. Sequencing reads were separated according to their individual tags, quality-filtered and individually assembled using 454 Newbler assembler software version 2.5.3. Automatic annotation was performed using DNA Master software v. 5.0.2 (JG Lawrence, http://cobamide2.bio.pitt.edu/computer.htm), followed by manual curation. Hypothetical protein genes were analyzed for the presence of conserved domains using online tools CDD-NCBI  and InterPro-EMBL-EBI . Genome maps were generated using Phamerator  with the database ‘Actinobacteriophage_554’. This database contains 554 phage genomes of phages infecting Actinobacterial hosts, consisting of 56,572 genes, organized into a total of 8032 phamilies (phams) composed of sequence-related proteins. Parameters for pham construction were as described previously . Network relationships were derived from the Phamerator Pham table, and represented using Splitstree selecting only the Cluster F1 phages.
BLAST, basic local alignment search tool; BLASTN, nucleotide basic local alignment search tool; bp, base-pair; CATG, center for advanced technologies in genomics; DNA, deoxyribonucleic acid; dsDNA, double-stranded deoxyribonucleic acid; FPZSP, Sao Paulo Zoo Park Foundation; g, force of gravity; g, gram; GC%, percent of guanine and cytosine; GTP, guanosine-5’-triphosphate; h, hour; HNH, histidine-asparagine-histidine motif; PCR, polymerase chain reaction; hsp65, heat shock protein 65 gene; m, meter; M, molar; M., Mycobacterium; m2, square meter; m3, cubic meter; Mbp, mega base-pair; MBTA, Middlebrook top agar; min, minute; mL, mililitre; mM, milimolar; nm, nanometer; nt, nucleotide; OADC, oleic acid, albumin, catalase and dextrose supplement; ORF, open reading frame; PANTA, antibiotic solution containing final concentrations of 40 U/mL polymyxin, 4 μg/mL amphotericin B, 16 μg/mL nalidixic acid, 4 μg/mL trimethoprim and 4 μg/mL azlocillin; PCR, polymerase chain reaction; pH, potential of hydrogen; PRA-hsp65, polymerase chain reaction restriction enzyme analysis of the hsp65 gene; rpoB, gene that encodes the beta subunit of bacterial ribonucleic acid polymerase; RT, room temperature; SDS, sodium dodecyl sulfate; USA, United States of America; w/v, weight/volume percent; WGS, whole genome sequencing; μL, microlitre; μm, micrometer
We acknowledge the staff from FPZSP for helping with the collections. We are very grateful to Drs Joao Batista da Cruz and Paulo Bressan, respectively Scientific Director and President of the São Paulo Zoo Park Foundation (FPZSP) for their continued support to this project. Prof. Mario Hiroyuki Hirata is acknowledged for providing phage D29 for this study. We thank Charles Bowman and Steve Cresawn for assistance in database construction.
This study received financial support from São Paulo Research Foundation (FAPESP) (grants 2011/18326-4 and 2011/50870-6), and from the Howard Hughes Medical Institute grant 54308198 to GFH. JDL-Jr and GEM received fellowships from FAPESP (2012/04911-5 and 2014/01825-6). DVCO received undergraduate fellowship from National Council for Scientific and Technological Development (CNPq). JM-Jr is supported by a fellowship from Coordination for the Improvement of Higher Education Personnel (CAPES). LAD, AMDS, JCS and SCL are partially supported by Research Fellowship Awards from National Council for Scientific and Technological Development (CNPq). Funding agencies had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Availability of data and materials
The mycobacterial isolates and mycobacteriophages were deposited in the collection of the Sao Paulo Zoo Park Foundation (FPZSP) and in the collection of the Laboratory of Mycobacteria, Departamento de Microbiologia, Imunologia e Parasitologia, Escola Paulista de Medicina, Universidade Federal de São Paulo, São Paulo, Brazil. The material is available for non-commercial purposes upon request.
JDL-Jr carried out mycobacteriophage isolation, identification of mycobacteria and helped to draft the manuscript, CV-N carried out mycobacterial isolation and helped to draft the manuscript, DVCO carried out identification of mycobacteria and experiments of isolation of phages using mycobacteria from compost, GEM helped in the identification of mycobacteria, MCSR helped to standardize the experiments of phage isolation in the laboratory, JM-Jr participated in sequence assembly, LFM carried out sequencing of the mycobacteriophages, LAD participated in sequence assembly, AMDS participated in coordination of phage sequencing and helped to draft the manuscript, JCS coordinated phage assembly and sequence analysis and helped to draft the manuscript, DAR carried out sequence review and analysis, DJ-S carried out annotation reviews and prepared figures, WHP reviewed the analysis of the mycobacteriophages, GFH coordinated the review of the analysis and helped to draft the manuscript, SCL conceived the study and participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.
JDL-Jr, DVCO, MCSR and GEM were students in the Laboratory of Mycobacteria, Universidade Federal de São Paulo; CV-N and SCL are Associate Professors in the same University and have a long-lasting collaboration in projects related to Molecular Biology of Mycobacteria; JM-Jr, LFM, LAD, AMDS and JCS belong to a group that work at Instituto de Quimica, Universidade de São Paulo in projects related to WGS and Bioinformatics. DAR, DJ-S, WHP and GFH are members of the Pittsburgh Bacteriophage Institute and have great expertise in mycobacteriophage characterization and analysis.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
This study has been approved by the Ethics Committee from Universidade Federal de Sao Paulo in May 11, 2012, under number 0068/12.
Nucleotide sequence accession numbers
Complete genome sequences for eight phages obtained in this work were deposited in GenBank under accession numbers: [GenBank:KR997929, GenBank:KR997930, GenBank:KR997931, GenBank:KR997932, GenBank:KR997933, GenBank:KR997934, GenBank:KR997967 and GenBank:KR997968].
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Hambly E, Suttle CA. The viriosphere, diversity, and genetic exchange within phage communities. Curr Opin Microbiol. 2005;8(4):444–50. doi:10.1016/j.mib.2005.06.005.View ArticlePubMedGoogle Scholar
- Rohwer F, Thurber RV. Viruses manipulate the marine environment. Nature. 2009;459(7244):207–12. doi:10.1038/nature08060.View ArticlePubMedGoogle Scholar
- Ackermann HW. Bacteriophage observations and evolution. Res Microbiol. 2003;154(4):245–51. doi:10.1016/S0923-2508(03)00067-6.View ArticlePubMedGoogle Scholar
- Maniloff J, Ackermann HW. Taxonomy of bacterial viruses: establishment of tailed virus genera and the order Caudovirales. Arch Virol. 1998;143(10):2051–63.View ArticlePubMedGoogle Scholar
- Gardner GM, Weiser RS. A bacteriophage for Mycobacterium smegmatis. Proc Soc Exp Biol Med. 1947;66:205–6.View ArticlePubMedGoogle Scholar
- Whittaker E. Two bacteriophages for Mycobacterium smegmatis. Can J Public Health. 1950;41:431–6.PubMedGoogle Scholar
- Bardarov S, Kriakov J, Carriere C, Yu S, Vaamonde C, McAdam RA, et al. Conditionally replicating mycobacteriophages: a system for transposon delivery to Mycobacterium tuberculosis. Proc Natl Acad Sci U S A. 1997;94(20):10961–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Jacobs Jr WR, Kalpana GV, Cirillo JD, Pascopella L, Snapper SB, Udani RA, et al. Genetic systems for mycobacteria. Methods Enzymol. 1991;204:537–55.View ArticlePubMedGoogle Scholar
- Kubica GP. Phage typing of Mycobacterium tuberculosis: a time for standardization. Am Rev Respir Dis. 1982;126(1):3–4.PubMedGoogle Scholar
- Rado TA, Bates JH, Engel HW, Mankiewicz E, Murohashi T, Mizuguchi Y, et al. World Health Organization studies on bacteriophage typing of mycobacteria. Subdivision of the species Mycobacterium tuberculosis. Am Rev Respir Dis. 1975;111(4):459–68.PubMedGoogle Scholar
- Rondon L, Piuri M, Jacobs Jr WR, de Waard J, Hatfull GF, Takiff HE. Evaluation of fluoromycobacteriophages for detecting drug resistance in Mycobacterium tuberculosis. J Clin Microbiol. 2011;49(5):1838–42. doi:10.1128/JCM.02476-10.View ArticlePubMedPubMed CentralGoogle Scholar
- Hatfull GF. Mycobacteriophages: genes and genomes. Annu Rev Microbiol. 2010;64:331–56. doi:10.1146/annurev.micro.112408.134233.View ArticlePubMedGoogle Scholar
- Hatfull GF, Science Education Alliance Phage Hunters Advancing G, Evolutionary Science P, KwaZulu-Natal Research Institute for T, Students HIVMGC, Phage Hunters Integrating R, et al. Complete genome sequences of 138 mycobacteriophages. J Virol. 2012;86(4):2382–4. doi:10.1128/JVI.06870-11.View ArticlePubMedPubMed CentralGoogle Scholar
- Mankiewicz E. Mycobacteriophages isolated from persons with tuberculous and non-tuberculous conditions. Nature. 1961;191:1416–7.View ArticlePubMedGoogle Scholar
- Mankiewicz E, Liivak M. Mycobacteriophages isolated from human sources. Nature. 1967;216(5114):485–6.View ArticlePubMedGoogle Scholar
- Cambir S. Some data concerning the distribution of mycobacteriophages in nature (soil, manure, excrements of birds and of bovines, etc.). Arch Roum Pathol Exp Microbiol. 1966;25(2):396–401.PubMedGoogle Scholar
- Sassi M, Bebeacua C, Drancourt M, Cambillau C. The first structure of a mycobacteriophage, the Mycobacterium abscessus subsp. bolletii phage Araucaria. J Virol. 2013;87(14):8099–109. doi:10.1128/JVI.01209-13.View ArticlePubMedPubMed CentralGoogle Scholar
- Sula L, Sulova J, Matejka M, Malkova J. Isolation of mycobacterial phage from the laboratory strain Mycobacterium leprae murium "Douglas". Czech Med. 1984;7(3):174–80.PubMedGoogle Scholar
- Fan X, Teng T, Wang H, Xie J. Biology of a novel mycobacteriophage, SWU1, isolated from Chinese soil as revealed by genomic characteristics. J Virol. 2012;86(18):10230–1. doi:10.1128/JVI.01568-1286/18/10230.View ArticlePubMedPubMed CentralGoogle Scholar
- Henry M, O'Sullivan O, Sleator RD, Coffey A, Ross RP, McAuliffe O, et al. In silico analysis of Ardmore, a novel mycobacteriophage isolated from soil. Gene. 2010;453(1–2):9–23. doi:10.1016/j.gene.2009.12.007.View ArticlePubMedGoogle Scholar
- Mageeney C, Pope WH, Harrison M, Moran D, Cross T, Jacobs-Sera D, et al. Mycobacteriophage Marvin: a new singleton phage with an unusual genome organization. J Virol. 2012;86(9):4762–75. doi:10.1128/JVI.00075-12.View ArticlePubMedPubMed CentralGoogle Scholar
- Pope WH, Bowman CA, Russell DA, Jacobs-Sera D, Asai DJ, Cresawn SG, et al. Whole genome comparison of a large collection of mycobacteriophages reveals a continuum of phage genetic diversity. eLife. 2015;4:e06416. doi:10.7554/eLife.06416.View ArticlePubMedPubMed CentralGoogle Scholar
- Hatfull GF, Pedulla ML, Jacobs-Sera D, Cichon PM, Foley A, Ford ME, et al. Exploring the mycobacteriophage metaproteome: phage genomics as an educational platform. PLoS Genet. 2006;2(6):e92.View ArticlePubMedPubMed CentralGoogle Scholar
- Ryckeboer J, Mergaert J, Coosemans J, Deprins K, Swings J. Microbiological aspects of biowaste during composting in a monitored compost bin. J Appl Microbiol. 2003;94(1):127–37.View ArticlePubMedGoogle Scholar
- Kumar S. Composting of municipal solid waste. Crit Rev Biotechnol. 2011;31(2):112–36. doi:10.3109/07388551.2010.492207.View ArticlePubMedGoogle Scholar
- Martins LF, Antunes LP, Pascon RC, de Oliveira JC, Digiampietri LA, Barbosa D, et al. Metagenomic analysis of a tropical composting operation at the sao paulo zoo park reveals diversity of biomass degradation functions and organisms. PLoS One. 2013;8(4):e61928. doi:10.1371/journal.pone.0061928.View ArticlePubMedPubMed CentralGoogle Scholar
- Kitamura R, Ishii K, Maeda I, Kozaki T, Iwabuchi K, Saito T. Evaluation of bacterial communities by bacteriome analysis targeting 16S rRNA genes and quantitative analysis of ammonia monooxygenase gene in different types of compost. J Biosci Bioeng. 2015. doi:10.1016/j.jbiosc.2015.05.005.Google Scholar
- Partanen P, Hultman J, Paulin L, Auvinen P, Romantschuk M. Bacterial diversity at different stages of the composting process. BMC Microbiol. 2010;10:94. doi:10.1186/1471-2180-10-94.View ArticlePubMedPubMed CentralGoogle Scholar
- Adekambi T, Stein A, Carvajal J, Raoult D, Drancourt M. Description of Mycobacterium conceptionense sp. nov., a Mycobacterium fortuitum group organism isolated from a posttraumatic osteitis inflammation. J Clin Microbiol. 2006;44(4):1268–73. doi:10.1128/JCM.44.4.1268-1273.2006.View ArticlePubMedPubMed CentralGoogle Scholar
- Brown-Elliott BA, Wallace Jr RJ. Clinical and taxonomic status of pathogenic nonpigmented or late-pigmenting rapidly growing mycobacteria. Clin Microbiol Rev. 2002;15(4):716–46.View ArticlePubMedPubMed CentralGoogle Scholar
- Lamy B, Marchandin H, Hamitouche K, Laurent F. Mycobacterium setense sp. nov., a Mycobacterium fortuitum-group organism isolated from a patient with soft tissue infection and osteitis. Int J Syst Evol Microbiol. 2008;58(Pt 2):486–90. doi:10.1099/ijs.0.65222-0.View ArticlePubMedGoogle Scholar
- Schinsky MF, Morey RE, Steigerwalt AG, Douglas MP, Wilson RW, Floyd MM, et al. Taxonomic variation in the Mycobacterium fortuitum third biovariant complex: description of Mycobacterium boenickei sp. nov. Mycobacterium houstonense sp. nov. Mycobacterium neworleansense sp. nov. and Mycobacterium brisbanense sp. nov. and recognition of Mycobacterium porcinum from human clinical isolates. Int J Syst Evol Microbiol. 2004;54(Pt 5):1653–67.Google Scholar
- Sahraoui N, Ballif M, Zelleg S, Yousfi N, Ritter C, Friedel U, et al. Mycobacterium algericum sp. nov., a novel rapidly growing species related to the Mycobacterium terrae complex and associated with goat lung lesions. Int J Syst Evol Microbiol. 2011;61(Pt 8):1870–4. doi:10.1099/ijs.0.024851-0.View ArticlePubMedGoogle Scholar
- Tortoli E, Gitti Z, Klenk HP, Lauria S, Mannino R, Mantegani P, et al. Survey of 150 strains belonging to the Mycobacterium terrae complex and description of Mycobacterium engbaekii sp. nov., Mycobacterium heraklionense sp. nov. and Mycobacterium longobardum sp. nov. Int J Syst Evol Microbiol. 2013;63(Pt 2):401–11. doi:10.1099/ijs.0.038737-0.View ArticlePubMedGoogle Scholar
- Tortoli E, Baruzzo S, Heijdra Y, Klenk HP, Lauria S, Mariottini A, et al. Mycobacterium insubricum sp. nov. Int J Syst Evol Microbiol. 2009;59(Pt 6):1518–23. doi:10.1099/ijs.0.003459-0.View ArticlePubMedGoogle Scholar
- Slany M, Mrlik V, Kriz P, Pavlik I. First isolation of a newly described Mycobacterium insubricum from freshwater fish. Vet Microbiol. 2010;144(1–2):254–5. doi:10.1016/j.vetmic.2009.12.015.View ArticlePubMedGoogle Scholar
- Waterbury JB, Valois FW. Resistance to co-occurring phages enables marine synechococcus communities to coexist with cyanophages abundant in seawater. Appl Environ Microbiol. 1993;59(10):3393–9.PubMedPubMed CentralGoogle Scholar
- Staley JT, Konopka A. Measurement of in situ activities of nonphotosynthetic microorganisms in aquatic and terrestrial habitats. Annu Rev Microbiol. 1985;39:321–46. doi:10.1146/annurev.mi.39.100185.001541.View ArticlePubMedGoogle Scholar
- Jones Jr WD, Greenberg J. Use of phage F-phi WJ-1 of Mycobacterium fortuitum to discern more phage types of Mycobacterium tuberculosis. J Clin Microbiol. 1976;3(3):324–6.PubMedPubMed CentralGoogle Scholar
- Hatfull GF, Jacobs-Sera D, Lawrence JG, Pope WH, Russell DA, Ko CC, et al. Comparative genomic analysis of 60 Mycobacteriophage genomes: genome clustering, gene acquisition, and gene size. J Mol Biol. 2010;397(1):119–43. doi:10.1016/j.jmb.2010.01.011.View ArticlePubMedPubMed CentralGoogle Scholar
- Pope WH, Jacobs-Sera D, Russell DA, Rubin DH, Kajee A, Msibi ZN, et al. Genomics and proteomics of mycobacteriophage Patience, an accidental tourist in the Mycobacterium neighborhood. mBio. 2014;5(6):e02145. doi:10.1128/mBio.02145-14.View ArticlePubMedPubMed CentralGoogle Scholar
- Hatfull GF, Sarkis GJ. DNA sequence, structure and gene expression of mycobacteriophage L5: a phage system for mycobacterial genetics. Mol Microbiol. 1993;7(3):395–405.View ArticlePubMedGoogle Scholar
- Mediavilla J, Jain S, Kriakov J, Ford ME, Duda RL, Jacobs Jr WR, et al. Genome organization and characterization of mycobacteriophage Bxb1. Mol Microbiol. 2000;38(5):955–70.View ArticlePubMedGoogle Scholar
- Jain S, Hatfull GF. Transcriptional regulation and immunity in mycobacteriophage Bxb1. Mol Microbiol. 2000;38(5):971–85.View ArticlePubMedGoogle Scholar
- Ghosh P, Wasil LR, Hatfull GF. Control of Phage Bxb1 Excision by a Novel Recombination Directionality Factor. PLoS Biol. 2006;4(6):e186.View ArticlePubMedPubMed CentralGoogle Scholar
- Pedulla ML, Ford ME, Houtz JM, Karthikeyan T, Wadsworth C, Lewis JA, et al. Origins of highly mosaic mycobacteriophage genomes. Cell. 2003;113(2):171–82.View ArticlePubMedGoogle Scholar
- Gil F, Grzegorzewicz AE, Catalao MJ, Vital J, McNeil MR, Pimentel M. Mycobacteriophage Ms6 LysB specifically targets the outer membrane of Mycobacterium smegmatis. Microbiology. 2010;156(Pt 5):1497–504. doi:10.1099/mic.0.032821-0.View ArticlePubMedPubMed CentralGoogle Scholar
- Henry M, Coffey A, O' Mahony J, Sleator RD. Comparative modelling of LysB from the mycobacterial bacteriophage Ardmore. Bioengineered bugs. 2011;2(2):88–95. doi:10.4161/bbug.2.2.14138.View ArticlePubMedGoogle Scholar
- Payne K, Sun Q, Sacchettini J, Hatfull GF. Mycobacteriophage Lysin B is a novel mycolylarabinogalactan esterase. Mol Microbiol. 2009;73(3):367–81. doi:10.1111/j.1365-2958.2009.06775.x.View ArticlePubMedPubMed CentralGoogle Scholar
- Bitencourt AL, Vallim MA, Maia D, Spinelli R, Angeloni R, Principal L, et al. Core sampling test in large-scale compost cells for microorganism isolation. Afr J Microbiol Res. 2010;4(15):1631–4.Google Scholar
- Parashar D, Chauhan DS, Sharma VD, Chauhan A, Chauhan SV, Katoch VM. Optimization of procedures for isolation of mycobacteria from soil and water samples obtained in northern India. Appl Environ Microbiol. 2004;70(6):3751–3. doi:10.1128/AEM.70.6.3751-3753.2004.View ArticlePubMedPubMed CentralGoogle Scholar
- Radomski N, Cambau E, Moulin L, Haenn S, Moilleron R, Lucas FS. Comparison of culture methods for isolation of nontuberculous mycobacteria from surface waters. Appl Environ Microbiol. 2010;76(11):3514–20. doi:10.1128/AEM.02659-09.View ArticlePubMedPubMed CentralGoogle Scholar
- Chimara E, Ferrazoli L, Ueky SY, Martins MC, Durham AM, Arbeit RD, et al. Reliable identification of mycobacterial species by PCR-restriction enzyme analysis (PRA)-hsp65 in a reference laboratory and elaboration of a sequence-based extended algorithm of PRA-hsp65 patterns. BMC Microbiol. 2008;8:48.View ArticlePubMedPubMed CentralGoogle Scholar
- Adékambi T, Colson P, Drancourt M. rpoB-based identification of nonpigmented and late-pigmenting rapidly growing mycobacteria. J Clin Microbiol. 2003;41(12):5699–708.View ArticlePubMedPubMed CentralGoogle Scholar
- Selvaraju SB, Khan IU, Yadav JS. A new method for species identification and differentiation of Mycobacterium chelonae complex based on amplified hsp65 restriction analysis (AHSPRA). Mol Cell Probes. 2005;19(2):93–9.View ArticlePubMedGoogle Scholar
- Marchler-Bauer A, Lu S, Anderson JB, Chitsaz F, Derbyshire MK, DeWeese-Scott C, et al. CDD: a Conserved Domain Database for the functional annotation of proteins. Nucleic Acids Res. 2011;39(Database issue):D225–9. doi:10.1093/nar/gkq1189.View ArticlePubMedGoogle Scholar
- Quevillon E, Silventoinen V, Pillai S, Harte N, Mulder N, Apweiler R, et al. InterProScan: protein domains identifier. Nucleic Acids Res. 2005;33(Web Server issue):W116–20. doi:10.1093/nar/gki442.View ArticlePubMedPubMed CentralGoogle Scholar
- Cresawn SG, Bogel M, Day N, Jacobs-Sera D, Hendrix RW, Hatfull GF. Phamerator: a bioinformatic tool for comparative bacteriophage genomics. BMC bioinformatics. 2011;12:395. doi:10.1186/1471-2105-12-395.View ArticlePubMedPubMed CentralGoogle Scholar
- Huson DH. SplitsTree: analyzing and visualizing evolutionary data. Bioinformatics. 1998;14(1):68–73.View ArticlePubMedGoogle Scholar