High taxonomic level fingerprint of the human intestinal microbiota by Ligase Detection Reaction - Universal Array approach
© Candela et al; licensee BioMed Central Ltd. 2010
Received: 16 September 2009
Accepted: 19 April 2010
Published: 19 April 2010
Affecting the core functional microbiome, peculiar high level taxonomic unbalances of the human intestinal microbiota have been recently associated with specific diseases, such as obesity, inflammatory bowel diseases, and intestinal inflammation.
In order to specifically monitor microbiota unbalances that impact human physiology, here we develop and validate an original DNA-microarray (HTF-Microbi.Array) for the high taxonomic level fingerprint of the human intestinal microbiota. Based on the Ligase Detection Reaction-Universal Array (LDR-UA) approach, the HTF-Microbi.Array enables specific detection and approximate relative quantification of 16S rRNAs from 30 phylogenetically related groups of the human intestinal microbiota. The HTF-Microbi.Array was used in a pilot study of the faecal microbiota of eight young adults. Cluster analysis revealed the good reproducibility of the high level taxonomic microbiota fingerprint obtained for each of the subject.
The HTF-Microbi.Array is a fast and sensitive tool for the high taxonomic level fingerprint of the human intestinal microbiota in terms of presence/absence of the principal groups. Moreover, analysis of the relative fluorescence intensity for each probe pair of our LDR-UA platform can provide estimation of the relative abundance of the microbial target groups within each samples. Focusing the phylogenetic resolution at division, order and cluster levels, the HTF-Microbi.Array is blind with respect to the inter-individual variability at the species level.
Human beings have been recently reconsidered as superorganisms in co-evolution with an immense microbial community living in the gastrointestinal tract (GIT), the human intestinal microbiota [1, 2]. Providing important metabolic functions that we have not evolved by our own , the intestinal microbiota has a fundamental role for the human health and well being [4, 5]. Several of our physiological features, such as nutrient processing, maturation of the immune system, pathogen resistance, and development of the intestinal architecture, strictly depend on the mutualistic symbiotic relationship with the intestinal microbiota . On the basis of its global impact on human physiology, the intestinal microbiota has been considered an essential organ of the human body .
The composition of the adult intestinal microbiota has been determined in three large scale 16S rRNA sequences surveys [7–11]. The phylogenetic analysis of a total of 45,000 bacterial 16S rRNA data from 139 adults revealed that, at the phylum level, only a small fraction of the known bacterial diversity is represented in our GIT. The vast majority of bacteria in the human intestinal microbiota (>99%) belongs to six bacterial phyla: Firmicutes, Bacteroidetes, Actinobacteria, Proteobacteria, Fusobacteria and Verrucomicrobia. The two dominant divisions are Firmicutes and Bacteroidetes, which represent together up to 90% of the total microbiota, with a relative abundance of 65% and 25%, respectively. Actinobacteria, Proteobacteria, Verrucomicrobia and Fusobacteria are the subdominants phyla with a relative abundance up to 5, 8, 2 and 1%, respectively. On the contrary, at lower taxonomic levels, we assist to a real explosion of the bacterial diversity in the human GIT. At least 1,800 genera [≥ 90% of sequence identity (ID)] and 16,000 phylotypes at the species level (≥ 97% ID) have been identified until now, predicting even a greater diversity at the species level . Since 70% of these phylotypes are subject-specific, and no phylotype is present at more than 0.5% abundance in all subjects , the intestinal microbiota of each individual has been shown to consist in a subject specific complement of hundreds of genera and thousands of species. However, the large degree of functional redundancy between species and genera allowed identifying a core microbiome at the gene level which is shared between all individuals . Coding for genes involved in important metabolic functions, this core functional microbiome is fundamental to support the mutualistic symbiotic relationship with the human host.
Recently, 16S rRNA sequences studies have been carried out with the attempt to describe disease-associated unbalances of the human intestinal microbiota. Even though species variability was associated with inter-individual variability, phylum-level changes of the intestinal microbiota were associated with specific diseases. In particular, obesity was characterized by a higher proportion of Firmicutes and Actinobacteria with respect to Bacteroidetes and an overall reduced bacterial diversity [12, 13]. Differently, inflammatory bowel diseases (IBD) were characterized by a marked reduction of bacterial diversity in the Clostridium cluster IV and XIVa belonging to Firmicutes, a decline in Bacteroidetes biodiversity, and a correspondent increase in Proteobacteria and Bacillus[14, 15]. Analogously, intestinal inflammation has been generally related with a marked increase in Enterobacteriaceae and a correspondent decrease in members of the resident colonic bacteria [16, 17]. In the light of these findings, it has been recently hypothesized that these high level taxonomic unbalances of the human intestinal microbiota can cause deviations from the core functional microbiome with a final impact on the host physiological state [12, 18, 19].
Since more than 75% of the phylotypes detected in the human GIT does not correspond to cultured species , phylogenetic DNA-microarrays have been recognized as a valuable tool for a high-throughput, quantitative and systematic analysis of the human intestinal microbiota . Recently, three different small ribosomal subunit RNA (SSU rRNA) based high-density phylogenetic microarrays for studying the human microbiota have been developed [22–24]. Targeting thousands bacterial phylo-types, these DNA-microarrays have been successfully applied in studies for the deep phylogenetic characterization of the human intestinal microbiota.
In order to specifically monitor the microbiota unbalances that impact on human physiology independently of the inter-individual variability, here we developed an original DNA-microarray for the high taxonomic level fingerprint of the human intestinal microbiota, called HTF-Microbi.Array (High Taxonomic Fingerprint Microbiota Array). The relatively low number of targets allowed implementing the Ligase Detection Reaction (LDR) technology [25, 26] for the development of the HTF-Microbi.Array. This enzymatic in vitro reaction, based on the discriminative properties of the DNA ligation enzyme, requires the design of a pair of two adjacent oligonucleotides specific for each target sequence: a probe specific for the variation (called "Discriminating Probe", or DS) which carries a 5'-fluorescent label, and a second probe, named "Common Probe" (or CP), starting one base 3'-downstream of the DS that carries a 5'-phosphate group and a unique sequence named cZipCode at its 3'-end. The oligonucleotide probe pairs and a thermostable DNA ligase are used in a LDR reaction with previously PCR-amplified DNA fragments. This reaction is cycled to increase product yield. The LDR products, obtained only in presence of a perfectly matching template by action of the DNA ligase, are addressed to a precise location onto a Universal Array (UA), where a set of artificial sequences, called Zip-codes are arranged. These products carry both the fluorescent label and a unique cZipCode sequence and can be detected by laser scanning and identified according to their location within the array. The LDR approach is a highly specific and sensitive assay for detecting single nucleotide variations; thus, differences of a single base along the 16S rRNA gene can be employed to distinguish among different microbial lineages. The HTF-Microbi.Array was successfully tested in a pilot study for the characterization of the faecal microbiota of eight healthy young adults.
Target selection and probe design
Probe set of the HTF-Microbi.Array.
H.G. AB %
Species (et rel)
Species (et rel)
Species (et rel)
Species (et rel)
Species (et rel)
Species (et rel)
Species (et rel)
Species (et rel)
Species (et rel)
Clostridium I and II
Cl I and II
Species (et rel)
Cl I and II
Species (et rel)
Species (et rel)
Species (et rel)
Species (et rel)
Specificity and coverage of each candidate probe was assessed by using the tool Probe Match of the RDP database. The probe pairs selected for the HTF-Microbi.Array were required to perfectly match the sequences of the positive set and to possess at least a mismatch at the 3' end of the discriminating probe respect to the entire negative set. The designed probes pairs had an average melting temperature (Tm) of 67.8 ± 0.9°C (n = 60) and an average length of 35.6 ± 4.9 nucleotides. Sixteen out of the 30 probe pairs were characterized by having no degenerated bases, whereas only one probe pair (i.e. the one for Clostridium cluster I and II) had 4 and 3 ambiguous bases on DS and CP, respectively (Additional file 2).
Validation of the HTF-Microbi.Array
LDR probe pair specificity
B. fragilis ATCC25285
B. thetaiotaomicrom ATCC29143
L. gasseri DSM20243
P. melaninogenica ATCC25845
B. subtilis DSM704
E. coli ATCC11105
P. mirabilis DSM4479
B. bifidum DSM20456
L. casei DSM20011
Lactobacillaceae, L. casei
Y. enterocolitica (faecal isolate)
Yersinia enterocolitica, Enterobacteriaceae
B. cereus DSM31
B. adolescentis ATCC15703
L. ramnosus DSM20021
Lactobacillaceae, L. casei
L. delbrueckii DSM20074
L. pentosus DSM20314
L. acidophilus DSM20079
L. reuteri DSM20016
L. plantarum DSM21074
Lactobacillaceae, L. plantarum
C. difficile ATCCBAA1382
Clostridium XI, Clostridium difficile
C. jejuni ATCC33292
V. parvula ATCC10790
Veillonella, Clostridium IX
B. breve DSM20091
B. longum ATCC15707
Bifidobacteriaceae, B. longum
R. productus ATCC 23340
L. salivarius SV2
Lactobacillaceae, L. salivarius
E. faecalis ATCC700802
C. leptum DSM73
Ruminocuccos bromii Clostridium IV
R. albus DSM20455
Ruminocuccos albus Clostridium IV
Evaluation of the LDR sensitivity and relative abundance detection level
Characterization of the faecal microbiota of eight healthy young adults
In these last years, 16S rRNA microarrays emerged as a sensitive and efficient way to screen complex bacterial communities. Here we describe and validate the HTF-Microbi.Array, a new phylogenetic DNA microarray designed for the high taxonomic level fingerprint of the human intestinal microbial community. The HTF-Microbi.Array is based on the LDR-UA approach, which is a fast and sensitive tool for the characterization of complex microbial communities with high sensitivity and specificity [25, 26]. The use of this molecular technique allows overcoming the major limitations of DNA microarrays whose discriminative power is based on hybridization. In fact, a) optimization of the hybridization conditions for each probe set is not required; b) problems due to the secondary structures of the target DNA are minimized, c) steric hindrances of differentially sized nucleic acid hybrids formed on the array after the hybridization are decreased . The final probe set of the HTF-Microbi.Array allows a high taxonomic level fingerprint of the human intestinal microbiota, with a good coverage of the major and minor components, as well as some of the most important pathogens and opportunistic bacteria . The LDR probes were designed by choosing DS oligonucleotides whose 3'end allowed the perfect discrimination of the target species from the non-target ones on the basis of our 16S rRNA sequence database. Definition of accurate and specific negative sets of gut microbiota sequences by ORMA tool  allowed the selection of maximally discriminative probe pairs. Probe specificity was confirmed on the entire known 16S rRNA gene sequences environment by the RDP Probe Match tool. This requirement is fundamental, since the primer set used for the PCR amplification was the "universal" 16S rRNA primer set designed by Edwards and co-workers .
The HTF-Microbi.Array recognized without ambiguity the 16S rRNA amplicons obtained from 28 members of the intestinal microbiota belonging to Bacteroides/Prevotella, Clostridium clusters IV, IX, XIVa, XI, I and II, Bifidobacteriaceae, Lactobacillaceae, Bacillus, Enterococcus, Enterobacteriaceae and Campylobacter, demonstrating the specificity of all the probe pairs. The sensitivity of the HTF-Microbi.Array was evaluated by using different concentrations of an artificial mix of 16S rRNA amplicons obtained from 6 microorganisms members of the human intestinal microbiota. To compensate the eventual drop in the signal due to a very low target concentrations, lower than 0.7 fmol (i.e. a percentage lower than 1.5% of the commonly used quantity of 50 fmol), a slightly relaxed criteria for significance of the t-test to α = 0.05 was chosen. All PCR products were specifically recognized in a concentration range from 75 to 0.7 fmol, showing high array sensitivity. The efficiency of the HTF-Microbi.Array in the detection of a particular target in a complex DNA environment was also determined. According to our data, the array is able to detect a specific DNA target down to 0.02% of the total 16S rRNA, which is comparable to the values obtained by Rajilic-Stojanovic et al. and Palmer et al.. Thus the HTF-Microbi.Array shows the potentiality to sense low abundant species of the gastrointestinal microbiota, enabling the detection of the 16S rRNA of a peculiar target group present at a fractional abundance <0.1% in an artificial mixture.
The HTF-Microbi.Array was used in a pilot study to characterize the faecal microbiota of eight young adults. Faecal microbiota was chosen as DNA source since sample collection is not invasive, samples contain large amount of microbes, and, most important, it is representative of interpersonal differences in distal gut microbial ecology . In order to have a good representation of the less abundant species of the intestinal microbial community, LDR reactions were performed starting from 50 fmol of PCR product. Cluster analysis of the presence-absence probes profiles enabled the identification of a reproducible high taxonomic level microbiota fingerprint for each subject. As expected, the intestinal microbial community of the voluntaries in the study resembled the typical fingerprint of healthy adults . According to our data, the faecal microbiota of the enrolled subjects was dominated by major mutualistic symbionts. In fact, members of Bacteroidetes, Clostridium clusters IV, IX and XIVa were all represented in 100% of the subjects. On the other hand, minor mutualistic symbionts, such as Lactobacillaceae, B. subtilis et re., Fusobacterium and Cyanobacteria, were detected in 55, 37, 50, and 63% of the subjects, respectively. Opportunistic pathogens, such as E. faecalis et rel., members of the Clostridium cluster I and II and Enterobacteriaceae, were represented only in 43, 25 and 12% of the subjects, respectively. Most importantly, enteropathogens such as, C. difficile, C. perfringens, E. faecium et rel., B. cereus et rel., and Campylobacter were never detected. A discrepancy between our data and the literature is the relatively low prevalence of the health promoting Bifidobacteriaceae in our samples (only 13% of samples). However, the low prevalence of bifidobacteria is a typical bias for several phylogenetic DNA microarrays [22, 23]. Probably this is due to the intrinsic low efficiency of amplification of the bifidobacterial genome with universal primer sets for the 16S rRNA gene . Surprisingly, a high prevalence was obtained for the minor mutualistic symbiont B. clausii et rel., 100% of samples, and the opportunistic pathogen Proteus, 50% of samples. For each subject the relative IF contributions of the probes were calculated, obtaining an approximate evaluation of the relative abundance of the principal microbial groups of the faecal microbiota. In general agreement with previous metagenomic studies [7–11] and SSU rRNA phylogenetic microarray investigations [22, 23], mutualistic symbionts such as Bacteroidetes, Clostridium clusters IV, IX and XIVa largely dominated the faecal microbiota, contributing for the 65 to 80% of total microbiota, depending on the subject. Differently, with an overall contribution ranging from 10 to 30%, minor mutualistic symbionts such as B. clausii et rel., Bifidobacteriaceae, Lactobacillaceae, B. subtilis et rel., Fusobacterium, and Cyanobacteria were largely subdominant. Opportunistic pathogens represented only a small fraction of the intestinal microbiota. Even if subjects under study show a common trend when the ratio between the relative IF of major, minor and opportunistic components were considered, differences in the relative IF contribution of single probes were detectable and subject specific profiles were identified. For instance, subject n. 1 showed a higher relative fluorescence for probes targeting major mutualistic symbionts and a lower relative fluorescence for minor mutualistic symbionts and opportunistic pathogens than subjects n. 4 and 15. On the other hand subjects n. 15 and 17 were characterized by a lower ratio Bacteroidetes/Firmicutes with respect to all the other subjects. It is tempting to hypothesize that differences in relative IF contribution within samples could represent an approximation of differences in relative abundances of the targeted groups in the faecal microbiota. However, caution must be taken when microarray based methods for the relative quantification of bacterial groups in complex microbial communities are used. In fact, biases are introduced at several levels of the experimental procedure: DNA extraction and purification, PCR amplification of the 16S rRNA gene, and interspecies variation of the rRNA gene copy number .
The HTF-Microbi.Array has been revealed a fast and sensitive tool for the high taxonomic level fingerprint of the human intestinal microbiota in terms of presence/absence of the principal groups. Since the flexibility of the universal array platform allow the addition of new probe pairs without a further optimization of the hybridization conditions [25, 26], the HTF-Microbi.Array can be easy implemented with the addition of new probe pairs targeting emerging microbial groups of the human intestinal microbiota, such as, for instance, the mucin degrading bacterium Akkermansia muciniphila. The evaluation of the relative abundance of the target groups on the bases of the relative IF probes response still has some hindrances. However, considered all the possible biases (i.e. DNA extraction/purification, PCR, copy number variations, etc.) typical of the microarray technology, analysis of IFs from our LDR-UA platform can be useful in the estimation of the relative abundance of the targets groups within each sample. Focusing the phylogenetic resolution at division, order and cluster levels, the HTF-Microbi.Array results blind with respect to the inter-individual variability at the species level. Its potential to characterize the high order taxonomic unbalances of the human intestinal microbiota associated with specific diseases will be assessed in further studies.
Eight healthy Italian individuals of 30 years old were enrolled for the study. None of the subjects had dietary restrictions except for antibiotics, probiotics and functional foods for at least 4 weeks prior to sampling. None of the selected subjects had a history of gastrointestinal disorders at the time of sampling. The study protocol was approved by the Ethical committee of Sant'Orsola-Malpighi Hospital (Bologna, Italy) and an informed consent was obtained from each enrolled subject. Faeces were collected for each subject and stored at -20°C.
Bacterial strains and culture conditions
The bifidobacterial strains used in this study were Bifidobacterium adolescentis ATCC15703, B. bifidum DSM20456, B. breve DSM20091, B. longum ATCC15707. The Lactobacillus strains were Lactobacillus plantarum DSM21074, L. casei DSM20011, L. ramnosus DSM20021, L. salivarius SV2 (strain from our collection), L. delbrueckii DSM 20314, L. gasseri DSM20243, L. reuteri DSM20016, L. pentosus DSM20134, L. acidophilus DSM20079. All bifidobacteria and Lactobacillus strains were grown on De Man-Rogosa-Sharpe (MRS) broth with cysteine (0.5 g/l) at 37°C under an anaerobic atmosphere (Anaerocult, Merck, Darmstadt, Germany). Escherichia coli ATCC11105 was cultivated at 37°C aerobically on TY-broth. Salmonella cholerasuis typhimurium and Yersinia enterocolitica-type, kindly provided by A. Essig, Dept. of Medical Microbiology, University of Ulm, Germany, were cultivated aerobically at 30°C on BHI-broth.
Target selection and consensus extraction
A database of 16S rRNA sequences was created by integration of the 16S rRNA database of the ARB Project (release February, 2005) (http://www.arb-home.de; ) with the database of the Ribosomal Database Project (RDP; release September, 2007) (http://rdp.cme.msu.edu/; [36, 37]). A phylogenetic tree was obtained in the ARB software, by using the neighbour-joining algorithm for the sequence alignment. The tree was used for the rational selection of phylogenetically related groups of bacteria belonging to the human intestinal microbiota which correspond to nodes of the phylogenetic tree (Additional file 1). Group specific consensus sequences were extracted, with a cut-off of 75% for base calling. Nucleotides which occurred at lower frequencies were replaced by the appropriate IUPAC ambiguity code.
Multiple alignment step of the selected sequences was performed in ClustalW . Since the taxonomic classification of the 30 groups selected for the probe design varied from species to phylum level, careful grouping of the sequences was performed for the multiple alignment step: (a) for higher level probes, only family/phylum consensus sequences were used as a negative set for probe design; (b) for genus/species level probes, only sequences belonging to other families/phyla were selected. All the LDR probe pairs were designed using ORMA . Both DS and CP were required to be between 25 and 60 bases pair, with a Tm of 68 ± 1°C, and with maximum 4 degenerated bases. In-silico check versus a publicly available database (i.e.: RDP) was then performed for assessing probe pair specificity.
Total DNA was extracted from 109 bacterial cells by using the DNeasy Tissue Kit 50 (Quiagen, Düsseldorf, Germany) following the manufacturer instructions. Bacterial DNA was also extracted from lyophilized bacterial cells of the following DSMZ (Braunschweig, Germany) collection strains: Clostridium leptum DSM73, Ruminococcus albus DSM20455, Eubacterium siraeum DSM15700, C. viride DSM6836, Megasphera micrinuciformis DSM17226, Bacillus clausii DSM2515, B. subtilis DSM704, B. cereus DSM21, and Proteus mirabilis DSM4479. Lyophilized bacterial cells were suspended in 1 ml of lysis buffer (500 mM NaCl, 50 mM Tris-HCl pH 8, 50 mM EDTA, 4% SDS) and DNA extraction was carried out by employing the same procedure used for the extraction of genomic DNA from faecal samples, according to the following procedure. Total DNA from faecal material was extracted using QIAamp DNA Stool Min Kit (Qiagen) with a modified protocol. 250 mg of faeces were suspended in 1 ml of lysis buffer. Four 3 mm glass beads and 0.5 g of 0.1 mm zirconia beads were added, and the samples were treated in FastPrep (MP Biomedical, Irvine, CA, USA) at 5.5 ms for 3 min. Samples were heated at 95°C for 15 minutes, then centrifuged for 5 min at full speed to pellet stool particles. Supernatants were collected and 260 μl of 10 M ammonium acetate were added, followed by incubation in ice for 5 min and centrifugation at full speed for 10 min. One volume of isopropanol was added to each supernatant and incubated in ice for 30 min. The precipitated nucleic acids were collected by centrifugation for 15 min at full speed and washed with 70% ethanol. Pellets were resuspended in 100 μl of TE buffer and treated with 2 μl of DNase-free RNase (10 mg/ml) at 37°C for 15 min. Protein removal by Proteinase K treatment and DNA purification with QIAamp Mini Spin columns were performed following the kit protocol. 200 μl of TE buffer were used for DNA elution. Final DNA concentration was determined by using NanoDrop ND-1000 (NanoDrop Technologies, Wilmington, DE). The bacterial DNA from the following 11 ATCC strains was directly obtained from the ATCC: Bacteroides fragilis ATCC25285, B. thetaiotaomicron ATCC29148, Prevotella melaninogenica ATCC25845, Veilonella parvula ATCC10790, C. difficile ATCCBAA1382, C. acetobutilicum ATCC824, C. perfringens ATCC13124, Enterococcus faecalis ATCC700802, E. faecium ATCC51559, Campylobacter jejuni ATCC33292, R. productus 23340.
Polymerase Chain Reaction (PCR)
All the oligonucleotide primers and probe pairs were synthesized by Thermo Electron (Ulm, Germany). PCR amplifications were performed with Biometra Thermal Cycler II and Biometra Thermal Cycler T Gradient (Biometra, Germany). PCR products were purified by using a Wizard SV gel and PCR clean-up System purification kit (Promega Italia, Milan, Italy), according to the manufacturer's instructions, eluted in 20 μl of sterile water, and quantified with the DNA 7500 LabChip Assay kit and BioAnalyzer 2100 (Agilent Technologies, Palo Alto, CA, USA). 16S rRNA was amplified using universal forward primer 16S27F (5'-AGAGTTTGATCMTGGCTCAG-3') and reverse primer r1492 (5'-TACGGYTACCTTGTTACGACTT-3'), following the protocol described in Castiglioni et al. except for using 50 ng of starting DNA and 0.5 U of DNAzyme DNA polymerase II (Finnzymes, Espoo, Finland).
LDR/Universal Array approach
Phenylen-diisothiocyanate (PDITC) activated chitosan glass slides were used as surfaces for the preparation of universal arrays , comprising a total of 49 Zip-codes. Hybridization controls (cZip 66 oligonucleotide, complementary to zip 66, 5'-Cy3-GTTACCGCTGGTGCTGCCGCCGGTA-3') were used to locate the submatrixes during the scanning. The entire experimental procedure for both the chemical treatment and the spotting is described in detail in Consolandi et al.. An overview of the Universal Array layout and ZipCodes is provided as Additional file 6. Ligase Detection Reaction and hybridization of the products on the universal arrays were performed according to the protocol described in Castiglioni et al., except for the probe annealing temperature, set at 60°C.
The LDRs were carried out in a final volume of 20 μl with different quantities of purified PCR products: a) all LDRs for specificity tests were performed on 50 fmol of initial PCR product, for having no issues related to target; b) sensitivity tests were performed with decreasing PCR product concentration from 75 to 0.7 fmol; c) relative abundance tests were performed on 1 fmol E. coli PCR amplicon, mixed with human genomic DNA extracted from whole blood, at decreasing concentrations, from 4%, down to 0.02%; d) LDR experiments on the eight faecal samples were performed on 50 fmol of PCR product.
All arrays were scanned with ScanArray 5000 scanner (Perkin Elmer Life Sciences, Boston, MA, USA), at 10 μm resolution. In the experiments, the fluorescent images were obtained with different acquisition parameters on both laser power and photo-multiplier gain, in order to avoid saturation. IF were quantitated by ScanArray Express 3.0 software, using the "Adaptive circle" option, letting diameters vary from 60 to 300 μm. No normalization procedures on the IFs have been performed. To assess whether a probe pair was significantly above the background (i.e. was "present" or not), we performed a one-sided t-test (α = 0.01). The criteria was relaxed to α = 0.05 for sensitivity tests. The null distribution was set as the population of "Blank" spots (e.g. with no oligonucleotide spotted, n = 6). Two times the standard deviation of pixel intensities of the same spots was added to obtain a conservative estimate. For each zip-code, we considered the population of the IFs of all the replicates (n = 4) and tested it for being significantly above the null-distribution (H0: μtest = μnull; H1: μtest>μnull). In case one replicate in the test population was below 2.5 times the distribution mean, this was considered an outlier and was discarded from the analyses. We calculated the ratio between the signal intensities of the specific probes on the blank intensity (SNRs) and the ratio between all the other probes and the blank intensity (SNRns).
Hierarchical clustering of HTF-Microbi.Array profiles was carried out using the statistical software R http://www.r-project.org. The Euclidean distance among sample profiles was calculated and Ward's method was used for agglomeration.
This work was funded by the Micro(bi)array project of the University of Bologna, Italy. Our thanks to Maria Vurchio for help with administrative issues and to Giada Caredda for the support in the experimental phase.
- Gill SR, Pop M, Deboy RT, Eckburg PB, Turnbaugh PJ, Samuel BS, Gordon JI, Relman DA, Fraser-Liggett CM, Nelson KE: Metagenomic analysis of the human distal gut microbiome. Science. 2006, 312 (5778): 1355-9. 10.1126/science.1124234.PubMed CentralView ArticlePubMedGoogle Scholar
- Ley RE, Hamady M, Lozupone C, Turnbaugh PJ, Ramey RR, Bircher JS, Schlegel ML, Tucker TA, Schrenzel MD, Knight R, Gordon JI: Evolution of mammals and their gut microbes. Science. 2008, 320 (5883): 1647-51. 10.1126/science.1155725.PubMed CentralView ArticlePubMedGoogle Scholar
- Turnbaugh PJ, Ley RE, Hamady M, Fraser-Liggett CM, Knight R, Gordon JI: The human microbiome project. Nature. 2007, 449 (7164): 804-10. 10.1038/nature06244.PubMed CentralView ArticlePubMedGoogle Scholar
- Egert M, de Graaf AA, Smidt H, de Vos WM, Venema K: Beyond diversity: functional microbiomics of the human colon. Trends Microbiol. 2006, 14 (2): 86-91. 10.1016/j.tim.2005.12.007.View ArticlePubMedGoogle Scholar
- Neish AS: Microbes in gastrointestinal health and disease. Gastroenterology. 2009, 136 (1): 65-80. 10.1053/j.gastro.2008.10.080.PubMed CentralView ArticlePubMedGoogle Scholar
- Round JL, Mazmanian SK: The gut microbiota shapes intestinal immune responses during health and disease. Nat Rev Immunol. 2009, 9 (5): 313-23. 10.1038/nri2515.PubMed CentralView ArticlePubMedGoogle Scholar
- Peterson DA, Frank DN, Pace NR, Gordon JI: Metagenomic approaches for defining the pathogenesis of inflammatory bowel diseases. Cell Host Microbe. 2008, 3 (6): 417-27. 10.1016/j.chom.2008.05.001.PubMed CentralView ArticlePubMedGoogle Scholar
- Hattori M, Taylor TD: The human intestinal microbiome: a new frontier of human biology. DNA Res. 2009, 16 (1): 1-12. 10.1093/dnares/dsn033.PubMed CentralView ArticlePubMedGoogle Scholar
- Eckburg PB, Bik EM, Bernstein CN, Purdom E, Dethlefsen L, Sargent M, Gill SR, Nelson KE, Relman DA: Diversity of the human intestinal flora. Science. 2005, 308: 1635-1638. 10.1126/science.1110591.PubMed CentralView ArticlePubMedGoogle Scholar
- Andersson AF, Lindberg M, Jakobsson H, Backhed F, Nyrén P, Engstrand L: Comparative analysis of human gut microbiota by barcoded pyrosequencing. PloS ONE. 2008, 3: e2836-10.1371/journal.pone.0002836.PubMed CentralView ArticlePubMedGoogle Scholar
- Claesson MJ, O'Sullivan O, Wang Q, Nikkila J, Marchesi JR, Smidt H, de Vos WM, O'Toole PW: Comparative analysis of pyrosequencing and a phylogenetic microarray for exploring microbial community structures in the human distal intestine. PLoS ONE. 2009, 4: e6669-10.1371/journal.pone.0006669.PubMed CentralView ArticlePubMedGoogle Scholar
- Turnbaugh PJ, Hamady M, Yatsunenko T, Cantarel BL, Duncan A, Ley RE, Sogin ML, Jones WJ, Roe BA, Affourtit JP, Egholm M, Henrissat B, Heath AC, Knight R, Gordon JI: A core gut microbiome in obese and lean twins. Nature. 2009, 457 (7228): 480-4. 10.1038/nature07540.PubMed CentralView ArticlePubMedGoogle Scholar
- Turnbaugh PJ, Ley RE, Mahowald MA, Magrini V, Mardis ER, Gordon JI: An obesity-associated gut microbiome with increased capacity for energy harvest. Nature. 2006, 444 (7122): 1027-31. 10.1038/nature05414.View ArticlePubMedGoogle Scholar
- Frank DN, St Amand AL, Feldman RA, Boedeker EC, Harpaz N, Pace NR: Molecular-phylogenetic characterization of microbial community imbalances in human inflammatory bowel diseases. Proc Natl Acad Sci USA. 2007, 104 (34): 13780-5. 10.1073/pnas.0706625104.PubMed CentralView ArticlePubMedGoogle Scholar
- Sokol H, Pigneur B, Watterlot L, Lakhdari O, Bermúdez-Humarán LG, Gratadoux JJ, Blugeon S, Bridonneau C, Furet JP, Corthier G, Grangette C, Vasquez N, Pochart P, Trugnan G, Thomas G, Blottière HM, Doré J, Marteau P, Seksik P, Langella P: Faecalibacterium prausnitzii is an anti-inflammatory commensal bacterium identified by gut microbiota analysis of Crohn disease patients. Proc Natl Acad Sci USA. 2008, 105 (43): 16731-6. 10.1073/pnas.0804812105.PubMed CentralView ArticlePubMedGoogle Scholar
- Stecher B, Robbiani R, Walker AW, Westendorf AM, Barthel M, Kremer M, Chaffron S, Macpherson AJ, Buer J, Parkhill J, Dougan G, von Mering C, Hardt WD: Salmonella enterica serovar typhimurium exploits inflammation to compete with the intestinal microbiota. PLoS Biol. 2007, 5 (10): 2177-89. 10.1371/journal.pbio.0050244.View ArticlePubMedGoogle Scholar
- Pédron T, Sansonetti P: Commensals, bacterial pathogens and intestinal inflammation: an intriguing ménage à trois. Cell Host Microbe. 2008, 3 (6): 344-7. 10.1016/j.chom.2008.05.010.View ArticlePubMedGoogle Scholar
- Mazmanian SK, Round JL, Kasper DL: A microbial symbiosis factor prevents intestinal inflammatory disease. Nature. 2008, 453 (7195): 620-5. 10.1038/nature07008.View ArticlePubMedGoogle Scholar
- Hamady M, Knight R: Microbial community profiling for human microbiome projects: Tools, techniques, and challenges. Genome Res. 2009, 19 (7): 1141-52. 10.1101/gr.085464.108.PubMed CentralView ArticlePubMedGoogle Scholar
- Flint HJ, Duncan SH, Scott KP, Louis P: Interactions and competition within the microbial community of the human colon: links between diet and health. Environ Microbiol. 2007, 9 (5): 1101-11. 10.1111/j.1462-2920.2007.01281.x.View ArticlePubMedGoogle Scholar
- Palmer C, Bik EM, Eisen MB, Eckburg PB, Sana TR, Wolber PK, Relman DA, Brown PO: Rapid quantitative profiling of complex microbial populations. Nucleic Acids Res. 2006, 34 (1): e5-10.1093/nar/gnj007.PubMed CentralView ArticlePubMedGoogle Scholar
- Palmer C, Bik EM, DiGiulio DB, Relman DA, Brown PO: Development of the human infant intestinal microbiota. PLoS Biol. 2007, 5 (7): e177-10.1371/journal.pbio.0050177.PubMed CentralView ArticlePubMedGoogle Scholar
- Rajilić-Stojanović M, Heilig HG, Molenaar D, Kajander K, Surakka A, Smidt H, de Vos WM: Development and application of the human intestinal tract chip, a phylogenetic microarray: analysis of universally conserved phylotypes in the abundant microbiota of young and elderly adults. Environ Microbiol. 2009,Google Scholar
- Paliy O, Kenche H, Abernathy F, Michail S: High-throughput quantitative analysis of the human intestinal microbiota with a phylogenetic microarray. Appl Environ Microbiol. 2009, 75 (11): 3572-9. 10.1128/AEM.02764-08.PubMed CentralView ArticlePubMedGoogle Scholar
- Castiglioni B, Rizzi E, Frosini A, Sivonen K, Rajaniemi P, Rantala A, Mugnai MA, Ventura S, Wilmotte A, Boutte C, Grubisic S, Balthasart P, Consolandi C, Bordoni R, Mezzelani A, Battaglia C, De Bellis G: Development of a universal microarray based on the ligation detection reaction and 16S rrna gene polymorphism to target diversity of cyanobacteria. Appl Environ Microbiol. 2004, 70 (12): 7161-72. 10.1128/AEM.70.12.7161-7172.2004.PubMed CentralView ArticlePubMedGoogle Scholar
- Hultman J, Ritari J, Romantschuk M, Paulin L, Auvinen P: Universal ligation-detection-reaction microarray applied for compost microbes. BMC Microbiol. 2008, 30 (8): 237-10.1186/1471-2180-8-237.View ArticleGoogle Scholar
- Collins MD, Lawson PA, Willems A, Cordoba JJ, Fernandez-Garayzabal J, Garcia P, Cai J, Hippe H, Farrow JA: The phylogeny of the genus Clostridium: proposal of five new genera and eleven new species combinations. Int J Syst Bacteriol. 1994, 44 (4): 812-26. 10.1099/00207713-44-4-812.View ArticlePubMedGoogle Scholar
- Rajilić-Stojanović M, Smidt H, de Vos WM: Diversity of the human gastrointestinal tract microbiota revisited. Environ Microbiol. 2007, 9 (9): 2125-36. 10.1111/j.1462-2920.2007.01369.x.View ArticlePubMedGoogle Scholar
- Peplies J, Glöckner FO, Amann R: Optimization strategies for DNA microarray-based detection of bacteria with 16S rRNA-targeting oligonucleotide probes. Appl Environ Microbiol. 2003, 69 (3): 1397-407. 10.1128/AEM.69.3.1397-1407.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Jin LQ, Li JW, Wang SQ, Chao FH, Wang XW, Yuan ZQ: Detection and identification of intestinal pathogenic bacteria by hybridization to oligonucleotide microarrays. World J Gastroenterol. 2005, 11 (48): 7615-9.PubMedGoogle Scholar
- Severgnini M, Cremonesi P, Consolandi C, Caredda G, De Bellis G, Castiglioni B: ORMA: a tool for identification of species-specific variations in 16S rRNA gene and oligonucleotides design. Nucleic Acids Res. 2009, 37 (16): e109-10.1093/nar/gkp499.PubMed CentralView ArticlePubMedGoogle Scholar
- Edwards U, Rogall T, Blöcker H, Emde M, Böttger EC: Isolation and direct complete nucleotide determination of entire genes. Characterization of a gene coding for 16S ribosomal RNA. Nucleic Acids Res. 1989, 19: 7843-53. 10.1093/nar/17.19.7843.View ArticleGoogle Scholar
- Eckburg PB, Bik EM, Bernstein CN, Purdom E, Dethlefsen L, Sargent M, Gill SR, Nelson KE, Relman DA: Diversity of the human intestinal microbial flora. Science. 2005, 308 (5728): 1635-8. 10.1126/science.1110591.PubMed CentralView ArticlePubMedGoogle Scholar
- Derrien M, Collado MC, Ben-Amor K, Salminem S, de Vos WM: The mucin degrader Akkermansia muciniphila is an abundant resident of the human intestinal tract. Appl Environ Microbiol. 2008, 74 (5): 1646-48. 10.1128/AEM.01226-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Ludwig W, Strunk O, Westram R, Richter L, Meier H, Yadhukumar , Buchner A, Lai T, Steppi S, Jobb G, Förster W, Brettske I, Gerber S, Ginhart AW, Gross O, Grumann S, Hermann S, Jost R, König A, Liss T, Lüssmann R, May M, Nonhoff B, Reichel B, Strehlow R, Stamatakis A, Stuckmann N, Vilbig A, Lenke M, Ludwig T, Bode A, Schleifer KH: ARB: a software environment for sequence data. Nucleic Acids Res. 2004, 32 (4): 1363-71. 10.1093/nar/gkh293.PubMed CentralView ArticlePubMedGoogle Scholar
- Cole JR, Chai B, Farris RJ, Wang Q, Kulam-Syed-Mohideen AS, McGarrell DM, Bandela AM, Cardenas E, Garrity GM, Tiedje JM: The ribosomal database project (RDP-II): introducing myRDP space and quality controlled public data. Nucleic Acids Res. 2007, D169-72. 10.1093/nar/gkl889. 35 Database
- Wang Q, Garrity GM, Tiedje JM, Cole JR: Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007, 73 (16): 5261-7. 10.1128/AEM.00062-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Chenna R, Sugawara H, Koike T, Lopez R, Gibson TJ, Higgins DG, Thompson JD: Multiple sequence alignment with the Clustal series of programs. Nucleic Acids Res. 2003, 31 (13): 3497-500. 10.1093/nar/gkg500.PubMed CentralView ArticlePubMedGoogle Scholar
- Gerry NP, Witowski NE, Day J, Hammer RP, Barany G, Barany F: Universal DNA microarray method for multiplex detection of low abundance point mutations. J Mol Biol. 1999, 292 (2): 251-62. 10.1006/jmbi.1999.3063.View ArticlePubMedGoogle Scholar
- Consolandi C, Severgnini M, Castiglioni B, Bordoni R, Frosini A, Battaglia C, Rossi Bernardi L, De Bellis G: A structured chitosan-based platform for biomolecule attachment to solid surfaces: application to DNA microarray preparation. Bioconjug Chem. 2006, 17 (2): 371-77. 10.1021/bc050285a.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.