Abstract
Introduction:
Autism spectrum disorder (ASD) is a complex neurodevelopmental disorder, which brings a great burden to the family and society. Gut microbiota is considered to be an important factor in ASD that easily affects function and development of the immune, metabolic, and nervous systems. However, most available studies have mainly focused on the altered gut bacteria, our knowledge of gut viruses in ASD children remains limited.
Methods:
In this study, we collected fecal samples from ASD children and healthy controls, then analyzed and compared the differences of the gut viral communities between the two groups by viral metagenomic techniques.
Results:
The alpha diversity of the ASD virome was lower than that of the healthy virome, and the beta diversity had a significant difference between ASD and healthy children. Podoviridae accounted for the highest proportion of viruses in ASD patients, while Alphaflexiviridae was dominant in healthy controls. There was a statistical difference in the abundance of Microviridae between the two groups. Additionally, human astrovirus, picobirnavirus, and norovirus were detected by phylogenetic analysis.
Discussion:
This study revealed that alpha diversity was reduced in children with ASD, and different compositions in gut viral communities were observed between ASD patients and healthy controls. Changes in viral diversity and composition deepen our understanding of the differences in the gut viral communities between ASD and healthy children, and also provides a perspective for further exploration of viruses related to ASD children.
1 Introduction
Autism spectrum disorder (ASD) is a complex neurodevelopmental disorder, which is characterized by impairments in social communication and interaction, restricted patterns of interest, as well as repetitive and stereotyped behaviours (). The incidence and prevalence of ASD have increased dramatically worldwide over the past decade (Zeidan et al., 2022). Today, ASD is one of the most prevalent neurodevelopmental disorders that affect children (). It is expected that the number of children with ASD will exceed the number of children with cancer, juvenile diabetes, and pediatric acquired immune deficiency syndrome (AIDS) combined, which can result in serious medical and social problems (). Frustratingly, little is understood about the etiology of ASD, but it is known that one of the main pathophysiology features of ASD is chronic inflammation in intestinal mucosa (), and increasing evidences indicate that gut microbiota that impacts development and function of the immune, metabolic, and nervous systems may play an important role in ASD (; ; Zhang et al., 2020; Willyard, 2021).
The gut microbiota includes bacteria, archaea, viruses, fungi, and helminths (Vemuri et al., 2020). Collectively, their respective genomes are called the gut microbiome (Vemuri et al., 2018). Bacteria have gained the most attention in these microorganisms. There is considerable evidence that gut bacteria is closely associated with ASD (; ). Patients with ASD were shown to have decreased abundance of Bifidobacterium, Blautia, Dialister, Prevotella, Veillonella, and Turicibacter (). Futhermore, decreased Bacteroidetes/Firmicutes ratio, increased Lactobacillus and Desulfovibrio species were also present in ASD patients (Tomova et al., 2015). However, merging findings have suggested that as part of the gut microecology, gut virome also appears associated with ASD (; Wan et al., 2024b).
Viruses are abundant in the gut. It is estimated that each gram of human intestinal content includes at least 108–1010 virus-like particles(VLPs) (). In addition to affecting the digestive system, enteroviruses can also impact the neurological system. At present, several viruses, such as those from the families Astroviridae, Picobirnaviridae, and Caliciviridae, have been reported in some cases of neurological damage (; ; ; ). Some enteroviruses can directly infect neurons or glial cells, causing cell damage and death through viral replication and lysis infection, thereby disrupting neural function (). Furthermore, viral infection can also activate the immune system, leading to an increase in the level of pro-inflammatory cytokines and inducing a systemic inflammatory response, which could affect brain development and thereby cause neurodevelopmental deficits (). A recent report found that there are significant differences in gut-DNA virome abundance and heterogeneity between ASD and typically developing (TD) children (Wan et al., 2024b). In addition, it has been reported that alpha diversity was significantly lower in ASD children especially those who have gastrointestinal symptoms and elevated Skunavirus abundance was also observed in children with ASD compared to healthy controls (). However, these studies mainly focus on DNA viruses and do not cover RNA viruses or include analysis of the detected viruses to determine their relationship with previously identified viruses. Overall, in the existing microbiota studies of ASD, the exploration of gut viruses is relatively insufficient.
In recent years, with the advancement of next-generation sequencing technologies, viral metagenomics has been increasingly applied in the study of microbial communities. Different from traditional methods, this technology can rapidly identify all viral sequences in a sample, providing a more comprehensive acquisition of information about the viral community within the sample. It not only expands the information reserve of known viruses but also reveals the characteristics of viruses (Wen et al., 2024).
Therefore, in this study, we further investigated the composition and differences of the gut viral community between ASD patients and healthy controls using viral metagenomic techniques, which provided new insights for us to better understand the impact of gut microecological environment on ASD.
2 Materials and methods
2.1 Sample collection and preparation
To analyze the gut viral communities between ASD and healthy children, we collected 11 fecal samples from ASD children admitted to the Pediatric Department, Affiliated Hospital of Jiangsu University (Jiangsu Province, China) and 11 fecal samples from healthy children who were in the same hospital for physical examination from September 2023 to March 2024. About 10g of feces were taken from each sample. All ASD patients met the diagnosis standards stipulated in the American Diagnostic and Statistical Manual of Mental Disorders, 5th edition (). Meantime, all participants were required to be no use of antibiotics or antiviral drugs within 3 weeks, no infectious diseases or psychiatric disorders. Subjects with complaints of gastrointestinal symptoms through parent interviews were also excluded. All samples were stored in sterile covered containers and transported to the laboratory via dry ice. In advance of viral metagenomic analysis, each sample was resuspended in 10 volumes of sterile-filtered Dulbecco’s phosphate-buffered saline (DPBS) and vortexed for 5 min vigorously. Then after centrifugation (10 min, 15,000 g, 4°C), the supernatants were collected and stored at -80°C (Wang et al., 2019). Sample collection and all experiments in this study obtained approval from the Ethics Committee of Jiangsu University and the reference number is No.KY2023K0805.
2.2 Viral metagenomic analysis
About five hundred microliters of supernatant was drawn from each sample and pooled into the separate sample pool. Sample pools were centrifuged (5 min, 12,000 g, 4°C), and then the supernatant was filtered through a 0.45 μm filter (Millipore) to remove eukaryotic and bacterial cell-sized particles (). The filtrate rich in viral particles was treated with DNase and RNase to break down unprotected nucleic acid at 37°C for 60min (Zhang et al., 2014, 2016, 2017). The treatment efficiency of DNase and RNase was validated by quantifying the abundance of the nucleic acid gene via qPCR before and after digestion, which showed a significant reduction. According to the manufacturer’s protocol, total viral nucleic acids were extracted using the QIAamp Viral RNA Mini Kit (Qiagen). Reverse transcription reactions were carried out on these nucleic acid samples composed of DNA and RNA viral sequences with SuperScript III reverse transcriptase and random hexamer primers, followed by a single round of DNA synthesis using Klenow fragment polymerase. All dsDNA products were used for library construction. Then, 22 libraries were constructed using Nextera XT DNA Sample Preparation Kit (Illumina). For bioinformatics analysis, 150-bp paired-end reads generated by NovaSeq were debarcoded using vendor software provided by Illumina (). Clonal reads were removed and low-quality sequencing tails were trimmed using Phred quality score 20 as the threshold. The cleaned reads were then de novo assembled by SOAPdenovo2 version r240 using Kmer size 63 with default settings. The assembled contigs and singlets were matched against an in-house viral proteome database using BLASTx with a cut-off E-value of <10-5, where the virus BLASTx database was compiled using NCBI virus reference proteome (ftp://ftp.ncbi.nih.gov/refseq/release/viral/) and viral protein sequences from NCBI nr FASTA files (based on annotation taxonomy in Virus Kingdom) (). The selection of the BLASTx E-value cutoff of <10-5 could ensured annotation accuracy while effectively controlled the false discovery rate (). Candidate viral sequences were further compared to the non-virus non-redundant protein database to eliminate false-positive viral sequences (). Contigs without significant BLASTx similarity was compared to viral protein families in vFam database using HMMER3 to identify remote viral protein similarities (; ; ). These BLASTx results generated by DIAMOND (DAA format) were imported into Megan6 software for generation of rma6 format files, which could be further used for subsequent analysis including species accumulation curve, and species rarefaction curve.
2.3 Statistical analysis
Utilizing Megan v6.25.7 and R v4.3.2, the statistical analysis related to this study was performed. The composition analysis from 22 libraries was standardized and compared by Megan (). The viral community structure and richness were visualized using the pheatmap and vegan package. Besides, the differences in viral communities were demonstrated using the ggplot2 package and Statistical Analysis of Metagenomic Profiles (STAMP) between the ASD patients group and the healthy controls group. For the multiple-testing involved in the analysis, Benjamini-Hochberg False Discovery Rate (FDR) correction method was adopted. When P <Â 0.05, it was considered to be statistically significant.
2.4 Phylogenetic analysis
Phylogenetic analysis was based on predicted nucleotide and protein sequences in this study along with their closest viral relatives that were downloaded from the NCBI GenBank database. Depending on the distance of the evolutionary relationship, for viruses with close kinship, we chose nucleic acid sequences to construct phylogenetic trees, while amino acid sequences were selected for viruses with distant genetic relationship. Related nucleotide and protein sequences were aligned using MUSCLE in MEGA v11.0.13 with default settings (). Sites which contained more than 50% gaps were temporarily removed from alignments. Bayesian inference trees were then built using MrBayes v.3.2.7 (). Moreover, maximum likelihood trees were constructed to confirm all the Bayesian inference trees by MEGA software with 1000 bootstrap replicates (). Adobe Illustrator and iTOL (https://iTOL.embl.de/) were used to represent the phylogenetic trees.
2.5 Quality control
In order to prevent cross-sample contamination and nucleic acid degradation, a robust viral metagenomics pipeline integrated with multiple contamination controls was followed throughout the process. Specifically, sterile ddH2O (Sangon Biotech) was prepared and further processed under the same conditions as a blank control to rule out the possibility of nucleic acid contamination. Additionally, all materials that came into direct contact with nucleic acid samples were free of DNase and RNase. DEPC-treated water and RNase inhibitors were used to dissolve nucleic acid samples to preserve nucleic acid integrity.
3 Results
3.1 Baseline characteristics of ASD and healthy children
The study participants included 11 ASD children and 11 healthy children. The mean age was 5.2 and 4.6 years, respectively. There was no significant statistical difference in age, gender or mode of birth between the two groups (P > 0.05) (Supplementary Table 1).
3.2 Overview of the virome
In the study, after library construction and next generation sequencing on the Illumina NovaSeq platform, the 22 libraries yielded a total of 52,261,771 raw reads with an average percentage of GC content (GC%) of 44.8%. 37,554,145 reads were related to virus in these libraries according to the de novo assembled and compared to the GenBank non-redundant protein database using BLASTx (Supplementary Table 2). The species richness of all collected samples was revealed by the species rarefaction and accumulation curves (Supplementary Figure 1). In majority of 22 libraries, the observed virus species tended to level out at the end, which indicates that the sequencing depth was sufficient to cover all species present in the samples. Even if the sequencing data is increased, the diversity of virus species cannot increase. And the accumulation curves revealed that there may be more than 800 different viruses in 22 libraries.
3.3 Diversity and richness of gut virome in ASD and healthy children
To investigate the differences in diversity and richness between the ASD group and the healthy group, alpha diversity analysis and beta diversity analysis were performed on the viral communities from the two groups. In alpha diversity analysis, we could see that α diversity was reduced in children with ASD. It was significantly lower than that in the healthy controls group at the phylum, class, and order levels, but not at the family or genus levels (Figure 1). In terms of beta diversity analysis, principal coordinate analysis (PCoA) based on unweighted UniFrac metrics revealed that the composition of the viral communities exhibited statistically significant differences across the two groups at the phylum, class, order, family, and genus levels (Figure 2).
Figure 1
Figure 2
3.4 Composition and comparison in the viral communities
Heatmap at the family level revealed that 51 viral families consisting of 29 double-stranded DNA (dsDNA) viral families, 8 single-stranded DNA (ssDNA) viral families, 2 dsRNA viral families, and 12 ssRNA viral families were present in the 22 libraries (Figure 3A). Among them, a total of 18 families showed different compositions (≥2-fold change) between ASD patients and healthy controls (Supplementary Table 3). In ASD patients, Podoviridae accounted for the highest proportion of viruses, reaching 58.81%, followed by Caliciviridae (31.06%), Virgaviridae (5.33%), and Siphoviridae (2.51%). Each of the family Microviridae, Myoviridae, and Salasmaviridae accounted for more than 1% of viral community compositions. On the other hand, Alphaflexiviridae was dominant in the healthy controls group with 50.82%, which was a drastic increase from 0.08% in ASD patients. The next highest proportions of viruses were Podoviridae (31.64%) and Microviridae (14.46%). For the top ten most abundant viral families between the two groups, we noticed that variations in their proportion in different libraries (Figure 3B). To further investigate the differences in the composition of gut viral communities between the ASD group and the healthy group, Statistical Analysis of Metagenomic Profiles analysis was performed on the viral communities from the two groups. According to the results of STAMP analysis, Phixviricota, Malgrandaviricetes, Petitvirales, and Microviridae made the greatest contribution to the differences in the healthy control group and the ASD patients group (P < 0.01), followed by Cressdnaviricota, Repensiviricetes, and Geplafuvirales (P < 0.05) (Figure 3C).
Figure 3
3.5 Astroviridae
Astroviruses are non-enveloped and single-stranded positive RNA viruses of the Astroviridae family with a size of about 6.8-7.9 kb that infect a variety of mammals and avian species (; ). The viral genome contains three overlapping open reading frames (ORFs): ORF1a, ORF1b and ORF2. In astroviruses, ORF1a and ORF1b encode the viral protease and RNA-dependent RNA polymerase (RdRp) respectively, while ORF2 encodes the capsid protein (). In this study, we obtained a nearly complete astrovirus genome from the library Humanfe016. According to the result of the search in GenBank, it was 99.73% identical to Human astrovirus 1 (GenBank no.PQ510859) collected in Yunnan, China, in 2024. Then, we performed a phylogenetic tree based on the RdRp protein, which indicated AstroviridaeHumanfe016 clustered inside the HAstV-1 genus (Figure 4).
Figure 4
3.6 Picobirnaviridae
Picobirnaviruses, the only genus within the family Picobirnaviridae, are non-enveloped icosahedral viruses with a diameter of 33–37 nm consisting of two segments (). The large segment encodes capsid protein while the small segment encodes RdRp (Xiao et al., 2020). Picobirnaviruses can infect humans and a wide range of animals (). We identified two strains belonging to Picobirnaviridae in this study. In order to analyze the relationship between the newly discovered genome and other known picobirnaviruses, the respective phylogenetic analysis trees were constructed based on the capsid protein (Figure 5A) and RdRp (Figure 5B), respectively. The results indicated that our sequences clustered with others throughout the phylogenetic trees.
Figure 5
3.7 Caliciviridae
Caliciviruses are non-enveloped icosahedral viruses with a single-stranded, positive-sense, polyadenylated RNA genome of 7.4-8.3 kb arranged in two or three major ORFs (Vinjé et al., 2019). Caliciviruses can cause a wide range of diseases in humans and animals, including gastroenteritis, neurological disorders, respiratory infections, and vesicular lesions (). The family Caliciviridae comprises eleven virus genera, mainly norovirus, lagovirus, nebovirus, sapovirus, and vesivirus (). The most important clinical representative is noroviruses.
We identified and assembled 6,971,747 reads obtained from sequencing, obtained the genomes of 10 nearly complete caliciviruses through alignment with the NCBI GenBank database using BLASTx data, and then constructed a phylogenetic tree based on the RdRp protein. The phylogenetic tree indicated that caliciviruses in this study had high diversity and were clustered into groups (Figure 6). By sequence analysis and phylogeny, nine of these caliciviruses were closely related to genus Norovirus GI.3 which shared more than 99% similarities of their nucleic acid sequences, while Humanfe021C2 and Norovirus GI.7 (Genbank No.LC646334) shared a branch.
Figure 6
4 Discussion
In recent years, the worldwide prevalence of ASD has been on the rise (Zeidan et al., 2022), which may be due to the increase in people’s awareness of the disease, the changes in diagnostic criteria, and the increase in the availability of diagnosis (). The exact etiopathogenesis of ASD is poorly understood and appear to involve a complex interplay of genetic and environmental factors, of which the gut microbiome is an important environmental factor related to a broad range of neurological and psychiatric disorders and diseases (; ; ; Vemuri et al., 2020). It is reported that the microbiota could shape ASD development and progression by impacting the host defense and inflammatory response (). However, most of the available research literature has focused on describing the altered gut bacteria, there are few studies on the gut virome in the feces of ASD.
In order to compare the gut viral communities between ASD children and healthy controls, we collected fecal samples from children with ASD and healthy controls, then used viral metagenomic analysis to investigate the composition and differences of the gut virome between the two groups in this study. The results showed that there were no differences between ASD and healthy children in terms of age, gender or mode of birth. By depicting the species rarefaction curve and species accumulation curve, the rationality of the sequencing data was illustrated through the species rarefaction curve and species accumulation curve. It should be acknowledged that the inclusion of over 800 different viral species contained in the 22 libraries appears a bit high. This may be because the reads belonging to the same viral genome have different best matches in BLASTx (E-value cutoff of <10−5), leading to an overestimation of the species of viruses. To mitigate this, CD-HIT could be used to cluster the assembled contigs and singlets, thereby avoiding redundant annotations from cross-homology.
The human gut virome is characterized by high diversity (). However, compare to healthy children, studies have shown a decrease in the diversity of ASD subjects (Tong et al., 2022; ). In our study, alpha diversity was reduced in ASD patients, which was consistent with former studies, and it was significantly lower than that of the healthy controls group in phylum to order level. It is reported that a decrease in microbial diversity is considered to represent in states of disease (). Although the etiology of ASD remains unclear, experimental and clinical studies suggest that ASD is associated with gut microbiota dysbiosis and dysregulation (Zhang et al., 2020; ). Reduced gut viral diversity in children with ASD may create an environment conducive to the overgrowth of pathogenic bacteria and reduced colonization of beneficial bacterial, thereby exacerbating intestinal inflammation (; Wan et al., 2024a). Previous studies have linked the loss of gut viral diversity to dysregulation of bacterial hosts and disrupted intestinal barrier function, both of which are pathways implicated in neurodevelopmental alterations associated with ASD via gut-brain axis signaling (; ). Thus, our findings, together with existing literature, indicate that reduced gut viral alpha diversity is potentially correlated with ASD susceptibility. Besides, the composition of the viral communities exhibited statistically significant differences in the two groups based on beta diversity.
We also observed that there were differences in the abundance of gut virome between ASD patients and healthy controls, and Phixviricota–Malgrandaviricetes–Petitvirale–Microviridae made the greatest contribution to the differences. Recent studies indicated that alterations in the abundance of gut virome may drive gut viral dysbiosis in children with ASD (Wan et al., 2024b; Zeng et al., 2024). Microviridae is a small, single-stranded DNA virus that primarily infects bacteria. As one of the common bacteriophages in human intestine, the family Microviridae can influence microbial abundance and the composition of microbial communities. By infecting and lysing specific bacterial hosts, Microviridae plays a role in regulating the quantity and structure of bacterial populations in the intestinal microecology (). The infecting hosts of this viral family contain Bacteroides and Prevotella (). We found a decreased abundance of Microviridae in ASD children in our study. Meanwhile, recent research has shown that the amount of Bacteroides and Prevotella increased in children with ASD compared to the controls (; Zou et al., 2020). This implies that the weakened lytic regulation effect of such bacteriophages on their host bacteria may reduce the growth inhibition pressure on potential host bacteria, which creates conditions for the increased abundance of the latter. Furthermore, in this study, DNA bacteriophage Podoviridae of Caudovirales was dominated in the gut virome of ASD children and its abundance is higher compared to healthy children. It is important to note that this observed dominance of Podoviridae could be influenced by confounding factors not measured in our study, such as dietary habits or environmental exposures, which may affect the gut virome composition. Therefore, the association between Podoviridae and ASD identified here may reflect a broader ecological shift. DNA nature bacteriophages are the most predominant members of the gut virome, which account for 90% of the intestinal viruses (). Phages have been reported to shape human microbiome diversity and functions through their ability to lyse and kill host bacteria, thereby affecting brain function, mood, and memory (; ; ; Wu et al., 2022). Wan et al. found an altered gut viral community accompanied by the higher abundance of Clostridium phage, Bacillus phage,and Enterobacteria phage in children with ASD than in healthy individuals (Wan et al., 2024b). Changes in the abundance of these phages are also considered a triggering factor for gut viral ecological system dysbiosis. However, this study used shotgun sequencing, a method that often yields data mixed with a large number of non-viral impurity sequences. This makes it difficult to efficiently and accurately isolate gut viral sequences from the dataset, and the subsequent data analysis work is also extremely cumbersome. The team of Tong reported that Streptococcal phages were significantly enriched in the oral phage communities of ASD (Tong et al., 2022). However, the underlying reasons for the enrichment of Streptococcal phages in the oral communities of ASD remain unknown. Moreover, many factors, such as narrow food ranges caused by allergies and intolerances, toxic exposures, and poor oral hygiene, may influence the oral microbial environment. Overall, these studies indicate that bacteriophages are potentially associated with the development of ASD.
Moreover, in this study, phylogenetic analysis trees of possible pathogenic viruses were also constructed. Human astroviruses (HAstVs) mainly affect children, elderly people and immunocompromised individuals (). As neurotropic enteroviruses, astroviruses may lead to neuroinflammation and neuronal damage, thereby affecting brain development and function (). Notably, some studies have found that besides enteric diseases, HAstV infections are associated with neurological disorders such as acute encephalopathy (; ; ). We acquired one strain of astrovirus in the fecal sample of a patient with ASD which shared 99.73% nucleic acid identity to their best match. Based on phylogenetic analysis, the human astrovirus was identified as Human astrovirus 1 (HAstV1). Another virus we found is picobirnavirus which was considered opportunistic enteric pathogens in many species, including humans (). The phylogentic analyses showed that two strains from library Humanfe013 clustered together with other sequences, respectively.
As the major viral enteric pathogen, noroviruses also affect the neurological system in addition to affecting the digestive system (). The presence of norovirus in the central nervous system has been confirmed in the last few years (; ). Noroviruses can infect brain endothelial cells and increase the expression of matrix metalloproteinases, changing the expression of tight junction proteins and destroying the blood brain barrier (). A recent study revealed that they detected human norovirus in stool samples from a patient with seizures and white matter lesions (). All calicivirus sequences found in this study were norovirus. Moreover, the noroviruses identified in this study were detected in the ASD group and accounted for a high proportion of the virome, but not in the healthy control group. While this suggests a potential link between norovirus and ASD, this association should be interpreted cautiously. Many confounding factors, including environmental exposure, dietary patterns, genetic background, and immune status, may contribute to norovirus persistence and altered neurological outcomes. Further investigation is needed to determine whether norovirus participates in ASD pathogenesis.
It is worth noting that this study has several limitations. Firstly, the relatively small sample size may not only limit the statistical power and generalizability of the findings but also increase potential false-positive risks. Although sequencing depth was sufficient, the diversity estimates may still be influenced by sample heterogeneity. Secondly, we did not conduct grouping studies on children with ASD in different disease severity and stages, which may also affect gut virome composition. Thirdly, the study only focused on gut viral communities, we did not explore interactions between viruses and other microbiota components such as bacteria, fungi, or archaea. These cross-domain interactions may play a critical role in shaping gut ecology and influencing ASD pathophysiology. The lack of integrated bacteriome data limits the interpretation of virus-bacteria interactions, such as phage-bacterial host dynamics and potential changes in microbial function. In the future, the sample size can be expanded through multi-center research, and group studies can be carried out on children with ASD at different disease degrees and periods to reduce sample deviation and improve the reliability, accuracy and clinical applicability of the research conclusions. Meanwhile, studies on the gut virome and intestinal flora can also be conducted to fill the shortcomings of a single study of viruses or intestinal flora, thereby providing more insights for understanding the mechanisms of ASD and developing new therapies.
In conclusion, this study compared the gut viral communities between ASD patients and healthy controls. Our study revealed a significant difference in the abundance of Microviridae between the two groups. The study also confirmed that the alpha diversity of the ASD patients group was lower than that of the healthy controls group, and the beta diversity of the two groups was significantly different. Since there are some limitations, our study deepens the understanding of the differences in the gut viral communities between between ASD and healthy children and contributes to further research on the characteristics and roles of the gut viruses in children with ASD.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics statement
The studies involving humans were approved by the Ethics Committee of the Affiliated Hospital of Jiangsu University. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the participants’ legal guardians/next of kin.
Author contributions
MY: Writing – original draft, Investigation, Visualization, Data curation, Methodology. QW: Writing – original draft, Visualization, Data curation, Investigation, Methodology. YL: Writing – original draft, Data curation. PX: Data curation, Writing – original draft. CP: Data curation, Writing – original draft. WZ: Conceptualization, Resources, Supervision, Writing – review & editing, Funding acquisition.. HL: Conceptualization, Resources, Supervision, Writing – review & editing, Funding acquisition.
Funding
The author(s) declare financial support was received for the research and/or publication of this article. This research was supported by National Key Research and Development Programs of China (No. 2022YFC2603801), Natural Science Foundation of China (No. 82341106), and Scientific Research Project of Health Commission of Jiangsu Province (M2022043).
Acknowledgments
We thank all authors for their contribution to this study and all participants for participation.
Conflict of interest
The authors affirm that this study was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2025.1660970/full#supplementary-material
Supplementary Table 1Baseline characteristics of ASD and healthy children.
Supplementary Table 2Library information.
Supplementary Table 3The summary of the percentage information of viruses at the family level.
Supplementary Figure 1The diversity of viral species in the 22 libraries. (A) Species rarefaction curve drawn in Megan6 software after log-scale transformation. Legends are presented at the bottom of the figure. Each legend represents a library whose color corresponds to the color of the curve. (B) Accumulation curve of viral species in viral metagenomes. The X axis represents the number of sampled libraries, and the Y axis represents the cumulative number of species found. The light blue area represents the 95% confidence interval.
References
1
AlmandE. A.MooreM. D.JaykusL. A. (2017). Virus-bacteria interactions: an emerging topic in human infection. Viruses9. doi:Â 10.3390/v9030058
2
Al-ObaidiM. M. J.BahadoranA.WangS. M.ManikamR.RajuC. S.SekaranS. D. (2018). Disruption of the blood brain barrier is vital property of neurotropic viral infection of the central nervous system. Acta Virol.62, 16–27. doi: 10.4149/av_2018_102
3
AriasC. F.DuBoisR. M. (2017). The astrovirus capsid: A review. Viruses9. doi:Â 10.3390/v9010015
4
BoschA.PintóR. M.GuixS. (2014). Human astroviruses. Clin. Microbiol. Rev.27, 1048–1074. doi: 10.1128/cmr.00013-14
5
Conceição-NetoN.ZellerM.LefrèreH.De BruynP.BellerL.DeboutteW.et al. (2015). Modular approach to customise sample preparation procedures for viral metagenomics: a reproducible protocol for virome analysis. Sci. Rep.5, 16532. doi: 10.1038/srep16532
6
CorettiL.PaparoL.RiccioM. P.AmatoF.CuomoM.NataleA.et al. (2018). Gut microbiota features in young children with autism spectrum disorders. Front. Microbiol.9. doi:Â 10.3389/fmicb.2018.03146
7
CryanJ. F.O'RiordanK. J.SandhuK.PetersonV.DinanT. G. (2020). The gut microbiome in neurological disorders. Lancet Neurol.19, 179–194. doi: 10.1016/s1474-4422(19)30356-4
8
DebS.MondalR.LahiriD.ShomeG.RoyA. G.SarkarV.et al. (2023). Norovirus-associated neurological manifestations: summarizing the evidence. J. Neurovirol29, 492–506. doi: 10.1007/s13365-023-01152-0
9
De BenedictisP.Schultz-CherryS.BurnhamA.CattoliG. (2011). Astrovirus infections in humans and animals - molecular biology, genetic diversity, and interspecies transmissions. Infect. Genet. Evol.11, 1529–1544. doi: 10.1016/j.meegid.2011.07.024
10
DelmasB.AttouiH.GhoshS.MalikY. S.MundtE.VakhariaV. N.et al. (2019). ICTV virus taxonomy profile: Picobirnaviridae. J. Gen. Virol.100, 133–134. doi: 10.1099/jgv.0.001186
11
DengX.NaccacheS. N.NgT.FedermanS.LiL.ChiuC. Y.et al. (2015). An ensemble strategy that significantly improves de novo assembly of microbial genomes from metagenomic next-generation sequencing data. Nucleic Acids Res.43, e46. doi:Â 10.1093/nar/gkv002
12
FinnR. D.ClementsJ.EddyS. R. (2011). HMMER web server: interactive sequence similarity searching. Nucleic Acids Res.39, W29–W37. doi: 10.1093/nar/gkr367
13
FulciV.StronatiL.CucchiaraS.LaudadioI.CarissimiC. (2021). Emerging roles of gut virome in pediatric diseases. Int. J. Mol. Sci.22. doi:Â 10.3390/ijms22084127
14
GaneshB.MasachessiG.MladenovaZ. (2014). Animal picobirnavirus. Virusdisease25, 223–238. doi: 10.1007/s13337-014-0207-y
15
GhoshA.Delgado-CunninghamK.LópezT.GreenK.AriasC. F.DuBoisR. M. (2024). Structure and antigenicity of the divergent human astrovirus VA1 capsid spike. PloS Pathog.20, e1012028. doi: 10.1371/journal.ppat.1012028
16
GreenK. Y.AndoT.BalayanM. S.BerkeT.ClarkeI. N.EstesM. K.et al. (2000). Taxonomy of the caliciviruses. J. Infect. Dis.181 Suppl 2, S322–S330. doi: 10.1086/315591
17
GregoryA. C.ZablockiO.ZayedA. A.HowellA.BolducB.SullivanM. B. (2020). The gut virome database reveals age-dependent patterns of virome diversity in the human gut. Cell Host Microbe28, 724–740.e728. doi: 10.1016/j.chom.2020.08.003
18
Gutiérrez-CamusA.ValverdeE.Del RosalT.UtrillaC.Ruiz-CarrascosoG.Falces-RomeroI.et al. (2022). Norovirus-associated white matter injury in a term newborn with seizures. Pediatr. Infect. Dis. J.41, 917–918. doi: 10.1097/inf.0000000000003677
19
HughesH. K.MorenoR. J.AshwoodP. (2023). Innate immune dysfunction and neuroinflammation in autism spectrum disorder (ASD). Brain Behav. Immun.108, 245–254. doi: 10.1016/j.bbi.2022.12.001
20
HusonD. H.AuchA. F.QiJ.SchusterS. C. (2007). MEGAN analysis of metagenomic data. Genome Res.17, 377–386. doi: 10.1101/gr.5969107
21
Iglesias-VázquezL.Van Ginkel RibaG.ArijaV.CanalsJ. (2020). Composition of gut microbiota in children with autism spectrum disorder: A systematic review and meta-analysis. Nutrients12. doi: 10.3390/nu12030792
22
Integrative HMP (iHMP) Research Network Consortium. (2019). The integrative human microbiome project. Nature569, 641–648. doi: 10.1038/s41586-019-1238-8
23
JohnsonL. S.EddyS. R.PortugalyE. (2010). Hidden Markov model speed heuristic and iterative HMM search procedure. BMC Bioinf.11, 431. doi:Â 10.1186/1471-2105-11-431
24
KashnikovA. Y.EpifanovaN. V.NovikovaN. A. (2023). On the nature of picobirnaviruses. Vavilovskii Zhurnal Genet. Selektsii27, 264–275. doi: 10.18699/vjgb-23-32
25
KearseM.MoirR.WilsonA.Stones-HavasS.CheungM.SturrockS.et al. (2012). Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics28, 1647–1649. doi: 10.1093/bioinformatics/bts199
26
KimY. S.LeventhalB. L. (2015). Genetic epidemiology and insights into interactive genetic and environmental effects in autism spectrum disorders. Biol. Psychiatry77, 66–74. doi: 10.1016/j.biopsych.2014.11.001
27
KimuraE.GotoH.MigitaA.HaradaS.YamashitaS.HiranoT.et al. (2010). An adult norovirus-related encephalitis/encephalopathy with mild clinical manifestation. BMJ Case Rep.2010. doi:Â 10.1136/bcr.03.2010.2784
28
KumarS.StecherG.LiM.KnyazC.TamuraK. (2018). MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol.35, 1547–1549. doi: 10.1093/molbev/msy096
29
LiuS.LiE.SunZ.FuD.DuanG.JiangM.et al. (2019b). Altered gut microbiota and short chain fatty acids in Chinese children with autism spectrum disorder. Sci. Rep.9, 287. doi:Â 10.1038/s41598-018-36430-z
30
LiuF.LiJ.WuF.ZhengH.PengQ.ZhouH. (2019a). Altered composition and function of intestinal microbiota in autism spectrum disorders: a systematic review. Transl. Psychiatry9, 43. doi:Â 10.1038/s41398-019-0389-6
31
LiuZ.YangS.WangY.ShenQ.YangY.DengX.et al. (2016). Identification of a novel human papillomavirus by metagenomic analysis of vaginal swab samples from pregnant women. Virol. J.13, 122. doi:Â 10.1186/s12985-016-0583-6
32
LordC.ElsabbaghM.BairdG.Veenstra-VanderweeleJ. (2018). Autism spectrum disorder. Lancet392, 508–520. doi: 10.1016/s0140-6736(18)31129-2
33
LuX.HuaX.WangY.ZhangD.JiangS.YangS.et al. (2021b). Comparison of gut viral communities in diarrhoea and healthy dairy calves. J. Gen. Virol.102. doi:Â 10.1099/jgv.0.001663
34
LuL.ZhongH.XuM.SuL.CaoL.JiaR.et al. (2021a). Molecular and epidemiological characterization of human adenovirus and classic human astrovirus in children with acute diarrhea in Shanghai 2017-2018. BMC Infect. Dis.21, 713. doi:Â 10.1186/s12879-021-06403-1
35
LuoZ.SuR.WangW.LiangY.ZengX.ShereenM. A.et al. (2019). EV71 infection induces neurodegeneration via activating TLR7 signaling and IL-6 production. PloS Pathog.15, e1008142. doi:Â 10.1371/journal.ppat.1008142
36
MartinC. R.OsadchiyV.KalaniA.MayerE. A. (2018). The brain-gut-microbiome axis. Cell Mol. Gastroenterol. Hepatol.6, 133–148. doi: 10.1016/j.jcmgh.2018.04.003
37
MarvinS. A.HuertaC. T.SharpB.FreidenP.ClineT. D.Schultz-CherryS. (2016). Type I interferon response limits astrovirus replication and protects against increased barrier permeability in vitro and in vivo. J. Virol.90, 1988–1996. doi: 10.1128/jvi.02367-15
38
Mayneris-PerxachsJ.Castells-NobauA.Arnoriaga-RodrÃguezM.Garre-OlmoJ.PuigJ.RamosR.et al. (2022). Caudovirales bacteriophages are associated with improved executive function and memory in flies, mice, and humans. Cell Host Microbe30, 340–356.e348. doi: 10.1016/j.chom.2022.01.013
39
MeltzerA.Van de WaterJ. (2017). The role of the immune system in autism spectrum disorder. Neuropsychopharmacology42, 284–298. doi: 10.1038/npp.2016.158
40
MilaniC.DurantiS.BottaciniF.CaseyE.TurroniF.MahonyJ.et al. (2017). The first microbial colonizers of the human gut: composition, activities, and health implications of the infant gut microbiota. Microbiol. Mol. Biol. Rev.81. doi:Â 10.1128/mmbr.00036-17
41
MorellaN. M.GomezA. L.WangG.LeungM. S.KoskellaB. (2018). The impact of bacteriophages on phyllosphere bacterial abundance and composition. Mol. Ecol.27, 2025–2038. doi: 10.1111/mec.14542
42
Morris-RosendahlD. J.CrocqM. A. (2020). Neurodevelopmental disorders-the history and future of a diagnostic concept. Dialogues Clin. Neurosci.22, 65–72. doi: 10.31887/DCNS.2020.22.1/macrocq
43
PapeK.TamouzaR.LeboyerM.ZippF. (2019). Immunoneuropsychiatry - novel perspectives on brain disorders. Nat. Rev. Neurol.15, 317–328. doi: 10.1038/s41582-019-0174-4
44
PortoP. S.RiveraA.MoonrintaR.WobusC. E. (2023). Entry and egress of human astroviruses. Adv. Virus Res.117, 81–119. doi: 10.1016/bs.aivir.2023.08.001
45
PosarA.RescaF.ViscontiP. (2015). Autism according to diagnostic and statistical manual of mental disorders 5(th) edition: The need for further improvements. J. Pediatr. Neurosci.10, 146–148. doi: 10.4103/1817-1745.159195
46
QuA.DuanB.WangY.CuiZ.ZhangN.WuD. (2023). Children with autism show differences in the gut DNA virome compared to non-autistic children: a case control study. BMC Pediatr.23, 174. doi:Â 10.1186/s12887-023-03981-8
47
ReyesA.SemenkovichN. P.WhitesonK.RohwerF.GordonJ. I. (2012). Going viral: next-generation sequencing applied to phage populations in the human gut. Nat. Rev. Microbiol.10, 607–617. doi: 10.1038/nrmicro2853
48
RochaL. C. D.EstofoleteC. F.MilhimB.AugustoM. T.ZiniN.SilvaG.et al. (2021). Enteric viruses circulating in undiagnosed central nervous system infections at tertiary hospital in São José do Rio Preto, São Paulo, Brazil. J. Med. Virol.93, 3539–3548. doi: 10.1002/jmv.26216
49
RoggeN.JanssenJ. (2019). The economic costs of autism spectrum disorder: A literature review. J. Autism Dev. Disord.49, 2873–2900. doi: 10.1007/s10803-019-04014-z
50
RonquistF.TeslenkoM.van der MarkP.AyresD. L.DarlingA.HöhnaS.et al. (2012). MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst. Biol.61, 539–542. doi: 10.1093/sysbio/sys029
51
RutschA.KantsjöJ. B.RonchiF. (2020). The gut-brain axis: how microbiota and host inflammasome influence brain physiology and pathology. Front. Immunol.11. doi: 10.3389/fimmu.2020.604179
52
SatoM.KurodaM.KasaiM.MatsuiH.FukuyamaT.KatanoH.et al. (2016). Acute encephalopathy in an immunocompromised boy with astrovirus-MLB1 infection detected by next generation sequencing. J. Clin. Virol.78, 66–70. doi: 10.1016/j.jcv.2016.03.010
53
SharonG.CruzN. J.KangD. W.GandalM. J.WangB.KimY. M.et al. (2019). Human gut microbiota from autism spectrum disorder promote behavioral symptoms in mice. Cell177, 1600–1618.e1617. doi: 10.1016/j.cell.2019.05.004
54
ShkoporovA. N.ClooneyA. G.SuttonT. D. S.RyanF. J.DalyK. M.NolanJ. A.et al. (2019). The human gut virome is highly diverse, stable, and individual specific. Cell Host Microbe26, 527–541.e525. doi: 10.1016/j.chom.2019.09.009
55
ShkoporovA. N.HillC. (2019). Bacteriophages of the human gut: the "Known unknown" of the microbiome. Cell Host Microbe25, 195–209. doi: 10.1016/j.chom.2019.01.017
56
ShuidA. N.JayusmanP. A.ShuidN.IsmailJ.Kamal NorN.MohamedI. N. (2021). Association between viral infections and risk of autistic disorder: an overview. Int. J. Environ. Res. Public Health18. doi:Â 10.3390/ijerph18062817
57
Skewes-CoxP.SharptonT. J.PollardK. S.DeRisiJ. L. (2014). Profile hidden Markov models for the detection of viruses within metagenomic sequence data. PloS One9, e105067. doi:Â 10.1371/journal.pone.0105067
58
SuQ.WongO. W. H.LuW.WanY.ZhangL.XuW.et al. (2024). Multikingdom and functional gut microbiota markers for autism spectrum disorder. Nat. Microbiol.9, 2344–2355. doi: 10.1038/s41564-024-01739-1
59
TomovaA.HusarovaV.LakatosovaS.BakosJ.VlkovaB.BabinskaK.et al. (2015). Gastrointestinal microbiota in children with autism in Slovakia. Physiol. Behav.138, 179–187. doi: 10.1016/j.physbeh.2014.10.033
60
TongZ.ZhouX.ChuY.ZhangT.ZhangJ.ZhaoX.et al. (2022). Implications of oral streptococcal bacteriophages in autism spectrum disorder. NPJ Biofilms Microbiomes8, 91. doi:Â 10.1038/s41522-022-00355-3
61
VemuriR.GundamarajuR.ShastriM. D.ShuklaS. D.KalpurathK.BallM.et al. (2018). Gut microbial changes, interactions, and their implications on human lifecycle: an ageing perspective. BioMed. Res. Int.2018, 4178607. doi:Â 10.1155/2018/4178607
62
VemuriR.ShankarE. M.ChieppaM.EriR.KavanaghK. (2020). Beyond just bacteria: functional biomes in the gut ecosystem including virome, mycobiome, archaeome and helminths. Microorganisms8. doi:Â 10.3390/microorganisms8040483
63
VinjéJ.EstesM. K.EstevesP.GreenK. Y.KatayamaK.KnowlesN. J.et al. (2019). ICTV virus taxonomy profile: caliciviridae. J. Gen. Virol.100, 1469–1470. doi: 10.1099/jgv.0.001332
64
WanY.WongO. W. H.TunH. M.SuQ.XuZ.TangW.et al. (2024a). Fecal microbial marker panel for aiding diagnosis of autism spectrum disorders. Gut Microbes16, 2418984. doi:Â 10.1080/19490976.2024.2418984
65
WanY.ZhangL.XuZ.SuQ.LeungT. F.ChanD.et al. (2024b). Alterations in fecal virome and bacteriome virome interplay in children with autism spectrum disorder. Cell Rep. Med.5, 101409. doi:Â 10.1016/j.xcrm.2024.101409
66
WangH.LingY.ShanT.YangS.XuH.DengX.et al. (2019). Gut virome of mammals and birds reveals high genetic diversity of the family Microviridae. Virus Evol.5, vez013. doi:Â 10.1093/ve/vez013
67
WenQ.WangJ.DaiL.MomingA.FanZ.HuangY.et al. (2024). Viromes of three phyla of aquatic organisms in the south China sea. J. Zoonoses4. doi:Â 10.15212/zoonoses-2024-0041
68
WillyardC. (2021). How gut microbes could drive brain disorders. Nature590, 22–25. doi: 10.1038/d41586-021-00260-3
69
WuJ.ChaiT.ZhangH.HuangY.PerryS. W.LiY.et al. (2022). Changes in gut viral and bacterial species correlate with altered 1,2-diacylglyceride levels and structure in the prefrontal cortex in a depression-like non-human primate model. Transl. Psychiatry12, 74. doi:Â 10.1038/s41398-022-01836-x
70
XiaoY.WangH.FengL.PanJ.ChenZ.WangH.et al. (2020). Fecal, oral, blood and skin virome of laboratory rabbits. Arch. Virol.165, 2847–2856. doi: 10.1007/s00705-020-04808-y
71
ZeidanJ.FombonneE.ScorahJ.IbrahimA.DurkinM. S.SaxenaS.et al. (2022). Global prevalence of autism: A systematic review update. Autism Res.15, 778–790. doi: 10.1002/aur.2696
72
ZengC.WanS. R.GuoM.TanX. Z.ZengY.WuQ.et al. (2024). Fecal virome transplantation: A promising strategy for the treatment of metabolic diseases. BioMed. Pharmacother177, 117065. doi:Â 10.1016/j.biopha.2024.117065
73
ZhangM.ChuY.MengQ.DingR.ShiX.WangZ.et al. (2020). A quasi-paired cohort strategy reveals the impaired detoxifying function of microbes in the gut of autistic children. Sci. Adv.6. doi:Â 10.1126/sciadv.aba3760
74
ZhangW.LiL.DengX.BlümelJ.NüblingC. M.HunfeldA.et al. (2016). Viral nucleic acids in human plasma pools. Transfusion56, 2248–2255. doi: 10.1111/trf.13692
75
ZhangW.LiL.DengX.KapusinszkyB.PesaventoP. A.DelwartE. (2014). Faecal virome of cats in an animal shelter. J. Gen. Virol.95, 2553–2564. doi: 10.1099/vir.0.069674-0
76
ZhangW.YangS.ShanT.HouR.LiuZ.LiW.et al. (2017). Virome comparisons in wild-diseased and healthy captive giant pandas. Microbiome5, 90. doi:Â 10.1186/s40168-017-0308-0
77
ZouR.XuF.WangY.DuanM.GuoM.ZhangQ.et al. (2020). Changes in the gut microbiota of children with autism spectrum disorder. Autism Res.13, 1614–1625. doi: 10.1002/aur.2358
Summary
Keywords
autism spectrum disorder, gut viral communities, virome, diversity, viral metagenomics
Citation
Yuan M, Wang Q, Lu Y, Xu P, Pan C, Zhang W and Lu H (2025) Comparison of gut viral communities between autism spectrum disorder and healthy children. Front. Cell. Infect. Microbiol. 15:1660970. doi: 10.3389/fcimb.2025.1660970
Received
07 July 2025
Accepted
26 September 2025
Published
14 October 2025
Volume
15 - 2025
Edited by
Stefano Marini, Independent Researcher, Termoli, Italy
Updates
Copyright
© 2025 Yuan, Wang, Lu, Xu, Pan, Zhang and Lu.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hongyan Lu, hylu@ujs.edu.cn; Wen Zhang, z0216wen@yahoo.com
†These authors have contributed equally to this work and share first authorship
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.