- 1Florey Institute of Neuroscience and Mental Health, The University of Melbourne, Parkville, VIC, Australia
- 2Department of Anatomy and Neuroscience, The University of Melbourne, Parkville, VIC, Australia
Schizophrenia (SZ) is a psychiatric disorder that constitutes one of the top 10 global causes of disability. More recently, a potential pathogenic role for the gut microbial community (microbiota) has been highlighted, with numerous studies describing dysregulated microbial profiles in SZ patients when compared to healthy controls. However, no animal model of SZ has previously recapitulated the gut dysbiosis observed clinically. Since the metabotropic glutamate receptor 5 (mGlu5) knockout mice provide a preclinical model of SZ with strong face and predictive validity, in the present study we performed gut microbiome profiling of mGlu5 knockout (KO) and wild-type (WT) mice by 16S rRNA sequencing of bacterial genomic DNA from fecal samples, analyzing bacterial diversity and taxonomic composition, as well as gastrointestinal parameters as indicators of gut function. We found a significant genotype difference in microbial beta diversity. Analysis of composition of microbiomes (ANCOM) models were performed to evaluate microbiota compositions, which identified a decreased relative abundance of the Erysipelotrichaceae family and Allobaculum genus in this mouse model of SZ. We also identified a signature of bacteria discriminating between the genotypes (KO and WT), consisting of the Erysipelotrichales, Bacteroidales, and Clostridiales orders and macroscopic gut differences. We thus uncovered global differential community composition in the gut microbiota profile between mGlu5 KO and WT mice, outlining the first evidence for gut dysbiosis in a genetic animal model of SZ. Our findings suggest that this widely used preclinical model of SZ also has substantial utility for investigations of gut dysbiosis and associated signaling via the microbiota–gut–brain axis, as potential modulators of SZ pathogenesis. Our discovery opens up new avenues to explore gut dysbiosis and its proposed links to brain dysfunction in SZ, as well as novel therapeutic approaches to this devastating disorder.
Introduction
Schizophrenia (SZ) is a devastating psychiatric disorder characterized by positive (e.g., hallucinations) and negative (e.g., reduced motivation) symptoms, and cognitive deficits, which existing pharmacological treatments generally fail to comprehensively treat (Leung et al., 2019). While it affects both men and women, there are sex differences in SZ, with men showing an earlier age at onset together with higher propensity to negative symptoms, co-morbid substance abuse and lower social functioning (reviewed in Li et al., 2016). Importantly, chronic gastrointestinal (GI) tract issues exist as comorbid symptoms of SZ, including gut inflammation (Severance et al., 2015; McCutcheon et al., 2020).
More recently a potential role for the microbial community that resides in the gut (called gut microbiota) in SZ pathogenesis has been highlighted, with multiple groups describing dysregulated microbial profiles in SZ patients when compared to healthy controls (Castro-Nallar et al., 2015; Schwarz et al., 2018; Shen et al., 2018; Xu et al., 2019; Szeligowski et al., 2020). Disruption of microbiota can lead to gut microbial imbalance, called gut dysbiosis, and this condition has been shown to occur in various brain disorders, including SZ (Cryan and Dinan, 2012). Gut dysbiosis not only means a differential gut microbial population constituting a pathogenic profile but is also associated with impairment in gut integrity, functionality, intestinal permeability, as well as gut inflammation, and thus an aberrant gut environment. This profile generates a milieu of signaling molecules that ultimately can communicate with the brain through neural communication (via the vagus nerve), endocrine signaling, the immune system and microbial metabolites, thus modulating brain function, and, most remarkably, cognition (reviewed in Gubert et al. (2020)).
Due to discrepancies between clinical studies regarding specific taxa abundance and enrichment, a unique gut microbiome profile for SZ patients is still not determined (Macedo e Cordeiro et al., 2020). However, several clinical studies using different approaches in order to determine the SZ microbiota profiling or signature by determining the population diversity and taxonomic composition or bacterial abundance have been performed, revealing a potential new etiological aspect of SZ associated with an overall gut dysbiosis (reviewed by Rodrigues-Amorim et al., 2018). In fact, recent studies suggest that the disruption of the microbiota–gut–brain axis may promote the development of SZ and additionally, a few bacterial taxa including Veillonellaceae and Lachnospiraceae may associate with SZ severity (Zheng et al., 2019). Further evidence that SZ microbiota may be disease-causing is provided by a metagenome-wide association study showing differences in short-chain fatty acids (SCFAs), and neurotransmitter metabolism, synthesis and degradation in patient microbiota (Zhu et al., 2020). Given the emerging evidence supporting a role for gut dysbiosis in SZ pathogenesis, there is great potential to exploit this as an intervention target, as microbiota can be modulated via a range of therapeutic approaches.
Mice are frequently used to probe mechanistic questions relating to complex brain disorders, including in gut microbiota research (Nguyen et al., 2015). In several studies, fecal matter transplant of gut microbiota from SZ patients has been reported to result in behavioral changes in mice reminiscent of SZ symptoms (Zheng et al., 2019; Zhu et al., 2019, 2020). Interestingly, germ-free mice receiving SZ microbiome fecal transplants demonstrated lower glutamate and higher glutamine and GABA while SZ-relevant behavior in mice (Zheng et al., 2019). However, despite clinical evidence, no animal model of SZ has previously modeled gastrointestinal dysfunction or dysbiosis, although gut dysbiosis has been characterized in animal models of other psychiatric disorders, such as depression (Cheung et al., 2019; Sun et al., 2019).
In the present study, we aimed to profile the gut microbiome of a well characterized mouse model of relevance to SZ, metabotropic glutamate receptor 5 (mGlu5) knockout mice and their wild-type (WT) controls (Gray et al., 2009; Burrows et al., 2015; Zeleznikow-Johnston et al., 2018). Mouse models of schizophrenia are scrutinized against construct, face and predictive validity. While construct validity, or how well an animal model recapitulates the etiology, is difficult to achieve given the complex genetic and environmental drivers of SZ, mGlu5 KO mice model glutamatergic disruption hypothesized to contribute to the disorder (Devon et al., 2001; Gupta et al., 2005) and modulators of the receptor have been of interest in the context of treatment (Balu et al., 2016; Maksymetz et al., 2017; Stansley and Conn, 2018; Uno and Coyle, 2019). The mGlu5 KO mouse model exhibits abnormal brain maturation (Hannan et al., 2001) and a behavioral phenotype that mimics the positive and cognitive symptoms of SZ. Specifically this includes impairments in pre-pulse inhibition (a translational measure of sensorimotor gating in SZ patients and mice), hyperlocomotion (proxy for psychomotor agitation in psychosis) and short-term and long-term spatial memory deficits (Brody and Geyer, 2004; Brody et al., 2004; Gray et al., 2009; Xu et al., 2009; Bird et al., 2014; Burrows et al., 2015; Zeleznikow-Johnston et al., 2018; Lim et al., 2019; Aguilar et al., 2020). Predictive validity refers to response to a clinically effective treatment and beneficial effects on both positive and cognitive symptoms have been shown in mGlu5 mice following treatment with the antipsychotic drug clozapine (Gray et al., 2009).
To date, no study has scrutinized gut health in mGlu5 KO mice or any other mouse model of relevance to SZ, and we hypothesized that differences in the gut microbiome would exist between the two genotypes. To characterize the general gut health, we have performed microbiome profiling using 16S rRNA sequencing, analyzing bacterial diversity and taxonomic composition, as well as gastrointestinal parameters as indicators of gut function. The overall goal of this study was to establish whether this SZ mouse model could be utilized as an appropriate tool for preclinical study of gut dysbiosis.
Materials and Methods
Animal Husbandry
Wild-type (WT) and mGlu5 KO male mice (Grm5tm1Rod; Lu et al., 1997) were generated from heterozygous breeding pairs that had been maintained past generation F10 on a C57Bl/6 background. Genotypes were determined by PCR, from a tail biopsy. The mice were co-housed according to genotype after weaning at 4 weeks of age, in order to avoid potential sharing of microbiota due to the coprophagic nature of mice, unless indicated. The housing condition consisted of open-top standard mouse cages (34 × 16 × 16 cm3; 2–4 mice/box) with basic bedding and nesting materials. Cage changes were performed weekly.
For 16S rRNA sequencing and body weight assessment the same cohort of mice was used, with a total of 6 mGlu5 KO and 6 WT mice, both genotypes split between two cages. All mice had ad libitum access to food and water and were housed in a controlled room of 22°C of temperature and 45% of humidity on a 12:12 h light/dark cycle. All procedures were approved by The Florey Institute of Neuroscience and Mental Health Animal Ethics Committee and were performed in accordance with the relevant guidelines and regulations of the National Health and Medical Research Council Code of Practice for the Use of Animals for Scientific Purposes.
Food and Water Intake
Two extra cohorts were generated and used to assess food and water intake. While for the water intake measures, those mice were still separated by genotype, for food intake measures the respective mice were housed in mixed genotype groups (2–4 mice per box). Cumulative food and water intake were assessed, both over 4 weeks in grouped-housed mice, with the intake normalized to body weight to account for variability. These results are represented as mg of food per g of body weight and as g of water per g of body weight, respectively, to cumulative food and water intake.
Gastrointestinal Measurements
A separate cohort of 7 WT and 6 mGlu5 KO mice at 26 weeks of age was used in order to determine gastrointestinal parameters.
Fecal Output and Fecal Water Content
Fecal output and fecal water content were assessed by single caging the animals for 1 h and counting the number of fecal pellets expelled during this period. All the fecal pellets were collected, and the total weight was measured before being dried at 95 °C for > 3 h. The percentage of difference between the initial total feces weight and the dry weight is taken as fecal water content.
Gastrointestinal Transit Time
Gastrointestinal transit time was determined using the non-absorbable red dye Carmine red (Sigma-Aldrich), prepared as a 6% (w/v) solution in 0.5% methylcellulose (Sigma-Aldrich) filtered and autoclaved prior to use. Non-fasted mice were gavaged with 150 μL of the carmine solution and housed individually. The time taken from gavage to the first appearance of carmine in the feces was recorded as the total transit time.
FITC-Dextran Intestinal Permeability
Intestinal epithelial barrier permeability was assessed by 4 kDa fluorescein isothiocyanate (FITC)−dextran (Sigma-Aldrich). Mice were fasted for 4 h prior to oral gavage with 150 μL FITC-dextran (dissolved in PBS to a concentration of 100 mg/mL). Blood was collected 4 h post administration, and immediately transferred to an EDTA collection tube and further centrifugated at 1,000 × g for 10 min. Plasma was collected and fluorescence was quantified at an excitation wavelength of 485 nm and 528 nm emission wavelength (PHERAstar FSX, Millipore). Serial dilutions of FITC-dextran in PBS were used to calculate a standard curve.
Macroscopic Assessment
Mice were euthanized by cervical dislocation, and their intestines were excised. The intestine was placed in a non-absorbent surface and the length of cecum and colon was measured using a ruler. The weight of cecum with their contents was also evaluated. All of these parameters were normalized to the body weight (g) of the animals, unless indicated.
Fecal Sample Collection and DNA Extraction for 16S rRNA Sequencing
At 12 weeks of age, fecal samples were collected. Mice were individually placed in clean cages for up to 10 min. Fresh pellets were collected and immediately frozen in dry ice and stored at −80°C until further processing.
The fecal genomic DNA was extracted using the PowerSoil HTP kit (Qiagen). The number of sequences per sample was determined as the standard quality control metric (Supplementary Table 1), which showed a mean of approximately 75,500 sequences across all the samples. The extracted genomic DNA and amplified using universal prokaryotic 515F and 806R primers targeting the V4 hypervariable region of the 16S small subunit rRNA gene (Caporaso et al., 2011). Amplicon 16S rRNA gene sequences were generated using paired-end 150 bp sequencing on the Illumina MiSeq platform.
16S rRNA Sequencing, Bioinformatics, and Statistical Analysis
Illumina MiSeq raw FASTQ data were processed using Qiita for quality control, demultiplexing sequences, trimming, and resolving exact sequence variants (ESVs) or Amplicon Sequence Variants (ASVs) with Deblur (Amir et al., 2017). The representative sequences were mapped onto Silva_132 99% to obtain the taxonomic identity of ASVs (Quast et al., 2012). Downstream analyses were computed using R software version 3.5.2.
We first measured the alpha-diversity, which summarizes the diversity of microbial structure within a sample. Reads were rarefied to 1100 reads (Supplementary Figure 1), which is the global minimum of the number of reads in the sequencing data, to compute several alpha-diversity metrics, including species richness (Observed), Shannon, Inverse Simpson and Fisher metrics, using the “Phyloseq” R package (McMurdie and Holmes, 2013). Species richness (Observed) is the number of operational taxonomic units (OTU) observed in a given sample. Shannon, Inverse Simpson and Fisher diversity metrics are a composite measure of richness (number of OTUs present) and evenness (relative abundance of OTUs). Kruskal–Wallis test was used to compare the species richness and alpha-diversity measurements between the genotypes.
To estimate beta-diversity, which summarizes the diversity between samples, the counts were normalized to their relative abundance by dividing raw counts from a particular sample by the total number of reads in each sample. We applied a relative abundance cut-off of 0.01% on the data (including only the relative abundance > 0.01% of the data detected), resulting in a total of 202 ASVs for the subsequent analysis. Two measures of beta-diversity, Bray-Curtis and unweighted UniFrac distances (Lozupone and Knight, 2005), were calculated and used in principal coordinates analysis (PCoA). The unweighted UniFrac distance accounts for the phylogenetic relationship between the OTUs, whereas the Bray-Curtis distance accounts for the abundance of the OTUs. To determine whether the visually observed differences were statistically significant, Adonis (Permutation multivariate ANOVA—PERMANOVA) from the “vegan” R package was performed with 999 permutations (Anderson, 2001; Dixon, 2003). The R2 value reported by Adonis PERMANOVA indicates the amount of variance, on a scale of 0 to 1, in the data which can be explained by the factors tested.
Before testing for differential abundance in the various taxonomic levels, centered log-ratio (CLR) transformation was applied to account for compositionality in microbiome data. Analysis of composition of microbiomes (ANCOM) was used to identify differential relative abundance of each bacterial family and genus between genotypes (Mandal et al., 2015). ANCOM accounts for the compositional nature of the taxa relative abundances and reside on the analysis of difference in pairwise log-ratio while controlling for false discoveries. We applied ANCOM with FDR correction of 0.05. A high “w score” generated by this test indicates the greater likelihood that the null hypothesis can be rejected, indicating the number of times a parameter is significantly different between groups.
Moreover, sparse PLS discriminant analysis (sPLS-DA) from the “mixOmics” package in R was used to identify a signature of discriminative ASVs associated with genotype (Lê Cao et al., 2016; Rohart et al., 2017). The performance of the model is evaluated using the leave-one-out cross validation method.
Co-occurrence Network Analysis
Given that microbiome data is compositional and transformation from absolute to relative abundances introduces spurious correlations, Sparse Correlations for Compositional data (SparCC) was employed to calculate correlations between ASVs to build the co-occurrence network (Friedman and Alm, 2012). SparCC and calculation of two-sided pseudo p values were performed based on bootstrapping of 100 repetitions. Correlations with p-values of < 0.05 and magnitude of = 0.6 (indicating strong co-abundance relationships) or =−0.6 (indicating strong exclusion relationships) were used to construct a network using qgraph (Epskamp et al., 2012). Clusters within the network were determined using the walktrap algorithm with 4 steps implemented in igraph (Csardi, 2006). The network was then visualized and plotted in Cytoscape 3.8 using a customized layout (Shannon, 2003).
General Statistical Analysis
In addition to the gut microbiome analyses described above, to verify the effects of genotype on the body weight, food intake, and water intake, unpaired t-tests were performed. These analyses and respective graphs were constructed using Graph® Pad Prism 8 software. In all cases, significance level was set to p < 0.05.
Results
Gut Dysbiosis in mGlu5 KO Mice
To characterize the gut microbiome, the genomic DNA from feces collected at 12 weeks of age was extracted for 16S amplicon sequencing. Estimates of alpha and beta diversity, both indicative of differences in microbial communities, were obtained using a bioinformatic approach to cluster genetic data from mGlu5 KO and WT mice samples.
The alpha diversity indices (Observed, Shannon, Inverse Simpson, and Fisher) were calculated and plotted to visualize the difference between the two genotypes. Further statistical testing using Kruskal–Wallis method revealed no significance in any of the alpha diversity indices measured [Figure 1A, Kruskal Wallis test, Observed (H = 2.08, df = 1, p = 0.1488); Shannon (H = 0.1, df = 1, p = 0.7488); Inverse Simpson (H = 0.1, df = 1, p = 0.7488); Fisher (H = 2.08, df = 1, p = 0.1488)].
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/582320/fcell-08-582320-HTML-r1/image_m/fcell-08-582320-g001.jpg)
Figure 1. Differential microbiota diversity profiling in mGlu5 KO mice. (A) Alpha diversity metrics including Observed, Shannon, Inverse Simpson (InvSimpson), and Fisher indices. There was no difference in all four indices between WT and mGlu5 KO mice (Kruskal–Wallis test). Beta diversity principal coordinate analysis of WT and mGlu5 KO data using (B) Bray-Curtis distance and (C) unweighted UniFrac distances. The samples were clustered according to the genotype, and there were significant differences between the two genotypes in both Bray-Curtis and unweighted UniFrac distances. Permutation multivariate ANOVA is (n = 6) for both WT and mGlu5 KO mice groups).
For beta diversity, both Bray-Curtis and unweighted UniFrac distance were calculated and the samples were ordinated based on those distances. When the entire population was examined together on a PCoA plot, using both the Bray-Curtis and unweighted UniFrac distance, samples tended to cluster according to the genotype. For Bray-Curtis PCoA, the stratification of samples according to groups was mostly on the second component (Axis 2), accounting for approximately 17% of variation in the data (Figure 1B; Bray Curtis dissimilarity distance PermanovaR2 = 0.19, p = 0.021). Whereas for unweighted UniFrac PCoA, the stratification of samples according to groups was mostly on the first component (Axis 1), accounting for approximately 32% of variation in the data (Figure 1C; unweighted UniFrac distance PermanovaR2 = 0.23, p = 0.003).
Across all mice, the most abundant Phyla was Bacteroidetes (74.9%), followed by Verrucomicrobia (14.6%), and then the Firmicutes (9.0%), which together made up approximately 98.5% of total abundance (Figure 2A). At Class level, the most abundant were Bacteroidia (74.8%), Verrucomicrobiae (14.5%), Clostridia (5.5%) and Erysipelotrichi (3.5%), which together are 98.3% of total abundance (Figure 2B). In addition, the most abundant Order were Bacteriodales (74.8%), Verrucomicrobiales (14.5%), Clostridiales (5.5%) followed by Erysipelotrichales (2.6%) (Figure 2C). A total of 27 bacterial families were detected in genomic data, and 23 remained after passing through a relative abundance cut-off of 0.01% (Figure 2D). The most abundant bacterial Family in its turn was S24.7 (from Bacteriodales order), comprising 51.5% of the bacterial population, followed by Verrucomicrobiaceae (14.6%), Prevotellaceae (8.5%), Bacteroidaceae (8.5%), and Erysipelotrichaceae (2.6%). Finally, a total of 39 genus were detected and 32 of these passed the cut off relative abundance of 0.01% (Figure 2E). The most abundant genus was Unclassified (from S24.7 family) with 51.4% of the bacterial population, followed by Akkermansia (14.5%), Bacteroides (8.5%), and Prevotella (8.4%).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/582320/fcell-08-582320-HTML-r1/image_m/fcell-08-582320-g002.jpg)
Figure 2. Taxonomic composition distribution histograms of each mGlu5 KO and WT sample. At (A) Phylum level, (B) class level, (C) order level, (D) family level, and (E) genus level. n = 6 for both WT and mGlu5 KO mice groups.
Next, we tested differential abundance of bacteria between mGlu5 KO and WT, using analysis of composition of microbiomes (ANCOM), which identified the family Erysipelotrichaceae (Figure 3A) and the genus Allobaculum (Figure 3B) as having different levels between genotypes. The relative abundance of these bacteria was identified as being decreased in both family (Figure 3C) and genus (Figure 3D) in mGlu5 KO mice when compared to WT littermate controls.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/582320/fcell-08-582320-HTML-r1/image_m/fcell-08-582320-g003.jpg)
Figure 3. Differential microbiota abundance in mGlu5 KO mice. Volcano plot for the analysis of composition of microbiomes (ANCOM) test for (A) family and (B) genus level. Only significant bacterial taxa are labeled and colored in blue. Significant taxa are showing high w-stats. Boxplots comparing the relative abundance of the family and genus whose abundance was found to differ significantly between mGlu5 KO and WT (C) Erysipelotrichaceae and (D) Allobaculum. The plots show mean ± SEM (n = 6 for both WT and mGlu5 KO mice groups).
Signature of Bacteria Discriminating Both Genotypes (KO and WT)
In order to detect specific ASVs which could contribute to the stratification of samples according to their genotype (as seen above on the PCoA), sPLS-DA was performed (Figure 4A). sPLS-DA is a multivariate method performed on the clr-transformed microbiome data to identify microbial drivers discriminating particular phenotype groups. The method identified a signature of bacteria discriminating WT and mGlu5 KO mice, detecting 7 ASVs whose respective contribution is shown on the loading plot (classification error rate = 0.11, Figure 4B). For WT mice, this signature consisted of the Erysipelotrichales, Bacteroidales and Clostridiales order. For mGlu5 KO, this signature consisted of bacteria from the Bacteroidales order.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/582320/fcell-08-582320-HTML-r1/image_m/fcell-08-582320-g004.jpg)
Figure 4. Identification of a bacterial signature discriminating the genotypes with sPLS-DA. (A) Sample plots with 0.95 confidence ellipse plots for WT and mGlu5 KO show a strong discrimination between genotypes. The top 7 resulting bacterial signature contributing to component 1 is displayed in (B). The length of the bars indicates the importance of each amplicon sequence variants (ASVs) (at the Order taxonomic rank) in the signature (from bottom to top: decreasing importance) with color indicating the phenotype group with maximum median abundance (n = 6 for both WT and mGlu5 KO mice groups, classification error rate = 0.11).
SparCC-Derived Co-abundance Network Analysis
To obtain a measure of association between ASVs while incorporating their abundance, we used SparCC correlation coefficients that are robust for analyzing compositional microbiome data. A total of 20,301 associations were found, and 78 associations with a p value of less than 0.05 were observed in the network: 46 (59%) of those were positive (r = −0.6) and 32 (41%) were negative (r =−0.6).
The top 10 keystone taxa (taxa with the most interactions with other taxa) detected in this network are ASV_377, ASV_361, ASV_252, ASV_107, ASV_181, ASV_223, ASV_159, ASV_105, ASV_287 and ASV_99 (none of which are detected by the sPLS-DA). Articulation points are defined as nodes in which when removed with its associated edges, makes the graph disconnected. Some of the articulation points in the network include the ASVs also detected by sPLS-DA: including ASV_557 (genus Ileibacterium), ASV_138 (family Muribaculaceae), ASV_134 (family Muribaculaceae), and ASV_275 (family Muribaculaceae) (Figure 5).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/582320/fcell-08-582320-HTML-r1/image_m/fcell-08-582320-g005.jpg)
Figure 5. SparCC-derived co-abundance network analysis. Nodes belonging to the same clusters calculated using the walktrap algorithm are depicted by their shapes (all nodes with circles belong to cluster (A), all nodes with triangles belong to cluster (B)). Nodes that were detected by sPLS-DA and their first neighbors were labeled. Nodes colored in magenta are amplicon sequence variants (ASVs) detected by sPLS-DA to be increased in KO mice, while the nodes colored in dark green are ASVs detected by sPLS-DA to be decreased in KO mice. ASVs are shown at the Genus level, in the case where Genus is not assigned, the next highest taxonomy level is shown. The edges (lines connecting two circles) are colored according to their correlation: blue for negative correlation, red for positive correlation. Isolated nodes are omitted from the network.
Reduced Body Weight in mGlu5 KO Mice Despite No Changes in Food and Water Intake
mGlu5 KO mice presented a reduction in body weight when compared to WT controls (Figure 6A; t = 5.72, df = 10, p = 0.0002). No differences were found in food (Figure 6B; t = 0.04, df = 11, p = 0.65) or water (Figure 6C; t = 1.21, df = 14, p = 0.2) intake between mGlu5 KO mice and WT controls.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/582320/fcell-08-582320-HTML-r1/image_m/fcell-08-582320-g006.jpg)
Figure 6. Reduced body weight in mGlu5 KO mice despite no changes in food and water intake. (A) Decreased body weight in mGlu5 KO mice compared to WT at 12 weeks of age, unpaired t-test. (B) Cumulative food intake (represented as mg of food per g of body weight) and (C) cumulative water intake (represented as g of water per g of body weight), both over 4 weeks between mGlu5 KO mice and WT controls, unpaired t-test. The plots show mean ± SEM (n = 6–8 for both WT and mGlu5 KO mice groups). ***p < 0.001.
Gastrointestinal Parameters
We did not observe any genotype difference in fecal water content (Figure 7A; t = 1.002, df = 11, p = 0.33), fecal output (Figure 7B; t = 0.91, df = 11, p = 0.37), gastrointestinal transit time (Figure 7C; t = 1.06, df = 11, p = 0.3) or gut permeability (Figure 7D; t = 0.72, df = 11, p = 0.48). On the other hand, we observed macroscopic differences between mGlu5 KO and WT mice (Figure 7E), including an increase in cecum length (Figure 7F; t = 4.02, df = 11, p = 0.002) and weight (Figure 7G); t = 2.52, df = 11, p = 0.02), together with an increase in colon length (Figure 7H; t = 8.12, df = 11, p < 0.0001), all relative to body weight. We did not find any macroscopic difference before normalizing to the body weight (Supplementary Figure 2).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/582320/fcell-08-582320-HTML-r1/image_m/fcell-08-582320-g007.jpg)
Figure 7. Gastrointestinal parameters of mGlu5 KO mice. (A) Fecal water content as a percentage of difference between the initial total feces weight and the dry weight. (B) Fecal output, counting the number of fecal pellets expelled during 1 h. (C) Gastrointestinal transit time, using the non- absorbable red dye Carmine red and recording first appearance in the feces. (D) Gut permeability measured by 4 kDa fluorescein isothiocyanate (FITC)-dextran assay. (E) Macroscopic evaluation of intestine. (F) Cecum length, (G) cecum weight, and (H) colon length normalized by body weight, unpaired t-test. The plots show mean ± SEM (n = 6–7 for both WT and mGlu5 KO mice groups). *p < 0.05, **p < 0.01, ***p < 0.001, and ns p > 0.05.
Discussion
This study provides the first evidence that an animal model of SZ replicates the gut dysbiosis observed in patients. In this study, we performed gut microbiome profiling of the mGlu5 KO and WT mice, by 16S rRNA sequencing of fecal samples. We found a significant genotype difference in beta diversity, indicating that mGlu5 KO and WT mice show differences in microbial communities. Specifically, using ANCOM, we uncovered differences in the relative abundance of Erysipelotrichaceae family and Allobaculum genus. We also identified a signature of bacteria discriminating both genotypes, consisting of the bacterial Order Erysipelotrichales, Bacteroidales, and Clostridiales. This new evidence for gut dysbiosis in mGlu5 KO mice provides unique opportunities to explore new approaches to improve treatment for SZ.
There are interesting overlaps between our findings in mice and clinical findings. Similar to our results, previous microbiome profiles of SZ patients have shown no difference between the most abundant phyla Firmicutes and Bacteroidetes (reviewed by Epskamp et al., 2012). Another parallel between our data and clinical SZ findings is that there is no difference in alpha diversity, but a significant difference in beta diversity between SZ patients and controls (Shen et al., 2018; Nguyen et al., 2019). However, it is important to note that while we found a significantly different beta diversity between mGlu5 KO and WT mice, the averaged distance-based coefficient of determination (R2 of the Bray Curtis and UniFrac results) was approximately 0.2, indicating that only around 20% of the microbiome variation is associated with the genotype. Still, a highly diverse gut microbiome is considered an advantage, being related to a healthy environment and lifestyle, usually associated with low or absent pathogenic bacterial species and also associated with improved cognition (e.g., learning, memory and behavioral flexibility; reviewed by Davidson et al., 2018). Therefore, our finding showing a subtle but significant differential beta diversity in mGlu5 KO mice relative to their WT littermate controls indicates a global differential community composition, and thus gut dysbiosis in this mouse model of SZ.
A highly diverse, as well as a balanced composition, of the microbial community is crucial to the homeostasis of the gut microbiome. Gut bacteria regulate basic processes such as digestion, maintenance of intestinal epithelium integrity and, importantly, commensal bacteria promote the first immune response against pathogenic bacteria (Rinninella et al., 2019). Furthermore, gut bacterial diversity and community composition determine the abundance of SCFAs, end-products of microbial fermentation which are vital for gut health (Campbell et al., 2016). One limitation of our study is that we had a relatively low sequencing depth at the alpha rarefaction curve, which could have influenced the alpha diversity analysis specifically. However, when we conducted an analysis investigating the microbial diversity for each sample at various sampling depths, we observed that in all cases the alpha diversity between the two genotypes was not significantly different (Supplementary Table 2).
When comparing humans and mice, there are extensive similarities in anatomy, physiology and genetics which have allowed numerous inferences about human health and pathologies. In gut microbiota research, mouse models have been increasingly used to study the role of the gut microbiota and its association with diseases (Nguyen et al., 2015). However, there are intrinsic differences between these two mammalian systems that make absolute comparisons challenging, particularly at the species level (Hugenholtz and de Vos, 2018). When care is taken in drawing direct parallels between murine and human gut, with regard to microbiota composition, mouse models provide valuable tools to study microbiota imbalance or dysfunction. Thus, our main goal with this study was to establish whether our mouse model of SZ provided an appropriate tool for preclinical studies of SZ. Our findings in the present study demonstrate that these mGlu5 KO mice may model the gut microbiota differential community composition observed in patients.
Using ANCOM we observed a decrease in the relative abundance of the family Erysipelotrichaceae in mGlu5 KO mice while the order Erysipelotrichales was one of the associated bacteria-signature discriminating mGlu5 KO and WT mice. In contrast, the same bacterial family was increased in adult mice born from a maternal immune activation (MIA) animal model of gestational influenza virus infection (Saunders et al., 2020) and a species from this family, Catenibacterium mitsuokai, was found relatively more abundant in SZ oral samples than in controls (Castro-Nallar et al., 2015). ANCOM analysis also indicated a decrease in the genus Allobaculum in mGlu5 KO mice when compared to WT. This genus is a member of the family Erysipelotrichaceae and is associated with a healthy microbiome, suggested to have beneficial gut effects related to mucus properties and formation (Everard et al., 2014). A decrease in the abundance of Allobaculum has been observed in the APP/PS1 mouse model of Alzheimer’s disease (Harach et al., 2017) and absent in restrained mice as a depression-like animal model, while high in the control group (Koo et al., 2010). Importantly, the treatment with the antipsychotic risperidone in mice was shown to increase the abundance of Allobaculum (Bahr et al., 2015).
The order Bacteriodales was one of our signature bacteria responsible for discriminating our genotypes (mGlu5 KO versus WT mice) and it was previously associated with individuals with ultra-high risk for psychosis, when compared to healthy controls (He et al., 2018). In addition, the genus of the same taxa Bacteroides, was associated with SZ in a recent gene set enrichment analysis (GSEA) study that used published psychiatric disorder GWAS data, as well as GWAS of gut microbiota summary data (Cheng et al., 2020). Clostridiales by its turn, another signature bacteria responsible to discriminate our SZ mouse model from control WT mice in our study, was also already associated with ultra-high risk for psychosis individuals when compared to healthy controls (He et al., 2018). Families from this order were demonstrated to be decreased in SZ patients, suggesting an impairment in the maintenance of gut health (Zheng et al., 2019). Other studies also demonstrated an association between lower levels of the Clostridiales order in SZ patients (Shen et al., 2018). Furthermore, a negative correlation between Clostridium_sensu_stricto_1 and the positive and negative syndrome scale (PANSS) score, which is a medical scale used for measuring symptom severity of patients with SZ, was recently reported, suggesting that changes in intestinal microbiota may modulate the prognosis of the disorder (Pan et al., 2020). We also performed a co-occurrence network analysis using SparCC to identify associations between the microbes. We observed that the main hubs (or keystone taxa) of the network are not microbes detected by sPLS-DA. The microbes identified by sPLS-DA are clustered together and their interactions are mainly local (within their own neighborhood). Some of these sPLS-DA identified microbes are also articulation points, in which they bridge the connection between two separate neighborhoods. However, we acknowledge that this approach may be insufficient in inferring microbial community structure (Hirano and Takemoto, 2019).
We have therefore identified interesting overlaps between our present findings and previous clinical SZ literature, with promising translational bacterial targets. However, it is important to note that, regarding specific taxa abundance and enrichment, the discrepancies between clinical studies are even more prominent, mainly due to methodological differences and inherent limitation of the clinical heterogeneity of SZ (Macedo e Cordeiro et al., 2020). Therefore, a simple gut microbiome signature for SZ patients is still not determined. It is essential in clinical studies to rigorously control for treatment, diet, physical activity and previous prebiotic/probiotic/antibiotic use, among other confounders; however, this is not always possible. Nevertheless, the fully controlled environment (as well as genetic control) that animal studies provide can overcome these limitations, and preclinical studies that replicate a general dysbiosis observed in patients, including this mGlu5 KO mouse model, represent a useful tool for the study of the gut microbiome in psychiatric disorders.
A few studies with animal models of relevance to SZ have previously suggested a promising potential of a microbial community role associated to these preclinical models. Sub-chronic N-methyl-D-aspartate receptor (NMDAR) antagonist phencyclidine (PCP) administration in male Lister Hooded rats resulted in cognitive impairment and differential beta diversity (Pyndt Jørgensen et al., 2015). A correlation between gut microbiota and memory performance was also demonstrated, in addition to the administration of antibiotic being able to reverse the cognitive deficits observed in sub-chronic PCP treated rats (Pyndt Jørgensen et al., 2015). However, this last result should be carefully interpreted since they did not investigate the direct central effect that ampicillin is able to perform, and thus this antibiotic may have acted directly on the brain rather than via the gut (Nau et al., 2010). Similarly, in an MIA model of maternal influenza viral infection and associated neurodevelopmental disorder, cognitive impairments were observed in the offspring and antibiotic administration during their prepuberal period was able to prevent that outcome (Saunders et al., 2020). Additionally, the same study showed that adult MIA offspring displayed altered gut microbiota, with differential relative abundance between components of the gut microbiota, including Ruminococcaceae. Likewise, in Wistar rats, an animal model combining MIA with polyI:C and adolescent cannabinoid exposure demonstrated a sex-specific effect, with adolescent female offspring exhibiting decreased fecal levels of Bifidobacterium longum, measured using qPCR (Katz-Barber et al., 2020). While very promising, these aforementioned studies should be interpreted with caution, since recapitulating the complexity of SZ in an animal model is still a challenge, and each study has its own limitations. In particular, the PCP model lacks construct validity (Sams-Dodd, 1999) and MIA animal models are associated with a broad range of disorders, including ASD, and thus MIA may be considered to be a non-specific primer of multiple neurodevelopmental disorders (Estes and McAllister, 2016).
The genetics of SZ is complex and heterogeneous and thus genome-wide association studies (GWAS) do not identify single mutations strongly associated with this psychiatric disorder. However, mGlu5 has been implicated in SZ pathogenesis (Devon et al., 2001; Gupta et al., 2005), along with many related proteins associated with synapse development and function. Decreased mGlu5 mRNA and protein levels in SZ patients have been reported in postmortem studies (Volk et al., 2010; Matosin et al., 2015). The mGlu5 KO mice demonstrate impaired prepulse inhibition (PPI) (Brody and Geyer, 2004; Burrows et al., 2015), baseline hyperactivity and a hypersensitivity to MK-801-induced hyperlocomotion (Lipina et al., 2007; Gray et al., 2009), impaired circadian process (Aguilar et al., 2020), clinically relevant cognitive deficits on touchscreen tasks (Zeleznikow-Johnston et al., 2018) and thus excellent face validity as a preclinical model of SZ. Lastly, chronic treatment with clozapine is able to reverse SZ-related behaviors (Gray et al., 2009), demonstrating strong predictive validity. Taken together, these studies indicate that mGlu5 KO mouse model has great utility for the study of SZ, including the altered microbiome associated with gut dysbiosis.
Recently, a study focused on the association between depressive behavior and the mGlu5 receptor reported 16S rRNA sequencing in feces from another line of mGlu5 KO mice, concluding that there were no gut microbiota changes between genotypes (Cai et al., 2020). Discrepancies between this study and ours are likely to result from multiple factors, including differences in the generation and genetic background of mGlu5 KO mice. In their study, Cai and colleagues generated mGlu5 KO mice by crossing Grm5flox/flox mice with B6.C-Tg (CMV-cre) mice but did not explore SZ-like endophenotypes. Our Grm5tm1Rod line of mice have a constitutive null mutation in the mGlu5 gene, with strong face and predictive validity as a preclinical model of SZ (Hannan et al., 2001; Gray et al., 2009; Burrows et al., 2015). In addition, Cai and collaborators have not extensively analyzed their microbiome data. In particular, similar to their study, we didn’t observe any difference between alpha diversity, but when we analyzed other levels of complexity, such as beta diversity and sPLS-DA (low classification error rate), we discovered a genotype effect. Therefore, following thorough and exhaustive microbiome profiling, our present findings are the first to demonstrate gut dysbiosis in the mGlu5 KO mouse model of SZ.
Consistent with previous findings (Bradbury et al., 2005), we also observed a reduced body weight in mGlu5 KO mice when compared to WT controls, despite no changes in food and water intake between groups. In fact, mGlu5 has been suggested to be a mediator of energy balance in rodents by decreasing caloric efficiency, suggesting increased energy expenditure in mGlu5 KO mice (Bradbury et al., 2005). Corroborating these findings, it was also demonstrated that activation of lateral hypothalamic mGlu5 receptors elicits feeding in rats (Ploj et al., 2010; Charles et al., 2014) and mGlu5 KO mice are resistant to diet-induced obesity (Bradbury et al., 2005). Interestingly, the abundance of the genus Allobaculum, which we have found decreased in mGlu5 KO mice, has been correlated with body weight and dietary-induced inflammation markers, including leptin and IL-22 (Ravussin et al., 2012; Zenewicz et al., 2013; Everard et al., 2014). Considering the highly explored and relevant role of gut microbiota composition in the pathogenesis of obesity (Leocádio et al., 2020), our novel findings showing a differential gut microbiome profile between mGlu5 KO and WT mice may have relevance to the study of metabolic and eating disorders, including obesity.
Glutamate is a central and peripheric modulator, with glutamatergic dysfunction associated with both central nervous system (CNS) disorders and GI diseases, which in many cases display intercorrelated co-morbidities (e. g. inflammatory bowel disease (IBD) and depression) (Keefer and Kane, 2017; Baj et al., 2019). Not surprisingly, SZ is also associated with GI issues, such as gut inflammation and gut cellular barrier defects (Severance et al., 2015). In fact, the glutamatergic system is directly implicated in GI modulation, with the mGlu5 receptor being also expressed in the GI tract (Ferrigno et al., 2017) and involved with the peripheral excitatory modulation of vagal afferent mechanosensitivity (Slattery et al., 2006). Interestingly, antipsychotics, the main treatment for SZ, induce weight gain, constipation and metabolic syndrome (Kanji et al., 2018; Zeng et al., 2020). It is possible that the negative side effects of antipsychotics may be at least partly related to an antibiotic-like side effect on gut microbiota, by altering its composition and decreasing diversity (reviewed in (Dinan and Cryan, 2018; Skonieczna-Źydecka et al., 2019)). Thus, microbiota interventions adjunctively with antipsychotics are promising to alleviate adverse side effects and mGlu5 KO mice may provide a useful preclinical tool to investigate these and other aspects of SZ pathogenesis and treatment.
Considering that the interactions between microbiota composition and gut status are bidirectional, with the bacteria population playing an important role in gut function and the gut function affecting the diversity of the microbiome in the GI tract (Barbara et al., 2005), we evaluated GI health parameters. We didn’t see any difference in fecal water content, fecal output or gastrointestinal transit time, which have been interrelated with gut microbiota composition (Kashyap et al., 2013; Vandeputte et al., 2016). Similarly, we didn’t see a difference in gut permeability, which has been linked to dysbiosis and increased inflammation (Chakaroun et al., 2020). One possible explanation is that, due to the fact that only around 20% of the microbiome variation revealed by beta diversity is associated to the genotype, dysbiosis can potentially not be strong enough to modulate or be reflected in the gastrointestinal parameters analyzed. On the other hand, after normalization by body weight we observed macroscopic differences between mGlu5 KO and WT mice, including an increase in cecum length and weight, together with an increased colon length. Since colon length can be considered a marker of colonic inflammation (Chassaing et al., 2014) and cecum measurements also associates with gut microbiota composition (Ge et al., 2017) further study is needed to explore whether mGlu5 KO mice present gastrointestinal dysfunction associated with inflammation. Therefore, clarifying the extent of the gut dysbiosis in this mGlu5 KO mouse model and establishing whether gut dysbiosis is causally associated with GI dysfunction, would be of great interest.
Conclusion
In conclusion, the microbiome phenotype observed in mGlu5 KO mice is to some extent in line with reports of gut dysbiosis in SZ patients, and thus this animal model provides a novel tool to explore the mechanistic understanding of how and when dysbiosis arises in SZ. Furthermore, our present findings can be used to explore the link between dysbiosis and SZ symptoms (i.e., behavioral impairments) as well as the potential utility of gut microbiome restoration as a therapeutic approach (e.g., fecal matter transplantation, probiotics and prebiotics), while clarifying the extent of microbiota–gut–brain axis dysfunction in this model. Moreover, considering that mGlu5 KO mice also present deficits in basic discrimination learning and cognitive flexibility (Zeleznikow-Johnston et al., 2018), this mouse model can be also used as a tool for the investigation of the role of the gut microbiome in generalized cognitive impairments. Thus, our present findings may also inform other cognitive disorders where mGlu5 receptors play a role and the microbiome association is still not well understood, such as addiction (Terbeck et al., 2015; Meckel and Kiraly, 2019).
Our results provide the first evidence that a genetic animal model of SZ, exhibiting both face and predictive validity, at least partly replicates the gut dysbiosis observed in patients. These new findings of gut dysbiosis in the mGlu5 KO mice provide opportunities to explore novel approaches focusing on the microbiota–gut–brain axis in SZ. We therefore propose that this preclinical model of SZ can be used as a tool to investigate how gut dysbiosis may contribute to SZ pathogenesis, via the microbiota–gut–brain axis. Our new findings may also inform the development of novel therapeutic approaches for SZ, a devastating disorder for which new treatments are urgently needed.
Data Availability Statement
The datasets and metadata related to this study have been deposited in the NCBI Sequence Read Archive under BioProject number PRJNA659149. Furthermore, the R code for the analysis has been uploaded to the following site: https://github.com/gkong1/mGLU5KO_microbiome_analysis.
Ethics Statement
The animal study was reviewed and approved by the Florey Institute of Neuroscience and Mental Health Animal Ethics Committee.
Author Contributions
CG was involved in the experimental design, data collection, data analysis, and manuscript writing. GK was involved in microbiome data analysis. VU, AZ-J, EB, and TR were involved in the data collection, analysis, and editing of the manuscript. AH was involved in the experimental design, project management and funding, data analysis, and drafting of the manuscript. All authors contributed to the article and approved the submitted version.
Funding
CG is a recipient of Melbourne Early Career Researcher Award. EB is an NHMRC-ARC Dementia Research Development Fellow. TR is an NHMRC Dementia Fellow. AH is an NHMRC Principal Research Fellow and is also supported by NHMRC Project Grants, an ARC Discovery Project and the DHB Foundation, Equity Trustees.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We thank Jamie Liew for helping with sample collection and Shanshan Li for laboratory management.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2020.582320/full#supplementary-material
References
Aguilar, D. D., Strecker, R. E., Basheer, R., and McNally, J. M. (2020). Alterations in sleep, sleep spindle, and EEG power in mGluR5 knockout mice. J. Neurophysiol. 123, 22–33. doi: 10.1152/jn.00532.2019
Amir, A., McDonald, D., Navas-Molina, J. A., Kopylova, E., Morton, J. T., Zech Xu, Z., et al. (2017). Deblur rapidly resolves single-nucleotide community sequence patterns. mSystems 2:e00191-16.
Anderson, M. J. (2001). A new method for non-parametric multivariate analysis of variance: non-Parametric Manova For Ecology. Austral Ecol. 26, 32–46. doi: 10.1111/j.1442-9993.2001.01070.pp.x
Bahr, S. M., Weidemann, B. J., Castro, A. N., Walsh, J. W., deLeon, O., Burnett, C. M. L., et al. (2015). Risperidone-induced weight gain is mediated through shifts in the gut microbiome and suppression of energy expenditure. EBioMedicine 2, 1725–1734. doi: 10.1016/j.ebiom.2015.10.018
Baj, A., Moro, E., Bistoletti, M., Orlandi, V., Crema, F., and Giaroni, C. (2019). Glutamatergic signaling along the microbiota-gut-brain Axis. Int. J. Mol. Sci. 20:1482. doi: 10.3390/ijms20061482
Balu, D. T., Li, Y., Takagi, S., Presti, K. T., Ramikie, T. S., Rook, J. M., et al. (2016). An mGlu5-positive allosteric modulator rescues the neuroplasticity deficits in a genetic model of NMDA receptor hypofunction in schizophrenia. Neuropsychopharmacology 41, 2052–2061. doi: 10.1038/npp.2016.2
Barbara, G., Stanghellini, V., Brandi, G., Cremon, C., Nardo, G. D., De Giorgio, R., et al. (2005). Interactions between commensal bacteria and gut sensorimotor function in health and disease. Am. J. Gastroenterol. 100, 2560–2568. doi: 10.1111/j.1572-0241.2005.00230.x
Bird, M. K., Lohmann, P., West, B., Brown, R. M., Kirchhoff, J., Raymond, C. R., et al. (2014). The mGlu5 receptor regulates extinction of cocaine-driven behaviours. Drug Alcohol Depend. 137, 83–89. doi: 10.1016/j.drugalcdep.2014.01.017
Bradbury, M. J., Campbell, U., Giracello, D., Chapman, D., King, C., Tehrani, L., et al. (2005). Metabotropic glutamate receptor mGlu5 is a mediator of appetite and energy balance in rats and mice. J. Pharmacol. Exp. Ther. 313, 395–402. doi: 10.1124/jpet.104.076406
Brody, S. A., Dulawa, S. C., Conquet, F., and Geyer, M. A. (2004). Assessment of a prepulse inhibition deficit in a mutant mouse lacking mGlu5 receptors. Mol. Psychiatry 9, 35–41. doi: 10.1038/sj.mp.4001404
Brody, S. A., and Geyer, M. A. (2004). Interactions of the mGluR5 gene with breeding and maternal factors on startle and prepulse inhibition in mice. Neurotox. Res. 6, 79–90. doi: 10.1007/BF03033300
Burrows, E. L., McOmish, C. E., Buret, L. S., Van den Buuse, M., and Hannan, A. J. (2015). Environmental enrichment ameliorates behavioral impairments modeling schizophrenia in mice lacking metabotropic glutamate receptor 5. Neuropsychopharmacology 40, 1947–1956. doi: 10.1038/npp.2015.44
Cai, G., Zhu, Y., Chen, J., Zhao, S., Wang, L., Wang, M., et al. (2020). Analysis of the gut microbiota and inflammatory factors in mGluR5-Knockout Mice. Front. Psychiatry 11:335. doi: 10.3389/fpsyt.2020.00335
Campbell, S. C., Wisniewski, P. J., Noji, M., McGuinness, L. R., Häggblom, M. M., Lightfoot, S. A., et al. (2016). The effect of diet and exercise on intestinal integrity and microbial diversity in mice. PLoS One 11:e0150502. doi: 10.1371/journal.pone.0150502
Caporaso, J. G., Lauber, C. L., Walters, W. A., Berg-Lyons, D., Lozupone, C. A., Turnbaugh, P. J., et al. (2011). “Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample,” in Proceedings of the National Academy of Sciences 108, 4516–4522. doi: 10.1073/pnas.1000080107
Castro-Nallar, E., Bendall, M. L., Pérez-Losada, M., Sabuncyan, S., Severance, E. G., Dickerson, F. B., et al. (2015). Composition, taxonomy and functional diversity of the oropharynx microbiome in individuals with schizophrenia and controls. PeerJ 3:e1140. doi: 10.7717/peerj.1140
Chakaroun, R. M., Massier, L., and Kovacs, P. (2020). Gut microbiome, intestinal permeability, and tissue bacteria in metabolic disease: Perpetrators or bystanders? Nutrients 12:1082. doi: 10.3390/nu12041082
Charles, J. R., Duva, M. A., Ramirez, G. J., Lara, R. L., Yang, C. R., and Stanley, B. G. (2014). Activation of lateral hypothalamic mGlu1 and mGlu5 receptors elicits feeding in rats. Neuropharmacology 79, 59–65. doi: 10.1016/j.neuropharm.2013.10.033
Chassaing, B., Aitken, J. D., Malleshappa, M., and Vijay-Kumar, M. (2014). Dextran sulfate sodium (DSS)-induced colitis in mice. Curr. Protoc. Immunol. 104, 15.25.1–15.25.14. doi: 10.1002/0471142735.im1525s104
Cheng, S., Han, B., Ding, M., Wen, Y., Ma, M., Zhang, L., et al. (2020). Identifying psychiatric disorder-associated gut microbiota using microbiota-related gene set enrichment analysis. Brief. Bioinformatics 21, 1016–1022. doi: 10.1093/bib/bbz034
Cheung, S. G., Goldenthal, A. R., Uhlemann, A.-C., Mann, J. J., Miller, J. M., and Sublette, M. E. (2019). Systematic review of gut microbiota and major depression. Front. Psychiatry 10:34. doi: 10.3389/fpsyt.2019.00034
Cryan, J. F., and Dinan, T. G. (2012). Mind-altering microorganisms: the impact of the gut microbiota on brain and behaviour. Nat. Rev. Neurosci. 13, 701–712. doi: 10.1038/nrn3346
Csardi, G. N. T. (2006). The igraph software package for complex network research. InterJ. Complex Syst. 1695, 1–9.
Davidson, G. L., Cooke, A. C., Johnson, C. N., and Quinn, J. L. (2018). The gut microbiome as a driver of individual variation in cognition and functional behaviour. Philos. Trans. R. Soc. B Biol. Sci. 373:20170286. doi: 10.1098/rstb.2017.0286
Devon, R. S., Anderson, S., Teague, P. W., Muir, W. J., Murray, V., Pelosi, A. J., et al. (2001). The genomic organisation of the metabotropic glutamate receptor subtype 5 gene, and its association with schizophrenia. Mol. Psychiatry 6, 311–314. doi: 10.1038/sj.mp.4000848
Dinan, T. G., and Cryan, J. F. (2018). Schizophrenia and the microbiome: time to focus on the impact of antipsychotic treatment on the gut microbiota. World J. Biol. Psychiatry 19, 568–570. doi: 10.1080/15622975.2018.1540793
Dixon, P. (2003). VEGAN, a package of R functions for community ecology. J. Veg. Sci. 14, 927–930. doi: 10.1111/j.1654-1103.2003.tb02228.x
Epskamp, S., Cramer, A. O. J., Waldorp, L. J., Schmittmann, V. D., and Borsboom, D. (2012). qgraph?: network visualizations of relationships in psychometric data. J. Stat. Softw. 48, 1–18. doi: 10.18637/jss.v048.i04
Estes, M. L., and McAllister, A. K. (2016). Maternal immune activation: implications for neuropsychiatric disorders. Science 353, 772–777. doi: 10.1126/science.aag3194
Everard, A., Lazarevic, V., Gaïa, N., Johansson, M., Ståhlman, M., Backhed, F., et al. (2014). Microbiome of prebiotic-treated mice reveals novel targets involved in host response during obesity. ISME J. 8, 2116–2130. doi: 10.1038/ismej.2014.45
Ferrigno, A., Berardo, C., Di Pasqua, L. G., Siciliano, V., Richelmi, P., and Vairetti, M. (2017). Localization and role of metabotropic glutamate receptors subtype 5 in the gastrointestinal tract. World J. Gastroenterol. 23, 4500–4507. doi: 10.3748/wjg.v23.i25.4500
Friedman, J., and Alm, E. J. (2012). Inferring correlation networks from genomic survey data. PLoS Comput. Biol. 8:e1002687. doi: 10.1371/journal.pcbi.1002687
Ge, X., Ding, C., Zhao, W., Xu, L., Tian, H., Gong, J., et al. (2017). Antibiotics-induced depletion of mice microbiota induces changes in host serotonin biosynthesis and intestinal motility. J. Transl. Med. 15:13.
Gray, L., van den Buuse, M., Scarr, E., Dean, B., and Hannan, A. J. (2009). Clozapine reverses schizophrenia-related behaviours in the metabotropic glutamate receptor 5 knockout mouse: association with N-methyl-d-aspartic acid receptor up-regulation. Int. J. Neuropsychopharmacol. 12, 45–60. doi: 10.1017/S1461145708009085
Gubert, C., Kong, G., Renoir, T., and Hannan, A. J. (2020). Exercise, diet and stress as modulators of gut microbiota: implications for neurodegenerative diseases. Neurobiol. Dis. 134:104621. doi: 10.1016/j.nbd.2019.104621
Gupta, D. S., McCullumsmith, R. E., Beneyto, M., Haroutunian, V., Davis, K. L., and Meador-Woodruff, J. H. (2005). Metabotropic glutamate receptor protein expression in the prefrontal cortex and striatum in schizophrenia. Synapse 57, 123–131. doi: 10.1002/syn.20164
Hannan, A. J., Blakemore, C., Katsnelson, A., Vitalis, T., Huber, K. M., Bear, M., et al. (2001). PLC-β1, activated via mGluRs, mediates activity-dependent differentiation in cerebral cortex. Nat. Neurosci. 4, 282–288. doi: 10.1038/85132
Harach, T., Marungruang, N., Duthilleul, N., Cheatham, V., Mc Coy, K. D., Frisoni, G., et al. (2017). Reduction of Abeta amyloid pathology in APPPS1 transgenic mice in the absence of gut microbiota. Sci. Rep. 7:41802. doi: 10.1038/srep41802
He, Y., Kosciolek, T., Tang, J., Zhou, Y., Li, Z., Ma, X., et al. (2018). Gut microbiome and magnetic resonance spectroscopy study of subjects at ultra-high risk for psychosis may support the membrane hypothesis. Eur. Psychiatry 53, 37–45. doi: 10.1016/j.eurpsy.2018.05.011
Hirano, H., and Takemoto, K. (2019). Difficulty in inferring microbial community structure based on co-occurrence network approaches. BMC Bioinformatics 20:329. doi: 10.1186/s12859-019-2915-1
Hugenholtz, F., and de Vos, W. M. (2018). Mouse models for human intestinal microbiota research: a critical evaluation. Cell. Mol. Life Sci. 75, 149–160. doi: 10.1007/s00018-017-2693-8
Kanji, S., Fonseka, T. M., Marshe, V. S., Sriretnakumar, V., Hahn, M. K., and Müller, D. J. (2018). The microbiome-gut-brain axis: implications for schizophrenia and antipsychotic induced weight gain. Eur. Arch. Psychiatry Clin. Neurosci. 268, 3–15. doi: 10.1007/s00406-017-0820-z
Kashyap, P. C., Marcobal, A., Ursell, L. K., Larauche, M., Duboc, H., Earle, K. A., et al. (2013). Complex interactions among diet, gastrointestinal transit, and gut microbiota in humanized Mice. Gastroenterology 144, 967–977. doi: 10.1053/j.gastro.2013.01.047
Katz-Barber, M. W., Hollins, S. L., Cuskelly, A., Leong, A. J. W., Dunn, A., Harms, L., et al. (2020). Investigating the gut-brain axis in a neurodevelopmental rodent model of schizophrenia. Brain Behav. Immun. Health 3:100048. doi: 10.1016/j.bbih.2020.100048
Keefer, L., and Kane, S. V. (2017). Considering the bidirectional pathways between depression and IBD: recommendations for comprehensive IBD care. Gastroenterol. Hepatol. 13, 164–169.
Koo, J. W., Russo, S. J., Ferguson, D., Nestler, E. J., and Duman, R. S. (2010). Nuclear factor- B is a critical mediator of stress-impaired neurogenesis and depressive behavior. Proc. Natl. Acad. Sci. U.S.A. 107, 2669–2674. doi: 10.1073/pnas.0910658107
Lê Cao, K.-A., Costello, M.-E., Lakis, V. A., Bartolo, F., Chua, X.-Y., Brazeilles, R., et al. (2016). MixMC: a multivariate statistical framework to gain insight into microbial communities. PLoS One 11:e0160169. doi: 10.1371/journal.pone.0160169
Leocádio, P. C. L., Oriá, R. B., Crespo-Lopez, M. E., and Alvarez-Leite, J. I. (2020). Obesity: more than an inflammatory, an infectious disease? Front. Immunol. 10:3092. doi: 10.3389/fimmu.2019.03092
Leung, C. C.-Y., Gadelrab, R., Ntephe, C. U., McGuire, P. K., and Demjaha, A. (2019). Clinical course, neurobiology and therapeutic approaches to treatment resistant schizophrenia. Toward an integrated view. Front. Psychiatry 10:601. doi: 10.3389/fpsyt.2019.00601
Li, R., Ma, X., Wang, G., Yang, J., and Wang, C. (2016). Why sex differences in schizophrenia? J. Transl. Neurosci. 1, 37–42.
Lim, J., Kim, E., Noh, H. J., Kang, S., Phillips, B. U., Kim, D. G., et al. (2019). Assessment of mGluR5 KO mice under conditions of low stress using a rodent touchscreen apparatus reveals impaired behavioural flexibility driven by perseverative responses. Mol. Brain 12:37.
Lipina, T., Weiss, K., and Roder, J. (2007). The Ampakine CX546 restores the prepulse inhibition and latent inhibition deficits in mGluR5-Deficient Mice. Neuropsychopharmacology 32, 745–756. doi: 10.1038/sj.npp.1301191
Lozupone, C., and Knight, R. (2005). UniFrac: a new phylogenetic method for comparing microbial communities. Appl. Environ. Microbiol. 71, 8228–8235. doi: 10.1128/AEM.71.12.8228-8235.2005
Lu, Y.-M., Jia, Z., Janus, C., Henderson, J. T., Gerlai, R., Wojtowicz, J. M., et al. (1997). Mice lacking metabotropic glutamate receptor 5 show impaired learning and reduced CA1 Long-Term Potentiation (LTP) But Normal CA3 LTP. J. Neurosci. 17, 5196–5205. doi: 10.1523/JNEUROSCI.17-13-05196.1997
Macedo e Cordeiro, T., Zhang, X., Graubics, K., Colwell, R., and Teixeira, A. L. (2020). Microbiome and schizophrenia: current evidence and future challenges. Curr. Behav. Neurosci. Rep. 7, 51–61. doi: 10.1007/s40473-020-00206-5
Maksymetz, J., Moran, S. P., and Conn, P. J. (2017). Targeting metabotropic glutamate receptors for novel treatments of schizophrenia. Mol. Brain 10:15. doi: 10.1186/s13041-017-0293-z
Mandal, S., Van Treuren, W., White, R. A., Eggesbø, M., Knight, R., and Peddada, S. D. (2015). Analysis of composition of microbiomes: a novel method for studying microbial composition. Microb. Ecol. Health Dis. 26:27663. doi: 10.3402/mehd.v26.27663
Matosin, N., Fernandez-Enright, F., Fung, S. J., Lum, J. S., Engel, M., Andrews, J. L., et al. (2015). Alterations of mGluR5 and its endogenous regulators Norbin, Tamalin and Preso1 in schizophrenia: towards a model of mGluR5 dysregulation. Acta Neuropathol. 130, 119–129. doi: 10.1007/s00401-015-1411-6
McCutcheon, R. A., Reis Marques, T., and Howes, O. D. (2020). Schizophrenia—An Overview. JAMA Psychiatry 77, 201–210. doi: 10.1001/jamapsychiatry.2019.3360
McMurdie, P. J., and Holmes, S. (2013). phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One 8:e61217. doi: 10.1371/journal.pone.0061217
Meckel, K. R., and Kiraly, D. D. (2019). A potential role for the gut microbiome in substance use disorders. Psychopharmacology 236, 1513–1530. doi: 10.1007/s00213-019-05232-0
Nau, R., Sörgel, F., and Eiffert, H. (2010). Penetration of drugs through the blood-cerebrospinal fluid/blood-brain barrier for treatment of central nervous system infections. Clin. Microbiol. Rev. 23, 858–883. doi: 10.1128/cmr.00007-10
Nguyen, T. L. A., Vieira-Silva, S., Liston, A., and Raes, J. (2015). How informative is the mouse for human gut microbiota research? Dis. Models Mech. 8, 1–16. doi: 10.1242/dmm.017400
Nguyen, T. T., Kosciolek, T., Maldonado, Y., Daly, R. E., Martin, A. S., McDonald, D., et al. (2019). Differences in gut microbiome composition between persons with chronic schizophrenia and healthy comparison subjects. Schizophr. Res. 204, 23–29. doi: 10.1016/j.schres.2018.09.014
Pan, R., Zhang, X., Gao, J., Yi, W., Wei, Q., and Su, H. (2020). Analysis of the diversity of intestinal microbiome and its potential value as a biomarker in patients with schizophrenia: a cohort study. Psychiatry Res. 291:113260. doi: 10.1016/j.psychres.2020.113260
Ploj, K., Albery-Larsdotter, S., Arlbrandt, S., Kjaer, M. B., Skantze, P. M. C., and Storlien, L. H. (2010). The metabotropic glutamate mGluR5 receptor agonist CHPG stimulates food intake. Neuroreport 21, 704–708. doi: 10.1097/WNR.0b013e32833b4fe7
Pyndt Jørgensen, B., Krych, L., Pedersen, T. B., Plath, N., Redrobe, J. P., Hansen, A. K., et al. (2015). Investigating the long-term effect of subchronic phencyclidine-treatment on novel object recognition and the association between the gut microbiota and behavior in the animal model of schizophrenia. Physiol. Behav. 141, 32–39. doi: 10.1016/j.physbeh.2014.12.042
Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2012). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596. doi: 10.1093/nar/gks1219
Ravussin, Y., Koren, O., Spor, A., LeDuc, C., Gutman, R., Stombaugh, J., et al. (2012). Responses of gut microbiota to diet composition and weight loss in lean and obese mice. Obesity 20, 738–747. doi: 10.1038/oby.2011.111
Rinninella, E., Raoul, P., Cintoni, M., Franceschi, F., Miggiano, G., Gasbarrini, A., et al. (2019). What is the healthy gut microbiota composition? A changing ecosystem across age, environment, diet, and diseases. Microorganisms 7:14. doi: 10.3390/microorganisms7010014
Rodrigues-Amorim, D., Rivera-Baltanás, T., Regueiro, B., Spuch, C., de las Heras, M. E., Vázquez-Noguerol Méndez, R., et al. (2018). The role of the gut microbiota in schizophrenia: current and future perspectives. World J. Biol. Psychiatry 19, 571–585. doi: 10.1080/15622975.2018.1433878
Rohart, F., Gautier, B., Singh, A., and Lê Cao, K.-A. (2017). mixOmics: an R package for ‘omics feature selection and multiple data integration. PLoS Comput. Biol. 13:e1005752. doi: 10.1371/journal.pcbi.1005752
Sams-Dodd, F. (1999). Phencyclidine in the social interaction test: an animal model of Schizophrenia with face and predictive validity. Rev. Neurosci. 10, 59–90. doi: 10.1515/REVNEURO.1999.10.1.59
Saunders, J. M., Moreno, J. L., Ibi, D., Sikaroodi, M., Kang, D. J., Muñoz-Moreno, R., et al. (2020). Gut microbiota manipulation during the prepubertal period shapes behavioral abnormalities in a mouse neurodevelopmental disorder model. Sci. Rep. 10:4697.
Schwarz, E., Maukonen, J., Hyytiäinen, T., Kieseppä, T., Orešič, M., Sabunciyan, S., et al. (2018). Analysis of microbiota in first episode psychosis identifies preliminary associations with symptom severity and treatment response. Schizophr. Res. 192, 398–403. doi: 10.1016/j.schres.2017.04.017
Severance, E. G., Prandovszky, E., Castiglione, J., and Yolken, R. H. (2015). Gastroenterology issues in Schizophrenia: why the gut matters. Curr. Psychiatry Rep. 17:27.
Shannon, P. (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498–2504. doi: 10.1101/gr.1239303
Shen, Y., Xu, J., Li, Z., Huang, Y., Yuan, Y., Wang, J., et al. (2018). Analysis of gut microbiota diversity and auxiliary diagnosis as a biomarker in patients with schizophrenia: a cross-sectional study. Schizophr. Res. 197, 470–477. doi: 10.1016/j.schres.2018.01.002
Skonieczna-Źydecka, K., Łoniewski, I., Misera, A., Stachowska, E., Maciejewska, D., Marlicz, W., et al. (2019). Second-generation antipsychotics and metabolism alterations: a systematic review of the role of the gut microbiome. Psychopharmacology 236, 1491–1512. doi: 10.1007/s00213-018-5102-6
Slattery, J. A., Page, A. J., Dorian, C. L., Brierley, S. M., and Blackshaw, L. A. (2006). Potentiation of mouse vagal afferent mechanosensitivity by ionotropic and metabotropic glutamate receptors: glutamatergic modulation of afferents. J. Physiol. 577, 295–306. doi: 10.1113/jphysiol.2006.117762
Stansley, B. J., and Conn, P. J. (2018). The therapeutic potential of metabotropic glutamate receptor modulation for Schizophrenia. Curr. Opin. Pharmacol. 38, 31–36. doi: 10.1016/j.coph.2018.02.003
Sun, L., Zhang, H., Cao, Y., Wang, C., Zhao, C., Wang, H., et al. (2019). Fluoxetine ameliorates dysbiosis in a depression model induced by chronic unpredicted mild stress in mice. Int. J. Med. Sci. 16, 1260–1270. doi: 10.7150/ijms.37322
Szeligowski, T., Yun, A. L., Lennox, B. R., and Burnet, P. W. J. (2020). The Gut Microbiome and Schizophrenia: the current state of the field and clinical applications. Front. Psychiatry 11:156. doi: 10.3389/fpsyt.2020.00156
Terbeck, S., Akkus, F., Chesterman, L. P., and Hasler, G. (2015). The role of metabotropic glutamate receptor 5 in the pathogenesis of mood disorders and addiction: combining preclinical evidence with human Positron Emission Tomography (PET) studies. Front. Neurosci. 9:86. doi: 10.3389/fnins.2015.00086
Uno, Y., and Coyle, J. T. (2019). Glutamate hypothesis in schizophrenia. Psychiatry Clin. Neurosci. 73, 204–215. doi: 10.1111/pcn.12823
Vandeputte, D., Falony, G., Vieira-Silva, S., Tito, R. Y., Joossens, M., and Raes, J. (2016). Stool consistency is strongly associated with gut microbiota richness and composition, enterotypes and bacterial growth rates. Gut 65, 57–62. doi: 10.1136/gutjnl-2015-309618
Volk, D. W., Eggan, S. M., and Lewis, D. A. (2010). Alterations in metabotropic glutamate receptor 1α and regulator of G protein signaling 4 in the prefrontal cortex in Schizophrenia. Am. J. Psychiatry 167, 1489–1498. doi: 10.1176/appi.ajp.2010.10030318
Xu, J., Zhu, Y., Contractor, A., and Heinemann, S. F. (2009). mGluR5 has a critical role in inhibitory learning. J. Neurosci. 29, 3676–3684. doi: 10.1523/JNEUROSCI.5716-08.2009
Xu, R., Wu, B., Liang, J., He, F., Gu, W., Li, K., et al. (2019). Altered gut microbiota and mucosal immunity in patients with schizophrenia. Brain Behav. Immun. 85, 120–127. doi: 10.1016/j.bbi.2019.06.039
Zeleznikow-Johnston, A. M., Renoir, T., Churilov, L., Li, S., Burrows, E. L., and Hannan, A. J. (2018). Touchscreen testing reveals clinically relevant cognitive abnormalities in a mouse model of schizophrenia lacking metabotropic glutamate receptor 5. Sci. Rep. 8:16412.
Zenewicz, L. A., Yin, X., Wang, G., Elinav, E., Hao, L., Zhao, L., et al. (2013). IL-22 deficiency alters colonic microbiota to be transmissible and colitogenic. J. Immunol. 190, 5306–5312. doi: 10.4049/jimmunol.1300016
Zeng, C., Yang, P., Cao, T., Gu, Y., Li, N., Zhang, B., et al. (2020). Gut microbiota: an intermediary between metabolic syndrome and cognitive deficits in schizophrenia. Prog. Neuropsychopharmacol. Biol. Psychiatry. doi: 10.1016/j.pnpbp.2020.110097 [Epub ahead of print].
Zheng, P., Zeng, B., Liu, M., Chen, J., Pan, J., Han, Y., et al. (2019). The gut microbiome from patients with schizophrenia modulates the glutamate-glutamine-GABA cycle and schizophrenia-relevant behaviors in mice. Sci. Adv. 5:eaau8317. doi: 10.1126/sciadv.aau8317
Zhu, F., Guo, R., Wang, W., Ju, Y., Wang, Q., Ma, Q., et al. (2019). Transplantation of microbiota from drug-free patients with schizophrenia causes schizophrenia-like abnormal behaviors and dysregulated kynurenine metabolism in mice. Mol. Psychiatry. doi: 10.1038/s41380-019-0475-4 [Epub ahead of print]
Keywords: gut dysbiosis, microbiota, microbiome, psychiatric disorders, schizophrenia, mGlu5 knockout mice
Citation: Gubert C, Kong G, Uzungil V, Zeleznikow-Johnston AM, Burrows EL, Renoir T and Hannan AJ (2020) Microbiome Profiling Reveals Gut Dysbiosis in the Metabotropic Glutamate Receptor 5 Knockout Mouse Model of Schizophrenia. Front. Cell Dev. Biol. 8:582320. doi: 10.3389/fcell.2020.582320
Received: 11 July 2020; Accepted: 08 October 2020;
Published: 29 October 2020.
Edited by:
Paola Brun, University of Padua, ItalyReviewed by:
Tomasz Kosciolek, Jagiellonian University, PolandZhihao Jia, Purdue University, United States
Copyright © 2020 Gubert, Kong, Uzungil, Zeleznikow-Johnston, Burrows, Renoir and Hannan. 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: Anthony J. Hannan, anthony.hannan@florey.edu.au