Terminase Genes Research Articles

DETECTION OF ELEPHANT ENDOTHELIOTROPIC HERPESVIRUS 1A IN ARCHIVAL TISSUE USING RNASCOPE® IN SITU HYBRIDIZATION.

Hemorrhagic disease due to elephant endotheliotropic herpesvirus infection (EEHV-HD) is an important cause of calf mortality in managed and free-ranging Asian (Elephas maximus) and African elephant (Loxodonta spp.) populations. Consequently, infection has profound implications for elephant population growth and sustainability. The mechanisms of disease caused by EEHV (i.e., infection, dissemination, shedding, latency) are relatively undefined, in part because of a lack of robust validated assays for detecting viral gene products in relevant samples. To address this issue, we used RNAscope® in situ hybridization (ISH) based on EEHV1A DNA polymerase and terminase genes to detect EEHV1A RNA in archival formalin-fixed, paraffin-embedded Asian elephant heart and tongue from PCR-confirmed cases (n = 4) of EEHV-HD and Asian elephants (n = 2) that died from other causes. EEHV1A-positive cases had positive hybridization signal in endothelial cell nuclei of both tissues for both DNA polymerase and terminase. EEHV-negative cases lacked signal. In positive cases, the number of positive nuclei was manually assessed to provide an estimate of the viral load and compare sensitivity of the two probes. In all cases, heart had greater signal than tongue for both probes (Wilcoxon rank test; P ≤ 0.01). Overall, terminase hybridization signal was greater than DNA polymerase signal (Wilcoxon rank test; P ≤ 0.01). Results indicate RNAscope ISH is a valuable tool for detection of EEHV in archival samples and for confirming infection. Additionally, the terminase gene is the optimal target and heart is preferable to tongue for detection in cases of EEHV-HD. Results will inform future investigations of viral tropism in EEHV-HD cases due to EEHV1A.

Muramidase, nuclease, or hypothetical protein genes intervene between paired genes encoding DNA packaging terminase and portal proteins in Wolbachia phages and prophages.

Genomes of the obligate intracellular alpha proteobacterium Wolbachia pipientis often encode prophage-like regions, and in a few cases, purified particles have been recovered. Because the structure of a conserved WO phage genome has been difficult to establish, we examined paired terminase and portal genes in Wolbachia phages and prophages, relative to those encoded by the gene transfer agent RcGTA from the free-living alpha proteobacterium Rhodobacter capsulatus. Terminase and portal proteins from Wolbachia have higher similarity to orthologs encoded by RcGTA than to orthologs encoded by bacteriophage lambda. In lambdoid phages, these proteins play key roles in assembly of mature phage particles, while in less well-studied gene transfer agents, terminase and portal proteins package random fragments of bacterial DNA, which could confound elucidation of WO phage genomes. In WO phages and prophages, terminase genes followed by a short gpW gene may be separated from the downstream portal gene by open-reading frames encoding a GH_25 hydrolase/muramidase, a PD-(D/E)XK nuclease, a hypothetical protein and/or a RelE/ParE toxin-antitoxin module. These aspects of gene organization, coupled with evidence for a low, non-inducible yield of WO phages, and the small size of WO phage particles described in the literature raise the possibility that Wolbachia prophage regions participate in processes that extend beyond conventional bacteriophage lysogeny and lytic replication. These intervening genes, and their possible relation to functions associated with GTAs, may contribute to variability among WO phage genomes recovered from physical particles and impact the ability of WO phages to act as transducing agents.

Development of in situ hybridization for detection of elephant endotheliotropic herpesvirus in Asian elephants.

Elephant endotheliotropic herpesvirus (EEHV) is one of the most important viral infectious diseases affecting the elephant population worldwide, especially juveniles and young adults. We developed a chromogenic in situ hybridization (ISH) test for detection of EEHV in Asian elephants ( Elephas maximus). Digoxigenin (DIG) DNA probes from the polymerase and terminase genes of EEHV were synthesized using a PCR DIG-labeling method, and detection of hybridized probe to target EEHV DNA was carried out by anti-DIG immunolabeling. Distribution of EEHV-1A and EEHV-4 genomes was found to be prominent in mononuclear phagocytic cells of spleen and endothelial cells of visceral organs. ISH enables the detection of EEHV infection and has applications in understanding pathogenesis of EEHV in Asian elephants.

Co-infection of Acipenserid herpesvirus 2 (AciHV-2) and Streptococcus iniae in cultured white sturgeon Acipenser transmontanus.

A mortality event in cultured white sturgeon Acipenser transmontanus (Richardson, 1836) sub-adults was investigated. After transfer between farms, high mortality was observed in fish, associated with back arching, abnormal swimming, and ulcerative skin lesions. Necropsy of moribund individuals revealed hemorrhagic ascites and petechial hemorrhages in the coelomic peritoneum and serosa of internal organs. Acipenserid herpesvirus 2 (AciHV-2) was isolated from external tissue samples, then identified and genotyped by sequencing of the terminase and polymerase genes. In addition, Streptococcus iniae was recovered from internal organs of affected fish. Histologic changes were limited to interstitial hematopoietic areas of the kidney and consisted of small foci of necrosis accompanied by fibrin deposition, minimal inflammatory response, and small numbers of bacterial cocci compatible with streptococci. Identity was confirmed by partial sequencing of the 16S rRNA, rpoB, and gyrB genes. Genetic fingerprinting demonstrated a genetic profile distinct from S. iniae isolates recovered from previous outbreaks in wild and cultured fish in North America, South America, and the Caribbean. Although the isolates were resistant to white sturgeon complement in serum killing assays, in vivo challenges failed to fulfill Koch's postulates. However, the clinical presentation, coupled with consistent recovery of S. iniae and AciHV-2 from moribund fish, suggests viral and bacterial co-infection were the proximate cause of death. To our knowledge, this represents the first report of AciHV-2 and S. iniae co-infection in cultured white sturgeon.

Inducible Prophage Mutant of Escherichia coli Can Lyse New Host and the Key Sites of Receptor Recognition Identification.

The use of bacteriophages as therapeutic agents is hindered by their narrow and specific host range, and by a lack of the knowledge concerning the molecular mechanism of receptor recognition. Two P2-like coliphages, named P88 and pro147, were induced from Escherichia coli strains K88 and DE147, respectively. A comparison of the genomes of these two and other P2-like coliphages obtained from GenBank showed that the tail fiber protein genes, which are the key genes for receptor recognition in other myoviridae phages, showed more diversity than the conserved lysin, replicase, and terminase genes. Firstly, replacing hypervariable region 2 (HR2: amino acids 716–746) of the tail fiber protein of P88 with that of pro147 changed the host range of P88. Then, replacing six amino acids in HR2 with the corresponding residues from pro147 altered the host range only in these mutants with changes at position 730 (leucine) and 744 (glutamic acid). Thus, we predicted that these amino acids are vital to establish the host range of P88. This study provided a vector of lysogenic bacteria that could be used to change or expand the phage host range of P88. These results illustrated that, in P2-like phage P88, the tail fiber protein determined the receptor recognition. Amino acids 716–746 and the amino acids at positions 730 and 744 were important for receptor recognition.

Genes essential for the morphogenesis of the Shiga toxin 2-transducing phage from Escherichia coli O157:H7.

Shiga toxin 2 (Stx2), one of the most important virulence factors of enterohaemorrhagic Escherichia coli (EHEC), is encoded by phages. These phages (Stx2 phages) are often called lambda-like. However, most Stx2 phages are short-tailed, thus belonging to the family Podoviridae, and the functions of many genes, especially those in the late region, are unknown. In this study, we performed a systematic genetic and morphological analysis of genes with unknown functions in Sp5, the Stx2 phage from EHEC O157:H7 strain Sakai. We identified nine essential genes, which, together with the terminase genes, determine Sp5 morphogenesis. Four of these genes most likely encoded portal, major capsid, scaffolding and tail fiber proteins. Although exact roles/functions of the other five genes are unknown, one was involved in head formation and four were required for tail formation. One of the four tail genes encoded an unusually large protein of 2,793 amino-acid residues. Two genes that are likely required to maintain the lysogenic state were also identified. Because the late regions of Stx2 phages from various origins are highly conserved, the present study provides an important basis for better understanding the biology of this unique and medically important group of bacteriophages.

Initial Detection and Molecular Characterization of Namaycush Herpesvirus (Salmonid Herpesvirus 5) in Lake Trout.

A novel herpesvirus was found by molecular methods in samples of Lake Trout Salvelinus namaycush from Lake Erie, Pennsylvania, and Lake Ontario, Keuka Lake, and Lake Otsego, New York. Based on PCR amplification and partial sequencing of polymerase, terminase, and glycoprotein genes, a number of isolates were identified as a novel virus, which we have named Namaycush herpesvirus (NamHV) salmonid herpesvirus 5 (SalHV5). Phylogenetic analyses of three NamHV genes indicated strong clustering with other members of the genus Salmonivirus, placing these isolates into family Alloherpesviridae. The NamHV isolates were identical in the three partially sequenced genes; however, they varied from other salmonid herpesviruses in nucleotide sequence identity. In all three of the genes sequenced, NamHV shared the highest sequence identity with Atlantic Salmon papillomatosis virus (ASPV; SalHV4) isolated from Atlantic Salmon Salmo salar in northern Europe, including northwestern Russia. These results lead one to believe that NamHV and ASPV have a common ancestor that may have made a relatively recent host jump from Atlantic Salmon to Lake Trout or vice versa. Partial nucleotide sequence comparisons between NamHV and ASPV for the polymerase and glycoprotein genes differ by >5% and >10%, respectively. Additional nucleotide sequence comparisons between NamHV and epizootic epitheliotropic disease virus (EEDV/SalHV3) in the terminase, glycoprotein, and polymerase genes differ by >5%, >20%, and >10%, respectively. Thus, NamHV and EEDV may be occupying discrete ecological niches in Lake Trout. Even though NamHV shared the highest genetic identity with ASPV, each of these viruses has a separate host species, which also implies speciation. Additionally, NamHV has been detected over the last 4 years in four separate water bodies across two states, which suggests that NamHV is a distinct, naturally replicating lineage. This, in combination with a divergence in nucleotide sequence from EEDV, indicates that NamHV is a new species in the genus Salmonivirus. Received April 20, 2015; accepted October 11, 2015.

Atlantic salmon papillomatosis in Russia and molecular characterization of the associated herpesvirus

Papillomatosis of Atlantic salmon Salmo salar has been reported for decades in Russia, Scandinavia and Scotland. The disease is typically benign although heavy losses have occasionally been reported. A herpesviral etiology has been suggested based on ultrastructural evidence; however, the virus has not been isolated or genetically characterized. In this study, we provide the first viral sequences detected in the papillomas from diseased Russian Atlantic salmon. Phylogenetic analyses, based on the partial sequences of the herpesviral polymerase and terminase genes, supported the virus as a novel member of the genus Salmonivirus within the family Alloherpesviridae. The sequences of the Atlantic salmon papillomatosis virus differ markedly from those of the 3 known salmoniviruses; therefore, the authors propose the species designation Salmonid herpesvirus 4 to be considered for approval by the International Committee on Taxonomy of Viruses.

Natural zinc ribbon HNH endonucleases and engineered zinc finger nicking endonuclease

Many bacteriophage and prophage genomes encode an HNH endonuclease (HNHE) next to their cohesive end site and terminase genes. The HNH catalytic domain contains the conserved catalytic residues His-Asn-His and a zinc-binding site [CxxC]2. An additional zinc ribbon (ZR) domain with one to two zinc-binding sites ([CxxxxC], [CxxxxH], [CxxxC], [HxxxH], [CxxC] or [CxxH]) is frequently found at the N-terminus or C-terminus of the HNHE or a ZR domain protein (ZRP) located adjacent to the HNHE. We expressed and purified 10 such HNHEs and characterized their cleavage sites. These HNHEs are site-specific and strand-specific nicking endonucleases (NEase or nickase) with 3- to 7-bp specificities. A minimal HNH nicking domain of 76 amino acid residues was identified from Bacillus phage γ HNHE and subsequently fused to a zinc finger protein to generate a chimeric NEase with a new specificity (12–13 bp). The identification of a large pool of previously unknown natural NEases and engineered NEases provides more ‘tools’ for DNA manipulation and molecular diagnostics. The small modular HNH nicking domain can be used to generate rare NEases applicable to targeted genome editing. In addition, the engineered ZF nickase is useful for evaluation of off-target sites in vitro before performing cell-based gene modification.

Molecular prophage typing of avian pathogenic Escherichia coli

Escherichia coli prophages confer virulence and resistance to physico-chemical, nutritional, and antibiotic stresses on their hosts, and they enhance the evolution of E. coli. Thus, studies on profiles of E. coli prophages are valuable to understand the population structure and evolution of E. coli pathogenicity. Large terminase genes participate in phage genome packaging and are one of the cornerstones for the identification of prophages. Thus, we designed primers to detect 16 types of large terminase genes and analyzed the genomes of 48 E. coli and Shigella reference strains for the prophage markers. We also investigated the distribution of the 16 prophage markers among 92 avian pathogenic E. coli (APEC) strains. APEC strains were classified into 61 prophage types (PPTs). Each strain was different from the reference strains as measured by the PPTs and from the frequency of each prophage marker. Investigation of the distribution of prophage-related serum resistance (bor), toxin (stx1 and cdtI), and T3SS effector (lom, espK, sopE, nleB, and ospG) genes revealed the presence of bor (44.1%), lom (95.5%) and cdtI (9.1%) in APEC strains with related prophages. Therefore, the molecular prophage typing method may be useful to understand population structure and evolution of E. coli pathogenicity, and further studies on the mobility of the prophages and the roles of virulence genes in APEC pathogenicity may be valuable.

Genetic diversity of equine gammaherpesviruses (γ-EHV) and isolation of a syncytium forming EHV-2 strain from a horse in Iceland

The horse population in Iceland is a special breed, isolated from other equines for at least one thousand years. This provides an exceptional opportunity to investigate old and new pathogens in a genetically closed herd. Both types of equine gammaherpesviruses, EHV-2 and EHV-5, are common in Iceland.Genetic variation was examined by sequencing four genes, glycoprotein B (gB), glycoprotein H (gH), DNA polymerase and DNA terminase for 12 Icelandic and seven foreign EHV-2 strains. One Icelandic virus isolate, gEHV-Dv, induced syncytium formation, an uncharacteristic cytopathy for EHV-2 in equine kidney cells. When sequenced, the glycoprotein genes were different from both EHV-2 and EHV-5, but the polymerase and terminase genes had 98–99% identity to EHV-2. Therefore the gEHV-Dv strain can be considered a variant of EHV-2.Substantial genetic variability was seen within the EHV-2 glycoprotein genes but limited in the polymerase and terminase genes. The Icelandic EHV-2 strains do not seem to differ phylogenetically from the foreign viruses, despite isolation for over a thousand years.

Virus genes in Arctic marine bacteria identified by metagenomic analysis

The relationship between viruses and prokaryotes in the Arctic appears to differ from that at lower latitudes because of low temperature and low abundances of viruses and bacteria, but the impact of viruses on Arctic bacterial genomes is unclear. In the present study, genome fragments of bacteria from the western Arctic Ocean were se quenced to explore the prevalence of virus genes in the genomes of bacteria inhabiting these perennially cold waters. Arctic bacterial genomes were sampled by cloning bacterial environmental DNA into a fosmid vector, and virus DNA within the bacterial genomes was identified using marine viral meta genomes to query the bacterial genomes using the basic local alignment search tool (BLAST). Virus DNA flanked by bacterial DNA was identified as being part of bacterial genomes. Virus DNA was ~3-fold more abundant in Arctic bacterial genomes than in bacterial genomes from Monterey Bay and ~10-fold more abundant than in bacterial genomes from Antarctic waters. Phage terminase genes involved in packaging DNA into phage capsids were the most abundant gene family identified in the Arctic bacterial genomes. Viruses appear to have a larger impact on prokaryotic communities in the Arctic than what might be inferred from the low bacterial and viral abundances in this high latitude ocean.

A Novel Herpesvirus in 3 Species of Pheasants: Mountain Peacock Pheasant (Polyplectron inopinatum), Malayan Peacock Pheasant (Polyplectron malacense), and Congo Peafowl (Afropavo congensis)

The mountain peacock pheasant (Polyplectron inopinatum), the Malayan peacock pheasant (Polyplectron malacense), and the Congo peafowl (Afropavo congensis) are all listed as vulnerable to extinction under the International Union for Conservation of Nature Red List of Threatened Species. Here the authors report fatal infection with a novel herpesvirus in all 3 species of birds. DNA was extracted from the livers of birds with hepatocellular necrosis and intranuclear eosinophilic inclusions consistent with herpesvirus infection. Based on degenerate herpesvirus primers and polymerase chain reaction, 220- and 519-base pair products of the herpes DNA polymerase and DNA terminase genes, respectively, were amplified. Sequence analysis revealed that all birds were likely infected with the same virus. At the nucleotide level, the pheasant herpesvirus had 92% identity with gallid herpesvirus 3 and 77.7% identity with gallid herpesvirus 2. At the amino acid level, the herpes virus had 93.8% identity with gallid herpesvirus 3 and 89.4% identity with gallid herpesvirus 2. These findings indicate that the closest relative to this novel herpesvirus is gallid herpesvirus 3, a nonpathogenic virus used widely in a vaccine against Marek's disease. In situ hybridization using probes specific to the peacock pheasant herpesvirus DNA polymerase revealed strong intranuclear staining in the necrotic liver lesions of an infected Malayan peacock pheasant but no staining in normal liver from an uninfected bird. The phasianid herpesvirus reported here is a novel member of the genus Mardivirus of the subfamily Alphaherpesvirinae and is distinct from other galliform herpesviruses.

Isolation and Partial Sequencing of Equid Herpesvirus 5 from a Horse in Iceland

Horses are hosts to 2 types of gammaherpesviruses, Equid herpesvirus 2 and 5 (EHV-2 and EHV-5, respectively). Both EHV-2 and EHV-5 are common in horses in Iceland. An Icelandic EHV-5 isolate was recovered by sequential culture in primary fetal horse kidney and rabbit kidney cells. Glycoprotein B, glycoprotein H, and DNA terminase genes of the isolate were fully sequenced, and the DNA polymerase gene was partly sequenced. To date, the glycoprotein B gene of EHV-5 was the only gene that has been reported to be completely sequenced in addition to small parts of the glycoprotein H, DNA polymerase, and DNA terminase genes. The present report, therefore, is a significant addition to previously reported EHV-5 sequences.

Isolation and characterization of bacteriophage ΦBP fromPaenibacillus polymyxaCCM 7400

A bacteriophage PhiBP infecting Paenibacillus polymyxa CCM 7400 was isolated from culture lysate. Electron microscopy of lysate samples revealed the presence of bacteriophage particles with polyhedral heads 56 nm in diameter and flexible noncontractile tails 144 nm in length. The profile of PhiBP structural proteins resembles that of other bacteriophages. The PhiBP genome consists of double-stranded DNA of 43-kbp size. Homology search of sequenced DNA fragments from EcoRI digest revealed regions with significant similarity to other known bacteriophage genes. Regions similar to phage terminase genes were identified within the 1.2-kbp fragment. Three lytic genes, two holin genes and one endolysin gene were identified within the 2.5-kbp fragment. We tested the isolates of P. polymyxa CCM 7400 for the presence of phage DNA on bacterial chromosome using PCR amplification with primers derived from proposed terminase and holin gene sequences. We confirmed the presence of PhiBP DNA on P. polymyxa chromosome by Southern hybridization. The bacteriophage PhiBP was capable of causing lysis of a P. polymyxaPhiBP lysogen despite the presence of the phage DNA on bacterial chromosome. Therefore, we concluded that PhiBP was a virulent mutant phage.

Genomic and biological analysis of phage Xfas53 and related prophages of Xylella fastidiosa.

We report the plaque propagation and genomic analysis of Xfas53, a temperate phage of Xylella fastidiosa. Xfas53 was isolated from supernatants of X. fastidiosa strain 53 and forms plaques on the sequenced strain Temecula. Xfas53 forms short-tailed virions, morphologically similar to podophage P22. The 36.7-kb genome is predicted to encode 45 proteins. The Xfas53 terminase and structural genes are related at a protein and gene order level to P22. The left arm of the Xfas53 genome has over 90% nucleotide identity to multiple prophage elements of the sequenced X. fastidiosa strains. This arm encodes proteins involved in DNA metabolism, integration, and lysogenic control. In contrast to Xfas53, each of these prophages encodes head and DNA packaging proteins related to the siphophage lambda and tail morphogenesis proteins related to those of myophage P2. Therefore, it appears that Xfas53 was formed by recombination between a widespread family of X. fastidiosa P2-related prophage elements and a podophage distantly related to phage P22. The lysis cassette of Xfas53 is predicted to encode a pinholin, a signal anchor and release (SAR) endolysin, and Rz and Rz1 equivalents. The holin gene encodes a pinholin and appears to be subject to an unprecedented degree of negative regulation at both the level of expression, with rho-independent transcriptional termination and RNA structure-dependent translational repression, and the level of holin function, with two upstream translational starts predicted to encode antiholin products. A notable feature of Xfas53 and related prophages is the presence of 220- to 390-nucleotide degenerate tandem direct repeats encoding putative DNA binding proteins. Additionally, each phage encodes at least two BroN domain-containing proteins possibly involved in lysogenic control. Xfas53 exhibits unusually slow adsorption kinetics, possibly an adaptation to the confined niche of its slow-growing host.

Phylogenetic relationships in the family Alloherpesviridae

Phylogenetic relationships among herpesviruses (HVs) of mammals, birds, and reptiles have been studied extensively, whereas those among other HVs are relatively unexplored. We have reconstructed the phylogenetic relationships among 13 fish and amphibian HVs using maximum likelihood and Bayesian analyses of amino acid sequences predicted from parts of the DNA polymerase and terminase genes. The relationships among 6 of these viruses were confirmed using the partial DNA polymerase data plus the complete sequences of the terminase, helicase, and triplex protein genes; the position of these viruses among all other sequenced HVs was also investigated using the complete terminase gene. The results established the monophyly of the fish and amphibian HVs (Alloherpesviridae) separate from the HVs of mammals, birds, and reptiles (Herpesviridae) and the single recognized HV of bivalve mollusks (Malacoherpesviridae) in the order Herpesvirales. Two major clades in the family Alloherpesviridae were recognized: one consisting of viruses from cyprinid and anguillid hosts and the other of viruses from ictalurid, salmonid, acipenserid, and ranid hosts. A comparison of virus and host phylogenies suggested that closely related HVs in this family may have coevolved with their hosts, whereas significant codiversification was not apparent for the more distantly related viruses.

The Temperate Marine Phage ΦHAP-1 of Halomonas aquamarina Possesses a Linear Plasmid-Like Prophage Genome

A myovirus-like temperate phage, PhiHAP-1, was induced with mitomycin C from a Halomonas aquamarina strain isolated from surface waters in the Gulf of Mexico. The induced cultures produced significantly more virus-like particles (VLPs) (3.73 x 10(10) VLP ml(-1)) than control cultures (3.83 x 10(7) VLP ml(-1)) when observed with epifluorescence microscopy. The induced phage was sequenced by using linker-amplified shotgun libraries and contained a genome 39,245 nucleotides in length with a G+C content of 59%. The PhiHAP-1 genome contained 46 putative open reading frames (ORFs), with 76% sharing significant similarity (E value of <10(-3)) at the protein level with other sequences in GenBank. Putative functional gene assignments included small and large terminase subunits, capsid and tail genes, an N6-DNA adenine methyltransferase, and lysogeny-related genes. Although no integrase was found, the PhiHAP-1 genome contained ORFs similar to protelomerase and parA genes found in linear plasmid-like phages with telomeric ends. Southern probing and PCR analysis of host genomic, plasmid, and PhiHAP-1 DNA indicated a lack of integration of the prophage with the host chromosome and a difference in genome arrangement between the prophage and virion forms. The linear plasmid prophage form of PhiHAP-1 begins with the protelomerase gene, presumably due to the activity of the protelomerase, while the induced phage particle has a circularly permuted genome that begins with the terminase genes. The PhiHAP-1 genome shares synteny and gene similarity with coliphage N15 and vibriophages VP882 and VHML, suggesting an evolutionary heritage from an N15-like linear plasmid prophage ancestor.

Genome classification by gene distribution: an overlapping subspace clustering approach.

BackgroundGenomes of lower organisms have been observed with a large amount of horizontal gene transfers, which cause difficulties in their evolutionary study. Bacteriophage genomes are a typical example. One recent approach that addresses this problem is the unsupervised clustering of genomes based on gene order and genome position, which helps to reveal species relationships that may not be apparent from traditional phylogenetic methods.ResultsWe propose the use of an overlapping subspace clustering algorithm for such genome classification problems. The advantage of subspace clustering over traditional clustering is that it can associate clusters with gene arrangement patterns, preserving genomic information in the clusters produced. Additionally, overlapping capability is desirable for the discovery of multiple conserved patterns within a single genome, such as those acquired from different species via horizontal gene transfers. The proposed method involves a novel strategy to vectorize genomes based on their gene distribution. A number of existing subspace clustering and biclustering algorithms were evaluated to identify the best framework upon which to develop our algorithm; we extended a generic subspace clustering algorithm called HARP to incorporate overlapping capability. The proposed algorithm was assessed and applied on bacteriophage genomes. The phage grouping results are consistent overall with the Phage Proteomic Tree and showed common genomic characteristics among the TP901-like, Sfi21-like and sk1-like phage groups. Among 441 phage genomes, we identified four significantly conserved distribution patterns structured by the terminase, portal, integrase, holin and lysin genes. We also observed a subgroup of Sfi21-like phages comprising a distinctive divergent genome organization and identified nine new phage members to the Sfi21-like genus: Staphylococcus 71, phiPVL108, Listeria A118, 2389, Lactobacillus phi AT3, A2, Clostridium phi3626, Geobacillus GBSV1, and Listeria monocytogenes PSA.ConclusionThe method described in this paper can assist evolutionary study through objectively classifying genomes based on their resemblance in gene order, gene content and gene positions. The method is suitable for application to genomes with high genetic exchange and various conserved gene arrangement, as demonstrated through our application on phages.

DNA recombination with a heterospecific Cre homolog identified from comparison of the pac-c1 regions of P1-related phages

Sequencing of the 7 kb immC region from four P1-related phages identified a novel DNA recombinase that exhibits many Cre-like characteristics, including recombination in mammalian cells, but which has a distinctly different DNA specificity. DNA sequence comparison to the P1 immC region showed that all phages had related DNA terminase, C1 repressor and DNA recombinase genes. Although these genes from phages P7, phi(w39) and p15B were highly similar to those from P1, those of phage D6 showed significant divergence. Moreover, the D6 sequence showed evidence of DNA deletion and substitution in this region relative to the other phages. Characterization of the D6 site-specific DNA recombinase (Dre) showed that it was a tyrosine recombinase closely related to the P1 Cre recombinase, but that it had a distinct DNA specificity for a 32 bp DNA site (rox). Cre and Dre are heterospecific: Cre did not catalyze recombination at rox sites and Dre did not catalyze recombination at lox sites. Like Cre, Dre catalyzed both integrative and excisive recombination and required no other phage-encoded proteins for recombination. Dre-mediated recombination in mammalian cells showed that, like Cre, no host bacterial proteins are required for efficient Dre-mediated site-specific DNA recombination.