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.

Host traits shape virome composition and virus transmission in wild small mammals.

Bats, rodents, and shrews are the most important animal sources of human infectious diseases. However, the evolution and transmission of viruses among them remain largely unexplored. Through the meta-transcriptomic sequencing of internal organ and fecal samples from 2,443 wild bats, rodents, and shrews sampled from four Chinese habitats, we identified 669 viruses, including 534 novel viruses, thereby greatly expanding the mammalian virome. Our analysis revealed high levels of phylogenetic diversity, identified cross-species virus transmission events, elucidated virus origins, and identified cases of invertebrate viruses in mammalian hosts. Host order and sample size were the most important factors impacting virome composition and patterns of virus spillover. Shrews harbored a high richness of viruses, including many invertebrate-associated viruses with multi-organ distributions, whereas rodents carried viruses with a greater capacity for host jumping. These data highlight the remarkable diversity of mammalian viruses in local habitats and their ability to emerge in new hosts.

Virome characterization of game animals in China reveals a spectrum of emerging pathogens.

Game animals are wildlife species traded and consumed as food and are potential reservoirs for SARS-CoV and SARS-CoV-2. We performed a meta-transcriptomic analysis of 1,941 game animals, representing 18 species and five mammalian orders, sampled across China. From this, we identified 102 mammalian-infecting viruses, with 65 described for the first time. Twenty-one viruses were considered as potentially high risk to humans and domestic animals. Civets (Paguma larvata) carried the highest number of potentially high-risk viruses. We inferred the transmission of bat-associated coronavirus from bats to civets, as well as cross-species jumps of coronaviruses from bats to hedgehogs, from birds to porcupines, and from dogs to raccoon dogs. Of note, we identified avian Influenza A virus H9N2 in civets and Asian badgers, with the latter displaying respiratory symptoms, as well as cases of likely human-to-wildlife virus transmission. These data highlight the importance of game animals as potential drivers of disease emergence.

Genomic Analyses from Non-invasive Prenatal Testing Reveal Genetic Associations, Patterns of Viral Infections, and Chinese Population History.

We analyze whole-genome sequencing data from 141,431 Chinese women generated for non-invasive prenatal testing (NIPT). We use these data to characterize the population genetic structure and to investigate genetic associations with maternal and infectious traits. We show that the present day distribution of alleles is a function of both ancient migration and very recent population movements. We reveal novel phenotype-genotype associations, including several replicated associations with height and BMI, an association between maternal age and EMB, and between twin pregnancy and NRG1. Finally, we identify a unique pattern of circulating viral DNA in plasma with high prevalence of hepatitis B and other clinically relevant maternal infections. A GWAS for viral infections identifies an exceptionally strong association between integrated herpesvirus 6 and MOV10L1, which affects piwi-interacting RNA (piRNA) processing and PIWI protein function. These findings demonstrate the great value and potential of accumulating NIPT data for worldwide medical and genetic analyses.

Factors Affecting Variation of the Human Gut Phageome.

The gut microbiome is a dense and metabolically active consortium of microorganisms and viruses located in the lower gastrointestinal tract of the human body. Bacteria and their viruses (phages) are the most abundant members of the gut microbiome. Investigating their biology and the interplay between the two is important if we are to understand their roles in human health and disease. In this review, we summarize recent advances in resolving the taxonomic structure and ecological functions of the complex community of phages in the human gut-the gut phageome. We discuss how age, diet, and geography can all have a significant impact on phageome composition. We note that alterations to the gut phageome have been observed in several diseases such as inflammatory bowel disease, irritable bowel syndrome, and colorectal cancer, and we evaluate whether these phageome changes can directly or indirectly contribute to disease etiology and pathogenesis. We also highlight how lack of standardization in studying the gut phageome has contributed to variation in reported results.

Immunological features of bats: resistance and tolerance to emerging viruses.

Bats are among the most diverse mammalian species, representing over 20% of mammalian diversity. The past two decades have witnessed a disproportionate spillover of viruses from bats to humans compared with other mammalian hosts, attributed to the viral richness within bats, their phylogenetic likeness to humans, and increased human contact with wildlife. Unique evolutionary adaptations in bat genomes, particularly in antiviral protection and immune tolerance genes, enable bats to serve as reservoirs for pandemic-inducing viruses. Here, we discuss current limitations and advances made in understanding the role of bats as drivers of pandemic zoonoses. We also discuss novel technologies that have revealed spatial, dynamic, and physiological factors driving virus and host coevolution.

Diversity, evolution, and emergence of fish viruses.

The production of aquatic animals has more than doubled over the last 50 years and is anticipated to continually increase. While fish are recognized as a valuable and sustainable source of nutrition, particularly in the context of human population growth and climate change, the rapid expansion of aquaculture coincides with the emergence of highly pathogenic viruses that often spread globally through aquacultural practices. Here, we provide an overview of the fish virome and its relevance for disease emergence, with a focus on the insights gained through metagenomic sequencing, noting potential areas for future study. In particular, we describe the diversity and evolution of fish viruses, for which the majority have no known disease associations, and demonstrate how viruses emerge in fish populations, most notably at an expanding domestic-wild interface. We also show how wild fish are a powerful and tractable model system to study virus ecology and evolution more broadly and can be used to identify the major factors that shape vertebrate viromes. Central to this is a process of virus-host co-divergence that proceeds over many millions of years, combined with ongoing cross-species virus transmission.

Wastewater-based epidemiology applied at the building-level reveals distinct virome profiles based on the age of the contributing individuals.

BACKGROUND: Human viruses released into the environment can be detected and characterized in wastewater. The study of wastewater virome offers a consolidated perspective on the circulation of viruses within a population. Because the occurrence and severity of viral infections can vary across a person's lifetime, studying the virome in wastewater samples contributed by various demographic segments can provide valuable insights into the prevalence of viral infections within these segments. In our study, targeted enrichment sequencing was employed to characterize the human virome in wastewater at a building-level scale. This was accomplished through passive sampling of wastewater in schools, university settings, and nursing homes in two cities in Catalonia. Additionally, sewage from a large urban wastewater treatment plant was analysed to serve as a reference for examining the collective excreted human virome. RESULTS: The virome obtained from influent wastewater treatment plant samples showcased the combined viral presence from individuals of varying ages, with astroviruses and human bocaviruses being the most prevalent, followed by human adenoviruses, polyomaviruses, and papillomaviruses. Significant variations in the viral profiles were observed among the different types of buildings studied. Mamastrovirus 1 was predominant in school samples, salivirus and human polyomaviruses JC and BK in the university settings while nursing homes showed a more balanced distribution of viral families presenting papillomavirus and picornaviruses and, interestingly, some viruses linked to immunosuppression. CONCLUSIONS: This study shows the utility of building-level wastewater-based epidemiology as an effective tool for monitoring the presence of viruses circulating within specific age groups. It provides valuable insights for public health monitoring and epidemiological studies.

The virota and its transkingdom interactions in the healthy infant gut.

SignificanceMicrobes colonizing the infant gut during the first year(s) of life play an important role in immune system development. We show that after birth the (nearly) sterile gut is rapidly colonized by bacteria and their viruses (phages), which often show a strong cooccurrence. Most viruses infecting the infant do not cause clinical signs and their numbers strongly increase after day-care entrance. The infant diet is clearly reflected by identification of plant-infecting viruses, whereas fungi and parasites are not part of a stable gut microbiota. These temporal high-resolution baseline data about the gut colonization process will be valuable for further investigations of pathogenic viruses, dynamics between phages and their bacterial host, as well as studies investigating infants with a disturbed microbiota.

Spatial turnover of soil viral populations and genotypes overlain by cohesive responses to moisture in grasslands.

Viruses shape microbial communities, food web dynamics, and carbon and nutrient cycling in diverse ecosystems. However, little is known about the patterns and drivers of viral community composition, particularly in soil, precluding a predictive understanding of viral impacts on terrestrial habitats. To investigate soil viral community assembly processes, here we analyzed 43 soil viromes from a rainfall manipulation experiment in a Mediterranean grassland in California. We identified 5,315 viral populations (viral operational taxonomic units [vOTUs] with a representative sequence ≥10 kbp) and found that viral community composition exhibited a highly significant distance-decay relationship within the 200-m2 field site. This pattern was recapitulated by the intrapopulation microheterogeneity trends of prevalent vOTUs (detected in ≥90% of the viromes), which tended to exhibit negative correlations between spatial distance and the genomic similarity of their predominant allelic variants. Although significant spatial structuring was also observed in the bacterial and archaeal communities, the signal was dampened relative to the viromes, suggesting differences in local assembly drivers for viruses and prokaryotes and/or differences in the temporal scales captured by viromes and total DNA. Despite the overwhelming spatial signal, evidence for environmental filtering was revealed in a protein-sharing network analysis, wherein a group of related vOTUs predicted to infect actinobacteria was shown to be significantly enriched in low-moisture samples distributed throughout the field. Overall, our results indicate a highly diverse, dynamic, active, and spatially structured soil virosphere capable of rapid responses to changing environmental conditions.

Reactivation of a transplant recipient's inherited human herpesvirus 6 and implications to the graft.

BACKGROUND: The implications of inherited chromosomally integrated human herpesvirus 6 (iciHHV-6) in solid organ transplantation remain uncertain. Although this trait has been linked to unfavorable clinical outcomes, an association between viral reactivation and complications has only been conclusively established in a few cases. In contrast to these studies, which followed donor-derived transmission, our investigation is the first to examine the pathogenicity of a recipient´s iciHHV-6B and its impact on the graft. METHODS: We used hybrid capture sequencing for in-depth analysis of the viral sequences reconstructed from sequential liver biopsies. Moreover, we investigated viral replication through in situ hybridization (U38-U94 genes), real-time PCR (U89/U90 genes), immunohistochemistry, and immunofluorescence (against viral lysate). We also performed whole transcriptome sequencing of the liver biopsies to profile the host immune response. RESULTS: We report a case of reactivation of a recipient´s iciHHV-6B and subsequent infection of the graft. Using a novel approach integrating the analysis of viral and mitochondrial DNAs, we located the iciHHV-6B intra-graft. We demonstrated active replication via the emergence of viral minor variants across time points, in addition to positive viral mRNAs and antigen stainings in tissue sections. Furthermore, we detected significant upregulation of cell surface molecules, transcription factors, and cytokines associated with antiviral immune responses, arguing against immunotolerance. CONCLUSIONS: Our analysis underscores the potential pathological impact of iciHHV-6B, emphasizing the need for close monitoring of reactivation in transplant recipients. Most crucially, it highlights the critical role that the host's virome can play in shaping the outcome of transplantation, urging further investigations.

Entourage: all-in-one sequence analysis software for genome assembly, virus detection, virus discovery, and intrasample variation profiling.

BACKGROUND: Pan-virus detection, and virome investigation in general, can be challenging, mainly due to the lack of universally conserved genetic elements in viruses. Metagenomic next-generation sequencing can offer a promising solution to this problem by providing an unbiased overview of the microbial community, enabling detection of any viruses without prior target selection. However, a major challenge in utilising metagenomic next-generation sequencing for virome investigation is that data analysis can be highly complex, involving numerous data processing steps. RESULTS: Here, we present Entourage to address this challenge. Entourage enables short-read sequence assembly, viral sequence search with or without reference virus targets using contig-based approaches, and intrasample sequence variation quantification. Several workflows are implemented in Entourage to facilitate end-to-end virus sequence detection analysis through a single command line, from read cleaning, sequence assembly, to virus sequence searching. The results generated are comprehensive, allowing for thorough quality control, reliability assessment, and interpretation. We illustrate Entourage's utility as a streamlined workflow for virus detection by employing it to comprehensively search for target virus sequences and beyond in raw sequence read data generated from HeLa cell culture samples spiked with viruses. Furthermore, we showcase its flexibility and performance on a real-world dataset by analysing a preassembled Tara Oceans dataset. Overall, our results show that Entourage performs well even with low virus sequencing depth in single digits, and it can be used to discover novel viruses effectively. Additionally, by using sequence data generated from a patient with chronic SARS-CoV-2 infection, we demonstrate Entourage's capability to quantify virus intrasample genetic variations, and generate publication-quality figures illustrating the results. CONCLUSIONS: Entourage is an all-in-one, versatile, and streamlined bioinformatics software for virome investigation, developed with a focus on ease of use. Entourage is available at https://codeberg.org/CENMIG/Entourage under the MIT license.

A metagenomic investigation of the faecal RNA virome structure of asymptomatic chickens obtained from a commercial farm in Durban, KwaZulu-Natal province, South Africa.

BACKGROUND: Virome studies on birds, including chickens are relatively scarce, particularly from the African continent. Despite the continuous evolution of RNA viruses and severe losses recorded in poultry from seasonal viral outbreaks, the information on RNA virome composition is even scantier as a result of their highly unstable nature, genetic diversity, and difficulties associated with characterization. Also, information on factors that may modulate the occurrence of some viruses in birds is limited, particularly for domesticated birds. Viral metagenomics through advancements in sequencing technologies, has enabled the characterization of the entire virome of diverse host species using various samples. METHODS: The complex RNA viral constituents present in 27 faecal samples of asymptomatic chickens from a South African farm collected at 3-time points from two independent seasons were determined, and the impact of the chicken's age and collection season on viral abundance and diversity was further investigated. The study utilized the non-invasive faecal sampling method, mRNA viral targeted enrichment steps, a whole transcriptome amplification strategy, Illumina sequencing, and bioinformatics tools. RESULTS: The results obtained revealed a total of 48 viral species spanning across 11 orders, 15 families and 21 genera. Viral RNA families such as Coronaviridae, Picornaviridae, Reoviridae, Astroviridae, Caliciviridae, Picorbirnaviridae and Retroviridae were abundant, among which picornaviruses, demonstrated a 100% prevalence across the three age groups (2, 4 and 7 weeks) and two seasons (summer and winter) of the 27 faecal samples investigated. A further probe into the extent of variation between the different chicken groups investigated indicated that viral diversity and abundance were significantly influenced by age (P = 0.01099) and season (P = 0.00099) between chicken groups, while there was no effect on viral shedding within samples in a group (alpha diversity) for age (P = 0.146) and season (P = 0.242). CONCLUSION: The presence of an exceedingly varied chicken RNA virome, encompassing avian, mammalian, fungal, and dietary-associated viruses, underscores the complexities inherent in comprehending the causation, dynamics, and interspecies transmission of RNA viruses within the investigated chicken population. Hence, chickens, even in the absence of discernible symptoms, can harbour viruses that may exhibit opportunistic, commensal, or pathogenic characteristics.

Benchmarking of virome metagenomic analysis approaches using a large, 60+ members, viral synthetic community.

IMPORTANCE: We report here efforts to benchmark performance of two widespread approaches for virome analysis, which target either virion-associated nucleic acids (VANA) or highly purified double-stranded RNAs (dsRNAs). This was achieved using synthetic communities of varying complexity levels, up to a highly complex community of 72 viral agents (115 viral molecules) comprising isolates from 21 families and 61 genera of plant viruses. The results obtained confirm that the dsRNA-based approach provides a more complete representation of the RNA virome, in particular, for high complexity ones. However, for viromes of low to medium complexity, VANA appears a reasonable alternative and would be the preferred choice if analysis of DNA viruses is of importance. Several parameters impacting performance were identified as well as a direct relationship between the completeness of virome description and sample sequencing depth. The strategy, results, and tools used here should prove useful in a range of virome analysis efforts.

Distinct characteristics of the gut virome in patients with osteoarthritis and gouty arthritis.

BACKGROUND/PURPOSE(S): The gut microbiota and its metabolites play crucial roles in pathogenesis of arthritis, highlighting gut microbiota as a promising avenue for modulating autoimmunity. However, the characterization of the gut virome in arthritis patients, including osteoarthritis (OA) and gouty arthritis (GA), requires further investigation. METHODS: We employed virus-like particle (VLP)-based metagenomic sequencing to analyze gut viral community in 20 OA patients, 26 GA patients, and 31 healthy controls, encompassing a total of 77 fecal samples. RESULTS: Our analysis generated 6819 vOTUs, with a considerable proportion of viral genomes differing from existing catalogs. The gut virome in OA and GA patients differed significantly from healthy controls, showing variations in diversity and viral family abundances. We identified 157 OA-associated and 94 GA-associated vOTUs, achieving high accuracy in patient-control discrimination with random forest models. OA-associated viruses were predicted to infect pro-inflammatory bacteria or bacteria associated with immunoglobulin A production, while GA-associated viruses were linked to Bacteroidaceae or Lachnospiraceae phages. Furthermore, several viral functional orthologs displayed significant differences in frequency between OA-enriched and GA-enriched vOTUs, suggesting potential functional roles of these viruses. Additionally, we trained classification models based on gut viral signatures to effectively discriminate OA or GA patients from healthy controls, yielding AUC values up to 0.97, indicating the clinical utility of the gut virome in diagnosing OA or GA. CONCLUSION: Our study highlights distinctive alterations in viral diversity and taxonomy within gut virome of OA and GA patients, offering insights into arthritis etiology and potential treatment and prevention strategies.

Integrative metagenomic analysis reveals distinct gut microbial signatures related to obesity.

Obesity is a metabolic disorder closely associated with profound alterations in gut microbial composition. However, the dynamics of species composition and functional changes in the gut microbiome in obesity remain to be comprehensively investigated. In this study, we conducted a meta-analysis of metagenomic sequencing data from both obese and non-obese individuals across multiple cohorts, totaling 1351 fecal metagenomes. Our results demonstrate a significant decrease in both the richness and diversity of the gut bacteriome and virome in obese patients. We identified 38 bacterial species including Eubacterium sp. CAG:274, Ruminococcus gnavus, Eubacterium eligens and Akkermansia muciniphila, and 1 archaeal species, Methanobrevibacter smithii, that were significantly altered in obesity. Additionally, we observed altered abundance of five viral families: Mesyanzhinovviridae, Chaseviridae, Salasmaviridae, Drexlerviridae, and Casjensviridae. Functional analysis of the gut microbiome indicated distinct signatures associated to obesity and identified Ruminococcus gnavus as the primary driver for function enrichment in obesity, and Methanobrevibacter smithii, Akkermansia muciniphila, Ruminococcus bicirculans, and Eubacterium siraeum as functional drivers in the healthy control group. Additionally, our results suggest that antibiotic resistance genes and bacterial virulence factors may influence the development of obesity. Finally, we demonstrated that gut vOTUs achieved a diagnostic accuracy with an optimal area under the curve of 0.766 for distinguishing obesity from healthy controls. Our findings offer comprehensive and generalizable insights into the gut bacteriome and virome features associated with obesity, with the potential to guide the development of microbiome-based diagnostics.

Metagenomic analysis of the gut virome in patients with irritable bowel syndrome.

Irritable bowel syndrome (IBS), a chronic functional gastrointestinal disorder, is recognized for its association with alterations in the gut microbiome and metabolome. This study delves into the largely unexplored domain of the gut virome in IBS patients. We conducted a comprehensive analysis of the fecal metagenomic data set from 277 IBS patients and 84 healthy controls to characterize the gut viral community. Our findings revealed a distinct gut virome in IBS patients compared to healthy individuals, marked by significant variances in between-sample diversity and altered abundances of 127 viral operational taxonomic units (vOTUs). Specifically, 111 vOTUs, predominantly belonging to crAss-like, Siphoviridae, Myoviridae, and Quimbyviridae families, were more abundant in IBS patients, whereas the healthy control group exhibited enrichment of 16 vOTUs from multiple families. We also investigated the interplay between the gut virome and bacteriome, identifying a correlation between IBS-enriched bacteria like Klebsiella pneumoniae, Fusobacterium varium, and Ruminococcus gnavus, and the IBS-associated vOTUs. Furthermore, we assessed the potential of gut viral signatures in predicting IBS, achieving a notable area under the receiver operator characteristic curve (AUC) of 0.834. These findings highlight significant shifts in the viral diversity, taxonomic distribution, and functional composition of the gut virome in IBS patients, suggesting the potential role of the gut virome in IBS pathogenesis and opening new avenues for diagnostic and therapeutic strategies targeting the gut virome in IBS management.

Metagenomic analysis reveals altered gut virome and diagnostic potential in pancreatic cancer.

Pancreatic cancer (PC) is a highly aggressive malignancy with a poor prognosis, making early diagnosis crucial for improving patient outcomes. While the gut microbiome, including bacteria and viruses, is believed to be essential in cancer pathogenicity, the potential contribution of the gut virome to PC remains largely unexplored. In this study, we conducted a comparative analysis of the gut viral compositional and functional profiles between PC patients and healthy controls, based on fecal metagenomes from two publicly available data sets comprising a total of 101 patients and 82 healthy controls. Our results revealed a decreasing trend in the gut virome diversity of PC patients with disease severity. We identified significant alterations in the overall viral structure of PC patients, with a meta-analysis revealing 219 viral operational taxonomic units (vOTUs) showing significant differences in relative abundance between patients and healthy controls. Among these, 65 vOTUs were enriched in PC patients, and 154 were reduced. Host prediction revealed that PC-enriched vOTUs preferentially infected bacterial members of Veillonellaceae, Enterobacteriaceae, Fusobacteriaceae, and Streptococcaceae, while PC-reduced vOTUs were more likely to infect Ruminococcaceae, Lachnospiraceae, Clostridiaceae, Oscillospiraceae, and Peptostreptococcaceae. Furthermore, we constructed random forest models based on the PC-associated vOTUs, achieving an optimal average area under the curve (AUC) of up to 0.879 for distinguishing patients from controls. Through additional 10 public cohorts, we demonstrated the reproducibility and high specificity of these viral signatures. Our study suggests that the gut virome may play a role in PC development and could serve as a promising target for PC diagnosis and therapeutic intervention. Future studies should further explore the underlying mechanisms of gut virus-bacteria interactions and validate the diagnostic models in larger and more diverse populations.

Identification and genetic characterization of a recently identified enterovirus C116 in China.

Enterovirus C116 (EV-C116) is a new member of the enterovirus C group which is closely associated with several infectious diseases. Although sporadic studies have detected EV-C116 in clinical samples worldwide, there is currently limited information available. In this study, two EV-C-positive fecal specimens were detected in apparently healthy children, which harbored low abundance, through meta-transcriptome sequencing. Based on the prototypes of several EV-Cs, two lineages were observed. Lineage 1 included many types that could not cause EV-like cytopathic effect in cell culture. Three genogroups of EV-C116 were divided in the maximum likelihood tree, and the two strains in this study (XZ2 and XZ113) formed two different lineages, suggesting that EV-C116 still diffuses worldwide. Obvious inter-type recombination events were observed in the XZ2 strain, with CVA22 identified as a minor donor. However, another strain (XZ113) underwent different recombination situations, highlighting the importance of recombination in the formation of EV-Cs biodiversity. The EV-C116 strains could propagate in rhabdomyosarcoma cell cultures at low titer; however, EV-like cytopathic effects were not observed. HEp-2, L20B, VERO, and 293T cell lines did not provide an appropriate environment for EV-C116 growth. These results challenge the traditional recognition of the uncultured nature of EV-C116 strains and explain the difficulty of clinical detection.

Multikingdom characterization of gut microbiota in patients with rheumatoid arthritis and rheumatoid arthritis-associated interstitial lung disease.

Rheumatoid arthritis-associated interstitial lung disease (RA-ILD) is a serious and common extra-articular disease manifestation. Patients with RA-ILD experience reduced bacterial diversity and gut bacteriome alterations. However, the gut mycobiome and virome in these patients have been largely neglected. In this study, we performed whole-metagenome shotgun sequencing on fecal samples from 30 patients with RA-ILD, and 30 with RA-non-ILD, and 40 matched healthy controls. The gut bacteriome and mycobiome were explored using a reference-based approach, while the gut virome was profiled based on a nonredundant viral operational taxonomic unit (vOTU) catalog. The results revealed significant alterations in the gut microbiomes of both RA-ILD and RA-non-ILD groups compared with healthy controls. These alterations encompassed changes in the relative abundances of 351 bacterial species, 65 fungal species, and 4,367 vOTUs. Bacteria such as Bifidobacterium longum, Dorea formicigenerans, and Collinsella aerofaciens were enriched in both patient groups. Ruminococcus gnavus (RA-ILD), Gemmiger formicilis, and Ruminococcus bromii (RA-non-ILD) were uniquely enriched. Conversely, Faecalibacterium prausnitzii, Bacteroides spp., and Roseburia inulinivorans showed depletion in both patient groups. Mycobiome analysis revealed depletion of certain fungi, including Saccharomyces cerevisiae and Candida albicans, in patients with RA compared with healthy subjects. Notably, gut virome alterations were characterized by an increase in Siphoviridae and a decrease in Myoviridae, Microviridae, and Autographiviridae in both patient groups. Hence, multikingdom gut microbial signatures showed promise as diagnostic indicators for both RA-ILD and RA-non-ILD. Overall, this study provides comprehensive insights into the fecal virome, bacteriome, and mycobiome landscapes of RA-ILD and RA-non-ILD gut microbiota, thereby offering potential biomarkers for further mechanistic and clinical research.