Alpha-synuclein alters the faecal viromes of rats in a gut-initiated model of Parkinson’s disease

Parkinson’s disease (PD) is a chronic neurological disorder associated with the misfolding of alpha-synuclein (α-syn) into aggregates within nerve cells that contribute to their neurodegeneration. Recent evidence suggests α-syn aggregation may begin in the gut and travel to the brain along the vagus nerve, with microbes potentially a trigger initiating α-syn misfolding. However, the effects α-syn alterations on the gut virome have not been investigated. In this study, we show longitudinal faecal virome changes in rats administered either monomeric or preformed fibrils (PFF) of α-syn directly into their enteric nervous system. Differential changes in rat viromes were observed when comparing monomeric and PFF α-syn, with alterations compounded by the addition of LPS. Changes in rat faecal viromes were observed after one month and did not resolve within the study’s five-month observational period. These results suggest that virome alterations may be reactive to host α-syn changes that are associated with PD development.

P arkinson's disease (PD) is the second most common progressive neurodegenerative disorder, affecting 0.1-0.2% of the population at any given time. This frequency increases with age, with more than 1% of individuals over 60 living with PD 1 . While the overall cost of PD varies from country to country, in the UK conservative and liberal estimates calculated the annual cost of the illness between £449 million and £3.3 billion, respectively 2 . As the mean age of populations increase, the social and economic burden of PD will rise. In fact, worldwide cases of PD are expected to double to 14 million by 2040 3 ; therefore, a greater understanding of the disease and new interventions are urgently required.
Genetic and epidemiological studies investigating PD have not demonstrated a specific aetiology that triggers PD. Instead, many hypothesise that PD is a multifactorial disease involving a complex combination of host genetics and environmental factors 4 . The pathophysiological onset of PD is associated with the loss of dopamine-producing neurons in the pars compacta subregion of the substantia nigra structure of the midbrain 5 . However, patients typically only display Parkinsonism, the motor features associated with PD, after 50-80% of dopaminergic neurons are lost 6 . The presynaptic neuronal protein alpha-synuclein (α-syn) is linked to the loss of dopaminergic neurons, whereby the monomeric form of α-syn aggregates within nerve cells into toxic Lewy body (LB) formations that disrupt cellular homoeostasis 7 .
Gastrointestinal (GI) comorbidities are frequently associated with PD, with cohorts from independent studies reporting that 70-78% of patients experience excessive saliva, 30-97% of individuals have difficulty swallowing and 20-89% report constipation 8 . Non-motor symptoms such as GI comorbidities associated with PD may also involve α-syn, as its aggregation is detected earlier and more frequently in the enteric (i.e. GI) nervous system of PD patients 9,10 . Recent animal models have also demonstrated the spread of pathologic α-syn from the gut to the brain through the vagus nerve [11][12][13][14] . This has led researchers to investigate the GI tract and its community of microbes as a potential environmental trigger for the aggregation of α-syn that may contribute to the development of PD.
There are 10 13 to 10 14 bacteria, archaea, fungi and viruses associated with the human body, known collectively as the microbiota. Non-culture-based approaches to study their genomic material (microbiome) have demonstrated the gut microbiome of PD patients differs from controls [15][16][17][18][19][20][21][22] . In addition, transplantation of human faecal microbiotas from either PD patients or non-PD controls to mice resulted in a motor function decline amongst the mice humanised with a PD microbiota 23 . Specific bacterial taxa, and even bacterial-infecting viruses (termed bacteriophages or phages), have been proposed as potential contributors to PD 19,24 . However, caution is often required when interpreting microbial changes associated with specific conditions, as alterations observed in the microbiome may be causal, or reactive, to human host-associated changes 25 .
This study was specifically performed in young rats to directly test if α-syn affects the virome at a younger age-span between 2 and 7 months, where 7 months is estimated to be equivalent to~18-20 years in humans 26 . Previous studies have shown that α-syn levels increased significantly in the gut of paediatric patients and correlated with infection by microbial pathogens 27 . Increases in α-syn levels at early ages is known to precede PD motor symptoms by decades. We thus rationalised that it would be most informative to understand how α-syn affects the virome in a younger healthier context without the confounding variable of age-related alterations in health and the microbiome. Furthermore, PD motor symptoms can present as young as 28 years due to triplication of the α-syn gene (SNCA) 28 . Thus, an increase in the expression of α-synuclein alone, as introduced in the gut in this study, is sufficient to cause PD with a younger age of onset. Finally, as we administered LPS as an adjunct treatment, we were concerned that this may have a detrimental impact on older animals compared to a younger cohort.
In this study, we investigated if the viral component of the rat microbiome (virome) changes in response to α-syn injected into rat enteric nervous systems. Both the monomeric and aggregated fibrillar forms of α-syn were injected directly into the myenteric plexus of surgically-prepared rat duodenal walls, with and without intraperitoneal (i.p.) injection of bacterial lipopolysaccharide (LPS) as a potential inflammatory adjunct. Virome diversity changes were observed between rats receiving specific α-syn treatments, and also longitudinally across individual rats. Specific viral contigs and viral clusters (VCs) were noted as differentially abundant across treatment groups, with microbiome changes reacting to α-syn and potentially PD inducing conditions.

Results
Defining the rat faecal virome. Viral contigs of the rat microbiome were identified using database-dependent and -independent approaches, applied to both the viral-enriched and wholegenome shotgun (WGS) sequencing data (see Methods). Viral sequences identified in the viral-enriched and WGS data were merged into a single database and made nonredundant by removing duplicate sequences and retaining the larger of any two homologous sequences (≥90% identity and ≥90% coverage). The majority of sequences in the final rat viral database originated from the viral-enriched fraction of the microbiome (3506/ 4009 sequences). The largest viral contig was 291,209 bp in length, while the mean and median length of viral contigs is 4123 and 1621 bp, respectively (Fig. 1a).
The viral nucleic acid extraction and sequencing protocols used in this study did not employ random DNA amplification, such as multiple displacement amplification (MDA), known to bias the viral composition towards small circular genetic elements. Therefore, the relative abundance of specific viral taxa such as small single-stranded DNA Microviridae may be more accurate than previous faecal virome studies 29 . The average raw sequencing reads per sample for rat faecal viromes was 2,926,309 (max. 3,616,794 and min. 2,221,657). The rat faecal viromes analysed in this study were dominated by bacterial-infecting phage. However, the term virome (and not 'phageome') is used throughout this study, as this study's viral-detecting protocol does not preclude the identification of rat-infecting viruses and these may be present amongst the 'unclassified' viruses. The absence of eukaryotic viruses is predicted to reflect the dominance of prokaryotetargeting phages in the faeces of healthy animals that likely requires greater sequencing depth to capture the total virome and the use of chloroform in the concentration of virus-like particles that would diminish membrane-containing viruses (see Methods).
Of the tailed Caudovirales phages assigned a putative familial taxonomic rank, Siphoviridae phages were the most abundant, with Myoviridae and Podoviridae phages detected at decreasing frequencies (Fig. 1b). The accumulative relative abundance of these tailed phages followed the same trend ( Supplementary  Fig. 1). All viral sequences identified within this study were clustered into pseudogenera termed VCs using the vConTACT2 gene sharing network-based algorithm. A total of 920 VCs were generated, with the largest cluster containing 20 sequences.
The majority of phage contigs identified in this study are linear (99.45%; Fig. 1c). While the rat faecal virome might have relatively few circular genomes, it is more likely that greater sequencing depth or an alternative strategy is required to complete and circularise more viral genome sequences. For instance, a hybrid assembly approach combining short-read and long-read sequencing technologies recently overcame viral population microdiversity to assemble longer and more circular contigs 30 . The viral sequences identified in this study were queried for the presence of phage integrases and/or recombinases indicative of lysogenic replication. Only 8.4% of contigs possessed genes encoding for lysogenic conversion that would facilitate integration and replication as a prophage (Fig. 1d).
Using a database-dependent and -independent approach to identify viral contigs, the average relative abundance of the identified rat faecal viruses accounts for 58.9% of the virome sequencing reads generated ( Fig. 1e; min. 47.6%, max. 70.2%). Therefore, a large proportion of the rat faecal virome data is not analysed and could be considered as 'viral dark matter' 31 . While significantly lower than within the whole shotgun sequencing, there are instances where the bacterial housekeeping gene, cpn60, is detected within the virome sequencing data (Fig. 1f). Therefore, potential bacterial contamination may have been inadvertently categorised as viral, and a greater understanding of the rat microbiome will be required to distinguish the origin of sequences.
The rat microbiome changes over time. The male Sprague-Dawley rats employed in this study were~2 months of age at the study's onset and 7 months old at its conclusion. Therefore, the potential impact of time in shaping the rat faecal microbiome was investigated. During the MetaPhlAn2 analysis of the WGS data, an average of 99.7 and 0.3% of the microbiome sequencing reads that could be assigned to microbial taxa were classified as bacterial and viral, respectively. No statistical difference was observed in the overall relative abundance of lysogenic phages over the course of the study ( Supplementary Fig. 2a). A longitudinal axis to the rat faecal virome was obtained by aligning the WGS sequencing data onto the final rat viral database. An average of 44,431 WGS sequencing reads aligned to viral sequences (min. 572; max. 256,910; median 33,298).
There was a statistically significant change in the faecal virome and total microbiome composition of the rats over the study's duration, with time explaining 14.6 and 13.3% of their variances (R 2 ), respectively (Fig. 2a, b). The two-dimensional separation of faecal samples by dissimilarity (β-diversity) was not driven by intra-sample composition (α-diversity). No statistical difference in α-diversities after Bonferroni correction are observed in rat faecal viromes by treatment group, examining either virome or WGS data (Supplementary Table 1). The PCoA axes 1 and 2 for the microbiome data explain more of the variation than for the virome data (28.5 and 15.5% versus 21.1 and 9.4%, respectively).
α-syn alters the rat faecal virome. As it is established that the rat faecal virome changed over the course of the study, any alterations occurring in α-syn injected rats need to be examined in light of the Sham and LPS control treatment groups. The Sham and Fig. 1 The rat faecal virome. a Contig length (bp) versus sequencing coverage of identifiable viruses present in the faeces of male rats. The shape aesthetic was used to highlight the topology of assembled viral contigs, while colour represents the proposed taxonomic assignment. b The breakdown of rat faecal viruses into their putative familial-level taxonomic assignments. c The topology of identified viral contigs, whereby circular sequences are most likely complete genomes. d Viral contigs encoding a phage integrase and/or recombinase were annotated as lysogenic phages. e The relative abundance of viral contigs within the virome sequencing data, aggregated by taxonomic assignment. Sequencing reads aligning to contigs not recognised as viral are omitted. f The relative abundance of the bacterial housekeeping gene cpn60 detected in the virome sequencing data versus the WGS data.
LPS controls experienced 22.6 and 19.5% variance, respectively, in their viromes over the study's duration (Fig. 3a, b). Less than a 1% difference in the variance of the rat faecal viral compositions are observed in the rats that received the α-syn monomer or α-syn monomer with LPS (28.0 and 27.6% variances explained, respectively; Fig. 3c, d). However, the effect of the α-syn monomer on rat viromes, with or without LPS, is far greater than differences attributable to either the Sham or LPS control treatments. The PFF form of α-syn appears to have had only a small effect on the rat virome relative to the Sham and LPS controls; however, the PFF plus LPS treatments caused the greatest virome alteration (variances explained 23.9 and 30.8%, respectively; Fig. 3e, f). There is an overall trend for the relative abundance of lysogenic phages to increase in rats treated with α-syn, monomer or fibrillar, yet only the monomer plus LPS and PFF treatments were statistically significant from the Sham control group (Supplementary Fig. 2b).
Changes in individual rat viromes. In order to investigate if changes in the rat faecal virome were directional, whereby a sample is not merely more or less similar to its baseline through random fluctuations, each individual rat's faecal virome was investigated relative to its baseline state. Using Bray-Curtis index values that indicate the dissimilarity of two given virome samples, it was observed that the greatest variation in the rat faecal virome occurs between Base and T1, compared to T1 and T5 (Fig. 4a). There was no statistical difference in the dissimilarity values between the 'Base versus T1' and 'Base versus T5'. Additionally, the virome at 5 months was more frequently dissimilar to the baseline compared to the 1-month time point (Fig. 4b).
For the individual rat viromes segregated by treatment conditions (Fig. 4c-h), the least variation in Bray-Curtis dissimilarity values between time points was observed between the average Sham, LPS, and α-syn monomer treatment baselines versus T1 (Fig. 4c-e). Whereas within the first month, the α-syn monomer plus LPS, the PFF α-syn and the PFF α-syn plus LPS treatments had the least similarity when averaged across rats ( Fig. 4f-h). The viromes of rats receiving the α-syn monomer plus LPS are the most dissimilar to their baseline after the study's 5 months (Fig. 4f).
Additionally, we investigated if virome changes associated with specific α-syn treatments were the result of cage effects. Each rat was assigned a numerical value and a code based on its cage. For example, code C21-R2 specifies cage number 21 and rat 2. Two cages, with two-three rats per cage, were used for each treatment. Looking at each individual rat, no specific trend across any of the treatment conditions was observed, such as all the rats within a specific cage dramatically differing at T1 or T5 relative to their baseline ( Fig. 4c-h). Also, the average Bray-Curtis dissimilarity values were similar for each specific treatment across the two corresponding cages.
Viral changes caused by α-syn treatments. If the α-syn treatments investigated in this study altered the rat microbiome and/or virome, it is expected that the abundance of viral contigs would be significantly altered in one or more time points (p value < 1E-05). The Sham and LPS treatment arms of this study resulted in the fewest differential viral and VC changes (Fig. 5a, b). There were only two viral contigs increased at the 1-and 5-month time points of the Sham-administered rats, while no viral contigs were consistently altered amongst LPS treated rats.
For rats receiving the monomeric form of α-syn, 16 and six viral contigs had decreased and increased differential abundance at the 5-month time point (Fig. 5c). The α-syn monomer treatment did not induce a change in a viral contig or virus cluster across all the time points investigated. The α-syn monomer plus LPS, the PFF and the PFF plus LPS treatment arms of the study all induced significant abundance changes in viral contigs or VCs consistent across T1 and T5 time points. These treatments caused decreased differential abundance in six, two and five viral contigs and VCs, respectively ( Fig. 5d-f). While the α-syn monomer plus LPS and the PFF treatments both increased two viral groups across both time points (Fig. 5d, e).
Common viral changes across PD inducing treatments. In order to discern if a common pattern of changes occurs across different treatment conditions, all viral contigs and VCs that were noted to have a differential abundance at any time point were compared across the treatments tested. A total of 78 viral contigs and VCs (36 and 42, respectively), predicted from the virome sequencing and the WGS sequencing data, were recorded as differentially abundant at one or more time points relative to their baseline (Fig. 6). Of these, the majority (48 of 78, 61.5%) were only altered in one treatment condition. These are most probably the result of random changes that have occurred in the rat microbiome over the course of the study.
Nine viral contigs or VCs were altered across both the T1 and T5 time points analysed, supporting the hypothesis that these are prolonged microbiome changes. Only four viral contigs are concordantly changed across the Sham and LPS controls, seven in  the α-syn monomer and α-syn monomer plus LPS treatments, and finally, nine in the α-syn PFF and α-syn PFF plus LPS treated rats. There was one viral contig with an increased differential abundance at T5 across all treatment groups (contig Rat_15_NO-DE_1829_length_1188_cov_3.973522). Whereas a second viral contig (C22 − R1 − T5_NODE_102_length_3426_cov_2.959953) has an increased differential abundance detected in all but the αsyn PFF treatment group. Therefore, these two contigs may be fragments of the same altered virus, with the latter sequence putative annotated as a Siphoviridae phage.

Discussion
There are a limited number of viral metagenomic analyses conducted on wild rats as they are a potential source of zoonotic infections [32][33][34] . However, despite their frequent use as models of human diseases, to our knowledge, the virome of lab rats has not been investigated (see Methods). In this study, the majority of faecal viruses analysed are putatively annotated as tailed phages belonging to the Siphoviridae, Myoviridae and Podoviridae families, with few Microviridae observed. A similar trend is observed in the accumulative relative abundance of viral taxa within microbiome sequencing data ( Supplementary Fig. 1). These top-level hierarchical viral taxa are similarly observed in the faeces of mice and humans 35,36 .
Unlike many mammalian faecal virome analyses, advances in sequencing technologies have enabled us to conduct our analysis without a nucleic acid amplification step. Therefore, the relative abundance of viral taxa and the overall diversity of viromes presented here is likely to be more accurate than previous studies. With a sequencing depth of~3 M reads per sample, fewer genomes than expected were circular. Without the random amplification bias towards small circular viral elements, and an overall better representation of the virome diversity within the microbiome, an even greater sequencing depth or a combined shortand long-read sequencing approach may be required to complete metagenomic viral genomes. On average, the abundance of  lysogenic phages did not statistically differ over the study's duration ( Supplementary Fig. 2). However, the relative abundance of phages without detectable lysogeny genes, annotated as replicating through an unknown strategy, decreased (p value < 0.05) in the monomeric α-syn plus LPS and PFF treatment groups relative to the study's Sham treatment group. This was investigated because changes in the ratio of virulent to temperate phages were observed in a comparison of healthy and inflammatory bowel disease human viromes 37 . The majority of the rat faecal virome can be considered as 'dark matter', whereby assembled sequences are unlike those deposited in reference databases. This is similar to studies of the human virome, where 40-90% of sequences were unidentifiable 38 . In this study, an average of 58.9% of the rat virome sequencing data could be aligned back to the viral database generated from the sequencing of their faecal viromes. Therefore, while a proportion of the virome sequencing data was not analysed, the remaining uncharacterised data was a cautionary measure to prevent unintentional interpretation of bacterial-contaminating sequences as viral.
The diversity of rat faecal viromes changed over the study's 5-month duration. As we are not aware of any published studies on the longitudinal rat faecal virome, it is difficult to ascertain the cause(s) for this variation. Interestingly, the magnitude of change observed for the whole microbiome and virome was similar using the same methodological procedures (variances of 14.6 and 13.3%, respectively). While this result may be coincidental, it would be logical that changes in bacterial populations would have a concomitant effect on the abundance of their viruses. One potential explanation for the shifts observed in the microbiome and virome of rats over 5 months would be their accelerated aging relative to humans. When comparing their lifespans, every 13.8 rat days approx. equates to a human year 26 . Therefore, rats were potentially first analysed in this study during puberty and reached social maturity at its conclusion. The virome changes associated with these developmental periods are poorly described across all mammalian hosts.
Knowing that the rat faecal virome changes over time, all fluctuations associated with inter-sample diversity must be interpreted cautiously. The variance (R 2 ) explained by viral βdiversities did not differ dramatically between the Sham, LPS or α-syn PFF treatments (22.6, 19.5, and 23.9%, respectively). However, large β-diversity differences occurred over the study's duration between control groups and rats administered the α-syn monomer, α-syn monomer plus LPS and α-syn PFF plus LPS treatments (28.0, 27.6, and 30.8%, respectively). Therefore, specific α-syn treatments affected the overall longitudinal β-diversity of rats differently, with the fibril form of α-syn with LPS injection as an inflammatory adjunct resulting in the greater variance.
When broken down into an investigation of individual rat viromes, the greatest variation was observed between the baseline and the first month, T1. However, the viromes of rats continued to change between T1 and T5, and more frequently the viromes were more dissimilar from their basal state at the 5-month time point, T5. The magnitude by which individual rat viromes transformed over the study was smallest for the Sham, LPS and αsyn monomer treatments. Once more, the greatest changes observed in individual rats occurred in those receiving the α-syn monomer plus LPS, the α-syn PFF and α-syn PFF plus LPS treatments. These results suggest the initial treatment of α-syn monomer plus LPS, or the fibrillar form of α-syn, caused an initial and irreversible change in the virome of rats.
The individual viral contigs and VCs responsible for the rat faecal virome changes were investigated. A similar trend to the βdiversity analysis was observed across all of the study's treatment groups. More viral contigs were differentially altered after 5 months, compared to 1 month. Amongst the Sham rats investigated, only two viral contigs of VCs were differentially altered at 1 month with the change maintained until the 5-month time point. The LPS treated control rats had no viral contig or VC altered across both time points. The abundance and magnitude of viral changes are amplified in rats treated with the α-syn monomer plus LPS, the α-syn PFF and the α-syn PFF plus LPS treated rats.
Changes in rat faecal viromes were examined for consistent patterns across treatment groups. Ideally, specific viruses would be altered across the various α-syn treated rats indicating a consistent effect on the viral composition. However, few viral contigs and VCs were differentially abundant across α-syn treatments, examining either the monomeric or fibrillar form of α-syn with and without LPS (Supplementary Fig. 3). Taken together, viral changes in the rat faecal virome may be indirectly associated with inflammatory or microbial changes associated with α-syn in the gut. However, no specific alteration in the virome could be definitively linked to the native or misfolded αsyn, despite several viral contigs differentially abundant across different treatments.
In conclusion, we present the first analysis of the lab rat faecal virome. The study focused on the effect caused by the addition of α-syn on the viral composition over 5 months. There is clear evidence that the virome of rats changed over time and following the injection of the monomer or fibrillar forms of α-syn into the enteric nervous system. The magnitude of change observed across faecal viromes was consistently greater amongst rats administer α-syn monomer plus LPS, α-syn PFF or α-syn PFF plus LPS. However, additional studies are required to clearly establish if αsyn induced virome alterations are associated with the pathogenesis or GI comorbidities of PD.

Methods
Generation of α-syn PFFs. Recombinant monomeric mouse α-syn was purified and used to generate α-syn PFFs, as previously described 39 . A Pierce LAL highcapacity endotoxin removal resin was used to minimise endotoxin. Endotoxin levels were 0.017 Unit/µg of protein as determined using a Pierce LAL endotoxin quantitation kit. Briefly, α-syn monomer, kindly provided by Laura Volpicelli-Daley, was thawed rapidly at 37°C in a 1.5 ml microcentrifuge tube, then centrifuged at 100,000 × g at 4°C for 60 min. The supernatant was removed and its protein concentration was determined using a Nanodrop spectrophotometer, using the mouse α-syn coefficient 7450 M −1 cm −1 and molecular weight~15 kDa. The supernatant was diluted in 50 mM Tris pH 7.5, 150 mM KCL to a final protein concentration of 5 mg/ml. A proportion of the solution was retained for subsequent ENS injection of α-syn monomer. The remainder was placed in a thermomixer at 37°C and shaken at 1000 rpm for 7 days. Following 7 days, the clear solution became turbid, which is indicative of PFF formation. The formation of PFFs was validated using the sedimentation assay and transmission electron microscopy as previously described 35 . Immediately prior to surgery, PFFs were sonicated at 30% power on the Misonix Sonicator 3000, 1 s on, 1 s off for 120 s total, in two 30 s sessions, in order to break the fibrils into small fragments. Transmission electron micrographs showed an aggregated protein structure, with filaments that are typical of PFFs. Measurement of fibril length after sonication revealed that >60% of the fibrils were less than 60 nm in diameter. A sedimentation assay showed that there were approximately equal amounts of α-syn PFFs in the supernatant and pellet, indicating efficient PFF formation 39 .
Animal housing and husbandry. Adult male Sprague-Dawley rats were purchased from Envigo, UK and maintained on a 12 h:12 h light:dark cycle (lights on at 08:00 h) at 21 ± 2°C and humidity (30−50%), with access to standard rat chow and water ad libitum. All experiments were conducted in accordance with the European Directive 2010/63/EU and under a project authorisation issued by the Health Products Regulatory Authority, Ireland (AE19130/P036). Animals were randomly divided into six treatment groups, as follows: (a) Sham: bovine serum albumin (BSA) injected into the duodenal wall (n = 12), (b) LPS: i.p. injection of LPS (n = 12), (c) monomer: α-syn monomer injected into the duodenal wall (n = 12), (d) monomer plus LPS: α-syn monomer injected into duodenal wall plus an immediate i.p. injection of LPS (n = 12), (e) PFF: PFF injected into the duodenal wall (n = 12), and (f) PFF plus LPS: PFF injected into duodenal wall plus an immediate i.p. injection of LPS (n = 12). Animals were housed three per cage, according to their treatment group. After 1 month, six animals from each group were culled and both gut and brain samples were retrieved for further analysis, the other six were culled at 5 months.
Administration of α-syn into the duodenal wall. Animals were anaesthetised using gaseous isoflurane (4-5% to induce anaesthesia; 2-3% to maintain anaesthesia). On the day before surgery, animals were pretreated systemically with the nonsteroidal anti-inflammatory analgesic, carprofen (Rimadyl ® , 1 ml/500 ml drinking water). Animals were maintained at a constant body temperature by the use of a heating pad. Each animal's abdomen was shaved and sterilised using betadine. The laparotomy procedure was conducted as follows. A 2.5-cm-incision was made vertically along the abdomen using a sterile blade, then scissors was used to cut through the muscle layer. Approximately 3 cm in length of the stomach and the duodenum was partially exposed to aid visualisation of the injection site. A Harvard automatic pump connected by tubing to a 50 µl Hamilton syringe with a 25-gauge needle was used to perform the injections. In order to validate the injection procedure and the transmission from duodenal myenteric neurons to the CNS, 2 µl of 1% GFP-non-toxic fluorescent cholera toxin (Thermo Fisher, catalogue number C34775) was injected as described previously 40,41 into three sites in the duodenal wall. Following this validation, 2 µl of either mouse α-syn monomer (5 µg/µl) or mouse PFF (5 µg/µl) was injected into each of these three sites. Animals in the LPS groups also received an i.p. injection of LPS (2.5 mg/kg, in 0.9% w/v sterile saline) immediately after the laparotomy (LPS: Sigma, catalogue number L2630, 0111:B4) 42 . Sham animals received an injection of BSA (2 µl per site, 5 µg/ µl) into the duodenum, as previously described 11,[43][44][45] and LPS animals only received 2.5 mg/kg LPS as i.p. injection. All animals were treated for 2 days following surgery with the nonsteroidal anti-inflammatory analgesic, carprofen (Rimadyl ® , 1 ml/500 ml drinking water) to ensure pain relief.
Faecal collection and virome shotgun sequencing. Faecal pellets were directly collected from animals at three time points; at the pre-surgical baseline, termed: baseline, 1-month post-surgery (T1) and 5 months post-surgery (T5). Samples were stored on dry ice and then frozen at −80°C. Nucleic corresponding to the viromes of each group of rats was purified from faecal samples, with approximately half a pellet per rat included in the extraction. Faecal samples were homogenised in 10 ml SM buffer followed by centrifugation twice at 5000 × g at 10°C for 10 min and filtration through a 0.45 µm syringe filter to remove particulates and bacterial cells. NaCl (0.5 M final concentration; Sigma) and 10% w/v polyethylene glycol (PEG-8000; Sigma) were added to the resulting filtrate and incubated at 4°C overnight. Following centrifugation at 5000 × g at 4°C for 20 min, the pellet was resuspended in 400 µl SM buffer. An equal volume of chloroform (Fisher) was added and following 30 sec of vortexing the sample was centrifuged at 2500 × g for 5 min at RT. The aqueous top layer was retained and subjected to RNase I (10 U final concentration; Ambion) and DNase (20 U final concentration; TURBO DNA-free™ Kit, Invitrogen) treatment in accordance with the manufacturer's guidelines.
To isolate the nucleic acids, virus-like particles were incubated with 20 μL of 10% SDS and 2 μL of proteinase K (Sigma, 20 mg/mL) for 20 min at 56°C, prior to lysis by the addition of 100 µL of Phage Lysis Buffer (4.5 M guanidine thiocyanate; 45 mM sodium citrate; 250 mM sodium lauroyl sarcosinate; 562.5 mM βmercaptoethanol; pH 7.0) with incubation at 65°C for 10 min. Viral nucleic acid was purified by two treatments with an equal volume of phenol:chloroform:isoamyl alcohol (25:24:1), finally passing it through a QIAGEN Blood and Tissue Purification Kit and eluting samples in 30 µL of AE Buffer. In order to include any RNA viruses in the downstream analysis, 22 µl of the eluent underwent reverse transcriptase (RT) treatment (Superscript IV, Thermo Fisher) as per the manufacturer's instructions, albeit doubling all volumes to increase the volume of the template used. Prior to library preparation, this RT product was sonicated following adjustment of the volume to 52.5 µl with low-EDTA TE buffer. Shearing of unamplified DNA/cDNA mixture was performed on M220 Focused-Ultrasonicator (Covaris) with the following settings: peak power of 50 W, the duty factor of 20%, 200 cycles per burst, total duration of 35 sec. All steps were performed in accordance with the manufacturer's protocol.
In order to concentrate the DNA/cDNA fragments entering the library, it was passed through the Zymo Genomic DNA Clean & Concentrator™ Kit, eluting in 17 µl DNA Elution Buffer supplied in the kit. Library preparation was carried out using Accel-NGS 1 S Plus kit (Swift Biosciences) according to the manufacturer's instructions, with the addition of a final purification step using a 1:1 DNA/AMPure beads ratio to eliminate any small fragment contamination prior to sequencing. A single-indexed pooled library was sequenced using a 2 × 150 nt paired-end sequencing run on an Illumina Novaseq platform at Genewiz, Germany.
Virome data processing. Illumina sequencing reads were downloaded and processed as follows. Read quality was assessed pre-and post-processing using FastQC version 0.11.3 46 . Adaptors and sequenced nucleotides demonstrating low quality were pruned using Trimmomatic version 0.36 in paired-end mode when the quality dropped below a Phred score of 30 for a 4 bp sliding window 47 . While processed reads less than 70 bp were dumped, the surviving paired and unpaired reads were used in all subsequent analyses. The level of bacterial contamination within the virome sequencing data was estimated by aligning reads against a cpn60 database 48 .
Reads were assembled into contigs using metaSPAdes version 3.11.1 49 . All contigs less than 1000 bp were discarded. Within the virome sequencing, viruses were predicted through several different procedures. The assembled sequences were queried against the NCBI viral RefSeq database version 95 (E-value 1E-05) using BLAST version 2.6.0+ 50 . Circular contigs were detected by analysing their sequence termini. All circular contigs in the virome sequencing data were considered viral. Assembled nucleotide sequences were run through the VirSorter pipeline, where category 1 and 2 viruses were examined further. Viral dark matter was detected by a BLASTn search against NCBI's NT database, with assembled sequences from the viral-enriched shotgun data with no representative in the NT database included in downstream analyses. Viral proteins were predicted using Prodigal version 2.6.3 51 with the 'meta' option enabled and Shine-Dalgarno training bypassed. Contig-encoded proteins were queried against the Prokaryotic Viral Orthologous Groups database (pVOGs) using HMMER version 3.1.b2 52 . The following cut-offs were employed to detect sequences rich in viral proteins: contigs <5 kb needed ≥3 pVOG hits; ≥5 and <10 kb, 4 pVOGs; ≥10 and <20 kb, 5 pVOGs; ≥20 and <40 kb, 6 pVOGs; ≥40 and <60 kb, 7 pVOGs; and ≥60 kb, 8 pVOGs. Sequences identified through the different approaches were pooled together and made nonredundant, keeping the larger of two sequences when the BLAST identity and coverage between sequences exceeded 90%.
Subsequent to the detection of putative viruses in the rat faecal virome sequencing data, potential contaminant sequences were removed as follows. The encoded proteins of all putative viruses were queried against a database of ribosomal proteins. All contigs encoding ribosomal proteins were removed from further analysis. Additionally, proteins were queried for similarity to the Pfam plasmid replication proteins: PF01051, PF01446, PF01719, PF04796, PF05732 and PF06970. Contigs containing an HMMER hit with a score ≥15 were discarded.
Virome sequencing reads were aligned to the final viral contigs database using the WGS sequencing data (kindly supplied by Teagasc, Moorepark). While the virome sequencing data was generated on pooled rat faecal samples, the WGS sequencing was conducted on the total microbiome and performed on multiple time points. Reads were aligned using Bowtie2 version 2.3.4.1 in end-to-end mode 53 . The breadth of coverage of reads mapping across each viral contig was calculated using the BEDTools version 2.26.0 coverage function 54 . In order to determine if a viral sequence was truly present, and not the result of multiple reads stacking onto a single conserved element, the following breadth of coverage cut-offs were applied: minimum three sequencing reads; for contigs <5 kb, 50% of the sequence needed to be covered in aligned reads; contigs ≥5 and <20 kb, 30% of the sequence needed to be covered, and for contigs ≥20 kb, 10% of the sequence needs to be covered. If the breadth of coverage filter was not met for a specific contig within a sample, zero reads were recorded. Viral nucleotide sequences were grouped into VCs, a pseudogenera taxonomic rank, using vConTACT2 version 0.9.8 55 . VCs were generated using default settings with the inclusion of known viruses from the Bacterial and Archaeal Viral RefSeq version 85, with ICTV and NCBI taxonomy. The sequencing reads of viruses belonging to specific VCs were aggregated. Viruses that grouped into ambiguous vConTACT2 clusters (i.e. outliers, overlap) were treated as singletons. Viruses not grouping into VCs were nonetheless retained during the analysis. MetaPhlAn2 was used to profile the WGS sequencing data 56 . The putative familial-level taxonomic compositions of viral sequences were predicted using the Demovir script (https:// github.com/feargalr/Demovir) with default settings. Integrases and recombinases were identified amongst viral-encoded proteins using a probabilistic hidden Markov model (HMM) generated from Pfam protein families; PF00239, PF00589, PF02899, PF07508, PF09003, PF09299, PF10136, PF13356, PFF14659, PF16795 and PF18644.
Virome analysis. Analyses were conducted in R version 3.6.1 57 , implemented through R studio. Images were drawn using 'ggplot2' 58 , utilising the 'RColor-Brewer' and 'jcolors' packages for colour schemes 59,60 . Statistical comparisons were performed using the 'ggpubr' package 61 , with two and multiple group comparisons executed through Wilcoxon and Kruskal-Wallis tests, respectively. Boxplots display the median, hinges (25th and 75th percentiles), whiskers (1.5× inter-quartile range from the hinge), and outliers as individual points. αand β-diversities of samples were calculated with the 'vegan' and 'phyloseq' packages, respectively 62,63 . Differences in sample β-diversities displayed in this study were calculated using the Bray-Curtis dissimilarity index with principal coordinate analysis (PCoA) twodimensional ordination. Permutational multivariate analysis of variance (PER-MANOVA) statistical comparisons of PCoA-grouped samples was conducted with the adonis function from the vegan package. Differential abundance changes in viruses and VCs were calculated using the 'DESeq2' package 64 . Reported p values were adjusted using Bonferroni correction. The heatmap annotation bar was generated using the 'ComplexHeatmap' package 65 . Literature searches were conducted in PubMed using the search term: 'rat' AND 'phage' OR 'rat' AND 'virome'.
Reporting Summary. Further information on research design is available in the Nature Research Reporting Summary linked to this article.

Data availability
The supplementary data and scripts required to generate the images and interpret the results of this study are openly available in Figshare at https://doi.org/10.6084/ m9.figshare.14332985. Virome and WGS sequence data reported are available in the Sequence Read Archive (SRA) database under the BioProject accession numbers PRJNA753493 and PRJEB32682, respectively.