ABSTRACT Dystonia is a complex neurological movement disorder characterized by involuntary muscle contractions. Increasing studies implicate the microbiome as a possible key susceptibility factor for neurological disorders, but the relationship between the gut microbiota and dystonia remains poorly explored. Here, the gut microbiota of 57 patients with isolated dystonia and 27 age- and environment-matched healthy controls was analyzed by 16S rRNA gene amplicon sequencing. Further, integrative analysis of the gut microbiome and serum metabolome measured by high-performance liquid chromatography-mass spectrometry was performed. No difference in α-diversity was found, while β-diversity was significantly different, with a more heterogeneous community structure among dystonia patients than among controls. The most significant changes in dystonia highlighted an increase in Clostridiales, including Blautia obeum, Dorea longicatena, and Eubacterium hallii, and a reduction in Bacteroides vulgatus and Bacteroides plebeius. The functional analysis revealed that genes related to tryptophan and purine biosynthesis were more abundant in gut microbiota from patients with dystonia, while genes linked to citrate cycle, vitamin B[6], and glycan metabolism were less abundant. The evaluation of serum metabolites revealed altered levels of l-glutamic acid, taurine, and d-tyrosine, suggesting changes in neurotransmitter metabolism. The most modified metabolites strongly inversely correlated with the abundance of members belonging to the Clostridiales, revealing the effect of the gut microbiota on neurometabolic activity. This study is the first to reveal gut microbial dysbiosis in patients with isolated dystonia and identified potential links between gut microbiota and serum neurotransmitters, providing new insight into the pathogenesis of isolated dystonia. IMPORTANCE Dystonia is the third most common movement disorder after essential tremor and Parkinson’s disease. However, the cause for the majority of cases is not known. This is the first study so far that reveals significant alterations of gut microbiome and correlates the alteration of serum metabolites with gut dysbiosis in patients with isolated dystonia. We demonstrated a general overrepresentation of Clostridiales and underrepresentation of Bacteroidetes in patients with dystonia in comparison with healthy controls. The functional analysis found that genes related to the biosynthesis of tryptophan, which is the precursor of the neurotransmitter serotonin, were more active in isolated dystonia patients. Altered levels of several serum metabolites were found to be associated with microbial changes, such as d-tyrosine, taurine, and glutamate, indicating differences in neurotransmitter metabolism in isolated dystonia. Integrative analysis suggests that neurotransmitter system dysfunction may be a possible pathway by which the gut microbiome participates in the development of dystonia. The gut microbiome changes provide new insight into the pathogenesis of dystonia, suggesting new potential therapeutic directions. KEYWORDS: dystonia, metagenomics, 16S rRNA, gut microbiome, serum metabolome INTRODUCTION Dystonia is a complex, highly variable neurological movement disorder characterized by sustained or intermittent muscle contractions ([39]1). It can be the manifestation of many neurological disorders, either in isolation (isolated dystonia) or in combination with additional symptoms (combined dystonia) ([40]2). Isolated dystonia includes syndromes in which dystonia is the sole phenotypic manifestation, with the exception of tremor ([41]1). Until now, most patients with common adult-onset dystonia tend to arise sporadically with no clear etiology and pathogenesis, though several associated genes (TOR1A, THAP1, GCH1, and KMT2B) have been identified as the cause of childhood- and adolescent-onset dystonia ([42]3), and these mechanisms include abnormalities in transcriptional regulation, striatal dopaminergic signaling and synaptic plasticity, and a loss of inhibition at neuronal circuits. Clinical electrophysiological studies showed that the dystonic central nervous system (CNS) exhibits to some extent deficient inhibitory neurotransmission, and dystonia is more likely a disorder of motor circuits in a particular brain structure. Researchers believe that dystonia results from an abnormality in the basal ganglia circuits or other brain regions that control movement and, possibly, abnormalities in processing a group of chemicals called neurotransmitters, such as dopamine, serotonin (5-HT), gamma-aminobutyrate (GABA), acetylcholine, purine, and norepinephrine, which help cells in the brain communicate with each other ([43]4, [44]5). The gut microbiome contains a complex community of microbes that live within the gastrointestinal tract. The brain-gut-microbiome axis allows the brain to control gut function but also provides an opportunity for the gut microbiome to influence the brain ([45]6). Emerging evidence has implicated the gut microbiome as an important player in neurodevelopment, as well as in brain diseases such as Parkinson's disease (PD) ([46]7), Alzheimer’s disease (AD) ([47]8), and multiple sclerosis ([48]9). Gut microbes communicate with the central nervous system through nervous-, endocrine-, and immune-signaling mechanisms. As being increasingly recognized, gut permeability, bacterial metabolites (such as short-chain fatty acids [SCFA] and lipopolysaccharide [LPS]), neurotransmitters (e.g., GABA, dopamine, and serotonin), and the vagal nerve are perhaps the most important factors in initiating and mediating microbial interactions with the central nervous system of the body. It has been reported that anticlostridial agents were capable of reversing both gastrointestinal and neurological symptoms in one sporadic case of myoclonus-dystonia ([49]10), providing clues to the mechanisms driving dystonia and potentially other neurological conditions with a gut microbial component. However, the relationship between the gut microbiota and isolated dystonia remains poorly explored. Using 16S rRNA gene amplicon and metagenomic sequencing, we investigated gut microbial characteristics in isolated dystonia and found that dystonia patients carry an altered gut microbiota enriched with Eubacterium hallii, Blautia obeum, and Dorea longicatena. Functional analyses also suggested differences in microbial tyrosine and tryptophan biosynthesis, purine metabolism, and sulfur metabolism pathways. We further performed integrative analysis of gut metagenome and nontargeted serum metabolome profiling of patients with significant gut dysbiosis and revealed an association of the altered gut microbiome with serum metabolites (especially neurotransmitters) in the host. RESULTS Demographics and clinical data. Fifty-seven isolated dystonia patients and 27 healthy controls (HCs) were recruited, and an integrative analysis of the gut microbiome and serum metabolome was performed. Clinical details of the study cohort are shown in [50]Table 1. TABLE 1. Demographics of subjects Parameter[51]^a Value for: __________________________________________________________________ P value Dystonia patients (n = 57) HCs (n = 27) No. (%) of males 15 (26.32) 9 (33.33) 0.506 No. (%) of females 42 (73.68) 18 (66.67) BMI (mean ± SD) 22.71 ± 2.04 23.61 ± 1.80 0.341 Age (yr) (mean ± SD) 45.25 ± 13.00 41.1 ± 11.15 0.158 Age at onset (yr) (mean ± SD) 45.05 ± 12.86 Duration (mo) (mean ± SD) 5.32 ± 3.77 UDRS score (mean ± SD) 9.91 ± 14.43 GDS score (mean ± SD) 6.54 ± 4.06 No. (%) with family history (n; %) 0 HAMA score (mean ± SD) 8.51 ± 4.21 6.78 ± 3.11 0.060 HAMD score (mean ± SD) 8.14 ± 4.13 6.52 ± 2.67 0.066 [52]Open in a new tab ^a UDRS, Unified Dystonia Rating Scale; GDS, Global Dystonia Rating Scale; HAMA, Hamilton Anxiety Scale; HAMD, Hamilton Depression Scale. Data overview. In total, 4,397,348 high-quality 16S rRNA gene V4 sequences were obtained, with 52,349 ± 8,462 (mean ± standard deviation [SD]) per sample after demultiplexing and quality control (see [53]Table S1 in the supplemental material). A total of 997 operational taxonomic units (OTUs) (>97% similarity, excluding singletons and low-abundance taxa with a frequency of <3 for samples; range, 204 to 461; median, 306 taxa per sample) were identified, representing 22 taxonomic phyla. The dominant phyla were Firmicutes (66.4%), Bacteroidetes (17.5%), Proteobacteria (9.4%), and Actinobacteria (5.3%). At the genus level, Bacteroides (13%), Faecalibacterium (10%), unidentified Lachnospiraceae (7%), Agathobacter (6%), and Blautia (5%) were the most abundant, of which all but Bacteroides belong to Clostridiales in Firmicutes ([54]Fig. S1). TABLE S1 Detailed data information obtained by 16S rRNA gene amplicon sequencing for the studied dystonia patients and HCs. Download [55]Table S1, XLSX file, 0.6 MB^ (611.3KB, xlsx) . Copyright © 2021 Ma et al. This content is distributed under the terms of the [56]Creative Commons Attribution 4.0 International license. FIG S1 Bacterial composition of fecal samples collected from 57 isolated dystonia patients and 27 healthy controls and the relative abundance of the most abundant phyla and genera (representing 98.6% of the total reads). Download [57]FIG S1, TIF file, 0.6 MB^ (663.5KB, tif) . Copyright © 2021 Ma et al. This content is distributed under the terms of the [58]Creative Commons Attribution 4.0 International license. Gut microbiota revealed by 16S rRNA gene amplicon sequencing. Gut microbial α-diversity was not significantly different as calculated by Shannon and Simpson indices ([59]Fig. 1A), while nonmetric multidimensional scaling (NMDS) analysis demonstrated a significant difference between isolated dystonia patients and HCs ([60]Fig. 1B) (permutational multivariate analysis of variance [PERMANOVA], R^2 = 0.157, P = 0.002). We observed higher β-diversity in gut microbiota of dystonia patients based on analysis of similarity (ANOSIM), indicating a more heterogeneous community structure among dystonia patients than among HCs ([61]Fig. 1C). Linear discriminant analysis effect size (LEfSe) analysis showed that Clostridiales and Bacteroidetes were the key taxa distinguishing isolated dystonia patients from HCs. Dystonia patients were significantly enriched in Blautia, Bifidobacterium, and unidentified Lachnospiraceae while depleted in Bacteroidetes and unidentified Clostridiales in comparison with HCs ([62]Fig. 1D). Further, the relative abundances of the most abundant taxa (≥1%) at different taxonomic levels were further compared ([63]Table S2 and [64]Fig. S2). At the species level, Eubacterium hallii, Blautia obeum, and Dorea longicatena were enriched in dystonia subjects, while Bacteroides vulgatus, Bacteroides plebeius, Fusobacterium varium, and Clostridium clostridioforme were enriched in HCs ([65]Fig. 1E). FIG 1. [66]FIG 1 [67]Open in a new tab Gut microbiome differences between dystonia patients and HCs. (A) Alpha-diversity calculated by Shannon and Simpson indices. Significance was confirmed using Mann-Whitney U tests. (B) Beta-diversity of gut bacterial communities between isolated dystonia patients and HCs using unsupervised ordination (NMDS). Significant differences were confirmed by PERMANOVA (R^2 = 0.157, P = 0.002). Centroids are indicated by crosses. (C) Within-group and between-group dissimilarities analyzed by ANOSIM (R = 0.087, P = 0.028). (D) LEfSe analysis revealing significant taxonomic differences between isolated dystonia patients and HCs. LDA, linear discriminant analysis. (E) Comparison of the relative abundances of the most abundant taxa (≥1%) at the species level in 16S rRNA gene amplicon sequencing. (F) Comparison of the relative abundances of the most abundant taxa (≥1%) at the species level in metagenomic sequencing. TABLE S2 The significant differences in relative abundance of the most abundant taxa (≥1%) at different taxonomic levels were compared between 16S rRNA gene and metagenomic analyses. The red font denotes the bacteria with significantly higher relative abundance in dystonia patients, while the green font denotes the bacteria with significantly lower relative abundance in dystonia patients. The yellow fill color denotes the significant difference found in both 16S rRNA gene and metagenomic analyses; the light blue fill color denotes the significant difference found only in metagenomic analysis; the light orange fill color denotes the significant difference found only in 16S rRNA analysis. Download [68]Table S2, XLSX file, 0.01 MB^ (12.2KB, xlsx) . Copyright © 2021 Ma et al. This content is distributed under the terms of the [69]Creative Commons Attribution 4.0 International license. FIG S2 Average relative abundance comparisons at different taxonomic levels between the isolated dystonia patients and healthy controls revealed by 16S rRNA gene amplicon sequencing. P values were determined using the Mann-Whitney U test, and the Benjamini-Hochberg procedure (false discovery rate method correction) was applied to obtain adjusted P values for multiple comparisons between groups. Significant differences in taxa (adjusted P value < 0.05) are indicated by an asterisk. Boxes indicate 5th to 95th percentiles, with median relative abundances marked as lines and whiskers indicating the range (minimum/maximum) multiplied by the interquartile range (5th to 95th percentiles) from the boxes. Bacterial taxa are ranked by average relative abundances of the overall gut microbiome. Download [70]FIG S2, TIF file, 0.5 MB^ (512.4KB, tif) . Copyright © 2021 Ma et al. This content is distributed under the terms of the [71]Creative Commons Attribution 4.0 International license. Stratification of gut bacterial community in patients with dystonia. Given more microbial community heterogeneity, we applied Dirichlet multinomial mixture (DMM) analysis to our cohort to investigate possible gut microbial subgroups in dystonia and identified two significantly distinct microbial community states (MCSs) (n = 8 and n = 49; R^2 = 0.138; P < 0.001) using a Laplace approximation ([72]Fig. 2A). Specific cooccurring bacterial families were characteristically enriched in these groups; the MCS1 (n = 8) microbiota was characteristically enriched by Bacteroidaceae, Ruminococcaceae, and Lachnospiraceae. The second larger group, MCS2, which included MCS2A (n = 20) and MCS2B (n = 29), exhibited a sharply decreased abundance of Bacteroidaceae, especially in MCS2B ([73]Fig. 2B). NMDS analysis confirmed a strong and significant relationship among MCS classes and bacterial β-diversity (PERMANOVA, R^2 = 0.29, P < 0.001), corroborating the existence of compositionally distinct microbial states ([74]Fig. 2C). These distinct microbial states exhibited significant differences in diversity (Shannon, Kruskal-Wallis tests, P = 0.002) ([75]Fig. 2D), with MCS2B exhibiting the lowest index. We also investigated whether specific factors (illness duration, age at onset, body mass index [BMI], and severity evaluation scores) were differentially associated with microbial community states. However, no significant correlation was detected. FIG 2. [76]FIG 2 [77]Open in a new tab Two compositionally distinct gut microbial states exist in isolated dystonia patients. (A) Based on Laplace approximation, Dirichlet multinomial mixture (DMM) analysis identified two compositionally distinct bacterial communities (n = 8 and n = 49) in the gut of patients with isolated dystonia. (B) Mean community composition of each microbial state at the family level. (C) NMDS analysis illustrates that DMM-defined gut bacterial communities are compositionally distinct (PERMANOVA, R^2 = 0.29, P = 0.001). (D) Shannon diversity differs significantly across microbial states (ANOVA, P < 0.001). Functional analysis of isolated dystonia by metagenomic sequencing. As shown above, 29 of 57 (51%) of the isolated dystonia patients underwent a gut microbial dysbiosis. To further explore the microbial functional feature of dystonia, we performed metagenomic sequencing of a subgroup consisting of 13 patients from MCS2B and 13 gender- and age-matched HCs. An average of 10.7 ± 2.1 Gb of clean reads was generated per sample, except for one HC sample that was removed due to low depth of sequencing ([78]Table S3). Consistent with the 16S rRNA gene analysis, B. obeum, D. longicatena, and E. hallii were significantly increased in dystonia patients, while B. vulgatus and B. plebeius were decreased ([79]Fig. 1F and [80]Fig. S3). In addition, the genera Alistipes and Parabacteroides were also found to be decreased in patients with dystonia in comparison to HCs ([81]Table S2). TABLE S3 Detailed data information obtained by shotgun metagenomic sequencing for the studied dystonia patients and HCs. Download [82]Table S3, XLSX file, 0.01 MB^ (14.1KB, xlsx) . Copyright © 2021 Ma et al. This content is distributed under the terms of the [83]Creative Commons Attribution 4.0 International license. FIG S3 Average relative abundance comparisons at different taxonomic levels between the isolated dystonia patients and healthy controls revealed by shotgun metagenomic sequencing. P values were determined using the Mann-Whitney U test, and the Benjamini-Hochberg procedure (false discovery rate method correction) was applied to obtain adjusted P values for multiple comparisons between groups. Significant differences in taxa (adjusted P value < 0.05) are indicated by an asterisk. Boxes indicate 5th to 95th percentiles, with median relative abundances marked as lines and whiskers indicating the range (minimum/maximum) multiplied by the interquartile range (5th to 95th percentiles) from the boxes. Bacterial taxa are ranked by average relative abundances of the overall gut microbiome. Download [84]FIG S3, TIF file, 0.5 MB^ (483.6KB, tif) . Copyright © 2021 Ma et al. This content is distributed under the terms of the [85]Creative Commons Attribution 4.0 International license. We further estimated the abundance of metabolic pathways using metagenomic reads mapped to functional orthologs from the KEGG databases to explore differences in the metabolic potential of gut microbiomes between dystonia patients and HCs. Dystonia communities were functionally different from healthy communities and less closely clustered together among individuals, suggesting that interindividual functional variation was higher in dystonia patients than in HCs ([86]Fig. 3A). A total of 22 metabolic pathways (KEGG level 3) were found to be significantly different (P < 0.05, false discovery rate [q] < 0.02) in abundance (>0.01%) between dystonia patients and HCs, including those involved in nucleotide, amino acid, carbohydrate, and lipid metabolism ([87]Table S4). Interestingly, we identified several different pathways, including phenylalanine, tyrosine, and tryptophan biosynthesis (ko00350, ko00400), purine metabolism (ko00230), and sulfur metabolism (ko00920), which were mainly associated with neurotransmitter’s metabolism, that appeared to be more active in the microbiome of dystonia patients ([88]Fig. 3B). In contrast, genes related to the tricarboxylic acid (TCA) cycle (ko00020, ko00630, ko00720), glycan biosynthesis and metabolism pathways (ko00511, ko00531, ko00540, ko00600, ko00603), and vitamin B[6] metabolism (ko00750) were less abundant in dystonia patients ([89]Fig. 3C). In addition, we identified an increased gene abundance for acetyl coenzyme A (acetyl-CoA) (an important component in the biogenic synthesis of the neurotransmitter acetylcholine) biosynthesis in dystonia patients compared to controls (system: trans-cinnamate degradation, trans-cinnamate → acetyl-CoA; phenylacetate degradation, phenylacetate → acetyl-CoA/succinyl-CoA) ([90]Table S4, KEGG module number M0000360, P = 0.07, q = 0.04). We also traced the contributing genes and determined their likely taxonomic origin to determine which bacteria are involved in these pathways ([91]Fig. 3B and [92]C). FIG 3. [93]FIG 3 [94]Open in a new tab Functional characterization of genes related to dystonia. (A) Functional characterization of the dystonia microbiome. (B and C) NMDS plot of Bray-Curtis resemblance generated from square root-transformed KEGG pathway (level 3) relative abundances. Functional differences in dystonia patients are based on selected metabolic pathways via genus-level composition of 16 modules increased (B) or decreased (C) in isolated dystonia patients compared with HCs. TABLE S4 Differences in metabolic potential of gut microbiome between dystonia patients and HCs determined by estimating the abundance of metabolic pathways using our metagenomic reads mapped to functional orthologs from the KEGG databases. Download [95]Table S4, XLSX file, 0.03 MB^ (29.7KB, xlsx) . Copyright © 2021 Ma et al. This content is distributed under the terms of the [96]Creative Commons Attribution 4.0 International license. Associations between gut microbial species and serum metabolites. We performed nontargeted metabolomics profiling of serum samples from the 13 dystonia patients and the 13 HCs in metagenomic analysis and identified 1,543 features defined by retention time and mass/charge ratio, with 224 identified in both positive and negative ionization modes ([97]Table S5). Patients with isolated dystonia showed pronounced metabolic alterations compared with HCs. A supervised partial least-squares discriminant analysis (PLS-DA) using two components (R^2X = 0.525, R^2Y[cum] = 0.99, Q[cum]^2 = 0.975, P < 0.001) was performed, resulting in some separation tendencies between dystonia cases and HCs ([98]Fig. 4A). Based on the PLS-DA models of metabolite profiling data, 242 metabolites were found to be significantly different in abundance between dystonia cases and HCs, with 82 metabolites having a higher concentration and 160 metabolites having a lower concentration in patients with dystonia (P < 0.05; variable importance in projection [VIP] value of >1 and fold change [FC] of >2 or <0.5) ([99]Fig. 4B and [100]Table S5). Among these, 82 metabolites had KEGG identifiers (IDs) in the KEGG database. The dystonia patient serum samples clearly showed the most dramatically decreased levels of several metabolites, including taurine and l-glutamic acid, the most important neurotransmitters, vitamin C, l-serine, ginkgoic acid, aminomethylphosphonic acid (AMPA), and hypoxanthine. In contrast, l-aspartic acid, phosphoric acid, sulfoserine, (S)-anti-citrullinated protein antibodies (ACPA), d-tyrosine, and caffeine levels were significantly higher. The VIP scores of these top metabolites were an indication of their contribution to the PLS model ([101]Table S5). FIG 4. [102]FIG 4 [103]Open in a new tab Associations of gut microbial species with circulating metabolites. (A) Supervised partial least-squares discriminant analysis (PLS-DA) showing differences between dystonia patients and HCs. (B) Volcano plot of metabolic features detected in serum samples after background subtraction and removal of the features found in <30% of the data. Positive log[2] fold change (FC) indicates increased abundance in isolated dystonia patients; +, positive ionization mode; −, negative ionization mode. All P values were adjusted using the Bonferroni method. (C) The top 20 pathways of KEGG pathway enrichment. The x axis shows the ratio of the number of metabolites in the corresponding pathway to the total number of identified metabolites. The greater the value, the higher the degree of metabolite enrichment in the pathway; the color of the dots represents the (log) P value and the size of the dot represents the number of metabolites in the corresponding pathway. The larger number denotes more metabolites in the pathway. The color blocks in the right panel represent different KEGG pathways. The color blocks on the right side were assigned to different KEGG pathways. (D) Significant KEGG pathway enrichment of the significantly changed metabolites. The color blocks in the right panel represent different KEGG pathways. (E) Heatmap of Spearman’s rank correlations between the five bacterial species displaying altered abundance in dystonia patients and 242 different metabolites (only metabolites correlated with at least one species with an adjusted P value of <0.01 are shown; *, P < 0.05; **, P < 0.01; ***, P < 0.001). TABLE S5 Serum metabolites were profiled by HPLC-MS. Download [104]Table S5, XLSX file, 1 MB^ (1,010KB, xlsx) . Copyright © 2021 Ma et al. This content is distributed under the terms of the [105]Creative Commons Attribution 4.0 International license. To further identify pathways affected by dystonia disease, we performed KEGG pathway enrichment analysis on the 242 significantly different metabolites with known KEGG IDs and identified 14 metabolites (VIP > 1.5, P < 0.05) enriched in 9 KEGG pathways: taurine and hypotaurine metabolism, glyoxylate and dicarboxylate metabolism, porphyrin and chlorophyll metabolism, ovarian steroidogenesis, vitamin digestion and absorption, alanine, aspartate, and glutamate metabolism, purine and caffeine metabolism, and biosynthesis of alkaloids derived from histidine and purine ([106]Fig. 4C and [107]D; [108]Table S5). In addition, neurotransmitter-associated metabolites described in the literature, such as d-tyrosine, l-aspartic acid, acetylcarnitine, oleamide, sulfoserine, formiminoglutamic acid, arachidonic acid, and AMPA were also investigated manually. The analyses showed that cell signaling and environmental stress responses were the most relevant biological processes. The most represented metabolic pathways were taurine/hypotaurine, caffeine/purine, and cofactor/vitamin metabolism. Further, Spearman’s correlation coefficients were computed for relationships between the relative abundance of the five dystonia-associated species identified and the different 242 normalized individual metabolomic features. We identified dystonia-associated gut microbial species linked to changes in serum metabolites. Metabolites were grouped into two clusters depending on the correlations, and correlation coefficients with significant P values (<0.01) are shown in [109]Fig. 4E. The abundance of the Clostridiales species enriched in patients with dystonia were positively correlated with the first metabolite cluster, including l-aspartic acid, tyrosine, and sulfoserine, and were otherwise negatively correlated with the second cluster, including l-glutamic acid, taurine, pyruvic acid, and vitamin C ([110]Fig. 4E). Finally, we checked, using the Virtual Metabolic Human database, whether they were made by the species for which we found the most correlations. For l-glutamic acid, taurine, pyruvic acid, l-aspartic acid, tyrosine, and hypoxanthine, we found that the metabolites themselves and their precursors and degradation products had been present in bacteria. No bacterial associations were found for acetylcarnitine, arachidonic acid, and sulfoserine, etc., consequently suggesting that their origin was exclusively human. DISCUSSION Using 16S rRNA gene amplicon and metagenomic sequencing, we investigated gut microbial characteristics in isolated dystonia patients and found that these patients carry an altered gut microbiota enriched with E. hallii, B. obeum, and D. longicatena. Functional analyses also suggested differences in microbial phenylalanine, tyrosine and tryptophan biosynthesis, purine metabolism, and sulfur metabolism pathways. An integrative analysis of the gut microbiome and serum metabolome was employed to gain insight into the links between gut dysbiosis and dystonia, and we found that altered levels of several serum neurotransmitters, especially taurine, l-glutamic acid, and d-tyrosine, were associated with microbial changes, indicating an abnormality in neurotransmitter metabolism in patients with isolated dystonia. Evidence is accumulating for a potential role of the gut microbiota in the pathogenesis of neurological disorders ([111]7[112]–[113]10). One myoclonic dystonia patient was reported to present a 90% improvement in tremor and abnormal posture after taking anticlostridial antibiotics, hinting at a possible correlation between the pathogenesis of dystonia and the gut microbiome underlying the potential role of Clostridial species ([114]10). In this study, intestinal dysbiosis was characterized by significantly increased Clostridiales and decreased Bacteroidetes in dystonia patients, which might explain the neurological improvement in the above-mentioned patient after taking anticlostridial agents, suggesting that Clostridiales may have a potential role in the pathogenesis of dystonia. Blautia and unidentified Lachnospiraceae, belonging to the Lachnospiraceae family, contributed principally to the increase in Clostridiales in dystonia patients. The Lachnospiraceae family is known to participate in the breakdown of carbohydrates into SCFAs ([115]11), which are believed to play a key role in microbiota-gut-brain cross talk ([116]12, [117]13). SCFAs, with acetate, propionate, and butyrate being the most abundant components, may influence the brain by crossing the blood-brain barrier (BBB) ([118]14), modulating neurotransmission, influencing levels of neurotrophic factors, and promoting serotonin (5hydroxytryptamine [5HT]) biosynthesis ([119]13). In a cell model, acetate and butyrate promoted serotonin biosynthesis by inducing the expression of tryptophan 5-hydroxylase 1 (the rate-limiting enzyme for 5HT synthesis) transcription ([120]15). Another study found that butyrate and propionate play a critical role in promoting host 5HT biosynthesis and in regulating levels of 5HT in both colon and serum ([121]16). Furthermore, they also induce tyrosine hydroxylase gene transcription ([122]17), which catalyzes the ratelimiting step in the biosynthesis of dopamine, noradrenaline, and adrenaline ([123]18). A decrease in SCFA production can cause intestinal barrier dysfunction, but whether an increased dose of SCFAs could lead directly to adverse effects on locomotor behavior has not been elucidated and requires animal studies in the future. In the present study, the increased abundance of Clostridiales might lead to elevated SCFA levels and a consequent increase in 5HT and dopamine in the gastrointestinal tract. The abundance of Bacteroides, Parabacteroides, and Alistipes was significantly decreased in isolated dystonia patients in comparison with HCs. A decrease or low level of Bacteroides is believed to be associated with metabolic diseases such as obesity ([124]19) and diabetes ([125]20) and psychiatric disorders such as depression ([126]21). Given that no subjects were associated with obesity or diabetes and that there were no differences in body mass index (BMI) between the two groups, an association between Bacteroides and metabolic factors can be ruled out in our study. Growing evidence suggests that nonmotor symptoms are essential components of dystonia, with psychiatric comorbidity being most prevalent ([127]22), and a link between dystonia and depression has been confirmed ([128]23). In our cohort, patients were newly diagnosed and there was no statistical difference in Hamilton Depression Scale (HAMD) scores between the two groups, so it is unlikely that depression was a confounding factor in our study. Alistipes species are indole positive and may thus influence tryptophan availability ([129]24). Since tryptophan is also the precursor of serotonin, a decreased abundance of Alistipes might therefore disrupt the balance of the intestinal serotonergic system. Until now, the pathophysiological basis of dystonia remains poorly understood. Evidence indicates that the pathophysiological basis of dystonia involves dysfunction of the basal ganglia and cortico-striatal-thalamo-cortical circuits ([130]25), whose activity is modulated among other neurotransmitters by the dopaminergic system ([131]26). In addition to dopamine, a significant neurotransmitter in motor control, serotonin and serotonin-dopamine interactions have received attention ([132]27). Disturbance of the serotonergic system is known to be involved in extrapyramidal movement and psychiatric disorders. Most studies point toward a hypofunction of the serotonergic system in dystonia ([133]27). Analysis of 5-hydroxyindolacetic acid (5-HIAA; the breakdown product of serotonin) in cerebrospinal fluid (CSF) showed significantly reduced levels in patients with idiopathic adult-onset focal dystonia compared with controls ([134]28). A [^11C]DASB [3-amino-4-(2-dimethylaminomethylphenylsulfanyl)-benzonitrile] positron emission tomography study of cervical dystonia patients revealed a significant correlation between serotonin transporter (SERT) binding in the dorsal raphe nucleus and motor and nonmotor symptoms ([135]29). Higher nondisplaceable binding potential (BPND) was statistically correlated with motor symptom severity, pain, and sleep disturbance, with motor symptom severity being the most important predictor of SERT binding ([136]29). Since Lachnospiraceae and Alistipes are involved in the production and regulation of 5-HT biosynthesis, the increase in the Lachnospiraceae family and the decrease in Alistipes species observed in our study may lead to an imbalance in peripheral levels of 5-HT, which might, in turn, modulate brain function by affecting the immune system ([137]30) or signaling to the brain via 5HT receptors on vagal afferent fibers ([138]31). Functional analyses revealed that phenylalanine, tyrosine, and tryptophan biosynthesis tend to be more active in dystonia, as reported previously for PD ([139]32). These aromatic amino acids are important as precursors of 5-HT and catecholamines. Correspondingly, metabolome analysis found that the serum concentration of d-tyrosine, the precursor for dopamine, was significantly elevated in dystonia cases. Composition analysis, functional analysis of the gut microbiome, and serum metabolome analysis all point to dopamine and serotonin metabolic disturbance, suggesting that dopaminergic/serotonergic imbalances may be an important modulator of dystonia. Additionally, the sulfur metabolism pathway appeared to be more active in the microbiome of dystonia patients. In the brain, sulfur metabolism is intimately involved in metabolic cooperation and provides cells with four major reagents critical for methylation reactions (S-adenosylmethionine), antioxidant capacity (glutathione), signaling (H[2]S), and cell volume regulation (taurine). An important intermediate in sulfur metabolism is homocysteine (HCY), a neuronal excitotoxic amino acid, and hyperhomocysteinemia is considered an independent vascular risk factor. Increased HCY has also been observed in movement disorders such as PD ([140]33), Huntington’s disease (HD) ([141]34), and dystonia ([142]35, [143]36). In patients with isolated dystonia, hyperhomocysteinemia was observed independently from immune activation ([144]36). Our study indicates that hyperhomocysteinemia might, to some extent, be related to an increase in sulfur metabolism caused by disturbance of the gut microbiome. It should be noted that the diversity of the gut microbiome is influenced by several factors, including health status, age, diet, and antibiotics ([145]37). One highlight of our study is the enrollment of drug-naive patients. Evidence is accumulating on the effects that drugs can have on the composition and function of the gut microbiome ([146]38). A case-control study on cause-effect relationships could be performed in the future by analyzing drug-naive dystonia patients. Despite being, to our knowledge, the first study to explore the alterations of gut microbiome in dystonia, the present study has some limitations that should be acknowledged. First, as this is a hospital-based, case-control study, selection bias is unavoidable, and the subjects are not fully representative of the general population. However, our hospital has a big center for neurological diseases and appeals to patients across the country. Patients in our cohort came from 12 provinces and are, to some extent, representative. Second, the cross-sectional design did not allow us to determine the occurrence sequence of gut microbial dysbiosis and dystonia nor to elaborate whether functional abnormalities are directly related to the primary pathophysiology or represent compensatory changes. Third, we are aware that a single case-control study may not be sufficient to fully interpret the relationship between gut microbial dysbiosis and isolated dystonia because of the relatively small number of patients involved. Larger numbers of subjects in longitudinal studies or animal experiments are needed to confirm our findings in the future. In summary, this is the first study to reveal a significant gut dysbiosis in isolated dystonia patients, which may be associated with altered serum metabolites. Integrative analysis of the microbiome and metabolome suggest that neurotransmitter system dysfunction may be a possible pathway by which the gut microbiome participates in the development of dystonia. The findings provide new insight into the pathogenesis of dystonia, indicating the potential involvement of the enteric neurotransmitter system in dystonic patients. Further validation is needed in animal experiments and in large populations to explore a clearer role of the gut microbiome in the pathogenesis of dystonia. MATERIALS AND METHODS Participants. Fifty-seven patients with isolated dystonia (15 male, 42 female) and 27 age- and sex-matched healthy participants (9 male, 18 female) were enrolled. Among the patients, 44 had focal dystonia, 11 had segmental dystonia, and two had general dystonia. All patients, from 12 provinces across China, admitted to the Movement Disorders and Botulinum Toxin-A Treatment center underwent a standardized neurological examination by a specialist in movement disorders. All dystonia patients were newly diagnosed and did not receive oral medicine or botulinum toxin injection. None of the participants had a family history of dystonia. All patients were negative for the genetic testing of DYT1 (TOR1A), DYT4 (TUBB4), DYT6 (THAP1), DYT23 (CIZ1), DYT24 (ANO3), DYT25 (GNAL), and COL6A3. Subjects who had no neurological or psychiatric disease were enrolled as healthy controls. Subjects were excluded if they had gastrointestinal diseases, malignant tumors, autoimmune disorders, infectious diseases, or a history of gastrointestinal surgery or had been administered antibiotics or laxatives or immunosuppressive agents in the past 3 months. No subjects were or are on any antianxiety or antidepressant medication. No relevant constipation was present in any of the subjects; dietary and smoking habits did not differ between groups. No detailed dietary plan was requested prior to feces collection, and samples were collected as first bowel movement of the day. Demographic information (gender, age, age of disease onset, disease duration, height, and weight) and severity evaluation scores, by use of the Global Dystonia Rating Scale (GDS) and the Unified Dystonia Rating Scale (UDRS), were collected from all patients. Evaluation for anxiety and depression, by use of the Hamilton Anxiety Scale (HAMA) and the Hamilton Depression Scale (HAMD), was performed for all subjects. Written informed consent was obtained from all individuals. The study was approved by the ethics committees of Beijing Tiantan Hospital. Analysis of fecal microbiota by 16s rRNA gene amplificon sequencing. Stool samples freshly collected from each participant were immediately transported to the laboratory and frozen at −80°C. Bacterial DNA was extracted, and the V4 region of the 16S rRNA gene was amplified by PCR using barcoded universal fusion primers 515 and 806. PCR products were sequenced and used for taxonomic assignment. The number of reads in each sample was standardized to 37,666 for further analysis. Dirichlet multinomial mixtures (DMM) modeling examines taxon frequencies and determines how many microbial community states (MCSs) exist within a data set by using the DirichletMultinomial ([147]39). See [148]Text S1 in the supplemental material for further details on the laboratory protocols. TEXT S1 Supplemental materials and methods. Download [149]Text S1, DOCX file, 0.02 MB^ (20.5KB, docx) . Copyright © 2021 Ma et al. This content is distributed under the terms of the [150]Creative Commons Attribution 4.0 International license. Functional analysis of microbiome by shotgun metagenomic sequencing. Total DNA from the 26 fecal samples (13 cases and 13 controls) were subjected to shotgun metagenomic sequencing. Sheared DNA (350 bp) was used for Illumina library construction, and index codes were added to attribute sequences to each sample. Functional annotation was carried out by BLASTP searches against the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (E value of ≤1e−5 and high-scoring segment pair score of >60). For each functional feature (KO in the KEGG database), we estimated its abundance by accumulating the relative abundance of all genes belonging to that feature. See [151]Text S1 for further details. Untargeted metabolomics analysis of serum samples. All serum samples were thawed on ice, and a quality control (QC) sample, made by blending equal volumes (20 μl) of each serum sample, was used to estimate a mean profile representing all the analytes encountered during analysis. Metabolite extraction was performed by mixing 100 μl thawed serum and 400 μl ice-cold 80% methanol and 0.1% formic acid, vortexing the mixture for 30 s, and keeping it at −20°C for 1 h. After centrifugation at 14,000 × g and 4°C for 20 min, some of supernatant was diluted to a final concentration containing 53% methanol by liquid chromatography-mass spectrometry (LC-MS)-grade water. The samples were subsequently transferred to a fresh Eppendorf tube and then were centrifuged at 14,000 × g at 4°C for 20 min. Finally, the supernatants were subjected to metabolomics profiling by high-performance liquid chromatography (HPLC)-MS in both positive and negative ionization modes. The raw data files were processed using the Compound Discoverer 3.1 (CD3.1; ThermoFisher). Peak intensities were normalized to the total spectral intensity. All data were mean centered and unit variance (UV) scaled before multivariate statistical analysis, and partial least-squares discriminant analysis (PLS-DA) was performed. R^2X, R^2Y[cum], and Q[cum]^2 are diagnostics to assess the statistical significance of the model (“cum” = cumulative). R^2X and R^2Y[cum] represent the interpretation rate of the model for X matrix and Y matrix, respectively. The larger the Q[cum]^2, the better the prediction performance of the model. Significantly different metabolites were determined by a combination of a VIP value of >1 from the PLS-DA model and P values of <0.05 from two-tailed Student’s t tests. See [152]Text S1 for further details. Spearman’s multiomics correlation analysis. Spearman’s correlation coefficients were computed for relationships between the relative abundance of identified dystonia-associated species and normalized individual metabolomic features. A scaled heatmap was constructed for the correlation matrix, including cladogram classification of variables, using the default clustering method. Visual presentation of multiple omics correlations was performed using the Corrplot package in R (v3.5.3). Significant microbiota-metabolite correlations were determined based on an r value of less than −0.3 or greater than 0.3 and a false discovery rate (FDR)-adjusted P value of <0.05. Data availability. The sequencing data set has been deposited in the Genome Sequence Archive in the BIG Data Center, Beijing Institute of Genomics, Chinese Academy of Sciences, under BioProject accession no. [153]PRJCA003075. ACKNOWLEDGMENTS