Long-Read Nanopore-Based Sequencing of Anelloviruses.

Routinely used metagenomic next-generation sequencing (mNGS) techniques often fail to detect low-level viremia (<104 copies/mL) and appear biased towards viruses with linear genomes. These limitations hinder the capacity to comprehensively characterize viral infections, such as those attributed to the Anelloviridae family. These near ubiquitous non-pathogenic components of the human virome have circular single-stranded DNA genomes that vary in size from 2.0 to 3.9 kb and exhibit high genetic diversity. Hence, species identification using short reads can be challenging. Here, we introduce a rolling circle amplification (RCA)-based metagenomic sequencing protocol tailored for circular single-stranded DNA genomes, utilizing the long-read Oxford Nanopore platform. The approach was assessed by sequencing anelloviruses in plasma drawn from people who inject drugs (PWID) in two geographically distinct cohorts. We detail the methodological adjustments implemented to overcome difficulties inherent in sequencing circular genomes and describe a computational pipeline focused on anellovirus detection. We assessed our protocol across various sample dilutions and successfully differentiated anellovirus sequences in conditions simulating mixed infections. This method provides a robust framework for the comprehensive characterization of circular viruses within the human virome using the Oxford Nanopore.

Signatures of COVID-19 Severity and Immune Response in the Respiratory Tract Microbiome.

Viral infection of the respiratory tract can be associated with propagating effects on the airway microbiome, and microbiome dysbiosis may influence viral disease. Here, we investigated the respiratory tract microbiome in coronavirus disease 2019 (COVID-19) and its relationship to disease severity, systemic immunologic features, and outcomes. We examined 507 oropharyngeal, nasopharyngeal, and endotracheal samples from 83 hospitalized COVID-19 patients as well as non-COVID patients and healthy controls. Bacterial communities were interrogated using 16S rRNA gene sequencing, and the commensal DNA viruses Anelloviridae and Redondoviridae were quantified by qPCR. We found that COVID-19 patients had upper respiratory microbiome dysbiosis and greater change over time than critically ill patients without COVID-19. Oropharyngeal microbiome diversity at the first time point correlated inversely with disease severity during hospitalization. Microbiome composition was also associated with systemic immune parameters in blood, as measured by lymphocyte/neutrophil ratios and immune profiling of peripheral blood mononuclear cells. Intubated patients showed patient-specific lung microbiome communities that were frequently highly dynamic, with prominence of Staphylococcus. Anelloviridae and Redondoviridae showed more frequent colonization and higher titers in severe disease. Machine learning analysis demonstrated that integrated features of the microbiome at early sampling points had high power to discriminate ultimate level of COVID-19 severity. Thus, the respiratory tract microbiome and commensal viruses are disturbed in COVID-19 and correlate with systemic immune parameters, and early microbiome features discriminate disease severity. Future studies should address clinical consequences of airway dysbiosis in COVID-19, its possible use as biomarkers, and the role of bacterial and viral taxa identified here in COVID-19 pathogenesis. IMPORTANCE COVID-19, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection of the respiratory tract, results in highly variable outcomes ranging from minimal illness to death, but the reasons for this are not well understood. We investigated the respiratory tract bacterial microbiome and small commensal DNA viruses in hospitalized COVID-19 patients and found that each was markedly abnormal compared to that in healthy people and differed from that in critically ill patients without COVID-19. Early airway samples tracked with the level of COVID-19 illness reached during hospitalization, and the airway microbiome also correlated with immune parameters in blood. These findings raise questions about the mechanisms linking SARS-CoV-2 infection and other microbial inhabitants of the airway, including whether the microbiome might regulate severity of COVID-19 disease and/or whether early microbiome features might serve as biomarkers to discriminate disease severity.

A novel anellovirus from hospitalized neonates.

Here, using viral metagenomics combined with conventional PCR, the complete genome sequence of a novel anellovirus (named anel-ch-zj and GenBank no. MT157223) from nasopharynx secretion specimens from hospitalized neonates was determined, and the deduced amino acid sequence of its ODF1 protein was found to be only 33.19%-39.33% identical to those of related anelloviruses with sequences available in the GenBank database, suggesting that it represents a putative new genus within the family Anelloviridae. PCR screening of 135 samples (including 45 nasopharynx secretion, 45 blood, and 45 fecal specimens collected from 45 individual hospitalized neonates) indicated that two nasopharynx secretion, two blood, and four fecal samples were positive for anel-ch-zj. Further PCR screening of 50 blood samples, 115 fecal samples, and 396 nasopharynx secretions collected from hospitalized children 1-5 years old did not yield any positive results. Whether this novel anellovirus detected in neonates is associated with specific diseases needs future investigation.

Identification of novel circovirus and anelloviruses from wolverines using a non-invasive faecal sampling approach.

Viruses in the families Circoviridae and Anelloviridae have circular single-stranded DNA genomes and have been identified in various animal species. Some members of the Circoviridae family such as beak and feather disease and porcine circovirus have been found to cause disease in their host animals. Anelloviruses on the other hand have not been identified to cause disease in their hosts but are highly prevalent in mammalian species. Using a non-invasive sampling approach, we identified novel circovirus and anelloviruses from faecal samples of wolverines dwelling in Montana, USA. Wolverines are forest carnivores that feed on a wide variety of carrion and other prey species, and they occupy diverse habitats across northern Europe to North America. Little is known about viruses associated with wild wolverines. Our investigation of the faecal samples resulted in the identification of a novel circovirus from three out of four wolverine samples, two collected in 2018 and one in 2019. Comparison with other circoviruses shows it is most closely related to a porcine circovirus 3, sharing ~69% identity. Additionally, three anellovirus genomes were recovered from two wolverine faecal samples which share 68--69% ORF1 nucleotide similarity with an anellovirus from another mustelid species, pine martens. Here we identify novel single-stranded DNA viruses associated with wolverine and open up new avenues for research.

Deep viral blood metagenomics reveals extensive anellovirus diversity in healthy humans.

Human blood metagenomics has revealed the presence of different types of viruses in apparently healthy subjects. By far, anelloviruses constitute the viral family that is more frequently found in human blood, although amplification biases and contaminations pose a major challenge in this field. To investigate this further, we subjected pooled plasma samples from 120 healthy donors in Spain to high-speed centrifugation, RNA and DNA extraction, random amplification, and massive parallel sequencing. Our results confirm the extensive presence of anelloviruses in such samples, which represented nearly 97% of the total viral sequence reads obtained. We assembled 114 different viral genomes belonging to this family, revealing remarkable diversity. Phylogenetic analysis of ORF1 suggested 28 potentially novel anellovirus species, 24 of which were validated by Sanger sequencing to discard artifacts. These findings underscore the importance of implementing more efficient purification procedures that enrich the viral fraction as an essential step in virome studies and question the suggested pathological role of anelloviruses.

Diverse cressdnaviruses and an anellovirus identified in the fecal samples of yellow-bellied marmots.

Over that last decade, coupling multiple strand displacement approaches with high throughput sequencing have resulted in the identification of genomes of diverse groups of small circular DNA viruses. Using a similar approach but with recovery of complete genomes by PCR, we identified a diverse group of single-stranded viruses in yellow-bellied marmot (Marmota flaviventer) fecal samples. From 13 fecal samples we identified viruses in the family Genomoviridae (n = 7) and Anelloviridae (n = 1), and several others that ware part of the larger Cressdnaviricota phylum but not within established families (n = 19). There were also circular DNA molecules identified (n = 4) that appear to encode one viral-like gene and have genomes of <1545 nts. This study gives a snapshot of viruses associated with marmots based on fecal sampling.

Detection of severe fever with thrombocytopenia syndrome virus and other viruses in cats via unbiased next-generation sequencing.

We used unbiased next-generation sequencing (NGS) to detect unknown viruses in cats. Serum or plasma samples were obtained from clinically ill cats with suspected acute viral infections. Nucleic acid was extracted from serum or plasma samples to construct a complementary DNA library for NGS. Comprehensive nucleotide sequencing analyses enabled detection of the genomes of various viruses, including the severe fever with thrombocytopenia syndrome virus, feline immunodeficiency virus, feline morbillivirus, parvovirus, and Torque teno felis virus. Our findings indicate that comprehensive nucleotide analyses of serum or plasma samples can be used to detect infections with unknown viruses in cats.

Diversity and dynamic changes of anelloviruses in plasma following allogeneic hematopoietic stem cell transplantation.

Monitoring of alphatorquevirus (torque teno virus [TTV]) DNA in plasma may prove to be useful to assess the net state of immune competence following allogeneic hematopoietic stem cell transplantation (allo-HSCT). There are scarce data published on the prevalence of beta (torque teno mini virus [TTMV]) and gammatorqueviruses (torque teno midi virus [TTMDV]) and, in particular, on the dynamics of anelloviruses in allo-HSCT patients. Twenty-five allo-HSCT recipients with available plasma specimens obtained before conditioning and after engraftment were included. Degenerated primers targeting a highly conserved genomic sequence across all anelloviruses were designed for genomic amplification and high-throughput sequencing. Co-detection of TTV, TTMV, and TTMDV both in pre-transplant and post-engraftment plasma specimens was documented in more than two-thirds of patients. The use of quantitative real-time polymerase chain reaction (PCR) assays targeting TTMV and TTMDV in addition to TTV may add value to TTV-specific PCR assays in the inference of the net state of immunosuppresion or immune competence in this clinical setting.

Complex Virome in a Mesenteric Lymph Node from a Californian Sea Lion (Zalophus Californianus) with Polyserositis and Steatitis.

An emaciated subadult free-ranging California sea lion (Csl or Zalophus californianus) died following stranding with lesions similar to 11 other stranded animals characterized by chronic disseminated granulomatous inflammation with necrotizing steatitis and vasculitis, involving visceral adipose tissues in the thoracic and peritoneal cavities. Histologically, affected tissues had extensive accumulations of macrophages with perivascular lymphocytes, plasma cells, and fewer neutrophils. Using viral metagenomics on a mesenteric lymph node six mammalian viruses were identified consisting of novel parvovirus, polyomavirus, rotavirus, anellovirus, and previously described Csl adenovirus 1 and Csl bocavirus 4. The causal or contributory role of these viruses to the gross and histologic lesions of this sea lion remains to be determined.

Dynamics of the Stool Virome in Very Early-Onset Inflammatory Bowel Disease.

BACKGROUND AND AIMS: Dysbiosis of the gut microbiota is a well-known correlate of the pathogenesis of inflammatory bowel disease [IBD]. However, few studies have examined the microbiome in very early-onset [VEO] IBD, which is defined as onset of IBD before 6 years of age. Here we focus on the viral portion of the microbiome-the virome-to assess possible viral associations with disease processes, reasoning that any viruses potentially associated with IBD might grow more robustly in younger subjects, and so be more detectable. METHODS: Virus-like particles [VLPs] were purified from stool samples collected from patients with VEO-IBD [n = 54] and healthy controls [n = 23], and characterized by DNA and RNA sequencing and VLP particle counts. RESULTS: The total number of VLPs was not significantly different between VEO-IBD and healthy controls. For bacterial viruses, the VEO-IBD subjects were found to have a higher ratio of Caudovirales vs to Microviridae compared to healthy controls. An increase in Caudovirales was also associated with immunosuppressive therapy. For viruses infecting human cells, Anelloviridae showed higher prevalence in VEO-IBD compared to healthy controls. Within the VEO-IBD group, higher levels of Anelloviridae DNA were also positively associated with immunosuppressive treatment. To search for new viruses, short sequences enriched in VEO-IBD samples were identified, and some could be validated in an independent cohort, although none was clearly viral; this provides sequence tags to interrogate in future studies. CONCLUSIONS: These data thus document perturbations to normal viral populations associated with VEO-IBD, and provide a biomarker-Anelloviridae DNA levels-potentially useful for reporting the effectiveness of immunosuppression.

Semen virome of men with HIV on or off antiretroviral treatment.

OBJECTIVES: Improving immune status of people living with HIV through antiretroviral therapy (ART) may also reduce shedding of other viruses in semen. We characterized the seminal fluid virome of men with HIV and tested potential associations between viruses present and CD4 T-cell count, HIV viremia, and antiretroviral therapy (ART) status. DESIGN AND METHODS: Metagenomics was used to enrich and sequence viral nucleic acids from the seminal fluid of 55 semen samples from 42 men living with HIV from San Francisco with a median age of 33 (IQR, 28.7-45) and median CD4 T-cell counts of 837 cells/µl (IQR, 258-1571 cells/µl). All samples were collected between 2005 and 2015, and ART status was ascertained from medical records. RESULTS: Anelloviruses, cytomegalovirus (CMV), and multiple genotypes of human papillomaviruses were detected. Participants shed from 0 to 4 distinct human viruses. Longitudinally collected seminal fluid samples showed changes in the viruses shed. Viruses were more frequently shed by individuals with detectable HIV viremia (43.7 vs. 15.4%, P = 0.042). A trend was seen for increased shedding by individuals who were not on ART (42.8 vs. 17.8%, P = 0.082) or with CD4 T-cell count less than 350 cells/µl (35.3 vs. 20%, P = 0.27). CONCLUSION: Seminal fluid from men with HIV from San Francisco contains nucleic acids from three different DNA viral families. A greater number of viruses, particularly CMV, were shed by participants with detectable HIV viremia (18.9 vs. 0%, P = 0.022). Control of viremia through ART may lower shedding of other viruses in semen in addition to HIV.

In-depth serum virome analysis in patients with acute liver failure with indeterminate etiology.

In clinical virome research, whole-genome/transcriptome amplification is required when starting material is limited. An improved method, named "template-dependent multiple displacement amplification" (tdMDA), has recently been developed in our lab (Wang et al. in BioTechniques 63:21-25. https://doi.org/10.2144/000114566, 2017). In combination with Illumina sequencing and bioinformatics pipelines, its application in virome sequencing was explored using a serum sample from a patient with chronic hepatitis C virus (HCV) infection. In comparison to an amplification-free procedure, virome sequencing via tdMDA showed a 9.47-fold enrichment for HCV-mapped reads and, accordingly, an increase in HCV genome coverage from 28.5% to 70.1%. Eight serum samples from acute patients liver failure (ALF) with or without known etiology were then used for virome sequencing with an average depth at 94,913x. Both similarity-based (mapping, NCBI BLASTn, BLASTp, and profile hidden Markov model analysis) and similarity-independent methods (machine-learning algorithms) identified viruses from multiple families, including Herpesviridae, Picornaviridae, Myoviridae, and Anelloviridae. However, their commensal nature and cross-detection ruled out an etiological interpretation. Together with a lack of detection of novel viruses in a comprehensive analysis at a resolution of single reads, these data indicate that viral agents might be rare in ALF cases with indeterminate etiology.

High-Throughput Sequencing-Based Investigation of Viruses in Human Cancers by Multienrichment Approach.

BACKGROUND: Viruses and other infectious agents cause more than 15% of human cancer cases. High-throughput sequencing-based studies of virus-cancer associations have mainly focused on cancer transcriptome data. METHODS: In this study, we applied a diverse selection of presequencing enrichment methods targeting all major viral groups, to characterize the viruses present in 197 samples from 18 sample types of cancerous origin. Using high-throughput sequencing, we generated 710 datasets constituting 57 billion sequencing reads. RESULTS: Detailed in silico investigation of the viral content, including exclusion of viral artefacts, from de novo assembled contigs and individual sequencing reads yielded a map of the viruses detected. Our data reveal a virome dominated by papillomaviruses, anelloviruses, herpesviruses, and parvoviruses. More than half of the included samples contained 1 or more viruses; however, no link between specific viruses and cancer types were found. CONCLUSIONS: Our study sheds light on viral presence in cancers and provides highly relevant virome data for future reference.

Metagenomic Next-Generation Sequencing Reveals Individual Composition and Dynamics of Anelloviruses during Autologous Stem Cell Transplant Recipient Management.

Over recent years, there has been increasing interest in the use of the anelloviruses, the major component of the human virome, for the prediction of post-transplant complications such as severe infections. Due to an important diversity, the comprehensive characterization of this viral family over time has been poorly studied. To overcome this challenge, we used a metagenomic next-generation sequencing (mNGS) approach with the aim of determining the individual anellovirus profile of autologous stem cell transplant (ASCT) patients. We conducted a prospective pilot study on a homogeneous patient cohort regarding the chemotherapy regimens that included 10 ASCT recipients. A validated viral mNGS workflow was used on 108 plasma samples collected at 11 time points from diagnosis to 90 days post-transplantation. A complex interindividual variability in terms of abundance and composition was noticed. In particular, a strong sex effect was found and confirmed using quantitative PCR targeting torque teno virus, the most abundant anellovirus. Interestingly, an important turnover in the anellovirus composition was observed during the course of the disease revealing a strong intra-individual variability. Although more studies are needed to better understand anellovirus dynamics, these findings are of prime importance for their future use as biomarkers of immune competence.

Torquetenovirus detection in exosomes enriched vesicles circulating in human plasma samples.

BACKGROUND: Torquetenovirus (TTV) belongs to Anelloviridae family, infects nearly all people indefinitely without causing overt disease establishing a fine and successful interaction with the host. Increasing evidence have shown some human viruses exploit extracellular vesicles thereby helping viral persistence in the host. Here, the presence of TTV in extracellular vesicles circulating in human plasma was investigated. METHODS: TTV DNA was quantified in plasma-derived exosomes from 122 samples collected from 97 diseased patients and 25 healthy donors. Exosomes enriched vesicles (EEVs) were extracted from plasma and characterized by Nanoparticle tracking analysis, by western blot for presence of tetraspanin CD63, CD81 and annexin II protein and, finally, by electron microscopy (EM). Presence and quantitation of TTV DNA were assessed with an universal single step real-time TaqMan PCR assay. RESULTS: Preliminary investigation showed that the human plasma extracted extracellular vesicles exhibited a main size of 70 nm, had concentration of 2.5 × 109/ml, and scored positive for tetraspanin CD63, CD81 and annexin II, typical characteristic of the exosomes vesicles. EEVs extracted from pooled plasma with TTV DNA viremia of 9.7 × 104 copies/ml showed to contain 6.3 × 102 TTV copies/ml, corresponding to 0.65% of total viral load. Important, TTV yield changed significantly following freezing/thawing, detergents and DNAse treatment of plasma before EEVs extraction. EEVs purified by sucrose-density gradient centrifugation and analysis of gradient fraction positive for exosomes marker CD63 harbored 102 TTV copies/ml. Moreover, EM evidenced the presence of TTV-like particles in EEVs. Successive investigation of plasma EEVs from 122 subjects (37 HIV-positive, 20 HCV infected, 20 HBV infected, 20 kidney transplant recipients, and 25 healthy) reported TTV DNA detection in 42 (34%) of the viremic samples (37 were from diseased patients and 5 from healthy people) at a mean level of 4.8 × 103 copies/ml. The examination of EEVs selected samples reported the presence of TTV genogroup 1, 3, 4 and 5, with genogroup 3 highly observed. CONCLUSIONS: Collectively, although these observations should be confirmed by further studies, circulation of TTV particles in EEVs opens new avenues and mechanistic insights on the molecular strategies adopted by anelloviruses to persist in the host.

Identification and whole genome characterization of novel anelloviruses in masked palm civets (Paguma larvata): Segregation into four distinct clades.

The members of the family Anelloviridae are small and single-stranded DNA viruses with marked diversity in sequence and length, which ubiquitously infect many vertebrates, including mammals, birds and reptiles. The anelloviruses isolated from mammals are currently classified into 11 assigned and four proposed genera; some anelloviruses remain unassigned. The present study was conducted to identify anelloviruses in wild-caught masked palm civets (Paguma larvata) in Japan using a rolling-circle amplification method. Thirteen novel anellovirus strains were identified from 8 of 10 masked palm civets and their entire genomic sequences (2039-2535 nucleotides) were determined; they were classifiable into four distinct clades. Comparative analyses of all reported anelloviruses for which the entire or near-entire genomic sequences have been determined, including the 13 strains obtained in the present study, revealed that anelloviruses can provisionally be classified into 20 clades, which may correspond to 20 genera (including 11 assigned and four proposed genera) by a >70% amino acid sequence difference in open reading frame 1 (ORF1). This study suggested that novel anelloviruses of marked diversity are circulating in animals worldwide, and that the rolling-circle amplification method would be useful for identifying novel anelloviruses and other viruses with a circular DNA genome.

Novel anelloviruses identified in buccal swabs of Antarctic fur seals.

Viral diversity associated with Antarctic wildlife remains poorly studied. Nonetheless, over the past 5 years, there has been a concerted effort using viral metagenomics approaches to identify and characterize viruses associated with Antarctic pinniped and avian species. Here we used a viral metagenomics approach to identify circular DNA viruses in buccal swab samples from Antarctic fur seals (Arctocephalus gazella) breeding on Livingston Island, Antarctica during the 2016/2017 field season. We identified two new lineages of anelloviruses, torque teno Arctocephalus gazella virus (TTAgV)-1 (2083 nts) and -2 (2127-2143 nts), which are most closely related to and cluster with a previously identified anellovirus associated with California sea lions (Zalophus californianus) sharing ~ 60% genome-wide pairwise identity. The ORF1 of TTAgVs share 26-41% amino acid similarity to the ORF1 of other previously identified pinniped-associated anelloviruses.

Prevalence of human anelloviruses in Romanian healthy subjects and patients with common pathologies.

BACKGROUND: Human anelloviruses (TTV, TTMDV and TTMV) are at high prevalence all across the globe, having also a controversial disease-inducing potential. This study aimed to estimate the prevalence of anelloviral DNA in the Romanian human population and to investigate the association of infections with common pathologies in Romanian population. METHODS: After informed consent, blood samples were collected from 2000 subjects represented by: clinically healthy individuals (n = 701) and a group of patients with pathologies linked to low grade inflammation or alteration of carbohydrate metabolism (n = 1299). All samples were analysed for the presence of TTV, TTMDV and TTMV DNA by hemi-nested PCR. RESULTS: The prevalence of TTV, TTMDV and TTMV in the studied population was 68.2, 54.4%, respectively 40.1%, lower than the recent reports from other geographic regions. The three viral species were significantly more frequent in the group of patients compared to the healthy subjects and were associated with type 2 diabetes mellitus. The presence of anelloviral DNA was also associated with medical procedures (e.g. haemodialysis/transfusions, surgical procedures) and previous hepatitis A virus infection. Lifestyle choices related to alcohol consumption, smoking, physical activity and living environment were not associated with differences in distribution of the three viruses. CONCLUSION: Further evidence is needed to establish a correlation between infection with human anelloviruses and a pathology or group of pathologies.

Viremia preceding multiple sclerosis: Two nested case-control studies.

Infections have been suggested to be involved in Multiple Sclerosis (MS). We used metagenomic sequencing to detect both known and yet unknown microorganisms in 2 nested case control studies of MS. Two different cohorts were followed for MS using registry linkages. Serum samples taken before diagnosis as well as samples from matched control subjects were selected. In cohort1 with 75 cases and 75 controls, most viral reads were Anelloviridae-related and >95% detected among the cases. Among samples taken up to 2 years before MS diagnosis, Anellovirus species TTMV1, TTMV6 and TTV27 were significantly more common among cases. In cohort2, 93 cases and 93 controls were tested under the pre-specified hypothesis that the same association would be found. Although most viral reads were again related to Anelloviridae, no significant case-control differences were seen. We conclude that the Anelloviridae-MS association may be due to multiple hypothesis testing, but other explanations are possible.