Viral Metagenomics Research Articles

Viral Metagenomic Analysis of Japanese Quail (Coturnix japonica) with Enteritis in the Republic of Korea.

We performed viral metagenomics analysis of Japanese quail affected with enteritis to elucidate the viral etiology. Metagenomics generated 21,066,442 sequence reads via high-throughput sequencing, with a mean length of 136 nt. Enrichment in viral sequences suggested that at least three viruses were present in quail samples. Coronavirus and picornavirus were identified and are known as pathogens causing quail enteritis that match the observed morphology. Abundant reads of coronavirus from quail samples yielded four fragment sequences exhibiting six genomes of avian coronavirus. Sequence analysis showed that this quail coronavirus was related to turkey coronavirus and chicken infectious bronchitis virus. Quail picornavirus 8177 bp in size was identified and was similar to the QPV1/HUN/01 virus detected in quails without clinical symptoms in Hungary with 84.6% nucleotide and 94.6% amino acid identity. Our results are useful for understanding the genetic diversity of quail viruses. Further studies must be performed to determine whether quail coronavirus and quail picornavirus are pathogens of the digestive tract of quails.

Mysterious syndrome causing high mortality in wild brown trout in Eastern Switzerland, pathology and search for a possible cause.

Since 2016, annually occurring species-specific die-offs of brown trout (Salmo trutta fario) occurred in the Thur River, situated in the Eastern part of Switzerland. These events lead to drastically reduced population densities in the impacted river regions. Clinical signs in brown trout and mortality were restricted to few weeks in August/September. To characterize the syndrome and to find possible causes, from end of March to November 2018, one-year-old brown trout (Salmo trutta fario) and rainbow trout (Oncorhynchus mykiss) were exposed to water from Thur River, fish were sampled regularly and screened for infectious agents, including viral metagenomics, and pathology was described. Starting approximately four months post-exposure, brown trout showed severe lymphohistiocytic pancarditis and necrotizing and haemorrhagic hepatitis. These lesions were recorded until the end of the experiment in November. Rainbow trout were not affected at any point in time. No infectious agents could be identified so far as cause of disease, especially no viral aetiology. Even if pathogenesis and pathology point in the direction of an infectious agent, a causative relationship could not be confirmed and aetiology remains unclear.

Viral metagenomics in Brazilian multiply transfused patients with sickle cell disease as an indicator for blood transfusion safety.

Patients with sickle cell disease (SCD) are submitted to multiple transfusions in order to increase the oxygen capacity of the blood, decrease blood viscosity, and suppress the sickling of the cells. Multiply transfused patients with SCD represent significant risk of acquiring parenterally transmitted infections. The analysis of the virome profile of high-risk multiply transfused patients with SCD can reveal the presence of parenterally transmitted viruses and therefore be used an indirect approach for evaluation of blood transfusion safety. Blood samples were collected from 45 patients with SCD receiving multiple transfusions and analyzed by metagenomic analyses. The samples were assembled in pools f which were submitted to nucleic acids extraction and sequencing by Illumina NextSeq 550 equipment. For bioinformatic analysis, we used a specific in-house developed pipeline specialized in identification of emerging viruses. The virome composition of SCD patients revealed the presence of commensal viruses represented by anelloviruses and Human Pegivirus-1 (HPgV-1, GB virus C). Contaminant viral sequences belonging to human lentiviruses (rev, env genes), cytomegalovirus and murine leukemia virus were also identified and are attributed to vectors used in the laboratory practice. No novel or unsuspected pathogenic viruses were identified. This study evaluates for the first time the virome of multiply transfused patients with SCD. Exclusively genetic material of commensal viruses was annotated. Therefore, we believe that viral metagenomics applied in patients with high risk for acquiring parenterally transmitted infections can serve as a direct indicator for evaluation of transfusion safety.

Improved diagnosis of viral encephalitis in adult and pediatric hematological patients using viral metagenomics.

Metagenomic sequencing is a powerful technique that enables detection of the full spectrum of pathogens present in any specimen in a single test. Hence, metagenomics is increasingly being applied for detection of viruses in clinical cases with suspected infections of unknown etiology and a large number of relevant potential causes. This is typically the case in patients presenting with encephalitis, in particular when immunity is impaired by underlying disorders. In this study, viral metagenomics has been applied to a cohort of hematological patients with encephalitis of unknown origin. Because viral loads in cerebrospinal fluid of patients with encephalitis are generally low, the technical performance of a metagenomic sequencing protocol with viral enrichment by capture probes targeting all known vertebrate viral sequences was studied. Subsequently, the optimized viral metagenomics protocol was applied to a cohort of hematological patients with encephalitis of unknown origin. Viral enrichment by capture probes increased the viral sequence read count of metagenomics on cerebrospinal fluid samples 100 - 10.000 fold, compared to unenriched metagenomic sequencing. In five out of 41 (12%) hematological patients with encephalitis, a virus was detected by viral metagenomics which had not been detected by current routine diagnostics. BK polyomavirus, hepatitis E virus, human herpes virus-6 and Epstein Barr virus were identified by this unbiased metagenomic approach. This study demonstrated that hematological patients with encephalitis of unknown origin may benefit from early viral metagenomics testing as a single step approach.

Viral metagenomics analysis of kidney donors and recipients: Torque teno virus genotyping and prevalence.

The viral load of the ubiquitous and nonpathogenic torque teno virus (TTV) is associated with the grade of immunosuppression of its host. The development of next-generation sequencing (NGS) allowed to describe the human virome of the blood compartment in transplanted patients, and showed that TTV is the most important part of the virome. This study is a descriptive retrospective pilot study of sequencing plasma samples from 15 matched donors and recipients. After sample processing, nucleic acids were amplified by rolling circle amplification and submitted to NGS by ion proton sequencing technology. Results were analyzed after mapping of reads on the 29 TTV reference genomes and de novo assembling of the reads with MIRA software. The number of TTV species present in donors and recipients was, on average, 12 in donors and 33 in recipients. TTV species predominantly present in donors were TTV-13 and TTV-18; and in recipients were TTV-P15-2, TTV-27, TTV-HD14a, and TTV-22. We highlighted a significant variability in abundance and composition in sequential samples from recipients. Temporal evolution of TTV populations was clearly observed in recipients, but no preferential transmission of a species from donor to recipient was evidenced. Diversity and population expansion were observed in kidney recipients. Further study of TTV species could help assess the potential impact of each species of this virus.

Viral Metagenomic Content Reflects Seawater Ecological Quality in the Coastal Zone

Viruses interfere with their host’s metabolism through the expression of auxiliary metabolic genes (AMGs) that, until now, are mostly studied under large physicochemical gradients. Here, we focus on coastal marine ecosystems and we sequence the viral metagenome (virome) of samples with discrete levels of human-driven disturbances. We aim to describe the relevance of viromics with respect to ecological quality status, defined by the classic seawater trophic index (TRIX). Neither viral (family level) nor bacterial (family level, based on 16S rRNA sequencing) community structure correlated with TRIX. AMGs involved in the Calvin and tricarboxylic acid cycles were found at stations with poor ecological quality, supporting viral lysis by modifying the host’s energy supply. AMGs involved in “non-traditional” energy-production pathways (3HP, sulfur oxidation) were found irrespective of ecological quality, highlighting the importance of recognizing the prevalent metabolic paths and their intermediate byproducts. Various AMGs explained the variability between stations with poor vs. good ecological quality. Our study confirms the pivotal role of the virome content in ecosystem functioning, acting as a “pool” of available functions that may be transferred to the hosts. Further, it suggests that AMGs could be used as an ultra-sensitive metric of energy-production pathways with relevance in the vulnerable coastal zone and its ecological quality.

Viral metagenomics in blood donations with post-donation illness reports from Brazil.

Post-donation illness can be described as appearance of clinical symptoms in blood donors after donation. The consequent call back of the donor to report these symptoms to the blood collection institution is considered a post-donation illness report (PDIR). The most suitable way to examine whether PDIR is related to infection is to apply next-generation sequencing (NGS) and viral metagenomics. Investigation into a PDIR can reveal its importance for transfusion safety and help elaborate strategies for donor education in order to prevent the transfusion transmission of infections which are not routinely tested by the blood collection services. We applied NGS and viral metagenomics on blood donations which were deferred due to a PDIR. Thirty-three PDIR donations obtained in the Blood Center of Ribeirão Preto, Southeast Brazil, were evaluated. Sequencing was performed using Illumina NextSeq 550 (Illumina Inc, San Diego, CA, USA) equipment and the reads obtained for each sample were analysed by specific bioinformatic pipeline for the classification and discovery of emerging viruses. The identified viral agents by metagenomics were directly confirmed by molecular methods. In all PDIR donations, we found abundant reads of commensal viruses belonging to the Anelloviridae family as well as human pegivirus-1. However, we were also able to identify blood donations positive for clinically important viruses like dengue serotype-2 (DENV-2) of the Asian-American genotype and parvovirus B19 (B19V). Both viruses were also confirmed by real-time polymerase chain reaction, detecting DENV-2 RNA in a significant number of cases (7 samples, 21.2%), compared to B19V which was confirmed in 1 case (3.0%). Our study applies for the first time viral metagenomics to evaluate the significance of PDIRs. We confirm the crucial importance of the donor providing a timely PDIR for the prevention of transfusion transmission of viral infections which are not routinely tested in the blood banks worldwide.

Unusual mono-reassortant of a Wa-like G1P[8] species A rotavirus containing a DS-1-like (genotype 2) NSP4 gene.

Species A rotaviruses are a major cause of acute gastroenteritis in infants and young children worldwide. Reassortment is a common phenomenon due to the segmented nature of the rotavirus genome. The complete coding sequences of a species A rotavirus strain isolated from the feces of a child with acute gastroenteritis in Japan in 2018 were determined using an unbiased viral metagenomics approach. The genetic analysis revealed that the rotavirus strain had an unusual genomic constellation (G1-P[8]-I1-R1-C1-M1-A1-N1-T1-E2-H1), suggesting reassortment of a genotype 1 with a genotype 2 rotavirus, from which the NSP4-encoding gene was acquired.

Depicting SARS-CoV-2 faecal viral activity in association with gut microbiota composition in patients with COVID-19

ObjectiveAlthough severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA was detected in faeces of patients with COVID-19, the activity and infectivity of the virus in the GI tract during disease...

PuMA: A papillomavirus genome annotation tool.

High-throughput sequencing technologies provide unprecedented power to identify novel viruses from a wide variety of (environmental) samples. The field of ‘viral metagenomics’ has dramatically expanded our understanding of viral diversity. Viral metagenomic approaches imply that many novel viruses will not be described by researchers who are experts on (the genomic organization of) that virus family. We have developed the papillomavirus annotation tool (PuMA) to provide researchers with a convenient and reproducible method to annotate and report novel papillomaviruses. PuMA currently correctly annotates 99% of the papillomavirus genes when benchmarked against the 655 reference genomes in the papillomavirus episteme. Compared to another viral annotation pipeline, PuMA annotates more viral features while being more accurate. To demonstrate its general applicability, we also developed a preliminary version of PuMA that can annotate polyomaviruses. PuMA is available on GitHub (https://github.com/KVD-lab/puma) and through the iMicrobe online environment (https://www.imicrobe.us/#/apps/puma).

Viral metagenomics in Brazilian Pekin ducks identifies two gyrovirus, including a new species, and the potentially pathogenic duck circovirus

Viral metagenomics coupled to high-throughput sequencing has provided a powerful tool for large-scale detection of known and unknown viruses associated to distinct hosts and environments. Using this approach, known and novel viruses have been characterized from sylvatic and commercial avian hosts, increasing our understanding of the viral diversity in these species. In the present work we applied an exploratory viral metagenomics on organs (spleen, liver and bursa of Fabricious) of Pekin ducks from Southern Brazil. The virome contained sequences related to a known duck pathogen (duck circovirus) and a number of other circular ssDNA viruses. Additionally, we detected avian gyrovirus 9 (to date detected only in human feces) and one new avian gyrovirus species, to which is proposed the name avian gyrovirus 13 (GyV13). This study is expected to contribute to the knowledge of the viral diversity in Pekin ducks.

Evaluation of a commercial exogenous internal process control for diagnostic RNA virus metagenomics from different animal clinical samples

Metagenomic next generation sequencing (mNGS) is increasingly recognized as an important complementary tool to targeted human and animal infectious disease diagnostics. It is, however, sensitive to biases and errors that are currently not systematically evaluated by the implementation of quality controls (QC) for the diagnostic use of mNGS. We evaluated a commercial reagent (Mengovirus extraction control kit, CeraamTools, bioMérieux) as an exogenous internal control for mNGS. It validates the integrity of reagents and workflow, the efficient isolation of viral nucleic acids and the absence of inhibitors in individual samples (verified using a specific qRT-PCR). Moreover, it validates the efficient generation of viral sequence data in individual samples (verified by normalized mengoviral read counts in the metagenomic analysis). We show that when using a completely random metagenomics workflow: (1) Mengovirus RNA can be reproducibly detected in different animal sample types (swine feces and sera, wild bird cloacal swabs), except for tissue samples (swine lung); (2) the Mengovirus control kit does not contain other contaminating viruses that may affect metagenomic experiments (using a cutoff of minimum 1 Kraken classified read per million (RPM)); (3) the addition of 2.17 × 106Mengovirus copies/mL of sample does not affect the virome composition of pig fecal samples or wild bird cloacal swab samples; (4) Mengovirus Cq values (using as cutoff the upper limit of the 99 % confidence interval of Cq values for a given sample matrix) allow the identification of samples with poor viral RNA extraction or high inhibitor load; (5) Mengovirus normalized read counts (cutoff RPM > 1) allow the identification of samples where the viral sequences are outcompeted by host or bacterial target sequences in the random metagenomic workflow. The implementation of two QC testing points, a first one after RNA extraction (Mengoviral qRT-PCR) and a second one after metagenomic data analysis provide valuable information for the validation of individual samples and results. Their implementation in addition to external controls validating runs or experiments should be carefully considered for a given sample type and workflow.

Freshwater viral metagenome reveals novel and functional phage-borne antibiotic resistance genes

BackgroundAntibiotic resistance developed by bacteria is a significant threat to global health. Antibiotic resistance genes (ARGs) spread across different bacterial populations through multiple dissemination routes, including horizontal gene transfer mediated by bacteriophages. ARGs carried by bacteriophages are considered especially threatening due to their prolonged persistence in the environment, fast replication rates, and ability to infect diverse bacterial hosts. Several studies employing qPCR and viral metagenomics have shown that viral fraction and viral sequence reads in clinical and environmental samples carry many ARGs. However, only a few ARGs have been found in viral contigs assembled from metagenome reads, with most of these genes lacking effective antibiotic resistance phenotypes. Owing to the wide application of viral metagenomics, nevertheless, different classes of ARGs are being continuously found in viral metagenomes acquired from diverse environments. As such, the presence and functionality of ARGs encoded by bacteriophages remain up for debate.ResultsWe evaluated ARGs excavated from viral contigs recovered from urban surface water viral metagenome data. In virome reads and contigs, diverse ARGs, including polymyxin resistance genes, multidrug efflux proteins, and β-lactamases, were identified. In particular, when a lenient threshold of e value of ≤ 1 × e−5 and query coverage of ≥ 60% were employed in the Resfams database, the novel β-lactamases blaHRV-1 and blaHRVM-1 were found. These genes had unique sequences, forming distinct clades of class A and subclass B3 β-lactamases, respectively. Minimum inhibitory concentration analyses for E. coli strains harboring blaHRV-1 and blaHRVM-1 and catalytic kinetics of purified HRV-1 and HRVM-1 showed reduced susceptibility to penicillin, narrow- and extended-spectrum cephalosporins, and carbapenems. These genes were also found in bacterial metagenomes, indicating that they were harbored by actively infecting phages.ConclusionOur results showed that viruses in the environment carry as-yet-unreported functional ARGs, albeit in small quantities. We thereby suggest that environmental bacteriophages could be reservoirs of widely variable, unknown ARGs that could be disseminated via virus-host interactions.4ZRryVp6ySP7-sdMyhXuoSVideo abstract.

The Role of Rhizosphere Bacteriophages in Plant Health.

Microbiomes and their hosts influence each other; for instance, the microbiome improves host fitness, whereas the host supports microbiome nutrition. Most studies on this topic have focused on the role of bacteria and fungi, although research on viruses that infect bacteria, known as 'bacteriophages' (phages), has gained importance due to the potential role bacteriophages play in the resilience and functionality of the gut microbiome. Like the gut microbiome, the rhizosphere harbors a complex microbiome, but little is known about the role of phages in this ecosystem. In this opinion, we extend the knowledge gained in human gut virus metagenomics (viromics) to disentangle the potential role of phages in driving the resilience and functionality of the rhizosphere microbiome. We propose that future comparative studies across environments are necessary to unravel the underlying mechanisms through which phages drive the composition and functionality of the rhizosphere microbiome and its interaction with the plant host. Importantly, such understanding might generate strategies to improve plant resistance and resilience in the context of climate change.

The Need of Non-traditional Techniques to Screen for the Virus

BACKGROUND: At the present moment, the etiological diagnosis of SARS-CoV-2 is based on the polymerase chain reaction (PCR). False negative cases are increasingly reported in several studies using reverse transcription-PCR (RT-PCR). For example, the positive rate of RT-PCR for throat swabs was reported to be about 60% in early stage of COVID-19.
 AIM: We aimed to present metagenomic next-generation sequencing (mNGS) as a potential tool to detect pathogens.
 METHODS: In the recent year, mNGS is shown the potential to detect pathogens without the need of hypothesis guided approach and is proven to be highly effective.
 RESULTS: A recent prospective study in the United States compared the diagnostic performance of routine diagnostic tests with mNGS and showed that mNGS detected a bacteria or virus in the CSF of 13 of 58 patients presenting with meningoencephalitis who were negative for or not assessed with routine diagnostic test including PCR. NGS also has the advantage to cover entire viral genomes.
 CONCLUSION: As viral metagenomics has significantly improved in recent years and become more cost effective, we think that a change in the approach toward a shot-gun metagenomic testing should be explored and could potentially aid the diagnosis of COVID-19 cases and the management of this pandemic.

Virome of a Feline Outbreak of Diarrhea and Vomiting Includes Bocaviruses and a Novel Chapparvovirus.

An unexplained outbreak of feline diarrhea and vomiting, negative for common enteric viral and bacterial pathogens, was subjected to viral metagenomics and PCR. We characterized from fecal samples the genome of a novel chapparvovirus we named fechavirus that was shed by 8/17 affected cats and identified three different feline bocaviruses shed by 9/17 cats. Also detected were nucleic acids from attenuated vaccine viruses, members of the normal feline virome, viruses found in only one or two cases, and viruses likely derived from ingested food products. Epidemiological investigation of disease signs, time of onset, and transfers of affected cats between three facilities support a possible role for this new chapparvovirus in a highly contagious feline diarrhea and vomiting disease.

Genomic Analyses of Potential Novel Recombinant Human Adenovirus C in Brazil.

Human Adenovirus species C (HAdV-C) is the most common etiologic agent of respiratory disease. In the present study, we characterized the nearly full-length genome of one potential new HAdV-C recombinant strain constituted by Penton and Fiber proteins belonging to type 89 and a chimeric Hexon protein of types 1 and 89. By using viral metagenomics techniques, we screened out, in the states of Tocantins and Pará, Northern and North regions of Brazil, from 2010 to 2016, 251 fecal samples of children between 0.5 to 2.5 years old. These children were presenting acute diarrhea not associated with common pathogens (i.e., rotavirus, norovirus). We identified two HAdV-C strains in two distinct patients. Phylogenetic analysis performed using all complete genomes available at GenBank database indicated that one strain (HAdV-C BR-245) belonged to type 1. The phylogenetic analysis also indicated that the second strain (HAdV-C BR-211) was located at the base of the clade formed by the newly HAdV-C strains type 89. Recombination analysis revealed that strain HAdV-C BR-211 is a chimera in which the variable regions of Hexon gene combined HAdV-C1 and HAdV-C89 sequences. Therefore, HAdV-C BR-211 strain possesses a genomic backbone of type HAdV-C89 and a unique insertion of HAdV-C1 in the Hexon sequence. Recombination may play an important driving force in HAdV-C diversity and evolution. Studies employing complete genomic sequencing on circulating HAdV-C strains in Brazil are needed to understand the clinical significance of the presented data.

Viral metagenomics performed in patients with acute febrile syndrome during Toxoplasma gondii outbreak in south Brazil

Toxoplasmosis is a zoonotic infection caused by the protozoan parasite Toxoplasma gondii. The infection is widely disseminated in the human population and is usually benign or asymptomatic. Systemic T. gondii infection presents risks for pregnant women and AIDS patients. Although rare, T. gondii can cause outbreaks in urban centers. The origin of these outbreaks is not completely understood but probably results from introduction of zoonotic T. gondii strains in the population. During such outbreaks other pathogens which mimic T. gondii acute febrile syndrome may also circulate; therefore, detailed investigation of the outbreak is of extreme importance. In this study we performed viral metagenomics next-generation sequencing (mNGS) in patient samples obtained during T. gondii outbreak in Santa Maria city, South Brazil. Specific bioinformatics pipelines specialized in virus discovery were applied in order to identify co-circulating vial agents. Epstein Barr virus and Parvovirus B19 contigs were assembled and these viruses can cause symptoms similar to toxoplasmosis. In conclusion, our findings show the importance of Metagenomics next generation sequencing (mNGS) use to help characterize the outbreak more completely and in the management of the affected patients.

Viral metagenomics revealed diverse CRESS-DNA virus genomes in faeces of forest musk deer

BackgroundMusk deer can produce musk which has high medicinal value and is closely related to human health. Viruses in forest musk deer both threaten the health of forest musk deer and human beings.MethodsUsing viral metagenomics we investigated the virome in 85 faeces samples collected from forest musk deer.ResultsIn this article, eight novel CRESS-DNA viruses were characterized, whole genomes were 2148 nt–3852 nt in length. Phylogenetic analysis indicated that some viral genomes were part of four different groups of CRESS-DNA virus belonging in the unclassified CRESS-DNA virus, Smacoviridae, pCPa-like virus and pPAPh2-like virus. UJSL001 (MN621482), UJSL003 (MN621469) and UJSL017 (MN621476) fall into the branch of unclassified CRESS-DNA virus (CRESSV1–2), UJSL002 (MN621468), UJSL004 (MN621481) and UJSL007 (MN621470) belong to the cluster of Smacoviridae, UJSL005 (MN604398) showing close relationship with pCPa-like (pCRESS4–8) clusters and UJSL006 (MN621480) clustered into the branch of pPAPh2-like (pCRESS9) virus, respectively.ConclusionThe virome in faeces samples of forest musk deer from Chengdu, Sichuan province, China was revealed, which further characterized the diversity of viruses in forest musk deer intestinal tract.

Genomic characterization of Changuinola viruses from Panama: evidence for multiple genome segment reassortment.

The complete coding sequences of five divergent strains of Changuinola virus (CGLV), collected over a 16-year period in Panama, were determined, using viral metagenomics. Each strain had 10 RNA segments that encoded structural and non-structural proteins with amino acid identities ranging from 33 to 99% with sequences of other 15 members of the Changuinola virus (Reoviridae: Orbivirus) species group. Genetic analyses of the five Panamanian virus strains revealed probable reassortment among multiple segments of the viruses.