- Open Access
Clinical detection of human probiotics and human pathogenic bacteria by using a novel high-throughput platform based on next generation sequencing
Journal of Clinical Bioinformatics volume 4, Article number: 1 (2014)
The human body plays host to a vast array of bacteria, found in oral cavities, skin, gastrointestinal tract and the vagina. Some bacteria are harmful while others are beneficial to the host. Despite the availability of many methods to identify bacteria, most of them are only applicable to specific and cultivable bacteria and are also tedious. Based on high throughput sequencing technology, this work derives 16S rRNA sequences of bacteria and analyzes probiotics and pathogens species.
We constructed a database that recorded the species of probiotics and pathogens from literature, along with a modified Smith-Waterman algorithm for assigning the taxonomy of the sequenced 16S rRNA sequences. We also constructed a bacteria disease risk model for seven diseases based on 98 samples. Applicability of the proposed platform is demonstrated by collecting the microbiome in human gut of 13 samples.
The proposed platform provides a relatively easy means of identifying a certain amount of bacteria and their species (including uncultivable pathogens) for clinical microbiology applications. That is, detecting how probiotics and pathogens inhabit humans and how affect their health can significantly contribute to develop a diagnosis and treatment method.
High throughput sequencing can analyze a large amount of sequences, enabling sequencing of 16S rRNA to identify complex bacteria species of pathogens and probiotic bacteria. Many naturally occurring bacteria form a complex population in the environment. The human body plays host to a vast array of bacteria, found in oral cavities, skin, gastrointestinal tract and the vagina. Some bacteria are harmful while others are beneficial to the host.
A pathogen is a microorganism that causes disease in its host. For example, bacterial pathogen include Corynebacterium diphtheria (causes diphtheria), Listeria monocytogenes (causes food poisons), and Legionella pneumophila (causes Legionnaires’ disease). Probiotics, another microorganism, benefit the host and has received considerable attention in recent years. A FAO report in 2001  cited the advantages of probiotics as increasing immunity [2, 3], reducing gastrointestinal discomfort [4, 5], and protecting the flora within urogenital tract . As is well known, probiotics can ameliorate symptoms of diseases  and reduce the risk of suffering from diseases [8, 9].
Despite the availability of many approaches to identify probiotics and pathogens, most of them are only applicable to specific and cultivable bacteria but time consuming. For instance, conventional methods detect growth of cultured bacteria in approximately two days, or an additional five days to obtain no-growth culture results , which is laborious. Besides, some bacteria cannot be cultured , subsequently increasing the difficulty of specifying pathogenic bacteria. Moreover, it is hard to determine whether an infection is caused by one or more bacteria types.
16S rRNA sequences, capable of identifying bacteria on a molecular level, can detect uncultivable bacteria . Use of 16S rRNA sequencing can overcome some problems of conventional culture method . Although 16S rRNA sequencing is a more effective means of identifying bacteria than conventional culture method, 16S rRNA sequencing takes a considerable amount of time in amplifying DNA sequences . Sanger sequencing known as “first-generation” or “conventional” sequencing has been used for DNA sequencing for almost two decades. Next generation sequencing (NGS) can analyze large-scale sequences quicker, enable massively parallel analysis, reduce reagent costs and the size of sample components, and perform high throughput . Thus NGS is more efficient than the Sanger method, which generates one read per sample. In addition, NGS of 16S rRNA more easily identify cultivable or uncultivable bacteria .
Because of the improvement of sequencing technology and Bioinformatics approaches, the accuracy in distinguishing bacteria with those methods has been increased. Based on high throughput sequencing technology, this work identifies 16S rRNA sequences of bacteria and analyzes bacteria species. High-throughput sequencing can sequence a large number of 16S rRNA sequence more efficiently; with high-throughput sequencing, researchers can acquire information to identify pathogens and probiotic bacteria [16–18].
Platform application: gut probiotics and pathogens detection
The read statistics of quality filtering and taxonomy assignment are demonstrated in Table 1. Figure 1A illustrated the percentage of probiotics detected by the proposed platform. Table 2 listed the quantities (matched sequenced reads) of probiotics identified in the samples in the case study. The top three identified probiotics in 12 samples are Lactococcus salivarius, Streptococcus thermophilus, and Bifidobacterium longum. Figure 1B and Table 3 listed the proportion and quantities of pathogens, of which top three pathogens are Escherichia coli, Salmonella enteric, and Haemophilus influenza.
Table 4 listed the results of disease risk evaluations. It showed that three diseases of two samples (B031 and B034) had similar distributions in the control group. The three diseases are obesity, colorectal cancer, and constipation. Sample B031 had reached the significance level with P-value 0.0333 and 0.0121 < 0.05 of distribution in constipation and colorectal cancer respectively compared to 98 samples control group using binomial test. Sample B034 had reached the significance level with P-value 0.00257 and 0.0121 < 0.05 of distribution in obesity and colorectal cancer. Evaluated by the association of bacterial risk markers and disease, the results suggested that these two samples had higher risk than 98 samples control group in constipation, colorectal cancer, and obesity. Their enterotypes of gut probiotics and pathogens may be one of risk factors which would cause disease.
Reproducibility and accuracy evaluation of proposed platform
Two replicated experiments of four samples were performed to estimate the reproducibility of the proposed platform. The results of repeated experiments were consistent. The similarity between two repeated experiments were calculated by using UniFrac . Results of each sample pair (replicate 1 and 2) closely resemble each other. The similarity of UniFrac distance of each sample pair is higher than 0.96 (0.9617 for B014, 0.9872 for B018, 0.9914 for B020, 0.9722 for B033). This implies that the analysis results are reproducible.
Next, accuracy of the platform is evaluated by adding Lactobacillus reuteri to a stool sample (B050). Sample B050 contains 24,408 assigned taxons, and Lactobacillus reuteri has no detected count. Whether the counts of this species in positive control sample (B050S_L) are elevated must be determined. Analysis results indicate that 27,113 taxons are detected in sample B050S_L. In fact, the detected counts of Lactobacillus reuteri in sample B050S_L are 1,430, and the percentage of Lactobacillus reuteri markedly increases from 0% to 5%.
In short, our platform is accurate and reproducible in terms of detecting the quantities of bacterial species of the proposed platform. The results evaluate the accuracy and feasibility of proposed platform in order to identify probiotics and pathogens. While requiring only about one day for detection, not limited in identifying certain bacteria, the proposed platform can detect and quantify multiple bacteria simultaneously.
Because of the constraint of costs and technical limitations, 16S rRNA sequences obtained in most databases are partial sequences. Many studies thus assign taxonomy by using partial 16S rRNA sequences. In our probiotics and pathogens 16S rRNA sequence database, 17,964 sequences are collected from NCBI nucleotide database, NCBI 16S microbial rRNA database, Greengenes database, and SILVA. Our probiotics and pathogens 16S rRNA database contain less than 39% of 16S rRNA sequences which are longer than 1400 bps. Only 9% of the sequences are close to full length.
This work extracts the V4 region from full length 16S rRNA of microbiome in the human gut as a platform application. Some 16S rRNA variable regions are more dependable than other regions in assigning taxonomy like V3 and V4 [20, 21]; in addition, some 16S rRNA variable regions are much conserved. The proportion and diversity of probiotics and pathogens may be made diverse by using different 16S rRNA variable regions. The proposed platform is also applicable to other 16S rRNA variable regions for taxonomy assignment. Importantly, a more appropriate region than others must be selected to produce an outcome that is close to full length 16S rRNA sequence.
This work further attempt is to collect common probiotics and pathogens from the literature. Although it may be incomplete, recent advances in sequencing technology make it possible to identify and define an increasing number of bacteria, implying an obvious increase in the number of identified probiotics and pathogens in the future. Efforts are underway in our laboratory to update the list of used probiotics and pathogens.
Previous studies [22–24] identified pathogen or probiotic bacteria by using antibody, 16S rRNA gene microarrays, fluorescence in situ hybridization (FISH), and proteomic methods. In this work, the proposed platform can detect various pathogens and probiotics based on 16S rRNA (rDNA) sequences of bacteria using NGS and Bioinformatics method. An average of 126,451 reads was acquired per sample in this work. It is doubt that the sequencing depth is enough to detect a small amount of probiotics and pathogens. Although increasing the coverage of sequencing can advance the sensitivity of detecting probiotics and pathogens, the sequencing cost will increase. It is important to work out an appropriate coverage of sequencing for detecting probiotics and pathogens.
The results of disease risk evaluations revealed that most of 12 samples did not have resembled distributions of bacteria markers with control group. Only two samples had reached the significance level of distributions. The reason for the phenomenon may be the overlapped bacteria markers between diseases. 28 markers are used in colorectal cancer, and 17 markers are used in irritable bowel syndrome. Six markers are overlapped. For sample B031, the significant distributions in colorectal cancer were partly contributed to the significance in irritable bowel syndrome owing to the overlapped markers. Similarly, two overlapped markers for sample B034 were in colorectal cancer and obesity. In this kind of speculation, the influence of colorectal cancer to irritable bowel syndrome would be six (overlapped markers of CC and IBS) over seventeen (markers of IBS), and the influence of colorectal cancer to obesity would be two (overlapped markers of CC and obesity) over nine (markers of obesity). In addition, the influence of colorectal cancer to constipation and ulcerative colitis would be one over six and two over ten, respectively.
In addition to that some bacteria markers in species level are belong to the marker of genus level and species level, both genus marker and species markers may have associated with affecting the distributions mutually. Continually, collecting more markers and evaluating the distributions with markers in the same level are required for constructing a global prediction model in Taiwanese.
This work constructed a bacterial disease risk evaluation model for seven diseases and developed a novel platform by using NGS and Bioinformatics approach. Compared with the traditional bacteria culture method, our proposed platform can reduce experiment time. Besides, the proportion of probiotics and pathogens (including uncultivable pathogens) in the human body can be detected rapidly with 16 s RNA database of probiotics and pathogens. Furthermore, the proposed platform provides further insight into the cause of disease based on the relation of probiotics, pathogens, and disease. For instance, the type of antibiotics can be adjusted if the pathogens of disease are identified from infected patients. In addition, the proposed platform allows researchers to determine whether the intake of probiotics impacts the human body [25–29]. In the future, this preliminary study will be continuously extended for more bacterial disease markers. For more comprehensive applications, this work will also collect bacteria from other parts of human body as control group data. In fact, a detective method of how the probiotics and pathogens inhabit human can provide new insight for human health. It could improve diagnosis and treatment method.
Figure 2 illustrates the bioinformatics system flow of the proposed platform, which includes analysis pipeline of NGS. The Figure 2 contains four parts: sequence quality filtering, construction of bacteria sequence database, taxonomy assignment, and disease risk model evaluation. The detailed components in the proposed platform are described below.
In this study, stool samples of 98 Taiwan volunteers were gathered. The samples were collected by Sigma-transwab (Medical Wire) into a tube with Liquid Amies Transport Medium, and stored at 4°C until processing.
In the case study, fresh faeces were obtained from participants. DNA was extracted directly on stool samples by using a QIAamp DNA Stool Mini Kit (Qiagen). A swab was vortexed vigorously and incubated at room temperature for 1 min. The sample was then transferred to microcentrifuge tubes containing 560 μl Buffer ASL, vortexed, and incubated at 37°C for 30 min. In addition, the suspension was incubated at 95°C for 15 min, vortexed, and centrifuged at 14,000 rpm for 1 min into pellet stool particles. Extraction was performed following the protocol of the QIAamp DNA Stool Mini Kit. The DNA was eluted with 50 μl Buffer AE, and centrifuged at 14,000 rpm for 1 min. Moreover, the DNA extract was stored at-20°C until further analysis. Finally, DNA extraction was performed, depending on the sample collected.
Library construction and sequencing for V4 region of 16S ribosomal DNA
The PCR primers, F515 (5′-GTGCCAGCMGCCGCGGTAA-3′) and R806 (5′-GGACTACHVGGGTWTCTAAT-3′), were designed to amplify the V4 domain of bacterial 16S ribosomal DNA as described previously . PCR amplification was performed in a 50 μl reaction volume containing 25 μl 2X Taq Master Mix (Thermo Scientific), 0.2 μM of each forward and reverse primer, and 20 ng DNA template. The reaction conditions consisted of an initial 95°C for 5 min, followed by 30 cycles of 95°C for 30 sec, 54°C for 1 min, and 72°C for 1 min, as well as a final extension of 72°C for 5 min. Next, amplified products were checked by 2% agarose gel electrophoresis and ethidium bromide staining. Amplicons were purified using the AMPure XP PCR Purification Kit (Agencourt), and quantified using Qubit dsDNA HS Assay Kit (Qubit) on Qubit 2.0 Fluorometer (Qubit)–all according to respective manufacturer instructions. For V4 library preparation, Illumina adapters were attached to the amplicons using the Illumina TruSeq DNA Sample Preparation v2 Kit. Purified libraries were applied for cluster generation and sequencing on the MiSeq system. The raw sequence files are available for download at http://clinic.mbc.nctu.edu.tw/.
16S rRNA (rDNA) sequence data quality filtering
The raw fastq files obtained by Illumina sequencing machine were quality-filtered using the FASTX-Toolkita. The paired-end 150 bp reads were performed using the minimum acceptable phred quality score of 20, as well as the 70% of bases that must exceed 20 phred quality score. Sequence shorter than 100 nucleotides would be omitted after quality trimming from reads tail. Notably, reads containing ambiguous characters were discarded.
Construction of probiotics and pathogens database
The list of probiotics and pathogens were obtained from literatures or the claims of official departments. Additional file 1: Table S1 lists species of probiotics which were adapted from both literatures [7, 9] and the claims of official departments, such as Taiwan Food and Drug Administration  and Health Canada . 99 bacterial pathogens were collected from literature [25, 26, 33–42] and Taiwan Food and Drug Administration  (Additional file 1: Table S2).
The 16S rRNA sequences of probiotics and pathogens used for taxonomy mapping were retrieved from the NCBI nucleotide database, NCBI 16S microbial rRNA database, Greengenes database  and SILVA . Following sequence data collection, we assemble partial sequences which used the same species classification and removed redundant sequences. Additionally, we also removed the unique sequence from only one research support with 3% similarity which shared the same species classification with other sequence.
To generate taxonomy assignments, the proposed platform invoked a modified Smith-Waterman algorithm from miRExpress , which can compare pairs of sequences in parallel, for mapping reads to taxons. miRExpress was designed for identifying the best similarity between sequencing reads and miRNA precursor sequences. In our model, it was modified for identifying multiple hits of 16S rRNA sequence mapping results with similarity threshold 0.97. In order to reduce the storage space of output, the SAM format  was used to replace the original miRExpress output format for storing alignment results. Furthermore, two kinds of output format were designed. One format records whole mapped sequencing reads based on taxons. The other one records which taxons could be assigned based on sequencing reads. These two kinds of output could support the important information for assigning sequencing reads to suitable taxon. miRExpress was originally designed for dealing with single-end sequencing data. Therefore, the additional program was added for processing paired-end sequencing data. In this part, both end sequencing reads need to be assigned to the same taxon. If paired-end sequencing reads were mapped to different taxons, this paired sequence would be dropped. The probiotics and pathogens 16S rRNA sequence from our database were built in FASTA format. Following quality filtering, all paired-end sequences were aligned to the probiotics and pathogens database with whole read aligned from one end to the other end. Reads were then truncated with an identity lower than 97%, according to previous research in order to achieve a better compromise between sequences from PCR sequencing errors and taxonomic relatedness .
The construction of Bacterial disease risk evaluation model (BDREM)
To study the associations between bacteria and diseases, we collected related information from literatures. We concerned bacteria that are associated with seven diseases: constipation [28, 47, 48], obesity [29, 49–52], irritable bowel syndrome (IBS) [28, 53–58], ulcerative colitis (UC) [53, 59–61], colon cancer (CC) [62–64], Atopic Dermatitis (AD) and Allergic rhinitis (AR), were collected positive correlation and negative correlation data, and the individual risk of disease was evaluated.
The association data were majorly collected from case–control studies which the quantities of bacteria were obtained from NGS data, and few well-known bacteria validated by multiple studies through cultural experiments were also included. We further eliminated some conflicted data with both positive and negative correlation between bacteria and disease in different studies.
Health Asians stool samples of 98 Taiwan volunteers were gathered. Following deep sequencing and sequencing data processing, the proportion of 78 bacteria from control group was applied as risk markers (constipation: 6, obesity: 9, IBS: 17, UC: 10, CC: 28, AD: 4, AR: 4) to predict disease risk to seven diseases in this study (Table 5).
The mathematical formula of BDREM in this study was developed as the following steps. Let λ be a N × S matrix, where N is the number of markers selected in the prediction model of constipation and S is the number of health subjects in 7 prediction models. T i was defined as one of the two notches of median for each row of λ. T i is a threshold to distinguish λ ij from normal proportion level to abnormal (fail to success in one trail of binomial distribution). Smaller notch was selected to T i when each marker was recorded as a negative association to the disease, and a success trail was identified when λ ij is smaller than T i . On the opposite, larger notch was selected when association was positive, and a success trail was identified when λ ij is larger than T i .
Let P j be the probability of successful trails in the jth column of λ. The meaning of P j is the personal probability that abnormal proportion level happened.
Let P h be the mean of P j . It represents how frequent the abnormal proportion level happened to all P j in average, regarded as the hypothesized probability of success in each P j .
Assume P j obey a binomial distribution, and let P h be the hypothesized probability (0.05051 for constipation, 0.07239 for obesity, 0.06952 for IBS, 0.05227 for UC, 0.09280 for CC, 0.04924 for AD, 0.05114 for AR). A binomial test was used to P j and P h . Alpha = 0.05 was choose to judge if a subject is significantly differently from the others in λ.
Figure 3 illustrated an example for evaluating the risk of obesity of B034 and B031. The model used lower and upper proportion bound of 9 markers from 98 control samples to define risk markers of these two samples following by using binomial test. Four markers of B034 exceed the lower bound and upper bound of obesity. The binomial test P-Value of B034 is 2.572e-03 < 0.05, Since P-Value < = hypothesized probability 0.07239, this case is specifically associated (significantly) with disease than random chance. There are two markers of B031 exceed lower bound of obesity. The P-Value of B031 is 0.1344 > 0.05, the case is no more associated with disease than random chance. As the results, we can assume that B034 had higher probability to cause Obesity.
FAO/WHO: Health and nutritional properties of probiotics in food including powder milk with live lactic acid bacteria. Book Health and nutritional properties of probiotics in food including powder milk with live lactic acid bacteria. http://www.who.int/foodsafety/publications/fs_management/en/probiotics.pdf. City; 2001
Kalliomaki M, Salminen S, Arvilommi H, Kero P, Koskinen P, Isolauri E: Probiotics in primary prevention of atopic disease: a randomised placebo-controlled trial. Lancet. 2001, 357: 1076-1079. 10.1016/S0140-6736(00)04259-8.
Isolauri E, Sutas Y, Kankaanpaa P, Arvilommi H, Salminen S: Probiotics: effects on immunity. Am J Clin Nutr. 2001, 73: 444S-450S.
Gopal PK, Prasad J, Smart J, Gill HS: In vitro adherence properties of lactobacillus rhamnosus DR20 and bifidobacterium lactis DR10 strains and their antagonistic activity against an enterotoxigenic Escherichia coli. Int J Food Microbiol. 2001, 67: 207-216. 10.1016/S0168-1605(01)00440-8.
Ogawa M, Shimizu K, Nomoto K, Takahashi M, Watanuki M, Tanaka R, Tanaka T, Hamabata T, Yamasaki S, Takeda Y: Protective effect of Lactobacillus casei strain Shirota on Shiga toxin-producing Escherichia coli O157:H7 infection in infant rabbits. Infect Immun. 2001, 69: 1101-1108. 10.1128/IAI.69.2.1101-1108.2001.
Reid G: Probiotic agents to protect the urogenital tract against infection. Am J Clin Nutr. 2001, 73: 437S-443S.
Wan YM, Zhu YQ, Xia B, Luo J: Treating TNBS-induced colitis in rats with probiotics. Turk J Gastroenterol. 2011, 22: 486-493.
Taipale T, Pienihakkinen K, Isolauri E, Larsen C, Brockmann E, Alanen P, Jokela J, Soderling E: Bifidobacterium animalis subsp. lactis BB-12 in reducing the risk of infections in infancy. Br J Nutr. 2011, 105: 409-416. 10.1017/S0007114510003685.
Wickens K, Black PN, Stanley TV, Mitchell E, Fitzharris P, Tannock GW, Purdie G, Crane J: A differential effect of 2 probiotics in the prevention of eczema and atopy: a double-blind, randomized, placebo-controlled trial. J Allergy Clin Immunol. 2008, 122: 788-794. 10.1016/j.jaci.2008.07.011.
Jordan JA, Durso MB: Real-time polymerase chain reaction for detecting bacterial DNA directly from blood of neonates being evaluated for sepsis. J Mol Diagn. 2005, 7: 575-581. 10.1016/S1525-1578(10)60590-9.
Wilson M: Bacteriology of Humans: An Ecological Perspective. 2008, Hoboken: Wiley-Blackwell
Goldenberger D, Kunzli A, Vogt P, Zbinden R, Altwegg M: Molecular diagnosis of bacterial endocarditis by broad-range PCR amplification and direct sequencing. J Clin Microbiol. 1997, 35: 2733-2739.
Woo PC, Lau SK, Woo GK, Fung AM, Ngan AH, Hui WT, Yuen KY: Seronegative bacteremic melioidosis caused by Burkholderia pseudomallei with ambiguous biochemical profile: clinical importance of accurate identification by 16S rRNA gene and groEL gene sequencing. J Clin Microbiol. 2003, 41: 3973-3977. 10.1128/JCM.41.8.3973-3977.2003.
Fox GE, Magrum LJ, Balch WE, Wolfe RS, Woese CR: Classification of methanogenic bacteria by 16S ribosomal RNA characterization. Proc Natl Acad Sci U S A. 1977, 74: 4537-4541. 10.1073/pnas.74.10.4537.
Metzker ML: Sequencing technologies–the next generation. Nat Rev Genet. 2010, 11: 31-46. 10.1038/nrg2626.
Kuczynski J, Lauber CL, Walters WA, Parfrey LW, Clemente JC, Gevers D, Knight R: Experimental and analytical tools for studying the human microbiome. Nat Rev Genet. 2012, 13: 47-58.
Woo PC, Lau SK, Teng JL, Tse H, Yuen KY: Then and now: use of 16S rDNA gene sequencing for bacterial identification and discovery of novel bacteria in clinical microbiology laboratories. Clin Microbiol Infect. 2008, 14: 908-934. 10.1111/j.1469-0691.2008.02070.x.
Lazarevic V, Whiteson K, Huse S, Hernandez D, Farinelli L, Osteras M, Schrenzel J, Francois P: Metagenomic study of the oral microbiota by Illumina high-throughput sequencing. J Microbiol Methods. 2009, 79: 266-271. 10.1016/j.mimet.2009.09.012.
Lozupone C, Hamady M, Knight R: UniFrac–an online tool for comparing microbial community diversity in a phylogenetic context. BMC Bioinformatics. 2006, 7: 371-10.1186/1471-2105-7-371.
Bokulich NA, Joseph CM, Allen G, Benson AK, Mills DA: Next-generation sequencing reveals significant bacterial diversity of botrytized wine. PloS one. 2012, 7: e36357-10.1371/journal.pone.0036357.
Claesson MJ, Wang Q, O’Sullivan O, Greene-Diniz R, Cole JR, Ross RP, O’Toole PW: Comparison of two next-generation sequencing technologies for resolving highly complex microbiota composition using tandem variable 16S rRNA gene regions. Nucleic Acids Res. 2010, 38: e200-10.1093/nar/gkq873.
Russmann H, Kempf VA, Koletzko S, Heesemann J, Autenrieth IB: Comparison of fluorescent in situ hybridization and conventional culturing for detection of Helicobacter pylori in gastric biopsy specimens. J Clin Microbiol. 2001, 39: 304-308. 10.1128/JCM.39.1.304-308.2001.
Barber JJ, Grichting WL: Australia’s media campaign against drug abuse. Int J Addict. 1990, 25: 693-708.
Rigsbee L, Agans R, Foy BD, Paliy O: Optimizing the analysis of human intestinal microbiota with phylogenetic microarray. FEMS Microbiol Ecol. 2011, 75: 332-342. 10.1111/j.1574-6941.2010.01009.x.
Maley MW, Kociuba K, Chan RC: Prevention of laboratory-acquired brucellosis: significant side effects of prophylaxis. Clin Infect Dis. 2006, 42: 433-434. 10.1086/499112.
Mantur BG, Amarnath SK, Shinde RS: Review of clinical and laboratory features of human brucellosis. Indian J Med Microbiol. 2007, 25: 188-202. 10.4103/0255-0857.34758.
Hummelen R, Fernandes AD, Macklaim JM, Dickson RJ, Changalucha J, Gloor GB, Reid G: Deep sequencing of the vaginal microbiota of women with HIV. PloS one. 2010, 5: e12078-10.1371/journal.pone.0012078.
Lyra A, Rinttila T, Nikkila J, Krogius-Kurikka L, Kajander K, Malinen E, Matto J, Makela L, Palva A: Diarrhoea-predominant irritable bowel syndrome distinguishable by 16S rRNA gene phylotype quantification. World J Gastroenterol. 2009, 15: 5936-5945. 10.3748/wjg.15.5936.
Furet JP, Kong LC, Tap J, Poitou C, Basdevant A, Bouillot JL, Mariat D, Corthier G, Dore J, Henegar C, et al: Differential adaptation of human gut microbiota to bariatric surgery-induced weight loss: links with metabolic and low-grade inflammation markers. Diabetes. 2010, 59: 3049-3057. 10.2337/db10-0253.
Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Lozupone CA, Turnbaugh PJ, Fierer N, Knight R: Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc Natl Acad Sci U S A. 2011, 108 (Suppl 1): 4516-4522.
Khamis A, Raoult D, La Scola B: Comparison between rpoB and 16S rRNA gene sequencing for molecular identification of 168 clinical isolates of corynebacterium. J Clin Microbiol. 2005, 43: 1934-1936. 10.1128/JCM.43.4.1934-1936.2005.
Lin J: Effect of antibiotic growth promoters on intestinal microbiota in food animals: a novel model for studying the relationship between gut microbiota and human obesity?. Front Microbiol. 2011, 2: 53-
Hegerle N, Guiso N: Epidemiology of whooping cough and typing of Bordetella pertussis. Future Microbiol. 2013, 8 (11): 1391-1403. 10.2217/fmb.13.111.
Seleem MN, Boyle SM, Sriranganathan N: Brucella: a pathogen without classic virulence genes. Vet Microbiol. 2008, 129: 1-14. 10.1016/j.vetmic.2007.11.023.
Zhan P, Suo LJ, Qian Q, Shen XK, Qiu LX, Yu LK, Song Y: Chlamydia pneumoniae infection and lung cancer risk: a meta-analysis. Eur J Cancer. 2011, 47: 742-747. 10.1016/j.ejca.2010.11.003.
Parkhill J, Wren BW, Mungall K, Ketley JM, Churcher C, Basham D, Chillingworth T, Davies RM, Feltwell T, Holroyd S, et al: The genome sequence of the food-borne pathogen campylobacter jejuni reveals hypervariable sequences. Nature. 2000, 403: 665-668. 10.1038/35001088.
Guidoboni M, Ferreri AJ, Ponzoni M, Doglioni C, Dolcetti R: Infectious agents in mucosa-associated lymphoid tissue-type lymphomas: pathogenic role and therapeutic perspectives. Clin Lymphoma Myeloma. 2006, 6: 289-300. 10.3816/CLM.2006.n.003.
Harvey RA, Champe PC, Fisher BD: Lippincott’s Illustrated Reviews: Microbiology. 2007, Philadelphia: Lippincott Williams & Wilkins, 2
Fredlund H, Falk L, Jurstrand M, Unemo M: Molecular genetic methods for diagnosis and characterisation of Chlamydia trachomatis and Neisseria gonorrhoeae: impact on epidemiological surveillance and interventions. APMIS. 2004, 112: 771-784. 10.1111/j.1600-0463.2004.apm11211-1205.x.
Berdichevski T, Barshack I, Bar-Meir S, Ben-Horin S: Pseudomembranes in a patient with flare-up of inflammatory bowel disease (IBD): is it only clostridium difficile or is it still an IBD exacerbation?. Endoscopy. 2010, 42 (Suppl 2): E131-
Satterfield BA, Stewart AF, Lew CS, Pickett DO, Cohen MN, Moore EA, Luedtke PF, O’Neill KL, Robison RA: A quadruplex real-time PCR assay for rapid detection and differentiation of the Clostridium botulinum toxin genes A, B, E and F. J Med Microbiol. 2010, 59: 55-64. 10.1099/jmm.0.012567-0.
Brook I: The role of anaerobic bacteria in cutaneous and soft tissue abscesses and infected cysts. Anaerobe. 2007, 13: 171-177. 10.1016/j.anaerobe.2007.08.004.
DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, Huber T, Dalevi D, Hu P, Andersen GL: Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006, 72: 5069-5072. 10.1128/AEM.03006-05.
Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, Peplies J, Glockner FO: The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013, 41: D590-D596. 10.1093/nar/gks1219.
Wang WC, Lin FM, Chang WC, Lin KY, Huang HD, Lin NS: miRExpress: analyzing high-throughput sequencing data for profiling microRNA expression. BMC Bioinformatics. 2009, 10: 328-10.1186/1471-2105-10-328.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R, Genome Project Data Processing S: The sequence alignment/map format and SAMtools. Bioinformatics. 2009, 25: 2078-2079. 10.1093/bioinformatics/btp352.
Chmielewska A, Szajewska H: Systematic review of randomised controlled trials: probiotics for functional constipation. World J Gastroenterol. 2010, 16: 69-75.
Chassard C, Dapoigny M, Scott KP, Crouzet L, Del’homme C, Marquet P, Martin JC, Pickering G, Ardid D, Eschalier A, et al: Functional dysbiosis within the gut microbiota of patients with constipated-irritable bowel syndrome. Aliment Pharmacol Ther. 2012, 35: 828-838. 10.1111/j.1365-2036.2012.05007.x.
Schwiertz A, Taras D, Schafer K, Beijer S, Bos NA, Donus C, Hardt PD: Microbiota and SCFA in lean and overweight healthy subjects. Obesity (Silver Spring). 2010, 18: 190-195. 10.1038/oby.2009.167.
Balamurugan R, George G, Kabeerdoss J, Hepsiba J, Chandragunasekaran AM, Ramakrishna BS: Quantitative differences in intestinal Faecalibacterium prausnitzii in obese Indian children. Br J Nutr. 2010, 103: 335-338. 10.1017/S0007114509992182.
Million M, Maraninchi M, Henry M, Armougom F, Richet H, Carrieri P, Valero R, Raccah D, Vialettes B, Raoult D: Obesity-associated gut microbiota is enriched in lactobacillus reuteri and depleted in Bifidobacterium animalis and methanobrevibacter smithii. Int J Obes (Lond). 2012, 36: 817-825. 10.1038/ijo.2011.153.
Armougom F, Henry M, Vialettes B, Raccah D, Raoult D: Monitoring bacterial community of human gut microbiota reveals an increase in lactobacillus in obese patients and methanogens in anorexic patients. PloS one. 2009, 4: e7125-10.1371/journal.pone.0007125.
Noor SO, Ridgway K, Scovell L, Kemsley EK, Lund EK, Jamieson C, Johnson IT, Narbad A: Ulcerative colitis and irritable bowel patients exhibit distinct abnormalities of the gut microbiota. BMC Gastroenterol. 2010, 10: 134-10.1186/1471-230X-10-134.
Carroll IM, Ringel-Kulka T, Siddle JP, Ringel Y: Alterations in composition and diversity of the intestinal microbiota in patients with diarrhea-predominant irritable bowel syndrome. Neurogastroenterol Motil. 2012, 24: 521-530. 10.1111/j.1365-2982.2012.01891.x. e248
Duboc H, Rainteau D, Rajca S, Humbert L, Farabos D, Maubert M, Grondin V, Jouet P, Bouhassira D, Seksik P, et al: Increase in fecal primary bile acids and dysbiosis in patients with diarrhea-predominant irritable bowel syndrome. Neurogastroenterol Motil. 2012, 24: 513-520. 10.1111/j.1365-2982.2012.01893.x. e246-517
Kerckhoffs AP, Ben-Amor K, Samsom M, van der Rest ME, de Vogel J, Knol J, Akkermans LM: Molecular analysis of faecal and duodenal samples reveals significantly higher prevalence and numbers of Pseudomonas aeruginosa in irritable bowel syndrome. J Med Microbiol. 2011, 60: 236-245. 10.1099/jmm.0.022848-0.
Saulnier DM, Riehle K, Mistretta TA, Diaz MA, Mandal D, Raza S, Weidler EM, Qin X, Coarfa C, Milosavljevic A, et al: Gastrointestinal microbiome signatures of pediatric patients with irritable bowel syndrome. Gastroenterology. 2011, 141: 1782-1791. 10.1053/j.gastro.2011.06.072.
Rajilic-Stojanovic M, Biagi E, Heilig HG, Kajander K, Kekkonen RA, Tims S, de Vos WM: Global and deep molecular analysis of microbiota signatures in fecal samples from patients with irritable bowel syndrome. Gastroenterology. 2011, 141: 1792-1801. 10.1053/j.gastro.2011.07.043.
Sokol H, Seksik P, Furet JP, Firmesse O, Nion-Larmurier I, Beaugerie L, Cosnes J, Corthier G, Marteau P, Dore J: Low counts of faecalibacterium prausnitzii in colitis microbiota. Inflamm Bowel Dis. 2009, 15: 1183-1189. 10.1002/ibd.20903.
Andoh A, Imaeda H, Aomatsu T, Inatomi O, Bamba S, Sasaki M, Saito Y, Tsujikawa T, Fujiyama Y: Comparison of the fecal microbiota profiles between ulcerative colitis and Crohn’s disease using terminal restriction fragment length polymorphism analysis. J Gastroenterol. 2011, 46: 479-486. 10.1007/s00535-010-0368-4.
Sokol H, Lepage P, Seksik P, Dore J, Marteau P: Temperature gradient gel electrophoresis of fecal 16S rRNA reveals active Escherichia coli in the microbiota of patients with ulcerative colitis. J Clin Microbiol. 2006, 44: 3172-3177. 10.1128/JCM.02600-05.
Wang T, Cai G, Qiu Y, Fei N, Zhang M, Pang X, Jia W, Cai S, Zhao L: Structural segregation of gut microbiota between colorectal cancer patients and healthy volunteers. ISME J. 2012, 6: 320-329. 10.1038/ismej.2011.109.
Moore WE, Moore LH: Intestinal floras of populations that have a high risk of colon cancer. Appl Environ Microbiol. 1995, 61: 3202-3207.
Shen XJ, Rawls JF, Randall T, Burcal L, Mpande CN, Jenkins N, Jovov B, Abdo Z, Sandler RS, Keku TO: Molecular characterization of mucosal adherent bacteria and associations with colorectal adenomas. Gut Microbes. 2010, 1: 138-147. 10.4161/gmic.1.3.12360.
Mcgill R, Tukey JW, Larsen WA: Variations of box plots. Am Stat. 1978, 32: 12-16.
The authors would like to thank the National Science Council of the Republic of China for financially supporting this research under Contract No. NSC 101-2311-B-009-003-MY3 and NSC 100-2627-B-009-002. This work was supported in part by the UST-UCSD International Center of Excellence in Advanced Bio-engineering sponsored by the Taiwan National Science Council I-RiCE Program under Grant Number: NSC 101-2911-I-009-101 and NSC 102-2911-I-009-101, and Veterans General Hospitals and University System of Taiwan (VGHUST) Joint Research Program under Grant Number: VGHUST101-G5-1-1. This work was also partially supported by MOE ATU.
The authors declare no competing interests.
HDH conceived and supervised the study. CMC were responsible for the design, computational analyses, implementation of the system, and drafting the manuscript. FML, THC, CL, TY, TLY, SLW, WCH and WYW were in charge of manuscript revision and data update. All authors read and approved the final manuscript.
Chih-Min Chiu, Feng-Mao Lin contributed equally to this work.