Analysis of gut microbiota in patients with cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL).

Cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL) is a major hereditary small vessel disease caused by mutations in NOTCH3. The variations in progression and severity among patients suggest that the CADASIL phenotype is modified by some genetic and environmental factors. Recent studies have shown the potential roles of gut microbiota in human diseases. We hypothesized that gut microbiota modifies the disease phenotype. We performed gut microbial meta 16S rRNA analysis of fecal samples from 15 CADASIL patients and 16 controls. The microbial α- and β-diversities and taxonomy were compared between CADASIL patients and controls and between CADASIL patients with and without an ischemic stroke history. No significant difference in α- or β-diversity was observed in either case-control or subgroup comparisons. In the taxonomic microbial analysis, there was a significant increase in abundance of 6 genera and significant decrease in 2 genera in CADASIL patients compared with controls. There was a significant decrease in abundance of 2 genera in CADASIL patients with compared with those without stroke. This is the first study on CADASIL focusing on gut microbiota. Our findings suggest that gut microbiota modifies the onset and progression of CADASIL.

Introduction

Drowing evidence supports the potential roles of gut microbiota in human health and diseases, including Parkinson’s disease, multiple sclerosis, and atherosclerotic disease.( It may be difficult to clarify the influence of gut microbiota on sporadic cerebral infarction because stroke patients show different conventional risk factors, including hypertension, smoking, diabetes, hyperlipidemia, and aging. In contrast, patients with cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL), a major hereditary small vessel disease caused by mutations in NOTCH3,( generally have no common vascular risk factors. Typical clinical characteristics of CADASIL patients are migraines, recurrent ischemic strokes, vascular dementia, and specific white matter lesions with multiple lacunar infarcts.( It is well-known that the progression and severity of CADASIL show both inter- and intra-familial variations among patients. Age at onset of the first stroke ranges from 20–70 years old, and moreover, some patients remain asymptomatic even after 70 years old.( These suggest that the CADASIL phenotype is modified by some genetic or environmental factors in addition to NOTCH3 mutation. To clarify the contribution of gut microbiota to cerebral infarction, we focused on the profile of fecal microbiota in CADASIL patients.

Materials and Methods

Participants

We recruited nineteen Japanese CADASIL patients who visited our outpatient clinic at University Hospital, Kyoto Prefectural University of Medicine from January to November 2017. Eighteen of their family members (15 spouses, 2 siblings, and a child) were recruited as controls. Nutrition survey during the one month before and up to fecal sampling, using the brief-type self-administered diet history questionnaire (BDHQ).( Four CADASIL patients and 2 controls were excluded from the analysis because of taking medication that may interfere with gut microbiota (antibiotics: n = 1, proton pump inhibitors: n = 4, Metformin: n = 1). We finally analyzed fecal samples from 15 CADASIL patients and 16 controls. CADASIL patients were divided into two subgroups: those with a history of symptomatic ischemic stroke (n = 7) and those without it (n = 8). The Ethics Committee of Kyoto Prefectural University of Medicine approved the research protocol (permission No. ERB-C-725-1), and all participants provided written informed consent prior to enrollment.

Meta16S analysis of gut microbiota

Meta16S analysis of gut microbiota was performed according to Takagi et al.( In brief, fecal sampling was performed using a sampling kit (Techno Suruga Lab., Shizuoka, Japan). Meta16S analysis of fecal DNA was performed at the Biomedical Center, Takara Bio. PCR for DNA extracted from fecal samples was comprised of two steps. In the first step, PCR was performed to amplify with primers of 341F (5'-TCGTCGGCAGCGTCAGATGTGTATAAGAGACAGCCTACGGGNGGCWGCAG-3') and 806R (5'-GTCTCGTGGGCTC GGAGATGTGTATAAGAGACAGGGACTACHVGGGTWTCTAAT-3') corresponding to the V3–V4 region of the 16S rRNA gene. In the second PCR step, the index sequences for Illumina sequencer were added. Then, the libraries were subjected to the sequencing of paired-end 300 bases on MiSeq (Illumina).

Processing of sequence data, operational taxonomic unit (OTU) definition, and taxonomy assignment were performed using QIIME ver. 1.9, USEARCH ver. 8.0, and UCHIME ver. 4.2.40, respectively. Singletons were removed in the present study. Taxonomy assignment of the resulting OTU was completed using RDP classifier ver. 2.10.2 with the Greengenes database (May 2013).

Microbiota diversity within a sample (α-diversity) was evaluated as the Chao1 index (OTU richness estimation) and Shannon index (OTU evenness estimation) calculated by the R “phyloseq” package. Differences in diversity between samples (β-diversity) was estimated using the UniFrac metric by quantitative (weighted) or qualitative (unweighted) analysis, and the distances between the samples were visualized by principal coordinate analysis (PCoA) plotting, using the software QIIME (ver. 1.9.0).

Statistical analysis

Data were compared between groups (CADASIL patients vs controls, and CADASIL patients with and without stroke). Nutrition intake profiles were compared by the Mann-Whitney U test using SPSS 17.0. The relative abundance of bacterial genera, OTU, and α-diversity were compared with a Wilcoxon rank sum test using JMP14.0.0. Clustering on PCoA plot (β-diversity) and was statistically analyzed using permutational multivariate analysis of variance (PERMANOVA). Values of p<0.05 were considered significant.

Results

Backgrounds and nutrition survey

There was no significant difference in backgrounds and vascular risk factors between CADASIL patients and controls, or CADASIL subgroups, except for weak difference in gender between CADASIL patients with stroke and those without stroke (p = 0.04) (Table 1).

Table 1

Background of the participants

	CADASIL patients (n = 15)	Normal controls (n = 16)	p value	CADASIL patients		p value
				Stroke (+) (n = 7)	Stroke (–) (n = 8)
Age (range)	56.9 (45–74)	53.7 (28–71)	0.33	57.0 (45–66)	56.9 (47–74)	0.98
BMI	22.9 (18.0–33.0)	22.0 (18.2–29.4)	0.44	23.9 (20.7–30.9)	22.0 (18.0–28.4)	0.3
Male/Female	8 (53%)/7 (47%)	8 (50%)/8 (50%)	1	6 (86%)/1 (14%)	2 (25%)/6 (75%)	0.04
Hypertension	2/15 (13%)	2/16 (13%)	1	2/7 (29%)	0/8 (0%)	0.2
Dyslipidemia	5/15 (33%)	2/16 (13%)	0.22	3/7 (43%)	2/8 (25%)	0.61
Diabetes mellitus	0/15 (0%)	1/16 (6 %)	1	0/7 (0%)	0/8 (0%)	—
Smoking	1/15 (7%)	2/16 (13%)	1	1/7 (14%)	0/8 (0%)	0.47

Mean of age and body mass index (BMI), number of male (%)/female (%), and number of observed/studied (%) are shown. For p value calculation, the t test (age and BMI), or Fisher’s exact test was employed.

The nutrition survey was performed using a BDHQ.( Between CADASIL patients and controls, there was no significant difference in 97 nutrition items. Between CADASIL subgroups, no significant difference was noted except for the daily intake of arachidonic acid (patients with stroke: 22.28 mg/day vs patients without it: 29.39 mg/day; p = 0.028).

Analysis of gut microbiota diversity

In total, seven hundred ninety OTUs were detected after the removal of singletons in this study. No significant difference in α-diversity was observed between the CADASIL patients and controls, or between the patients with and without stroke (Fig. 1). As for β-diversity, PCoA plots showed no significant difference in clusters between CADASIL patients and controls, or between CADASIL patients with and without stroke (Fig. 2).

Fig. 1

Case-control and subgroup comparison of α-diversity of gut microbiota. Box and whisker plots of the Chao1 index (upper) and Shannon index (lower). These α-diversity indices were compared between CADASIL patients and controls, and between CADASIL patients with previous ischemic stroke and those without it (Wilcoxon rank sum test).

Fig. 2

Case-control and subgroup comparisons of β-diversity of gut microbiota. Distribution of β-diversity was visualized by principle coordinate analysis (PCoA) plots from weighted (A) and unweighted (B) UniFrac metrics, and compared between CADASIL patients and controls (left), and between CADASIL patients with previous ischemic stroke and those without it (right) (permutational multivariate analysis of variance).

Taxonomic microbial analysis

The taxonomic change in the microbiota was assessed at genus and OTU levels. In CADASIL patients compared with controls, there was a significant increase in abundance of 6 genera (Lachnospira, Odoribacter, Parvimonas, unclassified genus belonging to Barnesiellaceae and Lachnospiraceae, unclassified genus belonging to order SHA-98), and a significant decrease in 2 genera (Megasphaera and Acidaminococcus) (Table 2). Regarding the OTU level, there was a significant increase in abundance of 24 OTUs and significant decrease in 4 OTUs in CADASIL patients compared with controls (Supplemental Table 1*).

Table 2

Genus level taxonomic analysis of gut microbiota

Case-control
Phylum	Class	Order	Family	Genus	Average ± STDEV (%)			Wilcoxon
					Patients (n = 15)	Controls (n = 16)		p value
Firmicutes	Clostridia	Clostridiales	Veillonellaceae	Megasphaera	0.000171 ± 0.000513	0.281 ± 0.597		0.001**
Bacteroidetes	Bacteroidia	Bacteroidales	[Barnesiellaceae]	Unclassified	0.479 ± 4.69	0.0880 ± 0.148		0.008**
Firmicutes	Clostridia	Clostridiales	Lachnospiraceae	Unclassified	0.132 ± 0.142	0.0326 ± 0.0525		0.011*
Firmicutes	Clostridia	Clostridiales	Lachnospiraceae	Lachnospira	1.44 ± 2.11	0.183 ± 0.348		0.011*
Firmicutes	Clostridia	Clostridiales	Veillonellaceae	Acidaminococcus	0.0814 ± 0.315	0.405 ± 0.754		0.027*
Firmicutes	Clostridia	SHA-98	Unclassified	Unclassified	0.00459 ± 0.00654	0.000699 ± 0.00176		0.032*
Bacteroidete	Bacteroidia	Bacteroidales	[Odoribacteraceae]	Odoribacter	0.234 ± 0.224	0.0981 ± 0.0931		0.0326*
Firmicutes	Clostridia	Clostridiales	[Tissierellaceae]	Parvimonas	0.00287 ± 0.00320	0.00153 ± 0.0038		0.036*

The Genera of significant differences (p<0.05, Wilcoxon rank test) between CADASIL patients and controls (case-control) and between the CADASIL patients with and without stroke (subgroup) are shown. *p<0.05, **p<0.01.

In the CADASIL patients with compared with those without stroke, there was a significant decrease in abundance of 2 genera, Phascolarctobacterium and Paraprevotella (Table 2).

Regarding the OTU level, a significant increase in abundance of 13 OTUs and significant decrease in 3 OTUs were observed in CADASIL patients with compared with those without stroke (Supplemental Table 1*).

Discussion

To our knowledge, this is the first study addressing the possibility of gut microbiota contributing to the development of CADASIL. The major finding of this study was identification of significant taxonomic differences in gut microbiota between CADASIL patients and controls, and also between patients with and without previous ischemic stroke. These suggest the role of gut microbiota as modifiers of the onset and progression of CADASIL.

Taxonomic analysis showed significant differences in abundances of some bacteria between CADASIL patients and controls. According to previous reports, some of the representative bacteria shown in Table 1 might be associated with stroke onset in CADASIL patients. Firstly, Goodrich et al.( reported that the abundance of SHA-98, a member of the Christensenellaceae consortium, was associated with a gene encoding aldehyde dehydrogenase (ALDH1L1). Interestingly, an SNP in ALDH1L1 was reported to be associated with ischemic stroke in the Framingham Heart Study.( Secondly, Parvimonas, a Gram-positive anaerobic coccus, was also of interest according to recent knowledge on the association between periodontal disease and ischemic stroke.( Fak et al.( reported that bacteria from the oral cavity and gut can be recovered from atherosclerotic plaque, and that the abundance of oral Parvimonas was positively associated with c-reactive protein (CRP) in atherosclerotic patients. Thirdly, the family Lachnospiraceae has been reported to facilitate Treg differentiation and stimulate TGF-β and IL-10 production by immune cells.( TGF-β is a key molecule in the pathogenesis of CARASIL, an autosomal recessive small vessel disease caused by HTRA1 mutations.( Dysregulation of TGF-β signaling may also be associated with CADASIL according to a report that latent TGF-β binding protein 1 was sequestrated into CADASIL-related NOTCH3 extracellular domain deposits.( Taken together, the increase of Lachnospiraceae in CADASIL patients is of interest.

In recent gut microbial taxonomic studies, other bacteria were reported to be associated with ischemic stroke in sporadic patients.( In this study, however, we did not replicate any of these findings.

On the other hand, on taxonomic analysis in association with stroke among CADASIL patients, two genera, Phascolarctobacterium and Paraprevotella, showed a significant difference. Based on our search, however, no previous report addressed the association of these genera with cerebrovascular dysfunction.

A major limitation of our study was the small number of participants, due to the fact that CADASIL is a rare hereditary disease. In this regard, it is of note that Inoue et al.( reported robust evaluation of gut dysbiosis in type 2 diabetic patients even though the sample size, 12 patients and 10 controls, was as small as our study. By using 16S rRNA metagenomic analysis in combination with functional prediction software Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt) and Kyoto Encyclopedia of Genes and Genomes (KEGG) database, they identified not only composition but also functional profiles of gut microbiota associated with diabetes pathophysiology.( Then, we additionally performed similar analysis. Although there was no significant difference between CADASIL patients and normal controls, we found 32 KEGG pathways of which abundances were significantly different between CADASIL patients with and without stroke (Supplemental Fig. 1*). To interpret these findings, further studies including fecal metabolome analysis should be necessary. However, of the 32 pathways, lipid metabolism seems of note, because Žitňanová et al.( recently reported that subfractions of plasma lipid were significantly altered in acute phase of stroke, another limitation was that we did not adjust background factors to compare gut microbiota between groups. Previous report showed that several changes in gut microbiota were associated with age and sex, by analyzing 277 healthy individuals.( Moreover, it cannot be concluded whether taxonomic changes between CADASIL and controls or between CADASIL patients with and without stroke were causes or results of the disease. To address these issues, increasing the number of participants and longitudinal analysis of each patient are necessary.

In conclusion, this is the first study on CADASIL focusing on gut microbiota. Significant differences in taxonomy were noted between CADASIL patients and controls, and between CADASIL patients with and without stroke, suggesting that gut microbiota modifies the onset and progression of CADASIL.

Author Contributions

TM and YN supervised and designed the study; JM and TM collected feces; JM, AW, TK, MM, IM and TM collected the clinical information; JM, RI, TT, IM, YN and TM analyzed and interpreted the gut microbiota data; JM, SW and TM analyzed and interpreted the nutrition data; JM, RI and SW performed the statistical analysis. All authors drafted and edited the manuscript, and discussed the results and commented on the manuscript.