Objective

To identify features of the gut microbiome associated with multiple sclerosis activity over time.

Methods

We used 16S ribosomal RNA sequencing from stool of 55 recently diagnosed pediatric-onset multiple sclerosis patients. Microbiome features included the abundance of individual microbes and networks identified from weighted genetic correlation network analyses. Prentice-Williams-Peterson Cox proportional hazards models estimated the associations between features and three disease activity outcomes: clinical relapses and both new/enlarging T2 lesions and new gadolinium-enhancing lesions on brain MRI. Analyses were adjusted for age, sex, and disease-modifying therapies.

Results

Participants were followed, on average, 2.1 years. Five microbes were nominally associated with all three disease activity outcomes after multiple testing correction. These included butyrate producers Odoribacter (relapse hazard ratio = 0.46, 95% confidence interval: 0.24, 0.88) and Butyricicoccus (relapse hazard ratio = 0.49, 95% confidence interval: 0.28, 0.88). Two networks of co-occurring gut microbes were significantly associated with a higher hazard of both MRI outcomes (gadolinium-enhancing lesion hazard ratios (95% confidence intervals) for Modules 32 and 33 were 1.29 (1.08, 1.54) and 1.42 (1.18, 1.71), respectively; T2 lesion hazard ratios (95% confidence intervals) for Modules 32 and 33 were 1.34 (1.15, 1.56) and 1.41 (1.21, 1.64), respectively). Metagenomic predictions of these networks demonstrated enrichment for amino acid biosynthesis pathways.

Interpretation

Both individual and networks of gut microbes were associated with longitudinal multiple sclerosis activity. Known functions and metagenomic predictions of these microbes suggest the important role of butyrate and amino acid biosynthesis pathways. This provides strong support for future development of personalized microbiome interventions to modify multiple sclerosis disease activity.

Background

Rare genetic variants are emerging as important contributors to the heritability of multiple sclerosis (MS). Whether rare variants also contribute to pediatric-onset multiple sclerosis (POMS) is unknown.

Objective

To test whether genes harboring rare variants associated with adult-onset MS risk (PRF1, PRKRA, NLRP8, and HDAC7) and 52 major histocompatibility complex (MHC) genes are associated with POMS.

Methods

We analyzed DNA samples from 330 POMS cases and 306 controls from the US Network of Pediatric MS Centers and Kaiser Permanente Northern California for which Illumina ExomeChip genotypes were available. Using the gene-based method "SKAT-O," we tested the association between candidate genes and POMS risk.

Results

After correction for multiple comparisons, one adult-onset MS gene (PRF1, p = 2.70 × 10^-3) and two MHC genes (BRD2, p = 5.89 × 10^-5 and AGER, p = 7.96 × 10^-5) were significantly associated with POMS. Results suggest these are independent of HLA-DRB1*1501.

Conclusion

Findings support a role for rare coding variants in POMS susceptibility. In particular, rare minor alleles within PRF1 were more common among individuals with POMS compared to controls while the opposite was true for rare variants within significant MHC genes, BRD2 and AGER. These genes would not have been identified by common variant studies, emphasizing the merits of investigating rare genetic variation in complex diseases.