Skip to main content
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mBio
    • Editor in Chief
    • Board of Editors
    • AAM Fellows
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
mBio
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mBio
    • Editor in Chief
    • Board of Editors
    • AAM Fellows
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
Research Article

Intestinal Dysbiosis Associated with Systemic Lupus Erythematosus

Arancha Hevia, Christian Milani, Patricia López, Adriana Cuervo, Silvia Arboleya, Sabrina Duranti, Francesca Turroni, Sonia González, Ana Suárez, Miguel Gueimonde, Marco Ventura, Borja Sánchez, Abelardo Margolles
Nicola Segata, Invited Editor, Maria Gloria Dominguez Bello, Editor
Arancha Hevia
aInstituto de Productos Lácteos de Asturias, Consejo Superior de Investigaciones Científicas, Villaviciosa, Asturias, Spain
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Christian Milani
bLaboratory of Probiogenomics, Department of Life Sciences, University of Parma, Parma, Italy
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Patricia López
cImmunology Area, Department of Functional Biology, University of Oviedo, Asturias, Spain
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Adriana Cuervo
dPhysiology Area, Department of Functional Biology, University of Oviedo, Asturias, Spain
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Silvia Arboleya
aInstituto de Productos Lácteos de Asturias, Consejo Superior de Investigaciones Científicas, Villaviciosa, Asturias, Spain
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Sabrina Duranti
bLaboratory of Probiogenomics, Department of Life Sciences, University of Parma, Parma, Italy
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Francesca Turroni
bLaboratory of Probiogenomics, Department of Life Sciences, University of Parma, Parma, Italy
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Sonia González
dPhysiology Area, Department of Functional Biology, University of Oviedo, Asturias, Spain
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ana Suárez
cImmunology Area, Department of Functional Biology, University of Oviedo, Asturias, Spain
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Miguel Gueimonde
aInstituto de Productos Lácteos de Asturias, Consejo Superior de Investigaciones Científicas, Villaviciosa, Asturias, Spain
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Marco Ventura
bLaboratory of Probiogenomics, Department of Life Sciences, University of Parma, Parma, Italy
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Borja Sánchez
aInstituto de Productos Lácteos de Asturias, Consejo Superior de Investigaciones Científicas, Villaviciosa, Asturias, Spain
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Abelardo Margolles
aInstituto de Productos Lácteos de Asturias, Consejo Superior de Investigaciones Científicas, Villaviciosa, Asturias, Spain
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Nicola Segata
University of Trento
Roles: Invited Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Maria Gloria Dominguez Bello
New York University School of Medicine
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/mBio.01548-14
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

Systemic lupus erythematosus (SLE) is the prototypical systemic autoimmune disease in humans and is characterized by the presence of hyperactive immune cells and aberrant antibody responses to nuclear and cytoplasmic antigens, including characteristic anti–double-stranded DNA antibodies. We performed a cross-sectional study in order to determine if an SLE-associated gut dysbiosis exists in patients without active disease. A group of 20 SLE patients in remission, for which there was strict inclusion and exclusion criteria, was recruited, and we used an optimized Ion Torrent 16S rRNA gene-based analysis protocol to decipher the fecal microbial profiles of these patients and compare them with those of 20 age- and sex-matched healthy control subjects. We found diversity to be comparable based on Shannon’s index. However, we saw a significantly lower Firmicutes/Bacteroidetes ratio in SLE individuals (median ratio, 1.97) than in healthy subjects (median ratio, 4.86; P < 0.002). A lower Firmicutes/Bacteroidetes ratio in SLE individuals was corroborated by quantitative PCR analysis. Notably, a decrease of some Firmicutes families was also detected. This dysbiosis is reflected, based on in silico functional inference, in an overrepresentation of oxidative phosphorylation and glycan utilization pathways in SLE patient microbiota.

IMPORTANCE Growing evidence suggests that the gut microbiota might impact symptoms and progression of some autoimmune diseases. However, how and why this microbial community influences SLE remains to be elucidated. This is the first report describing an SLE-associated intestinal dysbiosis, and it contributes to the understanding of the interplay between the intestinal microbiota and the host in autoimmune disorders.

INTRODUCTION

Metagenomic studies on gut microbiota burst onto the scientific scene during the last decade, due to the advent of next-generation sequencing techniques. In a very short period of time, microbiologists moved from the study of single, isolated, cultivable microorganisms, specifically, those able to grow under standard laboratory conditions, to the investigation of very complex microbial communities, mainly composed of uncultivable bacteria (1, 2). The first metagenomics reports enabled an overview of the complexity of our gut microbial communities (3, 4). Further studies focused on establishing the correlation between the human gut microbiome, the collective genomes of all microbes inhabiting the gut (5), and different physiological states, including those having an influence on health. Currently, we know that the gut microbiota might affect food and drug metabolism (6), influences human behavior (7), shifts during the course of pregnancy (8), displays age-associated changes (9–12), and possesses distinctive features depending on geographical location (12, 13), among other features. It is also becoming clear that there is a strong link between dietary patterns and the gut microbial profile (14, 15). Furthermore, some links have been established between some disorders (for example, obesity and metabolic syndrome) and an imbalance in the gut microbial ecology, also called dysbiosis (16–18). Remarkably, intestinal dysbiosis has also been associated with autoimmune diseases, such as rheumatoid arthritis, type 1 diabetes, and inflammatory bowel disease (IBD) (19–21).

Systemic lupus erythematous (SLE) is a prototypical autoimmune disease in humans that is characterized by the presence of hyperactive immune cells and aberrant antibody responses to nuclear and cytoplasmic antigens. Genetic, immunological, hormonal, and environmental factors contribute to disease susceptibility (22), and its prevalence varies greatly depending on the population under study, although a prevalence of 2 to 5 cases per 10,000 inhabitants is reportedly considered normal (23). Among the environmental factors, growing evidence suggests that molecular mimicry as a result of viral infection may contribute to the development of lupus (24). Also, some reports have highlighted intestinal infections that may ameliorate SLE symptoms (25), and a marked difference in the specificity of antibodies to bacterial DNA in healthy people and SLE patients has been indicated (26). In fact, there is early evidence of a different abundance of cultivable intestinal bacteria in SLE (27). Remarkably, it has been suggested that novel SLE biomarkers can be potentially found in the human microbiota (28). However, a study of the potential dysbiosis associated with SLE had not been tackled until now. Therefore, in this report we took advantage of next-generation sequencing techniques to explore the potential interplay of the human microbiome and SLE. We have proven, for the first time, that there is a gut microbial dysbiosis associated with SLE.

RESULTS AND DISCUSSION

Despite all the scientific knowledge generated in the last few years, and although few studies published so far support the dysbiosis theory as a key factor promoting chronic inflammation in autoimmune diseases (29–32), there is no scientific work that has taken advantage of next-generation sequencing techniques to explore the potential interplay of the human microbiome and SLE, the prototypical autoimmune disease in humans. Therefore, we designed our work with the aim of answering if there is an SLE-associated intestinal dysbiosis and, if so, which microbial population groups are related to the dysbiosis.

We defined the SLE population group by considering that there is a census of about 300 SLE patients in Asturias (from a total population of about 1,000,000 inhabitants). Thus, we were able to obtain a group of SLE patients from a well-defined geographical location to compare them with a similar group of healthy controls (HC), considering factors such as sex, age, medication (absence of antibiotic, steroid, and immunological treatments during the last 6 months), medical history (presenting a wide variety of clinical SLE manifestations), duration of the disease (2 to 24 years), and absence of flares of disease activity at the time of sampling (systemic lupus erythematosus disease activity index [SLEDAI] score of ≤8 at the time of sample collection). The group of SLE patients included individuals with a large variety of symptoms (Table 1), allowing us to establish correlations between the microbial profile and SLE, which are very likely independent of a specific pattern of symptoms. This variability in the phenotype of the disease is an intrinsic characteristic of SLE (22). We also selected patients with no active disease at the time of sampling, because the clinical manifestations of the disease in this population group are not biased by the pharmacological treatment necessary to treat SLE individuals during disease relapse. Furthermore, mean dietary intakes of energy, macronutrients, micronutrients, fiber, and phyto-compounds were recorded, both from patients with SLE and healthy subjects, and we found that there was no significant difference between the 2 groups (Table 2). Also, no significant difference was found between the 2 groups regarding lifestyle-related factors (smoking, alcohol consumption, physical activity, and use of vitamin and mineral supplements [data not shown]). This reduced the possibility that our analysis was affected by factors shown to have an influence on the gut microbial profile, such as age (9), diet (15), or phenolic compound intake (33, 34).

View this table:
  • View inline
  • View popup
TABLE 1 

Demographic, clinical, and immunological features of SLE patients

View this table:
  • View inline
  • View popup
TABLE 2 

General characteristics and mean dietary intake of energy, macronutrients, fiber forms, and phyto-compounds in patients with SLE and healthy controls

Our work is based on 16S rRNA gene-based data for fecal microbiota and the bioinformatic analysis of the results. In a previous work (35), we optimized protocols to study the human fecal microbial population by using an Ion Torrent PGM sequencing platform. This methodology was applied in the current study, and we obtained an average of 592,305 high-quality reads per fecal sample (see Table S1 in the supplemental material). Rarefaction curves obtained by plotting the Shannon, Chao1, and phylogenetic diversity indexes against the number of sequences (see Fig. S1 in the supplemental material) showed that a large part of the diversity of the samples was detected. The microbiota composition at the phylum and family levels was obtained (Fig. 1; see also Fig. S2 in the supplemental material). Remarkably, even considering the broad heterogeneity of the clinical manifestations of SLE in the individuals under study (Table 1), we obtained a particular type of microbiota for the SLE group. In this regard, the presence of anti–double-stranded DNA (dsDNA) antibodies and other clinical data were organized in a metadata file for all the microbiota profiles. A principal component analysis (PCoA) was performed with both metadata/microbiota profiles, using the variability of the 16S rRNA gene profiling at different taxonomic levels. Sample classification according to the metadata revealed no specific clustering of the samples or correlations with the different clinical features or anti-dsDNA antibodies (data not shown).

Figure S1

Rarefaction curves generated for 16S rRNA gene sequences obtained from fecal samples of control subjects (HC samples; gray lines) and lupus-affected subjects (SLE samples; red lines). (a) Rarefaction curves derived using the Shannon index. (b) Rarefaction curves derived using the Chao1 index. (C) Number of phylotypes identified versus the number of sequences per sample. Download Figure S1, DOC file, 0.2 MB.
Copyright © 2014 Hevia et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution-Noncommercial-ShareAlike 3.0 Unported license, which permits unrestricted noncommercial use, distribution, and reproduction in any medium, provided the original author and source are credited.

Table S1

Sequence data features Table S1, PDF file, 0.04 MB.
Copyright © 2014 Hevia et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution-Noncommercial-ShareAlike 3.0 Unported license, which permits unrestricted noncommercial use, distribution, and reproduction in any medium, provided the original author and source are credited.

FIG 1 
  • Open in new tab
  • Download powerpoint
FIG 1 

Aggregate microbiota composition in fecal samples from control (HC) and lupus-affected (SLE) subjects at the phylum level (a) and family level (b). In panel b, only taxonomic groups representing >0.5% are shown.

In silico analysis of the sequences highlighted the key findings of our work. A high-quality filtering approach was used in order to process the Ion Torrent-generated sequencing data (see Table S1 in the supplemental material); a total of 293,436 unfiltered operational taxonomic units (OTUs) were identified by using uclust for de novo OTU picking. Based on each of five alpha-diversity measures (Chao1, PD whole tree, observed species, Shannon, and Simpson indexes), patients and controls were not significantly different (data not shown). Notably, one of the main results was the identification of a clear dysbiosis between the two study groups which was characterized by a higher relative abundance of Bacteroidetes in the SLE group. Overall, we detected a significant decrease in the Firmicutes/Bacteroidetes ratio in SLE individuals: the microbiota of SLE patients, compared with controls, had an almost-2.5-fold-decreased ratio (Fig. 2A) (P < 0.002). Looking at the different individuals, a gradient from lower (SLE) to higher (HC) Firmicutes/Bacteroidetes ratios was observed (Fig. 2B). These 16S rRNA gene-based analyses were corroborated by quantitative PCR (qPCR) analysis. The levels (reported as the log cells/g, with the interquartile range [IQR] in parentheses) of total fecal bacteria were 10.62 (9.46 to 10.80) in SLE patients and 10.35 (10.07 to 10.59) in controls. Firmicutes levels in the SLE and control groups reached 9.69 (8.86 to 10.38) and 9.99 (9.68 to 10.31), respectively, while those of Bacteroidetes were 10.52 (9.56 to 10.83) and 9.89 (9.59 to 10.23), respectively. No statistically significant differences in these levels of microbial groups were found between SLE patients and controls. However, when the data were expressed as the relative percentages of Firmicutes and Bacteroidetes with respect to the total bacterial level, a significantly higher (P < 0.05) percentage of Bacteroidetes was observed in the SLE group. Moreover, when the Firmicutes/Bacteroidetes ratio was calculated, a statistically significant (P < 0.01) decrease in the SLE group with respect to the control group was found (0.94 [0.90 to 0.98] versus 1.01 [0.96 to 1.06], respectively). The differences in the ratios based on qPCR were less pronounced than the differences obtained with the 16S rRNA profiling, likely because the two techniques provide different kinds of information: a relative proportion of sequences (from the 16S rRNA gene-based analysis) versus an absolute quantification of sequences (via qPCR). Thus, the fact that we obtained clear evidence for a significantly lower Firmicutes/Bacteroidetes ratio in SLE patients when we used two different culture-independent techniques supports the soundness and reliability of our analyses. The phyla Bacteroidetes and Firmicutes include the most abundant components of the human gut microbiota (36). Dysbiosis between Firmicutes and Bacteroidetes in the human gut has been described in previous studies in association with some disorders. The ratio between Firmicutes and Bacteroidetes decreases in human type 2 diabetes compared with controls (37). Also, most studies of the microbiota in people with Crohn’s disease report a decrease in the abundance of Firmicutes and an increase in Bacteroidetes in association with the disease (38). An opposite situation is observed in obesity, in which the dysbiosis is characterized by an increase in the Firmicutes/Bacteroidetes ratio (17). Therefore, this specific microbial balance between the more abundant phyla in the human gut seems to be dependent on the physiology of the disorder. In relation to this, it has been reported that this equilibrium is susceptible to modification by shifts in the dietary pattern. Wu et al. (15) reported that Firmicutes levels were positively associated with a low-fat/high-fiber diet. Also, dietary interventions, including whole grain in the diet, increase the Firmicutes/Bacteroidetes ratio (39).

FIG 2 
  • Open in new tab
  • Download powerpoint
FIG 2 

(A) Box plot of Firmicutes/Bacteroidetes ratios (median ± IQR) in SLE patients versus healthy controls. (B) Percentages of 16S rRNA reads of Bacteroidetes (red bars) and Firmicutes (blue bars) in the DNA extracted from fecal samples of SLE patients (SLE codes) and healthy controls (HC codes). Ratios are significantly different (P < 0.002).

We prepared scatter plots (see Fig. S3A in the supplemental material), and they clearly highlighted a positive association between Firmicutes in healthy controls (P < 0.01) and Bacteroidetes in SLE patients (P < 0.001). This association was confirmed at lower taxonomic levels, and normalized abundances of the classes Bacteroidia (phylum Bacteroidetes) and Clostridia (phylum Firmicutes) and the orders Bacteroidales and Clostridiales differed significantly between the SLE patients and the healthy controls (see Fig. S3B). At the family level, Lachnospiraceae (P < 0.05) and Ruminococcaceae (P < 0.05) were positively associated with healthy controls (see Fig. S4 and S5 and Table S2 in the supplemental material).

Figure S3

(A) Scatter plot at the phylum level in which the relative abundances (means) of Firmicutes and Bacteroidetes in healthy controls (x axis) and SLE patients (y axis) are represented. (B) The normalized abundances of the classes Bacteroidia and Clostridia and the orders Bacteroidales and Clostridiales differed significantly between the lupus (SLE) and the healthy control (HC) 16S rRNA sequence data, as represented in the error plots. Tables with relative abundances and sequence metadata were elaborated using QIIME and analyzed using STAMP v 2.0.0. Two-sided Welch’s tests were run on every pair of means. The axes represent mean proportions of the taxa in HC (x axis) or SLE (y axis) groups. The vertical and horizontal lines in each taxa represent data spread (2nd and 98th percentiles) in the HC group (horizontal lines) and in the SLE group (vertical lines). The statistical fit of the points to the gray dashed y = x line is represented by the coefficient of determination in the upper left part of the plot (R2). The y = x line represents the boundary delimiting taxa that are enriched in the HC group (blue dots) or in the SLE group (yellow dots). Download Figure S3, DOC file, 0.3 MB.
Copyright © 2014 Hevia et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution-Noncommercial-ShareAlike 3.0 Unported license, which permits unrestricted noncommercial use, distribution, and reproduction in any medium, provided the original author and source are credited.

Figure S4

Scatter plot at the family level in which the relative abundances (means) of the bacterial families in healthy controls (x axis) and lupus patients (y axis) are represented. Families showing the greatest relative abundances and divergence from the diagonal are labeled with arrows. Download Figure S4, DOC file, 0.2 MB.
Copyright © 2014 Hevia et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution-Noncommercial-ShareAlike 3.0 Unported license, which permits unrestricted noncommercial use, distribution, and reproduction in any medium, provided the original author and source are credited.

Figure S5

Summaries of the communities present in the different samples analyzed in this study. Plots were obtained at four taxonomic levels: phylum (A), class (B), order (C), and family (D), and the most abundant taxons are labeled with arrows. Data were obtained using QIME v 1.7.0 with the scripts pick_de_novo_otus.py and summarize_taxa_through_plots.py. Download Figure S5, DOC file, 0.5 MB.
Copyright © 2014 Hevia et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution-Noncommercial-ShareAlike 3.0 Unported license, which permits unrestricted noncommercial use, distribution, and reproduction in any medium, provided the original author and source are credited.

Table S2

Statistically significant differences between healthy controls and SLE patients at the phylum and family levels Table S2, PDF file, 0.4 MB.
Copyright © 2014 Hevia et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution-Noncommercial-ShareAlike 3.0 Unported license, which permits unrestricted noncommercial use, distribution, and reproduction in any medium, provided the original author and source are credited.

Statistical differences between the two groups (HC and SLE) were calculated by a PERMANOVA test, with the distance data obtained after ordination analysis (PCoA) using PAST v 3.1 (40). In all cases, data from the relative taxa abundances were used, and distances were computed according to the Bray-Curtis similarity index. The groups of SLE and HC differed statistically whenever phylum- or family-level data were used (P < 0.01 or P < 0.02, respectively).

Unsupervised PCoA of the 16S rRNA sequence data identified the phyla Bacteroidetes and Firmicutes as the main gradients for SLE and healthy control groups, respectively (Fig. 3A). PCoA at the family level showed that the families Lachnospiraceae and Ruminococcaceae were located near the healthy controls (Fig. 3B). Interestingly, Lachnospiraceae was the most abundant family in the feces of both study groups (Fig. 1; see also Fig. S2 in the supplemental material). Depletion of Lachnospiraceae and Ruminococcaceae has been associated with Clostridium difficile infections and nosocomial diarrhea (41). Also, several studies showed a decrease of Lachnospiraceae in IBD patients, and this family has been suggested as a biomarker of disease activity (42–44). However, the functional consequences of the depletion of these bacteria in the previously mentioned intestinal diseases remain to be investigated. Although the amplicons of the 16S rRNA sequences could be relatively short to perform a totally reliable population structure analysis at the genus/species level, it is noteworthy that the relative abundance of sequences assigned to Bacteroides spp. were significantly higher in SLE samples (P < 0.02). In the control group, we also found a significantly higher (P < 0.05) relative abundance of sequences tentatively assigned to Desulfovibrio, the most common genus of sulfate-reducing bacteria in the human gut (45). It is worth mentioning that some authors have described how dysbiosis could affect mucosal barrier function and impair immunoregulatory mechanisms, leading to pathological effects in systemic immunity. Using animal models, a direct involvement of components of the microbiota in chronic intestinal inflammation (46) and the protective role of specific commensals in avoiding bacterial translocation (47) have been demonstrated. Finally, we should bear in mind that nonsteroidal anti-inflammatory drugs and the antimalarial treatment of SLE patients could have an influence on the observed dysbiosis.

Figure S2

Interindividual variation in fecal samples of SLE subjects (a) and control subjects (b) in the proportion of the major microbial phyla/families. The aggregate microbiota composition of control (HC codes) and lupus-affected (SLE codes) subjects at the levels of phylum and family are indicated at the bottom of each inset. Only taxonomic groups representing >1% are shown in the keys. Download Figure S2, DOC file, 0.1 MB.
Copyright © 2014 Hevia et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution-Noncommercial-ShareAlike 3.0 Unported license, which permits unrestricted noncommercial use, distribution, and reproduction in any medium, provided the original author and source are credited.

FIG 3 
  • Open in new tab
  • Download powerpoint
FIG 3 

PCoA results for the 16S rRNA profiles at the phylum (A) and family (B) level. The presence or absence of SLE was further included as metadata. HC, closed circles; SLE, open triangles.

In our study, we also determined the tentative metagenomes from phylogenetically associated reference genomes. Our aim was to highlight metabolic pathways and shifts associated with the SLE population compared with controls. We inferred the functionality of the different putative metagenomes by using PICRUSt software, which allows the prediction of metabolic pathways from the 16S rRNA reads (48). A functional analysis using the data obtained from the KEGG pathways at level 3 allowed us to detect certain processes potentially associated with either healthy controls or SLE. KEGG levels are the different hierarchical subdivisions in which the functions of a biological system (cell, organism, or ecosystem) are arranged according to information organized in the KEGG database (http://www.genome.jp/kegg/kegg1a.html). Pathways displaying a difference in mean proportions between healthy controls and SLE groups of at least 0.1% are represented in Fig. S6 in the supplemental material. Some glycan degradation pathways are slightly overrepresented in the microbiota of SLE patients, probably due to the higher abundance of Bacteroidetes in these samples. Bacteroidetes, and specifically the main genus of this phylum, Bacteroides, have been shown to display broad glycan-degrading abilities (49). The same occurs with lipopolysaccharide biosynthesis proteins, which is in direct relation to the higher abundance of Bacteroidetes in the SLE samples. Remarkably, oxidative phosphorylation processes seem to be associated with SLE patients. This finding indicates that some bacteria able to perform oxidative phosphorylation may be better adapted to the intestinal ecosystem of individuals with SLE, and this could be related to the imbalance in the oxidative stress environment at the intestinal level that has been linked to some autoimmune diseases (48).

Figure S6

KEGG pathways (hierarchical level 3) for which normalized abundances differed between SLE patients and healthy controls. Functional metagenomes were predicted using PICRUSt scripts, and only pathways reaching statistically significance are shown. On the y axis, the KEGG pathways are indicated. Mean proportions for each metadata group (healthy controls versus SLE) are shown together with the confidence intervals of the data. Variations due to every single metagenome on mean proportions are shown by horizontal error bars. Download Figure S6, DOCX file, 0.2 MB.
Copyright © 2014 Hevia et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution-Noncommercial-ShareAlike 3.0 Unported license, which permits unrestricted noncommercial use, distribution, and reproduction in any medium, provided the original author and source are credited.

In summary, understanding and potentially manipulating immune responses through the action of intestinal microbiota comprise one of the most active fields in probiotic and prebiotic research. Most likely, the dysbiosis defined in this work is the consequence of the altered immune function of SLE patients. At present, the treatment of SLE patients is exclusively performed with drugs, and our findings could indicate that challenging the immune system with the bacteria depleted in SLE, or substances promoting their growth, could influence SLE physiology. To the best of our knowledge, experimental papers about the relationship of the human gut microbiome and SLE have not been published, and this is the first report describing an SLE-associated intestinal dysbiosis. Thus, our results establish the basis to delve deeper into the understanding of the relationship between the human gut microbiota and autoimmune diseases.

MATERIALS AND METHODS

Ethics statement.Ethics approval for this study (reference code AGL2010-14952; grant title “Towards a better understanding of gut microbiota functionality in some immune disorders”) was obtained from the Bioethics Committee of CSIC (Consejo Superior de Investigaciones Científicas) and from the Regional Ethics Committee for Clinical Research (Servicio de Salud del Principado de Asturias) in compliance with the Declaration of Helsinki. All determinations were performed with fully informed written consent from all participants involved in the study.

Study subjects.The study sample comprised 20 patients with SLE (SLE codes) and 20 healthy controls (HC codes). SLE patients were selected from the updated Asturian Register of Lupus (Asociación Lúpicos de Asturias, Oviedo, Spain). All of them fulfilled at least four of the American College of Rheumatology criteria for SLE (50). The individuals recruited for the SLE group were in remission and had not had any immunotherapy or corticoid treatment during the previous months, since those are the treatments that could have the strongest influence on the physiology of the patients. All patients were women of Caucasian origin, 49.2 ± 10.7 years old (mean ± standard deviation [SD]), and had no active disease at the time of sampling (SLEDAI score of ≤8). Information on clinical manifestations was obtained by reviewing clinical records (Table 1). Patients were also asked precise questions regarding treatment received during the previous 6 months. Only those individuals who had not used antibiotics, glucocorticoids, immunosuppressive drugs, monoclonal antibodies, or other immunotherapies were recruited for the study. Eighteen patients were receiving antimalarial treatment, and all of them were regular consumers of nonsteroidal anti-inflammatory drugs. Twenty age-matched healthy women (46.9 ± 8.6 years old) from the same population were recruited as controls.

(i) Nutritional assessment.Variables of macro- and micronutrient intakes were collected by means of a semiquantitative food frequency questionnaire that included 160 items. During a personal interview, carried out within 7 days after the collection of the fecal sample, subjects were asked item by item whether they usually ate the food and, if so, how much they ate. For this purpose, three different serving sizes of each cooked food were presented in pictures to the participants so that they could choose from up to 7 serving sizes (from “less than the small one” to “more than the large one”). For some of the foods consumed, amounts were recorded in household units, by volume, or by measuring with a ruler. Special attention was paid to cooking practices and number and amount of ingredients used in each recipe, as well as questions concerning menu preparation (e.g., type of oil or milk used).

Food intake was analyzed for energy and macro- and micronutrient contents by using the nutrient food composition tables developed by CESNID (51). Total, soluble, and insoluble fiber intake was completed based on Marlett food composition tables (52), and polyphenol content was calculated from the U.S. Department of Agriculture (USDA) Nutrient Database (53).

(ii) Anthropometric measures.The body mass index (BMI) was calculated from the following formula: weight/(height)2 (in kg/m2). Height was registered by using a stadiometer with an accuracy of ±1 mm (Año-Sayol, Barcelona, Spain). Subjects were barefoot, in an upright position, and with the head positioned in the Frankfort horizontal plane. Weight was measured on a scale with an accuracy of ±100 g (Seca, Hamburg, Germany).

(iii) Lifestyle-related factors.During the interview, other factors associated with the lifestyle of the subject were registered. Smoking habit, physical activity, alcohol consumption, and supplements use were included in the questionnaire. Smoking status was categorized as nonsmoker (including exsmokers and occasional smokers) or current smoker. Those subjects who reported that they never exercised were categorized as physically inactive. “Regular alcohol consumer” refers to those subjects who declared a regular consumption of alcoholic drinks. Also, the use of vitamin and mineral supplements during the last month was self-reported.

Fecal sample collection and DNA extraction.Fresh fecal material (between 10 and 50 g per person) was collected in a sterile container and immediately manipulated and homogenized within a maximum of 3 h from defecation. During the waiting period, from defecation to homogenization, samples were kept at 4°C. Thirty milliliters of RNAlater solution (Applied Biosystems, Foster City, CA) was added to 10 g of sample, and the mixture was homogenized in a sterile bag, using a stomacher apparatus (IUL Instruments, Barcelona, Spain) with three cycles at high speed, 1 min per cycle. Homogenized samples were then stored at −80°C until use.

For DNA extraction, samples were thawed and the QIAamp DNA stool minikit (Qiagen Ltd., Strasse, Germany) was used as previously described (35).

16S rRNA gene amplification.Partial 16S rRNA gene sequences were amplified from extracted DNA by using the primer pair Probio_Uni (5′-CCTACGGGRSGCAGCAG-3')/Probio_Rev (5′-ATTACCGCGGCTGCT-3′) (35), which targets the V3 region of the 16S rRNA gene sequence. The PCR conditions used were 5 min at 95°C, 35 cycles of 30 s at 94°C, 30 s at 55°C, and 90 s at 72°C, followed by 10 min at 72°C. Amplification was carried out using a Verity thermocycler (Applied Biosystems). The integrity of the PCR amplicons was analyzed by electrophoresis on an Experion workstation (Bio-Rad, Hertfordshire, United Kingdom).

Ion Torrent PGM sequencing of 16S rRNA gene-based amplicons.The PCR products derived from amplification of specific 16S rRNA gene hypervariable regions were purified by electrophoretic separation on a 1.5% agarose gel and the use of a Wizard SV Gen PCR cleanup system (Promega, Madison, WI), followed by a further purification step involving Agencourt AMPure XP DNA purification beads (Beckman Coulter Genomics GmbH, Bernried, Germany) in order to remove primer dimers. The DNA concentration of the amplified sequence library was estimated through use of the Experion system (Bio-Rad). From the concentration and the average size of each amplicon library, the amount of DNA fragments per microliter was calculated and libraries for each run were diluted to 3E9 DNA molecules prior to clonal amplification. Emulsion PCR was carried out using the Ion OneTouch 200 template kit v2 DL (Life Technologies, Guilford, CA) according to the manufacturer’s instructions. Sequencing of the amplicon libraries was carried out on 316 chips by using the Ion Torrent PGM system and employing the Ion Sequencing 200 kit (Life Technologies) according to the supplier’s instructions at the DNA sequencing facility, GenProbio s.r.l. After sequencing, the individual sequence reads were filtered with the PGM software to remove low-quality and polyclonal sequences. Sequences matching the PGM 3′ adaptor were also automatically trimmed. All PGM quality-approved, trimmed, and filtered data were exported as SFF files.

Sequence-based microbiota analysis.The SFF files were processed using QIIME (54). Quality control retained sequences had lengths between 150 and 200 bp, a mean sequence quality score of >25, and with truncation of a sequence at the first base if a low-quality rolling 10-bp window was found. Presence of homopolymers of >7 bp and sequences with mismatched primers were omitted. In order to calculate downstream diversity measures (alpha and beta diversity indices, Unifrac analysis), 16S rRNA OTUs were defined at ≥97% sequence homology by using uclust (55). Chimeric sequences were removed using ChimeraSlayer (56). Furthermore, OTUs that included less than 10 sequences were filtered using QIIME (54). All reads were classified to the lowest possible taxonomic rank by using QIIME and a reference data set from the Ribosomal Database Project (57). The sequence data features of all the samples are included in Table S1 in the supplemental material.

Different alpha diversity indexes were calculated using QIIME and information from the OTU tables using the alpha_diversity.py script. The different diversity metrics were set by passing the option –s to the script. The following indexes were calculated for every sample and compared between groups by using a two-sided Student’s t test: Chao1, PD whole tree, observed species, Shannon, and Simpson.

Analysis by qPCR.Quantification of total fecal bacteria, Firmicutes, and Bacteroidetes by qPCR was performed by using previously described primers and conditions (58–60). Analyses were done in duplicate in a 7500 Fast real-time PCR system (Applied Biosystems) using Sybr green PCR master mix (Applied Biosystems). Standard curves were made with pure cultures, grown under anaerobic conditions at 37°C, of Escherichia coli LMG 2092 in Gifu anaerobic medium (GAM; Nissui Pharmaceutical Co., Tokyo, Japan), Faecalibacterium prausnitzii DMSZ 17677 in RCM formula (Oxoid Ltd., Basingstoke, Hampshire, United Kingdom), and Bacteroides thetaiotaomicron DSMZ 2079 in GAM.

Functional inference.The functionality of the different metagenomes, grouped by disease status (healthy control versus SLE), was predicted using the software PICRUSt 0.9.1 (http://picrust.github.io), which has been explained in detail elsewhere (48). In short, this software allows the prediction of functional pathways from the 16S rRNA reads. First, a collection of closed-reference OTUs was obtained from the filtered reads by using QIIME v 1.7.0 (54) and by querying the data against the IMG/GG reference collection (GreenGenes database, May 2013 version; http://greengenes.lbl.gov). Reverse-strand matching was enabled during the query, and OTUs were picked at a 97% identity. A BIOM-formatted table (biological observation matrix [61]) was obtained with the pick_closed_reference_otus.py script. This table, containing the relative abundances of the different reference OTUs in all the metagenomes, was normalized based on the predicted 16S rRNA copy number by using the script normalize_by_copy_number.py. Final functional predictions, inferred from the metagenomes, were created with the script predict_metagenomes.py. When necessary, tab-delimited tables were obtained with the script convert_biom.py.

Analysis of predicted metagenomes.PICRUSt and QIIME provide a number of scripts that can be useful for analyzing both 16S rRNA gene relative abundances and the predicted metabolic data. Predicted metagenomic contents were collapsed at KEGG pathway level 3 (http://www.genome.jp/kegg/pathway.html) with the categorize_by_function.py script, and the data were analyzed statistically by using STAMP v 2.0.0 (62). STAMP allows data filtering and the application of different statistical tests and corrections, including PCoA. It also generates different graphics, including box plots, error plots, and scatter plots. Data of the KEGG pathway distributions were plotted by using the script summarize_taxa_through_plots.py. Associations of different taxonomic categories to SLE were statistically analyzed with the script otu_category_significance.py.

Statistical analyses.Statistical analysis was performed using IBM-SPSS version 19.0 (IBM SPSS, Inc., Chicago, IL). For descriptive purposes, in Table 2 the mean values are presented as means ± SD on untransformed variables. Differences between SLE patients and controls were compared by using a multivariate linear model, including energy intake and BMI as covariates.

Individuals were ordered according to their sequence data composition by principal component analysis using the taxonomic data at the phylum and family levels. Patterns were extracted using all the variations from the taxonomic data via an indirect method as a model and SLE as metadata. To analyze the associations between inferred metabolic pathways and SLE, metabolic pathways with very low abundance levels (<0.001 in 50% of the samples) were excluded from all analyses. Association of KEGG pathways to SLE were identified by running two-sided Welch tests on every pair of means. This test is a variation of Student’s t test and is used when equal variance cannot be assumed in both groups. Confidence intervals (95%) were obtained by inverting the Welch’s tests. The false discovery rate (FDR) correction (48, 63) was finally applied in all cases, and significant differences between healthy controls and SLE patients were only considered when below a P value of 0.05 or a q value below 0.2. P and q values at the phylum and family levels are included in Table S2 in the supplemental material. In the particular case of the family Desulfovibrionaceae (with P < 0.05), further statistical analysis was carried out for sequences tentatively assigned to the genus Desulfovibirio by using a one-sided t test.

In relation to qPCR results, not all the bacterial groups showed normal distribution; therefore, differences in bacterial levels between groups of individuals were analyzed using a nonparametric test (Mann-Whitney U test).

Nucleotide sequence accession number.The raw sequences reported in this article have been deposited in the NCBI Short Read Archive (SRA; study accession number SRP028162).

ACKNOWLEDGMENTS

This study was financed by European Union FEDER funds and the Spanish Plan Nacional de I+D (grant AGL2010-14952). Arancha Hevia was the recipient of an FPI grant, and Borja Sánchez was the recipient of a Juan de la Cierva postdoctoral contract, both from the Spanish Ministerio de Ciencia e Innovación.

FOOTNOTES

    • Received 25 June 2014
    • Accepted 29 August 2014
    • Published 30 September 2014
  • Copyright © 2014 Hevia et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution-Noncommercial-ShareAlike 3.0 Unported license, which permits unrestricted noncommercial use, distribution, and reproduction in any medium, provided the original author and source are credited.

REFERENCES

  1. 1.↵
    1. Fraher MH,
    2. O’Toole PW,
    3. Quigley EM
    . 2012. Techniques used to characterize the gut microbiota: a guide for the clinician. Nat. Rev. Gastroenterol. Hepatol. 9:312–322. doi:10.1038/nrgastro.2012.44.
    OpenUrlCrossRefPubMed
  2. 2.↵
    1. Weinstock GM
    . 2012. Genomic approaches to studying the human microbiota. Nature 489:250–256. doi:10.1038/nature11553.
    OpenUrlCrossRefPubMedWeb of Science
  3. 3.↵
    1. Gill SR,
    2. Pop M,
    3. Deboy RT,
    4. Eckburg PB,
    5. Turnbaugh PJ,
    6. Samuel BS,
    7. Gordon JI,
    8. Relman DA,
    9. Fraser-Liggett CM,
    10. Nelson KE
    . 2006. Metagenomic analysis of the human distal gut microbiome. Science 312:1355–1359. doi:10.1126/science.1124234.
    OpenUrlAbstract/FREE Full Text
  4. 4.↵
    1. Xu J,
    2. Mahowald MA,
    3. Ley RE,
    4. Lozupone CA,
    5. Hamady M,
    6. Martens EC,
    7. Henrissat B,
    8. Coutinho PM,
    9. Minx P,
    10. Latreille P,
    11. Cordum H,
    12. Van Brunt A,
    13. Kim K,
    14. Fulton RS,
    15. Fulton LA,
    16. Clifton SW,
    17. Wilson RK,
    18. Knight RD,
    19. Gordon JI
    . 2007. Evolution of symbiotic bacteria in the distal human intestine. PLoS Biol. 5:e156. doi:10.1371/journal.pbio.0050156.
    OpenUrlCrossRefPubMed
  5. 5.↵
    1. Bäckhed F,
    2. Ley RE,
    3. Sonnenburg JL,
    4. Peterson DA,
    5. Gordon JI
    . 2005. Host-bacterial mutualism in the human intestine. Science 307:1915–1920. doi:10.1126/science.1104816.
    OpenUrlAbstract/FREE Full Text
  6. 6.↵
    1. Holmes E,
    2. Kinross J,
    3. Gibson GR,
    4. Burcelin R,
    5. Jia W,
    6. Pettersson S,
    7. Nicholson JK
    . 2012. Therapeutic modulation of microbiota-host metabolic interactions. Sci. Transl. Med. 4:137rv6. doi:10.1126/scitranslmed.3004244.
    OpenUrlFREE Full Text
  7. 7.↵
    1. Collins SM,
    2. Surette M,
    3. Bercik P
    . 2012. The interplay between the intestinal microbiota and the brain. Nat. Rev. Microbiol. 10:735–742. doi:10.1038/nrmicro2876.
    OpenUrlCrossRefPubMed
  8. 8.↵
    1. Koren O,
    2. Goodrich JK,
    3. Cullender TC,
    4. Spor A,
    5. Laitinen K,
    6. Bäckhed HK,
    7. Gonzalez A,
    8. Werner JJ,
    9. Angenent LT,
    10. Knight R,
    11. Bäckhed F,
    12. Isolauri E,
    13. Salminen S,
    14. Ley RE
    . 2012. Host remodeling of the gut microbiome and metabolic changes during pregnancy. Cell 150:470–480. doi:10.1016/j.cell.2012.07.008.
    OpenUrlCrossRefPubMedWeb of Science
  9. 9.↵
    1. Claesson MJ,
    2. Jeffery IB,
    3. Conde S,
    4. Power SE,
    5. O’Connor EM,
    6. Cusack S,
    7. Harris HM,
    8. Coakley M,
    9. Lakshminarayanan B,
    10. O’Sullivan O,
    11. Fitzgerald GF,
    12. Deane J,
    13. O’Connor M,
    14. Harnedy N,
    15. O’Connor K,
    16. O’Mahony D,
    17. van Sinderen D,
    18. Wallace M,
    19. Brennan L,
    20. Stanton C,
    21. Marchesi JR,
    22. Fitzgerald AP,
    23. Shanahan F,
    24. Hill C,
    25. Ross RP,
    26. O’Toole PW
    . 2012. Gut microbiota composition correlates with diet and health in the elderly. Nature 488:178–184. doi:10.1038/nature11319.
    OpenUrlCrossRefPubMedWeb of Science
  10. 10.↵
    1. Sharon I,
    2. Morowitz MJ,
    3. Thomas BC,
    4. Costello EK,
    5. Relman DA,
    6. Banfield JF
    . 2013. Time series community genomics analysis reveals rapid shifts in bacterial species, strains, and phage during infant gut colonization. Genome Res. 23:111–120. doi:10.1101/gr.142315.112.
    OpenUrlAbstract/FREE Full Text
  11. 11.↵
    1. Turroni F,
    2. Peano C,
    3. Pass DA,
    4. Foroni E,
    5. Severgnini M,
    6. Claesson MJ,
    7. Kerr C,
    8. Hourihane J,
    9. Murray D,
    10. Fuligni F,
    11. Gueimonde M,
    12. Margolles A,
    13. De Bellis G,
    14. O’Toole PW,
    15. van Sinderen D,
    16. Marchesi JR,
    17. Ventura M
    . 2012. Diversity of bifidobacteria within the infant gut microbiota. PLoS One 7:e36957. doi:10.1371/journal.pone.0036957.
    OpenUrlCrossRefPubMed
  12. 12.↵
    1. Yatsunenko T,
    2. Rey FE,
    3. Manary MJ,
    4. Trehan I,
    5. Dominguez-Bello MG,
    6. Contreras M,
    7. Magris M,
    8. Hidalgo G,
    9. Baldassano RN,
    10. Anokhin AP,
    11. Heath AC,
    12. Warner B,
    13. Reeder J,
    14. Kuczynski J,
    15. Caporaso JG,
    16. Lozupone CA,
    17. Lauber C,
    18. Clemente JC,
    19. Knights D,
    20. Knight R,
    21. Gordon JI
    . 2012. Human gut microbiome viewed across age and geography. Nature 486:222–227.
    OpenUrlCrossRefPubMedWeb of Science
  13. 13.↵
    1. De Filippo C,
    2. Cavalieri D,
    3. Di Paola M,
    4. Ramazzotti M,
    5. Poullet JB,
    6. Massart S,
    7. Collini S,
    8. Pieraccini G,
    9. Lionetti P
    . 2010. Impact of diet in shaping gut microbiota revealed by a comparative study in children from Europe and rural Africa. Proc. Natl. Acad. Sci. U. S. A. 107:14691–14696. doi:10.1073/pnas.1005963107.
    OpenUrlAbstract/FREE Full Text
  14. 14.↵
    1. Arumugam M,
    2. Raes J,
    3. Pelletier E,
    4. Le Paslier D,
    5. Yamada T,
    6. Mende DR,
    7. Fernandes GR,
    8. Tap J,
    9. Bruls T,
    10. Batto JM,
    11. Bertalan M,
    12. Borruel N,
    13. Casellas F,
    14. Fernandez L,
    15. Gautier L,
    16. Hansen T,
    17. Hattori M,
    18. Hayashi T,
    19. Kleerebezem M,
    20. Kurokawa K,
    21. Leclerc M,
    22. Levenez F,
    23. Manichanh C,
    24. Nielsen HB,
    25. Nielsen T,
    26. Pons N,
    27. Poulain J,
    28. Qin J,
    29. Sicheritz-Ponten T,
    30. Tims S,
    31. Torrents D,
    32. Ugarte E,
    33. Zoetendal EG,
    34. Wang J,
    35. Guarner F,
    36. Pedersen O,
    37. de Vos WM,
    38. Brunak S,
    39. Doré J,
    40. MetaHIT Consortium
    41. Antolín M,
    42. Artiguenave F,
    43. Blottiere HM,
    44. Almeida M,
    45. Brechot C,
    46. Cara C,
    47. Chervaux C,
    48. Cultrone A,
    49. Delorme C,
    50. Denariaz G,
    51. Dervyn R,
    52. Foerstner KU,
    53. Friss C,
    54. van de Guchte M,
    55. Guedon E,
    56. Haimet F,
    57. Huber W,
    58. van Hylckama-Vlieg J,
    59. Jamet A,
    60. Juste C,
    61. Kaci G,
    62. Knol J,
    63. Lakhdari O,
    64. Layec S,
    65. Le Roux K,
    66. Maguin E,
    67. Mérieux A,
    68. Melo Minardi R,
    69. M’rini C,
    70. Muller J,
    71. Oozeer R,
    72. Parkhill J,
    73. Renault P,
    74. Rescigno M,
    75. Sanchez N,
    76. Sunagawa S,
    77. Torrejon A,
    78. Turner K,
    79. Vandemeulebrouck G,
    80. Varela E,
    81. Winogradsky Y,
    82. Zeller G,
    83. Weissenbach J,
    84. Ehrlich SD,
    85. Bork P
    . 2011. Enterotypes of the human gut microbiome. Nature 473:174–180. doi:10.1038/nature09944.
    OpenUrlCrossRefPubMedWeb of Science
  15. 15.↵
    1. Wu GD,
    2. Chen J,
    3. Hoffmann C,
    4. Bittinger K,
    5. Chen YY,
    6. Keilbaugh SA,
    7. Bewtra M,
    8. Knights D,
    9. Walters WA,
    10. Knight R,
    11. Sinha R,
    12. Gilroy E,
    13. Gupta K,
    14. Baldassano R,
    15. Nessel L,
    16. Li H,
    17. Bushman FD,
    18. Lewis JD
    . 2011. Linking long-term dietary patterns with gut microbial enterotypes. Science 334:105–108. doi:10.1126/science.1208344.
    OpenUrlAbstract/FREE Full Text
  16. 16.↵
    1. Ley RE
    . 2010. Obesity and the human microbiome. Curr. Opin. Gastroenterol. 26:5–11. doi:10.1097/MOG.0b013e328333d751.
    OpenUrlCrossRefPubMedWeb of Science
  17. 17.↵
    1. Turnbaugh PJ,
    2. Ley RE,
    3. Mahowald MA,
    4. Magrini V,
    5. Mardis ER,
    6. Gordon JI
    . 2006. An obesity-associated gut microbiome with increased capacity for energy harvest. Nature 444:1027–1031. doi:10.1038/nature05414.
    OpenUrlCrossRefPubMedWeb of Science
  18. 18.↵
    1. Vijay-Kumar M,
    2. Aitken JD,
    3. Carvalho FA,
    4. Cullender TC,
    5. Mwangi S,
    6. Srinivasan S,
    7. Sitaraman SV,
    8. Knight R,
    9. Ley RE,
    10. Gewirtz AT
    . 2010. Metabolic syndrome and altered gut microbiota in mice lacking Toll-like receptor 5. Science 328:228–231. doi:10.1126/science.1179721.
    OpenUrlAbstract/FREE Full Text
  19. 19.↵
    1. Manichanh C,
    2. Borruel N,
    3. Casellas F,
    4. Guarner F
    . 2012. The gut microbiota in IBD. Nat. Rev. Gastroenterol. Hepatol. 9:599–608. doi:10.1038/nrgastro.2012.152.
    OpenUrlCrossRefPubMed
  20. 20.↵
    1. Scher JU,
    2. Sczesnak A,
    3. Longman RS,
    4. Segata N,
    5. Ubeda C,
    6. Bielski C,
    7. Rostron T,
    8. Cerundolo V,
    9. Pamer EG,
    10. Abramson SB,
    11. Huttenhower C,
    12. Littman DR
    . 2013. Expansion of intestinal Prevotella copri correlates with enhanced susceptibility to arthritis. Elife 2:e01202. doi:10.7554/eLife.01202.
    OpenUrlCrossRefPubMed
  21. 21.↵
    1. Wen L,
    2. Ley RE,
    3. Volchkov PY,
    4. Stranges PB,
    5. Avanesyan L,
    6. Stonebraker AC,
    7. Hu C,
    8. Wong FS,
    9. Szot GL,
    10. Bluestone JA,
    11. Gordon JI,
    12. Chervonsky AV
    . 2008. Innate immunity and intestinal microbiota in the development of type 1 diabetes. Nature 455:1109–1113. doi:10.1038/nature07336.
    OpenUrlCrossRefPubMedWeb of Science
  22. 22.↵
    1. Agmon-Levin N,
    2. Mosca M,
    3. Petri M,
    4. Shoenfeld Y
    . 2012. Systemic lupus erythematosus one disease or many? Autoimmun. Rev. 11:593–595. doi:10.1016/j.autrev.2011.10.020.
    OpenUrlCrossRefPubMed
  23. 23.↵
    1. López P,
    2. Mozo L,
    3. Gutiérrez C,
    4. Suárez A
    . 2003. Epidemiology of systemic lupus erythematosus in a northern Spanish population: gender and age influence on immunological features. Lupus 12:860–865. doi:10.1191/0961203303lu469xx.
    OpenUrlCrossRefPubMedWeb of Science
  24. 24.↵
    1. James JA,
    2. Robertson JM
    . 2012. Lupus and Epstein-Barr. Curr. Opin. Rheumatol. 24:383–388. doi:10.1097/BOR.0b013e3283535801.
    OpenUrlCrossRefPubMed
  25. 25.↵
    1. Francis L,
    2. Perl A
    . 2010. Infection in systemic lupus erythematosus: friend or foe? Int. J. Clin Rheumtol. 5:59–74. doi:10.2217/ijr.09.72.
    OpenUrlCrossRefPubMed
  26. 26.↵
    1. Pisetsky D,
    2. Drayton D,
    3. Wu ZQ
    . 1999. Specificity of antibodies to bacterial DNA in the sera of healthy human subjects and patients with systemic lupus erythematosus. J. Rheumatol. 26:1934–1938.
    OpenUrlPubMedWeb of Science
  27. 27.↵
    1. Gul’neva M,
    2. Romanov VA,
    3. Shilkina NP
    . 2007. Intestinal microecology in some systemic connective tissue diseases. Zh. Mikrobiol. Epidemiol. Immunobiol. 4:38–41. (In Russian.).
    OpenUrl
  28. 28.↵
    1. Proal AD,
    2. Albert PJ,
    3. Marshall T
    . 2009. Autoimmune disease in the era of the metagenome. Autoimmun. Rev. 8:677–681. doi:10.1016/j.autrev.2009.02.016.
    OpenUrlCrossRefPubMed
  29. 29.↵
    1. Berer K,
    2. Mues M,
    3. Koutrolos M,
    4. Rasbi ZA,
    5. Boziki M,
    6. Johner C,
    7. Wekerle H,
    8. Krishnamoorthy G
    . 2011. Commensal microbiota and myelin autoantigen cooperate to trigger autoimmune demyelination. Nature 479:538–541. doi:10.1038/nature10554.
    OpenUrlCrossRefPubMedWeb of Science
  30. 30.↵
    1. Frank DN
    , St. Amand AL, Feldman RA, Boedeker EC, Harpaz N, Pace NR. 2007. Molecular-phylogenetic characterization of microbial community imbalances in human inflammatory bowel diseases. Proc. Natl. Acad. Sci. U. S. A. 104:13780–13785. doi:10.1073/pnas.0706625104.
    OpenUrlAbstract/FREE Full Text
  31. 31.↵
    1. Markle JG,
    2. Frank DN,
    3. Mortin-Toth S,
    4. Robertson CE,
    5. Feazel LM,
    6. Rolle-Kampczyk U,
    7. von Bergen M,
    8. McCoy KD,
    9. Macpherson AJ,
    10. Danska JS
    . 2013. Sex differences in the gut microbiome drive hormone-dependent regulation of autoimmunity. Science 339:1084–1088. doi:10.1126/science.1233521.
    OpenUrlAbstract/FREE Full Text
  32. 32.↵
    1. Proal AD,
    2. Albert PJ,
    3. Marshall TG
    . 2013. The human microbiome and autoimmunity. Curr. Opin. Rheumatol. 25:234–240. doi:10.1097/BOR.0b013e32835cedbf.
    OpenUrlCrossRefPubMed
  33. 33.↵
    1. Moco S,
    2. Martin FP,
    3. Rezzi S
    . 2012. Metabolomics view on gut microbiome modulation by polyphenol-rich foods. J. Proteome Res. 11:4781–4790. doi:10.1021/pr300581s.
    OpenUrlCrossRefPubMed
  34. 34.↵
    1. Queipo-Ortuño MI,
    2. Boto-Ordóñez M,
    3. Murri M,
    4. Gomez-Zumaquero JM,
    5. Clemente-Postigo M,
    6. Estruch R,
    7. Cardona Diaz F,
    8. Andrés-Lacueva C,
    9. Tinahones FJ
    . 2012. Influence of red wine polyphenols and ethanol on the gut microbiota ecology and biochemical biomarkers. Am. J. Clin. Nutr. 95:1323–1334. doi:10.3945/ajcn.111.027847.
    OpenUrlAbstract/FREE Full Text
  35. 35.↵
    1. Milani C,
    2. Hevia A,
    3. Foroni E,
    4. Duranti S,
    5. Turroni F,
    6. Lugli GA,
    7. Sanchez B,
    8. Martín R,
    9. Gueimonde M,
    10. van Sinderen D,
    11. Margolles A,
    12. Ventura M
    . 2013. Assessing the fecal microbiota: an optimized Ion Torrent 16S rRNA gene-based analysis protocol. PLoS One 8:e68739. doi:10.1371/journal.pone.0068739.
    OpenUrlCrossRefPubMed
  36. 36.↵
    1. Qin J,
    2. Li R,
    3. Raes J,
    4. Arumugam M,
    5. Burgdorf KS,
    6. Manichanh C,
    7. Nielsen T,
    8. Pons N,
    9. Levenez F,
    10. Yamada T,
    11. Mende DR,
    12. Li J,
    13. Xu J,
    14. Li S,
    15. Li D,
    16. Cao J,
    17. Wang B,
    18. Liang H,
    19. Zheng H,
    20. Xie Y,
    21. Tap J,
    22. Lepage P,
    23. Bertalan M,
    24. Batto JM,
    25. Hansen T,
    26. Le Paslier D,
    27. Linneberg A,
    28. Nielsen HB,
    29. Pelletier E,
    30. Renault P,
    31. Sicheritz-Ponten T,
    32. Turner K,
    33. Zhu H,
    34. Yu C,
    35. Li S,
    36. Jian M,
    37. Zhou Y,
    38. Li Y,
    39. Zhang X,
    40. Li S,
    41. Qin N,
    42. Yang H,
    43. Wang J,
    44. Brunak S,
    45. Doré J,
    46. Guarner F,
    47. Kristiansen K,
    48. Pedersen O,
    49. Parkhill J,
    50. Weissenbach J,
    51. MetaHIT Consortium
    52. Bork P,
    53. Ehrlich SD,
    54. Wang J
    . 2010. A human gut microbial gene catalogue established by metagenomic sequencing. Nature 464:59–65. doi:10.1038/nature08821.
    OpenUrlCrossRefPubMedWeb of Science
  37. 37.↵
    1. Larsen N,
    2. Vogensen FK,
    3. van den Beg FW,
    4. Nielsen DS,
    5. Andreasen AS,
    6. Pedersen BK,
    7. Al-Soud WA,
    8. Sørensen SJ,
    9. Hansen LH,
    10. Jakobsen M
    . 2010. Gut microbiota in human adults with type 2 diabetes differs from non-diabetic adults. PLoS One 5:e9085. doi:10.1371/journal.pone.0009085.
    OpenUrlCrossRefPubMed
  38. 38.↵
    1. Man SM,
    2. Kaakoush NO,
    3. Mitchell HM
    . 2011. The role of bacteria and pattern-recognition receptors in Crohn’s disease. Nat. Rev. Gastroenterol. Hepatol. 8:152–168. doi:10.1038/nrgastro.2011.3.
    OpenUrlCrossRefPubMed
  39. 39.↵
    1. Martínez I,
    2. Lattimer JM,
    3. Hubach KL,
    4. Case JA,
    5. Yang J,
    6. Weber CG,
    7. Louk JA,
    8. Rose DJ,
    9. Kyureghian G,
    10. Peterson DA,
    11. Haub MD,
    12. Walter J
    . 2013. Gut microbiome composition is linked to whole grain-induced immunological improvements. ISME J. 7:269–280. doi:10.1038/ismej.2012.104.
    OpenUrlCrossRefPubMed
  40. 40.↵
    1. Hammer Ø,
    2. Harper DAT,
    3. Ryan PD
    . 2001. Past: paleontological statistics software package for education and data analysis. Palaeontol. Electron. 4:1–9.
    OpenUrl
  41. 41.↵
    1. Antharam VC,
    2. Li EC,
    3. Ishmael A,
    4. Sharma A,
    5. Mai V,
    6. Rand KH,
    7. Wang GP
    . 2013. Intestinal dysbiosis and depletion of butyrogenic bacteria in Clostridium difficile infection and nosocomial diarrhea. J. Clin. Microbiol. 51:2884–2992. doi:10.1128/JCM.00845-13.
    OpenUrlAbstract/FREE Full Text
  42. 42.↵
    1. Berry D,
    2. Reinisch W
    . 2013. Intestinal microbiota: a source of novel biomarkers in inflammatory bowel diseases? Best Pract. Res. Clin. Gastroenterol. 27:47–58. doi:10.1016/j.bpg.2013.03.005.
    OpenUrlCrossRef
  43. 43.↵
    1. Tannock GW,
    2. Lawley B,
    3. Munro K,
    4. Lay C,
    5. Taylor C,
    6. Daynes C,
    7. Baladjay L,
    8. Mcleod R,
    9. Thompson-Fawcett M
    . 2012. Comprehensive analysis of the bacterial content of stool from patients with chronic pouchitis, normal pouches, or familial adenomatous polyposis pouches. Inflamm. Bowel Dis. 18:925–934. doi:10.1002/ibd.21936.
    OpenUrlCrossRefPubMedWeb of Science
  44. 44.↵
    1. Ye J,
    2. Lee JW,
    3. Presley LL,
    4. Bent E,
    5. Wei B,
    6. Braun J,
    7. Schiller NL,
    8. Straus DS,
    9. Borneman J
    . 2008. Bacteria and bacterial rRNA genes associated with the development of colitis in IL-10(-/-) mice. Inflamm. Bowel Dis. 14:1041–1050. doi:10.1002/ibd.20442.
    OpenUrlCrossRefPubMedWeb of Science
  45. 45.↵
    1. Rey FE,
    2. Gonzalez MD,
    3. Cheng J,
    4. Wu M,
    5. Ahern PP,
    6. Gordon JI
    . 2013. Metabolic niche of a prominent sulfate-reducing human gut bacterium. Proc. Natl. Acad. Sci. U. S. A. 110:13582–13587. doi:10.1073/pnas.1312524110.
    OpenUrlAbstract/FREE Full Text
  46. 46.↵
    1. Tlaskalová-Hogenová H,
    2. Stěpánková R,
    3. Kozáková H,
    4. Hudcovic T,
    5. Vannucci L,
    6. Tučková L,
    7. Rossmann P,
    8. Hrnčíř T,
    9. Kverka M,
    10. Zákostelská Z,
    11. Klimešová K,
    12. Přibylová J,
    13. Bártová J,
    14. Sanchez D,
    15. Fundová P,
    16. Borovská D,
    17. Srůtková D,
    18. Zídek Z,
    19. Schwarzer M,
    20. Drastich P,
    21. Funda DP
    . 2011. The role of gut microbiota (commensal bacteria) and the mucosal barrier in the pathogenesis of inflammatory and autoimmune diseases and cancer: contribution of germ-free and gnotobiotic animal models of human diseases. Cell. Mol. Immunol. 8:110–120. doi:10.1038/cmi.2010.67.
    OpenUrlCrossRefPubMed
  47. 47.↵
    1. Othman M,
    2. Agüero R,
    3. Lin HC
    . 2008. Alterations in intestinal microbial flora and human disease. Curr. Opin. Gastroenterol. 24:11–16. doi:10.1097/MOG.0b013e3282f2b0d7.
    OpenUrlCrossRefPubMedWeb of Science
  48. 48.↵
    1. Morgan XC,
    2. Tickle TL,
    3. Sokol H,
    4. Gevers D,
    5. Devaney KL,
    6. Ward DV,
    7. Reyes JA,
    8. Shah SA,
    9. LeLeiko N,
    10. Snapper SB,
    11. Bousvaros A,
    12. Korzenik J,
    13. Sands BE,
    14. Xavier RJ,
    15. Huttenhower C
    . 2012. Dysfunction of the intestinal microbiome in inflammatory bowel disease and treatment. Genome Biol. 13:R79. doi:10.1186/gb-2012-13-9-r79.
    OpenUrlCrossRefPubMed
  49. 49.↵
    1. Martens EC,
    2. Lowe EC,
    3. Chiang H,
    4. Pudlo NA,
    5. Wu M,
    6. McNulty NP,
    7. Abbott DW,
    8. Henrissat B,
    9. Gilbert HJ,
    10. Bolam DN,
    11. Gordon JI
    . 2011. Recognition and degradation of plant cell wall polysaccharides by two human gut symbionts. PLoS Biol. 9:e1001221. doi:10.1371/journal.pbio.1001221.
    OpenUrlCrossRefPubMed
  50. 50.↵
    1. Tan EM,
    2. Cohen AS,
    3. Fries JF,
    4. Masi AT,
    5. McShane DJ,
    6. Rothfield NF,
    7. Schaller JG,
    8. Talal N,
    9. Winchester RJ
    . 1982. The 1982 revised criteria for the classification of systemic lupus erythematosus. Arthritis Rheum. 25:1271–1277. doi:10.1002/art.1780251101.
    OpenUrlCrossRefPubMedWeb of Science
  51. 51.↵
    1. Palma I,
    2. Farran A,
    3. Cantos D,
    4. CESNID
    (Centre d’Ensenyament Superior de Nutrició i Dietètica). 2008. Tablas de composición de alimentos por medidas caseras de consumo habitual en España . McGraw-Hill/Interamericana de España, Barcelona, Spain.
  52. 52.↵
    1. Marlett JA,
    2. Cheung TF
    . 1997. Database and quick methods of assessing typical dietary fiber intakes using data for 228 commonly consumed foods. J. Am. Diet. Assoc. 97:1139–1148. doi:10.1016/S0002-8223(97)00275-7.
    OpenUrlCrossRefPubMedWeb of Science
  53. 53.↵
    1. U.S. Department of Agriculture Agricultural Research Service
    . 2011. Database for the flavonoid content of select foods. USDA Agricultural Research Service, Washington, DC. http:/www.nal.usda.gov/fnic/foodcomp.
  54. 54.↵
    1. Caporaso JG,
    2. Kuczynski J,
    3. Stombaugh J,
    4. Bittinger K,
    5. Bushman FD,
    6. Costello EK,
    7. Fierer N,
    8. Peña AG,
    9. Goodrich JK,
    10. Gordon JI,
    11. Huttley GA,
    12. Kelley ST,
    13. Knights D,
    14. Koenig JE,
    15. Ley RE,
    16. Lozupone CA,
    17. McDonald D,
    18. Muegge BD,
    19. Pirrung M,
    20. Reeder J,
    21. Sevinsky JR,
    22. Turnbaugh PJ,
    23. Walters WA,
    24. Widmann J,
    25. Yatsunenko T,
    26. Zaneveld J,
    27. Knight R
    . 2010. QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7:335–336. doi:10.1038/nmeth.f.303.
    OpenUrlCrossRefPubMedWeb of Science
  55. 55.↵
    1. Edgar RC
    . 2010. Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26:2460–2461. doi:10.1093/bioinformatics/btq461.
    OpenUrlCrossRefPubMedWeb of Science
  56. 56.↵
    1. Haas BJ,
    2. Gevers D,
    3. Earl AM,
    4. Feldgarden M,
    5. Ward DV,
    6. Giannoukos G,
    7. Ciulla D,
    8. Tabbaa D,
    9. Highlander SK,
    10. Sodergren E,
    11. Methé B,
    12. DeSantis TZ,
    13. Microbiome Human
    Consortium, Petrosino JF, Knight R, Birren BW. 2011. Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons. Genome Res. 21:494–504. doi:10.1101/gr.112730.110. PubMed.
  57. 57.↵
    1. Cole JR,
    2. Wang Q,
    3. Cardenas E,
    4. Fish J,
    5. Chai B,
    6. Farris RJ,
    7. Kulam-Syed-Mohideen AS,
    8. McGarrell DM,
    9. Marsh T,
    10. Garrity GM,
    11. Tiedje JM
    . 2009. The Ribosomal Database Project: improved alignments and new tools for rRNA analysis. Nucleic Acids Res. 37 (Database issue):D141–D145. doi:10.1093/nar/gkn879.
    OpenUrlCrossRefPubMedWeb of Science
  58. 58.↵
    1. Arboleya S,
    2. Binetti A,
    3. Salazar N,
    4. Fernández N,
    5. Solís G,
    6. Hernández-Barranco A,
    7. Margolles A,
    8. de los Reyes-Gavilán CG,
    9. Gueimonde M
    . 2012. Establishment and development of intestinal microbiota in preterm neonates. FEMS Microbiol. Ecol. 79:763–772. doi:10.1111/j.1574-6941.2011.01261.x.
    OpenUrlCrossRefPubMedWeb of Science
  59. 59.↵
    1. Denman SE,
    2. McSweeney CS
    . 2006. Development of a real-time PCR assay for monitoring anaerobic fungal and cellulolytic bacterial populations within the rumen. FEMS Microbiol. Ecol. 58:572–582. doi:10.1111/j.1574-6941.2006.00190.x.
    OpenUrlCrossRefPubMedWeb of Science
  60. 60.↵
    1. Guo X,
    2. Xia X,
    3. Tang R,
    4. Zhou J,
    5. Zhao H,
    6. Wang K
    . 2008. Development of a real-time PCR method for Firmicutes and Bacteroidetes in faeces and its application to quantify intestinal population of obese and lean pigs. Lett. Appl. Microbiol. 47:367–373. doi:10.1111/j.1472-765X.2008.02408.x.
    OpenUrlCrossRefPubMedWeb of Science
  61. 61.
    1. McDonald D,
    2. Clemente JC,
    3. Kuczynski J,
    4. Rideout JR,
    5. Stombaugh J,
    6. Wendel D,
    7. Wilke A,
    8. Huse S,
    9. Hufnagle J,
    10. Meyer F,
    11. Knight R,
    12. Caporaso JG
    . 2012. The Biological Observation Matrix (BIOM) format or: how I learned to stop worrying and love the ome-ome. Gigacience 1:7. doi:10.1186/2047-217X-1-7.
    OpenUrlCrossRefPubMed
  62. 62.↵
    1. Parks DH,
    2. Beiko RG
    . 2010. Identifying biologically relevant differences between metagenomic communities. Bioinformatics 26:715–721. doi:10.1093/bioinformatics/btq041.
    OpenUrlCrossRefPubMedWeb of Science
  63. 63.↵
    1. Benjamini Y,
    2. Hochberg Y
    . 1995. Controlling the false discovery rate—a practical and powerful approach to multiple testing. J. R. Stat. Soc. B Stat. Methodol. 57:289–300.
    OpenUrl
PreviousNext
Back to top
Download PDF
Citation Tools
Intestinal Dysbiosis Associated with Systemic Lupus Erythematosus
Arancha Hevia, Christian Milani, Patricia López, Adriana Cuervo, Silvia Arboleya, Sabrina Duranti, Francesca Turroni, Sonia González, Ana Suárez, Miguel Gueimonde, Marco Ventura, Borja Sánchez, Abelardo Margolles
mBio Sep 2014, 5 (5) e01548-14; DOI: 10.1128/mBio.01548-14

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print

Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this mBio article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Intestinal Dysbiosis Associated with Systemic Lupus Erythematosus
(Your Name) has forwarded a page to you from mBio
(Your Name) thought you would be interested in this article in mBio.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Intestinal Dysbiosis Associated with Systemic Lupus Erythematosus
Arancha Hevia, Christian Milani, Patricia López, Adriana Cuervo, Silvia Arboleya, Sabrina Duranti, Francesca Turroni, Sonia González, Ana Suárez, Miguel Gueimonde, Marco Ventura, Borja Sánchez, Abelardo Margolles
mBio Sep 2014, 5 (5) e01548-14; DOI: 10.1128/mBio.01548-14
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS AND DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

Related Articles

Cited By...

About

  • About mBio
  • Editor in Chief
  • Board of Editors
  • AAM Fellows
  • Policies
  • For Reviewers
  • For the Media
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Author Warranty
  • Article Types
  • Ethics
  • Contact Us

Follow #mBio

@ASMicrobiology

       

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

 

American Society for Microbiology
1752 N St. NW
Washington, DC 20036
Phone: (202) 737-3600

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

Online ISSN: 2150-7511