Report of the second international conference on next generation sequencing for adventitious virus detection in biologics for humans and animals.

The IABS-EU, in association with PROVAXS and Ghent University, hosted the "2nd Conference on Next Generation Sequencing (NGS) for Adventitious Virus Detection in Human and Veterinary Biologics" held on November 13th and 14th 2019, in Ghent, Belgium. The meeting brought together international experts from regulatory agencies, the biotherapeutics and biologics industries, contract research organizations, and academia, with the goal to develop a scientific consensus on the readiness of NGS for detecting adventitious viruses, and on the use of this technology to supplement or replace/substitute the currently used assays. Participants discussed the progress on the standardization and validation of the technical and bioinformatics steps in NGS for characterization and safety evaluation of biologics, including human and animal vaccines. It was concluded that NGS can be used for the detection of a broad range of viruses, including novel viruses, and therefore can complement, supplement or even replace some of the conventional adventitious virus detection assays. Furthermore, the development of reference viral standards, complete and correctly annotated viral databases, and protocols for the validation and follow-up investigations of NGS signals is necessary to enable broader use of NGS. An international collaborative effort, involving regulatory authorities, industry, academia, and other stakeholders is ongoing toward this goal.

Influence of pig gut microbiota on Mycoplasma hyopneumoniae susceptibility.

This study investigated the influence of gut microbiome composition in modulating susceptibility to Mycoplasma hyopneumoniae in pigs. Thirty-two conventional M. hyopneumoniae free piglets were randomly selected from six different litters at 3 weeks of age and were experimentally inoculated with M. hyopneumoniae at 8 weeks of age. Lung lesion scores (LS) were recorded 4 weeks post-inoculation (12 weeks of age) from piglet lungs at necropsy. Fecal bacterial community composition of piglets at 3, 8 and 12 weeks of age were targeted by amplifying the V3-V4 region of the 16S rRNA gene. The LS ranged from 0.3 to 43% with an evident clustering of the scores observed in piglets within litters. There were significant differences in species richness and alpha diversity in fecal microbiomes among piglets within litters at different time points (p < 0.05). The dissimilarity matrices indicated that at 3 weeks of age, the fecal microbiota of piglets was more dissimilar compared to those from 8 to 12 weeks of age. Specific groups of bacteria in the gut that might predict the decreased severity of M. hyopneumoniae associated lesions were identified. The microbial shift at 3 weeks of age was observed to be driven by the increase in abundance of the indicator family, Ruminococcaceae in piglets with low LS (p < 0.05). The taxa, Ruminococcus_2 having the highest richness scores, correlated significantly between litters showing stronger associations with the lowest LS (r = -0.49, p = 0.005). These findings suggest that early life gut microbiota can be a potential determinant for M. hyopneumoniae susceptibility in pigs.

Identification of a Divergent Lineage Porcine Pestivirus in Nursing Piglets with Congenital Tremors and Reproduction of Disease following Experimental Inoculation.

Congenital tremors is a sporadic disease of neonatal pigs characterized by action-related repetitive myoclonus. A majority of outbreaks of congenital tremors have been attributed to an unidentified virus. The objectives of this project were to 1) detect potential pathogen(s) in samples from piglets with congenital tremors and 2) develop an infection model to reproduce disease. Using next-generation sequencing, a divergent lineage pestivirus was detected in piglets with congenital tremors. The virus was originally most closely related to a bat pestivirus but is now more closely related to a recently published novel porcine pestivirus provisionally named atypical porcine pestivirus. A quantitative real-time PCR detected the virus in samples from neonatal piglets with congenital tremors from two separate farms, but not in samples from unaffected piglets from the same farm. To fulfill the second objective, pregnant sows were inoculated with either serum containing the pestivirus or PBS (control) by intravenous and intranasal routes simultaneously with direct inoculation of fetal amniotic vesicles by ultrasound-guided surgical technique. Inoculations were performed at either 45 or 62 days of gestation. All sows inoculated with the novel pestivirus farrowed piglets affected with congenital tremors while PBS-inoculated control piglets were unaffected. Tremor severity for each piglet was scored from videos taken 0, 1 and 2 days post-farrowing. Tremor severity remained relatively constant from 0 to 2 days post-farrowing for a majority of piglets. The prevalence of congenital tremors in pestivirus-inoculated litters ranged from 57% (4 out of 7 affected piglets) to 100% (10 out of 10 affected piglets). The virus was consistently detected by PCR in tissues from piglets with congenital tremors but was not detected in control piglets. Samples positive by PCR in greater than 90% of piglets sampled included brainstem (37 out of 41), mesenteric lymph node (37 out of 41), tracheobronchial lymph node (37 out of 41), and whole blood (19 out of 20). Although the first description of congenital tremors was in 1922, this is the first reported reproduction of congenital tremors following experimental inoculation with a divergent lineage porcine pestivirus. Studies investigating disease mechanism, epidemiology, and diagnostic assay development are needed to better understand the pathophysiology of congenital tremors due to this pestivirus.

Genome sequence and comparative virulence of raccoonpox virus: the first North American poxvirus sequence.

We report here the complete genome sequence of raccoonpox virus (RCNV), a naturally occurring North American poxvirus. This is the first such North American sequence to the best of our knowledge, and the data showed that RCNV forms a new phylogenetic branch between orthopoxviruses and Yoka poxvirus. RCNV shared overall similarity in genome organization with orthopoxviruses, and the proteins in the central conserved region shared approximately 90  % amino acid identity with orthopoxviruses. RCNV proteins shared approximately 81  % amino acid identity with Yokapox virus proteins. RCNV is missing 10 genes normally conserved in orthopoxviruses, most of which are implicated in virulence. These gene deletions may explain the attenuated phenotype of RCNV in mammals. RCNV contained one unique genome region containing approximately 1 kb of DNA sequence that is not present in any reported poxvirus. It contained a unique ORF predicted to encode a protein with a transmembrane domain. RCNV replicates well in mammalian cells, is naturally attenuated and has been shown to be effective as a vaccine vector platform, so we further tested its safety. We showed here that RCNV is substantially more attenuated than even the highly attenuated VACV-A35Del mutant virus in pregnant, nude and severe combined immunodeficient (SCID) mouse models. RCNV was much safer in pregnant mice and was cleared rapidly from tissues, even in immunocompromised animals, whereas the VACV-A35Del mutant retains virulence and persists in tissues. Thus, RCNV is expected to be a superior vaccine vector for infectious diseases and cancer due to its excellent safety profile, reported vaccine efficacy and ability to replicate in mammalian cells.

Bat guano virome: predominance of dietary viruses from insects and plants plus novel mammalian viruses.

Bats are hosts to a variety of viruses capable of zoonotic transmissions. Because of increased contact between bats, humans, and other animal species, the possibility exists for further cross-species transmissions and ensuing disease outbreaks. We describe here full and partial viral genomes identified using metagenomics in the guano of bats from California and Texas. A total of 34% and 58% of 390,000 sequence reads from bat guano in California and Texas, respectively, were related to eukaryotic viruses, and the largest proportion of those infect insects, reflecting the diet of these insectivorous bats, including members of the viral families Dicistroviridae, Iflaviridae, Tetraviridae, and Nodaviridae and the subfamily Densovirinae. The second largest proportion of virus-related sequences infects plants and fungi, likely reflecting the diet of ingested insects, including members of the viral families Luteoviridae, Secoviridae, Tymoviridae, and Partitiviridae and the genus Sobemovirus. Bat guano viruses related to those infecting mammals comprised the third largest group, including members of the viral families Parvoviridae, Circoviridae, Picornaviridae, Adenoviridae, Poxviridae, Astroviridae, and Coronaviridae. No close relative of known human viral pathogens was identified in these bat populations. Phylogenetic analysis was used to clarify the relationship to known viral taxa of novel sequences detected in bat guano samples, showing that some guano viral sequences fall outside existing taxonomic groups. This initial characterization of the bat guano virome, the first metagenomic analysis of viruses in wild mammals using second-generation sequencing, therefore showed the presence of previously unidentified viral species, genera, and possibly families. Viral metagenomics is a useful tool for genetically characterizing viruses present in animals with the known capability of direct or indirect viral zoonosis to humans.

Viral nucleic acids in live-attenuated vaccines: detection of minority variants and an adventitious virus.

Metagenomics and a panmicrobial microarray were used to examine eight live-attenuated viral vaccines. Viral nucleic acids in trivalent oral poliovirus (OPV), rubella, measles, yellow fever, varicella-zoster, multivalent measles/mumps/rubella, and two rotavirus live vaccines were partially purified, randomly amplified, and pyrosequenced. Over half a million sequence reads were generated covering from 20 to 99% of the attenuated viral genomes at depths reaching up to 8,000 reads per nucleotides. Mutations and minority variants, relative to vaccine strains, not known to affect attenuation were detected in OPV, mumps virus, and varicella-zoster virus. The anticipated detection of endogenous retroviral sequences from the producer avian and primate cells was confirmed. Avian leukosis virus (ALV), previously shown to be noninfectious for humans, was present as RNA in viral particles, while simian retrovirus (SRV) was present as genetically defective DNA. Rotarix, an orally administered rotavirus vaccine, contained porcine circovirus-1 (PCV1), a highly prevalent nonpathogenic pig virus, which has not been shown to be infectious in humans. Hybridization of vaccine nucleic acids to a panmicrobial microarray confirmed the presence of endogenous retroviral and PCV1 nucleic acids. Deep sequencing and microarrays can therefore detect attenuated virus sequence changes, minority variants, and adventitious viruses and help maintain the current safety record of live-attenuated viral vaccines.

Novel circular DNA viruses in stool samples of wild-living chimpanzees.

Viral particles in stool samples from wild-living chimpanzees were analysed using random PCR amplification and sequencing. Sequences encoding proteins distantly related to the replicase protein of single-stranded circular DNA viruses were identified. Inverse PCR was used to amplify and sequence multiple small circular DNA viral genomes. The viral genomes were related in size and genome organization to vertebrate circoviruses and plant geminiviruses but with a different location for the stem-loop structure involved in rolling circle DNA replication. The replicase genes of these viruses were most closely related to those of the much smaller (approximately 1 kb) plant nanovirus circular DNA chromosomes. Because the viruses have characteristics of both animal and plant viruses, we named them chimpanzee stool-associated circular viruses (ChiSCV). Further metagenomic studies of animal samples will greatly increase our knowledge of viral diversity and evolution.

A novel picornavirus associated with gastroenteritis.

A novel picornavirus genome was sequenced, showing 42.6%, 35.2%, and 44.6% of deduced amino acid identities corresponding to the P1, P2, and P3 regions, respectively, of the Aichi virus. Divergent strains of this new virus, which we named salivirus, were detected in 18 stool samples from Nigeria, Tunisia, Nepal, and the United States. A statistical association was seen between virus shedding and unexplained cases of gastroenteritis in Nepal (P = 0.0056). Viruses with approximately 90% nucleotide similarity, named klassevirus, were also recently reported in three cases of unexplained diarrhea from the United States and Australia and in sewage from Spain, reflecting a global distribution and supporting a pathogenic role for this new group of picornaviruses.

Cardioviruses are genetically diverse and cause common enteric infections in South Asian children.

Cardioviruses cause enteric infections in mice and rats which when disseminated have been associated with myocarditis, type 1 diabetes, encephalitis, and multiple sclerosis-like symptoms. Cardioviruses have also been detected at lower frequencies in other mammals. The Cardiovirus genus within the Picornaviridae family is currently made up of two viral species, Theilovirus and Encephalomyocarditis virus. Until recently, only a single strain of cardioviruses (Vilyuisk virus within the Theilovirus species) associated with a geographically restricted and prevalent encephalitis-like condition had been reported to occur in humans. A second theilovirus-related cardiovirus (Saffold virus [SAFV]) was reported in 2007 and subsequently found in respiratory secretions from children with respiratory problems and in stools of both healthy and diarrheic children. Using viral metagenomics, we identified RNA fragments related to SAFV in the stools of Pakistani and Afghani children with nonpolio acute flaccid paralysis (AFP). We sequenced three near-full-length genomes, showing the presence of divergent strains of SAFV and preliminary evidence of a distant recombination event between the ancestors of the Theiler-like viruses of rats and those of human SAFV. Further VP1 sequencing showed the presence of five new SAFV genotypes, doubling the reported genetic diversity of human and animal theiloviruses combined. Both AFP patients and healthy children in Pakistan were found to be excreting SAFV at high frequencies of 9 and 12%, respectively. Further studies are needed to examine the roles of these highly common and diverse SAFV genotypes in nonpolio AFP and other human diseases.

Genomic characterization of novel human parechovirus type.

Using a simple metagenomic approach, we identified a divergent human parechovirus (HPeV) in the stool of a child in Pakistan. Genomic characterization showed this virus was distinct enough from reported HPeV types to qualify as candidate prototype for the seventh HPeV type.

Metagenomic analyses of viruses in stool samples from children with acute flaccid paralysis.

We analyzed viral nucleic acids in stool samples collected from 35 South Asian children with nonpolio acute flaccid paralysis (AFP). Sequence-independent reverse transcription and PCR amplification of capsid-protected, nuclease-resistant viral nucleic acids were followed by DNA sequencing and sequence similarity searches. Limited Sanger sequencing (35 to 240 subclones per sample) identified an average of 1.4 distinct eukaryotic viruses per sample, while pyrosequencing yielded 2.6 viruses per sample. In addition to bacteriophage and plant viruses, we detected known enteric viruses, including rotavirus, adenovirus, picobirnavirus, and human enterovirus species A (HEV-A) to HEV-C, as well as numerous other members of the Picornaviridae family, including parechovirus, Aichi virus, rhinovirus, and human cardiovirus. The viruses with the most divergent sequences relative to those of previously reported viruses included members of a novel Picornaviridae genus and four new viral species (members of the Dicistroviridae, Nodaviridae, and Circoviridae families and the Bocavirus genus). Samples from six healthy contacts of AFP patients were similarly analyzed and also contained numerous viruses, particularly HEV-C, including a potentially novel Enterovirus genotype. Determining the prevalences and pathogenicities of the novel genotypes, species, genera, and potential new viral families identified in this study in different demographic groups will require further studies with different demographic and patient groups, now facilitated by knowledge of these viral genomes.

A highly prevalent and genetically diversified Picornaviridae genus in South Asian children.

Viral metagenomics focused on particle-protected nucleic acids was used on the stools of South Asian children with nonpolio acute flaccid paralysis (AFP). We identified sequences distantly related to Seneca Valley virus and cardioviruses that were then used as genetic footholds to characterize multiple viral species within a previously unreported genus of the Picornaviridae family. The picornaviruses were detected in the stools of >40% of AFP and healthy Pakistani children. A genetically diverse and highly prevalent enteric viral infection, characteristics similar to the Enterovirus genus, was therefore identified substantially expanding the genetic diversity of the RNA viral flora commonly found in children.

Rapid identification of known and new RNA viruses from animal tissues.

Viral surveillance programs or diagnostic labs occasionally obtain infectious samples that fail to be typed by available cell culture, serological, or nucleic acid tests. Five such samples, originating from insect pools, skunk brain, human feces and sewer effluent, collected between 1955 and 1980, resulted in pathology when inoculated into suckling mice. In this study, sequence-independent amplification of partially purified viral nucleic acids and small scale shotgun sequencing was used on mouse brain and muscle tissues. A single viral agent was identified in each sample. For each virus, between 16% to 57% of the viral genome was acquired by sequencing only 42-108 plasmid inserts. Viruses derived from human feces or sewer effluent belonged to the Picornaviridae family and showed between 80% to 91% amino acid identities to known picornaviruses. The complete polyprotein sequence of one virus showed strong similarity to a simian picornavirus sequence in the provisional Sapelovirus genus. Insects and skunk derived viral sequences exhibited amino acid identities ranging from 25% to 98% to the segmented genomes of viruses within the Reoviridae family. Two isolates were highly divergent: one is potentially a new species within the orthoreovirus genus, and the other is a new species within the orbivirus genus. We demonstrate that a simple, inexpensive, and rapid metagenomics approach is effective for identifying known and highly divergent new viruses in homogenized tissues of acutely infected mice.

Disruption of the putative splice acceptor site for SIV(mac239)Vif reveals tight control of SIV splicing and impaired replication in Vif non-permissive cells.

Vif is dispensable for simian immunodeficiency virus (SIV) replication in some cells, termed permissive (i.e., CEM-SS), but not in others, termed non-permissive (i.e., H9, CEMx174, and peripheral blood lymphocytes). Non-permissive cells express the RNA editing enzyme, APOBEC3G. To determine whether vif mRNA could be alternatively spliced, a mutation altering the putative vif splice acceptor site (SA1) was introduced into SIV(mac239) (SIV(Deltavif-SA)). Despite three consensus splice acceptor sites nearby SA1, SIV(Deltavif-SA) did not efficiently generate alternatively spliced vif mRNA. SIV(Deltavif-SA) was growth attenuated in CEMx174 and H9 cells but not in CEM-SS cells. Following SIV(Deltavif-SA), but not SIV(mac239), infection in either H9 or CEMx174 cells viral cDNA contained numerous G to A mutations; no such differences were observed in CEM-SS cells. This pattern is consistent with mutations generated by APOBEC3G in the absence of Vif. Therefore, efficient splicing of SIV vif mRNA is tightly controlled and requires the SA1 site.

L-chicoric acid inhibits human immunodeficiency virus type 1 integration in vivo and is a noncompetitive but reversible inhibitor of HIV-1 integrase in vitro.

The human immunodeficiency virus (HIV) integrase (IN) must covalently join the viral cDNA into a host chromosome for productive HIV infection. l-Chicoric acid (l-CA) enters cells poorly but is a potent inhibitor of IN in vitro. Using quantitative real-time polymerase chain reaction (PCR), l-CA inhibits integration at concentrations from 500 nM to 10 microM but also inhibits entry at concentrations above 1 microM. Using recombinant HIV IN, steady-state kinetic analyses with l-CA were consistent with a noncompetitive or irreversible mechanism of inhibition. IN, in the presence or absence of l-CA, was successively washed. Inhibition of IN diminished, demonstrating that l-CA was reversibly bound to the protein. These data demonstrate that l-CA is a noncompetitive but reversible inhibitor of IN in vitro and of HIV integration in vivo. Thus, l-CA likely interacts with amino acids other than those which bind substrate.

Replication kinetics for divergent type 1 human immunodeficiency viruses using quantitative SYBR green I real-time polymerase chain reaction.

A quantitative and sensitive measure of human immunodeficiency virus type 1 (HIV-1) replication is quantitative real-time polymerase chain reaction (PCR). Real-time PCR using SYBR green I and oligonucleotide primers that amplify early, intermediate, and late products of reverse transcription were optimized to measure HIV-1 replication of clade A, B, C, and D HIV-1 isolates in peripheral blood lymphocytes and in both transformed and viral-transformed CD4+ lymphocyte cell lines. Real-time PCR can detect HIV-1 replication as early as 1 hr postinfection and demonstrates that in established cell lines cDNA can be detected as early as 4 hr postinfection. The first round of HIV-1 replication in established cell lines is complete between 12 and 24 hr postinfection. Furthermore, real-time PCR can detect HIV-1 replication in fewer than 0.1% of cells. Patient isolates replicated at different rates in peripheral blood lymphocytes, with viral cDNA peaking between 48 and 120 hr, depending on the virus being studied. Real-time PCR differentiated the mechanisms of action of drugs targeted at HIV-1 entry, reverse transcription, and proteolytic processing and identified differences in the kinetics of reverse transcription between zidovudine-sensitive and zidovudine-resistant HIV in the presence of zidovudine. In summary, real-time PCR using SYBR green I dye is a sensitive, quantitative, and reproducible measure of replication kinetics for a variety of group M HIV-1 isolates.

Human immunodeficiency virus type-1 integrase containing a glycine to serine mutation at position 140 is attenuated for catalysis and resistant to integrase inhibitors.

L-chicoric acid (L-CA) is a potent inhibitor of HIV integrase (IN) in vitro. In this report, the effects of a glycine to serine mutation at position 140 (G140S) on HIV IN and its effects on IN inhibitor resistance are described. HIV containing the G140S mutation showed a delay in replication. Using real-time polymerase chain reaction, the delay was secondary to a failure in integration. The mutant protein (IN(G140S)) was attenuated approximately four-fold for catalysis under equilibrium conditions compared to wild-type IN (IN(WT)) and attenuated five-fold in steady-state kinetic analysis of disintegration. Fifty percent inhibitory concentration assays were performed with IN inhibitors against both IN proteins in disintegration and strand transfer reactions. IN(G140S) was resistant to both L-CA and L-731,988, a diketoacid. HIV containing the mutation was resistant to both inhibitors as well. The G140S mutation attenuates IN activity and confers resistance to IN inhibitors, suggesting that diketoacids and L-CA interact with a similar binding site on HIV IN.

Dicaffeoyltartaric acid analogues inhibit human immunodeficiency virus type 1 (HIV-1) integrase and HIV-1 replication at nontoxic concentrations.

The human immunodeficiency virus type 1 (HIV-1) is a major health problem worldwide. In this study, 17 analogues of L-chicoric acid, a potent inhibitor of HIV integrase, were studied. Of these analogues, five submicromolar inhibitors of integrase were discovered and 13 compounds with activity against integrase at less than 10 microM were identified. Six demonstrated greater than 10-fold selectivity for HIV replication over cellular toxicity. Ten analogues inhibited HIV replication at nontoxic concentrations. Alteration of the linkages between the two bis-catechol rings, including the use of amides, mixed amide esters, cholate, and alkyl bridges, was explored. Amides were as active as esters but were more toxic in tissue culture. Alkyl and cholate bridges were significantly less potent against HIV-1 integrase in vitro and were inactive against HIV-1 replication. Two amino acid derivates and one digalloylderivative of L-chicoric acid (L-CA) showed improved selectivity over L-CA against integration in cell culture. These data suggest that in addition to the bis-catechols and free carboxylic acid groups reported previously, polar linkages are important constituents for optimal activity against HIV-1 integrase and that new derivatives can be developed with increased specificity for integration over HIV entry in vivo.