<title>Introduction</title> <p>The human gut microbiota has been shaped by the long co-evolutionary history of symbiotic host-microbe interaction ##REF##16497592##[1]## and has evolved to play an important role in maintaining human health by preventing colonization by pathogens, degrading dietary and <italic>in situ</italic>-produced compounds, producing nutrients, and shaping and maintaining the normal mucosal immunity ##UREF##0##[2]##–##REF##9849357##[4]##. Recently, other important functions of commensal microbiota became apparent including an anti-inflammatory effect on gut epithelium ##REF##10968793##[5]##, ##REF##14691478##[6]##, the influence on lipid metabolism of the host ##REF##15505215##[7]## and association with obesity ##REF##17183309##[8]##, ##REF##17183312##[9]## as well as the involvement in intestinal homeostasis, repair and angiogenesis ##REF##12432102##[10]##, ##REF##15260992##[11]##.</p> <p>Recently, the intricate molecular mechanisms behind this host-microbe cross-talk began to emerge. The innate immune system of the mammalian host senses bacteria through an impressive array of receptors, called Toll-like (TLR) and Nod-like (NLR) (also called CATERPILLER, NOD-LRR, and NACHT-LRR) receptors. The first mammalian TLR was discovered a decade ago ##REF##9237759##[12]## and the functionality of the first NLRs (then called NODs) as sensors of bacterial ligands was established in 2001 ##REF##11058605##[13]##. TLRs are predominantly but not exclusively located in the transmembrane region of host cells, while NLR proteins are mostly located in the cytoplasm ##REF##17680011##[14]##. Both groups of receptors contain leucine-rich repeats (LRR), which are involved in detection of danger signals and, in particular, of whole bacteria and bacterial products through the recognition of pathogen associated molecular patterns (PAMP). PAMP signalling through the innate immunity receptors appeared to be important not only for the recognition of, and launching defences against, the invading pathogens but also for the maintenance of intestinal epithelial homeostasis and for protection against gut injury and associated mortality ##REF##15260992##[11]##. Bacterial products are also crucial for the normal immune system development; for example, a single bacterial polysaccharide (PSA) from a commensal bacterium <italic>Bacteroides fragilis</italic> was shown is capable of directing the complete cellular and physical maturation of the developing immune system ##REF##16009137##[15]##. In murine experimental colitis models, signalling by probiotic bacterial DNA ligands through TLR9 exerted a strong anti-inflammatory effect ##REF##14762789##[16]## while signalling through TLR2 seems important for regulation of mucosal inflammation and maintenance of intestinal epithelial barrier integrity ##REF##15236188##[17]##, ##REF##17408640##[18]##.</p> <p>Another type of regulation of the host metabolism by microbes can be through the low molecular weight bacterial metabolites such as short chain fatty acids (SCFAs) and, in particular, butyrate. Butyrate is the preferred energy source for the colonic epithelium and may play an important role in colonic health and prevention of colorectal cancer and colitis ##REF##7496324##[19]##–##REF##10716678##[21]##.</p> <p>Because of limitations imposed by the extreme diversity of microbiota in the conventional gut, which does not allow monitoring the behaviour of a single bacterium, much less is known of how the bacterial counterpart responds to the host-generated molecular signals. Initially, the bacterial response was considered to be limited to various strategies used by gut bacteria to forage the endogenously produced substrates of the host ##REF##10449780##[22]##, ##REF##15790854##[23]##. Interestingly, however, one of the aspects of this interaction, namely the specific fucosylation, seems important for both sides involved in the host-bacterial dialogue. On one side, the intestinal bacteria induce mucin fucosylation in the gut epithelium thus creating a specific nutritional econiche ##REF##10449780##[22]##, ##REF##17673542##[24]##. On the other side, fucosylation of the cell surface in a mammalian-like fashion appeared to be essential for the gut bacteria to establish the successful mutualistic relationship with the host. The lack of this property may result in the failure of bacteria to colonize the mammalian intestine ##REF##15774760##[25]##. Yet another “language” of host-bacterial communication could be through the quorum sensing (QS), which utilizes hormone-like compounds and, in <italic>E. coli</italic>, has been implicated in a potential cross-communication between the luxS/AI-3 bacterial QS system and the epinephrine host signalling system ##REF##12847292##[26]##. Bacteria may also sense the immune status of the host and respond by the expression of a quorum-sensing dependent virulence determinant ##REF##16051797##[27]##. These two examples of quorum-sensing-mediated signalling between the host and bacteria are taken from the pathogenic bacteria but it cannot be excluded that this type of communication is common among gut symbiotic bacteria as well.</p> <p>What happens if this fine-tuned molecular cross-talk is compromised and the host fails to recognize bacterial components as, for example, in the case of inflammatory bowel disease (IBD) such as Crohn's disease (CD)? In a subset of CD patients, mutations in the NOD2/CARD15 protein have been associated with susceptibility to the disease ##REF##11385576##[28]##, ##REF##11385577##[29]##. Bacterial ligand for this receptor has been recently identified as bacterial muramyl dipeptide ##REF##12514169##[30]##, ##REF##12527755##[31]## and failure to recognise this ligand in CD patients leads to the breach of tolerance and to the launch of an aggressive Th1-skewed response against the antigens of normal diet and commensal bacteria as well as to production of autoantibodies ##REF##12198693##[32]##–##UREF##2##[34]##. The microbiota of these patients demonstrates diminished diversity, in particular, due to the depletion of bacteria belonging to the <italic>Firmicutes</italic> and <italic>Bacteroidetes</italic> phyla ##REF##16188921##[35]##, ##REF##17699621##[36]##. Interestingly, the loss of the transcription factor T-bet, which regulates the innate immune system, influences bacterial populations in the murine gut in a way that they become colitogenic and this colitis is communicable to genetically intact hosts ##REF##17923086##[37]##. Presently, it is not clear how the composition of the normal microbiota is restructured to become a “pathobiota” but, nevertheless, this is the first experimental proof of the view that the community as a whole may demonstrate pathogenic properties, not a single pathogenic bacterium.</p> <p>In our study, we investigated the effect of mutations in a single gene of the host on commensal gut microbiota. The gene under investigation, <italic>MEFV</italic> (for MEditerranean FeVer), encodes a protein called pyrin/marenostrin, which is involved in regulation of innate immunity ##UREF##3##[38]##. Structurally, this 781-residue protein consists of a 92-amino acid N-terminal PYRIN (PYD) domain, a B-box zinc finger, a coiled-coil region and a ∼200-amino acid C-terminal B30.2/rfp/SPRY domain ##REF##16724804##[39]##. Mutations in this gene lead to an autoinflammatory disorder, familial Mediterranean fever (FMF, MIM249100), which is characterised by short recurrent self-resolving attacks of fever and polyserositis, with no clinical signs of disease in remission periods. The carrier rate and occurrence of disease is high in populations originating from the Mediterranean basin, including Sephardic Jews, Armenians, Arabs, and Turks. Articles describing positional cloning of the gene by two independent consortia appeared in 1997 ##REF##9288758##[40]##, ##REF##9288094##[41]## and to date 166 <italic>MEFV</italic> mutations and single nucleotide polymorphisms (SNPs) are detected in the gene (<ext-link ext-link-type="uri" xlink:href="http://fmf.igh.cnrs.fr/infevers">http://fmf.igh.cnrs.fr/infevers</ext-link>). The hotspot of <italic>MEFV</italic> mutations is localized in exon 10, which encodes the B30.2/rfp/SPRY domain, and this single locus inheritance greatly facilitated our analysis of host genetics in comparison with other disorders with more complex genetics. Genome-wide scans for IBD susceptibility, for example, have resulted in identification of seven loci that are confirmed and replicated in several studies, but with many more that need independent confirmation and verification ##REF##17894847##[42]##.</p> <p>The present work is the first attempt to analyze simultaneously the genetics of both sides involved in the host-gut microbiota relationship. All study participants were analyzed for <italic>MEFV</italic> mutations and the composition of gut microbiomes was assessed by two culture-independent techniques, fluorescent in situ hybridization (FISH) and 16S rDNA gene library analyses. For the first time, our findings conclusively establish the key role of the host genotype in host-microbe interaction, which determines a highly specific composition of commensal microbiota in the human gut.</p>
<title>Materials and Methods</title> <title>Subjects and sampling</title> <p>A total of 19 Armenian FMF patients with clinically confirmed disease, seven healthy individuals of the same ethnic background, and one healthy subject of non-Armenian origin participated in this study (##TAB##0##Table 1##). All FMF cases were diagnosed based on Tel-Hashomer criteria ##REF##9336425##[74]##. Blood and fecal samples were collected from FMF patients (15 males, four females; aged from 15 to 64 years, mean age–27 years) and control subjects (one male, seven females; aged from 32 to 67 years, mean age–43.8 years) for genotyping and gut bacterial diversity analyses (see below). None of the FMF patients and healthy individuals had used antibiotics within the three months prior to sampling. In the beginning of this study, ten FMF patients were not receiving any drugs and nine were undergoing regular colchicine therapy. During the study, one patient from the colchicine-free group was prescribed this therapy. All subjects were informed about the aim of this study and gave their consent to participate in it. The study was approved by the local ethical committee at IMB.</p> <title>Genetic diagnostics</title> <p>Venous blood samples were obtained from the clinically confirmed 19 FMF patients and eight healthy controls to reveal the <italic>MEFV</italic> mutation carrier status. Genetic analysis was based on screening of two mutational hot spots in the <italic>MEFV</italic> gene (exons 2 and 10). For this, genomic DNA was isolated from the anticoagulated venous blood samples using the Wizard Genomics DNA Purification kit (Promega, UK), according to the manufacturer's instructions. Exon 2 of the <italic>MEFV</italic> gene was PCR amplified from genomic DNA using a newly designed forward and reverse primer set: Exon2F - <named-content content-type="gene">5′-ATTCTCTCTCCTCTGCCCTG-3′</named-content>, Exon2R - <named-content content-type="gene">5′-CCATTCTTTCTCTGCAGCCG-3′</named-content>, yielding a 839 bp amplicon. Exon 10 was amplified using the primers Exon10F <named-content content-type="gene">5′-CCAGAAGAACTACCCTGTCCC-3′</named-content> and Exon10R <named-content content-type="gene">5′-TCCTCCTCTGAAATCCATGG-3′</named-content>, yielding a product of 887 bp ##REF##16724804##[39]##. For exon 10, amplification was performed in a final volume of 50 µl, in a mix containing 10 mM Tris-HCl, 1.5 mM MgCl<sub>2</sub>, 50 mM KCl, 200 µM of each dNTP, 10.0 pmol of each primer, approximately 20 ng of genomic DNA, and 1 U of Red Taq polymerase (Bioline, UK). The cycling conditions included the initial denaturation step for 5 min at 94°C, followed by 35 cycles consisting of denaturation at 94°C for 30 sec, annealing at 61°C for 60 sec, and extension at 72°C for 60 sec, with a final extension at 72°C for 10 min. For exon 2, PCR amplification was performed similarly, except that the PCR mix contained 1 M betaine (Sigma-Aldrich, UK). Amplicons of the expected size were excised from 1.5% agarose gel and purified using a Wizard SV Gel and PCR Clean-up system (Promega, UK), according to the manufacturer's instructions. Purified PCR products were sequenced on both directions, using the same primer set as for generation of templates and the sequences were read on an automated 8-channel capillary sequencer (Beckman, USA) The wild-type sequences of exons 2 and 10 of the <italic>MEFV</italic> gene were retrieved from GenBank (accession number AF111163).</p> <title>Generation of 16S rDNA clone libraries</title> <p>In order to derive a detailed phylogenetic inventory of the gut microflora in FMF and in healthy controls, we constructed bacterial 16S rDNA gene libraries from fecal samples. For this, freshly voided feces were collected into sterile containers from 12 FMF patients in remission, three patients in the acute stage of FMF, six healthy Armenians, and one healthy subject of non-Armenian origin (##TAB##0##Table 1##). One patient, FMF 13-123, provided the samples both in remission and attack. Total DNA was extracted from fecal samples using QIAamp DNA Stool Mini Kit (Qiagen, UK), according to the manufacturer's instructions. DNA was additionally purified by phenol-chloroform-isoamyl alcohol (25∶24∶1) extraction, ethanol precipitated, and dissolved in TE buffer. Amplification of the 16S rRNA genes was carried out with the universal primer set–27F (<named-content content-type="gene">5′-AGAGTTTGATCMTGGCTCAG-3′</named-content>; positions 8 to 27 in the <italic>Escherichia coli</italic> 16S rRNA gene) and 1492R (<named-content content-type="gene">5′-ACGGCTACCTTGTTACGACTT-3′</named-content>; positions 1510 to 1492 in the <italic>E. coli</italic> 16S rRNA gene) ##REF##1987160##[75]##. The PCR reaction mixture contained 10 mM Tris-HCl, 2 mM MgCl<sub>2</sub>, 50 mM KCl, 200 µM of each dNTP, 10.0 pmol of each primer, 1U of Bio-X-Act Short DNA polymerase (Bioline, UK), and 100 ng of template DNA in a final volume of 50 µl. To reduce PCR bias caused by elevated number of PCR cycles ##REF##12054235##[57]##, amplifications were performed using the following conditions: initial denaturation of template DNA at 94°C for 5 min, followed by 10 cycles of denaturation at 94°C for 30 sec, annealing at 57°C for 30 sec, and extension at 72°C for 2 min, with a final extension at 72°C for 10 min. PCR products were visualized on 1% agarose gel in TBE buffer stained with GelStar (Cambrex, UK). 16S PCR products were purified using Wizard SV Gel and PCR Clean-up system (Promega, UK), according to the manufacturer's instructions. This was followed by ethanol precipitation and the final pellets were suspended in 5 µl of TE buffer (pH 8.0). The purified PCR products were directly cloned into pCR-4 cloning vector and then transformed into <italic>E. coli</italic> TOP10 chemically competent cells using the TOPO TA Cloning Kit (Invitrogen, UK), according to the manufacturer's protocols. Colonies were blue/white screened on Luria-Bertani (LB) agar plates supplemented with ampicillin (50 µg/ml) and X-Gal (80 µg/ml). No IPTG was necessary in this system because the strain and plasmid do not carry the <italic>lac</italic> repressor. White recombinant colonies were randomly picked up into 200 µl of LB medium with ampicillin (50 µg/ml) in a 96-well format, sealed with a gas-permeable membrane, and incubated overnight at 37°C in a horizontal shaker. Two insert-surrounding primers, T3 (<named-content content-type="gene">5′-ATTAACCCTCACTAAAGGGA-3′</named-content>) and T7 (<named-content content-type="gene">5′-TAATACGACTCACTATAGGG-3′</named-content>) were used to amplify inserts using a 35-cycle colony PCR. Products were visualized after electrophoresis on 1% agarose gel in TAE buffer and ethidium bromide staining.</p> <title>Sequencing and phylogenetic analysis</title> <p>Successfully amplified 16S rDNA inserts were purified and sequenced on the automated 8-channel capillary sequencers (Beckman, USA). The sequencing reactions were carried out with 926R primer (<named-content content-type="gene">5′-CCGTCAATTCCTTTGAGTTT-3′</named-content>; positions 926 to 907 in the <italic>E. coli</italic> 16S rRNA gene) ##REF##2413450##[76]##. Nucleotide sequences were aligned using the multiple sequence alignment program CLUSTALX v1.83 ##UREF##7##[77]##. Each sequence was manually edited in conjunction with its chromatogram. Sequences were examined for the possible chimeras using CHIMERA_CHECK online analysis at RDP-II website (<ext-link ext-link-type="uri" xlink:href="http://rdp.cme.msu.edu/cgis/chimera.cgisuSSU">http://rdp.cme.msu.edu/cgis/chimera.cgisuSSU</ext-link>). Similarity search of the 16S rDNA sequences against database entries was performed using online Basic Local Alignment Search Tool (BLAST) program at the NCBI website (<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/BLAST">http://www.ncbi.nlm.nih.gov/BLAST</ext-link>). Distance and similarity matrices were constructed with the DNADIST program v3.6 from PHYLIP (<ext-link ext-link-type="uri" xlink:href="http://evolution.genetics.washington.edu/phylip/phylip.html">http://evolution.genetics.washington.edu/phylip/phylip.html</ext-link>) according to the Jukes-Cantor model. Phylogenetic analysis was performed using the neighbour-joining method ##REF##3447015##[78]##. Statistical significance of branching was verified by bootstrap analysis involving the construction and analysis of 1000 trees from bootstrapped data sets. Operational taxonomic units (OTUs), or phylotypes, were defined at 99% level as recommended for gut microbiota analyses ##REF##15831718##[59]##. To assign sequences quickly and accurately to OTUs, additionally, Jukes-Cantor corrected distance matrices were analyzed with DOTUR (Distance-Based OTU and Richness) by using the furthest-neighbor method ##REF##15746353##[79]##.</p> <p>About the half of OTUs (165 out of 346) were sequenced to completion. For this, a set of additional three universal bacterial primers, 27F (<named-content content-type="gene">5′-AGAGTTTGATCMTGGCTCAG-3′</named-content>; positions 8 to 27 in the <italic>E. coli</italic> 16S rRNA gene), 519R (<named-content content-type="gene">5′-GTATTACCGCGGCTGCTG-3′</named-content>; positions 536 to 519 in the <italic>E. coli</italic> 16S rRNA gene) ##REF##2413450##[76]## and RP2 (<named-content content-type="gene">5′- ACGGCTACCTTGTTACGACTT -3′</named-content>; positions 1510 to 1492 in the <italic>E. coli</italic> 16S rRNA gene) was used. These partial rRNA gene sequences (the average single sequence read length was ∼700 bp), including the sequence previously generated by the 926R primer, were assembled with ChromasPro v1.33 program to produce a consensus sequence for each of the 165 clones. Initially, all sequence data were assembled according to their overlapping sequence fragments into a contiguous consensus sequence using default program settings or, in some cases, adjusting the program settings, if necessary, followed by verification by eye. Afterwards, all chromatogram contigs obtained were manually edited removing ambiguous chromatogram areas, which resulted in the nearly full-length 16S bacterial rDNA sequences of about 1500 bp.</p> <p>All full-length sequences were checked for possible chimeras using Chimera Check program, v2.7 (online analysis at Ribosomal Database Project II website, <ext-link ext-link-type="uri" xlink:href="http://rdp.cme.msu.edu/cgis/chimera.cgisuSSU">http://rdp.cme.msu.edu/cgis/chimera.cgisuSSU</ext-link>). Chimeras and sequences of poor quality were omitted from further phylogenetic analyses. Phylogenetic analysis was performed in two stages. First, the almost full-length sequences from our libraries were combined with the closest relatives, cultivated and uncultivated, from databases, aligned, and subjected to phylogenetic analysis using the neighbor-joining method ##REF##3447015##[78]##. In the second stage, the partial sequences were aligned to the full-length sequence profile using the guide tree generated during the first stage. A total of 1328 validated 16S rDNA sequences participated in the final tree computing. Validation of the tree branching was done with bootstrap analysis ##REF##28561359##[80]## involving 1,000 reiterations with re-sampled data sets.</p> <title>Community structure analysis</title> <p>Rarefaction and collector's curves of observed phylotypes, richness estimates and diversity indices were determined with DOTUR program using Jukes-Cantor corrected distance matrix. The bias-corrected Chao 1 richness estimator was calculated after 1000 randomizations of sampling without replacement. Collector's curves of observed and estimated (Chao 1 and the abundance-based coverage estimator, ACE) richness were constructed. Diversity was estimated by the Shannon and Simpson indices, the Simpson reciprocal index was calculated as 1/D, and another version of the Simpson diversity index–as 1-D. The evenness was calculated as E = e<sup>H</sup>/N, where H is the Shannon diversity index, N is the number of phylotypes. The Good's coverage percentage was calculated with the formula [1−(n/N)]×100, where n is the number of phylotypes in a sample represented by one clone (singletons) and N is the total number of sequences in that sample ##UREF##8##[81]##.</p> <title>Fluorescent in situ hybridization (FISH)</title> <p>FISH analysis was applied to quantify the predominant groups of bacteria in fecal samples using the universal bacterial and group-specific 16S rRNA-based Cy3-labelled oligonucleotide probes (##TAB##8##Table 9##). The probes hybridized with the next bacterial groups: Eub338 with the total bacteria ##REF##2200342##[82]##; Bac303 with the <italic>Bacteroides</italic> group ##REF##8704951##[83]##; Erec482 with the members of the <italic>Clostridium coccoides</italic> group (cluster XIVa) ##REF##9726880##[84]##; Rrec584 with the <italic>Eubacterium rectale-Roseburia</italic> cluster, which is a component of clostridial cluster XIVa ##REF##16957265##[85]##; Fprau645 for <italic>Faecalibacterium prausnitzii</italic>-related bacteria–a component of clostridial cluster IV ##REF##11403393##[86]##; the Rfla729 and the Rbro730 probes are specific for ruminococci in clostridial cluster IV ##REF##12039758##[87]##; Prop853 for the members of clostridial cluster IX ##REF##16000778##[88]##; Bif164 for the genus <italic>Bifidobacterium</italic>
##REF##7487040##[89]##; Ato291 for the <italic>Atopobium</italic> group ##REF##11010909##[90]##; EnterobactD for the genus <italic>Enterobacteriaceae</italic>
##REF##12067375##[91]##; and Lab158 for lactobacilli and enterococci ##UREF##9##[92]## (##TAB##8##Table 9##). Together, the group-specific probes detected the vast majority of the total bacterial cells in the normal human gut ##REF##9726880##[84]##. The Rrec584 probe was excluded from this summation, since it detects the <italic>E. rectale-Roseburia</italic> cluster within the group already covered by the probe Erec482 ##REF##16957265##[85]##. The probe Rfla729 was used in conjunction with the probe Rbro730 as the target bacterial groups are overlapping ##REF##12039758##[87]##.</p> <p>Freshly voided feces from six FMF patients in remission, seven in the acute phase, six healthy Armenians and one healthy subject of non-Armenian descent were collected (##TAB##0##Table 1##). Two patients, FMF 13-123 and FMF 87, provided fecal samples at both disease stages. The samples in sterile plastic bags were kept at 4°C for no longer than 12 h prior processing. The specimens were kneaded mechanically for 5 min at 4°C to distribute the sample evenly. 0.5 g of each sample was suspended in 4.5 ml of filtered (0.2 µm pore-size filter) ice-cold phosphate-buffer-saline (PBS) and vortexed with a dozen of glass beads for at least 3 min to homogenize the sample and dislodge the bacteria from feed particles. The suspension was centrifuged at 700×g for 1 min to remove debris. One ml of supernatant was added to three ml of 4% paraformaldehyde (PFA) in PBS, fixed at 4°C for 16 hours, washed with ice-cold PBS, and stored in 50% (vol/vol) ethanol-PBS at −20°C until analysed.</p> <p>Depending on the expected number of target cells, the samples were diluted 40- to 1600-fold, and 10 µl of a diluted bacterial suspension was applied to gelatin-coated slides, air-dried and fixed in 100% ethanol. 10 µl of 50 ng/µl solution of oligonucleotide probes in 100 µl of hybridization buffer was added to the slides and the samples were hybridized overnight (except for the Bac303 probe, which was incubated for 2 h) in a precision incubator at different temperatures, depending on the probe (##TAB##8##Table 9##). If more stringent conditions were required, formamide was added to the hybridization buffer. To improve bacterial cell permeability samples were incubated before hybridization in 10 µl of lysozyme (final concentration 1 mg/ml) in 100 mM Tris-HCl (pH 8.5) at 37°C. The full panel of FISH oligonucleotide probes with their target bacterial groups and hybridization conditions used in this study is shown in ##TAB##8##Table 9##.</p> <p>After hybridization, the slides were soaked in a washing buffer for 20 min at 50°C, rinsed with Milli-Q water, and rapidly dried with compressed air. Slides then were covered with 50 µl of Vectashield (Vector Laboratories, Burlington, California) to prevent the fading of fluorescence and a coverslip. The enumeration was performed when the hybridization signals were clear and morphologically distinguishable. Bacterial cells were counted automatically using image analysis software Quantimet HR600 with a Leica DMRXA epifluorescence microscope (Wetzlar, Germany). Depending on the number of labelled cells, from 20 to 50 fields were counted to bring the coefficient of variance value below 10.</p> <title>Statistical analyses</title> <p>Statistical analyses were performed using the Statsoft Statistica package (<ext-link ext-link-type="uri" xlink:href="http://www.statsoft.com">www.statsoft.com</ext-link>). Comparisons were made among groups of subjects using the Mann-Whitney <italic>U</italic> test and Student's t-test. Multiple comparisons were carried out using Kruskal-Wallis test. In order to classify cases data were subjected to discriminant function analysis (DA). The P-values below 0.05 were considered statistically significant.</p> <p>Statistical comparison of 16S rDNA clone libraries was performed using online library compare tool at RDP-II website (<ext-link ext-link-type="uri" xlink:href="http://rdp.cme.msu.edu/comparison/comp.jsp">http://rdp.cme.msu.edu/comparison/comp.jsp</ext-link>). This tool uses the RDP Naive Bayesian classifier to provide rapid classification of library sequences into the bacterial taxonomy proposed by the Bergey's Trust. The P values &lt;0.01 were considered as having significant differences.</p> <p>
<bold>Data deposition:</bold> The sequences reported in this paper have been deposited in GenBank with accession numbers EU531867-EU532010 and EU532765-EU533948.</p>
<title>Results</title> <title>MEFV Genotyping</title> <p>Nineteen FMF patients with clinically confirmed disease and eight healthy individuals were subjected to genotyping (##TAB##0##Table 1##). We performed genotyping of the exon regions of the <italic>MEFV</italic> gene, which is thought to be a single genetic determinant of FMF, in particular exons 2 and 10, where the majority of disease-associated mutations were identified (<ext-link ext-link-type="uri" xlink:href="http://fmf.igh.cnrs.fr/infevers">http://fmf.igh.cnrs.fr/infevers</ext-link>). FMF-associated mutations were found in all FMF patients displaying the clinical signs of the disease, as well as in one healthy Armenian control subject (in a simple heterozygous state) without disease-related symptoms (##TAB##0##Table 1##). All mutations were located in exon 10, only synonymous substitutions were detected in exon 2, irrespectively of the health status. The majority of mutations in FMF patients were M694V, which were present in the compounded form with other disease alleles in eight patients and in a simple heterozygous form in six patients. In the latter group, no other disease-related mutations in exon 2, except neutral SNPs, were detected. Homozygotes were represented exclusively by the M680I mutation (##TAB##0##Table 1##). Interestingly, the SNPs in the control group were all in the heterozygous state, while the homozygous forms were encountered only among FMF patients.</p> <title>Phylogenetic and comparative analysis of 16S libraries</title> <p>Three 16S rDNA libraries from the fecal samples of healthy controls as well as of FMF patients in remission and disease attack were generated (designated as S, SR and SA in ##TAB##0##Table 1##). After the quality control a total of 1328 validated sequences (572 for healthy controls, 629 for FMF remission and 127 for FMF attack) were analyzed and phylogenetic analysis was performed to establish taxonomic positioning of sequences obtained (##TAB##1##Table 2##). Among 1328 clones analyzed, there were 268 distinct OTUs and the vast majority (95.95%) fell into the two major phyla, <italic>Bacteroidetes</italic> with 99 OTUs (38.15%) and 714 sequences (53.76%) and <italic>Firmicutes</italic> with 158 OTUs (57.8%) and 580 sequences (43.67%). <italic>Proteobacteria</italic> were represented by 9 OTUs (3.47%) and 30 sequences (2.26%). Only four <italic>Actinobacteria</italic> sequences (2 OTUs) were detected, which were all derived from one FMF patient in remission.</p> <p>In all three libraries, bacteria belonging to the <italic>Bacteroidetes</italic> phylum were the most abundant, containing more than the half of the sequences in each library (from 50.72% to 56.64%) but a relatively smaller proportion of OTUs (from 35.19% to 40.91%), whereas the <italic>Firmicutes</italic> phylum had the larger proportion of phylotypes (from 55.30% to 60.49%) with the less copious sequence numbers (from 41.43% to 45.95%) (##TAB##1##Table 2##). At the phylum level, the relative proportions of <italic>Bacteroidetes</italic> and <italic>Firmicutes</italic> were not significantly different among the three libraries, however, pairwise comparisons of each 16S rDNA library to every other library by using the online library compare tool at RDP-II (Naive Bayesian rRNA Classifier) revealed significant differences within the lower taxa levels (##TAB##2##Table 3##). In particular, in asymptomatic FMF patients in remission, the proportion of <italic>Enterobacteriaceae</italic>, <italic>Acidaminococcaceae</italic>, <italic>Ruminococcus</italic> and <italic>Megasphaera</italic> was significantly increased in comparison with control subjects, while <italic>Roseburia</italic> was significantly reduced (##TAB##2##Table 3##). In the active disease, the proportion of <italic>Prevotellaceae, Dialister</italic> and <italic>Prevotella</italic> was significantly lower than in healthy controls while the <italic>Porphyromonadaceae</italic>, <italic>Phascolarctobacterium</italic>, <italic>Faecalibacterium</italic>, and <italic>Parabacteroides</italic> were significantly increased (##TAB##2##Table 3##). Interestingly, there were also differences due to the disease state (remission or attack) and the affected bacterial groups were <italic>Acidaminococcaceae</italic>, <italic>Porphyromonadaceae</italic>, <italic>Megasphaera</italic>, <italic>Dialister</italic>, <italic>Faecalibacterium</italic>, and <italic>Parabacteroides</italic> (##TAB##2##Table 3##).</p> <title>Community analysis of 16S rDNA libraries</title> <p>The estimates of diversity, richness, evenness and library coverage for the three 16S rDNA clone libraries studied are shown in ##TAB##3##Table 4##. At the phylotype level, the most diverse gut bacterial community was observed in FMF patients in remission. The slightly lower indices were found for the healthy microbiota, while the bacterial diversity during disease attacks was consistently lower. The similar trend was observed for another community parameter, species richness, which was the highest in the gut of FMF patients in remission and in healthy volunteers but much lower during disease attacks. On the contrary, evenness index, which describes how evenly the individual sequences are represented in phylotypes, was the highest in the disease attack community (##TAB##3##Table 4##). Good's coverage was 89.35–90.73% for the sequence sets from FMF patients in remission and healthy controls but only 77.2% for sequences from FMF patients in attack reflecting the lower number of subjects and sequences collected.</p> <p>Rarefaction analysis of clone libraries essentially confirmed these findings and suggested that FMF patients in remission possess the most diverse gut bacterial community, which is even higher than in healthy subjects (##FIG##0##Fig. 1##). Because of limited number of sequences in the clone library from the acute disease, we extrapolated the curve using an empiric regression function (for details, see the footnote to ##FIG##0##Fig. 1##). The extrapolated curve confirmed the initial findings suggesting diminished diversity of gut bacteria in FMF attack in comparison with disease remission and control.</p> <p>Collector's curves of the observed and estimated phylotype richness are shown in ##FIG##1##Fig. 2##. In healthy subjects, both estimates, Chao1 and abundance-based coverage estimator (ACE), were highly similar and, up to the point of collection of ca. 320 clones, were in a steady and congruent increase together with the observed phylotype number (##FIG##1##Fig. 2A##). In this sampling area, almost any new clone sequence added represented a novel phylotype. Adding new sequences after this point, however, led to the gap expansion between the estimated and observed phylotype numbers suggesting repeated sampling of the phylotypes already present in the sample. In FMF remission, the gap between the observed and estimated phylotype richness was constant after the sampling point of 290 clones (##FIG##1##Fig. 2B##). Confirmatory to the previous analyses, the phylotype number estimates at the sampling point with 500 clones gave the higher numbers of phylotypes in FMF remission than in the norm (220 vs. 160). The difference between the estimates and observed phylotype richness is highest in acute FMF possibly reflecting the insufficient coverage in this group (77.2%).</p> <title>FISH analysis</title> <p>Twelve Cy3-labelled oligonucleotide probes that target the total and specific predominant groups of human gut bacteria were used in this study to enumerate bacteria in the feces of FMF patients and healthy controls (for details, see <xref ref-type="sec" rid="s4">Materials and Methods</xref>). Hybridization signals were detected for all bacterial groups, indicating their presence as dominant microbiota in all fecal samples analyzed.</p> <p>FISH analysis revealed that the composition of fecal bacteria in FMF patients in both inactive and active phases of the disease is different from that of the healthy controls. The results of enumeration of bacterial cells in each group and their relative proportions versus the universal bacterial probe Eub338 are summarized in ##TAB##4##Table 5##. All fecal samples were processed exactly in the same way and the count of cells hybridized with Eub338 can be calculated to give the absolute numbers of bacteria per gram of feces. No significant difference in concentration of total fecal bacteria was found between the healthy controls and FMF patients in remission, the mean values 5.47×10<sup>10</sup>/g and 5.12×10<sup>10</sup>/g, respectively (##TAB##4##Table 5##). In the acute stage of the disease, however, the concentration of bacteria in fecal material was substantially lower, at 3.56×10<sup>10</sup>/g.</p> <p>Among the taxonomic groups tested, the highest proportion was detected with the Bac303 probe, which targets the <italic>Bacteroides</italic> group. It gave the concentration of these bacteria at 1.51×10<sup>10</sup>/g (27.3% of total bacteria) in healthy controls, with elevation in FMF patients to 1.86×10<sup>10</sup>/g (39.1%) in remission and to 1.35×10<sup>10</sup>/g (38.6%) in the attack period (##TAB##4##Table 5##). Cluster XIVa bacteria (<italic>Clostridium coccoides</italic> group) was the second most abundant group enumerated with the Erec482 probe comprising 1.09×10<sup>10</sup>/g (21.1%) in the healthy controls and approximately 1.01×10<sup>10</sup>/g in both, FMF remission and attack phases (18% and 28.9% of total bacteria, respectively). The <italic>E. rectale-Roseburia</italic> probe Rrec584, which is a nested probe within the <italic>C. coccoides</italic> group (enumerated with Erec482), detected 6.5%, 8.3% and 18.0% within the total bacteria probe, Eub338, and 41.5%, 43.8% and 59.8% within the Erec482 probe in healthy controls, FMF remission and FMF attacks, respectively. There were significant populations of <italic>Fecalibacterium prausnitzii</italic> related bacteria (Fprau645), the second important group of butyrate producers within clostridial cluster IV, in all studied groups, with the tendency to increase in acute FMF. These bacteria comprised 0.54×10<sup>10</sup>/g (11.3%) in healthy controls, 0.5×10<sup>10</sup>/g (10.4%) in FMF remissions, and 0.53×10<sup>10</sup>/g (13.8%) in acute stage of FMF. Another group within clostridial cluster IV, ruminococci, which was detected with the combination of probes Rfla729 and Rbro730, demonstrated lower levels in FMF patients, especially in the acute phase of the disease, compared to healthy controls. In particular, the members of the <italic>R. flavefaciens</italic>/<italic>R. bromii</italic> group accounted for 0.59×10<sup>10</sup>/g (9.9%) in healthy controls, 0.29×10<sup>10</sup>/g (5.3%) in FMF remissions, and 0.09×10<sup>10</sup>/g (2.6%) in acute FMF. Similarly, clostridial cluster IX representatives, detected with the probe Prop853, were estimated to be at 0.64×10<sup>10</sup>/g (11%) in healthy group, whereas in FMF remission these bacteria were reduced to 0.44×10<sup>10</sup>/g (7.1%), with even more significant reduction in the acute stage of FMF, accounting for 0.13×10<sup>10</sup>/g (3.5%). <italic>Atopobium</italic> group also tended to be less represented in gut microflora of FMF patients compared with healthy subjects, steadily decreasing from 0.42×10<sup>10</sup>/g (7.8%) in healthy controls to 0.29×10<sup>10</sup>/g (6.3%) in remission and 0.17×10<sup>10</sup>/g (4.5%) in attack. Members of the <italic>Bifidobacterium</italic>, <italic>Enterobacteriaceae</italic> and <italic>Lactobacillus-Enterococcus</italic> groups, which were detected by the probes Bif164, EnterobacD and Lab158, respectively, did not reveal significant differences between the studied groups because they were found in low numbers and demonstrated high inter-individual variability (##TAB##4##Table 5##).</p> <p>We also attempted to evaluate the difference in gut microbiota composition in remission and attack periods of disease in the “pure” form, when the genetic, phenotypic and environmental variables are removed from the analysis. For this, we monitored the intra-individual shifts in total and main gut bacterial populations during the remission and attack phases of the disease in two patients and these observations confirmed and reinforced the tendencies found between the larger combined cohorts of FMF sufferers in the remission and acute stages (##TAB##5##Table 6##). First, the total numbers of bacteria in the gut enumerated with the Eub338 probe decreased about two-fold. Second, in a number of bacterial groups, the tendencies detected earlier became more profound and statistically significant. In particular, these changes included a dramatic increase in concentration of cluster XIVa bacteria and particularly of <italic>E. rectale</italic>-<italic>Roseburia</italic> cluster during the acute phase of the disease, with the reduction in concentration of cluster IX bacteria and ruminococci. The fall of proportion of populations of lactic acid bacteria enumerated with the Bif164 and Lab158 probes during the disease attack was significant in one patient, FMF 13-123, while in the other, FMF 87, the differences in this group were not significant. Thus, the monitoring of intra-individual bacterial shifts during the remission and attack phases validated and confirmed the findings involving the larger cohorts (##TAB##4##Tables 5## and ##TAB##5##6##).</p> <p>The part of FMF patients was undergoing colchicine therapy and we compared the composition of gut bacteria in respect to this treatment (##TAB##6##Table 7##). There was no difference in acute disease but in remission a certain tendency to normalization was observed in the <italic>Bacteroides</italic> group and in clostridial cluster IX bacteria. Nevertheless, these differences were not statistically significant and the data sets of colchicine-treated and untreated patients were combined together for the microbiota analyses.</p> <p>Multiple comparative analyses of relative proportions of bacterial populations in the gut of the three cohorts were performed using the Kruskal-Wallis test and the results are presented in ##FIG##2##Fig. 3##. The proportion of <italic>Bacteroides</italic> was significantly higher in both diseased cohorts in comparison with the healthy controls, while the proportion of cluster XIVa bacteria was only increased in the acute phase of the disease in comparison with the control group and remission (##FIG##2##Fig. 3A##). The fraction of <italic>E. rectale-Roseburia</italic> cluster among total bacteria was substantially increased, especially in active disease (##FIG##2##Fig. 3B##), while bacteria of cluster IX and <italic>Atopobium</italic> decreased in the attack phase of the disease in comparison with controls and remission (##FIG##2##Fig. 3C##). The proportion of ruminococci was significantly lower in FMF patients (##FIG##2##Fig. 3D##), while no significant differences were detected among <italic>F. prausnitzii</italic>, bifidobacteria, lactic acid bacteria, and enterobacteria (##FIG##2##Fig. 3B, 3D and 3E##).</p> <p>Discriminant function analysis (DA) is a powerful tool to extract specific patterns in multiple datasets and we used this approach to analyze the specific gut microbiota distribution patterns in our three cohorts. The summary of discriminant function analysis with 20 subjects and 10 variables is shown in ##FIG##3##Fig. 4##. It revealed three clusters in bacterial distribution that are highly specific, well separated and distinct for healthy controls, FMF remission, and FMF attack. As shown in ##TAB##7##Table 8##, the total percentage of correct classification was 95%, thus 95% of total subjects were classified into the correct groups (100% in healthy controls, 83.3% in remission and 100% in attack). The validity of the variables in differentiating the groups was supported by the Wilks's lambda coefficient, which may vary from 0 to 1 (the smaller the lambda, the more the variable differentiates the groups), and in this case it had a value of 0.045. The F-test of Wilks's lambda, which is used to test if the discriminant model as a whole is significant, produced (20,16) = 2.95; p&lt;0.016. There was only one case of misclassification from a total of 20 when the remission stage was assessed as acute. This patient refused the colchicine therapy and was suffering from very frequent FMF attacks and there is a possibility that at the time of sampling the microbiota was not restructured back to the composition specific for the remission stage. Thus the model displayed good predictive accuracy and suggested the specific patterns of gut microbiota distribution in the healthy and diseased states.</p>
<title>Discussion</title> <p>This is the first work reporting the specific changes in the human gut microbiota due to the mutations in a single gene of the host. Suggestions that the composition of gut microbiota is host genotype-dependant have been proposed earlier ##UREF##4##[43]##, ##REF##16278440##[44]##, however, these works were based on twin studies and the host genotype was not assessed. Gut inflammatory diseases such as IBD have been long known as having a substantial genetic predisposition component and a recent review of genome-wide scans reported on the seven confirmed and replicated cases of genetic susceptibility loci, with many more, however, that need independent confirmation and verification ##REF##17894847##[42]##. Although the gut microbiota of IBD patients was subjected to several diversity analyses ##REF##16188921##[35]##, ##REF##17699621##[36]##, ##REF##12524406##[45]##, the host genetics was not assessed, most probably because of complicated genetics of IBD. Thus, uncovering the genetic mechanisms that are involved in regulation of composition of human gut microbiota remains a challenging task. Evidence obtained in a study of bacterium-derived cellular fatty acids in the stool samples of six mouse strains congenic for the major histocompatibility complex (MHC) suggested the importance of host immune system in this process ##REF##11254595##[46]##. The obvious advantage of our model is that the genetics of FMF is well defined and confined to a single gene thus facilitating the host genotype analysis. Thus, for the first time, we had the opportunity for direct comparison of the host and gut community genotypes.</p> <p>Initially, we genotyped the FMF patients and healthy controls for the presence of disease-specific mutations in exons 2 and 10 of the <italic>MEFV</italic> gene. In all clinically confirmed cases we were able to detect the presence of known mutations, all of them in exon 10, which encodes the SPRY/B30.2 domain. Among our very limited cohort of seven healthy controls of Armenian origin, one heterozygote was detected, suggesting that the mutation carrier rate in this population is at least 14%, which is within the range of carrier rates found in populations of other ethnic origin ##UREF##5##[47]##–##REF##11464248##[50]##. The high carrier frequency in affected populations suggests the selective advantage conferred by the heterozygous state but no perceptible biological advantage of the carrier state is presently identified ##REF##11484206##[49]##. Positive selection to maintain the high frequency of the heterozygotes should also overcome the negative selection imposed by the increased morbidity and mortality rates among the homozygotes and compounded heterozygotes. Clinical manifestations during the disease attack are consistent with a strong inflammatory response, with fever, massive influx of polymorphonuclear leukocytes into the affected tissues, neutrophilia, and acute-phase response. In some patients this leads to systemic amyloidosis resulting in renal failure and death.</p> <p>The penetrance of mutations that are all missense mutations, varied widely. Apart from the combined heterozygotes, which support the recessive nature of mutations, there were also simple heterozygotes but with different phenotypes. The simple heterozygosity in clinically confirmed disease was exclusively exemplified by a M694V mutation and, similarly to another study ##REF##10787449##[51]##, suggested the autosomal dominant inheritance probably due to the importance of this residue in the interaction of pyrin with caspase-1 ##REF##16785446##[52]##. On the other end of the mutation spectrum is the simple M680I heterozygosity found in one of the healthy volunteers, C6b, which did not result in any clinical or subclinical (cytokine profile) FMF manifestation. Thus, depending on the locus involved, the outcome could be either an asymptomatic carrier state or a full-blown disease. In the case of M694V mutation, a 50% complementation by the normal allele may be not sufficient to provide the protection against clinical manifestations of the disease. This is the most abundant mutation in our FMF cohort and is present in 70% of the patients. Taking into the consideration the high penetrance and associated fitness cost of carrying this allele even in the heterozygous form it seems counterintuitive to suggest that it may confer any selective advantage in the heterozygous carrier state under the present day living conditions. Phenotypically, the asymptomatic heterozygous carriers of FMF mutations are characterized by the elevated levels of markers of inflammation such as CRP, SAA, and INF-γ ##UREF##6##[53]##, ##REF##12966607##[54]## and by specific glycosylation of AGP ##REF##11454642##[55]##. Whether this subclinical inflammation in asymptomatic patients may be protective against specific pathogens, probably endemic for the Mediterranean basin, remains to be seen since the previous studies failed to identify such agents ##REF##11484206##[49]##.</p> <p>In a parallel study, we measured the levels of systemic markers of inflammation in acute phase and remission of FMF ##REF##18440310##[56]##. The study confirmed earlier observations that even in remission, where the clinical signs of disease are absent, the level of certain markers of inflammation remains elevated suggesting chronic subclinical inflammation in FMF. How this inflammation may affect the gut microbiota? To answer this question, we used two independent techniques to assess gut molecular bacterial diversity in health and disease, the sequence analysis of 16S rRNA gene libraries and FISH. Both these techniques have their own advantages and disadvantages but in combination they are complementary and allow obtaining a more balanced overview of gut microbiota. In particular, 16S libraries may suffer from PCR biases introduced by elevated number of cycles ##REF##12054235##[57]## or poor amplification of certain templates ##REF##11375193##[58]## while FISH analysis is less sensitive since it targets higher taxonomic entities and no community analysis involving individuals and species/phylotypes, unless they very are numerous, is possible. On the other hand, FISH technique consistently detects higher levels of <italic>Bifidobacteria</italic> than clone library sequencing, it is less expensive and can be easily automated in conjunction with flow cytometry or computer-aided automatic count under the microscope for analysis of large cohorts. In order to minimize inter-individual variations, 16S rRNA gene libraries were constructed from fecal samples of 12 FMF patients in remission, three patients in the acute stage of FMF, six healthy Armenians, and one healthy subject of non-Armenian origin. The number of amplification cycles was maintained low (10 cycles) to prevent a potential PCR bias.</p> <p>Both analyses of gut microbiota produced the similar results in terms of general distribution of major phyla representatives in the human gut. Similarly to the previous studies ##REF##17183309##[8]##, ##REF##17699621##[36]##, ##REF##15831718##[59]## the vast majority of bacterial sequences and cells hybridized to specific probes in the gut fell into the <italic>Bacteroidetes</italic> and <italic>Firmicutes</italic> phyla. The <italic>Bacteroidetes</italic> contained the largest number of sequences but fewer phylotypes while the diverse <italic>Firmicutes</italic> were dominant in terms of the numbers and the relative proportion of phylotypes detected. The <italic>Bacteroidetes</italic> was also the largest group detected by the Bac303 probe accounting for 27–39% of the total bacterial count, which is lower than the clone library estimates. On the contrary, the combination of various FISH probes targeting the <italic>Firmicutes</italic> produced a higher estimate of bacteria belonging to this phylum than the clone libraries. The most profound difference between the two techniques in detection of bacteria in the human gut, however, concerned the enumeration of high G+C Gram-positive bacteria such as <italic>Atopobium</italic> spp. and <italic>Bifidobacterium</italic> spp.: FISH technique consistently detected a sizable population of these genera in the gut while they were hardly represented in PCR-generated libraries. Taking into the consideration an important role played by, for example, <italic>Bifidobacteria</italic> in maintaining the gastrointestinal health, consistent underestimates of these bacteria in clone library analyses may obscure the effect of this group in various gut diseases such as IBD as well as in the healthy gut.</p> <p>Examination of gut bacterial diversity in healthy and diseased cohorts by the two methods arrived essentially at the same conclusion suggesting the specific shifts in composition of bacteria imposed by mutations in the <italic>MEFV</italic> gene of the host. The gut microbiota in active disease was characterized by diminished bacterial diversity and major population restructuring within the phyla. In particular, the proportion of phylotypes belonging to of <italic>Porphyromonadaceae</italic>, <italic>Phascolarctobacterium</italic>, <italic>Faecalibacterium</italic>, and <italic>Parabacteroides</italic> were significantly increased while <italic>Prevotellaceae, Dialister</italic> and <italic>Prevotella</italic> were significantly lower in acute disease in comparison with control. Both approaches detected the fall in proportion of clostridial cluster IX in FMF attack: FISH demonstrated it for the whole group, while the sequence libraries elaborated this general trend as involving <italic>Acidaminococcaceae, Megasphaera,</italic> and <italic>Dialister.</italic> Evenness values for microbiota in active disease were higher reflecting the lower number of dominant phylotypes, although this could be a consequence of the smaller number of sequences collected. Also, in a parallel FISH approach, the fall in absolute numbers of bacteria during the inflammation attacks was detected. In most cases, this analysis is less sensitive than sequencing because it targets the larger taxonomic entities such as, for example, the whole <italic>Bacteroides</italic> group and cannot detect the differences at the lower taxa levels, e.g., in the <italic>Prevotellaceae</italic> and <italic>Porphyromonadaceae</italic> families and the <italic>Prevotella</italic> and <italic>Parabacteroides</italic> genera while sequence analysis operating at the individual sequence level is capable to do. In some cases, however, the situation may be opposite: Rfla729/Rbro730 FISH probe, for example, detects a subset of ruminococci, which is significantly reduced in attack (##TAB##4##Table 5##), while at the genus level the differences are not significant (##TAB##2##Table 3##). Thus, FISH served as a complementary tool and allowed analyzing bacterial populations that are routinely underestimated in sequence analyses. In particular, FISH detected the fall in absolute numbers of high G+C bacteria such as <italic>Bifidobacterium</italic> and <italic>Atopobium</italic>, which are not represented well in our and others SSU rRNA clone libraries ##REF##15831718##[59]##. Thus, in disease attack, the drop of total bacterial count and restructuring of bacterial composition in the gut may be the consequence of acute inflammation, which is accompanied by fever and influx of polymorphonuclear neutrophils (PMNs) into the affected areas including peritoneum. Elevated body temperatures (up to 40°C) and the localized respiratory burst of PMNs producing reactive oxygen species could be responsible for gut microbiota selection.</p> <p>As mentioned earlier, despite the absence of any clinical signs of disease in remission, certain systemic markers of inflammation are still elevated suggesting subclinical chronic inflammation ##REF##18440310##[56]##. Interestingly, the diversity values in remission were higher than in the healthy gut suggesting that the low-level subclinical inflammation in FMF may favour a greater diversity of intestinal bacteria. DA analysis of FISH data also confirmed the specific composition of bacterial populations in remission, which is different from control and forms a separate cluster. An additional internal control subject of different ethnic background was included in our study and the classification of the subject was consistently in the healthy group. The specific shifts of bacterial composition were also confirmed in two patients that were investigated during the remission and attack periods to exclude the influence of other genetic or environmental factors and these analyses confirmed the specific bacterial composition shifts due to disease state such as the increase of butyrate producers and the fall of clostridial cluster IX bacteria during disease attacks. While these shifts in attack may be explained by factors such as acute inflammation, with fever and PMN infiltration, the factors that govern the specific bacterial shifts in remission, which is free of disease symptoms, are harder to suggest. Nevertheless, our DA model consistently classified the patients in remission as a separate group, which suggests the specific patterns of gut microbiota in this cohort.</p> <p>Thus the question is: what are the molecular mechanisms responsible for these microbiota shifts imposed by the mutated version of pyrin? In a broad sense, the functional role of this protein is the regulation of innate immunity ##UREF##3##[38]##. The protein is expressed in granulocytes, activated monocytes, and serosal and synovial fibroblasts ##REF##10807793##[60]##–##REF##15529356##[62]##. Recently it becomes evident that the regulatory functions of pyrin are carried out through the interaction of two domains of the protein, PYRIN and B30.2/rfp/SPRY, with the proinflammatory caspase-1 activating complex, called the inflammasome. The N-terminal PYRIN domain of the protein binds and competes for ASC (apoptosis-associated speck-like protein containing a caspase-recruitment domain), an inflammasome component, thereby decreasing its availability for the cryopyrin/NALP3/CIAS1 inflammasome. This results in reduction of caspase-1 activation and pro-interleukin-1beta processing and secretion thus halting the inflammatory response ##REF##12667444##[63]##. The mutational hotspot in FMF, however, is located in the C-terminal B30.2/rfp/SPRY domain, which modulates the inflammatory response through the interaction with several components of the inflammasome. In addition to interaction with cryopyrin/NALP3/CIAS1, it interacts directly with caspase-1 and its substrate pro-IL-1β ##REF##16785446##[52]##, ##REF##17431422##[64]##. The net effect of these interactions is the suppression of IL-1β activation and block of inflammation. It was suggested that the heightened IL-1β responsiveness could be a factor selecting for mutations in the MEFV gene ##REF##16785446##[52]##. Taking into consideration the overall role of pyrin as a modulator/suppressor of the inflammatory response, another factor contributing to the autoinflammatory nature of the disease could be the reduced <italic>MEFV</italic> messenger RNA expression in patients with familial Mediterranean fever ##REF##12384939##[65]##. First, it may contribute to the lower titre of pyrin and its PYRIN domain molecules in the cell thus making more ASC molecules available to initiate caspase-1 activation. Second, the reduced concentration of pyrin and therefore of its B30.2/rfp/SPRY domain, which, in addition, is mutated to the loss of caspase-1 suppressor function in most of FMF cases, may provoke easier triggering the inflammation cascade through caspase-1 activation. Thus, both consequences of <italic>MEFV</italic> mutations may lead to the heightened responsiveness of cryopyrin/NALP3/CIAS1, which can be oligomerized and activated in response to a very diverse range of ligands such as bacterial muramyl dipeptide, ATP, toxins, bacterial and viral RNA, small antiviral compounds, <italic>Staphylococcus, Listeria</italic>, and uric acid crystals as well as by low intracellualr potassium concentration ##REF##15530394##[66]##–##REF##17599094##[71]##. While these exo- and endogenous stimuli are normally not perceived as danger signals or are efficiently suppressed, the facilitated oligomerization of the cryopyrin/NALP3/CIAS1, which is not adequately suppressed by the mutated pyrin in FMF, may be the underlying cause in this disease, resulting in excessive pro-caspase-1 and pro-IL1β processing. Indeed, monocytes from FMF patients in remission fail to induce LPS homologous tolerance and exhibit heightened sensitivity to bacterial endotoxin ##REF##16785765##[72]##, one of the important activators of the cryopyrin/NALP3/CIAS1 inflammasome ##REF##17186029##[73]##.</p> <p>In conclusion, we established in this work that mutations in a single host gene lead to specific restructuring of commensal gut microbiota. Although the exact molecular mechanisms behind this process are not fully understood, at our current genetic stage of analysis the work points to the host genotype as a key factor in symbiotic host-microbe interaction. Further developments in this area may focus on details of this interaction at the biochemical and physiological levels. Other models of gut inflammation such as IBD may help to establish the role of other host genes affecting the composition of commensal gut microbiota. Comparative community analyses from genetically different hosts are particularly interesting in identifying gut microbiota alterations specific for a particular disease genotype. At the same time, it cannot be excluded that the changes in the community structure may share some common features reflecting the similarities in disease phenotypes such as inflammation. And finally, how the gut microbiota is transformed in disease to become a “pathobiota” ##REF##17923086##[37]## remains an interesting question to answer.</p>
<p>Conceived and designed the experiments: ZK RA. Performed the experiments: ZK GM KG. Analyzed the data: ZK GM KG RA. Contributed reagents/materials/analysis tools: DK RA. Wrote the paper: ZK RA.</p> <title>Background</title> <p>The human gastrointestinal tract is inhabited by a very diverse symbiotic microbiota, the composition of which depends on host genetics and the environment. Several studies suggested that the host genetics may influence the composition of gut microbiota but no genes involved in host control were proposed. We investigated the effects of the wild type and mutated alleles of the gene, which encodes the protein called pyrin, one of the regulators of innate immunity, on the composition of gut commensal bacteria. Mutations in <italic>MEFV</italic> lead to the autoinflammatory disorder, familial Mediterranean fever (FMF, MIM249100), which is characterized by recurrent self-resolving attacks of fever and polyserositis, with no clinical signs of disease in remission.</p> <title>Methodology/Principal Findings</title> <p>A total of 19 FMF patients and eight healthy individuals were genotyped for mutations in the <italic>MEFV</italic> gene and gut bacterial diversity was assessed by sequencing 16S rRNA gene libraries and FISH analysis. These analyses demonstrated significant changes in bacterial community structure in FMF characterized by depletion of total numbers of bacteria, loss of diversity, and major shifts in bacterial populations within the <italic>Bacteroidetes, Firmicutes</italic> and <italic>Proteobacteria</italic> phyla in attack. In remission with no clinical signs of disease, bacterial diversity values were comparable with control but still, the bacterial composition was substantially deviant from the norm. Discriminant function analyses of gut bacterial diversity revealed highly specific, well-separated and distinct grouping, which depended on the allele carrier status of the host.</p> <title>Conclusions/Significance</title> <p>This is the first report that clearly establishes the link between the host genotype and the corresponding shifts in the gut microbiota (the latter confirmed by two independent techniques). It suggests that the host genetics is a key factor in host-microbe interaction determining a specific profile of commensal microbiota in the human gut.</p>
<p>We thank Sylvia Duncan, Alan Walker and Harry Flint for help with automated FISH.</p>
<fig id="pone-0003064-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003064.g001</object-id><label>Figure 1</label><caption><title>Rarefaction curves obtained for clone libraries from healthy controls and FMF patients.</title><p>The rarefaction curves were generated by plotting the number of phylotypes (defined at 99% level) against the number of clones sequenced. The shape of the curves of observed phylotype richness indicates a trend of diminishing chance of finding new phylotypes as the sampling continues. Data for library in attack were extrapolated by empiric regression function a×x<sup>2</sup>b×exp(−c×x). Coefficient of regression determination (R<sup>2</sup>) was 99.9%. Regression statistics: adjusted R<sup>2</sup> – 0.9991; standard error – 0.1720; a = 1.8690; b = 0.7104; c = 8.0466E-04. Healthy controls are shown in dark blue; FMF patients in remission – in red and attack – in green, the extrapolated curve is shown in light blue.</p></caption></fig> <fig id="pone-0003064-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003064.g002</object-id><label>Figure 2</label><caption><title>Collector's curves of the observed and estimated phylotype richness of 16S libraries.</title><p>Collector's curves of the observed (shown in dark blue) and estimated (Chao1 (light blue) and ACE (red)) phylotype richness calculated for healthy subjects (A), FMF remission (B), and FMF attack (C) with the phylogenetic depth defined at 99% level.</p></caption></fig> <fig id="pone-0003064-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003064.g003</object-id><label>Figure 3</label><caption><title>Comparison of the main bacterial groups identified by FISH between healthy and diseased groups.</title><p>Comparative analysis of relative proportions of bacterial groups detected by different oligonucleotide probes in fecal samples between FMF patients in remission and attack periods and healthy control group (by Kruskal-Wallis test for multiple comparisons). On box-and-whisker plots the points within boxes indicate the mean values of bacterial proportions, the boxes represent standard errors (SE) and the vertical bars represent standard deviations (SD).</p></caption></fig> <fig id="pone-0003064-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003064.g004</object-id><label>Figure 4</label><caption><title>Discriminant function analysis of FISH data.</title><p>Number of variables in model: 10; grouping - 3 groups. Wilks's lambda coefficient – 0.045; F-test - (20,16) = 2.95, p&lt;0.016. Root 1, 2 – discriminant functions 1 and 2 (1<sup>st</sup> and 2<sup>nd</sup> canonical roots). Regions of group distribution are outlined by ellipses; healthy is shown in blue triangles, remission – in red circles, and attack – in green squares.</p></caption></fig>
<table-wrap id="pone-0003064-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003064.t001</object-id><label>Table 1</label><caption><title>
<italic>MEFV</italic> genotypes of FMF patients and control subjects and gut microbiota analysis methods applied.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Subjects</td><td align="left" rowspan="1" colspan="1">Exon 10 FMF mutations</td><td colspan="2" align="left" rowspan="1">Exon 2 non FMF mutations</td><td align="left" rowspan="1" colspan="1">Gut microbiota analyses</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">D102D</td><td align="left" rowspan="1" colspan="1">G138G</td><td align="left" rowspan="1" colspan="1"/></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">C1b</td><td align="left" rowspan="1" colspan="1">N/D<xref ref-type="table-fn" rid="nt103">*</xref>
</td><td align="left" rowspan="1" colspan="1">HT</td><td align="left" rowspan="1" colspan="1">HT</td><td align="left" rowspan="1" colspan="1">S</td></tr><tr><td align="left" rowspan="1" colspan="1">C4b</td><td align="left" rowspan="1" colspan="1">N/D</td><td align="left" rowspan="1" colspan="1">N/D</td><td align="left" rowspan="1" colspan="1">N/D</td><td align="left" rowspan="1" colspan="1">S; F</td></tr><tr><td align="left" rowspan="1" colspan="1">C6b</td><td align="left" rowspan="1" colspan="1">M680I</td><td align="left" rowspan="1" colspan="1">HT</td><td align="left" rowspan="1" colspan="1">HT</td><td align="left" rowspan="1" colspan="1">S; F</td></tr><tr><td align="left" rowspan="1" colspan="1">C11b</td><td align="left" rowspan="1" colspan="1">N/D</td><td align="left" rowspan="1" colspan="1">N/D</td><td align="left" rowspan="1" colspan="1">N/D</td><td align="left" rowspan="1" colspan="1">S; F</td></tr><tr><td align="left" rowspan="1" colspan="1">C15b</td><td align="left" rowspan="1" colspan="1">N/D</td><td align="left" rowspan="1" colspan="1">HT</td><td align="left" rowspan="1" colspan="1">HT</td><td align="left" rowspan="1" colspan="1">S; F</td></tr><tr><td align="left" rowspan="1" colspan="1">C16b</td><td align="left" rowspan="1" colspan="1">N/D</td><td align="left" rowspan="1" colspan="1">N/D</td><td align="left" rowspan="1" colspan="1">N/D</td><td align="left" rowspan="1" colspan="1">F</td></tr><tr><td align="left" rowspan="1" colspan="1">C18b</td><td align="left" rowspan="1" colspan="1">N/D</td><td align="left" rowspan="1" colspan="1">HT</td><td align="left" rowspan="1" colspan="1">HT</td><td align="left" rowspan="1" colspan="1">S; F</td></tr><tr><td align="left" rowspan="1" colspan="1">Cm13b</td><td align="left" rowspan="1" colspan="1">N/D</td><td align="left" rowspan="1" colspan="1">N/D</td><td align="left" rowspan="1" colspan="1">N/D</td><td align="left" rowspan="1" colspan="1">S; F</td></tr><tr><td align="left" rowspan="1" colspan="1">FMF 0</td><td align="left" rowspan="1" colspan="1">M680I/V726A</td><td align="left" rowspan="1" colspan="1">N/A<xref ref-type="table-fn" rid="nt104">**</xref>
</td><td align="left" rowspan="1" colspan="1">N/A</td><td align="left" rowspan="1" colspan="1">SR</td></tr><tr><td align="left" rowspan="1" colspan="1">FMF 1</td><td align="left" rowspan="1" colspan="1">M694V/A761H</td><td align="left" rowspan="1" colspan="1">N/A</td><td align="left" rowspan="1" colspan="1">N/A</td><td align="left" rowspan="1" colspan="1">SR</td></tr><tr><td align="left" rowspan="1" colspan="1">FMF 2</td><td align="left" rowspan="1" colspan="1">M694V/A761H</td><td align="left" rowspan="1" colspan="1">N/A</td><td align="left" rowspan="1" colspan="1">N/A</td><td align="left" rowspan="1" colspan="1">SR</td></tr><tr><td align="left" rowspan="1" colspan="1">FMF 3</td><td align="left" rowspan="1" colspan="1">M694V/A761H</td><td align="left" rowspan="1" colspan="1">N/A</td><td align="left" rowspan="1" colspan="1">N/A</td><td align="left" rowspan="1" colspan="1">SR</td></tr><tr><td align="left" rowspan="1" colspan="1">FMF 5</td><td align="left" rowspan="1" colspan="1">M680I/M694V</td><td align="left" rowspan="1" colspan="1">N/A</td><td align="left" rowspan="1" colspan="1">N/A</td><td align="left" rowspan="1" colspan="1">SR</td></tr><tr><td align="left" rowspan="1" colspan="1">FMF 7</td><td align="left" rowspan="1" colspan="1">M694V</td><td align="left" rowspan="1" colspan="1">N/D</td><td align="left" rowspan="1" colspan="1">N/D</td><td align="left" rowspan="1" colspan="1">SR</td></tr><tr><td align="left" rowspan="1" colspan="1">FMF 8</td><td align="left" rowspan="1" colspan="1">M694V</td><td align="left" rowspan="1" colspan="1">HM</td><td align="left" rowspan="1" colspan="1">HM</td><td align="left" rowspan="1" colspan="1">SR</td></tr><tr><td align="left" rowspan="1" colspan="1">FMF 10</td><td align="left" rowspan="1" colspan="1">M694V</td><td align="left" rowspan="1" colspan="1">HT</td><td align="left" rowspan="1" colspan="1">HT</td><td align="left" rowspan="1" colspan="1">SR; FR</td></tr><tr><td align="left" rowspan="1" colspan="1">FMF 11</td><td align="left" rowspan="1" colspan="1">M694V</td><td align="left" rowspan="1" colspan="1">HT</td><td align="left" rowspan="1" colspan="1">HT</td><td align="left" rowspan="1" colspan="1">SR</td></tr><tr><td align="left" rowspan="1" colspan="1">FMF 12</td><td align="left" rowspan="1" colspan="1">M694V</td><td align="left" rowspan="1" colspan="1">HM</td><td align="left" rowspan="1" colspan="1">HM</td><td align="left" rowspan="1" colspan="1">SR; FR</td></tr><tr><td align="left" rowspan="1" colspan="1">FMF 13-123</td><td align="left" rowspan="1" colspan="1">M680I/V726A</td><td align="left" rowspan="1" colspan="1">N/A</td><td align="left" rowspan="1" colspan="1">N/A</td><td align="left" rowspan="1" colspan="1">SA; SR; FA; FR</td></tr><tr><td align="left" rowspan="1" colspan="1">FMF 14</td><td align="left" rowspan="1" colspan="1">M694V/A761H</td><td align="left" rowspan="1" colspan="1">N/A</td><td align="left" rowspan="1" colspan="1">N/A</td><td align="left" rowspan="1" colspan="1">SR; FR</td></tr><tr><td align="left" rowspan="1" colspan="1">FMF 87</td><td align="left" rowspan="1" colspan="1">M694V/A761H</td><td align="left" rowspan="1" colspan="1">HT</td><td align="left" rowspan="1" colspan="1">HT</td><td align="left" rowspan="1" colspan="1">FR; FA</td></tr><tr><td align="left" rowspan="1" colspan="1">FMF 103</td><td align="left" rowspan="1" colspan="1">M680I/M680I</td><td align="left" rowspan="1" colspan="1">N/A</td><td align="left" rowspan="1" colspan="1">N/A</td><td align="left" rowspan="1" colspan="1">SA; FA</td></tr><tr><td align="left" rowspan="1" colspan="1">FMF 109</td><td align="left" rowspan="1" colspan="1">M694V</td><td align="left" rowspan="1" colspan="1">HT</td><td align="left" rowspan="1" colspan="1">HT</td><td align="left" rowspan="1" colspan="1">FA</td></tr><tr><td align="left" rowspan="1" colspan="1">FMF 110</td><td align="left" rowspan="1" colspan="1">M694V/V726A</td><td align="left" rowspan="1" colspan="1">N/D</td><td align="left" rowspan="1" colspan="1">N/D</td><td align="left" rowspan="1" colspan="1">FA</td></tr><tr><td align="left" rowspan="1" colspan="1">FMF 111</td><td align="left" rowspan="1" colspan="1">M680I/M694V</td><td align="left" rowspan="1" colspan="1">HM</td><td align="left" rowspan="1" colspan="1">HM</td><td align="left" rowspan="1" colspan="1">FA</td></tr><tr><td align="left" rowspan="1" colspan="1">FMF 114</td><td align="left" rowspan="1" colspan="1">M680I/M680I</td><td align="left" rowspan="1" colspan="1">N/A</td><td align="left" rowspan="1" colspan="1">N/A</td><td align="left" rowspan="1" colspan="1">FR</td></tr><tr><td align="left" rowspan="1" colspan="1">FMF 124</td><td align="left" rowspan="1" colspan="1">M680I/M680I</td><td align="left" rowspan="1" colspan="1">N/A</td><td align="left" rowspan="1" colspan="1">N/A</td><td align="left" rowspan="1" colspan="1">SA; FA</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0003064-t002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003064.t002</object-id><label>Table 2</label><caption><title>Phylogenetic distribution of 16S rRNA gene sequences and phylotypes retrieved from healthy individuals and FMF patients in remission and attack.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Phylum</td><td colspan="2" align="left" rowspan="1">Healthy</td><td colspan="2" align="left" rowspan="1">FMF remission</td><td colspan="2" align="left" rowspan="1">FMF attack</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">OTUs<xref ref-type="table-fn" rid="nt105">*</xref>
</td><td align="left" rowspan="1" colspan="1">Clones</td><td align="left" rowspan="1" colspan="1">OTUs</td><td align="left" rowspan="1" colspan="1">Clones</td><td align="left" rowspan="1" colspan="1">OTUs</td><td align="left" rowspan="1" colspan="1">Clones</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">
<italic>Bacteroidetes</italic>
</td><td align="left" rowspan="1" colspan="1">54</td><td align="left" rowspan="1" colspan="1">324</td><td align="left" rowspan="1" colspan="1">57</td><td align="left" rowspan="1" colspan="1">319</td><td align="left" rowspan="1" colspan="1">21</td><td align="left" rowspan="1" colspan="1">71</td></tr><tr><td align="left" rowspan="1" colspan="1">% <italic>Bacteroidetes</italic>
</td><td align="left" rowspan="1" colspan="1">40.91</td><td align="left" rowspan="1" colspan="1">56.64</td><td align="left" rowspan="1" colspan="1">35.19</td><td align="left" rowspan="1" colspan="1">50.72</td><td align="left" rowspan="1" colspan="1">40.38</td><td align="left" rowspan="1" colspan="1">55.91</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>Firmicutes</italic>
</td><td align="left" rowspan="1" colspan="1">73</td><td align="left" rowspan="1" colspan="1">237</td><td align="left" rowspan="1" colspan="1">98</td><td align="left" rowspan="1" colspan="1">289</td><td align="left" rowspan="1" colspan="1">29</td><td align="left" rowspan="1" colspan="1">54</td></tr><tr><td align="left" rowspan="1" colspan="1">% <italic>Firmicutes</italic>
</td><td align="left" rowspan="1" colspan="1">55.30</td><td align="left" rowspan="1" colspan="1">41.43</td><td align="left" rowspan="1" colspan="1">60.49</td><td align="left" rowspan="1" colspan="1">45.95</td><td align="left" rowspan="1" colspan="1">55.77</td><td align="left" rowspan="1" colspan="1">42.52</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>Proteobacteria</italic>
</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">11</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">17</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">2</td></tr><tr><td align="left" rowspan="1" colspan="1">% <italic>Proteobacteria</italic>
</td><td align="left" rowspan="1" colspan="1">3.79</td><td align="left" rowspan="1" colspan="1">1.92</td><td align="left" rowspan="1" colspan="1">3.09</td><td align="left" rowspan="1" colspan="1">2.70</td><td align="left" rowspan="1" colspan="1">3.85</td><td align="left" rowspan="1" colspan="1">1.57</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>Actinobacteria</italic>
</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">% <italic>Actinobacteria</italic>
</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">1.23</td><td align="left" rowspan="1" colspan="1">0.63</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">Total</td><td align="left" rowspan="1" colspan="1">132</td><td align="left" rowspan="1" colspan="1">572</td><td align="left" rowspan="1" colspan="1">162</td><td align="left" rowspan="1" colspan="1">629</td><td align="left" rowspan="1" colspan="1">52</td><td align="left" rowspan="1" colspan="1">127</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0003064-t003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003064.t003</object-id><label>Table 3</label><caption><title>Comparison of 16S rRNA gene libraries derived from healthy controls and FMF patients in remission and attack.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Bacterial taxa</td><td align="left" rowspan="1" colspan="1">Control</td><td align="left" rowspan="1" colspan="1">Remission</td><td align="left" rowspan="1" colspan="1">Attack</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Class <italic>Gammaproteobacteria</italic>
<xref ref-type="table-fn" rid="nt106">*</xref>
</td><td align="left" rowspan="1" colspan="1">0.2%</td><td align="left" rowspan="1" colspan="1">2.1%</td><td align="left" rowspan="1" colspan="1">0%</td></tr><tr><td align="left" rowspan="1" colspan="1">Order <italic>Enterobacteriales</italic>
<xref ref-type="table-fn" rid="nt106">*</xref>
</td><td align="left" rowspan="1" colspan="1">0%</td><td align="left" rowspan="1" colspan="1">1.6%</td><td align="left" rowspan="1" colspan="1">0%</td></tr><tr><td align="left" rowspan="1" colspan="1">Family <italic>Enterobacteriaceae</italic>
<xref ref-type="table-fn" rid="nt106">*</xref>
</td><td align="left" rowspan="1" colspan="1">0%</td><td align="left" rowspan="1" colspan="1">1.6%</td><td align="left" rowspan="1" colspan="1">0%</td></tr><tr><td align="left" rowspan="1" colspan="1">Family <italic>Acidaminococcaceae</italic>
<xref ref-type="table-fn" rid="nt106">*</xref>
<sup>,</sup>
<xref ref-type="table-fn" rid="nt108">***</xref>
</td><td align="left" rowspan="1" colspan="1">9.2%</td><td align="left" rowspan="1" colspan="1">14.8%</td><td align="left" rowspan="1" colspan="1">3.1%</td></tr><tr><td align="left" rowspan="1" colspan="1">Family <italic>Prevotellaceae</italic>
<xref ref-type="table-fn" rid="nt107">**</xref>
</td><td align="left" rowspan="1" colspan="1">28.8%</td><td align="left" rowspan="1" colspan="1">24.5%</td><td align="left" rowspan="1" colspan="1">17.3%</td></tr><tr><td align="left" rowspan="1" colspan="1">Family <italic>Porphyromonadaceae</italic>
<xref ref-type="table-fn" rid="nt107">**</xref>
<sup>,</sup>
<xref ref-type="table-fn" rid="nt108">***</xref>
</td><td align="left" rowspan="1" colspan="1">2.3%</td><td align="left" rowspan="1" colspan="1">3.3%</td><td align="left" rowspan="1" colspan="1">9.4%</td></tr><tr><td align="left" rowspan="1" colspan="1">Genus <italic>Ruminococcus</italic>
<xref ref-type="table-fn" rid="nt106">*</xref>
</td><td align="left" rowspan="1" colspan="1">0.2%</td><td align="left" rowspan="1" colspan="1">2.8%</td><td align="left" rowspan="1" colspan="1">1.6%</td></tr><tr><td align="left" rowspan="1" colspan="1">Genus <italic>Roseburia</italic>
<xref ref-type="table-fn" rid="nt106">*</xref>
</td><td align="left" rowspan="1" colspan="1">2.8%</td><td align="left" rowspan="1" colspan="1">0.8%</td><td align="left" rowspan="1" colspan="1">2.4%</td></tr><tr><td align="left" rowspan="1" colspan="1">Genus <italic>Megasphaera</italic>
<xref ref-type="table-fn" rid="nt106">*</xref>
<sup>,</sup>
<xref ref-type="table-fn" rid="nt108">***</xref>
</td><td align="left" rowspan="1" colspan="1">0.2%</td><td align="left" rowspan="1" colspan="1">4.5%</td><td align="left" rowspan="1" colspan="1">0%</td></tr><tr><td align="left" rowspan="1" colspan="1">Genus <italic>Dialister</italic>
<xref ref-type="table-fn" rid="nt107">**</xref>
<sup>,</sup>
<xref ref-type="table-fn" rid="nt108">***</xref>
</td><td align="left" rowspan="1" colspan="1">7.9%</td><td align="left" rowspan="1" colspan="1">6.7%</td><td align="left" rowspan="1" colspan="1">0%</td></tr><tr><td align="left" rowspan="1" colspan="1">Genus <italic>Phascolarctobacterium</italic>
<xref ref-type="table-fn" rid="nt107">**</xref>
</td><td align="left" rowspan="1" colspan="1">0%</td><td align="left" rowspan="1" colspan="1">0.6%</td><td align="left" rowspan="1" colspan="1">3.1%</td></tr><tr><td align="left" rowspan="1" colspan="1">Genus <italic>Faecalibacterium</italic>
<xref ref-type="table-fn" rid="nt107">**</xref>
<sup>,</sup>
<xref ref-type="table-fn" rid="nt108">***</xref>
</td><td align="left" rowspan="1" colspan="1">6.5%</td><td align="left" rowspan="1" colspan="1">6.5%</td><td align="left" rowspan="1" colspan="1">14.2%</td></tr><tr><td align="left" rowspan="1" colspan="1">Genus <italic>Prevotella</italic>
<xref ref-type="table-fn" rid="nt107">**</xref>
</td><td align="left" rowspan="1" colspan="1">28.4%</td><td align="left" rowspan="1" colspan="1">23.1%</td><td align="left" rowspan="1" colspan="1">16.5%</td></tr><tr><td align="left" rowspan="1" colspan="1">Genus <italic>Parabacteroides</italic>
<xref ref-type="table-fn" rid="nt107">**</xref>
<sup>,</sup>
<xref ref-type="table-fn" rid="nt108">***</xref>
</td><td align="left" rowspan="1" colspan="1">1.2%</td><td align="left" rowspan="1" colspan="1">2.1%</td><td align="left" rowspan="1" colspan="1">5.9%</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0003064-t004" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003064.t004</object-id><label>Table 4</label><caption><title>Indices of diversity, richness and evenness and library coverage for 16S rRNA gene libraries.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Measurement</td><td align="left" rowspan="1" colspan="1">Healthy (n = 7)</td><td align="left" rowspan="1" colspan="1">FMF remission (n = 12)</td><td align="left" rowspan="1" colspan="1">FMF attack (n = 3)</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">No. of clones</td><td align="left" rowspan="1" colspan="1">572</td><td align="left" rowspan="1" colspan="1">629</td><td align="left" rowspan="1" colspan="1">127</td></tr><tr><td align="left" rowspan="1" colspan="1">No. of OTUs</td><td align="left" rowspan="1" colspan="1">132</td><td align="left" rowspan="1" colspan="1">162</td><td align="left" rowspan="1" colspan="1">52</td></tr><tr><td align="left" rowspan="1" colspan="1">Singletons</td><td align="left" rowspan="1" colspan="1">53</td><td align="left" rowspan="1" colspan="1">67</td><td align="left" rowspan="1" colspan="1">29</td></tr><tr><td align="left" rowspan="1" colspan="1">Doubletons</td><td align="left" rowspan="1" colspan="1">18</td><td align="left" rowspan="1" colspan="1">33</td><td align="left" rowspan="1" colspan="1">9</td></tr><tr><td align="left" rowspan="1" colspan="1">Chao1 estimator of species richness<xref ref-type="table-fn" rid="nt110">*</xref>
</td><td align="left" rowspan="1" colspan="1">204.5 (168, 280)</td><td align="left" rowspan="1" colspan="1">227 (197, 283)</td><td align="left" rowspan="1" colspan="1">92.5 (68, 156)</td></tr><tr><td align="left" rowspan="1" colspan="1">Shannon diversity index (H)</td><td align="left" rowspan="1" colspan="1">4.3</td><td align="left" rowspan="1" colspan="1">4.5</td><td align="left" rowspan="1" colspan="1">3.5</td></tr><tr><td align="left" rowspan="1" colspan="1">Simpson diversity index (1-D)</td><td align="left" rowspan="1" colspan="1">0.978</td><td align="left" rowspan="1" colspan="1">0.984</td><td align="left" rowspan="1" colspan="1">0.96</td></tr><tr><td align="left" rowspan="1" colspan="1">Simpson reciprocal index (1/D)</td><td align="left" rowspan="1" colspan="1">45.45</td><td align="left" rowspan="1" colspan="1">62.5</td><td align="left" rowspan="1" colspan="1">25</td></tr><tr><td align="left" rowspan="1" colspan="1">Evenness index</td><td align="left" rowspan="1" colspan="1">0.558</td><td align="left" rowspan="1" colspan="1">0.556</td><td align="left" rowspan="1" colspan="1">0.637</td></tr><tr><td align="left" rowspan="1" colspan="1">Good's estimator of coverage, %</td><td align="left" rowspan="1" colspan="1">90.73</td><td align="left" rowspan="1" colspan="1">89.35</td><td align="left" rowspan="1" colspan="1">77.2</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0003064-t005" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003064.t005</object-id><label>Table 5</label><caption><title>Mean numbers of predominant bacteria in fecal samples of FMF patients in both stages of the disease and healthy controls detected by FISH.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Probes</td><td colspan="2" align="left" rowspan="1">Controls (n = 7)</td><td colspan="4" align="left" rowspan="1">FMF patients</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td colspan="2" align="left" rowspan="1"/><td colspan="2" align="left" rowspan="1">Remission (n = 6)</td><td colspan="2" align="left" rowspan="1">Attack (n = 7)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Bacterial number/g<xref ref-type="table-fn" rid="nt112">*</xref>
</td><td align="left" rowspan="1" colspan="1">% bacteria<xref ref-type="table-fn" rid="nt113">**</xref>
</td><td align="left" rowspan="1" colspan="1">Bacterial number/g</td><td align="left" rowspan="1" colspan="1">% bacteria</td><td align="left" rowspan="1" colspan="1">Bacterial number/g</td><td align="left" rowspan="1" colspan="1">% bacteria</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Eub338</td><td align="left" rowspan="1" colspan="1">547.2±315.5</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">512.4±250.4</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">356.1±83.5</td><td align="left" rowspan="1" colspan="1">100</td></tr><tr><td align="left" rowspan="1" colspan="1">Bac303</td><td align="left" rowspan="1" colspan="1">151.3±106.0</td><td align="left" rowspan="1" colspan="1">27.3±9.5</td><td align="left" rowspan="1" colspan="1">185.9±136.0</td><td align="left" rowspan="1" colspan="1">39.1±18.7</td><td align="left" rowspan="1" colspan="1">135.4±55.6</td><td align="left" rowspan="1" colspan="1">38.6±15.7</td></tr><tr><td align="left" rowspan="1" colspan="1">Erec482</td><td align="left" rowspan="1" colspan="1">108.5±60.6</td><td align="left" rowspan="1" colspan="1">21.1±5.3</td><td align="left" rowspan="1" colspan="1">101.0±69.8</td><td align="left" rowspan="1" colspan="1">18.0±8.0</td><td align="left" rowspan="1" colspan="1">100.8±37.0</td><td align="left" rowspan="1" colspan="1">28.9±9.6</td></tr><tr><td align="left" rowspan="1" colspan="1">Rrec584</td><td align="left" rowspan="1" colspan="1">45.1±54.2</td><td align="left" rowspan="1" colspan="1">6.5±4.1</td><td align="left" rowspan="1" colspan="1">44.3±35.5</td><td align="left" rowspan="1" colspan="1">8.3±4.5</td><td align="left" rowspan="1" colspan="1">60.3±53.0</td><td align="left" rowspan="1" colspan="1">18.0±14.1</td></tr><tr><td align="left" rowspan="1" colspan="1">Fprau645</td><td align="left" rowspan="1" colspan="1">54.4±19.6</td><td align="left" rowspan="1" colspan="1">11.3±3.2</td><td align="left" rowspan="1" colspan="1">49.8±38.3</td><td align="left" rowspan="1" colspan="1">10.4±6.1</td><td align="left" rowspan="1" colspan="1">53.2±47.2</td><td align="left" rowspan="1" colspan="1">13.8±8.6</td></tr><tr><td align="left" rowspan="1" colspan="1">Rfla729/Rbro730</td><td align="left" rowspan="1" colspan="1">59.0±45.0</td><td align="left" rowspan="1" colspan="1">9.9±3.2</td><td align="left" rowspan="1" colspan="1">28.5±35.7</td><td align="left" rowspan="1" colspan="1">5.3±4.7<xref ref-type="table-fn" rid="nt114">†</xref>
</td><td align="left" rowspan="1" colspan="1">8.5±9.3<xref ref-type="table-fn" rid="nt115">††</xref>
</td><td align="left" rowspan="1" colspan="1">2.6±2.8<xref ref-type="table-fn" rid="nt115">††</xref>
</td></tr><tr><td align="left" rowspan="1" colspan="1">Prop853</td><td align="left" rowspan="1" colspan="1">64.2±51.7</td><td align="left" rowspan="1" colspan="1">10.8±3.2</td><td align="left" rowspan="1" colspan="1">44.3±56.2</td><td align="left" rowspan="1" colspan="1">7.1±7.6</td><td align="left" rowspan="1" colspan="1">12.8±10.4<xref ref-type="table-fn" rid="nt115">††</xref>
</td><td align="left" rowspan="1" colspan="1">3.5±2.5<xref ref-type="table-fn" rid="nt115">††</xref>
</td></tr><tr><td align="left" rowspan="1" colspan="1">Bif164</td><td align="left" rowspan="1" colspan="1">20.0±27.3</td><td align="left" rowspan="1" colspan="1">2.7±3.1</td><td align="left" rowspan="1" colspan="1">27.0±44.6</td><td align="left" rowspan="1" colspan="1">4.6±6.1</td><td align="left" rowspan="1" colspan="1">5.0±5.9</td><td align="left" rowspan="1" colspan="1">1.7±2.3</td></tr><tr><td align="left" rowspan="1" colspan="1">Ato291</td><td align="left" rowspan="1" colspan="1">42.2±24.9</td><td align="left" rowspan="1" colspan="1">7.8±4.1</td><td align="left" rowspan="1" colspan="1">28.9±20.6</td><td align="left" rowspan="1" colspan="1">6.3±3.0</td><td align="left" rowspan="1" colspan="1">16.9±12.7<xref ref-type="table-fn" rid="nt114">†</xref>
</td><td align="left" rowspan="1" colspan="1">4.5±2.9</td></tr><tr><td align="left" rowspan="1" colspan="1">EnterobactD</td><td align="left" rowspan="1" colspan="1">5.4±3.0</td><td align="left" rowspan="1" colspan="1">1.2±0.6</td><td align="left" rowspan="1" colspan="1">4.6±3.6</td><td align="left" rowspan="1" colspan="1">1.2±1.5</td><td align="left" rowspan="1" colspan="1">7.9±12.0</td><td align="left" rowspan="1" colspan="1">1.9±2.9</td></tr><tr><td align="left" rowspan="1" colspan="1">Lab158</td><td align="left" rowspan="1" colspan="1">2.0±3.1</td><td align="left" rowspan="1" colspan="1">0.4±0.4</td><td align="left" rowspan="1" colspan="1">1.9±2.4</td><td align="left" rowspan="1" colspan="1">0.3±0.3</td><td align="left" rowspan="1" colspan="1">1.3±2.0</td><td align="left" rowspan="1" colspan="1">0.4±0.8</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0003064-t006" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003064.t006</object-id><label>Table 6</label><caption><title>Distribution of main bacterial groups in two FMF patients in remission and attack.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Probes</td><td colspan="2" align="left" rowspan="1">Patient FMF 87</td><td align="left" rowspan="1" colspan="1">P-value</td><td colspan="2" align="left" rowspan="1">Patient FMF 13-123</td><td align="left" rowspan="1" colspan="1">P-value</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Remission</td><td align="left" rowspan="1" colspan="1">Attack</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Remission</td><td align="left" rowspan="1" colspan="1">Attack</td><td align="left" rowspan="1" colspan="1"/></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Eub338<xref ref-type="table-fn" rid="nt117">*</xref>
</td><td align="left" rowspan="1" colspan="1">100% (684.08±94.01)</td><td align="left" rowspan="1" colspan="1">100% (318.17±69.16)</td><td align="left" rowspan="1" colspan="1">0.0001</td><td align="left" rowspan="1" colspan="1">100% (656.03±144.49)</td><td align="left" rowspan="1" colspan="1">100% (389.53±86.23)</td><td align="left" rowspan="1" colspan="1">0.0001</td></tr><tr><td align="left" rowspan="1" colspan="1">Bac303</td><td align="left" rowspan="1" colspan="1">65.21±14.05</td><td align="left" rowspan="1" colspan="1">58.23±11.55</td><td align="left" rowspan="1" colspan="1">0.0037</td><td align="left" rowspan="1" colspan="1">26.02±7.11</td><td align="left" rowspan="1" colspan="1">25.01±6.17</td><td align="left" rowspan="1" colspan="1">0.2858</td></tr><tr><td align="left" rowspan="1" colspan="1">Erec482</td><td align="left" rowspan="1" colspan="1">12.32±3.5</td><td align="left" rowspan="1" colspan="1">26.63±5.9</td><td align="left" rowspan="1" colspan="1">0.0001</td><td align="left" rowspan="1" colspan="1">27.91±5.21</td><td align="left" rowspan="1" colspan="1">45.12±14.64</td><td align="left" rowspan="1" colspan="1">0.0001</td></tr><tr><td align="left" rowspan="1" colspan="1">Rrec584</td><td align="left" rowspan="1" colspan="1">4.29±1.02</td><td align="left" rowspan="1" colspan="1">18.04±4.08</td><td align="left" rowspan="1" colspan="1">0.0001</td><td align="left" rowspan="1" colspan="1">13.63±3.03</td><td align="left" rowspan="1" colspan="1">43.91±9.28</td><td align="left" rowspan="1" colspan="1">0.0001</td></tr><tr><td align="left" rowspan="1" colspan="1">Fprau645</td><td align="left" rowspan="1" colspan="1">8.51±2.53</td><td align="left" rowspan="1" colspan="1">12.0±3.24</td><td align="left" rowspan="1" colspan="1">0.0001</td><td align="left" rowspan="1" colspan="1">18.46±3.67</td><td align="left" rowspan="1" colspan="1">16.80±3.69</td><td align="left" rowspan="1" colspan="1">0.0092</td></tr><tr><td align="left" rowspan="1" colspan="1">Rfla729/Rbro730</td><td align="left" rowspan="1" colspan="1">0.57±0.25</td><td align="left" rowspan="1" colspan="1">0.31±0.13</td><td align="left" rowspan="1" colspan="1">0.0012</td><td align="left" rowspan="1" colspan="1">3.90±1.32</td><td align="left" rowspan="1" colspan="1">0.75±0.37</td><td align="left" rowspan="1" colspan="1">0.0001</td></tr><tr><td align="left" rowspan="1" colspan="1">Prop853</td><td align="left" rowspan="1" colspan="1">2.64±0.60</td><td align="left" rowspan="1" colspan="1">1.97±0.55</td><td align="left" rowspan="1" colspan="1">0.0010</td><td align="left" rowspan="1" colspan="1">5.22±1.16</td><td align="left" rowspan="1" colspan="1">0.37±0.22</td><td align="left" rowspan="1" colspan="1">0.0001</td></tr><tr><td align="left" rowspan="1" colspan="1">Bif164</td><td align="left" rowspan="1" colspan="1">1.82±0.42</td><td align="left" rowspan="1" colspan="1">0.59±0.22</td><td align="left" rowspan="1" colspan="1">0.0001</td><td align="left" rowspan="1" colspan="1">3.31±0.73</td><td align="left" rowspan="1" colspan="1">0.16±0.11</td><td align="left" rowspan="1" colspan="1">0.0001</td></tr><tr><td align="left" rowspan="1" colspan="1">Ato291</td><td align="left" rowspan="1" colspan="1">2.30±0.65</td><td align="left" rowspan="1" colspan="1">3.91±1.23</td><td align="left" rowspan="1" colspan="1">0.0001</td><td align="left" rowspan="1" colspan="1">4.06±0.86</td><td align="left" rowspan="1" colspan="1">3.59±1.33</td><td align="left" rowspan="1" colspan="1">0.0262</td></tr><tr><td align="left" rowspan="1" colspan="1">EnterobactD</td><td align="left" rowspan="1" colspan="1">0.50±0.1</td><td align="left" rowspan="1" colspan="1">0.32±0.12</td><td align="left" rowspan="1" colspan="1">0.0006</td><td align="left" rowspan="1" colspan="1">0.65±0.18</td><td align="left" rowspan="1" colspan="1">1.01±0.72</td><td align="left" rowspan="1" colspan="1">0.0013</td></tr><tr><td align="left" rowspan="1" colspan="1">Lab158</td><td align="left" rowspan="1" colspan="1">0.81±0.23</td><td align="left" rowspan="1" colspan="1">0.67±0.24</td><td align="left" rowspan="1" colspan="1">0.1617</td><td align="left" rowspan="1" colspan="1">0.64±0.16</td><td align="left" rowspan="1" colspan="1">N/D<xref ref-type="table-fn" rid="nt118">†</xref>
</td><td align="left" rowspan="1" colspan="1">–</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0003064-t007" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003064.t007</object-id><label>Table 7</label><caption><title>Distribution of main bacterial groups in respect to colchicine treatment.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Probes</td><td colspan="4" align="left" rowspan="1">FMF remission</td><td colspan="4" align="left" rowspan="1">FMF attack</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td colspan="2" align="left" rowspan="1">Colchicine treatment (n = 3)</td><td colspan="2" align="left" rowspan="1">No colchicine treatment (n = 3)</td><td colspan="2" align="left" rowspan="1">Colchicine treatment (n = 3)</td><td colspan="2" align="left" rowspan="1">No colchicine treatment (n = 4)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Bacterial number/g<xref ref-type="table-fn" rid="nt120">*</xref>
</td><td align="left" rowspan="1" colspan="1">% Bacteria<xref ref-type="table-fn" rid="nt121">**</xref>
</td><td align="left" rowspan="1" colspan="1">Bacterial number/g</td><td align="left" rowspan="1" colspan="1">% bacteria</td><td align="left" rowspan="1" colspan="1">Bacterial number/g</td><td align="left" rowspan="1" colspan="1">% bacteria</td><td align="left" rowspan="1" colspan="1">Bacterial number/g</td><td align="left" rowspan="1" colspan="1">% bacteria</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Eub338</td><td align="left" rowspan="1" colspan="1">668.7±21.6</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">356.0±288.0</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">385.9±22.8</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">333.7±109.7</td><td align="left" rowspan="1" colspan="1">100</td></tr><tr><td align="left" rowspan="1" colspan="1">Bac303</td><td align="left" rowspan="1" colspan="1">166.1±41.1</td><td align="left" rowspan="1" colspan="1">25±6.8</td><td align="left" rowspan="1" colspan="1">205.7±208.2</td><td align="left" rowspan="1" colspan="1">53.3±15.0</td><td align="left" rowspan="1" colspan="1">146.5±68.9</td><td align="left" rowspan="1" colspan="1">38.7±20.5</td><td align="left" rowspan="1" colspan="1">127.1±53.0</td><td align="left" rowspan="1" colspan="1">38.5±14.5</td></tr><tr><td align="left" rowspan="1" colspan="1">Erec482</td><td align="left" rowspan="1" colspan="1">156.0±42.0</td><td align="left" rowspan="1" colspan="1">23.5±6.9</td><td align="left" rowspan="1" colspan="1">46.0±36.5</td><td align="left" rowspan="1" colspan="1">12.6±4.7</td><td align="left" rowspan="1" colspan="1">110.9±56.6</td><td align="left" rowspan="1" colspan="1">28.8±14.5</td><td align="left" rowspan="1" colspan="1">93.2±20.8</td><td align="left" rowspan="1" colspan="1">28.9±6.8</td></tr><tr><td align="left" rowspan="1" colspan="1">Rrec584</td><td align="left" rowspan="1" colspan="1">69.0±34.6</td><td align="left" rowspan="1" colspan="1">10.4±5.4</td><td align="left" rowspan="1" colspan="1">19.6±11.3</td><td align="left" rowspan="1" colspan="1">6.2±2.7</td><td align="left" rowspan="1" colspan="1">78.1±83.0</td><td align="left" rowspan="1" colspan="1">20.4±21.2</td><td align="left" rowspan="1" colspan="1">46.9±21.6</td><td align="left" rowspan="1" colspan="1">16.2±9.5</td></tr><tr><td align="left" rowspan="1" colspan="1">Fprau645</td><td align="left" rowspan="1" colspan="1">62.3±51.2</td><td align="left" rowspan="1" colspan="1">9.4±7.9</td><td align="left" rowspan="1" colspan="1">37.3±24.0</td><td align="left" rowspan="1" colspan="1">11.3±5.4</td><td align="left" rowspan="1" colspan="1">42.5±20.1</td><td align="left" rowspan="1" colspan="1">11.0±5.2</td><td align="left" rowspan="1" colspan="1">61.2±63.1</td><td align="left" rowspan="1" colspan="1">15.8±10.8</td></tr><tr><td align="left" rowspan="1" colspan="1">Rfla729/Rbro730</td><td align="left" rowspan="1" colspan="1">50.4±41.8</td><td align="left" rowspan="1" colspan="1">7.4±5.9</td><td align="left" rowspan="1" colspan="1">6.6±2.3</td><td align="left" rowspan="1" colspan="1">3.2±2.5</td><td align="left" rowspan="1" colspan="1">11.0±12.7</td><td align="left" rowspan="1" colspan="1">2.8±3.1</td><td align="left" rowspan="1" colspan="1">6.6±7.3</td><td align="left" rowspan="1" colspan="1">2.4±3.0</td></tr><tr><td align="left" rowspan="1" colspan="1">Prop853</td><td align="left" rowspan="1" colspan="1">79.3±64.5</td><td align="left" rowspan="1" colspan="1">11.7±9.1</td><td align="left" rowspan="1" colspan="1">9.3±7.8</td><td align="left" rowspan="1" colspan="1">2.5±0.2</td><td align="left" rowspan="1" colspan="1">14±15.8</td><td align="left" rowspan="1" colspan="1">3.5±3.8</td><td align="left" rowspan="1" colspan="1">11.9±7.1</td><td align="left" rowspan="1" colspan="1">3.5±1.7</td></tr><tr><td align="left" rowspan="1" colspan="1">Bif164</td><td align="left" rowspan="1" colspan="1">46.4±61.8</td><td align="left" rowspan="1" colspan="1">6.8±8.9</td><td align="left" rowspan="1" colspan="1">7.5±4.4</td><td align="left" rowspan="1" colspan="1">2.5±1.2</td><td align="left" rowspan="1" colspan="1">3.4±5.4</td><td align="left" rowspan="1" colspan="1">0.9±1.5</td><td align="left" rowspan="1" colspan="1">6.2±6.8</td><td align="left" rowspan="1" colspan="1">2.3±2.8</td></tr><tr><td align="left" rowspan="1" colspan="1">Ato291</td><td align="left" rowspan="1" colspan="1">41.8±23.4</td><td align="left" rowspan="1" colspan="1">6.2±3.2</td><td align="left" rowspan="1" colspan="1">16.0±4.4</td><td align="left" rowspan="1" colspan="1">6.3±3.5</td><td align="left" rowspan="1" colspan="1">24.7±16.2</td><td align="left" rowspan="1" colspan="1">6.3±3.8</td><td align="left" rowspan="1" colspan="1">11.0±6.5</td><td align="left" rowspan="1" colspan="1">3.1±1.1</td></tr><tr><td align="left" rowspan="1" colspan="1">EnterobactD</td><td align="left" rowspan="1" colspan="1">4.5±3.2</td><td align="left" rowspan="1" colspan="1">0.7±0.5</td><td align="left" rowspan="1" colspan="1">4.8±4.7</td><td align="left" rowspan="1" colspan="1">1.8±2.0</td><td align="left" rowspan="1" colspan="1">13.2±17.3</td><td align="left" rowspan="1" colspan="1">3.3±4.2</td><td align="left" rowspan="1" colspan="1">4.0±6.3</td><td align="left" rowspan="1" colspan="1">0.9±1.2</td></tr><tr><td align="left" rowspan="1" colspan="1">Lab158</td><td align="left" rowspan="1" colspan="1">1.9±1.9</td><td align="left" rowspan="1" colspan="1">0.3±0.3</td><td align="left" rowspan="1" colspan="1">1.8±3.2</td><td align="left" rowspan="1" colspan="1">0.3±0.5</td><td align="left" rowspan="1" colspan="1">0.2±0.2</td><td align="left" rowspan="1" colspan="1">0.1±0.1</td><td align="left" rowspan="1" colspan="1">2.0±2.5</td><td align="left" rowspan="1" colspan="1">0.7±1.0</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0003064-t008" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003064.t008</object-id><label>Table 8</label><caption><title>Classification table of correct estimates based on FISH variables.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Studied groups</td><td align="left" rowspan="1" colspan="1">The percentage of correct classification</td><td align="left" rowspan="1" colspan="1">Healthy controls</td><td colspan="2" align="left" rowspan="1">Disease phases</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Remission</td><td align="left" rowspan="1" colspan="1">Attack</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Healthy controls</td><td align="left" rowspan="1" colspan="1">100.0% (35.0%)<xref ref-type="table-fn" rid="nt123">*</xref>
</td><td align="left" rowspan="1" colspan="1">7</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">FMF remission</td><td align="left" rowspan="1" colspan="1">83.3% (30.0%)</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">1</td></tr><tr><td align="left" rowspan="1" colspan="1">FMF attack</td><td align="left" rowspan="1" colspan="1">100.0% (35.0%)</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">7</td></tr><tr><td align="left" rowspan="1" colspan="1">Total</td><td align="left" rowspan="1" colspan="1">95.0% (33.3%)</td><td align="left" rowspan="1" colspan="1">7</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">8</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0003064-t009" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003064.t009</object-id><label>Table 9</label><caption><title>Oligonucleotide probes and hybridization conditions used in this study.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Probe</td><td align="left" rowspan="1" colspan="1">Sequence (5′-3′)</td><td align="left" rowspan="1" colspan="1">Targeted bacteria</td><td colspan="3" align="left" rowspan="1">Hybridization conditions<xref ref-type="table-fn" rid="nt124">*</xref>
</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">T</td><td align="left" rowspan="1" colspan="1">F</td><td align="left" rowspan="1" colspan="1">L</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Eub338</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">GCTGCCTCCCGTAGGAGT</named-content>
</td><td align="left" rowspan="1" colspan="1">
<italic>Bacteria</italic>
</td><td align="left" rowspan="1" colspan="1">50</td><td align="left" rowspan="1" colspan="1">-</td><td align="left" rowspan="1" colspan="1">-</td></tr><tr><td align="left" rowspan="1" colspan="1">Bac303</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">CCAATGTGGGGGACCTT</named-content>
</td><td align="left" rowspan="1" colspan="1">
<italic>Bacteroides</italic> group</td><td align="left" rowspan="1" colspan="1">47</td><td align="left" rowspan="1" colspan="1">-</td><td align="left" rowspan="1" colspan="1">-</td></tr><tr><td align="left" rowspan="1" colspan="1">Erec482</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">GCTTCTTAGTCAGGTACCG</named-content>
</td><td align="left" rowspan="1" colspan="1">Clostridial cluster XIVa</td><td align="left" rowspan="1" colspan="1">47</td><td align="left" rowspan="1" colspan="1">-</td><td align="left" rowspan="1" colspan="1">-</td></tr><tr><td align="left" rowspan="1" colspan="1">Rrec584</td><td align="left" rowspan="1" colspan="1">TCAGACTTGCCG(C/T)ACCGC</td><td align="left" rowspan="1" colspan="1">
<italic>Roseburia-E. rectale</italic> group</td><td align="left" rowspan="1" colspan="1">50</td><td align="left" rowspan="1" colspan="1">-</td><td align="left" rowspan="1" colspan="1">-</td></tr><tr><td align="left" rowspan="1" colspan="1">Fprau645</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">CCTCTGCACTACTCAAGAAAAAC</named-content>
</td><td align="left" rowspan="1" colspan="1">
<italic>F. prausnitzii</italic> group</td><td align="left" rowspan="1" colspan="1">50</td><td align="left" rowspan="1" colspan="1">-</td><td align="left" rowspan="1" colspan="1">-</td></tr><tr><td align="left" rowspan="1" colspan="1">Rfla729</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">AAAGCCCAGTAAGCCGCC</named-content>
</td><td align="left" rowspan="1" colspan="1">
<italic>R. flavefaciens</italic> subcluster</td><td align="left" rowspan="1" colspan="1">50</td><td align="left" rowspan="1" colspan="1">20</td><td align="left" rowspan="1" colspan="1">10</td></tr><tr><td align="left" rowspan="1" colspan="1">Rbro730</td><td align="left" rowspan="1" colspan="1">TAAAGCCCAG(C/T)AGGCCGC</td><td align="left" rowspan="1" colspan="1">
<italic>R. bromii</italic> sublcuster</td><td align="left" rowspan="1" colspan="1">50</td><td align="left" rowspan="1" colspan="1">20</td><td align="left" rowspan="1" colspan="1">10</td></tr><tr><td align="left" rowspan="1" colspan="1">Prop853</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">ATTGCGTTAACTCCGGCAC</named-content>
</td><td align="left" rowspan="1" colspan="1">Clostridial cluster IX</td><td align="left" rowspan="1" colspan="1">50</td><td align="left" rowspan="1" colspan="1">-</td><td align="left" rowspan="1" colspan="1">-</td></tr><tr><td align="left" rowspan="1" colspan="1">Bif164</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">CATCCGGCATTACCACCC</named-content>
</td><td align="left" rowspan="1" colspan="1">
<italic>Bifidobacterium</italic> genus</td><td align="left" rowspan="1" colspan="1">50</td><td align="left" rowspan="1" colspan="1">-</td><td align="left" rowspan="1" colspan="1">-</td></tr><tr><td align="left" rowspan="1" colspan="1">Ato291</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">GGTCGGTCTCTCAACCC</named-content>
</td><td align="left" rowspan="1" colspan="1">
<italic>Atopobium</italic> group</td><td align="left" rowspan="1" colspan="1">50</td><td align="left" rowspan="1" colspan="1">-</td><td align="left" rowspan="1" colspan="1">-</td></tr><tr><td align="left" rowspan="1" colspan="1">EnterobactD</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">TGCTCTCGCGAGGTCGCTTCTCTT</named-content>
</td><td align="left" rowspan="1" colspan="1">
<italic>Enterobacteriaceae</italic> genus</td><td align="left" rowspan="1" colspan="1">50</td><td align="left" rowspan="1" colspan="1">-</td><td align="left" rowspan="1" colspan="1">-</td></tr><tr><td align="left" rowspan="1" colspan="1">Lab158</td><td align="left" rowspan="1" colspan="1">GGTATTAGCA(C/T)CTGTTTCCA</td><td align="left" rowspan="1" colspan="1">
<italic>Lactobacillus-Enterococcus</italic>
</td><td align="left" rowspan="1" colspan="1">45</td><td align="left" rowspan="1" colspan="1">20</td><td align="left" rowspan="1" colspan="1">10</td></tr></tbody></table></alternatives></table-wrap>
<table-wrap-foot><fn id="nt101"><p>Control subjects denoted with prefix C; FMF patients denoted with prefix FMF.</p></fn><fn id="nt102"><p>Abbreviations: S, 16S library; F, FISH analysis; SR, 16S library in FMF remission; SA, 16S library in FMF attack; FR, FISH analysis in FMF remission; FA, FISH analysis in FMF attack; HT, heterozygote; HM, homozygote.</p></fn><fn id="nt103"><label>*</label><p>N/D, not detected.</p></fn><fn id="nt104"><label>**</label><p>N/A, not analyzed.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt105"><label>*</label><p>OTUs were defined at 99% similarity level.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt106"><label>*</label><p>p&lt;0.01, comparison of remission vs. control libraries.</p></fn><fn id="nt107"><label>**</label><p>p&lt;0.01, comparison of attack vs. control libraries.</p></fn><fn id="nt108"><label>***</label><p>p&lt;0.01, comparison of remission vs. attack libraries.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt109"><p>Calculations were made based on OTU definition at 99% sequence identity.</p></fn><fn id="nt110"><label>*</label><p>The values in parentheses are the 95% confidence intervals, which assess the precision of the richness estimates.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt111"><p>Mean count±SD.</p></fn><fn id="nt112"><label>*</label><p>Counts 10<sup>8</sup> cells per gram of fecal specimen.</p></fn><fn id="nt113"><label>**</label><p>The percentage relative to the total number of cells determined by Eub338.</p></fn><fn id="nt114"><label>†</label><p>p&lt;0.05 as compared to healthy controls.</p></fn><fn id="nt115"><label>††</label><p>p&lt;0.01 as compared to healthy controls.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt116"><p>Mean count of the percentage relative to the total number of cells determined by Eub338 ± SD.</p></fn><fn id="nt117"><label>*</label><p>The values in parentheses are the absolute numbers of cells determined by Eub338 ± SD (counts 10<sup>8</sup> cells per gram of fecal specimen).</p></fn><fn id="nt118"><label>†</label><p>N/D, not detected.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt119"><p>Mean count±SD.</p></fn><fn id="nt120"><label>*</label><p>Counts 10<sup>8</sup> cells per gram of fecal specimen.</p></fn><fn id="nt121"><label>**</label><p>The percentage relative to the total number of cells determined by Eub338.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt122"><p>Statistical parameters of the discriminant model: Wilks's lambda coefficient – 0.045. Significance of the model: F (20,16) = 2.95; p&lt;0.016.</p></fn><fn id="nt123"><label>*</label><p>
<italic>apriori</italic> probability of correct classification.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt124"><label>*</label><p>T, Hybridization temperature (°C); F, Formamide concentration (%); L, Lysozyme treatment (min).</p></fn></table-wrap-foot> <fn-group><fn fn-type="COI-statement"><p><bold>Competing Interests: </bold>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p><bold>Funding: </bold>This work was funded by the ISTC (project A-1055). RIA is supported by the RERAD.</p></fn></fn-group>
<graphic id="pone-0003064-t001-1" xlink:href="pone.0003064.t001"/> <graphic id="pone-0003064-t002-2" xlink:href="pone.0003064.t002"/> <graphic id="pone-0003064-t003-3" xlink:href="pone.0003064.t003"/> <graphic id="pone-0003064-t004-4" xlink:href="pone.0003064.t004"/> <graphic xlink:href="pone.0003064.g001"/> <graphic xlink:href="pone.0003064.g002"/> <graphic id="pone-0003064-t005-5" xlink:href="pone.0003064.t005"/> <graphic id="pone-0003064-t006-6" xlink:href="pone.0003064.t006"/> <graphic xlink:href="pone.0003064.g003"/> <graphic id="pone-0003064-t007-7" xlink:href="pone.0003064.t007"/> <graphic xlink:href="pone.0003064.g004"/> <graphic id="pone-0003064-t008-8" xlink:href="pone.0003064.t008"/> <graphic id="pone-0003064-t009-9" xlink:href="pone.0003064.t009"/>
[{"label": ["2"], "element-citation": ["\n"], "surname": ["Finegold", "Sutter", "Mathisen", "Hentges"], "given-names": ["SM", "VL", "GE", "DJ"], "year": ["1983"], "article-title": ["Human Intestinal Microflora in Health and Disease."], "publisher-loc": ["New York"], "publisher-name": ["Academic Press"], "fpage": ["3"], "lpage": ["31"]}, {"label": ["20"], "element-citation": ["\n"], "surname": ["Csordas"], "given-names": ["A"], "year": ["1996"], "article-title": ["Butyrate, aspirin and colorectal cancer."], "source": ["Eur J Cancer Prev"], "volume": ["4"], "fpage": ["221"], "lpage": ["231"]}, {"label": ["34"], "element-citation": ["\n"], "surname": ["Adams", "Heazlewood", "Gilshenan", "O'brien", "McGuckin"], "given-names": ["RJ", "SP", "KS", "M", "MA"], "year": ["2007"], "article-title": ["IgG antibodies against common gut bacteria are more diagnostic for Crohn's disease than IgG against mannan or flagellin."], "source": ["Am J Gastroenterol"], "comment": ["Oct 6 [Epub ahead of print]"]}, {"label": ["38"], "element-citation": ["\n"], "surname": ["Ting", "Kastner", "Hoffman"], "given-names": ["JP", "DL", "HM"], "year": ["2006"], "article-title": ["CATERPILLERs, pyrin and hereditary immunological disorders."], "source": ["Nat Rev Immunology"], "volume": ["3"], "fpage": ["183"], "lpage": ["195"]}, {"label": ["43"], "element-citation": ["\n"], "surname": ["Zoetendal", "Akkermans", "Akkermans van-Vliet", "de Visser", "de Vos"], "given-names": ["EG", "ADL", "WM", "JA", "WM"], "year": ["2001"], "article-title": ["The host genotype affects the bacterial community in the human gastrointestinal tract."], "source": ["Microbial Ecol Health Dis"], "volume": ["13"], "fpage": ["129"], "lpage": ["134"]}, {"label": ["47"], "element-citation": ["\n"], "surname": ["Stoffman", "Magal", "Shohat", "Lotan", "Koman"], "given-names": ["N", "N", "T", "R", "S"], "year": ["2000"], "article-title": ["Higher than expected carrier rates for familial Mediterranean fever in various Jewish ethnic groups."], "source": ["Eur J Hum Gen"], "volume": ["8"], "fpage": ["307"], "lpage": ["310"]}, {"label": ["53"], "element-citation": ["\n"], "surname": ["Tunca", "Kirkali", "Soyturk", "Akar", "Pepys"], "given-names": ["M", "G", "M", "S", "MB"], "year": ["1999"], "article-title": ["Acute phase response and evolution of familial Mediterranean fever."], "source": ["Lancet"], "volume": ["353"], "fpage": ["1415"]}, {"label": ["77"], "element-citation": ["\n"], "surname": ["Thompson", "Gibson", "Plewniak", "Jeanmougin", "Higgins"], "given-names": ["JD", "TJ", "F", "F", "DG"], "year": ["1997"], "article-title": ["The ClustalX windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools."], "source": ["Nucleic Acids Research"], "volume": ["24"], "fpage": ["4876"], "lpage": ["4882"]}, {"label": ["81"], "element-citation": ["\n"], "surname": ["Good"], "given-names": ["IJ"], "year": ["1953"], "article-title": ["The population frequencies of species and the estimation of population parameters."], "source": ["Biometrika"], "volume": ["40"], "fpage": ["237"], "lpage": ["264"]}, {"label": ["92"], "element-citation": ["\n"], "surname": ["Harmsen", "Elfferich", "Schut", "Welling"], "given-names": ["HJ", "P", "F", "GW"], "year": ["1999"], "article-title": ["A 16S rRNA-targeted probe for detection of lactobacilli and enterococci in fecal samples by fluorescent "], "italic": ["in situ"], "source": ["Microb Ecol Health Dis"], "volume": ["11"], "fpage": ["3"], "lpage": ["12"]}]
PLoS One. 2008 Aug 26; 3(8):e3064