Abstract
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.