Multi-Gene Detection and Identification of Mosquito-Borne RNA Viruses Using an Oligonucleotide Microarray

Background

Arthropod-borne viruses are important emerging pathogens world-wide. Viruses transmitted by mosquitoes, such as dengue, yellow fever, and Japanese encephalitis viruses, infect hundreds of millions of people and animals each year. Global surveillance of these viruses in mosquito vectors using molecular based assays is critical for prevention and control of the associated diseases. Here, we report an oligonucleotide DNA microarray design, termed ArboChip5.1, for multi-gene detection and identification of mosquito-borne RNA viruses from the genera Flavivirus (family Flaviviridae), Alphavirus (Togaviridae), Orthobunyavirus (Bunyaviridae), and Phlebovirus (Bunyaviridae).

Methodology/Principal Findings

The assay utilizes targeted PCR amplification of three genes from each virus genus for electrochemical detection on a portable, field-tested microarray platform. Fifty-two viruses propagated in cell-culture were used to evaluate the specificity of the PCR primer sets and the ArboChip5.1 microarray capture probes. The microarray detected all of the tested viruses and differentiated between many closely related viruses such as members of the dengue, Japanese encephalitis, and Semliki Forest virus clades. Laboratory infected mosquitoes were used to simulate field samples and to determine the limits of detection. Additionally, we identified dengue virus type 3, Japanese encephalitis virus, Tembusu virus, Culex flavivirus, and a Quang Binh-like virus from mosquitoes collected in Thailand in 2011 and 2012.

Conclusions/Significance

We demonstrated that the described assay can be utilized in a comprehensive field surveillance program by the broad-range amplification and specific identification of arboviruses from infected mosquitoes. Furthermore, the microarray platform can be deployed in the field and viral RNA extraction to data analysis can occur in as little as 12 h. The information derived from the ArboChip5.1 microarray can help to establish public health priorities, detect disease outbreaks, and evaluate control programs.     

Genomic Characterization of Novel Circular ssDNA Viruses from Insectivorous Bats in Southern Brazil

Circoviruses are highly prevalent porcine and avian pathogens. In recent years, novel circular ssDNA genomes have recently been detected in a variety of fecal and environmental samples using deep sequencing approaches. In this study the identification of genomes of novel circoviruses and cycloviruses in feces of insectivorous bats is reported. Pan-reactive primers were used targeting the conserved rep region of circoviruses and cycloviruses to screen DNA bat fecal samples. Using this approach, partial rep sequences were detected which formed five phylogenetic groups distributed among the Circovirus and the recently proposed Cyclovirus genera of the Circoviridae. Further analysis using inverse PCR and Sanger sequencing led to the characterization of four new putative members of the family Circoviridae with genome size ranging from 1,608 to 1,790 nt, two inversely arranged ORFs, and canonical nonamer sequences atop a stem loop.     

Profiling of the bacteria responsible for pyogenic liver abscess by 16S rRNA gene pyrosequencing

Pyogenic liver abscess (PLA) is a severe disease with considerable mortality and is often polymicrobial. Understanding the pathogens that cause PLA is the basis for PLA treatment. Here, we profiled the bacterial composition in PLA fluid by pyrosequencing the 16S ribosomal RNA (rRNA) gene based on next-generation sequencing (NGS) technology to identify etiological agents of PLA and to provide information of their 16S rRNA sequences for application to DNA-based techniques in the hospital. Twenty patients with PLA who underwent percutaneous catheter drainage, abscess culture, and blood culture for isolates were included. Genomic DNAs from abscess fluids were subjected to polymerase chain reaction and pyrosequencing of the 16S rRNA gene with a 454 GS Junior System. The abscess and blood cultures were positive in nine (45%) and four (20%) patients, respectively. Pyrosequencing of 16S rRNA gene showed that 90% of the PLA fluid samples contained single or multiple genera of known bacteria such as Klebsiella, Fusobacterium, Streptococcus, Bacteroides, Prevotella, Peptostreptococcus, unassigned Enterobacteriaceae, and Dialister. Klebsiella was predominantly found in the PLA fluid samples. All samples that carried unassigned bacteria had 26.8% reads on average. We demonstrated that the occurrence of PLA was associated with eight known bacterial genera as well as unassigned bacteria and that 16S rRNA gene sequencing was more useful than conventional culture methods for accurate identification of bacterial pathogens from PLA.     

Carrot yellow leaf virus Is Associated with Carrot Internal Necrosis

Internal necrosis of carrot has been observed in UK carrots for at least 10 years, and has been anecdotally linked to virus infection. In the 2009 growing season some growers had up to 10% of yield with these symptoms. Traditional diagnostic methods are targeted towards specific pathogens. By using a metagenomic approach with high throughput sequencing technology, other, as yet unidentified causes of root necrosis were investigated. Additionally a statistical analysis has shown which viruses are most closely associated with disease symptoms. Carrot samples were collected from a crop exhibiting root necrosis (102 Affected: 99 Unaffected) and tested for the presence of the established carrot viruses: Carrot red leaf virus (CtRLV), Carrot mottle virus (CMoV), Carrot red leaf associated viral RNA (CtRLVaRNA) and Parsnip yellow fleck virus (PYFV). The presence of these viruses was not associated with symptomatic carrot roots either as single viruses or in combinations. A sub-sample of carrots of mixed symptom status was subjected to MiSeq sequencing. The results from these tests suggested Carrot yellow leaf virus (CYLV) was associated with symptomatic roots. Additionally a novel Torradovirus, a novel Closterovirus and two novel Betaflexiviradae related plant viruses were detected. A specific diagnostic test was designed for CYLV. Of the 102 affected carrots, 98% were positive for CYLV compared to 22% of the unaffected carrots. From these data we conclude that although we have yet to practically demonstrate a causal link, CYLV appears to be strongly associated with the presence of necrosis of carrots.     

Viral population analysis and minority-variant detection using short read next-generation sequencing

RNA viruses within infected individuals exist as a population of evolutionary-related variants. Owing to evolutionary change affecting the constitution of this population, the frequency and/or occurrence of individual viral variants can show marked or subtle fluctuations. Since the development of massively parallel sequencing platforms, such viral populations can now be investigated to unprecedented resolution. A critical problem with such analyses is the presence of sequencing-related errors that obscure the identification of true biological variants present at low frequency. Here, we report the development and assessment of the Quality Assessment of Short Read (QUASR) Pipeline (http://sourceforge.net/projects/quasr) specific for virus genome short read analysis that minimizes sequencing errors from multiple deep-sequencing platforms, and enables post-mapping analysis of the minority variants within the viral population. QUASR significantly reduces the error-related noise in deep-sequencing datasets, resulting in increased mapping accuracy and reduction of erroneous mutations. Using QUASR, we have determined influenza virus genome dynamics in sequential samples from an in vitro evolution of 2009 pandemic H1N1 (A/H1N1/09) influenza from samples sequenced on both the Roche 454 GSFLX and Illumina GAIIx platforms. Importantly, concordance between the 454 and Illumina sequencing allowed unambiguous minority-variant detection and accurate determination of virus population turnover in vitro.     

Arbovirus and insect‐specific virus discovery in Kenya by novel six genera multiplex high‐resolution melting analysis

A broad diversity of arthropod‐borne viruses (arboviruses) of global health concern are endemic to East Africa, yet most surveillance efforts are limited to just a few key viral pathogens. Additionally, estimates of arbovirus diversity in the tropics are likely to be underestimated as their discovery has lagged significantly over past decades due to limitations in fast and sensitive arbovirus identification methods. Here, we developed a nearly pan‐arbovirus detection assay that uses high‐resolution melting (HRM) analysis of RT–PCR products from highly multiplexed assays to differentiate broad diversities of arboviruses. We differentiated 15 viral culture controls and seven additional synthetic viral DNA sequence controls, within Flavivirus, Alphavirus, Nairovirus, Phlebovirus, Orthobunyavirus and Thogotovirus genera. Among Bunyamwera, sindbis, dengue and Thogoto virus serial dilutions, detection by multiplex RT–PCR‐HRM was comparable to the gold standard Vero cell plaque assays. We applied our low‐cost method for enhanced broad‐range pathogen surveillance from mosquito samples collected in Kenya and identified diverse insect‐specific viruses, including a new clade in anopheline mosquitoes, and Wesselsbron virus, an arbovirus that can cause viral haemorrhagic fever in humans and has not previously been isolated in Kenya, in Culex spp. and Anopheles coustani mosquitoes. Our findings demonstrate how multiplex RT–PCR‐HRM can identify novel viral diversities and potential disease threats that may not be included in pathogen detection panels of routine surveillance efforts. This approach can be adapted to other pathogens to enhance disease surveillance and pathogen discovery efforts, as well as the study of pathogen diversity and viral evolutionary ecology.     

Virome Profiling of Bats from Myanmar by Metagenomic Analysis of Tissue Samples Reveals More Novel Mammalian Viruses

Bats are reservoir animals harboring many important pathogenic viruses and with the capability of transmitting these to humans and other animals. To establish an effective surveillance to monitor transboundary spread of bat viruses between Myanmar and China, complete organs from the thorax and abdomen from 853 bats of six species from two Myanmar counties close to Yunnan province, China, were collected and tested for their virome through metagenomics by Solexa sequencing and bioinformatic analysis. In total, 3,742,314 reads of 114 bases were generated, and over 86% were assembled into 1,649,512 contigs with an average length of 114 bp, of which 26,698 (2%) contigs were recognizable viral sequences belonging to 24 viral families. Of the viral contigs 45% (12,086/26,698) were related to vertebrate viruses, 28% (7,443/26,698) to insect viruses, 27% (7,074/26,698) to phages and 95 contigs to plant viruses. The metagenomic results were confirmed by PCR of selected viruses in all bat samples followed by phylogenetic analysis, which has led to the discovery of some novel bat viruses of the genera Mamastrovirus, Bocavirus, Circovirus, Iflavirus and Orthohepadnavirus and to their prevalence rates in two bat species. In conclusion, the present study aims to present the bat virome in Myanmar, and the results obtained further expand the spectrum of viruses harbored by bats.     

Bacterial Pathogens and Community Composition in Advanced Sewage Treatment Systems Revealed by Metagenomics Analysis Based on High-Throughput Sequencing

This study used 454 pyrosequencing, Illumina high-throughput sequencing and metagenomic analysis to investigate bacterial pathogens and their potential virulence in a sewage treatment plant (STP) applying both conventional and advanced treatment processes. Pyrosequencing and Illumina sequencing consistently demonstrated that Arcobacter genus occupied over 43.42% of total abundance of potential pathogens in the STP. At species level, potential pathogens Arcobacter butzleri, Aeromonas hydrophila and Klebsiella pneumonia dominated in raw sewage, which was also confirmed by quantitative real time PCR. Illumina sequencing also revealed prevalence of various types of pathogenicity islands and virulence proteins in the STP. Most of the potential pathogens and virulence factors were eliminated in the STP, and the removal efficiency mainly depended on oxidation ditch. Compared with sand filtration, magnetic resin seemed to have higher removals in most of the potential pathogens and virulence factors. However, presence of the residual A. butzleri in the final effluent still deserves more concerns. The findings indicate that sewage acts as an important source of environmental pathogens, but STPs can effectively control their spread in the environment. Joint use of the high-throughput sequencing technologies is considered a reliable method for deep and comprehensive overview of environmental bacterial virulence.     

The Genome Sequence of Lone Star Virus,a Highly Divergent Bunyavirus Found in the Amblyomma americanum Tick

Viruses in the family Bunyaviridae infect a wide range of plant, insect, and animal hosts. Tick-borne bunyaviruses in the Phlebovirus genus, including Severe Fever with Thrombocytopenia Syndrome virus (SFTSV) in China, Heartland virus (HRTV) in the United States, and Bhanja virus in Eurasia and Africa have been associated with acute febrile illness in humans. Here we sought to characterize the growth characteristics and genome of Lone Star virus (LSV), an unclassified bunyavirus originally isolated from the lone star tick Amblyomma americanum. LSV was able to infect both human (HeLa) and monkey (Vero) cells. Cytopathic effects were seen within 72 h in both cell lines; vacuolization was observed in infected Vero, but not HeLa, cells. Viral culture supernatants were examined by unbiased deep sequencing and analysis using an in-house developed rapid computational pipeline for viral discovery, which definitively identified LSV as a phlebovirus. De novo assembly of the full genome revealed that LSV is highly divergent, sharing <61% overall amino acid identity with any other bunyavirus. Despite this sequence diversity, LSV was found by phylogenetic analysis to be part of a well-supported clade that includes members of the Bhanja group viruses, which are most closely related to SFSTV/HRTV. The genome sequencing of LSV is a critical first step in developing diagnostic tools to determine the risk of arbovirus transmission by A. americanum, a tick of growing importance given its expanding geographic range and competence as a disease vector. This study also underscores the power of deep sequencing analysis in rapidly identifying and sequencing the genomes of viruses of potential clinical and public health significance.     

Prospecting for Novel Biocatalysts in a Soil Metagenome

The metagenomes of complex microbial communities are rich sources of novel biocatalysts. We exploited the metagenome of a mixed microbial population for isolation of more than 15 different genes encoding novel biocatalysts by using a combined cultivation and direct cloning strategy. A 16S rRNA sequence analysis revealed the presence of hitherto uncultured microbes closely related to the genera Pseudomonas, Agrobacterium, Xanthomonas, Microbulbifer, and Janthinobacterium. Total genomic DNA from this bacterial community was used to construct cosmid DNA libraries, which were functionally searched for novel enzymes of biotechnological value. Our searches in combination with cosmid sequencing resulted in identification of four clones encoding 12 putative agarase genes, most of which were organized in clusters consisting of two or three genes. Interestingly, nine of these agarase genes probably originated from gene duplications. Furthermore, we identified by DNA sequencing several other biocatalyst-encoding genes, including genes encoding a putative stereoselective amidase (amiA), two cellulases (gnuB and uvs080), an α-amylase (amyA), a 1,4-α-glucan branching enzyme (amyB), and two pectate lyases (pelA and uvs119). Also, a conserved cluster of two lipase genes was identified, which was linked to genes encoding a type I secretion system. The novel gene aguB was overexpressed in Escherichia coli, and the enzyme activities were determined. Finally, we describe more than 162 kb of DNA sequence that provides a strong platform for further characterization of this microbial consortium.     

BAC-Pool Sequencing and Analysis of Large Segments of A12 and D12 Homoeologous Chromosomes in Upland Cotton

Although new and emerging next-generation sequencing (NGS) technologies have reduced sequencing costs significantly, much work remains to implement them for de novo sequencing of complex and highly repetitive genomes such as the tetraploid genome of Upland cotton (Gossypium hirsutum L.). Herein we report the results from implementing a novel, hybrid Sanger/454-based BAC-pool sequencing strategy using minimum tiling path (MTP) BACs from Ctg-3301 and Ctg-465, two large genomic segments in A12 and D12 homoeologous chromosomes (Ctg). To enable generation of longer contig sequences in assembly, we implemented a hybrid assembly method to process ~35x data from 454 technology and 2.8-3x data from Sanger method. Hybrid assemblies offered higher sequence coverage and better sequence assemblies. Homology studies revealed the presence of retrotransposon regions like Copia and Gypsy elements in these contigs and also helped in identifying new genomic SSRs. Unigenes were anchored to the sequences in Ctg-3301 and Ctg-465 to support the physical map. Gene density, gene structure and protein sequence information derived from protein prediction programs were used to obtain the functional annotation of these genes. Comparative analysis of both contigs with Arabidopsis genome exhibited synteny and microcollinearity with a conserved gene order in both genomes. This study provides insight about use of MTP-based BAC-pool sequencing approach for sequencing complex polyploid genomes with limited constraints in generating better sequence assemblies to build reference scaffold sequences. Combining the utilities of MTP-based BAC-pool sequencing with current longer and short read NGS technologies in multiplexed format would provide a new direction to cost-effectively and precisely sequence complex plant genomes.     

Targeted next-generation sequencing of DNA regions proximal to a conserved GXGXXG signaling motif enables systematic discovery of tyrosine kinase fusions in cancer

Tyrosine kinase (TK) fusions are attractive drug targets in cancers. However, rapid identification of these lesions has been hampered by experimental limitations. Our in silico analysis of known cancer-derived TK fusions revealed that most breakpoints occur within a defined region upstream of a conserved GXGXXG kinase motif. We therefore designed a novel DNA-based targeted sequencing approach to screen systematically for fusions within the 90 human TKs; it should detect 92% of known TK fusions. We deliberately paired ‘in-solution’ DNA capture with 454 sequencing to minimize starting material requirements, take advantage of long sequence reads, and facilitate mapping of fusions. To validate this platform, we analyzed genomic DNA from thyroid cancer cells (TPC-1) and leukemia cells (KG-1) with fusions known only at the mRNA level. We readily identified for the first time the genomic fusion sequences of CCDC6-RET in TPC-1 cells and FGFR1OP2-FGFR1 in KG-1 cells. These data demonstrate the feasibility of this approach to identify TK fusions across multiple human cancers in a high-throughput, unbiased manner. This method is distinct from other similar efforts, because it focuses specifically on targets with therapeutic potential, uses only 1.5 µg of DNA, and circumvents the need for complex computational sequence analysis.     

Enhanced De Novo Assembly of High Throughput Pyrosequencing Data Using Whole Genome Mapping

Despite major advances in next-generation sequencing, assembly of sequencing data, especially data from novel microorganisms or re-emerging pathogens, remains constrained by the lack of suitable reference sequences. De novo assembly is the best approach to achieve an accurate finished sequence, but multiple sequencing platforms or paired-end libraries are often required to achieve full genome coverage. In this study, we demonstrated a method to assemble complete bacterial genome sequences by integrating shotgun Roche 454 pyrosequencing with optical whole genome mapping (WGM). The whole genome restriction map (WGRM) was used as the reference to scaffold de novo assembled sequence contigs through a stepwise process. Large de novo contigs were placed in the correct order and orientation through alignment to the WGRM. De novo contigs that were not aligned to WGRM were merged into scaffolds using contig branching structure information. These extended scaffolds were then aligned to the WGRM to identify the overlaps to be eliminated and the gaps and mismatches to be resolved with unused contigs. The process was repeated until a sequence with full coverage and alignment with the whole genome map was achieved. Using this method we were able to achieved 100% WGRM coverage without a paired-end library. We assembled complete sequences for three distinct genetic components of a clinical isolate of Providencia stuartii: a bacterial chromosome, a novel bla _NDM-1 plasmid, and a novel bacteriophage, without separately purifying them to homogeneity.     

Full Genome Characterization of the Culicoides-Borne Marsupial Orbiviruses: Wallal Virus,Mudjinbarry Virus and Warrego Viruses

Viruses belonging to the species Wallal virus and Warrego virus of the genus Orbivirus were identified as causative agents of blindness in marsupials in Australia during 1994/5. Recent comparisons of nucleotide (nt) and amino acid (aa) sequences have provided a basis for the grouping and classification of orbivirus isolates. However, full-genome sequence data are not available for representatives of all Orbivirus species. We report full-genome sequence data for three additional orbiviruses: Wallal virus (WALV); Mudjinabarry virus (MUDV) and Warrego virus (WARV). Comparisons of conserved polymerase (Pol), sub-core-shell ‘T2’ and core-surface ‘T13’ proteins show that these viruses group with other Culicoides borne orbiviruses, clustering with Eubenangee virus (EUBV), another orbivirus infecting marsupials. WARV shares <70% aa identity in all three conserved proteins (Pol, T2 and T13) with other orbiviruses, consistent with its classification within a distinct Orbivirus species. Although WALV and MUDV share <72.86%/67.93% aa/nt identity with other orbiviruses in Pol, T2 and T13, they share >99%/90% aa/nt identities with each other (consistent with membership of the same virus species - Wallal virus). However, WALV and MUDV share <68% aa identity in their larger outer capsid protein VP2(OC1), consistent with membership of different serotypes within the species - WALV-1 and WALV-2 respectively.     

Sequence analysis of the human virome in febrile and afebrile children

Unexplained fever (UF) is a common problem in children under 3 years old. Although virus infection is suspected to be the cause of most of these fevers, a comprehensive analysis of viruses in samples from children with fever and healthy controls is important for establishing a relationship between viruses and UF. We used unbiased, deep sequencing to analyze 176 nasopharyngeal swabs (NP) and plasma samples from children with UF and afebrile controls, generating an average of 4.6 million sequences per sample. An analysis pipeline was developed to detect viral sequences, which resulted in the identification of sequences from 25 viral genera. These genera included expected pathogens, such as adenoviruses, enteroviruses, and roseoloviruses, plus viruses with unknown pathogenicity. Viruses that were unexpected in NP and plasma samples, such as the astrovirus MLB-2, were also detected. Sequencing allowed identification of virus subtype for some viruses, including roseoloviruses. Highly sensitive PCR assays detected low levels of viruses that were not detected in approximately 5 million sequences, but greater sequencing depth improved sensitivity. On average NP and plasma samples from febrile children contained 1.5- to 5-fold more viral sequences, respectively, than samples from afebrile children. Samples from febrile children contained a broader range of viral genera and contained multiple viral genera more frequently than samples from children without fever. Differences between febrile and afebrile groups were most striking in the plasma samples, where detection of viral sequence may be associated with a disseminated infection. These data indicate that virus infection is associated with UF. Further studies are important in order to establish the range of viral pathogens associated with fever and to understand of the role of viral infection in fever. Ultimately these studies may improve the medical treatment of children with UF by helping avoid antibiotic therapy for children with viral infections.