This protocol provides a simple, cost-effective DNA isolation method for the analysis of murine gut microbiota alterations during the development of autoimmune disease.
Gut microbiota has an important role in educating the immune system. This relationship is extremely important for understanding autoimmune diseases that are not only driven by genetic factors, but also environmental factors that can trigger the onset and/or worsen the disease course. A previously published study on the dynamics of the gut microbiota in lupus-prone MRL/lpr female mice showed how changes of the gut microbiota can alter disease progression. Here, a protocol is described for extracting representative samples from the gut microbiota for studies of autoimmunity. Microbiota samples are collected from the anus and processed, from which the DNA is extracted using a phenol-chloroform method and purified by alcohol precipitation. After PCR is performed, purified amplicons are sequenced using a Next Generation Sequencing platform at Argonne National Laboratory. Finally, the 16S ribosomal RNA gene sequencing data is analyzed. As an example, data obtained from gut microbiota comparisons of MRL/lpr mice with or without CX3CR1 are shown. Results showed significant differences in genera containing pathogenic bacteria such as those in the phylum Proteobacteria, as well as the genus Bifidobacterium, which is considered part of the healthy commensal microbiota. In summary, this simple, cost-effective DNA isolation method is reliable and can help the investigation of gut microbiota changes associated with autoimmune diseases.
Humans and bacteria have coexisted for a long time. They have established a codependent relationship with mutual beneficial effects that influences host immune responses in quantitative and qualitative ways1. Recent studies suggest an association between the gut microbiota composition and the pathogenesis of autoimmune diseases that include multiple sclerosis2, rheumatoid arthritis3, type 2 diabetes4, Inflammatory bowel disease5, and systemic lupus erythematosus (SLE)6. However, whether the gut microbiota is the main cause or a secondary effect of these autoimmune diseases is still unclear7. Potentially, gut microbiota could exacerbate the disease during the effector phase of autoimmune disorders or play a role in regulating the induction of these diseases8.
Intestinal dysbiosis has been reported in female lupus-prone MRL/Mp-Faslpr (MRL/lpr) mice, and gut microbiota changes with a significant depletion of Lactobacilli were observed9. When a mixture of five Lactobacillus strains was orally administered, lupus-like symptoms were largely attenuated in these mice, suggesting an essential role of microbiota in regulating SLE pathogenesis.
The following technique of DNA extraction allows to follow microbiota fluctuations and analyze them qualitatively and quantitively during the course of murine SLE-like disease in lupus-prone mice. Whether to examine the healthy gut microbiota or define dysbiosis, it is important to critically examine how data is collected and whether it is accurate and reproducible10. Every step is critical in this process. An appropriate methodology must be used to extract microbial DNA as any possible problem introducing biases during the DNA extraction process could result in inaccurate microbial representation. While the phenol-chloroform method is described here, there are commercially available kits to extract DNA from bacteria that work well in particular cases11. However, their usability is limited by the cost and required sample quantity.
The protocol presented here is cost-effective and requires only a small amount of sample. It works fine with any kind of stool sample and is useful in studying the dynamics of the gut microbiota over time as well as comparisons among groups. DNA is isolated with a method of alcohol purification, which uses phenol, chloroform, and isoamyl alcohol. Alcohol-based extraction helps to clean and remove the sample of proteins and lipids, where DNA is precipitated in the final step. The proposed method has significantly high efficiency and quality and has been proven to be accurate in identifying bacterial populations. One critical note during the procedure is that DNA contamination can occur, and thus appropriate sample handling is required12.
The DNA is then analyzed by the Next Generation Sequencing platforms for the 16S rRNA gene, such as the Illumina MiSeq. In particular, the V4 hyper-variable region is analyzed to provide a better quantification for high-rank taxa13. The subsequent bioinformatics analysis is outsourced, followed by in-house analysis using standard statistical methods. There are numerous open-source bioinformatics software programs available for downstream sequencing, and the type of analyses performed depends heavily on the specific biological question of interest14. This protocol focuses specifically on the experimental steps prior to sequencing and provides a more versatile, cost-effective, comparable, and efficient method to obtain DNA from fecal samples.
The Cx3cr1gfp/gfp locus of B6.129P2(Cg)-Cx3cr1tm1Litt/J mice was backcrossed to MRL/MpJ-Faslpr/J (MRL/lpr) for 10 generations to generate MRL/lpr-CX3CR1gfp/gfp mice. Genome screening using single nucleotide polymorphism (SNP) panels confirmed that the genetic background of newly generated mice was more than 97% identical to that of MRL/lpr. After that, mice were bred and maintained in a specific pathogen-free environment following the specific requirements of the Institutional Animal Care and Use Committee (IACUC) at Virginia Tech. Samples were collected when the mice were 13, 14, and 15 weeks old.
1. Collection of microbiota samples from mice
2. DNA extraction
3. Sample preparation for 16S rRNA gene sequencing
The results from Argonne National Laboratory are analyzed by a qualified bioinformatician, followed by an in-house analysis of the data using standard statistical methods. Typical microbiome analyses involve clustering of similar sequences into operational taxonomic units (OTUs) or amplicon sequence variants (ASVs) as a proxy for different microorganisms in a sample. Counts of OTUs or ASVs across samples can then be used to test for changes in within-sample (alpha) diversity and between-sample (beta) diversity. They can also be used to test for taxa that are differentially abundant across metadata categories of interest.
As shown in Figure 1, the mouse strain lacking the receptor CX3CR1 has a different gut microbiota composition. Compared to wildtype MRL/lpr, MRL/lpr-CX3CR1gfp/gfp mice (expressing the green fluorescent protein instead of CX3CR1) show significant differences in certain genera relevant to potentially pathogenic bacteria such as those in the phylum Proteobacteria. In addition, some differences were noted for "healthy" commensal bacteria such as Bifidobacterium, but not Lactobacillus.
Figure 1: Comparisons of gut microbiota at the genus level for MRL/lpr and MRL/lpr-CX3CR1 mice. The abundance of different genera at 13, 14, and 15 weeks of age (n = 5 female mice per group). Two-way ANOVA was performed. No significance ('ns'), *P < 0.05, **P < 0.01. Please click here to view a larger version of this figure.
Balanced gut microbiota can protect the human body from diseases. Once this balance is disrupted by external or internal triggers, consequences can be devastating. This method presents a way to analyze the dynamics of gut microbiota in murine models. The method is suitable not only for comparisons among groups but also for tracking the gut microbiota over time to better identify time-dependent factors disrupting the gut microbiota.
All the mice in the experiment must be handled in the same environment. One essential matter to consider while collecting fecal samples is to be consistent with the time, day, and place. Mice are coprophagic. This means that they eat their feces or the ones around them to obtain essential nutrients. Therefore, feces cannot be collected from the cage. Samples need to be collected fresh and directly from the anus. The best method is to massage the mouse while holding it or place the mouse in a container. In addition, it is recommended that samples are processed at the same time and frozen at -80 °C while waiting for all the samples to be ready.
DNA extraction is extremely sensitive to pipetting and handling, so all the samples should be processed at the same time. Another consideration is to avoid environmental contamination; therefore, every single step not necessarily performed in the chemical hood due to the toxic reagents should be done in a biological safety cabinet. Even though the chances of environmental contamination are low and there is no step to allow bacteria to grow, this method is extremely sensitive, and avoiding any possible contamination is absolutely recommended. Additionally, weighing the pellets will help to increase consistency among samples. As an advantage over commercially available kits, this method is particularly good for achieving high concentrations of DNA. The DNA yield obtained with this protocol ranges between 0.03 g and 0.07 g depending on the weight of the fecal pellet. In comparison, commercial kits can yield between 0.005 g to 0.05 g, but with a higher chance of saturation as limited by the column retention step.
The bead-beating step is important in breaking the fecal pellet into small pieces so the lysis agent can reach all the bacteria. If this step is not well performed, pellets can be partially disintegrated, and fewer representative bacteria can be reached. In addition, when transferring supernatants, correct pipetting is important, since carrying part of the residual buffer from the previous step would alter the yield and quality of DNA at the end.
When diluting the samples for the 16S ribosomal RNA gene sequencing, it is important to act quickly before the samples are put into the freezer. It is also important to mix samples well prior to performing dilutions. Avoid leaving samples at room temperature for long periods of time. Since the volumes are small, freezing the samples in the plate before shipping them will help to minimize loss.
The Argonne National Laboratory analyses the V4 region of the 16S rRNA gene. The V4 region is better conserved than the other regions and has a better estimation for high-rank taxa13. While the other regions are better for the analysis of quickly evolving species and can help identify and better discriminate genus or species, those regions accumulate mutations faster than V4.
Once results are back and analyzed by a qualified bioinformatician, statistical software can be used to plot the results and compare groups. It is important to understand the limitations of this technique and the results after bioinformatics analysis. The lower the taxonomic level, the lesser the accuracy or greater the variations. From phylum to family, and even the genus level can be analyzed with reliability. However, annotations at the species level are not accurate with this method. Shotgun sequencing, on the other hand, can reach the species level; however, while it can provide more insight into the biological functions of a community, 16S rRNA gene sequencing is still an accurate and cost-efficient method to quantify taxonomical distributions of a microbial community15.
In summary, a cost-efficient protocol has been described to analyze the gut microbiota and study its dynamics with the maximal accuracy that can help to reveal longitudinal changes over time as well as differences among groups.
The authors have nothing to disclose.
We appreciate the help from the Argonne National Laboratory and our collaborating bioinformaticians. This work is supported by various NIH and internal grants.
0.1 mm glass beads | BioSpec Products | 11079101 | |
2 mL screw cap tubes | Thermo Fisher Scientific | 3488 | |
20% SDS | FisherScientific | BP1311-1 | SDS 20% |
96% Ethanol, Molecular Biology Grade | Thermo Fisher Scientific | T032021000CS | |
Ammonium Acetate (5 M) | Thermo Fisher Scientific | AM9071 | NH4AC 5M |
B6.129P2(Cg)-Cx3cr1tm1Litt/J | Jackson Laboratory | 005582 | |
Bullet Blender storm 24 | Next Advance | 4116-BBY24M | Homogenizer |
Chloroform | FisherScientific | C298-500 | |
DEPC-Treated Water | Thermo Fisher Scientific | AM9916 | |
Ethylenediamine Tetraacetic Acid | FisherScientific | BP118-500 | EDTA |
Foil plate seal | FisherScientific | NC0302491 | |
Kimwipes-Kimtech 34256 | FisherScientific | 06-666C | |
MRL/MpJ-Faslpr/J (MRL/lpr) mice | Jackson Laboratory | 000485 | |
Nanodrop 2000 spectrophotomer | Thermo Fisher Scientific | ND2000CLAPTOP | |
Phenol: chloroform: isoamylalchohol (25:24:1) | FisherScientific | BP1752I-400 | PCI |
Scale with 4 decimals | Mettler Toledo | MS205DU | |
Skirted 96-well plates | Thermo Fisher Scientific | AB-0800 | |
Sodium chloride | FisherScientific | 15528154 | NaCl |
Tris Hydrochloride | FisherScientific | BP1757-100 | |
Vortex | Scientific Industries | SI-0236 |