- Research article
- Open Access
Whole genome sequencing identifies an allele responsible for clear vs. turbid plaque morphology in a Mycobacteriophage
BMC Microbiology volume 20, Article number: 148 (2020)
Whole genome sequencing promises to revolutionize our ability to link genotypic and phenotypic variation in a wide range of model and non-model species.
Here we describe the isolation and characterization of a novel mycobacteriophage named BGlluviae that grows on Mycobacterium smegmatis mc2155. BGlluviae normally produces turbid plaques but a spontaneous clear plaque was also recovered. The genomic DNA from pure populations of the BGlluviae phage and the clear plaque mutant were sequenced. A single substitution, at amino acid 54 (I to T), in the immunity repressor protein resulted in a clear plaque phenotype.
This substitution is predicted to be located at the subunit interaction interface of the repressor protein, and thus prevents the establishment of lysogeny.
Bacteriophages such as λ and T4 have long served as model systems in genetic research since they are easy to culture, have simple genomes, and vary in numerous aspects of plaque morphology . Lambda phage is the most well studied bacteriophage, especially with respect to life cycle and plaque morphology . Wild-type λ phage is a temperate phage, and therefore has the ability to enter into either lytic or lysogenic life cycles . During the lytic life cycle, the phage utilizes the host’s proteins to transcribe and translate phage genes necessary for replication and the construction of new phages . The new phage genomes are packaged into new phage bodies, which eventually erupt from the cell, killing it, to infect nearby cells . Lytic phages leave clear plaques on a bacterial lawn, since they lyse or kill all bacterial cells they infect. In contrast, during the lysogenic life cycle, the phage genome integrates into the host genome . The phage genome is copied along with that of the host, but no new phages are produced and none leave the cell to infect other nearby cells . Temperate phages like λ leave turbid plaques on a bacterial lawn, since they do not lyse all bacterial cells they infect .
Kaiser  isolated 40 different mutant strains of λ phage that produce clear plaques owing to their inability to lysogenize (i.e. all bacterial cells are lysed). Through a series of complementation crosses, Kaiser mapped the mutations responsible for clear plaques to three tightly-linked genes, later termed cI, cII, and cIII . Later, Ptashne et. al.  isolated the first of these genes, cI, now called the CI repressor, which is responsible for repressing the lytic life cycle in favor of the lysogenic life cycle in bacteriophages. In addition to illuminating the workings of the bacteriophage life cycle, the CI repressor has continued to serve as an important model for how Repressor proteins interact with DNA and repress transcription .
More recently, mycobacteriophages have become important model systems for understanding phage biology and bacterial pathogens, especially since this group contains members that infect the bacterial species Mycobacterium tuberculosis, the pathogen that causes tuberculosis . For example, Donelly-Wu et al.  identified Repressor protein of mycobacteriophage L5 as a selectable marker for M. smegmatis growth but not for M. tuberculosis. Selectable markers can be used as tools to study the pathogenicity of a bacterium. Later, Petrova et al.  mapped the repressor gene of lytic mutants of the temperate M. smegmatis mycobacteriophage Adephagia. In this strain, the repressor gene is a selectable marker for slow and fast growing bacteria, as it confers immunity to Adephagia superinfection . The identification of this gene is important because it can be useful to manipulate the mycobacterial host and provide improved recombinant vaccines against tuberculosis. Specifically, transformants could carry the prophage, which would promote plasmid maintenance . Mycobacteriophages have even become important educational tools for understanding DNA and simple bacteriophage genomes .
In this study, we isolated novel mycobacteriophages specific to the host Mycobacterium smegmatis mc2155. We then used whole genome sequencing to identify the sequence of both a novel clear and turbid plaque producing mycobacteriophage in an effort to identify the mutations responsible for observed differences in plaque morphology.
Phage isolation and characterization
We isolated a mycobacteriophage that produced two plaque morphologies: turbid and clear plaques. Both plaque types are circular, approximately 2.5 mm in diameter with well-defined borders (Fig. 1a, b). We named the mycobacteriophage BGlluviae after the location in which it was discovered (Bowling Green, KY, from a drainage ditch; lluvia is the Spanish word for rain). Transmission electron microscopy of phages isolated from the turbid plaques revealed that the BGlluviae phage exhibits the siphoviral morphology that is typical of Mycobacteriophages, including an isometric head and long, flexible tail . In the case of BGlluviae, the capsid head is only 43 nm in diameter, with a length to width ratio of 0.887. This capsid size is relatively small, accounting for just under 25% of the total length of the phage (187 nm) (Fig. 1c). For comparison, mycobacteriophages Che8 and Fruitloop (Cluster F1) have similar genome sizes but capsid heads of 59 and 60 nm in diameter, respectively .
Genome sequencing and assembly
The assembly of short 100–200 bp next-generation sequencing reads for both BGlluviae and its clear plaque derivative produced three contigs containing approximately 98,792 reads. In both cases, the longest contig produced final genome assemblies that were both 59,308 bp in length, with sequencing coverage of at least 120-fold depth. A BLAST search with default options revealed that the BGlluviae turbid plaque genome is 99% identical to the mycobacteriophage Belladonna (score = 1.043e05, Identities = 52,745/52,793, Gaps = 2/52,793, E-value = 0.0,), differing primarily in the last 7000 bp where one gene of unknown function is absent in the BGlluviae genome. This result indicates that, like Belladonna, BGlluviae is a member of the K1 cluster of mycobacteriophages. K1 phage genomes are characterized by the 11 bp 3′ overhang 5′-CTCGTAGGCAT-3′. Automated and manual annotation of the BGlluviae genomes isolated from both turbid and clear plaques revealed the presence of 95 protein coding genes and 1 tRNA, with an overall GC content of 66.5%. The DNA sequence and annotation of the BGlluviae turbid plaque strain is published as GenBank accession number MN908692.
The genomes of the BGlluviae phage and the clear plaque mutant were aligned across their entire length using LAGAN . This alignment revealed the presence of just a single difference between the two genomes: a single nucleotide polymorphism (SNP) at position 34,936. This SNP, an T → C transition, is unlikely to be the result of sequencing error. The T genotype of the turbid phage strain is supported by 272/276 (98.6%) of sequencing reads at this position, while the C genotype of the clear phage strain is supported by 278/278 (100%) of reads. This SNP at position 34,936 is located at amino acid position 54 of BGlluviae Gene_43, where it produces a missense mutation that changes Ile to Thr.
Sequence homology of BGlluviae Gene_43
Manual annotation of BGlluviae Gene_43 using sequence and structural similarity to numerous extrinsic databases, including PhagesDB, BLASTP, and HHPRED, identified BGlluviae Gene_43 as a potential Immunity Repressor protein. For example, 50 out of 81 (62%) significant local alignments of BGlluviae Gene_43 to other mycobacteriophage proteins in PhagesDB are annotated as putative Immunity Repressors. The same is true of 36 out of the top 50 (72%) BLASTP hits for BGlluviae Gene_43 in the NCBI non-redundant protein database (Query coverage range from 65 to 100%; Percent identity range from 51.16 to 99.21%; E-value range from 4e-26 to 1e-82). In addition, HHPRED, which attempts to assign protein homology based on predicted secondary structure, revealed a significant structural alignment to the Repressor protein of the model bacteriophage Enterobacteria phage Lambda (PDB ID: 1LLI_B; Probability = 98.4; Query coverage = 53.17%; alignment from amino acids 3–70; E-value = 4.6e-8). The sequence similarity of BGlluviae Gene_43 to the proteins of viral genomes outside of the mycobacteriophages, such as Bacteriophage 434 (NCBI Taxid 10,712), Lactococcus phage TP901–1 (NCBI Taxid 35,345), and Salmonella phage P22 (NCBI Taxid 10,754) was generally much lower, although still significant (Query coverage range from 34 to 57%; Percent identity range from 10.77 to 18.06%, E-value range from 2e-10 to 2e-19). This low sequence similarity to more distantly related viral taxa was reflected in a maximum likelihood phylogeny of these proteins, which shows BGlluviae Gene_43 clearly related to other mycobacteriophage immunity repressors, but more distantly related to the repressors of other phages (Fig. 2). However, BGlluviae Gene_43 is more closely related to all other repressors than to the outgroup tail assembly chaperone proteins. Thus, based on its sequence similarity to these annotated proteins, as well as its placement in a monophyletic group with predicted or established Repressor proteins, we identified BGlluviae Gene_43 as a potential Immunity Repressor protein.
Structural modeling of BGlluviae immunity repressor
To infer the impact of the Ile54Thr substitution on BGlluviae Immunity Repressor function, we modeled the secondary, tertiary, and quaternary structure of this protein using various methods based on sequence similarity and de novo prediction (see Methods). Based on alignment to similar proteins in the Conserved Domain Database , the Ile54Thr substitution site is predicted to occur within a broad cro/C1-type helix-turn-helix domain found in many bacteriophage transcriptional repressors (Domain HTH_19, pfam12844; from amino acids 20–73; E-value = 9.61e-3), but outside the range of the narrower two-member helix-turn-helix motif used to bind DNA (Fig. 3). Instead, the Ile54Thr substitution occurs within helix 4, near the beginning of the interface range between amino acid residues 50–86 used for protein dimerization (Fig. 3, red triangle).
Modeling of the tertiary and quaternary structures of the BGlluviae turbid and clear Immunity Repressor proteins via threading in I-TASSER [13, 14] gives C-scores of − 4.05 and − 4.24, respectively. The two models of the BGlluviae Repressor proteins with the 54Ile and 54Thr substitutions are shown in Fig. 4a. The relative locations of the first five helices are fairly conserved. When these tertiary models are then superimposed onto the N-terminal arms of the lambda Repressor protein (PDB ID 1LMB) to generate a quaternary model of the Repressor dimer interacting with DNA, the two amino acid 54 residues can be seen pointing into the interface region (Fig. 4b). Figure 5 shows the interface from the top and reveals that amino acid 54 points towards helix 6 of the opposite paired subunit. Differences in this interactive interface could impact the ability of the protein to dimerize or alter the distance between the DNA binding domains such that stable DNA binding is not possible. The disruption of repressor dimerization or DNA binding would be consistent with the potential inhibition of lysogeny and the clear plaque phenotype found in 54Thr mutants.
We have isolated a novel mycobacteriophage from Bowling Green, Kentucky specific to the host Mycobacterium smegmatis mc2155. This bacteriophage, BGlluviae, is a member of the K1 cluster of temperate mycobacteriophages that normally produces small, turbid plaques. In the course of purification, a clear plaque mutant was recovered (Fig. 1). We sequenced the genomes of both the turbid and clear plaque forms of BGlluviae and have published the annotated genome in GenBank as accession MN908692. A single nucleotide polymorphism (SNP) was identified within BGlluviae Gene_43, a predicted Immunity Repressor protein (Fig. 3). This SNP, a T → C transition, results in a Ile to Thr missense mutation at position 54, in a helix motif, located at the interface between the subunits (Figs. 4 and 5).
The results of this study are consistent with the role of the Repressor protein in regulating the switch between the lytic and lysogenic life cycle in bacteriophages. In the most well-characterized example, Lambda phage, the Repressor (cI) works in conjunction with cro to determine lysogenic versus lytic life cycles (reviewed in ). During the lysogenic life cycle, the repressor cI is transcribed. Its protein product, CI, dimerizes with itself and cooperatively binds to two operators via a helix-turn-helix DNA binding domain [5, 15]. One of these operators, OR, is found directly upstream of the adjacent cro gene. When the operator is bound by CI, the expression of the cro gene is repressed, preventing the transcription of genes necessary for lysis . With the lytic life cycle repressed, Lambda phage integrates into the host genome. On a bacterial lawn, a small amount of background lysis from these otherwise lysogenic phages produces turbid plaques made up of a majority of living bacterial cells. This process of repression is further reinforced by the products of two additional genes, cII and cIII [17, 18].
In contrast, during the lytic life cycle, cro is transcribed. Its protein product, Cro, also binds to the operator OR, but does so in such a way that transcription of the upstream gene cI is blocked while that of the downstream cro is permitted. When the operator is bound by Cro, repression by CI is itself repressed, permitting the transcription of genes necessary for lysis. With the lysogenic life cycle repressed, the λ phages replicate themselves and express proteins necessary for phage head and tail formation, assembly, and lysis. On a bacterial lawn, the complete lysis of all infected cells produces clear plaques devoid of any living bacteria . Thus, we propose that the Ile to Thr substitution, found in the product of BGlluviae Gene_43, results in clear plaque formation by inhibiting the ability of the BGlluviae Immunity Repressor protein to properly bind between subunits, thus inhibiting its ability to repress the BGlluviae cro-like gene and prevent its transcription, leading to lysis. In support of this hypothesis, BGlluviae Gene_43 and 44, annotated as Repressor and Cro-like proteins, respectively, are located adjacent to each other in reverse and forward orientations similar to the orthologous genes in Lambda phage .
This conclusion, of course, depends on the accuracy of modeling the BGlluviae Gene_43 repressor protein and its comparison to the tertiary structure of the Lambda CI Repressor N-terminal structure (PDB ID 1LMB). Although BGlluviae Gene 43 and the Lambda CI Repressor share only 12% sequence identity, it is not unusual for distantly related homologs to diverge in sequence identity while still retaining overall structural similarity through the use of biochemically similar amino acids . Additionally, a pairwise alignment between these two proteins using DELTA-BLAST reveals that their sequences are most similar between amino acids 19 and 67, while the alignment of the predicted BGlluviae Repressor tertiary structure with lambda shows closely corresponding helices and similar location in 3D space over amino acids 9 and 87 (Fig. 4).
The Immunity Repressor has already been implicated in the formation of different plaque mutants in mycobacteriophages. In particular, Petrova et al.  found four independent protein mutations from seven clear plaques in lysates from the mycobacteriophage Adephagia. Adephagia has three important similarities to BGlluviae: (1) Adephagia is a member of the K1 subcluster of mycobacteriophages, (2) the Adephagia clear plaque mutations were all located within the Immunity Repressor protein, and (3) these mutations impact the cro/C1-type HTH domain (Fig. 2). In Adephagia, one mutation entirely deletes the first and second helices used for DNA binding (not shown), while the remaining three fall within the second, fourth, and fifth helices (Fig. 2). Two of these mutations also occur in the helices necessary for protein dimerization, and one, Arg55Gly, is located just a single amino acid away from the Ile54Thr mutation reported here in BGlluviae (Fig. 2). Interestingly, three of the four mutations found thus far in Adephagia exchange large amino acids (I, R, and Y) for smaller ones (T, G, and C respectively). Importantly, the Adephagia Immunity Repressor protein is almost certainly orthologous to the BGlluviae Immunity Repressor protein identified in this study based on the high sequence identity shown between the Immunity Repressor proteins of the BGlluviae turbid-plaque strain and Adephagia (BLASTP Identity = 100%, Coverage = 100%, E-value = 8 × 10− 93). Thus, the results found here for BGlluviae further support previous observations from the Adephagia mycobacteriophage that mutations within the Immunity Repressor protein contribute to clear plaque formation in mycobacteriophages.
We have isolated and characterized the genome of a mycobacteriophage from Bowling Green, Kentucky, BGlluviae (Genbank Accession Number MN908692), that produces clear and turbid plaques. Through the use of whole genome sequencing, sequence alignment, and protein modeling, we determined that the mutation responsible for clear plaque formation occurs in an important dimer interface helix of the Immunity Repressor protein, a protein that regulates the lytic vs. lysogenic life cycle in bacteriophages. This result confirms the utility of bacteriophages in general – and mycobacteriophages in particular – to model genetic processes, and confirms the utility of whole genome sequencing and protein modeling to quickly link genotype and phenotype.
Soil was collected from a storm drain in Bowling Green, KY (36° 57′ 25″N and 86° 27′ 50″ W) and mycobacteriophages were isolated from this sample using previously described protocols “(https://phagesdb.org/workflow/)”. Briefly, phages that use M. smegmatis mc2155 as their hosts were enriched by combining 1 g of soil with 20 ml water, 2.5 ml Middlebrook 7H9 broth base + glycerol, 2.5 ml albumin + dextrose supplement, 250 μl of 100 mM CaCl2, and 2.5 ml of an 48 h culture of M. smegmatis mc2155. The culture was shaken at 250 rpm for 24 h at 37 °C. After incubation, larger particles (soil and M. smegmatis cells) were removed by centrifugation at 3000 rpm for 10 min. The cleared supernatant was then filtered through a 0.2 μm syringe filter to remove any remaining cells.
Serial dilutions (10− 1 to 10− 5) of the filtrate were mixed with 0.5 ml aliquots of M. smegmatis cells for 30 min to allow for phage attachment. After incubation, 4.5 ml of Middlebrook Top Agar was added, and the entire mixture was transferred onto Luria agar plates supplemented with carbenicillin and cycloheximide. After the plates solidified, they were incubated at 37 °C for 24 h. A pure population of a turbid plaque forming phage was obtained using three rounds of plaque streak-plating. During the purification process, a clear plaque mutant appeared on the plates, and was separately purified using three rounds of plaque streak-plating. This mutant was spontaneous, and was not obtained by induced mutagenesis. Both the turbid plaque phage and the clear plaque phage were then grown to produce high-titer lysates for DNA isolation.
We pelleted a 1.0 ml sample of this high-titer lysate, removed the supernatant, and re-suspended the pellet in 100 μl fresh phage buffer. We then pipetted 10 μl of the sample onto an EM grid, washed the grid twice with sterile water, and stained the grid with 1% uranyl acetate. Bacteriophages were viewed and photographed used a JEOL 1500-Plus transmission electron microscope at Western Kentucky University’s Ogden College Electron Microscopy Facility.
DNA isolation and sequencing
Genomic DNA was isolated from high titer lysates using a Promega DNA Clean Up Kit (Madison, WI) following the protocol available here: https://phagesdb.org/workflow/Extraction/. Sequencing libraries of the turbid and the clear plaque derivative were prepared using the Ion Shear Kit V2 and the Ion Xpress Library Kit. Libraries were processed on an Ion Torrent PGM using a 314 chip for single-end reads of 100–200 bp.
Genome assembly and annotation
Using the FASTX genome analysis toolkit , we assembled 100,000 genomic reads into contigs using the overlap-layout-consensus algorithm implemented in the short read assembler Newbler v2.6  using default parameters. BLASTN was used to confirm that orientation coincided with other mycobacteriophages in the PhagesDB . Phage genes were identified using Glimmer , GeneMark , Aragorn , and tRNAscan-SE , and by manual inspection and annotation revision using DNA Master “(http://cobamide2.bio.pitt.edu)” and PECAAN . The genome sequence and annotation of the turbid plaque morphology phage is available at GenBank (https://www.ncbi.nlm.nih.gov) accession number MN908692. To compare the sequence of BGlluviae phages isolated from turbid and clear plaques, we aligned the entire genome sequence of each using LAGAN .
We identified potential homologs to BGlluviae proteins by using BLASTP and DELTA-BLAST  to search the PhageDB , NCBI RefSeq , and Protein Data Bank  databases for similar protein sequences, restricting all results to viral taxa (NCBI TaxID: 10239). From each search, we chose significant alignments (E < 0.05) representing different viral lineages, then aligned all sequences using the Expresso structural multiple alignment algorithm . Finally, we used maximum likelihood to infer the phylogenetic relationship of these homologs using W-IQ-TREE under all default parameters necessary to infer an appropriate model of sequence evolution and assess bootstrap support .
We characterized the impact of any mutations on protein secondary and tertiary structure using several methods to infer protein structure de novo or with reference to previously-established tertiary structures from homologous proteins using the programs DELTA-BLAST , JPred , PSIPRED [28, 35] and I_TASSER [13, 14]. In addition to local alignment, DELTA-BLAST can also be used to identify potential domains in protein sequences by searching for similar sequences in the Conserved Domain Database. PSIPRED and JPred are programs that predict protein secondary structure, including helices, from primary sequence information using position-specific scoring matrices similar to those used in PSI- and DELTA-BLAST . I- TASSER [13, 14] predicted protein tertiary structure by threading together structural information from alignments to multiple different homologs across multiple parts of a protein’s primary sequence. Finally, we used the program UCSF Chimera  to superimpose these predicted tertiary structures with experimentally-determined structures of homologous proteins from Protein Data Bank.
Availability of data and materials
Laboratory protocols and software documentation are available at phagesdb.org/documents. Details relating to the habitat and geographic location in which the phage were discovered, genome sequencing information, taxonomic and morphological characterization of the phage, plaque photographs, and TEM photographs of the phage are all available at phagesdb.org/phages/BGlluviae. The annotated sequence is available at GenBank (https://www.ncbi.nlm.nih.gov) as accession number MN908692.
Weaver RF. Molecular Biology. 5th ed. New York: McGraw-Hill; 2011.
Kaiser AD. Mutations in a temperate bacteriophage affecting its ability to lysogenize Escherichia coli. Virology. 1957;3(1):42–61.
Jacob F, Monod J. Genetic regulatory mechanisms in the synthesis of proteins. J Mol Biol. 1961;3(3):318–56.
Ptashne M. Isolation of the λ phage repressor. Proc Natl Acad Sci U S A. 1967;57(2):306.
Stayrook S, Jaru-Ampornpan P, Ni J, Hochschild A, Lewis M. Crystal structure of the λ repressor and a model for pairwise cooperative operator binding. Nature. 2008;452(7190):1022.
Hatfull GF. Mycobacteriophages: windows into tuberculosis. PLoS Pathog. 2014;10(3):e1003953.
Donnelly-Wu MK, Jacobs WR Jr, Hatfull GF. Superinfection immunity of mycobacteriophage L5: applications for genetic transformation of mycobacteria. Mol Microbiol. 1993;7(3):407–17.
Petrova ZO, Broussard GW, Hatfull GF. Mycobacteriophage-repressor-mediated immunity as a selectable genetic marker: Adephagia and BPs repressor selection. Microbiology. 2015;161(Pt 8):1539.
Hatfull GF, Pedulla ML, Jacobs-Sera D, Cichon PM, Foley A, Ford ME, Gonda RM, Houtz JM, Hryckowian AJ, Kelchner VA, Namburi S. Exploring the mycobacteriophage metaproteome: phage genomics as an educational platform. PLoS Genet. 2006;2(6):e92.
Hatfull GF, Jacobs-Sera D, Lawrence JG, Pope WH, Russell DA, Ko CC, Weber RJ, Patel MC, Germane KL, Edgar RH, Hoyte NN. Comparative genomic analysis of 60 mycobacteriophage genomes: genome clustering, gene acquisition, and gene size. J Mol Biol. 2010;397(1):119–43.
Brudno M, Do CB, Cooper GM, Kim MF, Davydov E, Green ED, Sidow A, Batzoglou S. NISC comparative sequencing program, LAGAN and multi-LAGAN: efficient tools for large-scale multiple alignment of genomic DNA. Genome Res. 2003;13(4):721–31.
Lu S, Wang J, Chitsaz F, Derbyshire MK, Geer RC, Gonzales NR, Gwadz M, Hurwitz DI, Marchler GH, Song JS, Thanki N, Yamashita RA, Yang M, Zhang D, Zheng C, Lanczycki CJ, Marchler-Bauer A. CDD/SPARCLE: the conserved domain database in 2020. Nucleic Acids Res. 2020;48(D1):D265–8.
Yang J, Yan R, Roy A, Xu D, Poisson J, Zhang Y. The I-TASSER suite: protein structure and function prediction. Nat Methods. 2015;12:7–8.
Zhang Y. I-TASSER server for protein 3D structure prediction. BMC Bioinformatics. 2008;9:40.
Pabo CO. Structure of the lambda complex at 2.5 a resolution: details of the repressor-operator interactions. Science. 1988;242(4880):893–9.
Meyer BJ, Kleid DG, Ptashne M. Lambda repressor turns off transcription of its own gene. Proc Natl Acad Sci. 1975;72(12):4785–9.
Shimatake H, Rosenberg M. Purified λ regulatory protein cII positively activates promoters for lysogenic development. Nature. 1981;292(5819):128–32.
Echols H, Green L. Establishment and maintenance of repression by bacteriophage lambda: the role of the cI, cII, and cIII proteins. Proc Natl Acad Sci. 1971;68(9):2190–4.
Chothia C, Lesk AM. The relation between the divergence of sequence and structure in proteins. EMBO J. 1986;5(4):823–6.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.
Chaisson MJ, Pevzner PA. Short read fragment assembly of bacterial genomes. Genome Res. 2008;18(2):324–30.
Russell DA, Hatfull GF. PhagesDB: the actinobacteriophage database. Bioinformatics. 2016;33(5):784–6.
Delcher AL, Harmon D, Kasif S, White O, Salzberg SL. Improved microbial gene identification with GLIMMER. Nucleic Acids Res. 1999;27(23):4636–41.
Lukashin AV, Borodovsky M. GeneMark. Hmm: new solutions for gene finding. Nucleic Acids Res. 1998;26(4):1107–15.
Laslett D, Canback B. ARAGORN, a program to detect tRNA genes and tmRNA genes in nucleotide sequences. Nucleic Acids Res. 2004;32(1):11–6.
Lowe TM, Eddy SR. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997;25(5):955–64.
Rinehart CA, Gaffney BL, Smith JR, Wood JD. "PECAAN, Phage Evidence Collection And Annotation Network.". 2016. https://discover.kbrinsgd.org/login. Accessed 23 Feb 2018.
Jones DT. Protein secondary structure prediction based on position-specific scoring matrices. J Mol Biol. 1999;292(2):195–202.
Brister JR, Ako-Adjei D, Bao Y, Blinkova O. NCBI viral genomes resource. Nucleic Acids Res. 2015;43(Database issue):D571–7.
Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE. The Protein Data Bank. Nucleic Acids Res. 2000;28:235–42.
Armougom F, Moretti S, Poirot O, Audic S, Dumas P, Schaeli B, Keduas V, Notredame C. Expresso: automatic incorporation of structural information in multiple sequence alignments using 3D-Coffee. Nucleic Acids Res. 2006;34(Web Server issue):W604–W608.
Trifinopoulos J, Nguyen L-T, von Haeseler A, Minh BQ. W-IQ-TREE: a fast online phylogenetic tool for maximum likelihood analysis. Nucleic Acids Res. 2016;44(W1):W232–5.
Boratyn GM, Schäffer AA, Agarwala R, Altschul SF, Lipman DJ, Madden TL. Domain enhanced lookup time accelerated BLAST. Biol Direct. 2012;7(1):12.
Cole C, Barber JD, Barton GJ. The Jpred 3 secondary structure prediction server. Nucleic Acids Res. 2008;36:W197–201.
Buchan DW, Minneci F, Nugent TC, Bryson K, Jones DT. Scalable web services for the PSIPRED protein analysis workbench. Nucleic Acids Res. 2013;41(W1):W349–57.
Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, Ferrin TE. UCSF Chimera: a visualization system for exploratory research and analysis. J Comput Chem. 2004;25(13):1605–12.
We would like to thank the Hatfull Laboratory at the University of Pittsburgh for invaluable advice and protocols.
The authors acknowledge support for this project from an Institutional Development Award (IDeA) from National Institute of General Medical Sciences, National Institutes of Health Grant 2P20GM103436. The National Institutes of Health FWO played no role in the design of the study, nor in the collection, analysis, and interpretation of data, nor in writing the manuscript. Molecular graphics and analyses performed with UCSF Chimera, developed by the Resource for Biocomputing, Visualization, and Informatics at the University of California, San Francisco, with support from NIH P41-GM103311.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Gudlavalleti, B.S., Phung, T., Barton, C.L. et al. Whole genome sequencing identifies an allele responsible for clear vs. turbid plaque morphology in a Mycobacteriophage. BMC Microbiol 20, 148 (2020). https://doi.org/10.1186/s12866-020-01833-4
- Mycobacterium smegmatis
- Immunity repressor protein