Optimizing microbiome sequencing for small intestinal aspirates: validation of novel techniques through the REIMAGINE study
BMC Microbiology volume 19, Article number: 239 (2019)
The human small intestine plays a central role in the processes of digestion and nutrient absorption. However, characterizations of the human gut microbiome have largely relied on stool samples, and the associated methodologies are ill-suited for the viscosity and low microbial biomass of small intestine samples. As part of the REIMAGINE study to examine the specific roles of the small bowel microbiome in human health and disease, this study aimed to develop and validate methodologies to optimize microbial analysis of the small intestine.
Subjects undergoing esophagogastroduodenoscopy without colon preparation for standard of care were prospectively recruited, and ~ 2 ml samples of luminal fluid were obtained from the duodenum using a custom sterile aspiration catheter. Samples of duodenal aspirates were either untreated (DA-U, N = 127) or pretreated with dithiothreitol (DA-DTT, N = 101), then cultured on MacConkey agar for quantitation of aerobic gram-negative bacteria, typically from the class Gammaproteobacteria, and on blood agar for quantitation of anaerobic microorganisms. DA-DTT exhibited 2.86-fold greater anaerobic bacterial counts compared to DA-U (P = 0.0101), but were not statistically different on MacConkey agar.
DNA isolation from DA-U (N = 112) and DA-DTT (N = 43) samples and library preparation for 16S rRNA gene sequencing were also performed using modified protocols. DA-DTT samples exhibited 3.81-fold higher DNA concentrations (P = 0.0014) and 4.18-fold higher 16S library concentrations (P < 0.0001) then DA-U samples. 16S rRNA gene sequencing revealed increases in the detected relative abundances of obligate and facultative anaerobes in DA-DTT samples, including increases in the genera Clostridium (false discovery rate (FDR) P = 4.38E-6), Enterococcus (FDR P = 2.57E-8), Fusobacterium (FDR P = 0.02) and Bacteroides (FDR P = 5.43E-9). Detected levels of Gram-negative enteropathogens from the phylum Proteobacteria, such as Klebsiella (FDR P = 2.73E-6) and Providencia (FDR P < 0.0001) (family Enterobacteriaceae) and Pseudomonas (family Pseudomonadaceae) (FDR P = 0.04), were also increased in DA-DTT samples.
This study validates novel DTT-based methodology which optimizes microbial culture and 16S rRNA gene sequencing for the study of the small bowel microbiome. The microbial analyses indicate increased isolation of facultative and obligate anaerobes from the mucus layer using these novel techniques.
The Human Microbiome Project  was a groundbreaking introduction to the understanding of the microbiome of the human body. In this effort, many areas of the human microbiome were sampled and characterized, such as the mouth, nose, integument and vaginal tract . The gastrointestinal tract was specifically characterized using stool as an easily available surrogate. Although stool sequencing may reveal the microbial signature of the distal colon, it is well known that stool does not adequately represent the entire gastrointestinal tract, given the multiple environments that exist as one travels from the stomach to the small bowel and then the colon . For example, conditions such as acidity and transit time vary tremendously along the intestine with likely significant effects on microbes by area of exam .
In contrast to the colon, the small intestine, which is divided into the duodenum, jejunum and ileum, is of central importance to digestion and nutrient absorption. Of these, the duodenum has great importance as the site of convergence of chyme from the stomach, enzymes from the pancreas and bile salts from the gall bladder. Clearly, characterizing the microbial populations of the small intestine is of central importance, but efforts to date have been hampered both by the difficulty of obtaining samples, and by challenges associated with adapting sample processing for DNA isolation techniques that were designed for stool.
The REIMAGINE (Revealing the Entire Intestinal Microbiota and its Associations with the Genetic, Immunologic, and Neuroendocrine Ecosystem) study is a large-scale initiative to examine the specific importance of the small bowel microbiome in human health and disease (https://www.cedars-sinai.org/programs/digestive-liver-diseases/clinical/small-bowel-diseases-nutrition/clinical-trials/small-intestinal-sampling-study.html). However, in order to adequately assess the small bowel, techniques need to be assessed and optimized for this location of the intestinal tract. In addition to maintaining sterility and preventing cross-contamination with oral and stomach microbes, particular challenges associated with processing and isolation of DNA from small intestinal samples include the viscosity, small sample volumes, and low microbial biomass. Recently, techniques and adaptations have been described aimed at optimizing DNA sequencing from low-biomass , but not high viscosity samples. Viscosity of the small bowel mucous could impede or affect DNA isolation. One possible remedy would be to reduce viscosity prior to DNA isolation by treating the aspirates with dithiothreitol (DTT), which is commonly used to reduce the disulfide bonds between cysteine residues of proteins and can also reduce the disulfide bonds linking mucin subunits in mucus, improving bacterial recovery and DNA extraction methods . DTT has previously been used to liquefy sputum samples for DNA extraction .
In this study we attempt to develop and validate methodologies to optimize microbial analysis of the small intestine as part of the REIMAGINE study. Optimization included utilization of a new catheter technique for aspiration, steps to improve DNA recovery using DTT, and a new DNA library preparation technique in comparison to conventional DNA isolation and sequencing.
Samples and treatment
A total of 228 subjects had DA samples collected and analyzed as shown in Fig. 1. Of these, 127 DA were not pretreated with DTT prior to microbial culture (DA-U, the untreated group), and 101 were pretreated with DTT prior to microbial culture (the DA-DTT group).
DTT effect on microbial cultures
No growth was observed in MacConkey agar plated with 1x DTT only (negative control). The CFU on MacConkey agar obtained from DA-U subjects ranged from 0 to 240 × 103 CFU/mL (Mean = 10.6 × 103 CFU/mL, Median = 0 CFU/mL, 25th percentile = 0 CFU/mL and 75th percentile = 0 CFU/mL). In the DA-DTT group, the CFU on MacConkey agar ranged from 0 to 1035 × 103 CFU/mL (Mean = 28.22 × 103 CFU/mL, Med = 0 CFU/mL, 25th percentile = 0 CFU/mL and 75th percentile = 2.5 × 103 CFU/mL) (Additional file 1). For the purposes of statistical analysis only, no growth was designated as 1 bacterial CFU/mL of aspirate. DA-DTT exhibited 2.6-fold greater bacterial colonies on MacConkey agar when compared to DA-U, but the p-value did not reach statistical significance (P = 0.14).
No growth was observed on blood agar cultured with 1x DTT only (negative control). The CFU on blood agar obtained from DA-U subjects ranged from 0 to 800 × 103 CFU/mL (Mean = 31.3 × 103 CFU/mL, Median = 0 CFU/mL, 25th percentile = 0 CFU/mL and 75th percentile = 20 × 103 CFU/mL) (Additional file 1). On blood agar, DA-DTT exhibited 2.86-fold greater anaerobic bacterial colonies when compared to DA-U (P = 0.0101). For the purposes of statistical analysis only, no growth was designated as 1 bacterial CFU/mL of aspirate. In the DA-DTT group, CFU on blood agar ranged from 0 to 2070 × 103 CFU/ml (Mean = 89.58 × 103 CFU/mL, Med = 6x103CFU/mL, 25th percentile = 0 CFU/mL and 75th percentile = 98.5 × 103 CFU/mL).
Immediate post aspiration DTT improves DNA extraction and 16S metagenomic library preparation for DA
A total of 155 subjects had their DA samples sequenced and analyzed as shown in Fig. 2 (Additional file 1). The concentrations of DNAs obtained from negative controls (DTT only) were undetectable. The concentrations of DNAs obtained from DA-U subjects ranged from undetectable levels (lower than 10 pg/μL) (n = 18) to 24.6 ng/μL (Med = 0.0908 ng/μL, 25th percentile = 0.02365 ng/μL and 75th percentile = 0.6875 ng/μL) (Additional file 1). Treatment with DTT improved DNA isolation. In the DA-DTT group, DNA concentrations ranged from undetectable levels (n = 3) to 68.8 ng/μL (Med = 0.346 ng/μL, 25th percentile = 0.0906 ng/μL and 75th percentile = 1.91 ng/μL) and were 3.81-fold higher than those from DA-U (Mann Whitney P = 0.0014).
The concentrations of the final 16S libraries amplified from negative controls were undetectable. The concentrations of the final 16S libraries amplified from DA-U samples (i.e. those for which DTT was added only for the removal of the All Protect reagent) ranged from 0.14 ng/μL to 136 ng/μL (median = 21.8 ng/μL, 25th percentile = 5.1 ng/μL and 75th percentile = 69.8 ng/μL) and correlated with the initial DNA concentrations (Spearman r = 0.316, P = 0.001) (Additional file 1).
The concentrations of the final 16S libraries amplified from DA-DTT samples (i.e. those for which DTT was added both before microbial culture and for removal of the All Protect reagent) were 4.18x higher than those of libraries from DA-U samples (P < 0.0001) (see Fig. 3). The library concentrations ranged from 1.69 ng/μL to 302 ng/μL (Med = 91.2 ng/μL, 25th percentile = 36.6 ng/μL and 75th percentile = 117 ng/μL) and correlated with the initial DNA concentrations (Spearman r = 0.443, P = 0.003) (Additional file 1).
All samples had at least 9000 sequences and no exclusions were performed. A total of 112 DA-U and 43 DA-DTT samples were sequenced. The difference in average library sizes between the groups was less than 2-fold (Additional file 2). Predictions for significant differentially abundant Operational Taxonomic Units (OTUs) were performed following recommendations from McMurdie and Holmes , and from Weiss et al. , used when the average library size for each group is approximately equal and/or the fold difference between groups is not high (> 2-3x on average) (Additional file 1).
Considering observations regarding contamination of DNA extraction kits with traces of bacterial DNA ,16S sequencing was also performed on negative control samples (DTT only). Less than 0.03% of the total sequences generated in each MiSeq run was assigned to negative control samples, 4433 sequences on average. 27.63% of the sequences assigned to negative controls were identified as bacterial DNA, mostly belong to the Pseudomonas genus OTU 646549 (63.5%), and Bacteroides genus OTUs 1,749,079, 193,591 and 359,538 (12%).
All OTUs observed in negative controls were also detected in both groups analyzed, DA-U and DA-DTT. The OTUs assigned to Bacteroides genus observed in negative controls represented less than 3% of all OTUs assigned to this same genus in DA-U and DA-DTT, thus these OTUs were not excluded during downstream analysis. The OTU assigned to Pseudomonas genus (646549) observed in negative controls represented 67% of the OTUs assigned to this same genus in DA-U and DA-DTT, and considering the high risk of bias during analysis the OTU 646549 was excluded during comparisons between DA-U and DA-DTT groups (Additional file 1).
DTT increases the detected relative abundance of anaerobic bacteria in DA
The main two dominant phyla observed in DA-DTT and DA-U were Firmicutes and Proteobacteria, followed by smaller proportions of Actinobacteria, Fusobacteria, Bacteroidetes and TM7 (Fig. 4, Table 1). After pretreatment with DTT, DA showed increased relative abundance of the phyla Proteobacteria (FC = 6.22, FDR P = 7.71E-7), Bacteroidetes (FC = 2.19, FDR P = 0.03) and Fusobacteria (FC = 1.96, FDR P = 0.03), when compared to DA-U (Table 1). There were also smaller changes in the relative abundances of Actinobacteria and TM7 that did not reach significance (Table 1).
Although no changes were seen in the detected relative abundances of Clostridia and Bacilli, the two main classes from the phylum Firmicutes, in DA-DTT vs. DA-U, specific increases were observed in families from both of these classes. Specifically, DA-DTT exhibited increased detected relative abundances of the family Clostridiaceae (FC = 5.10, FDR P < 0.0001) and genus Clostridium (FC = 4.06, FDR P = 4.38E-6), which are Gram-positive obligate anaerobes, and of the family Enterococcaceae (FC = 76.22, FDR P = 2.62E-11) and genus Enterococcus (FC = 42.18, FDR P = 2.57E-8), which are Gram-positive facultative anaerobic lactic acid bacteria (Table 2, Fig. 4).
The detected relative abundances of several obligate anaerobic bacteria were increased in DA-DTT vs. DA-U, including Fusobacterium (phylum Fusobacteria), which are Gram-negative bacilli (FC = 2.29, FDR P = 0.02), and Bacteroides (phylum Bacteroidetes) (FC = 28.08, FDR P = 5.43E-9) (Table 2, Fig. 4).
Pretreatment with DTT increased the detected relative abundance of gram-negative enteropathogens from the phylum Proteobacteria
The relative abundance of the phylum Proteobacteria, a major phylum of Gram-negative bacteria, detected in DA-DTT was increased compared to that detected in DA-U (FC = 6.22, FDR P = 7.71E-7). The detected relative abundances of three of the five most important classes from this phylum were significantly increased in DA-DTT compared to DA-U - class Gammaproteobacteria, which comprises several enteropathogens (FC = 8.44, FDR P = 4.25E-8) , class Alphaproteobacteria, which includes mainly phototrophic bacteria (FC = 7.94, FDR P = 2.60E-8), and class Deltaproteobacteria, which includes sulfate- and sulfur-reducing bacteria (FC = 6.35, FDR P = 9.7E-5) (Table 3). Smaller changes in classes Betaproteobacteria and Epsilonproteobacteria did not reach significance (Table 3).
The increase in the detected relative abundance of Gammaproteobacteria in DA-DTT was partially driven by higher detected relative abundances of Enterobacteriaceae family members (FC = 5.46, FDR P = 1.47E-3), including important enteropathogens and pathogens that cause infection in several parts of the human body, such as Klebsiella and Providencia (see Table 4) . The detected relative abundances of other members of the class Gammaproteobacteria were also increased in DA-DTT vs. DA-U, including the family Aeromonadaceae (FC = 63.61, FDR P = 1.18E-13) and the genus Pseudomonas (family Pseudomonadaceae) (FC = 2.65, FDR P = 0.04).
The increase in detected relative abundances of Alphaproteobacteria and Deltaproteobacteria in DA-DTT (Table 3) was driven by increases in detection of the order Rhizobiales (FC = 19.03, FDR P = 4.33E-13), and of sulfur-producing bacteria from the orders Desulfobacterales (FC = 42.61, FDR P < 0.0001) and Desulfovibrionales (FC = 6.41, FDR P = 6.71E-3), respectively.
Pretreatment with DTT does not affect microbial diversity in DA
Sample rarefaction curves showed a similar pattern, which verified that most of the species present in each sample from DA-DTT and DA-U groups were observed (Fig. 5) . DA-DTT exhibited the same alpha diversity as DA-U, as demonstrated by Simpson’s index (P = 0.9287) and Shannon entropy (P = 0.8066) (Fig. 6). Beta diversity of the DA-U and DA-DTT microbiome was analyzed based on the weighted UniFrac metric. Principal Coordinate Analysis plot showed no clustering of the DA-DTT (n = 43) and DA-U groups (n = 112) (see Fig. 7).
In this paper, we develop and validate a novel methodological approach based on the use of the reducing agent dithiothreitol (DTT) that resolves issues related to low microbial biomass from luminal duodenal aspirates (DA). The use of DTT clearly increases the number of bacteria detected on culture plates, and also increases DNA yields and the concentration of V3/V4 libraries for sequencing, which in turn results in important differences in the microbial populations detected in DA.
Given the central role of the small intestine in the processes of digestion and nutrient absorption, accurate characterization of the human small intestinal microbiome is an important future consideration. The small intestine is not as heavily colonized as the large intestine, ranging under healthy conditions from 103 to 104 bacteria per mL of intestinal content in the duodenum and jejunum to 108 bacteria per mL in the ileum, compared to 1011 bacteria per gram of wet stool in the colon [12, 13]. In addition to having low bacterial biomass, duodenal luminal contents are viscous due to the mucus layer present in the small intestine [14, 15], and require special handling during sample collection and processing prior to culture and DNA extraction, in order to increase the likelihood of the assessment of all microbial communities, including those associated with the mucus layer.
In the assessment of the microbiome, it is essential to accurately and completely assess the microbial components in a sample. Although standards have been set for the assessment of the stool microbiome, these standards have not been assessed for small intestinal fluid assessment. Mucous in general is a viscous fluid that can trap bacteria in its matrix and previous studies performed with sputum samples have shown that treating this viscosity has an impact on the microbial assessment [5, 16]. However, until now, no studies have investigated the impact on microbial assessment and DNA recovery in aspirates collected from small bowel. There are agents known to safely and effectively improve microbial assessment and DNA yield in viscous samples, and one such agent is DTT, which can reduce the disulfide bonds between mucin subunits.
In this study we established a methodology to improve microbial DNA recovery from small bowel aspirates, which includes different sample processing steps when compared to conventionally published methods for extracting DNA for microbiome assessment of gut materials such as stool [17,18,19]. The concentrations of DNAs extracted from DA ranged from very low levels to up to 70 ng/ml when samples were pretreated with DTT, more than 3-fold higher than those from samples which were not pretreated with the reducing agent. Initial DNA concentrations exhibited a higher positive correlation with those of the final V3/V4 libraries for DTT-pretreated DA compared to non-pretreated DA, which may indicate a specific increase in the isolation of bacterial DNAs. The use of a fixed initial DNA concentration during the preparation of sequencing libraries from DA should be carefully analyzed, and the addition of DTT during sample processing and for removal of the All Protect reagent is highly recommended as this increases the initial yield of microbial DNA prior sequencing library preparation.
In addition to increases in DNA yields and library concentrations, DA processed with DTT prior to microbial culture exhibited higher numbers of bacterial colonies on blood agar plates incubated under anaerobic conditions. The relative abundance of specific obligate and facultative anaerobes detected in DA-DTT was increased compared to DA-U samples, which may reflect the increased detection of microbes associated with the mucus layer. For example, the relative abundance of the genus Enterococcus (phylum Firmicutes) detected in DA-DTT was significantly increased. This genus comprises over 50 Gram-positive facultative anaerobic lactic acid cocci species isolated from numerous environments, including the human GI tract . Enterococcus species constitute up to 1% of the gut microbiota, and most species can grow on blood agar plates under anaerobic conditions. The relative abundance of the genus Clostridium (phylum Firmicutes), which is comprised of obligate anaerobes and some aerotolerant species, detected in DA-DTT samples was also increased.
The relative abundance of the genus Bacteroides, comprised of Gram-negative obligately anaerobic bacilli, detected in DA-DTT was also significantly increased compared to DA-U samples. Species from this genus can grow on blood agar and are well-adapted to the gastrointestinal tracts of mammals, including humans [21, 22]. The human large intestine is densely colonized with species from the genus Bacteroides (phylum Bacteroidetes) , many of which perform essential metabolic functions for the host, including the metabolism of proteins and complex sugars. In contrast, the small intestine is not as heavily colonized by members of the phylum Bacteroidetes, which comprised less than 4% of the total microbes detected in DA-U. With the addition of the reducing agent DTT, which breaks the disulfide bonds linking mucin subunits in mucus prior to microbial culture and DNA extraction, the relative abundance of Bacteroidetes detected in DA increased significantly from 4 to 7%, indicating a possible role for species from this phylum in mucus metabolism. Mucus is a dynamic matrix, consisting of mucin glycoproteins secreted by intestinal goblet cells, which lubricates the transit of intestinal contents, amongst other functions. Mucus glycoproteins can be used as a carbon source by many asaccharolytic microorganisms, and the low oxygen levels at atmospheric pressures allow the colonization and growth of anaerobes in mucus [24, 25].
The phylum Proteobacteria also includes aerotolerant asaccharolytic microorganisms that require proteinaceous substrates as carbon and energy sources, such as Campylobacter , as well as facultative anaerobes from the family Enterobacteriaceae included in the “Mucosally Associated Consortium” in the colon described by Albenberg et al. . Pretreatment of DA with DTT increases the detected relative abundance of many Enterobacteriaceae members, including the clinically important genera Klebsiella, Providencia and Salmonella as well as unknown members. Providencia and Salmonella include motile species that can adhere to mucus and epithelial cells and actively invade the host epithelium [27,28,29]. The relative abundance of the genus Pseudomonas, detected in DA-DTT was also increased compared to non-pretreated DA. Members of this genus, including the most studied species P. aeruginosa, are also motile and can be part of the normal human microflora, but are also important clinically as they are known to cause hospital-acquired infections such as pneumonia and urinary tract infections .
A limitation of this study was that DTT treatment vs. non-treatment could not be tested in different portions of the same samples. This was partly due to the small sizes of the individual samples, and partly due to the fact that their viscosity made it impossible to divide them evenly. We have attempted to compensate for this through the number of samples tested.
Surprisingly, the changes in several microbial taxa in DA samples after the addition of DTT did not affect the overall microbial diversity. These findings further suggest that the addition of the reducing agent DTT improves microbial assessment and DNA recovery without causing a dramatic change in the microbial balance in the aspirate samples.
This study validates methodology to optimize yield for culture, and for DNA extraction for analysis of the small bowel microbiome. Culture totals, microbial DNA and microbiome analysis demonstrate marked differences with this new technique. This suggests that conventional techniques for DNA isolation provide an incomplete picture of the microbial environment in the small bowel. Thus, this new technique appears ideal for small bowel microbiome assessment.
The REIMAGINE study is a large-scale study designed to examine the relationship between the small bowel microbiome in human health and disease. In brief, the study involves collecting data in consecutive patients undergoing routine upper gastrointestinal endoscopy (esophagogastroduodenoscopy). After a battery of questionnaires for comprehensive collection of health information, and obtaining serum and genetic samples, the endoscopy entailed the collection of small bowel aspirates and two mucosal biopsies. The first phase of this study was to develop, validate and optimize small intestinal aspiration techniques, and microbial sample acquisition.
Male and female subjects aged 18–85 undergoing esophagogastroduodenoscopy (EGD) without colon preparation for standard of care purposes were prospectively recruited for this study. Potential participants were identified by study staff and their eligibility was verified by co-investigators or the PI. Although there are no exclusion criteria for the REIMAGINE study, small bowel biopsies are not collected from subjects with bleeding disorders or advanced cirrhosis of the liver with coagulopathy and intestinal varices where the international normalized ratio (INR) was greater than 1.5, in order to minimize the risk of bleeding from the biopsy site. The REIMAGINE study protocol was approved by the Institutional Review Board at Cedars-Sinai Medical Center, and all subjects provided informed written consent prior to participating in the study. Samples from subjects taking antibiotics were not included in the present study.
Small intestinal sample procurement
During EGD, samples of luminal fluid (~ 2 ml) were obtained from the duodenum using a custom sterile aspiration catheter (Hobbs Medical, Inc.). The custom catheter consisted of a newly designed double lumen sterile catheter, with the inner lumen maintaining sterility during insertion by applying sterile bone wax into the open tip of the external catheter. During endoscopy, the endoscopist is instructed to immediately enter the duodenum and insert the aspiration catheter. The inner catheter then dislodges the bone wax, exposing the sterile inner catheter. This inner catheter is used to aspirate duodenal fluid through lasered side holes to acquire a volume of 2 mL. These precautions eliminated the risk of oral and gastric contamination.
Immediately after aspiration, samples from all duodenal aspirates were cultured on MacConkey agar (Becton Dickinson, Franklin Lakes, NJ, EUA) and blood agar (Becton Dickinson) for determination of the number of colony-forming units (CFU) per mL of aspirate. To assess the effect of viscosity on culture, a subset of aspirates were not pretreated and simply cultured (the DA-U group for “untreated”) and another subset were first pretreated with the reducing agent Dithiothreitol (DTT) (Sputolysin® Reagent, Cat. 560,000-1SET, EMD Millipore Corp, Billerica, MA, USA) (the DA-DTT group) (see Fig. 1).
For DA-DTT samples, 1x DTT (6.5 mM dithiothreitol in 100 mM phosphate buffer, pH 7.0) was added to an aliquot of the DA in a 1:1 ratio and the resulting mixture was vortexed until the sample was liquified (typically 30 s). 100 μL of the liquified mixture was serially diluted with 900 μL sterile 1x PBS and samples of the 1:10 and 1:100 dilutions were plated in duplicate on MacConkey agar under aerobic conditions for the quantitation of Gram-negative bacilli, and on blood agar under anaerobic conditions for the quantitation of total anaerobes. For DA-U samples, 100 μL of DA was diluted directly with 900 μL sterile 1x PBS and samples of the 1:10 and 1:100 dilutions were plated as described above. All plates were incubated at 37 °C for 16–18 h, after which colonies were electronically counted using the Scan 500 (Interscience, Paris, France). As a negative control, 100 μL of 1x DTT was also cultured aerobically on MacConkey agar and anaerobically on blood agar. All dilution factors were taken into account for final determination of microbial burden.
After aliquots for microbial culture were taken, remaining DA-U and DA-DTT samples were centrifuged at high speed (17,136 x g) for 10 min and the supernatant was carefully removed and stored at − 80 °C for future metabolomic analyses. 500 μL of sterile All Protect reagent (Qiagen, Hilden, Germany) was added to each pellet for stabilization of DNA, RNA and proteins, and the pellets were stored at − 80 °C prior to DNA isolation and analysis of the DA microbiome.
DNA extraction and quantification from aspirates
DA-U and DA-DTT samples were thawed on ice and 1x DTT was added in a 1:1 ratio, after which the samples were vortexed until the All Protect reagent was fully liquefied (around 30 s). DNA extraction was then performed for both groups using the MagAttract PowerSoil DNA KF Kit (Qiagen, cat. No. 27000–4-KF) with some modifications. DNA extraction was also performed on negative control samples (1x DTT) as a control.
The lysis step was carried out by adding garnet beads (Qiagen, cat. No. 13123–50) and 746 μL PowerBead Solution to each pellet-containing tube, followed by 4 μL RNase A and 60 μL SL Solution (Lysis buffer) in this specific order. Tubes were sealed with parafilm, vortexed horizontally for 15 min, and then centrifuged for 6 min at 4500 x g. The supernatants were transferred to new tubes containing 450 μL IR Solution, vortexed for 3 s, incubated at 4 °C for 10 min, and then centrifuged for 6 min at 4500 x g. The supernatants were transferred to new tubes and centrifuged for a further 6 min at 4500 x g. 450 μL of the resulting supernatants were added to deep 96-well KingFisher plates containing magnetic beads and DNA extraction was performed using the KingFisher Duo (Thermo Fisher Scientific, Waltham, MA, USA). The final DNA volume was 100 μL. DNAs were then quantified using Qubit ds DNA BR Assay kits (Invitrogen by Thermo Fisher Scientific, Waltham, MA, USA) on a Qubit 4 Fluorometer (Invitrogen).
Library preparation and 16S rRNA gene sequencing
16S library preparation for DNAs from all groups was performed according to the Illumina (Illumina, San Diego, CA, USA) protocol https://support.illumina.com/documents/documentation/chemistry_documentation/16s/16s-metagenomic-library-prep-guide-15044223-b.pdf, with some modifications. The V3 and V4 regions were amplified using the gene-specific primers S-D-Bact-0341-b-S-17 and S-D-Bact-0785-a-A-21 published and validated by Klindworth et al. . The primers were modified in accordance with the protocol by adding the Illumina sequencing adapters to each one.
The full-length primer sequences used were:
16S amplicon PCR forward primer: 5’TCGTCGGCAGCGTCAGATGTGTATAAGAGACAGCCTACGGGNGGCWGCAG
16S amplicon PCR reverse primer: 5’GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAGGACTACVHGGGTATCTAATCC
The 16S library preparation protocol was modified as follows: 5 μL of DNA was added to a Master Mix (0.5 μL of 10 μM 16S Amplicon PCR Forward primer, 0.5 μL of 10 μM 16S Amplicon PCR Reverse primer, 12.5 μL 2x KAPA HiFi HotStart ReadyMix and 6.5 μL of molecular grade PCR H2O) and the PCR was performed as follows:
Initial denaturation step at 95 °C for 3 min
27 cycles of: 95 °C for 30 s, 55 °C for 30 s and 72 °C for 30 s
72 °C for 5 min
Hold at 4 °C
An optimized Clean-Up step was performed with Agencourt AMPure XP beads using the modifications proposed by Quail et al. . After adding the beads to each Amplicon PCR 96-well plate on the magnetic stand, samples were incubated for five minutes followed by two wash steps with 80% ethanol. The beads were air dried for five minutes. After removing the plate from magnetic stand, beads were incubated with EB Buffer (Qiagen) for five minutes to elute the DNA. The plate was placed back on the magnetic stand and after 2–3 min the supernatant was transferred to an empty clean well, preventing the transfer of the beads with the supernatant.
Five μL of the final Amplicon PCR product was used for the Index PCR, which was performed using the Nextera XT Index kit and 2x KAPA HiFi HotStart ReadyMix, following the Illumina protocol for 8 cycles. After a second modified Clean-Up step, the final product was quantified using Qubit ds DNA BR Assay kits and Qubit 1X dsDNA HS Assay kits on a Qubit 4 Fluorometer and analyzed using Agilent DNA 1000 chips (Agilent Technologies, Santa Clara, CA) and Agilent HS DNA chips (Agilent) on an Agilent 2100 Bioanalyzer System.
16S rRNA gene sequencing and analysis
The V3 and V4 libraries prepared using DNAs from DA-DTT and DA-U groups were sequenced using a MiSeq Reagent Kit v3 (600-cycles) on a MiSeq System (Illumina, San Diego, California). 2 × 301 cycles paired-end sequencing was performed according to manufacturer’s protocol and 5% Phix (Illumina) was added to each library pool.
Operational Taxonomic Unit (OTU) clustering and taxonomic analyses were performed using CLC Genomics Workbench v. 10.1.1 and CLC Microbial Genomics Module v. 2.5 (Qiagen). Sequences were first trimmed to remove 13 bases at the 5′ terminal position and merged considering the alignment scores as follows: mismatch cost of 2, gap cost of 2, zero maximum unaligned end mismatches and minimum score of 30. After merging, sequences were clustered into OTUs at 97% sequence similarity level using the Amplicon-Based OTU clustering tool. The creation of new OTUs was allowed considering 97% taxonomic similarity. The most abundant sequences were selected as representative of each cluster, and then assigned to a taxonomy level using CLC Microbial Genomics default values and the Greengenes Database 2013 release. Alpha diversity indexes (Chao1, Simpson and Shannon) were calculated using the Abundance Analysis tool. The weighted Unifrac metric was used to calculate inter-sample diversity (beta diversity).
Multiple comparisons and statistical analyses were performed using CLC Genomics Workbench v. 10.1.1 and CLC Microbial Genomics Module v. 2.5 (Qiagen). A Negative Binomial Generalized Linear Model (GLM) model was used to obtain maximum likelihood estimates for an OTU’s log-fold change between two conditions, and the Wald test was used to determine significance, as part of the CLC package available at https://www.qiagenbioinformatics.com/products/clc-genomics-workbench/. False Discovery Rate (FDR) was performed to correct P-values. Fold changes are calculated from the GLM, which corrects for differences in library size between the samples and the effects of confounding factors. Again, these calculations were performed using the CLC package. It is therefore not possible to derive these fold changes from the original counts by simple algebraic calculations. Two-tailed Spearman r correlations, Mann-Whitney tests and graph construction were performed using GraphPad Prism 7.02 (GraphPad Software, La Jolla, CA, USA). For statistical analysis purposes only, no growth on blood agar and MacConkey agar (CFU/ml = 0) was assigned as 1 CFU/ml.
Availability of data and materials
The datasets generated during the current study are available in the National Center for Biotechnology Information (NCBI) BioProject Repository https://www.ncbi.nlm.nih.gov/bioproject under BioProject ID PRJNA520899.
Duodenal aspirates pretreated with dithiothreitol
Duodenal aspirates untreated
False discovery rate
Generalized linear model
International normalized ratio
Operational taxonomic unit
Polymerase chain reaction
- REIMAGINE study:
Revealing the Entire Intestinal Microbiota and its Associations with the Genetic, Immunologic, and Neuroendocrine Ecosystem study
ribosomal ribonucleic acid
Turnbaugh PJ, Ley RE, Hamady M, Fraser-Liggett CM, Knight R, Gordon JI. The human microbiome project. Nature. 2007;449:804–10.
Costello EK, Lauber CL, Hamady M, Fierer N, Gordon JI, Knight R. Bacterial community variation in human body habitats across space and time. Science. 2009;326:1694–7.
Hillman ET, Lu H, Yao T, Nakatsu CH. Microbial ecology along the gastrointestinal tract. Microbes Environ. 2017;32:300–13.
Minich JJ, Zhu Q, Janssen S, Hendrickson R, Amir A, Vetter R, et al. KatharoSeq Enables High-Throughput Microbiome Analysis from Low-Biomass Samples. mSystems. 2018;3.
Oriano M, Terranova L, Teri A, Sottotetti S, Ruggiero L, Tafuro C, et al. Comparison of different conditions for DNA extraction in sputum - a pilot study. Multidiscip Respir Med. 2019;14:6.
McMurdie PJ, Holmes S. Waste not, want not: why rarefying microbiome data is inadmissible. PLoS Comput Biol. 2014;10:e1003531.
Weiss S, Xu ZZ, Peddada S, Amir A, Bittinger K, Gonzalez A, et al. Normalization and microbial differential abundance strategies depend upon data characteristics. Microbiome. 2017;5:27.
Glassing A, Dowd SE, Galandiuk S, Davis B, Chiodini RJ. Inherent bacterial DNA contamination of extraction and sequencing reagents may affect interpretation of microbiota in low bacterial biomass samples. Gut Pathog. 2016;8:24.
Williams KP, Gillespie JJ, Sobral BW, Nordberg EK, Snyder EE, Shallom JM, et al. Phylogeny of gammaproteobacteria. J Bacteriol. 2010;192:2305–14.
Kolling G, Wu M, Guerrant RL. Enteric pathogens through life stages. Front Cell Infect Microbiol. 2012;2:114.
Cox KD, Black MJ, Filip N, Miller MR, Mohns K, Mortimor J, et al. Community assessment techniques and the implications for rarefaction and extrapolation with hill numbers. Ecol Evol. 2017;7:11213–26.
Leser TD, Molbak L. Better living through microbial action: the benefits of the mammalian gastrointestinal microbiota on the host. Environ Microbiol. 2009;11:2194–206.
Sender R, Fuchs S, Milo R. Revised estimates for the number of human and Bacteria cells in the body. PLoS Biol. 2016;14:e1002533.
Johansson ME, Sjovall H, Hansson GC. The gastrointestinal mucus system in health and disease. Nat Rev Gastroenterol Hepatol. 2013;10:352–61.
Ermund A, Schutte A, Johansson ME, Gustafsson JK, Hansson GC. Studies of mucus in mouse stomach, small intestine, and colon. I. Gastrointestinal mucus layers have different properties depending on location as well as over the Peyer's patches. Am J Physiol Gastrointest Liver Physiol. 2013;305:G341–7.
Terranova L, Oriano M, Teri A, Ruggiero L, Tafuro C, Marchisio P, et al. How to process sputum samples and extract bacterial DNA for microbiota analysis. Int J Mol Sci. 2018;19.
Jovel J, Patterson J, Wang W, Hotte N, O'Keefe S, Mitchel T, et al. Characterization of the gut microbiome using 16S or shotgun Metagenomics. Front Microbiol. 2016;7:459.
Wang Z, Zolnik CP, Qiu Y, Usyk M, Wang T, Strickler HD, et al. Comparison of fecal collection methods for microbiome and metabolomics studies. Front Cell Infect Microbiol. 2018;8:301.
Jones MB, Highlander SK, Anderson EL, Li W, Dayrit M, Klitgord N, et al. Library preparation methodology can influence genomic and functional predictions in human microbiome research. Proc Natl Acad Sci U S A. 2015;112:14024–9.
Dubin K, Pamer EG. Enterococci and their interactions with the intestinal microbiome. Microbiol Spectr. 2014;5.
Ley RE, Hamady M, Lozupone C, Turnbaugh PJ, Ramey RR, Bircher JS, et al. Evolution of mammals and their gut microbes. Science. 2008;320:1647–51.
Eckburg PB, Bik EM, Bernstein CN, Purdom E, Dethlefsen L, Sargent M, et al. Diversity of the human intestinal microbial flora. Science. 2005;308:1635–8.
Rajilic-Stojanovic M, de Vos WM. The first 1000 cultured species of the human gastrointestinal microbiota. FEMS Microbiol Rev. 2014;38:996–1047.
Albenberg L, Esipova TV, Judge CP, Bittinger K, Chen J, Laughlin A, et al. Correlation between intraluminal oxygen gradient and radial partitioning of intestinal microbiota. Gastroenterology. 2014;147:1055–63 e1058.
Li H, Limenitakis JP, Fuhrer T, Geuking MB, Lawson MA, Wyss M, et al. The outer mucus layer hosts a distinct intestinal microbial niche. Nat Commun. 2015;6:8292.
van der Stel AX, van de Lest CHA, Huynh S, Parker CT, van Putten JPM, Wosten M. Catabolite repression in campylobacter jejuni correlates with intracellular succinate levels. Environ Microbiol. 2018;20:1374–88.
O'Hara CM, Brenner FW, Miller JM. Classification, identification, and clinical significance of Proteus, Providencia, and Morganella. Clin Microbiol Rev. 2000;13:534–46.
Kurmasheva N, Vorobiev V, Sharipova M, Efremova T, Mardanova A. The potential virulence factors of Providencia stuartii: motility, adherence, and invasion. Biomed Res Int. 2018;2018:3589135.
Elhadad D, Desai P, Rahav G, McClelland M, Gal-Mor O. Flagellin is required for host cell invasion and Normal Salmonella Pathogenicity Island 1 expression by Salmonella enterica Serovar Paratyphi a. Infect Immun. 2015;83:3355–68.
Lister PD, Wolter DJ, Hanson ND. Antibacterial-resistant Pseudomonas aeruginosa: clinical impact and complex regulation of chromosomally encoded resistance mechanisms. Clin Microbiol Rev. 2009;22:582–610.
Klindworth A, Pruesse E, Schweer T, Peplies J, Quast C, Horn M, et al. Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. 2013;41:e1.
Quail MA, Swerdlow H, Turner DJ. Improved protocols for the illumina genome analyzer sequencing system. Curr Protoc Hum Genet. 2009; Chapter 18: Unit 18 12.
This study was supported in part by a grant from the Monica Lester Charitable Trust to Dr. Ruchi Mathur. This funding body played no role in the design of the study, in the collection, analysis, and interpretation of data, or the writing of the manuscript.
Ethics approval and consent to participate
The study protocol was approved by the Institutional Review Board at Cedars-Sinai Medical Center (Pro00035192), and all subjects provided informed written consent prior to participating in the study.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1. Title of data: Culture data, sequenced samples information and OTU table. Description of the data: the culture data section (sheet 1) shows the number of bacterial colonies observed on MacConkey agar (incubated aerobically) and blood agar (incubated anaerobically) for all DA-DTT (N = 101) and DA-U (N = 127) samples. Data are presented as CFU/mL of duodenal aspirate. The sequenced samples information section (sheet 2) provides data for the sequenced DA-DTT (N = 43) and DA-U (N = 112) samples, including the number of bacterial colonies observed on MacConkey agar and blood agar (in CFU/mL of duodenal aspirate), initial DNA concentration, library concentration and library size (number of sequences). The OTU table section (sheet 3) contains the OTU table generated from the sequenced DA-DTT (N = 43) and DA-U (N = 112) samples. The table provides the OTU numbers, taxonomic classification, number of reads for each OTU and the OTU sequence.
Additional file 2. Title of data: V3-V4 16S rRNA library size of DA-DTT and DA-U samples. Description of the data: library sizes (number of sequences in the OTU table) for DA-DTT and DA-U samples. The table provides the library size mean, standard deviation, standard error of the mean, median, 25 and 75% percentiles.
About this article
Cite this article
Leite, G.G.S., Morales, W., Weitsman, S. et al. Optimizing microbiome sequencing for small intestinal aspirates: validation of novel techniques through the REIMAGINE study. BMC Microbiol 19, 239 (2019). https://doi.org/10.1186/s12866-019-1617-1
- Small intestine
- Mucus layer
- Methodology optimization
- Microbial culture
- 16S rRNA gene sequencing