BACKGROUND
viral integration is a ubiquitous biological process where a virus inserts its genetic material into a host genome.
integration sometimes results in the entire viral genome being incorporated into the host genome, such as occurs with lysogenic bacteriophage or retroviruses  <cit> .
in other cases, only a portion of the viral genome is integrated.
for example, sub-genomic segments of human papilloma virus or hepatitis b virus  genomes are sometimes integrated in cervical and hepatocellular cancers  respectively  <cit> .
by identifying the lost and gained sub-genomic segments of host and virus, potential expression of viral and host genes can be estimated.
there is a growing amount of next generation sequencing  data available for cancer genomes  <cit> .
these datasets allow for massive viral integration analysis in which single virus and host chimeras have been identified and compared  <cit> .
some experimental studies have identified the complete viral integration by pcr  <cit> .
current viral integration detection software programs  <cit>  are designed to identify virus-host chimeras and report their genomic positions.
however, the complete mapping of viral integrations requires the association of two chimeric sequences representing the two virus-host junctions present in each integration event.
summonchimera is a perl program developed to detect chimeras from paired or single end ngs reads then associate two chimeras that describe the integration event.
mapping and association of both virus-host junctions allows the identification of viral and host sequences retained and lost during integration.
implementation
datasets
the single-end ion torrent sequencing reads of a purified salmonella enterica culture, contained  <dig> , <dig> reads with an average length of  <dig> residues, was obtained from russell et al.
<cit>  and is publicly available.
no irb approval is required for these studies.
hepatocellular carcinoma  paired-end whole genome sequencing  reads were generated on an illumina hiseq  <dig> sequencer platform by sung et al.
<cit> .
the reads were publicly available in the european nucleotide archive  .
samples t <dig> and t <dig> from the hcc wgs study, with  <dig> , <dig> and  <dig> , <dig> paired-end reads respectively, were chosen for two reasons.
first, subsets of chimeras were experimentally verified  <cit> .
second, these samples were used by wang et al.
<cit>  to compare current integration detection software  <cit> .
prinseq  <cit>  quality filtering was run on the detected chimeric reads from each dataset to eliminate reads with an average quality less than  <dig>  sequencing sets were run through the alignment subtraction pipeline outlined in figure  <dig> .figure  <dig>
chimera detection pipeline.
red arrows represent unaligned reads.
black arrows represent aligned reads and data.
this figure outlines the standard computational subtraction pipeline used to generate summonchimera input.
mapping of paired-end reads must occur to a database composed of both virus and humangenomes with bowtie <dig>  the blastn steps are ordered for optimal search time.
read counts at each step of the pipeline for each sample can be found in additional file  <dig>
paired-end mapping
default bowtie <dig>  <cit>  settings were used for mapping of discordant paired-end reads.
for both hcc samples, a single bowtie <dig> database was created with the ucsc hg <dig> human genome  and the hbv genome  from the ncbi genome database .
no mapping was run on the single end iontorrent salmonella dataset.
blastn pipeline
the blastn  <cit>  subtraction pipeline  was used for the identification of chimeric reads.
in cases where the reads were mapped, all unmapped reads were extracted from the resulting bam file with ‘samtools –f 4’.
blastn parameters ‘-use_index true’, ‘–outfmt 6’, ‘-word_size 16’, and ‘-perc_identity 95’ are used to minimize output file size and computation time of each step while maximizing the number of detected chimeras.
the input reads are first aligned to the virus database with blastn.
reads with a virus blast hit were then searched against the host database with blastn.
hcc blast databases used the same reference sequences that were used during the mapping step: host hg <dig> and virus hbv.
the salmonella blast databases were composed of the virus genome salmonella bacteriophage re- <dig>  and host genome salmonella enterica  from the ncbi genome database.
blast databases and indices were also created with default parameters.
summonchimera
sc was written such that sequence processing is left to the users.
blast and mapping are all done separate from sc and results are used as input.
figure  <dig> indicates the pipeline to be used for analysis.
sc was run with default parameters using the blast output and sam files as input.
sc’s first steps are to filter confounding results from the blast analysis.
blast alignments are filtered based on user defined parameters to determine which alignments are suitable.
from the remaining blast results, alignments with the highest bit score to virus and host are used to characterize the chimera.
chimeric reads with alignments having equivalent bit scores to multiple positions in virus or host are removed from analysis.
all mapped virus-host discordant pairs identified in the sam file are treated as a chimera.
chimeras were categorized in three ways: strand orientation, organism order, and ambiguous region .
strand orientation indicates whether the virus has integrated its genome parallel or anti-parallel to the host genome.
parallelism was determined by the virus and host reference strand as listed in the ncbi genome genbank record.
a parallel integration was defined as the viral reference strand covalently bonded to the host reference strand.
reference anti-parallel would be the virus non-reference strand linked to the host reference strand.
organism order is the 5’ to 3’ order of the organisms represented in the chimera based on the host reference strand.
integrations are inferred from strand orientation and organism order .
finally, an ambiguous region  is defined as the portion of the chimera that cannot be distinguished as virus or host.
chimeric reads were clustered based on these three characteristics.
blast chimeras were clustered if they shared all three characteristics and occurred at the same host and virus position.
mapped and blast chimeras were clustered if they had equivalent organism order, strand orientation, and occurred within a user defined paired-end read insert size on both virus and host.
insert size is used for clustering chimeric reads because it is hypothetically the farthest possible nucleotide distance between a mapped chimeric read and the chimeric junction as determined by blast.
the hcc dataset insert size was  <dig> nucleotides.
occasionally multiple blast chimera ambiguous regions were associated with a cluster.
in these cases, the most frequent blast ambiguous region was reported.figure  <dig>
chimeric cases.
there are four possible cases of chimeric reads.
arrows point from 5’ to 3’ on the reference strand of the host.
blue arrows represent host genome and red is virus.
the two columns indicate organism order cases, which as illustrated is always determined by comparison to the host reference strand.
there are two possible strand orientations which are determined by the orientation of the host and virus reference strand to each other.
three ambiguous cases for chimeric reads are possible.
blastn finds a single chimeric read and thus there is potential for overlap or void positions in the alignment.
mapped chimeras are two reads with an unknown nucleotide distance between them, as such their ambiguous region is defined as the closest two nucleotides mapped from virus and host.
reference parallel and anti-parallel integrations are reported if they are within the user defined host deletion size.
integration predictions were determined by the association of chimeras within the user defined nucleotide deletion size on the host genome .
nucleotide deletion size is the maximum expected portion of the host genome removed during viral integration.
further, chimeras involved in the integration must have equivalent strand orientation and host flanking virus organism order.
figure 2a rows  illustrate the two cases of chimera pairs that represent integration events.
all chimeras that are not paired are reported as orphan chimeras.
the user determines whether a reported orphan chimera or integration is real based on the context of their study.
RESULTS
sc was tested on two different data sets, a salmonella bacterium culture known to harbor an integrated bacteriophage and hbv integration in hcc.
the salmonella dataset was used to verify our integration inference method by comparison of sc results to a known integration of bacteriophage into a salmonella bacterium  <cit> .
second, two hcc tumor samples were used to show the application of sc on cancer ngs data and its ability to identify and associate experimentally verified chimeras.
bacteriophage integration
in some types of lysogenic bacteriophage integration into the host occurs within the homologous core of attp and attb sites.
attp occurs at the single stranded 5’ and 3’ ends of the linear bacteriophage genome.
when a bacteriophage genome enters the cell it circularizes and a double stranded attp is formed.
the enzymatic integration occurs without loss to the bacteriophage or bacterial genome  <cit> .
we tested sc on a bacterial strain known to harbor an integrated bacteriophage genome.
the genome of a salmonella strain carrying a prophage was sequenced by ion torrent using single end reads  <cit> .
thus, only blastn alignments were utilized for sc input.
sc identified a total of  <dig> chimeric reads and reported one integration with ten possible combinations of 5’ and 3’ ambiguous regions .
genome coverage of the  <dig> chimeric reads associated with the 5’ and 3’ ambiguous regions illustrates the integration prediction .
the remaining  <dig> reads  are individual chimeras and are likely artifacts of the sequencing process based on lack of chimeric junction coverage.
prinseq identified no chimeric reads with an average quality less than  <dig>  the 5’ chimeras were divided between two ambiguous regions and the 3’ chimeric reads were divided into five ambiguous regions.
each ambiguous region had single nucleotide variations which likely resulted from sequencing error neighboring the chimeric junction.
the integration with the most reads for the 5’ and 3’ chimera  was selected from sc output for further analysis.
the selected integration showed that the nucleotide compositions of the 5’ and 3’ chimera ambiguous regions were the same.
further, the chimeras occurred at the 5’ and 3’ ends of the bacteriophage genome.
analysis of the bacterial indices indicated that the bacterial genome remained intact .
sc predicted the identical integration as was determined by wgs dna sequencing  <cit> .figure  <dig>
summonchimera prediction of a salmonella bacteriophage integration.
these are segments taken from the aligned region of a 5’ and 3’ chimera.
blue is bacterial sequence , red is bacteriophage sequence , and green is homology to both.
the numbers indicate the reference sequence nucleotide position.
the ambiguous region  is  <dig> nucleotides long.
the blue  and red  arrows show the 5’ to 3’ direction of the host and virus reference strands.
the number of reads  at each nucleotide position in the 5’ and 3’ chimeras for host and virus are shown.
the ambiguous region of the chimera is counted twice  due to the unknown nature of where the integration occurs.
hepatocellular carcinoma hbv integration
in contrast to integration by lysogenic bacteriophage, the integration of dna tumor virus genomes occur at random positions with respect to both the viral and host genomes.
for example, when hbv is integrated in hcc there is a loss of both viral and human dna.
experimentally determined hbv integrations show no nucleotide specific human integration site  <cit> .
to test the ability of sc to detect hbv integrations we examined two hcc datasets, t <dig> and t <dig>  for which viral chimeras had been pcr verified  <cit> .
this data has been used to verify other chimera detection software  <cit> .
from these datasets sc identified  <dig> chimeric reads including all reported pcr verified chimeras .
prinseq was run on all detected chimeric reads and identified none with average quality less than  <dig>  additionally, sc identified a complete integration event in t <dig> and t <dig> by associating pcr verified chimeras .figure  <dig>
summonchimera detected recurrent hbv integration in hcc.
the human chromosome  <dig> region with the htert gene from the ucsc genome browser  view of the hg <dig> is shown  <cit> .
htert is a known recurrent hbv integration site.
the red vertical lines indicate the positions of hbv integration, as predicted by sc, that are illustrated in b and c.  the t <dig> hcc sample hbv integration has been shown to be reference parallel.
the red arrow points 5’ to 3’ on the virus reference strand and the numbers indicate the determined viral nucleotides included, keeping in mind hbvs circular genome.
green blocks represent the annotated gene coding regions within the integration.
the t <dig> hcc sample has an hbv integration that places the viral genes  in a transcription parallel orientation relative to htert.
both integrations detected by sc occurred within the htert gene, a known recurrent integration site for hbv in liver cancer  <cit>  .
host positions of these integrations are similar but, their viral genetic content is different .
the integration in sample t <dig>  is reference parallel and integrates the hbv non-coding strand into the coding strand of the htert gene.
the hbv integration in sample t <dig>  is reference anti-parallel, incorporating the hbv coding region into the htert coding strand.
thus while occurring in the same host gene these integrations have different biological implications.
this result illustrates that the identification of the complete integration event by sc provides additional genetic information beyond current chimeric detection software.
discussion
large ngs datasets containing genomic dna or rna sequences provides an opportunity to detect and examine viral integration events.
the precise mapping of both virus-host junctions in individual tumor samples allows determination of the host gene in which integration occurs, the estimation of viral gene expression, the presence of virus-host chimeric mrnas, and whether host and viral transcription have diverged from their known form.
sc was developed to map viral integration events from ngs.
sc is designed to keep all analysis under user control and for usability with any ngs dataset.
all chimeras and integrations detected by sc are output into a tab separated text document.
provided fields include ambiguous host and virus positions, aligned virus and host region, chimera strand orientation, organism order, host sequence name, and the nucleotide distance between the chimeras of a predicted integration.
with the provided junctions for the 5’ and 3’ chimeras and strand orientation the lost portion of viral and host genome can be inferred.
sc predictions were validated against a variety of experimentally pertinent data.
sc detected 100% of the verified chimeras from a hcc dataset and predicted two previously unknown integrations within the hcc genomes.
current reports suggest that an integration event may either activate or suppress neighboring host and viral genes  <cit> .
since sc determines the integrated viral genome and lost human genome, this information can be utilized to predict effects on host and viral transcripts.
ngs generates artificial chimeras that could be falsely interpreted as an integration by sc.
assuming artifact chimeras are generated at random the number of reads covering the same junction will be low relative to authentic chimeric junctions.
we developed a model to predict the probability of false-positive integration calls by sc.
the equation and calculations for the salmonella and hcc datasets are supplied   <cit> .
chimeric junctions and artifact chimeras are manually differentiated by read coverage.
it is incumbent upon the user to determine the authenticity of chimeras and precision of called integrations based on the context of their experiment.
recent studies of hpv integrations in cervical cancer have shown a positional correlation to genomic instability.
these genomic instability events included copy-number variations, translocations, and inversions  <cit> .
this genomic instability correlation implicates a potential problem for integration inference.
when a viral integration occurs near any of these events , then the chimeric cases  in sc are not sufficient for integration prediction.
further work is required to include genomic instability mapping into sc in order to describe integrations within these regions.figure  <dig>
summonchimera undetectable integrations.
viral integrations in cancer have been correlated with chromosomal inversion, translocation, and copy-number variation events.
when integrations occur as illustrated, summonchimera cannot predict the genetic content.
arrows point from 5’ to 3’ along the reference strand of each organism chromosome .
segment labels indicate a portion of the host genome.
in the inversion case one of the chimeras will be determined with opposite organism order and strand orientation, thus eluding integration prediction.
in case of translocation, the integration can occur between two separate host chromosomes.
the two chimeras will be detected on separate chromosomes and no integration would be predicted.
finally, the copy-number variation case, if the integration occurs between the duplicated segments the chimeras will be detected in the wrong chromosomal positions and the 3’ chimera will appear to be upstream of the 5’ chimera on the chromosome.
CONCLUSIONS
current viral integration software only detects and reports the position of individual chimeras.
sc reports chimera organism order, strand orientation and associates appropriate chimeras to accurately infer the genetic content of viral integrations .
with the identification of the strand orientation and lost portion of host and viral genome from an integration event, the molecular and genetic consequences of integration can be predicted.
availability and requirements
project name: summonchimera
project home page:http://pipaslab.webfactional.com/wp/wp-content/uploads/2014/01/summonchimera.tar.gz
operating system: linux, mac os x, windows
programming language: perl
other requirements: blast +  <dig> .
<dig> and bowtie <dig>  <dig> .0
license: gnu
restrictions for use by non-academics: none
additional files
additional file 1:
read counts for figure
1
computational subtraction pipeline.
description of dataset: read counts for aligned and unaligned at each step of the pipeline.
additional file 2:
summonchimera output for salmonella.
description of dataset: all reported integrations and orphan chimeras for the salmonella data.
additional file 3:
summonchimera output for hepatocellular carcinoma t <dig>  description of dataset: all reported integrations and orphan chimeras for the t <dig> data.
additional file 4:
summonchimera output for hepatocellular carcinoma t <dig>  description of dataset: all reported integrations and orphan chimeras for the t <dig> data.
additional file 5:
calculations of false-positive integration calls by summonchimera.
description of dataset: contains developed equation and calculations for all three datasets.
additional file 6:
summonchimera.
description of dataset: contains summonchimera source code and test data.
abbreviations
hpvhuman papilloma virus
hbvhepatitis b virus
scsummonchimera
hcchepatocellular carcinoma
wgswhole genome sequencing
ngsnext generation sequencing
competing interests
the authors declare that they have no competing interests.
authors’ contributions
jmp had the initial idea to associate two chimeras to identify integration and identified testable cases.
jpk developed the method and code.
jpk wrote the manuscript and jmp edited.
both authors read and approved the final manuscript.

 viral integration is a ubiquitous biological process where a virus inserts its genetic material into a host genome. integration sometimes results in the entire viral genome being incorporated into the host genome, such as occurs with lysogenic bacteriophage or retroviruses. in other cases, only a portion of the viral genome is integrated. for example, sub-genomic segments of human papilloma virus or hepatitis b virus genomes are sometimes integrated in cervical and hepatocellular cancers respectively. by identifying the lost and gained sub-genomic segments of host and virus, potential expression of viral and host genes can be estimated. there is a growing amount of next generation sequencing data available for cancer genomes. these datasets allow for massive viral integration analysis in which single virus and host chimeras have been identified and compared. some experimental studies have identified the complete viral integration by pcr. current viral integration detection software programs are designed to identify virus-host chimeras and report their genomic positions. however, the complete mapping of viral integrations requires the association of two chimeric sequences representing the two virus-host junctions present in each integration event. summonchimera is a perl program developed to detect chimeras from paired or single end ngs reads then associate two chimeras that describe the integration event. mapping and association of both virus-host junctions allows the identification of viral and host sequences retained and lost during integration. implementation datasets the single-end ion torrent sequencing reads of a purified salmonella enterica culture, contained, reads with an average length of residues, was obtained from russell et al. and is publicly available. no irb approval is required for these studies. hepatocellular carcinoma paired-end whole genome sequencing reads were generated on an illumina hiseq sequencer platform by sung et al.. the reads were publicly available in the european nucleotide archive. samples t and t from the hcc wgs study, with, and, paired-end reads respectively, were chosen for two reasons. first, subsets of chimeras were experimentally verified. second, these samples were used by wang et al. to compare current integration detection software. prinseq quality filtering was run on the detected chimeric reads from each dataset to eliminate reads with an average quality less than sequencing sets were run through the alignment subtraction pipeline outlined in figure .figure chimera detection pipeline. red arrows represent unaligned reads. black arrows represent aligned reads and data. this figure outlines the standard computational subtraction pipeline used to generate summonchimera input. mapping of paired-end reads must occur to a database composed of both virus and humangenomes with bowtie the blastn steps are ordered for optimal search time. read counts at each step of the pipeline for each sample can be found in additional file paired-end mapping default bowtie settings were used for mapping of discordant paired-end reads. for both hcc samples, a single bowtie database was created with the ucsc hg human genome and the hbv genome from the ncbi genome database. no mapping was run on the single end iontorrent salmonella dataset. blastn pipeline the blastn subtraction pipeline was used for the identification of chimeric reads. in cases where the reads were mapped, all unmapped reads were extracted from the resulting bam file with ‘samtools –f 4’. blastn parameters ‘-use_index true’, ‘–outfmt 6’, ‘-word_size 16’, and ‘-perc_identity 95’ are used to minimize output file size and computation time of each step while maximizing the number of detected chimeras. the input reads are first aligned to the virus database with blastn. reads with a virus blast hit were then searched against the host database with blastn. hcc blast databases used the same reference sequences that were used during the mapping step: host hg and virus hbv. the salmonella blast databases were composed of the virus genome salmonella bacteriophage re- and host genome salmonella enterica from the ncbi genome database. blast databases and indices were also created with default parameters. summonchimera sc was written such that sequence processing is left to the users. blast and mapping are all done separate from sc and results are used as input. figure  indicates the pipeline to be used for analysis. sc was run with default parameters using the blast output and sam files as input. sc’s first steps are to filter confounding results from the blast analysis. blast alignments are filtered based on user defined parameters to determine which alignments are suitable. from the remaining blast results, alignments with the highest bit score to virus and host are used to characterize the chimera. chimeric reads with alignments having equivalent bit scores to multiple positions in virus or host are removed from analysis. all mapped virus-host discordant pairs identified in the sam file are treated as a chimera. chimeras were categorized in three ways: strand orientation, organism order, and ambiguous region. strand orientation indicates whether the virus has integrated its genome parallel or anti-parallel to the host genome. parallelism was determined by the virus and host reference strand as listed in the ncbi genome genbank record. a parallel integration was defined as the viral reference strand covalently bonded to the host reference strand. reference anti-parallel would be the virus non-reference strand linked to the host reference strand. organism order is the 5’ to 3’ order of the organisms represented in the chimera based on the host reference strand. integrations are inferred from strand orientation and organism order. finally, an ambiguous region is defined as the portion of the chimera that cannot be distinguished as virus or host. chimeric reads were clustered based on these three characteristics. blast chimeras were clustered if they shared all three characteristics and occurred at the same host and virus position. mapped and blast chimeras were clustered if they had equivalent organism order, strand orientation, and occurred within a user defined paired-end read insert size on both virus and host. insert size is used for clustering chimeric reads because it is hypothetically the farthest possible nucleotide distance between a mapped chimeric read and the chimeric junction as determined by blast. the hcc dataset insert size was nucleotides. occasionally multiple blast chimera ambiguous regions were associated with a cluster. in these cases, the most frequent blast ambiguous region was reported.figure chimeric cases. there are four possible cases of chimeric reads. arrows point from 5’ to 3’ on the reference strand of the host. blue arrows represent host genome and red is virus. the two columns indicate organism order cases, which as illustrated is always determined by comparison to the host reference strand. there are two possible strand orientations which are determined by the orientation of the host and virus reference strand to each other. three ambiguous cases for chimeric reads are possible. blastn finds a single chimeric read and thus there is potential for overlap or void positions in the alignment. mapped chimeras are two reads with an unknown nucleotide distance between them, as such their ambiguous region is defined as the closest two nucleotides mapped from virus and host. reference parallel and anti-parallel integrations are reported if they are within the user defined host deletion size. integration predictions were determined by the association of chimeras within the user defined nucleotide deletion size on the host genome. nucleotide deletion size is the maximum expected portion of the host genome removed during viral integration. further, chimeras involved in the integration must have equivalent strand orientation and host flanking virus organism order. figure 2a rows illustrate the two cases of chimera pairs that represent integration events. all chimeras that are not paired are reported as orphan chimeras. the user determines whether a reported orphan chimera or integration is real based on the context of their study.  sc was tested on two different data sets, a salmonella bacterium culture known to harbor an integrated bacteriophage and hbv integration in hcc. the salmonella dataset was used to verify our integration inference method by comparison of sc results to a known integration of bacteriophage into a salmonella bacterium. second, two hcc tumor samples were used to show the application of sc on cancer ngs data and its ability to identify and associate experimentally verified chimeras. bacteriophage integration in some types of lysogenic bacteriophage integration into the host occurs within the homologous core of attp and attb sites. attp occurs at the single stranded 5’ and 3’ ends of the linear bacteriophage genome. when a bacteriophage genome enters the cell it circularizes and a double stranded attp is formed. the enzymatic integration occurs without loss to the bacteriophage or bacterial genome. we tested sc on a bacterial strain known to harbor an integrated bacteriophage genome. the genome of a salmonella strain carrying a prophage was sequenced by ion torrent using single end reads. thus, only blastn alignments were utilized for sc input. sc identified a total of chimeric reads and reported one integration with ten possible combinations of 5’ and 3’ ambiguous regions. genome coverage of the chimeric reads associated with the 5’ and 3’ ambiguous regions illustrates the integration prediction. the remaining reads are individual chimeras and are likely artifacts of the sequencing process based on lack of chimeric junction coverage. prinseq identified no chimeric reads with an average quality less than the 5’ chimeras were divided between two ambiguous regions and the 3’ chimeric reads were divided into five ambiguous regions. each ambiguous region had single nucleotide variations which likely resulted from sequencing error neighboring the chimeric junction. the integration with the most reads for the 5’ and 3’ chimera was selected from sc output for further analysis. the selected integration showed that the nucleotide compositions of the 5’ and 3’ chimera ambiguous regions were the same. further, the chimeras occurred at the 5’ and 3’ ends of the bacteriophage genome. analysis of the bacterial indices indicated that the bacterial genome remained intact. sc predicted the identical integration as was determined by wgs dna sequencing.figure summonchimera prediction of a salmonella bacteriophage integration. these are segments taken from the aligned region of a 5’ and 3’ chimera. blue is bacterial sequence, red is bacteriophage sequence, and green is homology to both. the numbers indicate the reference sequence nucleotide position. the ambiguous region is nucleotides long. the blue and red arrows show the 5’ to 3’ direction of the host and virus reference strands. the number of reads at each nucleotide position in the 5’ and 3’ chimeras for host and virus are shown. the ambiguous region of the chimera is counted twice due to the unknown nature of where the integration occurs. hepatocellular carcinoma hbv integration in contrast to integration by lysogenic bacteriophage, the integration of dna tumor virus genomes occur at random positions with respect to both the viral and host genomes. for example, when hbv is integrated in hcc there is a loss of both viral and human dna. experimentally determined hbv integrations show no nucleotide specific human integration site. to test the ability of sc to detect hbv integrations we examined two hcc datasets, t and t for which viral chimeras had been pcr verified. this data has been used to verify other chimera detection software. from these datasets sc identified chimeric reads including all reported pcr verified chimeras. prinseq was run on all detected chimeric reads and identified none with average quality less than additionally, sc identified a complete integration event in t and t by associating pcr verified chimeras.figure summonchimera detected recurrent hbv integration in hcc. the human chromosome region with the htert gene from the ucsc genome browser view of the hg is shown. htert is a known recurrent hbv integration site. the red vertical lines indicate the positions of hbv integration, as predicted by sc, that are illustrated in b and c. the t hcc sample hbv integration has been shown to be reference parallel. the red arrow points 5’ to 3’ on the virus reference strand and the numbers indicate the determined viral nucleotides included, keeping in mind hbvs circular genome. green blocks represent the annotated gene coding regions within the integration. the t hcc sample has an hbv integration that places the viral genes in a transcription parallel orientation relative to htert. both integrations detected by sc occurred within the htert gene, a known recurrent integration site for hbv in liver cancer. host positions of these integrations are similar but, their viral genetic content is different. the integration in sample t is reference parallel and integrates the hbv non-coding strand into the coding strand of the htert gene. the hbv integration in sample t is reference anti-parallel, incorporating the hbv coding region into the htert coding strand. thus while occurring in the same host gene these integrations have different biological implications. this result illustrates that the identification of the complete integration event by sc provides additional genetic information beyond current chimeric detection software. discussion large ngs datasets containing genomic dna or rna sequences provides an opportunity to detect and examine viral integration events. the precise mapping of both virus-host junctions in individual tumor samples allows determination of the host gene in which integration occurs, the estimation of viral gene expression, the presence of virus-host chimeric mrnas, and whether host and viral transcription have diverged from their known form. sc was developed to map viral integration events from ngs. sc is designed to keep all analysis under user control and for usability with any ngs dataset. all chimeras and integrations detected by sc are output into a tab separated text document. provided fields include ambiguous host and virus positions, aligned virus and host region, chimera strand orientation, organism order, host sequence name, and the nucleotide distance between the chimeras of a predicted integration. with the provided junctions for the 5’ and 3’ chimeras and strand orientation the lost portion of viral and host genome can be inferred. sc predictions were validated against a variety of experimentally pertinent data. sc detected 100% of the verified chimeras from a hcc dataset and predicted two previously unknown integrations within the hcc genomes. current reports suggest that an integration event may either activate or suppress neighboring host and viral genes. since sc determines the integrated viral genome and lost human genome, this information can be utilized to predict effects on host and viral transcripts. ngs generates artificial chimeras that could be falsely interpreted as an integration by sc. assuming artifact chimeras are generated at random the number of reads covering the same junction will be low relative to authentic chimeric junctions. we developed a model to predict the probability of false-positive integration calls by sc. the equation and calculations for the salmonella and hcc datasets are supplied. chimeric junctions and artifact chimeras are manually differentiated by read coverage. it is incumbent upon the user to determine the authenticity of chimeras and precision of called integrations based on the context of their experiment. recent studies of hpv integrations in cervical cancer have shown a positional correlation to genomic instability. these genomic instability events included copy-number variations, translocations, and inversions. this genomic instability correlation implicates a potential problem for integration inference. when a viral integration occurs near any of these events, then the chimeric cases in sc are not sufficient for integration prediction. further work is required to include genomic instability mapping into sc in order to describe integrations within these regions.figure summonchimera undetectable integrations. viral integrations in cancer have been correlated with chromosomal inversion, translocation, and copy-number variation events. when integrations occur as illustrated, summonchimera cannot predict the genetic content. arrows point from 5’ to 3’ along the reference strand of each organism chromosome. segment labels indicate a portion of the host genome. in the inversion case one of the chimeras will be determined with opposite organism order and strand orientation, thus eluding integration prediction. in case of translocation, the integration can occur between two separate host chromosomes. the two chimeras will be detected on separate chromosomes and no integration would be predicted. finally, the copy-number variation case, if the integration occurs between the duplicated segments the chimeras will be detected in the wrong chromosomal positions and the 3’ chimera will appear to be upstream of the 5’ chimera on the chromosome.  current viral integration software only detects and reports the position of individual chimeras. sc reports chimera organism order, strand orientation and associates appropriate chimeras to accurately infer the genetic content of viral integrations. with the identification of the strand orientation and lost portion of host and viral genome from an integration event, the molecular and genetic consequences of integration can be predicted. availability and requirements project name: summonchimera project home page:http://pipaslab.webfactional.com/wp/wp-content/uploads/2014/01/summonchimera.tar.gz operating system: linux, mac os x, windows programming language: perl other requirements: blast + . and bowtie.0 license: gnu restrictions for use by non-academics: none additional files additional file 1: read counts for figure  1 computational subtraction pipeline. description of dataset: read counts for aligned and unaligned at each step of the pipeline. additional file 2: summonchimera output for salmonella. description of dataset: all reported integrations and orphan chimeras for the salmonella data. additional file 3: summonchimera output for hepatocellular carcinoma t description of dataset: all reported integrations and orphan chimeras for the t data. additional file 4: summonchimera output for hepatocellular carcinoma t description of dataset: all reported integrations and orphan chimeras for the t data. additional file 5: calculations of false-positive integration calls by summonchimera. description of dataset: contains developed equation and calculations for all three datasets. additional file 6: summonchimera. description of dataset: contains summonchimera source code and test data. abbreviations hpvhuman papilloma virus hbvhepatitis b virus scsummonchimera hcchepatocellular carcinoma wgswhole genome sequencing ngsnext generation sequencing competing interests the authors declare that they have no competing interests. authors’ contributions jmp had the initial idea to associate two chimeras to identify integration and identified testable cases. jpk developed the method and code. jpk wrote the manuscript and jmp edited. both authors read and approved the final manuscript. 
without knowing the chimeric boundaries of an integration, its genetic content cannot be determined. these methods identify individual viral-host junctions but do not associate chimeric pairs as an integration event. summonchimera is a perl program that associates chimera pairs to infer the complete viral genomic integration event to the nucleotide level within single or paired-end ngs data. further, summonchimera integration inference precisely predicted bacteriophage integration. viral integration into a host genome is defined by two chimeric junctions that join viral and host dna. recently, computational tools have been developed that utilize ngs data to detect chimeric junctions. furthermore, summonchimera predicted integrations from experimentally verified hepatitis b virus chimeras within a paired-end whole genome sequencing hepatocellular carcinoma tumor database.
 viral integration into a host genome is defined by two chimeric junctions that join viral and host dna. recently, computational tools have been developed that utilize ngs data to detect chimeric junctions. these methods identify individual viral-host junctions but do not associate chimeric pairs as an integration event. without knowing the chimeric boundaries of an integration, its genetic content cannot be determined.  summonchimera is a perl program that associates chimera pairs to infer the complete viral genomic integration event to the nucleotide level within single or paired-end ngs data. summonchimera integration prediction was verified on a set of single-end iontorrent reads from a purified salmonella bacterium with an integrated bacteriophage. furthermore, summonchimera predicted integrations from experimentally verified hepatitis b virus chimeras within a paired-end whole genome sequencing hepatocellular carcinoma tumor database.  summonchimera identified all experimentally verified chimeras detected by current computational methods. further, summonchimera integration inference precisely predicted bacteriophage integration. the application of summonchimera to cancer ngs accurately identifies deletion of host and viral sequence during integration. the precise nucleotide determination of an integration allows prediction of viral and cellular gene transcription patterns. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords integrationviruscancerissue-copyright-statement© the author 2014 