De novo identification of viral pathogens from cell culture hologenomes
© Gore et al; licensee BioMed Central Ltd. 2011
Received: 27 August 2011
Accepted: 6 January 2012
Published: 6 January 2012
Fast, specific identification and surveillance of pathogens is the cornerstone of any outbreak response system, especially in the case of emerging infectious diseases and viral epidemics. This process is generally tedious and time-consuming thus making it ineffective in traditional settings. The added complexity in these situations is the non-availability of pure isolates of pathogens as they are present as mixed genomes or hologenomes. Next-generation sequencing approaches offer an attractive solution in this scenario as it provides adequate depth of sequencing at fast and affordable costs, apart from making it possible to decipher complex interactions between genomes at a scale that was not possible before. The widespread application of next-generation sequencing in this field has been limited by the non-availability of an efficient computational pipeline to systematically analyze data to delineate pathogen genomes from mixed population of genomes or hologenomes.
We applied next-generation sequencing on a sample containing mixed population of genomes from an epidemic with appropriate processing and enrichment. The data was analyzed using an extensive computational pipeline involving mapping to reference genome sets and de-novo assembly. In depth analysis of the data generated revealed the presence of sequences corresponding to Japanese encephalitis virus. The genome of the virus was also independently de-novo assembled. The presence of the virus was in addition, verified using standard molecular biology techniques.
Our approach can accurately identify causative pathogens from cell culture hologenome samples containing mixed population of genomes and in principle can be applied to patient hologenome samples without any background information. This methodology could be widely applied to identify and isolate pathogen genomes and understand their genomic variability during outbreaks.
KeywordsEpidemics Mixed Population Genomes Hologenome De novo assembly Japanese encephalitis Next generation sequencing
Viral pathogens have been a major cause of epidemics worldwide [1–3]. The regular surveillance system in many cases of viral epidemics is ineffective to directly identify the virus unless directed by telltale clinical features and clinical complications [4, 5]. In many cases the disease causing pathogen is not identified , and contributes to the inadequate and inappropriate management of the condition. Many of the disease outbreaks are also caused by changes in the environmental niche , and at least in some cases cause the emergence of new pathogens [6, 8–11]. It is estimated that at least 33 new pathogens have emerged during the last three decades . Identification of causative organisms during viral outbreak is a problem, primarily due to the tedious methods involving isolation and culturing of the pathogen [13, 14]. The relatively low concentration of the genetic material, especially in the case of RNA viruses and the heavy interference due to the host cell nucleic acid and other metagenomes makes sequencing based approaches ineffective, unless done at higher depths [15, 16].
The availability of next-generation sequencing (NGS) technology has enabled the scale and ease of addressing biological questions on a genomics perspective [17, 18]. The throughput of sequencing enables deep sequencing of nucleic acids, adequate to provide for enough reads of the pathogen, even while the interference of the host genetic material is very high. Metagenomics has been one of the major applications of NGS technology for understanding the composition and dynamics of mixed population of organisms . The field has now emerged to a vibrant area of genomics trying to understand a large spectrum of environmental niches right from human body in both disease as well as healthy states to natural geographical niches [20, 21]. Although NGS technology has been successfully used for addressing a large diversity of biological questions, its application to address questions pertaining to bio-surveillance and emerging infectious diseases on a large scale has been limited , in spite of the unprecedented opportunity provided by the scale and speed of operations .
Hologenome, a term borrowed from evolutionary biology is defined as the sum of the genetic information of the host and its microbiota . Hologenomics is an emerging field in genomics which deals with mixed population of genomes, as in the case of interacting populations in host-pathogen and commensals. Hologenome differs from the widely popular term Metagenome, which involves the study of communities of microbes directly in their natural environments .
Cultured population of viruses co-exists and interacts intricately with their host genomes . Although this presents a technical challenge in isolating individual genomes from mixed populations, it offers enormous possibility to understand the interactions and dynamics between the genomes in real-time.
Here we report the sequencing and analysis methodology, involving computational algorithms for reference mapping and de novo sequence assembly to accurately identify viral pathogens from mixed populations of genomes. The pipeline relies on the specificity of sequence mappings and the differential distribution of the mapped reads across genomes. As a proof of concept, we applied the methodology on a cell culture hologenome consisting of human, bacterial and viral genomes, and could specifically identify the viral pathogen. This methodology could potentially be applied for rapid and specific identification of viral pathogens during epidemic outbreaks.
Sample collection and RNA isolation
The sample was collected and isolated during an epidemic of acute viral encephalitis from an anonymous patient suffering from fever and acute encephalitis from Baba Raghav Das (BRD) Medical College and Nehru Hospital, Gorakhpur, India. Sample was procured and processed as per ethical procedures laid down by BRD Medical College, Gorakhpur, India and National Institute of Virology, Pune, India. Using standard virus isolation protocols samples were inoculated in human Rhabdosarcoma (RD) and Baby Hamster kidney (BHK) cell lines for virus isolation [25, 26]. Cells were observed for cytopathological effects (CPE) and passaged three times. The cell culture supernatant was filtered using 0.22 μm Millipore filters for every passage. RNA was isolated using Qiagen (QiaAMP viral RNA minikit) kit as per manufacturer's instructions. The RNA was eluted in 60 μlitre of AVE buffer.
Library preparation, sequencing and genome assembly
The RNA library was prepared according to the manufacturer's instructions using RNA Sample-prep kit (Illumina Inc, USA) for sequencing on Illumina sequencing platform. Two microgram of total RNA was fragmented using divalent cation. Cleaved RNA was converted to cDNA using reverse transcriptase (SSRT-II Invitrogen) and random primers. The fragments were further subjected to second strand cDNA synthesis using DNA polymerase as per manufacturer's instructions. End-repairing process followed by A-base addition and adapter ligation was further performed on the cDNA fragments. Approximately 350 base pair products were separated by gel excision and enriched with PCR to create the final library.
Clusters were generated on the flow cell using cBot Paired end cluster generation kit (Illumina Inc, USA) as per manufacturer's instructions. The sequencing runs were performed on Illumina Inc, USA) using 76 × 2 base reads. The sequence-quality files generated was transformed to Sanger quality scores using custom scripts.
Velvet , a popularly used de novo assembly algorithm based on de Bruijn graphs was used for the de novo assembly. The entire read data was partitioned into smaller subsets for analysis. De-novo assembly was attempted on the subsets with different k-mers. The data was compiled and compared using in-house scripts.
Experimental validation of the virus was performed using reverse transcription (RT)-polymerase chain reaction. Specific primer JEV PM1R (reverse primer): 5'-CGGARTCTCCTGCTTCGCTTGG-3' and JEV C1F (forward primer): 5'-GGCAGAAAGCAAAACAAAAGA-3' specific for Japanese encephalitis virus were used. RNA was initially reverse transcribed using JEV PM1R at 50°C using Superscript II reverse transcriptase (Invitrogen, Life Science Technologies). The cDNA was amplified using forward primer JEV C1F and reverse primer JEV PM1R. PCR amplification was carried out by denaturing 94°C for 5 min, followed by 35 cycles of 94°C for 30 sec 58°C for 30 sec, 72°C for 1 min and final extension of 3 min using Taq DNA polymerase (Fermentas).
Results and discussion
Reference mapping is an attractive strategy, as it is fast and frugal on compute and memory requirements, however is limited by the availability of the genome under consideration in the reference dataset. This could limit the widespread application of this methodology in identifying pathogens that have not been sequenced before, as in the case of new and emerging infectious diseases. This limitation could be potentially overcome using a de novo assembly strategy, which does not rely on prior knowledge of the genome sequence. We therefore attempted de novo assembly of the cell culture hologenome using Velvet. Different k-mer values and coverage cut-offs were used for the assembly. The largest contig was assembled at a k-mer of 27 and had a length of 10,758 bases [GenBank: JN644310]. The de novo assembled contig aligned with 99% identity to the Japanese encephalitis virus isolate 014178 genome in NCBI nr database, having a length of 10,976 bases, further confirming the identity of the virus. We also attempted to validate the identity of virus using specific reverse transcriptase PCR (RT-PCR). JE specific primers mapping in the C-prM region amplified approximately 400 bp amplicon as expected [25, 26], substantiating the identity of the virus. The data is represented in Additional file 2.
Furthermore the reads and alignments in the reference assembly that mapped to Japanese encephalitis genome were analyzed in depth to understand the genome organization and genomic variations of the isolate. We identified a total of 209 genomic variations in the genome with high confidence. The complete list of variations, their genomic loci, their type (i.e. synonymous and non-synonymous) and the genes that harbor them are summarized in Additional file 3. Fast and specific identification of viral pathogens would enable adequate and timely response to disease outbreaks. Next generation sequencing technology coupled with efficient computational algorithms could be effectively used for the identification of pathogens in such epidemics.
Here we have used next-generation sequencing approach and generated over 22 million sequence reads from a cell culture hologenome sample consisting of human, bacteria and virus. We also successfully applied a pipeline of computational algorithms including reference mapping and de novo assembly to identify the pathogen as Japanese encephalitis virus. Furthermore we also validated the identity of the virus using RT-PCR techniques. In summary we successfully demonstrate the utility of such approaches in identification of viral pathogens from mixed population of genomes. As cheaper and faster sequencing technology becomes widely available in research labs, this approach would open up immense opportunities in identifying viral pathogens in outbreaks, apart from being used in regular bio-surveillance of infectious diseases and screening for agents used in bioterrorism. This methodology could also find application in special situations like food processing, beverage, and water quality monitoring for identification of food borne and waterborne pathogens respectively.
Availability of supporting data
Raw sequence data is available at the NCBI Short Read Archive with ID: SRX099040. Sequence generated using de novo genome assembly tools is available at Genbank with ID: JN644310.
Authors acknowledge Dr. Bhupesh Taneja and Dr. Chetana Sachidanandan for reviewing the manuscript and Dr. Swati Subodh for discussions. Dr. Naresh Singh and Mr. Vikas Pandey are acknowledged for the technical support. This work was funded by the Council for Scientific and Industrial Research (CSIR), India through EMPOWER Grant EMP0005. AP acknowledges the Senior Research fellowship from CSIR, India. The Sequencing facility is supported through SIP006 and FAC002 Grant from CSIR, India and the Computational analysis was performed at the CSIR centre for in silico Biology at IGIB.
- Chan-Yeung M, Xu RH: SARS: epidemiology. Respirology. 2003, 8 (Suppl): S9-14.PubMedView ArticleGoogle Scholar
- Chua KB, Goh KJ, Wong KT, Kamarulzaman A, Tan PS, Ksiazek TG, et al: Fatal encephalitis due to Nipah virus among pig-farmers in Malaysia. Lancet. 1999, 354: 1257-1259. 10.1016/S0140-6736(99)04299-3.PubMedView ArticleGoogle Scholar
- World Health Organization: Chikungunya. Fact sheet N°327. 2008, [http://www.who.int/mediacentre/factsheets/fs327/en/]Google Scholar
- Petric M, Comanor L, Petti CA: Role of the laboratory in diagnosis of influenza during seasonal epidemics and potential pandemics. J Infect Dis. 2006, 194 (Suppl 2): S98-110.PubMedView ArticleGoogle Scholar
- Hatchette TF, Bastien N, Berry J, Booth TF, Chernesky M, Couillard M, et al: The limitations of point of care testing for pandemic influenza: what clinicians and public health professionals need to know. Can J Public Health. 2009, 100: 204-207.PubMedGoogle Scholar
- Feldmann H: Truly emerging-a new disease caused by a novel virus. N Engl J Med. 2011, 364: 1561-1563. 10.1056/NEJMe1102671.PubMedView ArticleGoogle Scholar
- Patz JA, Epstein PR, Burke TA, Balbus JM: Global climate change and emerging infectious diseases. JAMA. 1996, 275: 217-223. 10.1001/jama.1996.03530270057032.PubMedView ArticleGoogle Scholar
- Murphy FA: The evolution of viruses, the emergence of viral diseases: a synthesis that Martinus Beijerinck might enjoy. Arch Virol Suppl. 1999, 15: 73-85.PubMedGoogle Scholar
- Holmes EC, Drummond AJ: The evolutionary genetics of viral emergence. Curr Top Microbiol Immunol. 2007, 315: 51-66. 10.1007/978-3-540-70962-6_3.PubMedGoogle Scholar
- Fisher D, Hui D, Gao Z, Lee C, Oh MD, Bin C, et al: Pandemic response lessons from influenza H1N1 2009 in Asia. Respirology. 2011, 16: 876-882. 10.1111/j.1440-1843.2011.02003.x.PubMedView ArticleGoogle Scholar
- Fraser C, Donnelly CA, Cauchemez S, Hanage WP, Van K, Hollingsworth TD, et al: Pandemic potential of a strain of influenza A (H1N1): early findings. Sci. 2009, 19 (324): 1557-1561.View ArticleGoogle Scholar
- Heymann DL, Rodier GR: Hot spots in a wired world: WHO surveillance of emerging and re-emerging infectious diseases. Lancet Infect Dis. 2001, 1: 345-353. 10.1016/S1473-3099(01)00148-7.PubMedView ArticleGoogle Scholar
- Call DR: Challenges and opportunities for pathogen detection using DNA microarrays. Crit Rev Microbiol. 2005, 31: 91-99. 10.1080/10408410590921736.PubMedView ArticleGoogle Scholar
- Avarre JC, de LP, Bena G: Hybridization of genomic DNA to microarrays: a challenge for the analysis of environmental samples. J Microbiol Methods. 2007, 69: 242-248. 10.1016/j.mimet.2006.11.007.PubMedView ArticleGoogle Scholar
- Potgieter AC, Page NA, Liebenberg J, Wright IM, Landt O, van Dijk AA: Improved strategies for sequence-independent amplification and sequencing of viral double-stranded RNA genomes. J Gen Virol. 2009, 90: 1423-1432. 10.1099/vir.0.009381-0.PubMedView ArticleGoogle Scholar
- Hoper D, Hoffmann B, Beer M: A comprehensive deep sequencing strategy for full-length genomes of influenza A. PLoS One. 2011, 6: e19075-10.1371/journal.pone.0019075.PubMedPubMed CentralView ArticleGoogle Scholar
- Shendure J, Ji H: Next-generation DNA sequencing. Nat Biotechnol. 2008, 26: 1135-1145. 10.1038/nbt1486.PubMedView ArticleGoogle Scholar
- Mardis ER: The impact of next-generation sequencing technology on genetics. Trends Genet. 2008, 24: 133-141. 10.1016/j.tig.2007.12.007.PubMedView ArticleGoogle Scholar
- National Research Council (US) Committee: The new science of metagenomics: revealing the secrets of our microbial planet. 2007, National Academies Press (US)Google Scholar
- Singh J, Behal A, Singla N, Joshi A, Birbian N, Singh S, et al: Metagenomics: concept, methodology, ecological inference and recent advances. Biotechnol J. 2009, 4: 480-494. 10.1002/biot.200800201.PubMedView ArticleGoogle Scholar
- Tang P, Chiu C: Metagenomics for the discovery of novel human viruses. Future Microbiol. 2010, 5: 177-189. 10.2217/fmb.09.120.PubMedView ArticleGoogle Scholar
- Isakov O, Modai S, Shomron N: Pathogen Detection Using Short-RNA Deep Sequencing Subtraction and Assembly. Bioinformatics. 2011, 27: 2027-2030. 10.1093/bioinformatics/btr349.PubMedPubMed CentralView ArticleGoogle Scholar
- Zilber-Rosenberg I, Rosenberg E: Role of microorganisms in the evolution of animals and plants: the hologenome theory of evolution. FEMS Microbiol Rev. 2008, 32: 723-735. 10.1111/j.1574-6976.2008.00123.x.PubMedView ArticleGoogle Scholar
- Handelsman J: Metagenomics: application of genomics to uncultured microorganisms. Microbiol Mol Biol Rev. 2004, 68: 669-685. 10.1128/MMBR.68.4.669-685.2004.PubMedPubMed CentralView ArticleGoogle Scholar
- Sapkal GN, Bondre VP, Fulmali PV, Patil P, Gopalkrishna V, Dadhania V, et al: Enteroviruses in patients with acute encephalitis, uttar pradesh, India. Emerg Infect Dis. 2009, 15: 295-298. 10.3201/eid1502.080865.PubMedPubMed CentralView ArticleGoogle Scholar
- Fulmali PV, Sapkal GN, Athawale S, Gore MM, Mishra AC, Bondre VP: Introduction of Japanese encephalitis virus genotype I, India. Emerg Infect Dis. 2011, 17: 319-321.PubMedPubMed CentralView ArticleGoogle Scholar
- Li H, Ruan J, Durbin R: Mapping short DNA sequencing reads and calling variants using mapping quality scores. Genome Res. 2008, 18: 1851-1858. 10.1101/gr.078212.108.PubMedPubMed CentralView ArticleGoogle Scholar
- Wheeler DL, Chappey C, Lash AE, Leipe DD, Madden TL, Schuler GD, et al: Database resources of the National Center for Biotechnology Information. Nucleic Acids Res. 2000, 28: 10-14. 10.1093/nar/28.1.10.PubMedPubMed CentralView ArticleGoogle Scholar
- Zerbino DR, Birney E: Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008, 18: 821-829. 10.1101/gr.074492.107.PubMedPubMed CentralView ArticleGoogle Scholar