Dynamic changes in short- and long-term bacterial composition following fecal microbiota transplantation for recurrent Clostridium difficile infection
- Alexa Weingarden†1, 7,
- Antonio González2,
- Yoshiki Vázquez-Baeza†3,
- Sophie Weiss4,
- Gregory Humphry5,
- Donna Berg-Lyons5,
- Dan Knights6, 7,
- Tatsuya Unno7,
- Aleh Bobr8,
- Johnthomas Kang7, 8,
- Alexander Khoruts7, 8,
- Rob Knight2, 9, 10 and
- Michael J Sadowsky1, 7Email author
© Weingarden et al.; licensee BioMed Central. 2015
Received: 5 November 2014
Accepted: 29 January 2015
Published: 30 March 2015
Fecal microbiota transplantation (FMT) is an effective treatment for recurrent Clostridium difficile infection (CDI) that often fails standard antibiotic therapy. Despite its widespread recent use, however, little is known about the stability of the fecal microbiota following FMT.
Here we report on short- and long-term changes and provide kinetic visualization of fecal microbiota composition in patients with multiply recurrent CDI that were refractory to antibiotic therapy and treated using FMT. Fecal samples were collected from four patients before and up to 151 days after FMT, with daily collections until 28 days and weekly collections until 84 days post-FMT. The composition of fecal bacteria was characterized using high throughput 16S rRNA gene sequence analysis, compared to microbiota across body sites in the Human Microbiome Project (HMP) database, and visualized in a movie-like, kinetic format. FMT resulted in rapid normalization of bacterial fecal sample composition from a markedly dysbiotic state to one representative of normal fecal microbiota. While the microbiome appeared most similar to the donor implant material 1 day post-FMT, the composition diverged variably at later time points. The donor microbiota composition also varied over time. However, both post-FMT and donor samples remained within the larger cloud of fecal microbiota characterized as healthy by the HMP.
Dynamic behavior is an intrinsic property of normal fecal microbiota and should be accounted for in comparing microbial communities among normal individuals and those with disease states. This also suggests that more frequent sample analyses are needed in order to properly assess success of FMT procedures.
KeywordsShort- and long-term changes in microbiota following FMT
Fecal microbiota transplantation (FMT) has emerged in recent years as a highly effective treatment for refractory Clostridium difficile infection (CDI) that cannot be cured with antibiotics alone . The procedure leads to prompt engraftment of donor microbiota, attainment of donor-like bacterial diversity, and normalization of the overall microbial community structure [2-8]. However, existing data characterizing long-term stability of engrafted microbiota are limited. One recent study suggests the microbiota of patients after FMT may not fully recover until 16 weeks after the procedure . This type of analysis, however, is complicated by the fact that the microbial communities are intrinsically dynamic and affected by daily fluctuations in the host’s diet, activities, and health [10-12]. In addition, multiple fixed host factors, such as different states of immune competence, genetics, or gastrointestinal anatomy, likely also affect the composition, stability, or resilience of colonic microbiota [13-17]. Therefore, it is unclear whether divergence in post-FMT microbiota from that of donor implant material represents continued recovery, or whether these temporal changes are a general characteristic of host-associated gut microbiota in a changing host environment.
Here we describe both short- and long-term dynamic changes of fecal bacterial composition in four patients following FMT. All patients received microbiota from the same pre-qualified donor according to the standardized FMT protocol we described previously . Three patients received freshly prepared microbiota and one patient received microbiota that had previously been frozen. We compared pre- and post-FMT fecal microbial communities from these patients, as well as pre-FMT communities from 10 additional patients with multiply recurrent CDI (R-CDI), to the sequences of normal subjects described in the Human Microbiome Project . In addition, we compared temporal changes in fecal bacterial composition in recipients following FMT to temporal changes observed within samples from the donor.
Bacterial composition of fecal samples from patients with recurrent CDI becomes healthy and donor-like following FMT
Four patients (CD1 to CD4) with recurrent CDI were treated with FMT using material obtained from a single donor but from different time points, and fecal samples were collected from these patients before and after the procedure as well as from the donor at times of donation. Bacterial communities from these fecal samples were characterized by sequencing the V4 region of the 16S rRNA gene. Following trimming and quality filtering from a total of 12,536,492 sequences, we randomly subsampled to 5,000 sequences/sample in order to normalize read depth across all samples. All further analyses were performed using this rarefied read depth.
Fecal microbial communities remain dynamic following FMT
To determine whether this dynamic range of post-FMT microbial composition fits within the range seen across healthy individuals, we also compared communities in our samples to those in the HMP via weighted UniFrac and PCoA (Figure 2B). Again, fecal microbial communities prior to FMT were highly distinct from healthy fecal microbial communities, and following the procedure, these communities more closely resembled those of healthy individuals. Similar to the comparison with donor communities above, fecal microbial communities of recurrent CDI patients following FMT shifted within the cluster of communities from healthy individuals.
Rapid and substantial changes to Enterobacteriales in feces following FMT
We focused on these changes by examining the relative abundance of Enterobacteriales alone in each patient before and after FMT. The relative abundance of this taxon ranged from 44% to 82% in all four patient samples prior to FMT and rapidly dropped to undetectable levels within 1 week after the procedure. Moreover, abundance of this taxon remained low at 26 days after FMT, the latest time point shared by all four patients (Figure 4A), although other members of the Proteobacteria remain detectable if decreased in relative abundance (Figure 3). In addition, we generated individual value control charts based on the average abundance of this taxon in recurrent CDI patients. Compared to relative abundance, these control charts displayed the expected variation of the abundance of Enterobacteriales in these fecal samples. In all patients, the abundance of Enterobacteriales was above the expected variation (that is, more than three standard deviations above the mean relative abundance [the standard upper control limit, or UCL] of this order across all samples) prior to FMT, and rapidly fell below the upper control limit within 1 to 2 days after the procedure (Figure 4B). These results suggest that the relative abundance of Enterobacteriales significantly decreased in all patients soon after FMT to levels similar to donor samples and remained within a statistically expected range for the duration of sample collection (up to 151 days post-FMT).
Post-FMT communities are initially similar to donor samples but can later diverge
To further investigate how fecal microbial communities in these patients correlate to donor communities, we examined Pearson and Spearman correlations between donor and patient samples, which were common to each patient (pre-FMT samples and those up to 26 days post-FMT; Figure 5B,C and Additional file 4: Figure S2). While fecal microbial communities from patients before FMT were highly distinct from those in the donor, fecal microbial communities from samples 1 day after the procedure were highly correlated to donor communities via both Pearson and Spearman analyses in all patients. After the initial time point after FMT, the Pearson correlation values of patient to donor samples were highly variable within and across patients, although Spearman correlations remained high for three patients. To examine whether this variation is similar in healthy individuals, we determined Pearson and Spearman correlations within the four donor samples used in FMT, as well as eight additional donor samples from the same individual as a control. Results of this analysis revealed that donor microbiota also changed over time (Figure 5D). These findings suggested that the level of variability seen across patient post-FMT fecal microbial communities was within the range of normal microbiota behavior in a healthy individual.
Normalization and dynamic range of post-FMT patient fecal microbial communities are similar to donor communities
P values of normalization and dynamic range of patient samples sets versus donor set
It is now well understood that the fecal microbiota change substantially following FMT, typically shifting to fecal microbial communities more similar to those of the donor after transplant [2-8]. Here we show that these communities shift away from a dysbiotic state towards a composition that is representative of fecal microbial communities from hundreds of healthy individuals, collected in the HMP . Similarly to previous studies [4-8], the dysbiotic state in these patients with multiply recurrent CDI is characterized by a large expansion of Proteobacteria (primarily members of the order Enterobacteriales, which contains the family Enterobacteriaceae), and FMT is associated with reemergence of dominance by members of the Bacteroidetes and Firmicutes phyla.
Analysis of multiple donor and post-FMT samples demonstrates the dynamic behavior of fecal microbial communities over time. Both donor and recipient samples are characterized by highly dynamic shifts that nonetheless remain within the compositional range of normal fecal microbiota. This observation is consistent with known rapid responsiveness of the fecal microbiome to environmental inputs, such as dietary variations , and drifts in microbiota composition over time in healthy individuals .
The dynamic nature of intestinal microbiota is an intrinsic property, which should be taken into account when considering how therapeutic interventions, including FMT, impact its composition over time. In long-term post-FMT follow-up, Song and colleagues also noted dynamic changes in the fecal microbiome of R-CDI patients up to 16 weeks post-FMT . These investigators concluded that the fecal microbiome of post-FMT patients did not fully recover over this time, despite clinical recovery. Indeed, we observed divergence of microbiome in some of the patients away from the original implanted material over time. However, analysis of multiple donor samples showed that this movement is within the same dynamic range observed in the donor’s fecal microbiome. We therefore conclude that the dynamic behavior of microbiota needs to be taken into account in making comparisons between individuals, and should become an integral part of analysis of the success of FMT.
Three of the recipients in this study received freshly prepared microbiota, while one received frozen/thawed preparation. Use of frozen microbiota preparations is increasing in clinical practice , and its equivalency has not been rigorously established in randomized clinical trials. The ability to store microbiota allows the most up-to-date testing of the donor and fecal material for infectious pathogens, as some of the current tests may take several weeks to complete. Therefore, ability to preserve donor microbiota long-term is critical for its development as a therapeutic agent in clinical practice. Our results here, although limited in the number of patients, demonstrate indistinguishable behavior of fresh and frozen/thawed microbiota preparation.
The patients in this study did not have any significant gastrointestinal comorbidities. However, a significant proportion of patients with recurrent CDI have underlying inflammatory bowel disease, take potent immunosuppressive medications, or have multiple other medical problems [18,24]. The importance of these host factors in contributing to microbiota behavior is currently unknown, but is a subject of great interest . Understanding these influences will require analysis of multiple samples. Recently, Fuentes and colleagues  reported that some specific microbial groups and interactive networks are likely to be very important for the maintenance of microbiota in healthy individuals. However, although there is a great deal of effort focused on discovery of compositional differences in microbiota between normal subjects and individuals with different gastrointestinal and medical conditions, the dynamic behavior of fecal microbiota constitutes another dimension that may distinguish these cases. Thus, predictors of stable or dysbiotic intestinal microflora may also change over time. Further detailed studies of dynamic behavior of post-FMT microbiota may improve our understanding of causal connections between microbial communities and different disease states.
The fecal microbiota of patients with R-CDI continues to undergo change after FMT is performed, though these changes appear to fall within the range of normal variation of healthy individuals over time. Dynamic behavior is an intrinsic property of normal fecal microbiota and should be accounted for in comparing microbial communities among normal individuals and those with disease states.
Patients and donors
All patients suffered from multiply recurrent CDI refractory to standard antibiotic therapies. A single standard donor was used in the preparation of all fecal microbiota material as described previously . The Institutional Review Board at the University of Minnesota approved prospective collection of fecal specimens and their analysis. All patients satisfied the inclusion criteria for the FMT within our program, which included at least two spontaneous recurrences of CDI within a month of discontinuation of antibiotics and failure of at least one advanced antibiotic regimen such as a vancomycin pulse/taper protocol or vancomycin treatment followed by administration of rifaximin or fidaxomicin for 2 to 3 weeks. The specific clinical characteristics of patients involved in this study are summarized in Additional file 5: Table S2.
Fecal microbiota transplantation
FMT was done using a standardized preparation of concentrated fresh or frozen fecal bacteria via colonoscopy as previously described . All patients were treated with oral vancomycin, 125 mg four times daily, until 2 days prior to the procedure . The day before the procedure, patients received a polyethylene glycol-based colonoscopy prep (GoLYTELY® or MoviPrep®) to remove residual antibiotics and fecal material. Donor fecal microbiota was placed into the terminal ileum and/or cecum via the biopsy channel of the colonoscope. A total of 17 donor samples from the same individual were used in these studies. The CD1 to CD4 donor samples were given to patients CD1 to CD4, respectively. Patients CD1, CD3, and CD4 received freshly prepared fecal microbiota, while patient CD2 received a previously frozen preparation of fecal microbiota, all from the same standardized, anonymous donor.
Fecal samples were collected at home by the patients using swabs to sample feces deposited into a toilet hat immediately after production and stored frozen at approximately −20°C. Samples were subsequently transferred to the laboratory and stored at −80°C until used. Donor samples for DNA extraction were collected during processing of material for FMT and stored frozen at −80°C until used. Samples from patients CD1 to CD4 were obtained prior to FMT and between 1 to 151 days post-FMT, with daily collection until day 28, and weekly collection until day 84. Fecal material prior to FMT was obtained from patients CD5 to CD14.
DNA was extracted from donor and recipients’ pre- and post-FMT fecal samples using MOBIO PowerSoil DNA extraction kits (MOBIO, Carlsbad, CA, USA), according to the manufacturer’s instructions. Fecal DNA concentrations were measured using a QuBit DNA quantification system (Invitrogen, Carlsbad, CA, USA).
Extracted DNA was amplified using the EMP standard protocols at http://www.earthmicrobiome.org/ following the recommendations of Caporaso et al. . Briefly, F515/R806 primers were used, with 12-base Golay codes introduced on the 806 end to provide unique sample indices. Cycling and annealing conditions were as previously described .
DNA sequencing was performed as previously described  on an Illumina MiSeq platform using 2 × 150 bp paired-end reads and the Illumina v3 reagent chemistry.
Sequence processing and analysis
Sequence data was processed and analyzed using QIIME  according to the Illumina demultiplexing and processing protocol  and current quality-filtering recommendations , using the 1.8.0 pipeline and the default parameters in split_libraries_fastq.py. After quality control and demultiplexing, we picked close references at 97% similarity against the 97% similarity Greengenes database  version 13_8. All further analyses were performed at a rarefied depth of 5,000 reads/sample. EMPeror  was used for data visualization of BIOM-format  OTU tables. OTU analyses were performed by clustering at the 97% level with UCLUST , and data were integrated with the HMP dataset according to the protocols used for similar previous meta-analyses [15,32]. Sequences were analyzed by using both weighted and unweighted UniFrac , followed by principal coordinate analysis . Data were visualized using Phinch. The Phinch program provides an easy-to-use, browser-based, platform to visualize contingency tables along with their sample metadata (Bik et al., manuscript in preparation, https://github.com/PitchInteractiveInc/Phinch).
Analysis of microbiome stability and centrality
For each set of post-transplant patient samples, we assessed the similarity of that set to the set of reference samples from the donor (2,000 reads/sample). To reduce noise and compare patient samples along only relevant dimensions in UniFrac distance space, we applied PCoA to the unweighted UniFrac distance matrix containing only the post-transplant and donor samples for that donor-patient pair, then recalculated the distances using only the first n principal coordinates axes required to explain at least 80% of the variation in the distance matrix. An 80% cutoff was chosen to balance bias and overfitting. Distances were recalculated using Euclidean distances between points in PCoA space in order to convert PCoA coordinates to a distance matrix. The empirical P values for the ‘normality’ were obtained by comparing the mean distance between patient and donor samples to the histogram of within-donor distances (generated using all samples from a given donor by enumerating the pairwise distance between those samples). The empirical P values for the ‘dynamic range’ (stability) were obtained by comparing the mean distance within patient samples to the histogram of within-donor distances. These analyses were also performed using alternative parameters including, weighted UniFrac, Jensen-Shannon, root Jensen-Shannon, and Bray-Curtis.
Approval for this study was given by the University of Minnesota Institutional Review Board (Protocol Number: 0901M56962). All human subjects provided informed consent for participation in the study and collection and analysis of data. All human subjects gave their permission for their information to be published.
Clostridium difficile infection
fecal microbiota transplantation
Human Microbiome Project
lower control limit
operational taxonomic unit
principal coordinates analysis
recurrent Clostridium difficile infection
upper control limit
This research was supported, in part, by NIH Grant R21AI091907 (to AK and MJS). ARW was supported by a Doctoral Dissertation Fellowship from the University of Minnesota Graduate School and by the Dennis W. Watson Fellowship in Microbiology. Work done in the Knight lab was funded by the NIH, Crohn’s & Colitis Foundation of America, and the HHMI.
- Khoruts A, Weingarden AR. Emergence of fecal microbiota transplantation as an approach to repair disrupted microbial gut ecology. Immunology letters. 2014;162:77–81.View ArticlePubMedGoogle Scholar
- Khoruts A, Dicksved J, Jansson JK, Sadowsky MJ. Changes in the composition of the human fecal microbiome after bacteriotherapy for recurrent Clostridium difficile-associated diarrhea. J Clin Gastroenterol. 2010;44:354–60.PubMedGoogle Scholar
- Shahinas D, Silverman M, Sittler T, Chiu C, Kim P, Allen-Vercoe E, et al. Toward an understanding of changes in diversity associated with fecal microbiome transplantation based on 16S rRNA gene deep sequencing. mBio. 2012;3:e00338–12.View ArticlePubMed CentralPubMedGoogle Scholar
- Hamilton MJ, Weingarden AR, Unno T, Khoruts A, Sadowsky MJ. High-throughput DNA sequence analysis reveals stable engraftment of gut microbiota following transplantation of previously frozen fecal bacteria. Gut Microbes. 2013;4:125–35.View ArticlePubMed CentralPubMedGoogle Scholar
- van Nood E, Vrieze A, Nieuwdorp M, Fuentes S, Zoetendal EG, de Vos WM, et al. Duodenal infusion of donor feces for recurrent Clostridium difficile. N Engl J Med. 2013;368:407–15.View ArticlePubMedGoogle Scholar
- Weingarden AR, Chen C, Bobr A, Yao D, Lu Y, Nelson VM, et al. Microbiota transplantation restores normal fecal bile acid composition in recurrent Clostridium difficile infection. Am J Physiol Gastrointest Liver Physiol. 2014;306:G310–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Seekatz AM, Aas J, Gessert CE, Rubin TA, Saman DM, Bakken JS, et al. Recovery of the gut microbiome following fecal microbiota transplantation. mBio. 2014;5:e00893–00814.View ArticlePubMed CentralPubMedGoogle Scholar
- Fuentes S, van Nood E, Tims S, Heikamp-de Jong I, ter Braak CJ, Keller JJ, et al. Reset of a critically disturbed microbial ecosystem: faecal transplant in recurrent Clostridium difficile infection. ISME J. 2014;8:1621–33.View ArticlePubMedGoogle Scholar
- Song Y, Garg S, Girotra M, Maddox C, von Rosenvinge EC, Dutta A, et al. Microbiota dynamics in patients treated with fecal microbiota transplantation for recurrent Clostridium difficile infection. PLoS One. 2013;8:e81330.View ArticlePubMed CentralPubMedGoogle Scholar
- Dethlefsen L, Relman DA. Incomplete recovery and individualized responses of the human distal gut microbiota to repeated antibiotic perturbation. Proc Natl Acad Sci U S A. 2011;108 Suppl 1:4554–61.View ArticlePubMed CentralPubMedGoogle Scholar
- David LA, Maurice CF, Carmody RN, Gootenberg DB, Button JE, Wolfe BE, et al. Diet rapidly and reproducibly alters the human gut microbiome. Nature. 2014;505:559–63.View ArticlePubMed CentralPubMedGoogle Scholar
- Knights D, Ward TL, McKinlay CE, Miller H, Gonzalez A, McDonald D, et al. Rethinking "enterotypes". Cell Host Microbe. 2014;16:433–7.View ArticlePubMedGoogle Scholar
- Backhed F, Ley RE, Sonnenburg JL, Peterson DA, Gordon JI. Host-bacterial mutualism in the human intestine. Science. 2005;307:1915–20.View ArticlePubMedGoogle Scholar
- Lupp C, Robertson ML, Wickham ME, Sekirov I, Champion OL, Gaynor EC, et al. Host-mediated inflammation disrupts the intestinal microbiota and promotes the overgrowth of Enterobacteriaceae. Cell Host Microbe. 2007;2:119–29.View ArticlePubMedGoogle Scholar
- Koren O, Goodrich JK, Cullender TC, Spor A, Laitinen K, Backhed HK, et al. Host remodeling of the gut microbiome and metabolic changes during pregnancy. Cell. 2012;150:470–80.View ArticlePubMed CentralPubMedGoogle Scholar
- Lozupone CA, Li M, Campbell TB, Flores SC, Linderman D, Gebert MJ, et al. Alterations in the gut microbiota associated with HIV-1 infection. Cell Host Microbe. 2013;14:329–39.View ArticlePubMedGoogle Scholar
- Kawamoto S, Maruya M, Kato LM, Suda W, Atarashi K, Doi Y, et al. Foxp3(+) T cells regulate immunoglobulin a selection and facilitate diversification of bacterial species responsible for immune homeostasis. Immunity. 2014;41:152–65.View ArticlePubMedGoogle Scholar
- Hamilton MJ, Weingarden AR, Sadowsky MJ, Khoruts A. Standardized frozen preparation for transplantation of fecal microbiota for recurrent Clostridium difficile infection. Am J Gastroenterol. 2012;107:761–7.View ArticlePubMedGoogle Scholar
- Turnbaugh PJ, Ley RE, Hamady M, Fraser-Liggett CM, Knight R, Gordon JI. The human microbiome project. Nature. 2007;449:804–10.View ArticlePubMed CentralPubMedGoogle Scholar
- Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Huntley J, Fierer N, et al. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 2012;6:1621–4.View ArticlePubMed CentralPubMedGoogle Scholar
- Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.View ArticlePubMed CentralPubMedGoogle Scholar
- Bokulich NA, Subramanian S, Faith JJ, Gevers D, Gordon JI, Knight R, et al. Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing. Nat Methods. 2013;10:57–9.View ArticlePubMed CentralPubMedGoogle Scholar
- McDonald D, Price MN, Goodrich J, Nawrocki EP, DeSantis TZ, Probst A, et al. An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J. 2012;6:610–8.View ArticlePubMed CentralPubMedGoogle Scholar
- Vazquez-Baeza Y, Pirrung M, Gonzalez A, Knight R. EMPeror: a tool for visualizing high-throughput microbial community data. GigaSci. 2013;2:16.View ArticleGoogle Scholar
- McDonald D, Clemente JC, Kuczynski J, Rideout JR, Stombaugh J, Wendel D, et al. The Biological Observation Matrix (BIOM) format or: how I learned to stop worrying and love the ome-ome. GigaSci. 2012;1:7.View ArticleGoogle Scholar
- Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26:2460–1.View ArticlePubMedGoogle Scholar
- Lozupone CA, Stombaugh J, Gonzalez A, Ackermann G, Wendel D, Vazquez-Baeza Y, et al. Meta-analyses of studies of the human microbiota. Genome Res. 2013;23:1704–14.View ArticlePubMed CentralPubMedGoogle Scholar
- Lozupone C, Knight R. UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol. 2005;71:8228–35.View ArticlePubMed CentralPubMedGoogle Scholar
- Faith JJ, Guruge JL, Charbonneau M, Subramanian S, Seedorf H, Goodman AL, et al. The long-term stability of the human gut microbiota. Science. 2013;341:1237439.View ArticlePubMed CentralPubMedGoogle Scholar
- Youngster I, Sauk J, Pindar C, Wilson RG, Kaplan JL, Smith MB, et al. Fecal microbiota transplant for relapsing Clostridium difficile infection using a frozen inoculum from unrelated donors: a randomized, open-label, controlled pilot study. Clin Infect Dis. 2014;58:1515–22.View ArticlePubMedGoogle Scholar
- Kelly CR, Ihunnah C, Fischer M, Khoruts A, Surawicz C, Afzali A, et al. Fecal microbiota transplant for treatment of Clostridium difficile infection in immunocompromised patients. Am J Gastroenterol. 2014;109:1065–71.View ArticlePubMedGoogle Scholar
- Gevers D, Kugathasan S, Denson LA, Vazquez-Baeza Y, Van Treuren W, Ren B, et al. The treatment-naive microbiome in new-onset Crohn's disease. Cell Host Microbe. 2014;15:382–92.View ArticlePubMed CentralPubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.