Global genome analysis reveals a vast and dynamic anellovirus landscape within the human virome

The Blood Virome: A new frontier in biomedical science

Recent advances in metagenomic testing opened a new window into the mammalian blood virome. Comprised of well-known viruses like human immunodeficiency virus, hepatitis C virus, and hepatitis B virus, the virome also includes many other eukaryotic viruses and phages whose medical significance, lifecycle, epidemiology, and impact on human health are less well known and thus regarded as commensals. This review synthesizes available information for the so-called commensal virome members that circulate in the blood of humans considering their restriction to and interaction with the human host, their natural history, and their impact on human health and physiology.

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.

Hybrid capture shotgun sequencing detected unexpected viruses in the cerebrospinal fluid of children with acute meningitis and encephalitis

PURPOSE

Investigation of undiagnosed cases of infectious neurological diseases, especially in the paediatric population, remains a challenge. This study aimed to enhance understanding of viruses in CSF from children with clinically diagnosed meningitis and/or encephalitis (M/ME) of unknown aetiology using shotgun sequencing enhanced by hybrid capture (HCSS).

METHODS

A single-centre prospective study was conducted at Sant Joan de Déu University Hospital, Barcelona, involving 40 M/ME episodes of unknown aetiology, recruited from May 2021 to July 2022. All participants had previously tested negative with the FilmArray Meningitis/Encephalitis Panel. HCSS was used to detect viral nucleic acid in the patients' CSF. Sequencing was performed on Illumina NovaSeq platform. Raw sequence data were analysed using CZ ID metagenomics and PikaVirus bioinformatics pipelines.

RESULTS

Forty episodes of M/ME of unknown aetiology in 39 children were analysed by HCSS. A significant viral detection in 30 CSF samples was obtained, including six parechovirus A, three enterovirus ACD, four polyomavirus 5, three HHV-7, two BKV, one HSV-1, one VZV, two CMV, one EBV, one influenza A virus, one rhinovirus, and 13 HERV-K113 detections. Of these, one sample with BKV, three with HHV-7, one with EBV, and all HERV-K113 were confirmed by specific PCR. The requirement for Intensive Care Unit admission was associated with HCSS detections.

CONCLUSION

This study highlights HCSS as a powerful tool for the investigation of undiagnosed cases of M/ME. Data generated must be carefully analysed and reasonable precautions must be taken before establishing association of clinical features with unexpected or novel virus findings.

Torque Teno Virus DNA Load in Blood as an Immune Status Biomarker in Adult Hematological Patients: The State of the Art and Future Prospects

A solid body of scientific evidence supports the assumption that Torque teno virus (TTV) DNA load in the blood compartment may behave as a biomarker of immunosuppression in solid organ transplant recipients; in this clinical setting, high or increasing TTV DNA levels precede the occurrence of infectious complications, whereas the opposite anticipates the development of acute rejection. The potential clinical value of the TTV DNA load in blood to infer the risk of opportunistic viral infection or immune-related (i.e., graft vs. host disease) clinical events in the hematological patient, if any, remains to be determined. In fact, contradictory data have been published on this matter in the allo-SCT setting. Studies addressing this topic, which we review and discuss herein, are highly heterogeneous as regards design, patient characteristics, time points selected for TTV DNA load monitoring, and PCR assays used for TTV DNA quantification. Moreover, clinical outcomes are often poorly defined. Prospective, ideally multicenter, and sufficiently powered studies with well-defined clinical outcomes are warranted to elucidate whether TTV DNA load monitoring in blood may be of any clinical value in the management of hematological patients.

Case report of pediatric TTMV-related acute promyelocytic leukemia with central nervous system infiltration and rapid accumulation of RARA-LBD mutations

TTMV::RARA is a recently reported fusion gene associated with acute promyelocytic leukemia (APL), caused by the integration of torque teno mini virus (TTMV) genomic fragments into the second intron of the RARA gene. Currently, there have been only six documented cases, with clinical presentations showing significant variability. Although initial responses to all-trans retinoic acid (ATRA) treatment may be observed in patients with TTMV::RARA-APL, the overall prognosis remains unfavorable among infrequent reported cases. This article presents a pediatric case that manifested as PML::RARA-negative APL with central nervous system involvement at onset. The patient experienced both intramedullary and extramedullary relapse one year after undergoing allogeneic hematopoietic stem cell transplantation. Upon identification as TTMV::RARA-APL and subsequent administration of two rounds of ATRA-based treatment, the patient rapidly developed multiple RARA ligand-binding domain mutations and demonstrated extensive resistance to ATRA and various other therapeutic interventions. Additionally, the patient experienced ARID1A mutant clone expansion and progressed MYC-targeted gene activation. This case represents the first documentation of extramedullary involvement at both the initial diagnosis and relapse stages, emphasizing the intricate clinical features and challenges associated with the rapid accumulation of multiple ATRA-resistant mutations in TTMV::RARA-APL, characterizing it as a distinct and complex sub-entity of atypical APL.

The mysterious anelloviruses: investigating its role in human diseases

Anelloviruses (AVs) that infect the human population are members of the Anelloviridae family. They are widely distributed in human populations worldwide. Torque teno virus (TTV) was the first virus of this family to be identified and is estimated to be found in the serum of 80-90% of the human population. Sometime after the identification of TTV, Torque teno mini virus (TTMV) and Torque teno midi virus (TTMDV) were also identified and classified in this family. Since identifying these viruses, have been detected in various types of biological fluids of the human body, including blood and urine, as well as vital organs such as the liver and kidney. They can be transmitted from person to person through blood transfusions, fecal-oral contact, and possibly sexual intercourse. Recent studies on these newly introduced viruses show that although they are not directly related to human disease, they may be indirectly involved in initiating or exacerbating some human population-related diseases and viral infections. Among these diseases, we can mention various types of cancers, immune system diseases, viral infections, hepatitis, and AIDS. Also, they likely use the microRNAs (miRNAs) they encode to fulfill this cooperative role. Also, in recent years, the role of proliferation and their viral load, especially TTV, has been highlighted to indicate the immune system status of immunocompromised people or people who undergo organ transplants. Here, we review the possible role of these viruses in diseases that target humans and highlight them as important viruses that require further study. This review can provide new insights to researchers.

Microbial transmission in the social microbiome and host health and disease

Although social interactions are known to drive pathogen transmission, the contributions of socially transmissible host-associated mutualists and commensals to host health and disease remain poorly explored. We use the concept of the social microbiome-the microbial metacommunity of a social network of hosts-to analyze the implications of social microbial transmission for host health and disease. We investigate the contributions of socially transmissible microbes to both eco-evolutionary microbiome community processes (colonization resistance, the evolution of virulence, and reactions to ecological disturbance) and microbial transmission-based processes (transmission of microbes with metabolic and immune effects, inter-specific transmission, transmission of antibiotic-resistant microbes, and transmission of viruses). We consider the implications of social microbial transmission for communicable and non-communicable diseases and evaluate the importance of a socially transmissible component underlying canonically non-communicable diseases. The social transmission of mutualists and commensals may play a significant, under-appreciated role in the social determinants of health and may act as a hidden force in social evolution.

Characterization of a new anellovirus species infecting an ocelot (Leopardus pardalis) in Brazil

A complete genome of the first anellovirus infecting the wild felid Leopardus pardalis (ocelot) and a partial genome were assembled and annotated through high-throughput sequencing protocols followed by Sanger sequencing validation. The full-length virus obtained comprises 2,003 bp, while the partial genome comprises 1,224 bp. Phylogenetic analysis grouped these two sequences in two distinct clusters related to previously described Felidae anelloviruses. The ORF1 of the partial genome was identified as a new species provisionally called Torque teno ocelot virus, with 53.6% identity with its sister lineage. The complete genome was inferred as a new representative of the Torque teno felid virus 3 species, with 73.28% identity to the closest reference. This study expands known virus diversity and the host span of anelloviruses.

Respiratory torque teno virus load at emergency department visit predicts intensive care unit admission of SARS-CoV-2 infected patients

Accurate prediction of COVID-19 severity remains a challenge. Torque teno virus (TTV), recognized as a surrogate marker of functional immunity in solid organ transplant recipients, holds the potential for assessing infection outcomes. We investigated whether quantifying TTV in nasopharyngeal samples upon emergency department (ED) admission could serve as an early predictor of COVID-19 severity. Retrospective single-center study in the ED of Saint-Louis Hospital in Paris, France. TTV DNA was quantified in nasopharyngeal swab samples collected for SARS-CoV-2 testing. Among 295 SARS-CoV-2 infected patients, 92 returned home, 160 were admitted to medical wards, and 43 to the intensive care unit (ICU). Elevated TTV loads were observed in ICU patients (median: 3.02 log copies/mL, interquartile range [IQR]: 2.215-3.825), exceeding those in discharged (2.215, [0; 2.962]) or hospitalized patients (2.24, [0; 3.29]) (p = 0.006). Multivariate analysis identified diabetes, obesity, hepatitis, fever, dyspnea, oxygen requirement, and TTV load as predictors of ICU admission. A 2.91 log10 copies/mL TTV threshold independently predicted ICU admission. Nasopharyngeal TTV quantification in SARS-CoV-2 infected patients is linked to the likelihood of ICU admission and might reflect respiratory immunosuppression.

Exploring the relationship between anellovirus load and clinical variables in hospitalized COVID-19 patients: Implications for immune activation and inflammation

Objectives

Anelloviruses have been linked with host-immunocompetence and inflammation. Here, we studied the anellovirus load in hospitalized COVID-19 patients.

Methods

We collected samples of patients recruited in the DAWN-Plasma trial that received convalescent plasma (CP) therapy (four plasma units) combined with standard of care (SOC) or SOC of alone. Plasma samples were collected on day 0 and 6 of hospitalization and we quantified anellovirus load. With multivariate models, clinical variables were associated with changes in anellovirus load.

Results

Samples were collected on day 0 and 6 of 150 patients (103 CP + SOC and 47 SOC). Anellovirus load was higher on day 0 compared to day 6 and we found a significant drop in SOC patients. Patients receiving immunosuppressive drug had a lower anellovirus load (coefficient: 1.021, 95% confidence interval [CI] 0.270-1.772, P = 0.008), while patients admitted to the emergency room displayed a higher abundance on day 0 (1.308, 95% CI 0.443-2.173, P = 0.003). Unspecific markers of inflammation and organ damage, D-dimer (0.001, 95% CI <0.001-0.001, P = 0.001) and lactate dehydrogenase (0.002, 95% CI 0.001-0.004, P = 0.044), were positively associated with anellovirus load. Finally, anellovirus load on day 0 (-39.9, 95% CI -75.72 to -4.27, P = 0.029) was negatively associated with SARS-CoV-2 antibody response on day.

Conclusion

The results showed associations between clinical variables and anellovirus load in COVID-19 patients. Many variables share properties related to host immunocompetence or inflammation. Therefore, we expect that anellovirus abundance displays the net state of immune activation.

Natural recombination of the torque teno canis virus within the ORF1, -2, and -3 genes and the untranslated region

The torque teno canis virus (TTCaV) was first reported in 2001 and it shares similarities with the known Torque teno virus (TTV) in terms of genomic organization and putative transcriptional features. It is a single-stranded DNA virus characterized by high rates of recombination and nucleotide substitution, like RNA viruses. Studies reported recombination events in torque teno virus; however, there is limited reporting of TTCaV reorganization events. This study screened fecal samples from domestic dogs in Henan Province. There was a positivity rate of 16.5% (19/115) for TTCaV. Four nearly complete TTCaV genomes, namely Canine/HeNan/4, 5, 6, and 13/2019, were obtained from the 19 positive fecal samples, whose genome sequence was obtained using gap-filling PCR. Sequence analysis revealed two unique amino acid mutation sites in the TTCaV strains, K278Q (compared with the first isolate Cf-TTV10 in Japan) and V/L268I (compared with the TTCaV strain from southern China). Subsequently, 17 near full-length TTCaV genome sequences were subjected to phylogenetic and recombination detection program analyzes. We obtained evidence supporting recombination events in the Chinese TTCaV strains. These findings suggest that mutation and recombination occurred in the three individual gene segments (ORF1, ORF2, ORF3) and the untranslated region, an area of major recombination in the Chinese TTCaV strain GX265 genome. Interestingly, the TTCaV strain (Canine/HeNan/6/2019) was a major parent involved in the genetic recombination of the GX265 strain. This study provides insights into the genetic variability and evolution of TTCaV.

Characterization of Diverse Anelloviruses, Cressdnaviruses, and Bacteriophages in the Human Oral DNA Virome from North Carolina (USA)

The diversity of viruses identified from the various niches of the human oral cavity-from saliva to dental plaques to the surface of the tongue-has accelerated in the age of metagenomics. This rapid expansion demonstrates that our understanding of oral viral diversity is incomplete, with only a few studies utilizing passive drool collection in conjunction with metagenomic sequencing methods. For this pilot study, we obtained 14 samples from healthy staff members working at the Duke Lemur Center (Durham, NC, USA) to determine the viral diversity that can be identified in passive drool samples from humans. The complete genomes of 3 anelloviruses, 9 cressdnaviruses, 4 Caudoviricetes large bacteriophages, 29 microviruses, and 19 inoviruses were identified in this study using high-throughput sequencing and viral metagenomic workflows. The results presented here expand our understanding of the vertebrate-infecting and microbe-infecting viral diversity of the human oral virome in North Carolina (USA).

Spatial scale influences the distribution of viral diversity in the eukaryotic virome of the mosquito Culex pipiens

Our knowledge of the diversity of eukaryotic viruses has recently undergone a massive expansion. This diversity could influence host physiology through yet unknown phenomena of potential interest to the fields of health and food production. However, the assembly processes of this diversity remain elusive in the eukaryotic viromes of terrestrial animals. This situation hinders hypothesis-driven tests of virome influence on host physiology. Here, we compare taxonomic diversity between different spatial scales in the eukaryotic virome of the mosquito Culex pipiens. This mosquito is a vector of human pathogens worldwide. The experimental design involved sampling in five countries in Africa and Europe around the Mediterranean Sea and large mosquito numbers to ensure a thorough exploration of virus diversity. A group of viruses was found in all countries. This core group represented a relatively large and diverse fraction of the virome. However, certain core viruses were not shared by all host individuals in a given country, and their infection rates fluctuated between countries and years. Moreover, the distribution of coinfections in individual mosquitoes suggested random co-occurrence of those core viruses. Our results also suggested differences in viromes depending on geography, with viromes tending to cluster depending on the continent. Thus, our results unveil that the overlap in taxonomic diversity can decrease with spatial scale in the eukaryotic virome of C. pipiens. Furthermore, our results show that integrating contrasted spatial scales allows us to identify assembly patterns in the mosquito virome. Such patterns can guide future studies of virome influence on mosquito physiology.

Expansion of Betatorquevirus and/or Gammatorquevirus in Patients with Severe Clinical Outcomes of the Liver Diseases

Anellovirus (AV) is a ubiquitous virus in the human population. Individuals can be infected with multiple AV genera and species to form a heterogeneous repertoire, termed the anellome. Using advanced methods, we examined the anellomes from 12 paired serum and liver samples, as well as 2701 subjects with different clinical diagnoses. Overall, anellomes are remarkably individualized, with significant among-group differences (Kruskal-Wallis test p = 6.6 × 10-162 for richness and p = 7.48 × 10-162 for Shannon entropy). High dissimilarity scores (beta diversity) were observed between patient groups, except for paired serum and liver samples. At the population level, the relative abundance of combinational AV genus Betatorquevirus (torque teno mini viruses, TTMV), and Gammatorquevirus (torque teno midi viruses, TTMDV) exhibited an exponential distribution with a low bound point at 32%. Defined by this value, the AV TTMV/TTMDV-expanded anellome was significantly enriched among patients with acute liver failure (31.7%) and liver transplantation (40.7%), compared with other patient groups (χ2 test: p = 4.1 × 10-8-3.2 × 10-3). Therefore, anellome heterogeneity may be predictive of clinical outcomes in certain diseases, such as liver disease. The consistency of anellome between paired serum and liver samples indicates that a liquid biopsy approach would be suitable for longitudinal studies to clarify the causality of the AV TTMV/TTMDV-expanded anellome in the outcomes of liver disease.

SCANellome: Analysis of the Genomic Diversity of Human and Non-Human Primate Anelloviruses from Metagenomics Data

Anelloviruses are extremely prevalent in the human population and are considered to be commensal parts of the human virome. The best-known member in humans is the Torque teno virus. Recent metagenomic next-generation sequencing investigations have helped reveal the considerable number of species and genotypes from the same genus that can be co-detected within a single individual and that this diversity increases as a function of age during the first months/years of life. As a result, to date, the bioinformatics analysis of this genetic diversity remains complex and constraining for researchers. Here, we present SCANellome, a user-friendly tool to investigate the anellome composition at the genus, species, and genotype levels of samples from metagenomics data generated by the Illumina and Nanopore platforms. SCANellome is based on an in-house up-to-date database that includes all human and non-human primate anellovirus reference sequences available on GenBank and meets the latest classification criteria established by the International Committee on Taxonomy of Viruses.

Genomic Characteristics of Feline Anelloviruses Isolated from Domestic Cats in Shanghai, China

Viral metagenomics techniques allow the high-throughput discovery of possible pathogens carried by companion animals from their feces and other excreta. In this study, the viral metagenomics of 22 groups of fecal samples from domestic cats revealed a high prevalence of feline anelloviruses (FcTTV) infection in domestic cats in Shanghai, China. Serum samples from 30 cat individuals were further detected by polymerase chain reaction, and an average positive rate of 36.67% (11/30) of FcTTV infection was found. Next, the full-length sequences of five Shanghai FcTTV variants were obtained and submitted to GenBank with access numbers OP186140 to OP186144. Phylogenetic analysis indicates that the Shanghai FcTTV variants have relatively consistent genomic characteristics, with two variants from Zhejiang 2019 and one variant from the Czech Republic 2010. The recombination event analysis of the variants showed that one variant (OP186141_SH-02) had a primary parental sequence derived from a variant (KM229764) from the Czech Republic in 2010, while the secondary parental sequence was derived from OP186140_SH-01. The results revealed that FcTTV infection is prevalent in domestic cats and that the use of viral metagenomics to rapidly identify some infecting viruses whose hosts lack clinical features would be an effective approach.

Characterization of the Human Blood Virome in Iranian Multiple Transfused Patients

Blood transfusion safety is an essential element of public health. Current blood screening strategies rely on targeted techniques that could miss unknown or unexpected pathogens. Recent studies have demonstrated the presence of a viral community (virobiota/virome) in the blood of healthy individuals. Here, we characterized the blood virome in patients frequently exposed to blood transfusion by using Illumina metagenomic sequencing. The virome of these patients was compared to viruses present in healthy blood donors. A total number of 155 beta-thalassemia, 149 hemodialysis, and 100 healthy blood donors were pooled with five samples per pool. Members of the Anelloviridae and Flaviviridae family were most frequently observed. Interestingly, samples of healthy blood donors harbored traces of potentially pathogenic viruses, including adeno-, rota-, and Merkel cell polyomavirus. Viruses of the Anelloviridae family were most abundant in the blood of hemodialysis patients and displayed a higher anellovirus richness. Pegiviruses (Flaviviridae) were only observed in patient populations. An overall trend of higher eukaryotic read abundance in both patient groups was observed. This might be associated with increased exposure through blood transfusion. Overall, the findings in this study demonstrated the presence of various viruses in the blood of Iranian multiple-transfused patients and healthy blood donors.

Human Anelloviruses: Influence of Demographic Factors, Recombination, and Worldwide Diversity

Anelloviruses represent the major and most diverse component of the healthy human virome, referred to as the anellome. In this study, we determined the anellome of 50 blood donors, forming two sex- and age-matched groups. Anelloviruses were detected in 86% of the donors. The number of detected anelloviruses increased with age and was approximately twice as high in men as in women. A total of 349 complete or nearly complete genomes were classified as belonging to torque teno virus (TTV), torque teno mini virus (TTMV), and torque teno midi virus (TTMDV) anellovirus genera (197, 88, and 64 sequences, respectively). Most donors had intergenus (69.8%) or intragenus (72.1%) coinfections. Despite the limited number of sequences, intradonor recombination analysis showed 6 intragenus recombination events in ORF1. As thousands of anellovirus sequences have been described recently, we finally analyzed the global diversity of human anelloviruses. Species richness and diversity were close to saturation in each anellovirus genus. Recombination was found to be the main factor promoting diversity, although its effect was significantly lower in TTV than in TTMV and TTMDV. Overall, our results suggest that differences in diversity between genera may be caused by variations in the relative contribution of recombination. IMPORTANCE Anelloviruses are the most common human infectious viruses and are considered essentially harmless. Compared to other human viruses, they are characterized by enormous diversity, and recombination is suggested to play an important role in their diversification and evolution. Here, by analyzing the composition of the plasma anellome of 50 blood donors, we find that recombination is also a determinant of viral evolution at the intradonor level. On a larger scale, analysis of anellovirus sequences currently available in databases shows that their diversity is close to saturation and differs among the three human anellovirus genera and that recombination is the main factor explaining this intergenus variability. Global characterization of anellovirus diversity could provide clues about possible associations between certain virus variants and pathologies, as well as facilitate the implementation of unbiased PCR-based detection protocols, which may be relevant for using anelloviruses as endogenous markers of immune status.

Evolution of anelloviruses from a circovirus-like ancestor through gradual augmentation of the jelly-roll capsid protein

Anelloviruses are highly prevalent in diverse mammals, including humans, but so far have not been linked to any disease and are considered to be part of the 'healthy virome'. These viruses have small circular single-stranded DNA (ssDNA) genomes and encode several proteins with no detectable sequence similarity to proteins of other known viruses. Thus, anelloviruses are the only family of eukaryotic ssDNA viruses currently not included in the realm Monodnaviria. To gain insights into the provenance of these enigmatic viruses, we sequenced more than 250 complete genomes of anelloviruses from nasal and vaginal swab samples of Weddell seal (Leptonychotes weddellii) from Antarctica and a fecal sample of grizzly bear (Ursus arctos horribilis) from the USA and performed a comprehensive family-wide analysis of the signature anellovirus protein ORF1. Using state-of-the-art remote sequence similarity detection approaches and structural modeling with AlphaFold2, we show that ORF1 orthologs from all Anelloviridae genera adopt a jelly-roll fold typical of viral capsid proteins (CPs), establishing an evolutionary link to other eukaryotic ssDNA viruses, specifically, circoviruses. However, unlike CPs of other ssDNA viruses, ORF1 encoded by anelloviruses from different genera display remarkable variation in size, due to insertions into the jelly-roll domain. In particular, the insertion between β-strands H and I forms a projection domain predicted to face away from the capsid surface and function at the interface of virus-host interactions. Consistent with this prediction and supported by recent experimental evidence, the outermost region of the projection domain is a mutational hotspot, where rapid evolution was likely precipitated by the host immune system. Collectively, our findings further expand the known diversity of anelloviruses and explain how anellovirus ORF1 proteins likely diverged from canonical jelly-roll CPs through gradual augmentation of the projection domain. We suggest assigning Anelloviridae to a new phylum, 'Commensaviricota', and including it into the kingdom Shotokuvirae (realm Monodnaviria), alongside Cressdnaviricota and Cossaviricota.

Plasma virome dynamics in chronic hepatitis B virus infected patients

The virome remains an understudied domain of the human microbiome. The role of commensal viruses on the outcome of infections with known pathogens is not well characterized. In this study we aimed to characterize the longitudinal plasma virome dynamics in chronic hepatitis B virus (HBV) infected patients. Eighty-five longitudinal plasma samples were collected from 12 chronic HBV infected individuals that were classified in the four stages of HBV infection. The virome was characterized with an optimized viral extraction protocol and deep-sequenced on a NextSeq 2500 platform. The plasma virome was primarily composed of members of the Anello- Flavi-, and Hepadnaviridae (HBV) families. The virome structure and dynamics did not correlate with the different stages of chronic HBV infection nor with the administration of antiviral therapy. We observed a higher intrapersonal similarity of viral contigs. Genomic analysis of viruses observed in multiple timepoint demonstrated the presence of a dynamic community. This study comprehensively assessed the blood virome structure in chronic HBV infected individuals and provided insights in the longitudinal development of this viral community.

Second-Generation SARS-CoV-2 Recombinants: Lessons from Other Viruses

RNA viruses have developed notable strategies to evolve and escape host immunity [...].

Longitudinal Detection of Twenty DNA and RNA Viruses in Allogeneic Hematopoietic Stem Cell Transplant Recipients Plasma

Metagenomics revealed novel and routinely overlooked viruses, representing sources of unrecognized infections after allogeneic hematopoietic stem cell transplantation (allo-HSCT). We aim to describe DNA and RNA virus prevalence and kinetics in allo-HSCT recipients' plasma for one year post HSCT. We included 109 adult patients with first allo-HSCT from 1 March 2017 to 31 January 2019 in this observational cohort study. Seventeen DNA and three RNA viral species were screened with qualitative and/or quantitative r(RT)-PCR assays using plasma samples collected at 0, 1, 3, 6, and 12 months post HSCT. TTV infected 97% of patients, followed by HPgV-1 (prevalence: 26-36%). TTV (median 3.29 × 10⁵ copies/mL) and HPgV-1 (median 1.18 × 10⁶ copies/mL) viral loads peaked at month 3. At least one Polyomaviridae virus (BKPyV, JCPyV, MCPyV, HPyV6/7) was detected in >10% of patients. HPyV6 and HPyV7 prevalence reached 27% and 12% at month 3; CMV prevalence reached 27%. HSV, VZV, EBV, HHV-7, HAdV and B19V prevalence remained <5%. HPyV9, TSPyV, HBoV, EV and HPg-V2 were never detected. At month 3, 72% of patients had co-infections. TTV and HPgV-1 infections were highly prevalent. BKPyV, MCPyV and HPyV6/7 were frequently detected relative to classical culprits. Further investigation is needed into associations between these viral infections and immune reconstitution or clinical outcomes.

Herpesviruses, polyomaviruses, parvoviruses, papillomaviruses, and anelloviruses in vestibular schwannoma

Etiology of vestibular schwannoma (VS) is unknown. Viruses can infect and reside in neural tissues for decades, and new viruses with unknown tumorigenic potential have been discovered. The presence of herpesvirus, polyomavirus, parvovirus, and anellovirus DNA was analyzed by quantitative PCR in 46 formalin-fixed paraffin-embedded VS samples. Five samples were analyzed by targeted next-generation sequencing. Viral DNA was detected altogether in 24/46 (52%) tumor samples, mostly representing anelloviruses (46%). Our findings show frequent persistence of anelloviruses, considered normal virome, in VS. None of the other viruses showed an extensive presence, thereby suggesting insignificant role in VS.

Comprehensive analysis of circulating viral DNA in maternal plasma at population-scale using low-pass whole-genome sequencing

As the most readily adopted molecular screening test, low-pass WGS of maternal plasma cell-free DNA for aneuploidy detection generates a vast amount of genomic data. This large-scale method also allows for high-throughput virome screening. NIPT sequencing data, yielding 6.57 terabases of data from 187.8 billion reads, from 12,951 pregnant Turkish women was used to investigate the prevalence and abundance of viral DNA in plasma. Among the 22 virus sequences identified in 12% of participants were human papillomavirus, herpesvirus, betaherpesvirus and anellovirus. We observed a unique pattern of circulating viral DNA with a high prevalence of papillomaviruses. The prevalence of herpesviruses/anellovirus was similar among Turkish, European and Dutch populations. Hepatitis B prevalence was remarkably low in Dutch, European and Turkish populations, but higher in China. WGS data revealed that herpesvirus/anelloviruses are naturally found in European populations. This represents the first comprehensive research on the plasma virome of pregnant Turkish women.

Dynamics of Human Anelloviruses in Plasma and Clinical Outcomes Following Kidney Transplantation

BACKGROUND

Torque teno virus, the major member of the genus Alphatorquevirus , is an emerging biomarker of the net state of immunosuppression after kidney transplantation. Genetic diversity constitutes a main feature of the Anelloviridae family, although its posttransplant dynamics and clinical correlates are largely unknown.

METHODS

The relative abundance of Alphatorquevirus , Betatorquevirus , and Gammatorquevirus genera was investigated by high-throughput sequencing in plasma specimens obtained at various points during the first posttransplant year (n = 91 recipients). Total loads of all members of the Anelloviridae family were also quantified by an "in-house" polymerase chain reaction assay targeting conserved DNA sequences (n = 195 recipients). In addition to viral kinetics, clinical study outcomes included serious infection, immunosuppression-related adverse event (opportunistic infection and cancer)' and acute rejection.

RESULTS

Alphatorquevirus DNA was detected in all patients at every point, with an increase from pretransplantation to month 1. A variable proportion of recipients had detectable Betatorquevirus and Gammatorquevirus at lower frequencies. At least 1 change in the predominant genus (mainly as early transition to Alphatorquevirus predominance) was shown in 35.6% of evaluable patients. Total anelloviruses DNA levels increased from baseline to month 1, to peak by month 3 and decrease thereafter, and were higher in patients treated with T-cell depleting agents. There was a significant albeit weak-to-moderate correlation between total anelloviruses and TTV DNA levels. No associations were found between the predominant Anelloviridae genus or total anelloviruses DNA levels and clinical outcomes.

CONCLUSIONS

Our study provides novel insight into the evolution of the anellome after kidney transplantation.

Anellovirus evolution during long-term chronic infection

Human anelloviruses (AVs) are extremely genetically diverse, are widespread in the human population, and cause chronic infections. However, the evolutionary dynamics of AVs within single hosts is currently unknown, and it is unclear whether these changes have an implication on the long-term persistence of AVs in the host. Here, we assessed the evolutionary dynamics of six AV lineages during 30 years of chronic infection at single host resolution. The total number of substitutions and the number of variable sites increased over time. However, not all substitutions reached population fixation, showing that AV lineages form heterogeneous swarms within the host. Most substitutions occurred within a hypervariable region (HVR) located between nucleotide positions 800 and 1,300 of ORF1, which is known to be located within the spike domain. Different regions of the ORF1 gene undergo either positive or negative selection pressure. Sites under strong diversifying selection pressure were detected in the HVR, while the majority of the sites under purifying selection were detected outside this region. The HVR may play the role of an immunological decoy that prevents antibodies from binding to more vulnerable parts of ORF1. Moreover, the frequent substitutions in this region may increase the chances of AV particles escaping immune recognition.

Respiratory virome profiles reflect antiviral immune responses

BACKGROUND

From early life, respiratory viruses are implicated in the development, exacerbation and persistence of respiratory conditions such as asthma. Complex dynamics between microbial communities and host immune responses shape immune maturation and homeostasis, influencing health outcomes. We evaluated the hypothesis that the respiratory virome is linked to systemic immune responses, using peripheral blood and nasopharyngeal swab samples from preschool-age children in the PreDicta cohort.

METHODS

Peripheral blood mononuclear cells from 51 children (32 asthmatics and 19 healthy controls) participating in the 2-year multinational PreDicta cohort were cultured with bacterial (Bacterial-DNA, LPS) or viral (R848, Poly:IC, RV) stimuli. Supernatants were analysed by Luminex for the presence of 22 relevant cytokines. Virome composition was obtained using untargeted high throughput sequencing of nasopharyngeal samples. The metagenomic data were used for the characterization of virome profiles and the presence of key viral families (Picornaviridae, Anelloviridae, Siphoviridae). These were correlated to cytokine secretion patterns, identified through hierarchical clustering and principal component analysis.

RESULTS

High spontaneous cytokine release was associated with increased presence of Prokaryotic virome profiles and reduced presence of Eukaryotic and Anellovirus profiles. Antibacterial responses did not correlate with specific viral families or virome profile; however, low antiviral responders had more Prokaryotic and less Eukaryotic virome profiles. Anelloviruses and Anellovirus-dominated profiles were equally distributed among immune response clusters. The presence of Picornaviridae and Siphoviridae was associated with low interferon-λ responses. Asthma or allergy did not modify these correlations.

CONCLUSION

Antiviral cytokine responses at a systemic level reflect the upper airway virome composition. Individuals with low innate interferon responses have higher abundance of Picornaviruses (mostly Rhinoviruses) and bacteriophages. Bacteriophages, particularly Siphoviridae, appear to be sensitive sensors of host antimicrobial capacity, while Anelloviruses are not correlated with TLR-induced immune responses.

TTV and other anelloviruses: The astonishingly wide spread of a viral infection

The broad family of viruses known as anelloviruses (AV) infects both humans and numerous animal species. They have a tiny, covalently closed single-stranded DNA genome and the astonishing capacity to infect a very high percentage of healthy and ill people with chronic infections that could last a lifetime. AV, and particularly the prototype Torquetenovirus, have established a successful interaction with the host's immune system and the rate at which they replicate is a gauge to measure overall immune function, even though many aspects of their life cycle and pathogenesis are still poorly understood.

Control of Human Anelloviruses by Cytosine to Uracil Genome Editing

Anelloviruses are the most common viruses infecting humans. Every human carries a nonpathogenic personal anellovirus virome (anellome), yet it is unknown which mechanisms contribute to its stability. Here, we assessed the dynamics and impact of a host antiviral defense mechanism-cytidine deaminase activity leading to C to U editing in anelloviruses-on the stability of the anellome. We investigated anellome sequence data obtained from serum samples collected every 6 months from two healthy subjects followed for more than 30 years. The subjects were infected by a total of 64 anellovirus lineages. Minus-stranded C to U editing was observed in lineages belonging to the Alpha-, Beta-, and Gammatorquevirus genera. The edited genomes were present within virus particles, therefore editing must have occurred at the late stages of the virus life cycle. Editing was favored by 5'-TC contexts in the virus genome, indicating that apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like, catalytic subunit 3 or A3 (APOBEC3) proteins are involved. Within a lineage, mutational dynamics varied over time and few fixations of mutations were detected, indicating that C to U editing is a dead end for a virus genome. We detected an editing coldspot in the GC-rich regions, suggesting that the GC-rich region is crucial for genome packaging, since only packaged virus particles were included in the analysis. Finally, we noticed a lineage-specific reduced concentration after an editing event, yet no clearance. In conclusion, cytidine deaminase activity does not clear anelloviruses, nor does it play a major role in virus evolution, but it does contribute to the stability of the anellome. IMPORTANCE Despite significant attention on anellovirus research, the interaction between the anellovirus virome and the human host remains unknown. We show the dynamics of APOBEC3-mediated cytidine deaminase activity on anelloviruses during a 30-year period of chronic infection and postulate that this antiviral mechanism controls anelloviruses. These results expand our knowledge of anellovirus-host interactions, which may be important for the design of gene therapies.

Comprehensive profiling of antibody responses to the human anellome using programmable phage display

Anelloviruses represent a major constituent of the commensal human virome; however, little is known about their immunobiology. Here, we present "AnelloScan," a T7 phage library representing the open reading frame 1 (ORF1), ORF2, ORF3, and torque teno virus (TTV)-derived apoptosis-inducing protein (TAIP) sequences of more than 800 human anelloviruses and profile the antibody reactivities of serum samples from a cross-sectional cohort of 156 subjects by using phage-immunoprecipitation sequencing (PhIP-Seq). A majority of anellovirus peptides are not reactive in any of the subjects tested (n = ∼28,000; ∼85% of the library). Antibody-reactive peptides are largely restricted to the C-terminal region of the capsid protein ORF1. Moreover, using a longitudinal cohort of matched blood-transfusion donors and recipients, we find that most transmitted anelloviruses do not elicit a detectable antibody reactivity in the recipient and that the remainder elicit delayed responses appearing ∼100-150 days after transfusion.

A Unified Model of Age-Related Cardiovascular Disease

Despite progress in biomedical technologies, cardiovascular disease remains the main cause of mortality. This is at least in part because current clinical interventions do not adequately take into account aging as a driver and are hence aimed at suboptimal targets. To achieve progress, consideration needs to be given to the role of cell aging in disease pathogenesis. We propose a model unifying the fundamental processes underlying most age-associated cardiovascular pathologies. According to this model, cell aging, leading to cell senescence, is responsible for tissue changes leading to age-related cardiovascular disease. This process, occurring due to telomerase inactivation and telomere attrition, affects all components of the cardiovascular system, including cardiomyocytes, vascular endothelial cells, smooth muscle cells, cardiac fibroblasts, and immune cells. The unified model offers insights into the relationship between upstream risk factors and downstream clinical outcomes and explains why interventions aimed at either of these components have limited success. Potential therapeutic approaches are considered based on this model. Because telomerase activity can prevent and reverse cell senescence, telomerase gene therapy is discussed as a promising intervention. Telomerase gene therapy and similar systems interventions based on the unified model are expected to be transformational in cardiovascular medicine.

Impact of polymorphisms in genes orchestrating innate immune responses on replication kinetics of Torque teno virus after kidney transplantation

Background: Torque teno virus (TTV) DNAemia has been proposed as a surrogate marker of immunosuppression after kidney transplantation (KT), under the assumption that the control of viral replication is mainly exerted by T-cell-mediated immunity. However, Tthe impact on post-transplant TTV kinetics of single genetic polymorphisms (SNPs) in genes orchestrating innate responses remains unknown. We aimed to characterize the potential association between 14 of these SNPs and TTV DNA levels in a single-center cohort of KT recipients.

Methods: Plasma TTV DNAemia was quantified by real-time PCR in 221 KT recipients before transplantation (baseline) and regularly through the first 12 post-transplant months. We performed genotyping of the following SNPs: CTLA4 (rs5742909, rs231775), TLR3 (rs3775291), TLR9 (rs5743836, rs352139), CD209 (rs735240, rs4804803), IFNL3 (rs12979860, rs8099917), TNF (rs1800629), IL10 (rs1878672, rs1800872), IL12B (rs3212227) and IL17A (rs2275913).

Results: The presence of the minor G allele of CD209 (rs4804803) in the homozygous state was associated with undetectable TTV DNAemia at the pre-transplant assessment (adjusted odds ratio: 36.96; 95% confidence interval: 4.72–289.67; p-value = 0.001). After applying correction for multiple comparisons, no significant differences across SNP genotypes were observed for any of the variables of post-transplant TTV DNAemia analyzed (mean and peak values, areas under the curve during discrete periods, or absolute increments from baseline to day 15 and months 1, 3, 6 and 12 after transplantation).

Conclusion: The minor G allele of CD209 (rs4804803) seems to exert a recessive protective effect against TTV infection in non-immunocompromised patients. However, no associations were observed between the SNPs analyzed and post-transplant kinetics of TTV DNAemia. These negative results would suggest that post-transplant TTV replication is mainly influenced by immunosuppressive therapy rather than by underlying genetic predisposition, reinforcing its clinical application as a biomarker of adaptive immunity.

Integrated Immunologic Monitoring in Solid Organ Transplantation: The Road Toward Torque Teno Virus-guided Immunosuppression

Potent immunosuppressive drugs have been introduced into clinical care for solid organ transplant recipients. It is now time to guide these drugs on an individual level to optimize their efficacy. An ideal tool simultaneously detects overimmunosuppression and underimmunosuppression, is highly standardized, and is straightforward to implement into routine. Randomized controlled interventional trials are crucial to demonstrate clinical value. To date, proposed assays have mainly focused on the prediction of rejection and were based on the assessment of few immune compartments. Recently, novel tools have been introduced based on a more integrated approach to characterize the immune function and cover a broader spectrum of the immune system. In this respect, the quantification of the plasma load of a highly prevalent and apathogenic virus that might reflect the immune function of its host has been proposed: the torque teno virus (TTV). Although TTV control is driven by T cells, other major immune compartments might contribute to the hosts' response. A standardized in-house polymerase chain reaction and a conformité européenne-certified commercially available polymerase chain reaction are available for TTV quantification. TTV load is associated with rejection and infection in solid organ transplant recipients, and cutoff values for risk stratification of such events have been proposed for lung and kidney transplantation. Test performance of TTV load does not allow for the diagnosis of rejection and infection but is able to define at-risk patients. Hitherto TTV load has not been used in interventional settings, but two interventional randomized controlled trials are currently testing the safety and efficacy of TTV-guided immunosuppression.

Transmission of anelloviruses to HIV-1 infected children

Anelloviruses (AVs) are widespread in the population and infect humans at the early stage of life. The mode of transmission of AVs is still unknown, however, mother-to-child transmission, e.g., via breastfeeding, is one of the likely infection routes. To determine whether the mother-to-child transmission of AVs may still occur despite the absence of natural birth and breastfeeding, 29 serum samples from five HIV-1-positive mother and child pairs were Illumina-sequenced. The Illumina reads were mapped to an AV lineage database “Anellometrix” containing 502 distinct ORF1 sequences. Although the majority of lineages from the mother were not shared with the child, the mother and child anellomes did display a significant similarity. These findings suggest that AVs may be transmitted from mothers to their children via different routes than delivery or breastfeeding.

Discovery of two novel Torque Teno viruses in Callithrix penicillata provides insights on Anelloviridae diversification dynamics

The development of high-throughput sequencing (HTS) technologies and metagenomics protocols deeply impacted the discovery of viral diversity. Moreover, the characterization of novel viruses in the Neotropical primates (NP) is central for the comprehension of viral evolution dynamics in those hosts, due to their evolutionary proximity to Old World primates, including humans. In the present work, novel anelloviruses were detected and characterized through HTS protocols in the NP Callithrix penicillata, the common black-tufted marmoset. De novo assembly of generated sequences was carried out, and a total of 15 contigs were identified with complete Anelloviridae ORF1 gene, two of them including a flanking GC-rich region, confirming the presence of two whole novel genomes of ~3 kb. The identified viruses were monophyletic within the Epsilontorquevirus genus, a lineage harboring previously reported anelloviruses infecting hosts from the Cebidae family. The genetic divergence found in the new viruses characterized two novel species, named Epsilontorquevirus callithrichensis I and II. The phylogenetic pattern inferred for the Epsilontorquevirus genus was consistent with the topology of their host species tree, echoing a virus-host diversification model observed in other viral groups. This study expands the host span of Anelloviridae and provides insights into their diversification dynamics, highlighting the importance of sampling animal viral genomes to obtain a clearer depiction of their long-term evolutionary processes.

Genomic Diversity of Torque Teno Virus in Blood Samples from Febrile Paediatric Outpatients in Tanzania: A Descriptive Cohort Study

Torque teno virus (TTV) is considered to be an ubiquitous member of the commensal human blood virome commonly reported in mixed genotype co-infections. This study investigates the genomic diversity of TTV in blood samples from 816 febrile Tanzanian children. Metagenomic next-generation sequencing was used to screen for TTV in individual blood samples from a cohort of 816 febrile Tanzanian paediatric outpatients. For positive samples, the number of TTV species and genotypes present were evaluated. We investigate the linear relationship between individual TTV diversity and the patient age by linear regression. TTV was detected in 97.2% of sera. ORF1 analysis revealed the presence of 149 genotypes from 38 species, suggesting the presence of 13 new species. These genotypes were mostly present as co-infections with a median of 11 genotypes/subject (range: 1−71). In terms of species, we found a median of nine species/subject (range: 1−29). We further show a significant association between the diversity of co-detected TTV and the age of the subjects (p value < 0.0001). This study shows that significant TTV genomic diversity is acquired by the age of five and that this diversity tends to increase with age, which indicates a repetitive TTV acquisition during the first months/years of life.

The enigmatic roles of Anelloviridae and Redondoviridae in humans

Anelloviridae and Redondoviridae are virus families with small, circular, single-stranded DNA genomes that are common components of the human virome. Despite their small genome size of less than 5000 bases, they are remarkably successful - anelloviruses colonize over 90% of adult humans, while the recently discovered redondoviruses have been found at up to 80% prevalence in some populations. Anelloviruses are present in blood and many organs, while redondoviruses are found mainly in the ororespiratory tract. Despite their high prevalence, little is known about their biology or pathogenic potential. In this review, we discuss anelloviruses and redondoviruses and explore their enigmatic roles in human health and disease.

Plasma Virome Reveals Blooms and Transmission of Anellovirus in Intravenous Drug Users with HIV-1, HCV, and/or HBV Infections

Intravenous drug users (IDUs) are a high-risk group for HIV-1, hepatitis C virus (HCV), and hepatitis B virus (HBV) infections, which are the leading causes of death in IDUs. However, the plasma virome of IDUs and how it is influenced by above viral infections remain unclear. Using viral metagenomics, we determined the plasma virome of IDUs and its association with HIV-1, HCV, and/or HBV infections. Compared with healthy individuals, IDUs especially those with major viral infections had higher viral abundance and diversity. Anelloviridae dominated plasma virome. Coinfections of multiple anelloviruses were common, and anelloviruses from the same genus tended to coexist together. In this study, 4,487 anellovirus ORF1 sequences were identified, including 1,620 (36.1%) with less than 69% identity to any known sequences, which tripled the current number. Compared with healthy controls (HC), more anellovirus sequences were observed in neg-IDUs, and HIV-1, HCV, and/or HBV infections further expanded the sequence number in IDUs, which was characterized by the emergence of novel divergent taxons and blooms of resident anelloviruses. Pegivirus was mainly identified in infected IDUs. Five main pegivirus transmission clusters (TCs) were identified by phylogenetic analysis, suggesting a transmission link. Similar anellovirus profiles were observed in IDUs within the same TC, suggesting transmission of anellome among IDUs. Our data suggested that IDUs suffered higher plasma viral burden especially anelloviruses, which was associated with HIV-1, HCV, and/or HBV infections. Blooms in abundance and unprecedented diversity of anellovirus highlighted active evolution and replication of this virus in blood circulation, and an uncharacterized role it may engage with the host.

IMPORTANCE Virome is associated with immune status and determines or influences disease progression through both pathogenic and resident viruses. Increased viral burden in IDUs especially those with major viral infections indicated the suboptimal immune status and high infection risks of these population. Blooms in abundance and unprecedented diversity of anellovirus highlighted its active evolution and replication in the blood circulation, and sensitive response to other viral infections. In addition, transmission cluster analysis revealed the transmission link of pegivirus among IDUs, and the individuals with transmission links shared similar anellome profiles. In-depth monitoring of the plasma virome in high-risk populations is not only needed for surveillance for emerging viruses and transmission networks of major and neglected bloodborne viruses, but also important for a better understanding of commensal viruses and their role it may engage with immune system.

Torque Teno Virus Primary Infection Kinetics in Early Childhood

Human torque teno viruses (TTVs) are a diverse group of small nonenveloped viruses with circular, single-stranded DNA genomes. These elusive anelloviruses are harbored in the blood stream of most humans and have thus been considered part of the normal flora. Whether the primary infection as a rule take(s) place before or after birth has been debated. The aim of our study was to determine the time of TTV primary infection and the viral load and strain variations during infancy and follow-up for up to 7 years. TTV DNAs were quantified in serial serum samples from 102 children by a pan-TTV quantitative PCR, and the amplicons from representative time points were cloned and sequenced to disclose the TTV strain diversity. We detected an unequivocal rise in TTV-DNA prevalence, from 39% at 4 months of age to 93% at 2 years; all children but one, 99%, became TTV-DNA positive before age 4 years. The TTV-DNA quantities ranged from 5 × 10¹ to 4 × 10⁷ copies/mL, both within and between the children. In conclusion, TTV primary infections occur mainly after birth, and increase during the first two years with high intra- and interindividual variation in both DNA quantities and virus strains.

Diversity and Long-Term Dynamics of Human Blood Anelloviruses

Anelloviruses (AVs) are commensal members of the human blood virome. Even though it was estimated that over 90% of the human population carries AVs, the dynamics of the AV virome (“anellome”) are unknown. We investigated the dynamics of blood anellomes in two healthy people followed up for more than 30 years. Both subjects were positive for AVs in the majority of samples. Alphatorquevirus (torque teno virus [TTV]) was the most common genus in both subjects, followed by Betatorquevirus (torque teno minivirus [TTMV]) and Gammatorquevirus (torque teno midivirus [TTMDV]). Almost five times more lineages were found in subject 1 than in subject 2, and the anellomes differed phylogenetically. Both anellomes remained compositionally stable, and 9 out of 64 AV lineages were detected in over half of the time points. We confirmed the long-term and short-term persistence of 13 lineages by specific quantitative PCR (qPCR). AV lineages were detected in blood for over 30 years. Noticeable differences in anellome richness were found between the tested subjects, but both anellomes remained compositionally stable over time. These findings demonstrate that the human blood anellome is personal and that AV infection is chronic and potentially commensal.

IMPORTANCE Knowledge of the persistence of AVs in humans is crucial to our understanding of the nature of AV infection (chronic or acute) and the role of AV in the host. We therefore investigated the dynamics of anellovirus infection in two healthy people followed up for 30 years. Our findings suggest that the human blood anellovirus virome (anellome) remains stable and personal for decades.

Plasma Torquetenovirus (TTV) microRNAs and severity of COVID-19

Background

Torquetenovirus (TTV), a widespread anellovirus recognized as the main component of the healthy human virome, displays viremia that is highly susceptible to variations in immune competence. TTV possesses microRNA (miRNA)-coding sequences that might be involved in viral immune evasion. Among TTV-encoded miRNAs, miRNA t1a, t3b, and tth8 have been found in biological fluids. Here, the presence of TTV DNA and TTV miRNAs in the plasma of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-infected subjects was investigated to monitor the possible association with coronavirus disease 2019 (COVID-19) severity.

Methods

Detection of TTV DNA and miRNA t1a, t3b, and tth8 was investigated in plasma samples of 56 SARS-CoV-2-infected subjects with a spectrum of different COVID-19 outcomes. TTV DNA and TTV miRNAs were assessed with a universal single step real-time TaqMan PCR assay and miRNA quantitative RT-PCR miRNA assay, respectively.

Results

The TTV DNA prevalence was 59%, whereas at least one TTV miRNA was found in 94% of the patients tested. miRNA tth8 was detected in 91% of subjects, followed by miRNAs t3b (64%) and miRNAt1a (30%). Remarkably, although TTV DNA was unrelated to COVID-19 severity, miRNA tth8 was significantly associated with the degree of disease (adjusted incidence rate ratio (IRR) 2.04, 95% CI 1.14–3.63, for the subjects in the high severity group compared to those in the low severity group).

Conclusions

Our findings encourage further investigation to understand the potential role of TTV miRNAs in the different outcomes of COVID-19 at early and late stages.

Metagenomic Analysis Reveals High Abundance of Torque Teno Mini Virus in the Respiratory Tract of Children with Acute Respiratory Illness

Human Anelloviridae is a highly prevalent viral family, including three main genera—Alphatorquevirus (Torque teno virus, TTV), Betatorquevirus (Torque teno mini virus, TTMV), and Gammatorquevirus (Torque teno midi virus, TTMDV). To date, the characterization of Anelloviridae in the respiratory tract of children with acute respiratory infection (ARI) has been poorly reported and mainly focused on TTV. We performed a metagenomic analysis of eight respiratory samples collected from children with an ARI of unknown etiology (eight samples tested negative with a multiplex PCR assay, out of the 39 samples initially selected based on negative routine diagnostic testing). A total of 19 pediatric respiratory samples that tested positive for respiratory syncytial virus (RSV, n = 13) or influenza virus (n = 6) were also sequenced. Anelloviridae reads were detected in 16/27 samples, including 6/8 negative samples, 7/13 RSV samples and 3/6 influenza samples. For samples with a detection of at least one Anelloviridae genus, TTMV represented 87.1 (66.1−99.2)% of Anelloviridae reads, while TTV and TTMDV represented 0.8 (0.0−9.6)% and 0.7 (0.0−7.1)%, respectively (p < 0.001). Our findings highlight a high prevalence of TTMV in respiratory samples of children with an ARI of unknown etiology, as well as in samples with an RSV or influenza infection. Larger studies are needed to explore the role of TTMV in childhood respiratory diseases.

High-Content Screening and Computational Prediction Reveal Viral Genes That Suppress the Innate Immune Response

Suppression of the host innate immune response is a critical aspect of viral replication. Upon infection, viruses may introduce one or more proteins that inhibit key immune pathways, such as the type I interferon pathway. However, the ability to predict and evaluate viral protein bioactivity on targeted pathways remains challenging and is typically done on a single-virus or -gene basis. Here, we present a medium-throughput high-content cell-based assay to reveal the immunosuppressive effects of viral proteins. To test the predictive power of our approach, we developed a library of 800 genes encoding known, predicted, and uncharacterized human virus genes. We found that previously known immune suppressors from numerous viral families such as Picornaviridae and Flaviviridae recorded positive responses. These include a number of viral proteases for which we further confirmed that innate immune suppression depends on protease activity. A class of predicted inhibitors encoded by Rhabdoviridae viruses was demonstrated to block nuclear transport, and several previously uncharacterized proteins from uncultivated viruses were shown to inhibit nuclear transport of the transcription factors NF-κB and interferon regulatory factor 3 (IRF3). We propose that this pathway-based assay, together with early sequencing, gene synthesis, and viral infection studies, could partly serve as the basis for rapid in vitro characterization of novel viral proteins. IMPORTANCE Infectious diseases caused by viral pathogens exacerbate health care and economic burdens. Numerous viral biomolecules suppress the human innate immune system, enabling viruses to evade an immune response from the host. Despite our current understanding of viral replications and immune evasion, new viral proteins, including those encoded by uncultivated viruses or emerging viruses, are being unearthed at a rapid pace from large-scale sequencing and surveillance projects. The use of medium- and high-throughput functional assays to characterize immunosuppressive functions of viral proteins can advance our understanding of viral replication and possibly treatment of infections. In this study, we assembled a large viral-gene library from diverse viral families and developed a high-content assay to test for inhibition of innate immunity pathways. Our work expands the tools that can rapidly link sequence and protein function, representing a practical step toward early-stage evaluation of emerging and understudied viruses.

Early-Life Colonization by Anelloviruses in Infants

Anelloviruses (AVs) are found in the vast majority of the human population and are most probably part of a healthy virome. These viruses infect humans in the early stage of life, however, the characteristics of the first colonizing AVs are still unknown. We screened a collection of 107 blood samples from children between 0.4 and 64.8 months of age for the presence of three AV genera: the Alpha-, Beta- and Gammatorquevirus. The youngest child that was positive for AV was 1.2 months old, and a peak in prevalence (100% of samples positive) was reached between the twelfth and eighteenth months of life. Intriguingly, the beta- and gammatorqueviruses were detected most at the early stage of life (up to 12 months), whereas alphatorqueviruses, the most common AVs in adults, increased in prevalence in children older than 12 months. To determine whether that order of colonization may be related to oral transmission and unequal presence of AV genera in breast milk, we examined 63 breast milk samples. Thirty-two percent of the breast milk samples were positive in a qPCR detecting beta- and gammatorqueviruses, while alphatorqueviruses were detected in 10% of the samples, and this difference was significant (p = 0.00654). In conclusion, we show that beta- and gammatorqueviruses colonize humans in the first months of life and that breastfeeding could play a role in AV transmission.

Insidious Insights: Implications of viral vector engineering for pathogen enhancement

Optimizing viral vectors and their properties will be important for improving the effectiveness and safety of clinical gene therapy. However, such research may generate dual-use insights relevant to the enhancement of pandemic pathogens. In particular, reliable and generalizable methods of immune evasion could increase viral fitness sufficient to cause a new pandemic. High potential for misuse is associated with (1) the development of universal genetic elements for immune modulation, (2) specific insights on capsid engineering for antibody evasion applicable to viruses with pandemic potential, and (3) the development of computational methods to inform capsid engineering. These risks may be mitigated by prioritizing non-viral delivery systems, pharmacological immune modulation methods, non-genetic vector surface modifications, and engineering methods specific to AAV and other viruses incapable of unassisted human-to-human transmission. We recommend that computational vector engineering and the publication of associated code and data be limited to AAV until a technical solution for preventing malicious access to viral engineering tools has been established.

Within-host quantitation of anellovirus genome complexity from clinical samples

Anellovirus (AV) is a ubiquitous and diverse virus in the human population. An individual can be infected with multiple AV genera and species that form a heterogeneous repertoire, called the anellome. Due to its exceptional genetic diversity, efficient evaluation of anellome complexity remains a methodological challenge. In the current study, AV genome was first enriched from patient serum samples through two-phase rolling circle amplification. Following Illumina sequencing, anellome was analyzed with an advanced bioinformatics pipeline, including read extraction at three similarity levels, de novo assembly, species assignment, and determination of relative abundance among AV variants. The method was validated in the mock sample and then applied to 21 hepatitis C virus (HCV) patients with and without hepatocellular carcinoma (HCC). Overall, there was a large variance regarding AV richness, ranging from 2 to 51 AV species. In contrast to HCV patients without HCC, HCC incidence was associated with reduced richness (12.6 ± 14.4 vs. 35.4 ± 13.6, p = 0.001) and Shannon entropy (0.4 ± 0.34 vs. 0.61 ± 0.12, p = 0.095) at the AV species level. Interestingly, AV genus beta and gamma expanded in the anellome in 7 of 10 HCC patients. These observations shed light on the potential association between anellome and HCC incidence in patients with chronic HCV infection. The method presented here represents a valuable tool to investigate the role of anellome in human health and disease.

Plasma virome and the risk of blood-borne infection in persons with substance use disorder

There is an urgent need for innovative methods to reduce transmission of bloodborne pathogens like HIV and HCV among people who inject drugs (PWID). We investigate if PWID who acquire non-pathogenic bloodborne viruses like anelloviruses and pegiviruses might be at greater risk of acquiring a bloodborne pathogen. PWID who later acquire HCV accumulate more non-pathogenic viruses in plasma than matched controls who do not acquire HCV infection. Additionally, phylogenetic analysis of those non-pathogenic virus sequences reveals drug use networks. Here we find first in Baltimore and confirm in San Francisco that the accumulation of non-pathogenic viruses in PWID is a harbinger for subsequent acquisition of pathogenic viruses, knowledge that may guide the prioritization of the public health resources to combat HIV and HCV.

Spread of bloodborne infections, such as HCV and HIV, is a problem, particularly amongst people who inject drugs (PWID). Here, the authors describe and then confirm in observational PWID cohorts that those with more non-pathogenic viruses in plasma were more likely later to acquire HCV than PWID who had fewer of these non-pathogenic viruses.

Exploring the Diversity of the Human Blood Virome

Metagenomics is greatly improving our ability to discover new viruses, as well as their possible associations with disease. However, metagenomics has also changed our understanding of viruses in general. The vast expansion of currently known viral diversity has revealed a large fraction of non-pathogenic viruses, and offers a new perspective in which viruses function as important components of many ecosystems. In this vein, studies of the human blood virome are often motivated by the search for new viral diseases, especially those associated with blood transfusions. However, these studies have revealed the common presence of apparently non-pathogenic viruses in blood, particularly human anelloviruses and, to a lower extent, human pegiviruses (HPgV). To shed light on the diversity of the human blood virome, we subjected pooled plasma samples from 587 healthy donors in Spain to a viral enrichment protocol, followed by massive parallel sequencing. This showed that anelloviruses were clearly the major component of the blood virome and showed remarkable diversity. In total, we assembled 332 complete or near-complete anellovirus genomes, 50 of which could be considered new species. HPgV was much less frequent, but we, nevertheless, recovered 17 different isolates that we subsequently used for characterizing the diversity of this virus. In-depth investigation of the human blood virome should help to elucidate the ecology of these viruses, and to unveil potentially associated diseases.

Complex evolutionary history of felid anelloviruses

Anellovirus infections are highly prevalent in mammals, however, prior to this study only a handful of anellovirus genomes had been identified in members of the Felidae family. Here we characterise anelloviruses in pumas (Puma concolor), bobcats (Lynx rufus), Canada lynx (Lynx canadensis), caracals (Caracal caracal) and domestic cats (Felis catus). The complete anellovirus genomes (n = 220) recovered from 149 individuals were diverse. ORF1 protein sequence similarity network analysis coupled with phylogenetic analysis, revealed two distinct clusters that are populated by felid-derived anellovirus sequences, a pattern mirroring that observed for the porcine anelloviruses. Of the two-felid dominant anellovirus groups, one includes sequences from bobcats, pumas, domestic cats and an ocelot, and the other includes sequences from caracals, Canada lynx, domestic cats and pumas. Coinfections of diverse anelloviruses appear to be common among the felids. Evidence of recombination, both within and between felid-specific anellovirus groups, supports a long coevolution history between host and virus.