Supplementary Figures and Tables

metagenomics<br><br>Supplementary Figure 1. Comparison of Immune Cell Counts Between HIV-detected patients and HIV-undetected individuals. Comparison of CD4+ T-cell counts, CD8+ T-cell counts, and CD4+/CD8+ ratios between individuals with detectable HIV (HIV_detected) ad individuals with undetectable HIV (HIV_undetected). Each violin plot illustrates the distribution of the data, with the box plot representing the interquartile range (IQR), and the median indicated by the line within the box. (a) A highly significant difference in CD4+ T-cell counts is observed between the HIV_detected group and the HIV_undetected group (p&lt;0.0001). (b) CD8+ T-cell Count (cells/μL) showed no statistically significant difference (ns). (c) CD4+ T-cell /CD8+ T-cell Ratio showed significant differences between the HIV_detected group and the HIV_undetected group (p&lt;0.05).<br>Supplementary Figure 2. Evolutionary Analysis and Viral Structure. (a) Evolutionary analysis was conducted based on the Cap and RdRp genes of Picobirnavirus, with red and orange circles highlighting the viral sequences obtained in this study. (b) A phylogenetic tree was constructed based on the polymerase gene of HBV, with the viral sequences obtained in this study marked in red. On the right is the full genome structure and functional protein annotation of HBV obtained in this study. (c) Evolutionary analysis was performed based on the Cap and Rep genes of Brisavirus, with red and green circles indicating the viral sequences obtained in this study.Supplementary Figure 3. Gene Function Analysis. (a) Volcano plot of differentially expressed genes. The x-axis represents the log2-transformed fold change in gene expression (log2(Fold Change)), and the y-axis shows the negative logarithm of the p-value (-log10(p-value)). Red and blue dots indicate significantly up-regulated and down-regulated genes in the HIV patient group, respectively (filtering criteria: |log2FC| ≥ 1 and adjusted p-value (p.adjust) &lt; 0.05). (b) GSEA enrichment results plot.The most significantly enriched (up-regulated) InterPro entries are displayed in ascending order of adjusted p-value (p.adjust). Different colors represent distinct InterPro entries. (c) GSEA enrichment profiles for up- and down-regulated GO functional terms. The colored curves (red or blue) depict the running enrichment score (RES) as the analysis progresses through the ranked gene list (based on association strength with the phenotype). The peak of each curve indicates the point of maximal enrichment. Vertical tick marks represent the positions of the core enriched InterPro entries within the ranked list. The gray bars at the bottom of each plot illustrate the distribution of all genes associated with InterPro entries across the ranked list.Supplementary Table 1. Sequencing quality statisticsSupplementary Table 2. Viral read counts at the taxonomic family levelSupplementary Table 3. Viral read counts at the taxonomic genus levelSupplementary Table 4. Viral read counts at the taxonomic species levelSupplementary Table 5. Viral read counts in blood samples, stratified by taxonomic level, for LEfSe (LDA Effect Size) analysisSupplementary Table 6. Viral read counts in blood samples, stratified by taxonomic level, for LEfSe (LDA Effect Size) analysisSupplementary Table 7. Input counts for InterPro-based functional analysis (InterPro2GO; DIAMOND+MEGAN)<br><br>Supplementary Table 8. HIV viral load and immune cell levels<br><br>