Mycobacterium avium complex (MAC) causes a broad range of opportunistic infections in humans and animals  . A member of the MAC known as the Mycobacterium avium subspecies hominissuis (hereafter M. avium) is responsible for detrimental health conditions in certain human populations, such as the elderly, children, and AIDS patients   . M. avium is a highly prevalent environmental mycobacterium, with a widespread distribution in water and soil   . Using DNA fingerprinting, M. avium isolates from AIDS patients were detected in residential water, indicating that M. avium infections might be acquired from contaminated water sources  .
Bacteriophages (phages) infect bacterial hosts and are estimated to be the most abundant forms of life on the planet Earth  . While lytic phages lyse the host cells to release progeny phages, temperate phages enter a quiescent state upon viral DNA integration into the host chromosome and reproduce as prophages  . According to the selfish-gene concept, host-incorporated prophage genes are maintained because they contribute to the fitness of the host  . These genes encode for proteins capable of performing key functions, such as metabolism, adhesion, colonization, invasion, spreading, resistance to immune responses, antibiotic resistance, exotoxin production, and serum resistance   . The roles of prophages also include increasing survival or fitness of the host by introducing new fitness factors (lysogenic conversion and transduction), genome rearrangements, gene disruption, protection from lytic infection, and destruction of competitor strains  .
Prophages associated with host virulence have been described in a number of human pathogens, such as Vibrio cholerae, Staphylococcus aureus, Streptococcus sp. and others     . In fact, the public database has recorded large number of prophages identified in the various bacterial genomes   and indicates that phages represent a major driving force in the emergence and evolution of pathogenic bacteria through horizontal transfer of genes. Common examples of prophages-encoding virulence factors are diphtheria toxin, Shiga toxin, cholera toxin, type III secretion system effectors such as Salmonella-derived Sop E2.
Attempting to identify prophage sequences in the M. avium 104 genome, a genome- wide search of the National Center for Biotechnology Information (NCBI) database was performed. Three distinctive prophage genes MAV_0696 (hypothetical protein), MAV_2265 (putative prophage regulatory protein), and MAV_3971 (death on curing protein) were identified. Two other large regions, encompassing approximately 80 genes, have been more recently identified, but not studied in this report. The functions of these genes have yet to be defined. Therefore, biological experiments were designed to characterize M. avium prophage genes and to gain insight into their importance in bacterium pathogenesis.
2. Materials and Methods
2.1. Bacterial Strains and Growth Media
M. avium 104 is a virulent strain isolated from the blood of a patient with AIDS. Mycobacterium smegmatis mc2155 was a gift from Dr. William Jacobs Jr. (Albert Einstein College of Medicine of Yeshiva University, Bronx, NY). Both mycobacteria were grown in Middlebrook 7H9 broth (Difco Laboratories, Detroit, MI), enriched with 10% oleic acid, albumin, dextrose, and catalase (OADC, Difco) or plated on 7H10 agar (Difco) with OADC at 37˚C (pH 7.2). M. smegmatis transformant clones were plated onto 7H10 agar containing 50 mg/ml of kanamycin (Km). E. coli strain DH5B (Stratagene, La Jolla CA) was used as the host for plasmid constructions. Luria-Bertani (LB) broth or LB agar with 50 mg of Kanamycin/ml was employed for growing Escherichia coli transformants. To prepare M. avium 104 and M. smegmatis inocula, bacteria grown on 7H10 agar were resuspended in Hank’s buffered salt solution (HBSS) and adjusted to McFarland Standard #2 turbidity, which corresponds to a suspension of 3.0 ´ 108 colony forming unit (CFU) per ml. In addition, the inocula were serially diluted and plated to determine the number of CFU.
Identification of prophages in bacterial genome can be difficult because of the integration in the genome. Prophages harbor terminases, portal protein, head maturation protein, coat protein or still Tail tape protein. We used NCBI Gene database and PHAST  as the search tool to identify M. avium prophage sequences. M. avium 104 prophage genes (MAV_0696, MAV_2265, and MAV_3971) and encoding proteins were analyzed using the NCBI database (www.ncbi.nlm.nih.gov). The sequences were then matched again the data base.
2.3. Plasmids and PCR Settings
The promotorless mycobacterium shuttle vector pFJS3 (Table 1) was used for making the prophage gene constructs; it was propagated and purified from E. coli using the QIAprep Miniprep kit (QIAGEN, Valencia CA). M. avium strain 104 was used as a source of genomic DNA for prophage gene amplification by PCR. Three unique PCR primer sets with HindIII restriction site were designed to amplify the prophage genes with additional 150 bp sequence upstream of the start codon to include its native ribosomal binding site and promoter (Table 2). PCR amplification was performed in a 30 ml volume vial containing 13.5 ml of sterile H2O, 1 ml DMAO, 15 ml Fidelitaq mix (2X), 0.25 ml prophage-specific forward and reverse primers (100 mM). The PCR parameters were set as follows: 95˚C for 5 min for the initial denaturing step, followed by 35 cycles of 95˚C for 30 s, annealing at 55˚C for 30 s and 72˚C for 2 min, with a final extension at 72˚C for an additional 5 min, and then placed at 4˚C. Because MAC_3971 was not shown to have a role in any of the phenotype evaluated, we decided to use a couple va-
Table 1. Strains used in this study.
Table 2. The PCR primer sets for M. avium prophages containing genes.
riants of the promoter sequence, extended by 25 and 50 bp.
2.4. Construction of the Prophage Gene Clones
After PCR amplification of prophage gene sequences, the products were verified in 1% agarose gel electrophoresis, and sequenced at the Center for Gene Research and Biotechnology (CGRB), Oregon State University, Corvallis. The pFJS3 plasmid (a derivative of MV261 with the LS constitutive promoter downstream of the cloning site) was digested with HindIII restriction enzyme and CIP-treated for 1 h at 37˚C (Table 3). The plasmid DNA and digested prophage inserts were processed for ligation, along with a positive control (digested and CIP-treated pFJS3 without any insert). The resulted positive colonies were stored in 50% glycerol at −80˚C.
2.5. Preparation of M. smegmatis Competent Cells and Transformation
M. smegmatis mc2155 was obtained from glycerol stock and plated onto 7H10 agar. Bacteria were incubated at 37˚C for 5 days, and then transferred into 7H9 broth enriched with OADC until OD600 readings of 1.0 - 2.0 A. Electro-competent cells of M. smegmatis were prepared as the following: the bacterial pellet was washed three times with ice-chilled sterile 10% glycerol and 0.1% Tween 80 buffer. Each time, the pellet was vortex agitated in 10 ml solution before bringing the volume to 30 ml, and centrifuged at 3500 rpm for 20 min at 4˚C. The pellet obtained after the final wash was resuspended in 1 ml of the sterile 10% glycerol and placed on ice for immediate use and optimum transformation efficiency. The electroporation cuvettes were chilled at −20˚C for 1 h. To each of the 200 ml of competent cells in the Eppendorf tubes, duplicates of 5 ml of purified pFJS3 or 5 ml of each transformed plasmid DNA (pFJS3 + MAV_0696, pFJS3 + MAV_2265, pFJS3 + MAV_3971) were added. These contents were then
Table 3. Putative prophage genes found in the genome of M. avium 104.
transferred to chilled cuvettes and electroporated using the parameters set on the Bio- Rad GenePulser under the following conditions: capacitance 25 µF, resistance 1000 W, voltage 2.5 kv. Electroporated cells were immediately recovered by gently mixing 300 ml of sterile 7H9 broth and transfered into sterile tubes for 2 h incubation at 37˚C, with shaking. Each transformation was then plated onto 7H10 agar plates containing 50 mg of kanamycin/ml, and incubated for an additional 4 - 5 days. M. smegmatis positive clones for the three different genes were selected by PCR screening, employing the M. avium specific prophage primers.
To verify if the prophages genes were being expressed, five selected clones for each gene, were prepared and run in polyacrylamide gel. The observed overexpressed protein bands, using coomassie blue stain, were cut out of the gel, eluted and sent to mass spectrometry analysis at OSU Mass Spectrometry Facility, as previously reported (18). Clones in which the over expressed band correspond to the gene cloned were used in the subsequent assays.
2.6. Biofilm Assay
Biofilm formation was determined for M. smegmatis prophage clones (3 clones for each of the genes), and compared to M. smegmatis wild-type with or without pFJS3 vector and M. avium 104. Two hundred ml of 108 bacteria from each clone were inoculated into polyvinyl chloride (PVC) 96-well microplates (Becton-Dickinson Labware, Franklin Lakes NJ) at room temperature in an undisturbed drawer for 3, 5, 7, 11, and 14 days. To measure the biofilm formation, the supernatants were gently removed from each well by inverting the plate onto absorbent pads. Each well was washed with 200 ml of 1´ HBSS as reported  . The biofilms were fixed with 200 ml of methanol for 15 min with the plate lid on, followed by removal of the methanol, and incubation without the lid for an additional 15 min. Crystal violet dye can only stain the bacterial cells and not the PVC material. Fifty ml of a 1% crystal violet solution was added to each well, and the plates were incubated at room temperature for an additional 15 min. Then wells were rinsed three times with 200 ml of 1´ HBSS, and 200 ml of 95% ethanol was used to dissolve the crystal violet. Biofilm formation was analyzed at 570 nm using a microplate reader (Bio-Rad, Hercules, CA).
2.7. Effect of Temperature
To examine the temperature effects on bacteria growth, 100 m of 108 M. smegmatis wild-type, prophage gene clones, and wild-type M. avium 104 were seeded in PCR tubes and incubated for 3 h at different temperatures (20˚C, 30˚C, 50˚C) using a PCR temperature gradient machine. Bacterial colony counts were determined through serial dilution and plating onto 7H10 agar plates.
2.8. Acid Tolerance
Approximately 3.0 ´ 108 M. smegmatis containing the prophage gene clones and wild- type M. avium 104 were inoculated in HBSS at pH 2.0, pH 5.0, and pH 7.0. HBSS was adjusted to pH 2.0 and pH 5.0 with 5 M HCl and 1 M HCl, respectively, and the pH was determined using pH strips. The number of bacteria in the suspension was calculated by plating the suspension on agar plates after 2 h exposure to different pH. Acid inactivation prior to plating was carried out with NaOH.
2.9. Effect of UV Light
One hundred ml of 1´ HBSS, containing approximately 3.0 ´ 108 bacteria, was inoculated into PVC plastic 96-well microtiter plates and exposed to UV for different time intervals (5, 15, 30, 60 min) in the biological safety cabinet. The intensity of the UV lamp (1.d2; 40 microwatts/cm2) was emitted at the wavelength of 253.7 nm in the center of the work surface of the cabinet. Bacteria were placed at 25 cm from the bulb. Viable colonies from plate counts after UV exposure were compared to initial concentrations.
2.10. Prophage Gene Expression
To determine the expression of prophage genes upon biofilm formation, bacterial RNAs were obtained at 24h, days 3 and 5, after M. avium was placed in contact with a polyvinyl chloride plate surface. The Real-Time (RT) PCR was carried out using the conditions as previously described  . Briefly, total bacterial RNAs from broth grown bacteria (control) and from biofilms (experimental) were extracted with the combination of a guanidine thiocyanate-based buffer (Trisol) (Invitrogen, Carlsbad, CA) and rapid mechanical cell lysis of M.avium in a bead-beater. Prior to the real-time PCR, RNA was cleaned up with RNA clean kit (QIAGEN, Valencia, CA) and treated with DNase I. RNA quality was verified by ethidium bromide staining on the agarose gel and by OD260/280 nm absorption. Mycobacterial total RNA (1 mg) was reverse transcribed with 100U of Superscript II Plus RNase H− Reverse Transcriptase (Invitrogen, Carlsbad, CA), using RT primers according to the manufacturer’s instruction. M. avium gene expressions were quantified with SYBR Green I assay by Real-Time PCR detection system using gene-specific primers.
2.11. Macrophage Killing Assay
It was carried out as previously reported  . Briefly, THP-1 macrophages were infected with M. smegmatis controls and M. smegmatis expressing the three prophage genes for 1 h. Then, the extracellular bacteria were removed by washing as described. Intracellular bacteria were quantified at 1 h, 2 days and 3 days after infection. Monolayers were lysed in presence of sterile water and 0.05% SDS for 10 min and the lysate was serially diluted and plated onto 7H10 agar plates. The number of bacteria were determined after 4 days.
2.12. Statistical Analysis
Data represent the means ± standard deviations from three independent experiments. The results of experimental groups and controls were compared using the Student’s t-test (two groups) more than two groups were confirmed by using one-way analysis of variance or ANOVA accordingly. A p < 0.05 was considered significant.
3.1. M. avium Prophage Gene Analysis
Figure 1 shows the genomic sequence of the prophage containing genes. G-C content, conserved domains, and putative function identified in each prophage category of M. avium 104 prophages are summarized in Table 2. Growth rate of selected transformants of M. smegmatis was examined. Three different clones of each of the M. avium genes in M. smegmatis were then incubated in presence of 7H9 broth with OADC and their growth was monitored for up to 7 days. Table 4 shows that all strains had similar growth rates. Colony morphology of the clones was also observed and the strain M. smegmatis with MAV_0696 has a morphotype that is drier than the WT bacterium. The other two clones were indistinguishable from the WT bacterium (Figure 2). We cloned upstream of the gene, 150 bp sequence. All the upstream sequences cloned were examined, by using bioinformatics information, for the presence of a mycobacterial promoter sequence.
3.2. Biofilm Assay
To determine whether the prophages genes would have a role in biofilm formation, a common characteristic of mycobacteria in the environment, we seeded M. smegmatis wild-type and selected three clones overexpressing the prophage proteins in HBSS and compared for differences in the biofilm robustness over time. All three selected clones for each prophage gene showed similar results. We then chose one representative clone for each prophage, for which results are shown. As shown in Table 5, expression of the prophages MAV_2265 and MAV_0696 in M. smegmatis resulted in greater biofilm formation at each indicated time interval when compared to both wild-type and pFJS3 containing M. smegmatis control strains.
3.3. Effects of Temperature on M. smegmatis Survival
Since mycobacteria are exposed to a range of different temperatures in the environment, M. smegmatis expressing prophages genes and wild-type strain were incubated at
Figure 1. Sequences of the prophages MAV_0696, MAV_2265 and MAV_3971.
Figure 2. Morphotypes of M. smegmatis overexpressing the M. avium prophages MAV_0696, MAV_2265 and MAV_3971.
different temperatures (20˚C, 37˚C, 50˚C) to establish whether the prophages genes confer advantage to the bacterium to survive in extreme temperatures. In all three temperatures tested, the presence of prophages genes did not offer any significant advantage to M. smegmatis clones.
Table 4. Growth rate in 7H9 broth of M. smegmatis WT compared to M. smegmatis clones expressing MAC_2266, MAC_0691 and MAC_3971.
Table 5. Biofilm formation by M. smegmatis clones expressing the prophage containing genes.
*p < 0.05 compared with M. smegmatis wild-type and M. smegmatis + pFJS3. 3 different clones expressing the genes were analyzed. All three clones expressing the same gene showed similar results. The mean + SD of one clone is shown.
3.4. Effects of Acidic Environment on M. smegmatis Survival
Because M. avium may encounter acid in both the outside environment and when ingested by the host, we examined the effects of an acidic pH on bacterial growth (survival). M. smegmatis wild-type and prophage genes expressing clones from a five-day-old culture were resuspended in HBSS at acid pH 2.0, pH 5.0 and pH 7.0. Following 2 hours of exposure the acid pH was neutralized and the viability of M. smegmatis was determined by plating onto 7H10 plates. The average CFU is shown in Figure 3. Incubation of experimental and control strains of M. smegmatis in the neutral and highest acidic conditions (pH 2.0) did not result in any significant phenotype characteristic. However, MAV_2265 prophage gene containing M. smegmatis (1.5 × 107) had significantly greater survival rate at pH 5.0 than both the wild-type control (2.9 × 106) and the M. smegmatis containing control pFJS3 plasmid (3.2 × 106). The presence of the other of prophages did not increase in bacterial resistance to acidic conditions (Figure 3).
3.5. Effects of UV Exposure on M. smegmatis Survival
M. avium is subjected to UV exposure from the sunlight in the environment. To evaluate whether prophage genes play a role in UV protection, we exposed the transformants and wild-type strains to UV light at different time intervals and compared the bacterial CFU before and after exposure. It was observed that while 15-min UV exposure resulted in 2-log reduction among transformants clones and wild-type strains of
Figure 3. Viability of M. smegmatis prophage clones at different pHs: pH 2.0, pH 5.0 and pH 7.0. Bacteria were exposed to a range of different pH as described in Material and Methods. At pH 5.0, M. smegmatis overexpressing MAP_2265 was observed to have increased resistance to the environmental conditions. (*) p < 0.05 compared with the other genes.
M. smegmatis, 30-min UV exposure led to over 5-log decrease in bacterial CFUs as shown in Figure 4. However, there was not any significance observed in survival rate of M. smegmatis clones expressing prophage genes compared to both control strains.
3.6. Expression of Prophages under Biofilm Formation
To examine whether prophage genes were upregulated upon biofilm formation, we performed the quantitative Real-Time PCR using M. avium RNA. The expression levels of target genes from biofilms and from broth-grown bacteria were normalized to the expression level of the endogenous reference 16S rRNA in each sample. The M. avium prophage gene MAV_2265 showed the greatest level of expression in biofilm at all selected time points compared with the expression levels of MAV_0696 and MAV_3971 (Figure 5). While MAV_0696 was upregulated over 2-fold at day 3, the induction increased up to 3.5 at day 5. There were not any changes observed in MAV_3971 gene expression levels over time as shown in Figure 5.
3.7. Survival in Macrophages
M. smegmatis expressing M. avium prophages sequences were used to infect THP-1 macrophages. M. smegmatis with the empty plasmid was used as control. None of the prophages sequences had any impact on the ability of M. smegmatis to survive in macrophages. By day 2 after infection, there was a reduction of 36%, 40%, 38%, and 39% for the M. smegmatis with the empty plasmid and expressing MAV_2265, MAV_0696 and MAV_3971 respectively. By day 3 after infection, the bacterial load reduction was 67%, 71%, 74%, and 68% respectively.
Figure 4. Viability of M. smegmatis prophage clones after 5 min, 15 min and 30 min exposure to UV light. UV light is an important environmental factor. The fact that overexpression of prophages had no effect on the clones to resist to UV light suggest that the phages are not associated with environmental protection.
Figure 5. The Real Time PCR quantification of M. avium genes in biofilm. Bacterial RNA was purified, as described in materials and methods. Total RNA was used to determine the copy number of cDNAs for target prophage and reference 16S genes. Data were analyzed on the basis of Ct values of each sample and normalized with an internal housekeeping gene control, 16S rRNA. Values shown are representative of three different experiments with very similar results.
Both humans and animals are susceptible to M. avium infection. This pathogen has the ability to form biofilm or a biofilm-like structure, which has been associated with chronic bacterial infection   . Study of M. avium infection in the lung has suggested a possible association between biofilm formation and the difficulty in responding to therapy   . In the experiments reported here, we attempted to determine whether chromosomal prophage genes had any role in the adaptation of M. avium to the environment. We screened each prophage gene transformant in M. smegmatis (gain of function assays), and compared the results with wild-type strains of M. avium 104 and M. smegmatis. The biofilm assay revealed that both MAV_2265 and MAV_ 0696 transformed into M. smegmatis resulted in a significant increase in biofilm formation at each indicated time interval, when compared to other strains used in this study. The degree of biofilm formation induced in M. smegmatis by the phage was even greater than the ability of M. avium 104 to form biofilm, which suggested that the prophage gene may act in synergism with other biofilm related genes. In addition, it is more likely that either biofilm or biofilm-like structure establishment is more regulated in M. avium 104 than in M. smegmatis, since wild-type M. smegmatis and all M. smegmatis transformants form more biofilm than M. avium.
A number of different prophages genes have been shown to be associated with virulence in bacteria such as Staphylococcus aureus   , Streptococcus sp.    and Vibrio cholerae   , Shigella sp.  , Pseudomonas aeruginosa  , Salmonella enterica  , and others. In S. aureus for instance, prophages are expressed during animal infection and the absence of these phage genes results in virulence defects in a murine model of abscess formation  , thus revealing essential contributions of prophages to the pathogenesis. Cholera toxin (ctx AB) is encoded by genes present in the prophage ctx0  and is not produced in absence of the prophage. Prophages in Streptococcus and Staphylococcus are usually observed in “prophage regions” in the genome, in contrast to prophages in low CG content bacteria, such as gram-ne- gative bacteria  .
In the case of M. avium, we hypothesized that prophages would probably be linked to the ability of bacteria to interact with the environment. In fact, MAV_2265 and MAV_0696 were the prophage genes associated with the ability to form biofilm, and M. smegmatis MAV_2265 prophage clone had significant survival rate at pH 5.0 compared with control wild-type and other prophage clones. The MAV_2265 gene was upregulated in M. avium under conditions of biofilm formation by 6.8-fold at day 3. The expression of MAV_0696 resulted in a 3.5-fold increase in M. avium biofilms at day 5, while the evaluated MAV_3971 prophage gene showed no changes in expression over time.
The sequences of the three M. avium prophages in this study belong to XRE-family and DOC family. Two out of three sequences suggest regulatory function, which would explain the fact that they are located in isolated regions of the genome surrounded by transposases, phage intergrases and tRNA. Recently discovered prophages in Salmonella are linked to type III secretion systems such as SopE2  , demonstrating that prophage sequences can carry an array of functions associated with crucial phenotypes for pathogenic activity. Also recently, other prophages regions have been discovered in the M. avium genome. Since this finding happened after this paper had been submitted, the regions were not investigated.
In the case of MAV_2265, the discovery that the prophage gene is associated with an increase of biofilm formation raises important questions. The prophage GC content suggests acquisition of these genes from a different bacterial species, perhaps another environmental mycobacteria. The mechanism(s) involved of MAV_2265 and MAV_ 0696 participation in biofilm formation is still unknown at this point. Interestingly, the other prophage gene MAV_3971 has no significant role in any of the characteristics tested (acid resistance, temperature tolerance, biofilm formation, resistance to UV light), indicating that its function may be related with other aspects in the bacterial physiology. MAV_0696 is flanked by MAV_0697, which encodes for a protein containing a DNA metabolism/replication domain and may have a function associated with DNA replication under different stresses.
We have recently finished the sequence of 5 different M. avium-M. intracellulare complex strains (MAC 100, MAC 101, MAC A5, MAC 3388, MAC 3387). Search in the genome of the strains demonstrates that MAC 104 and MAC 101, two M. avium subsp. hominissuis, contain the MAV_0696 and MAV_2264 genes when their sequences are not present in the other strains. MAV_3971 in contrast is present in all strains. None of the sequences is present in Mycobacterium tuberculosis, Mycobacterium marinum, Mycobacterium leprae, which indicates that M. avium acquires them in a unique environment.
In summary, biofilm development has been shown to have an important impact on resistance to a number of stresses, including antibiotic exposure, while also having a significant impact on growth. Here we have identified three unique prophage sequences in the M. avium 104 genome and established that MAV_2265 and MAV_0697 genes participate in biofilm formation and low pH adaptation.
We thank Sandra Sommer, Michael McNamara, Sasha Rose, and Laura Meek for insightful discussion and suggestions; and Denny Weber and Beth Chamblin for preparing the manuscript.
 von Reyn, C.F., Maslow, J.N., Barber, T.W., Falkinham, J.O. and Arbeit, R.D. (1994) Persistent Colonisation of Potable Water as a Source of Mycobacterium avium Infection in AIDS. Lancet, 343, 1137-1141.
 Johansen, T.B., Olsen, I., Jensen, M.R., Dahle, U.R., Holstad, G. and Djonne, B. (2007) New Probes Used for IS1245 and IS1311 Restriction Fragment Length Polymorphism of Mycobacterium avium subsp. avium and Mycobacterium avium subsp. hominissuis Isolates of Human and Animal Origin in Norway. BMC Microbiology, 7, 14.
 Canchaya, C., Fournous, G., Chibani-Chennoufi, S., Dillmann, M.L. and Brussow, H. (2003) Phage as Agents of Lateral Gene Transfer. Current Opinion in Microbiology, 6, 417-424.
 Canchaya, C., Proux, C., Fournous, G., Bruttin, A. and Brussow, H. (2003) Prophage Genomics. Microbiology and Molecular Biology Reviews, 67, 238-276.
 Schuch, R. and Fischetti, V.A. (2006) Detailed Genomic Analysis of the W Beta and Gamma Phages Infecting Bacillus anthracis: Implications for Evolution of Environmental Fitness and Antibiotic Resistance. Journal of Bacteriology, 188, 3037-3051.
 Wagner, P.L., Livny, J., Neely, M.N., Acheson, D.W., Friedman, D.I. and Waldor, M.K. (2002) Bacteriophage Control of Shiga Toxin 1 Production and Release by Escherichia coli. Molecular Microbiology, 44, 957-970.
 Brussow, H., Canchaya, C. and Hardt, W.D. (2004) Phages and the Evolution of Bacterial Pathogens: From Genomic Rearrangements to Lysogenic Conversion. Microbiology and Molecular Biology Reviews, 68, 560-602.
 Bae, T., Baba, T., Hiramatsu, K. and Schneewind, O. (2006) Prophages of Staphylococcus aureus Newman and Their Contribution to Virulence. Molecular Microbiology, 62, 1035-1047.
 Boyd, E.F. and Brussow, H. (2002) Common Themes among Bacteriophage-Encoded Virulence Factors and Diversity among the Bacteriophages Involved. Trends in Microbiology, 10, 521-529.
 Sumby, P., Barbian, K.D., Gardner, D.J., Whitney, A.R., Welty, D.M., Long, R.D., Bailey, J.R., Parnell, M.J., Hoe, N.P., Adams, G.G., et al. (2005) Extracellular Deoxyribonuclease made by Group A Streptococcus Assists Pathogenesis by Enhancing Evasion of the Innate Immune Response. Proceedings of the National Academy of Sciences of the United States of America, 102, 1679-1684.
 Hatfull, G.F., Cresawn, S.G. and Hendrix, R.W. (2008) Comparative Genomics of the Mycobacteriophages: Insights into Bacteriophage Evolution. Research in Microbiology, 159, 332-339.
 Carter, G., Wu, M., Drummond, D.C. and Bermudez, L.E. (2003) Characterization of Biofilm Formation by Clinical Isolates of Mycobacterium avium. Journal of Medical Microbiology, 52, 747-752.
 Yamazaki, Y., Danelishvili, L., Wu, M., MacNab, M. and Bermudez, L.E. (2006) Mycobacterium avium Genes Associated with the Ability to Form a Biofilm. Applied and Environmental Microbiology, 72, 819-825.
 Danelishvili, L., Wu, M., Stang, B., Harriff, M., Cirillo, S.L., Cirillo, J.D., Bildfell, R., Arbogast, B. and Bermudez, L.E. (2007) Identification of Mycobacterium avium Pathogenicity Island Important for Macrophage and Amoeba Infection. Proceedings of the National Academy of Sciences of the United States of America, 104, 11038-11043.
 Olivier, K.N., Weber, D.J., Wallace Jr., R.J., Faiz, A.R., Lee, J.H., Zhang, Y., Brown-Elliot, B.A., Handler, A., Wilson, R.W., Schechter, M.S., et al. (2003) Nontuberculous Mycobacteria. I: Multicenter Prevalence Study in Cystic Fibrosis. American Journal of Respiratory and Critical Care Medicine, 167, 828-834.
 Yamazaki, Y., Danelishvili, L., Wu, M., Hidaka, E., Katsuyama, T., Stang, B., Petrofsky, M., Bildfell, R. and Bermudez, L.E. (2006) The Ability to form Biofilm Influences Mycobacterium avium Invasion and Translocation of Bronchial Epithelial Cells. Cellular Microbiology, 8, 806-814.
 Holden, M.T., Feil, E.J., Lindsay, J.A., Peacock, S.J., Day, N.P., Enright, M.C., Foster, T.J., Moore, C.E., Hurst, L., Atkin, R., et al. (2004) Complete Genomes of Two Clinical Staphylococcus aureus Strains: Evidence for the Rapid Evolution of Virulence and Drug Resistance. Proceedings of the National Academy of Sciences of the United States of America, 101, 9786-9791.
 Aziz, R.K., Edwards, R.A., Taylor, W.W., Low, D.E., McGeer, A. and Kotb, M. (2005) Mosaic Prophages with Horizontally Acquired Genes Account for the Emergence and Diversification of the Globally Disseminated M1T1 Clone of Streptococcus pyogenes. Journal of Bacteriology, 187, 3311-3318.
 Banks, D.J., Lei, B. and Musser, J.M. (2003) Prophage Induction and Expression of Prophage-Encoded Virulence Factors in Group A Streptococcus Serotype M3 Strain MGAS315. Infection and Immunity, 71, 7079-7086.
 Reidl, J. and Mekalanos, J.J. (1995) Characterization of Vibrio cholerae Bacteriophage K139 and Use of a Novel Mini-Transposon to Identify a Phage-Encoded Virulence Factor. Molecular Microbiology, 18, 685-701.
 Allison, G.E. and Verma, N.K. (2000) Serotype-Converting Bacteriophages and O-Antigen Modification in Shigella flexneri. Trends in Microbiology, 8, 17-23.
 Nakayama, K., Kanaya, S., Ohnishi, M., Terawaki, Y. and Hayashi, T. (1999) The Complete Nucleotide Sequence of phi CTX, a Cytotoxin-Converting Phage of Pseudomonas aeruginosa: Implications for Phage Evolution and Horizontal Gene Transfer via Bacteriophages. Molecular Microbiology, 31, 399-419.
 Figueroa-Bossi, N., Uzzau, S., Maloriol, D. and Bossi, L. (2001) Variable Assortment of Prophages Provides a Transferable Repertoire of Pathogenic Determinants in Salmonella. Molecular Microbiology, 39, 260-271.
 Boyd, E.F. and Brussow, H. (2002) Common Themes among Bacteriophage-Encoded Virulence Factors and Diversity among the Bacteriophages Involved. Trends in Microbiology, 10, 521-529.
 Mirold, S., Rabsch, W., Tschape, H. and Hardt, W.D. (2001) Transfer of the Salmonella Type III Effector sopE between Unrelated Phage Families. Journal of Molecular Biology, 312, 7-16.