Batch effect exerts a bigger influence on the rat urinary metabolome and gut microbiota than uraemia: a cautionary tale
Microbiome volume 7, Article number: 127 (2019)
Rodent models are invaluable for studying biological processes in the context of whole organisms. The reproducibility of such research is based on an assumption of metabolic similarity between experimental animals, controlled for by breeding and housing strategies that minimise genetic and environmental variation. Here, we set out to demonstrate the effect of experimental uraemia on the rat urinary metabolome and gut microbiome but found instead that the effect of vendor shipment batch was larger in both areas than that of uraemia.
Twenty four Wistar rats obtained from the same commercial supplier in two separate shipment batches underwent either subtotal nephrectomy or sham procedures. All animals undergoing subtotal nephrectomy developed an expected uraemic phenotype. The urinary metabolome was studied using 1H-NMR spectroscopy and found to vary significantly between animals from different batches, with substantial differences in concentrations of a broad range of substances including lactate, acetate, glucose, amino acids, amines and benzoate derivatives. In animals from one batch, there was a complete absence of the microbiome-associated urinary metabolite hippurate, which was present in significant concentrations in animals from the other batch. These differences were so prominent that we would have drawn quite different conclusions about the effect of uraemia on urinary phenotype depending on which batch of animals we had used. Corresponding differences were seen in the gut microbiota between animals in different batches when assessed by the sequencing of 16S rRNA gene amplicons, with higher alpha diversity and different distributions of Proteobacteria subtaxa and short-chain fatty acid producing bacteria in the second batch compared to the first. Whilst we also demonstrated differences in both the urinary metabolome and gut microbiota associated with uraemia, these effects were smaller in size than those associated with shipment batch.
These results challenge the assumption that experimental animals obtained from the same supplier are metabolically comparable, and provide metabolomic evidence that batch-to-batch variations in the microbiome of experimental animals are significant confounders in an experimental study. We discuss strategies for reducing such variability and the need for transparency in research publications about the supply of experimental animals.
The lack of reproducibility in pre-clinical animal research remains a major challenge in experimental biology  and is at least partially explained by variation between animal microbiomes . Animal research has been based on the assumption that whilst experimental animals in different facilities may have differences at species level between their gut microbiota , at a population level, in healthy laboratory animals on identical diets, these diverse collections of microorganisms achieve a shared set of basic metabolic functions—an assumption supported by evidence of significant functional redundancy within gut microbial communities .
A number of toxic molecules that accumulate in renal failure have been shown to be produced by bacterial metabolism of dietary protein in the large intestine [5, 6], leading to an interest in the gut microbiome as a potential therapeutic target to reduce the cardiovascular morbidity of patients with chronic kidney disease .
Based on an assumption of metabolic similarity between experimental animals, we sought to investigate this ‘gut-kidney axis’ in a rodent model of uraemia, by demonstrating the effect of experimental uraemia on the urinary metabolome and gut microbiota of rats, purchased from the same supplier in two separate shipment batches for logistical reasons. We actually found that the effect of shipment batch had a larger effect in both areas than uraemia and that conclusions drawn about the effect of uraemia on gut-derived metabolites would have been radically different depending on the batch of animals used.
We obtained 24 wild-type outbred Wistar International Genetic Standard (IGS) rats in two shipment batches, 3 weeks apart, from the same supplier (Charles Rivers, Kent, UK). Fourteen were rendered uraemic by undergoing a two-stage subtotal (five-sixth) nephrectomy (eight from batch 1, six from batch 2), whilst 10 underwent sham procedures (six from batch 1, four from batch 2, Fig. 1a). There were no differences in animal husbandry or diet between batches. At the time of sacrifice 8 weeks later, the urinary metabolome was assessed by untargeted proton nuclear magnetic resonance (1H-NMR) spectroscopy, and composition of the gut microbiota was assessed by sequencing 16S rRNA gene amplicons. All animals undergoing subtotal nephrectomy developed an expected uraemic phenotype, including elevations in serum urea and creatinine, weight loss, and polyuria compared to sham-operated controls, and there were no gross phenotypic differences between animals from different batches (Fig. 1b–f; Additional file 1).
Principal component analysis (PCA) of normalised and aligned urinary NMR spectral profiles identified that shipment batch was responsible for the largest source of variance in the biochemical data, seen chiefly in principal component 1, which accounted for 38% of variance. Surgical treatment accounted for a smaller but nonetheless definite source of variance, with these differences being seen chiefly in the second principal component, which accounted for 17.7% of total variance (Fig. 2a).
Separate orthogonal projection to latent structures discriminant analysis (OPLS-DA) models was constructed to elucidate biochemical variation associated with shipment batch and treatment class. The model built using shipment batch had a stronger predictive power (Q2Y = 0.66, p = 0.001) than the model built using treatment class (Q2Y = 0.48, p = 0.007). Discriminatory metabolites between the two shipment batches were identified from the OPLS-DA model (Fig. 2b), and their relative abundances were calculated from integration of the relevant regions of the aligned spectral profiles (Table 1).
Animals in batch 1 excreted significantly greater amounts of glycine (141.5 vs 68.5 relative units, Benjamini-Hochberg adjusted p < 0.001), alanine (29.3 vs 18.0 units, p < 0.001) and glucose (43.9 vs 19.7 units, p = 0.006) than animals in batch 2. They also excreted higher amounts of the potential gut bacterial products acetate (a short-chain fatty acid, 192.2 vs 105.2 units, p = 0.003), succinate (a bacterial metabolic product of dietary fibre digestion, 97.9 vs 72.6 units, p = 0.017) and lactate (571.7 vs 188.3 units, p = 0.001), compared with those in batch 2. Interestingly, hippurate was almost completely absent from the urine of batch 1 animals but present in urine from all animals in batch 2 (6.6 vs 34.5 units, p = 0.003). Correspondingly, benzoate, a gut microbially derived precursor of hippurate, was lower in the urine of batch 2 animals compared to those in batch 1 (111.0 vs 52.1 units, p < 0.001). Whilst a high degree of between-sample variation meant the batch effect did not reach the overall significance, on review of individual sample NMR spectra, it became clear that many animals had no detectable trimethylamine (TMA), a product of bacterial protein metabolism, including almost all of those in batch 1, whereas others (predominantly those in batch 2) had easily detectable concentrations.
To determine whether the substantial batch variations we had demonstrated could have led to erroneous conclusions about the effect of uraemia on the urinary metabolome, we built an OPLS-DA model for each shipment batch separately using surgical treatment class (subtotal nephrectomy vs sham) as the response variable. The model built on the batch 1 profiles was not found to be significant (Q2Y = 0.265, pQ2Y = 0.120), leading to the potential conclusion that the urinary metabolome is not influenced by uraemia. However, a significant predictive model was obtained using profiles from batch 2 (Q2Y = 0.543, pQ2Y = 0.049), despite small sample numbers, suggesting that uraemia does indeed determine urinary phenotype.
To assess whether differences in the gut microbiota between shipment batches and treatment classes might underlie these trends in the metabolomic data, sequencing of the V3 and V4 hypervariable regions of the amplified 16S rRNA gene in DNA extracted from caecal fluid was carried out. Sequence abundance data underwent isometric log-ratio transformation to allow compositional analysis of the different microbial communities.
Unsupervised PCA of the compositional data revealed that shipment batch had a larger impact on sample clustering than did treatment class (Fig. 3a). Consistent with this, permutational multivariate analysis of variance (PERMANOVA) was performed using an ADONIS analysis of a Euclidean distance matrix and confirmed that batch had a small but significant effect on the gut microbiome (R2 = 0.097, p = 0.001), whilst treatment class did not (R2 = 0.048, p = 0.227). This was further confirmed by showing that a valid predictive OPLS-DA model could be built using shipment batch as the response variable (Q2Y = 0.573, p < 0.05), but not when using treatment class (Q2Y = 0.206, p = 0.2).
The gut microbiotas of animals differed significantly in community structure between batches, with samples taken from animals in batch 2 displaying higher alpha diversity than those from animals in batch 1, across a range of measures including the inverse Simpson (40.7 vs 58.5, p = 0.043, Fig. 3b) and Shannon indices (4.53 vs 4.81, p = 0.046, Fig. 3c). Conversely, we did not demonstrate a difference in alpha diversity between uraemic and control animals.
To explore these differences more closely, populations were assessed on the basis of taxonomic assignments of OTUs at phylum, order, class, family and genus levels. Microbiotas in all animals were dominated by phyla Firmicutes (accounting for 83.1% of total reads) and Bacteroidetes (14.5%), with all other phyla (Verrucomicrobia, Tenericutes, Proteobacteria, Actinobacteria, Saccharibacteria and Deferribacteres) together representing less than 2.5% of total sequences when normalised across samples (Fig. 3d).
Differences in the abundances of OTUs and higher taxonomic groupings were analysed between shipment batches and treatment classes using the Analysis of Composition of Microbiomes (ANCOM) framework, based on isometrically log-ratio transformed abundance data and Benjamini-Hochberg adjustment for multiple hypothesis testing. Differential abundances between samples taken from animals in different shipment batches were apparent as high as at class level, with animals in batch 2 having higher relative abundances of Pseudomonadales in phylum Proteobacteria. No higher-order differences were demonstrated between uraemic and control animals.
On further analysis at OTU level, it became clear that it was primarily the less abundant OTUs which showed significant differences between batches, whilst OTUs differing significantly between uraemic and control animals were generally more abundant. Thus, whilst the relative abundance of 33/1110 OTUs (2.97% on the total) differed significantly between shipment batches, these represented only 3.80% of total sequences when analysed by the abundance of each OTU. However, the six OTUs which differed significantly between treatment classes (0.54% of the total) accounted for 5.13% of total sequences when adjusted for abundance.
These six OTUs showing significant abundance differences between uraemic and control animals were all from the family Lachnospiraceae: five from the NK4A136 group and one from the UCG-001 group. All but one showed significant decreases in relative abundance in uraemic animals, including the third most abundant OTU overall.
The 33 OTUs showing significant compositional differences between batches were drawn from five different phyla. In keeping with the higher alpha diversity seen in samples from batch 2 animals, 30/33 differentially abundant OTU between batches were seen in higher abundances in animals from this batch. Interestingly bacterial genera known to possess significant metabolic potential were prominently represented amongst these differentially abundant organisms, including a number of producers of short-chain fatty acids (Roseburia, Butyricicoccus, Butyrivibrio and Acetomaculum) and three from the phylum Proteobacteria.
The rodent gut microbiome is a complex community of several hundred different bacterial species that possess significant metabolic potential of immense relevance to the host organism. It has previously been demonstrated that this community differs according to a variety of factors including host age  and genetics [10, 11], caging arrangements [9, 11, 12], bedding material and water sterilisation technique  and vendor shipment batch . Xiao et al. generated a catalogue of the mouse metagenome by sequencing faecal material from 184 mice and found that vendor was a prime determinant in variation at a genetic and function level .
In light of these studies, we have demonstrated that predicted batch variations in gut microbiota are associated with multiple, major variations in a range of urinary metabolites, with the potential for significant downstream effects on wider areas of host phenotype. For example, circulating hippurate has recently been suggested as a biomarker for gut microbial diversity, associating with the risk of metabolic syndrome ; however, our results suggest it may be totally absent in the urine of experimental animals based on shipment batch. Likewise, the biological relevance of dietary amines has been demonstrated through the association of TMA and its metabolite trimethylamine N-oxide with cardiovascular disease [16, 17], including in patients with chronic kidney disease . However, our results suggest that rats purchased from the same supplier in different shipment batches may metabolise dietary amines in quite different ways, potentially questioning the generalisability of research based on individual batches of animal subjects.
Since the diet of animals in each group was identical, we conclude that differences in bacterial metabolic pathways are likely to underlie these differences in the urinary metabolome. We demonstrated batch differences in the relative abundances of a number of bacteria that are of known metabolic significance, including several that are major sources of short-chain fatty acids and associated with beneficial health outcomes [19,20,21], and several from the phylum Proteobacteria that has recently been shown to contribute significantly to functional variation between gut metagenomes .
These results challenge the assumption that in healthy organisms, different microbial communities achieve a common set of basic metabolic functions despite variation in the individual species present [23, 24]. It can no longer be assumed that healthy laboratory animals, purchased from the same supplier, are metabolically similar. The inherent microbial dissimilarity and associated metabolic differences between animals in different batches provide a significant source of experimental variation.
Such batch variations could easily lead to spurious positive results. For example, a group that demonstrates an effect in response to an experimental intervention with a small group of animals may decide to increase the number of animals in order to publish their findings; they purchase new animals from the same supplier, but fail to reproduce their earlier results because the new additions have significantly different microbial metabolic potential. Even worse, they may have carried out interventional procedures on one batch of animals, and then used animals from a different batch as controls, with exaggerated differences between groups reflecting underlying differences in microbiomes rather than any effect of the experimental procedure. The alternative in each case—to re-run the whole experiment with animals purchased in a new, single batch—may be prohibitively expensive, may fail to reproduce the initial results and seems to stand against the second of the ‘Three R’s’ governing ethical use of animals in research: the reduction of the number of animals used .
It is crucial that publishers maintain the requirement to document fully all aspects of animal use, including purchase details of the different batches of animals used in a study if these come from a commercial supplier. Furthermore, steps should be taken to reduce the amount of variation within batches, such as by using a standardised procedure for moving bedding between cages, which has been shown to reduce intra-batch variation . Many experimental groups breed their own animals, which may reduce intra-group variation, although potentially at the expense of generalisability with results from other laboratories. Statistical approaches including percentile normalisation have been suggested that would allow pooling of data between different batches on experimental subjects in different settings, although for this to be successful, large numbers of control subjects are required .
Finally, batch variation can be embraced as a reflection of real-world microbial variation. For this to be successful, it is important that researchers use unsupervised PCA plots of all experimental subjects, coloured according to batch, for quality control, using statistical methods to measure the effect of batch variation. Documenting whether the same observed changes were seen in animals from all batches, or whether different batches behaved differently, is helpful in assessing the generalisability of results, and many journals already have such requirements in place. An experiment showing the same effect in two or more smaller but separate batches of animals may be more striking than an experiment showing a larger effect size in a single animal batch.
Animal experiments were conducted in accordance with the UK Home Office Animals (Scientific Procedures) Act 1986, with local ethical committee approval. All animal works were carried out at the Biological Services Unit of Queen Mary University of London at Charterhouse Square, and complied fully with all relevant animal welfare guidance and legislation. The 24 male, outbred Wistar IGS rats were obtained from Charles Rivers (Kent, UK) in two shipment batches 3 weeks apart. All were housed in individually ventilated cages under 12-h light/dark cycles and were allowed unlimited access to water and chow (RM1 diet from Special Diet Services, Essex, UK).
After a week-long period of acclimatisation, rats underwent a two-stage surgical procedure involving either subtotal nephrectomy or a sham procedure. Subtotal nephrectomy involved exteriorisation of the left kidney with decapsulation and removal of the upper and lower poles and subsequent replacement of the middle pole only, followed by total right nephrectomy 2 weeks later. Sham procedures involved exteriorisation, decapsulation and replacement of the left kidney, followed by the same procedure on the right kidney 2 weeks later.
Following surgery, rats were weighed weekly. There were up to four rats per cage, and the animals were initially housed according to surgical intervention (subtotal nephrectomy or sham) for 2 weeks after the second stage surgery, before some were moved into mixed cages comprising both uraemic and control animals (this was in order to assess the role of individual cage variants, which on subsequent analysis, not presented here, was found not to be as important in explaining the key changes in urinary metabolome or gut microbiome as shipment batch or treatment class). There was no co-housing between batches. Each week, the animals were housed individually in metabolism cages to allow the collection of a 24-h urinary specimen which was frozen at − 80 °C until the time of analysis. Rats were killed by lethal injection of sodium thiopentone (LINK Pharmaceuticals, Horsham, UK), and caecal fluid was expressed, stored in foil and snap-frozen in liquid nitrogen and then at − 80 °C until the time of analysis. Blood samples were taken by cardiac puncture, and after centrifugation, the serum was frozen at − 80 °C until the time of analysis.
Quantification of serum urea and creatinine was done by IDEXX Bioresearch, Ludwigsberg, Germany.
Urine samples were randomised prior to dilution with buffer and running on the machine to remove potential for technical batch effects in processing and analysis, and prepared for 1H-NMR spectroscopy as described previously . All samples were analysed on an NMR spectrometer (Bruker) operating at 600.22 MHz 1H frequency.
Processing of NMR data
The NMR spectral profiles were digitised and imported into Matlab (Mathworks) using in-house scripts (Additional file 3). The raw spectra were adjusted for 24-h urine volumes by multiplying all NMR absorbance values by the urine volume in millilitres. The peaks for water and trimethylsilylpropanoic acid (TSP) were excised from the raw NMR spectra which were then aligned to adjust for variation in peak shift due to pH differences. Further normalisation was carried out using the probabilistic quotient method between samples in order to ensure comparable baselines between samples (Additional file 2).
Unsupervised PCA was used to identify sources of variation in the metabolic data. This was followed by supervised OPLS-DA analysis using both shipment batch and treatment class as the response variable. In-house-developed scripts were used to perform these multivariate statistical analyses. Peak integrals were calculated from metabolite peaks identified as discriminatory from the OPLS-DA models. Comparisons between these integrals were used to calculate differences in relative abundance according to shipment batch and treatment class using Microsoft Excel, with the Student’s t test and Welch’s correction used to assess significance. These p values were adjusted using the Benjamini-Hochberg method  and a false discovery rate of 0.15 using the q values  package in R (Additional file 4).
16S rRNA gene sequencing and analyses
DNA was extracted from samples of caecal fluid using the DNeasy PowerSoil kit from QIAGEN, used according to the manufacturer’s instructions. All samples were processed using the same kit, and a negative ‘kitome’ control was also included with samples . DNA diluted to 10 ng/μL (in 10 mM Tris HCl pH 8.5) was submitted to the Centre for Genomic Research at the University of Liverpool for library preparation and sequencing of the V3/V4 hypervariable region of the 16S rRNA gene. Sequence data were processed using QIIME v1.9 . Paired-end data were joined using join_paired_ends.py, and primer sequences removed from split library files using cutadapt . OTUs were picked using 99% BLAST identity using usearch; from these, a representative set of OTUs was selected. Sequences were aligned (PyNAST) against Silva v128 , and this database was also used to assign taxonomy. Singletons, mitochondria-, cyanobacteria- and control-associated OTUs were removed from the OTU table, as were OTUs unaffiliated with any taxonomic group. Data were then rarefied to 100,000 reads to account for differences in sequencing depth across samples and these relative abundances were used to calculate the overall abundances by phylum presented in Fig. 3d.
Raw (unrarefied) OTU abundance data were imported into R for analyses using Phyloseq  (Additional files 5 and 6). A phylogenetic tree was generated using MEGA v7.0  and rooted to a random node using the R package phytools . A pseudocount of 0.001 was added to all OTU abundances to avoid calculating log-ratios involving zeros, and then data was then made compositional through isometric log-ratio transformation using the R package philr . Ordination was carried out using the ‘ordinate’ function in Phyloseq, based on Euclidean distances in philr space. Permutational analysis of variance (PERMANOVA) was carried out using the ADONIS command in the R package vegan . OPLS-DA models were built using the ropls package in R . Alpha diversity was assessed using Phyloseq. Compositional analysis of the microbiota at six taxonomic levels was based on isometric log-ratio transformation of raw sequence abundances and adjusted for multiple testing using the Benjamini-Hochberg method, carried out using the ANCOM statistical framework  in R, with code obtained from the author’s webpage: https://sites.google.com/site/siddharthamandal1985/research.
Preparation of figures
In order to achieve uniformity, most figures except those demonstrating NMR spectral data were generated using GraphPad Prism 7 (GraphPad Software Inc., San Diego, California). NMR spectra and related figures were created using Matlab (Mathworks) with in-house scripts.
Availability of data and materials
Animal data (weight, urine output and serum biochemistry) are submitted as Additional file 1.
The raw NMR data is available as Additional file 2.
The rarefied 16S rRNA gene sequence data including the negative ‘kitome’ control have been deposited with links to BioProject accession number PRJNA525754 in the NCBI BioProject database (https://www.ncbi.nlm.nih.gov/bioproject/). Raw (non-rarefied) data are available as Additional file 5.
International Genetic Standard
Nuclear magnetic resonance
Orthogonal projection to latent structures discriminant analysis
Operational taxonomic unit
Principal component analysis
Collins FS, Tabak LA. NIH plans to enhance reproducibility. Nature. 2014;505:612–3.
Stappenbeck TS, Virgin HW. Accounting for reciprocal host-microbiome interactions in experimental science. Nature. 2016;534:191–9.
Parker KD, Albeke SE, Gigley JP, Goldstein AM, Ward NL. Microbiome composition in both wild-type and disease model mice is heavily influenced by mouse facility. Front Microbiol. 2018;9:1598.
Heintz-Buschart A, Wilmes P. Human gut microbiome: function matters. Trends Microbiol. 2018;26(7):563–74.
Aronov PA, Luo FJ, Plummer NS, Quan Z, Holmes S, Hostetter TH, Meyer TW. Colonic contribution to uremic solutes. J Am Soc Nephrol. 2011;22(9):1769–76.
Mishima E, Fukuda S, Mukawa C, Yuri A, Kanemitsu Y, Matsumoto Y, et al. Evaluation of the impact of gut microbiota on uremic solute accumulation by a CE-TOFMS-based metabolomics approach. Kidney Int. 2017;92(3):634–45.
Ramezani A, Raj DS. The gut microbiome, kidney disease, and targeted interventions. J Am Soc Nephrol. 2014;25(4):657–70.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Royal Stat Soc. 1995;57(1):289–300.
Lees H, Swann J, Poucher SM, Nicholson JK, Holmes E, Wilson ID, Marchesi JR. Age and microenvironment outweigh genetic influence on the Zucker rat microbiome. PLoS One. 2014;9(9):e100916.
Ericsson AC, Davis JW, Spollen W, Bivens N, Givan S, Hagan CE, et al. Effects of vendor and genetic background on the composition of the fecal microbiota of inbred mice. PLoS One. 2015;10(2):e0116704.
Hildebrand F, Nguyen TL, Brinkman B, Yunta RG, Cauwe B, Vandenabeele P, et al. Inflammation-associated enterotypes, host genotype, cage and inter-individual effects drive gut microbiota variation in common laboratory mice. Genome Biol. 2013;14(1):R4.
Swann JR, Tuohy KM, Lindfors P, Brown DT, Gibson GR, Wilson ID, et al. Variation in antibiotic-induced microbial recolonization impacts on the host metabolic phenotypes of rats. J Proteome Res. 2011;10(8):3590–603.
Bidot WA, Ericsson AC, Franklin CL. Effects of water decontamination methods and bedding material on the gut microbiota. PLoS One. 2018;13(10):e0198305.
Xiao L, Feng Q, Liang S, Sonne SB, Xia Z, Qiu X, et al. A catalog of the mouse gut metagenome. Nat Biotechnol. 2015;33(10):1103–8.
Pallister T, Jackson MA, Martin TC, Zierer J, Jennings A, Mohney RP, et al. Hippurate as a metabolomic marker of gut microbiome diversity: modulation by diet and relationship to metabolic syndrome. Sci Rep. 2017;7(1):13670.
Tang WH, Wang Z, Levison BS, Koeth RA, Britt EB, Fu X, et al. Intestinal microbial metabolism of phosphatidylcholine and cardiovascular risk. N Engl J Med. 2013;368(17):1575–84.
Qi J, You T, Li J, Pan T, Xiang L, Han Y, Zhu L. Circulating trimethylamine N-oxide and the risk of cardiovascular diseases: a systematic review and meta-analysis of 11 prospective cohort studies. J Cell Mol Med. 2017;22(1):185–94.
Stubbs JR, House JA, Ocque AJ, Zhang S, Johnson C, Kimber C, Schmidt K, et al. Serum trimethylamine-N-oxide is elevated in CKD and correlates with coronary atherosclerosis burden. J Am Soc Nephrol. 2016;27(1):305–13.
Tamanai-Shacoori Z, Smida I, Bousarghin L, Loreal O, Meuric V, Fong SB, et al. Roseburia spp.: a marker of health? Future Microbiol. 2017;12:157–70.
Ohkawara S, Furuya H, Nagashima K, Asanuma N, Hino T. Oral administration of Butyrivibrio fibrisolvens, a butyrate-producing bacterium, decreases the formation of aberrant crypt foci in the colon and rectum of mice. J Nutr. 2005;135(12):2878–83.
Eeckhaut V, Wang J, Van Parys A, Haesebrouck F, Joossens M, Falony G, et al. The probiotic Butyricicoccus pullicaecorum reduces feed conversion and protects from potentially harmful intestinal microorganisms and necrotic enteritis in broilers. Front Microbiol. 2016;7:1416.
Bradley PH, Pollard KS. Proteobacteria explain significant functional variability in the human gut microbiome. Microbiome. 2017;5(1):36.
Turnbaugh PJ, Hamady M, Yatsunenko T, Cantarel BL, Duncan A, Ley RE, et al. A core gut microbiome in obese and lean twins. Nature. 2009;457(7228):480–4.
Human Microbiome Project Consortium. Structure, function and diversity of the healthy human microbiome. Nature. 2012;486(7402):207–14.
Russell WMS, Burch RL. The principles of humane experimental technique. London: Methuen; 1959.
Miyoshi J, Leone V, Nobutani K, Musch MW, Martinez-Guryn K, Wang Y, Miyoshi S, et al. Minimizing confounders and increasing data quality in murine models for studies of the gut microbiome. Peer J. 2018;6:e5166.
Gibbons SM, Duvallet C, Alm EJ. Correcting for batch effects in case-control microbiome studies. PLoS Comput Biol. 2018;14(4):e1006102.
Dona AC, Jiménez B, Schäfer H, Humpfer E, Spraul M, Lewis MR, Pearce JT, et al. Precision high throughput proton NMR spectroscopy of human urine, serum, and plasma for large-scale metabolic phenotyping. Anal Chem. 2014;86:9887–94.
Storey JD. A direct approach to false discovery rates. Journal of the Royal Statistical Society, Series B, 64:479–498, 2002.
Salter SJ, Cox MJ, Turek EM, Calus ST, Cookson WO, Moffatt MF, et al. Reagent and laboratory contamination can critically impact sequence-based microbiome analyses. BMC Biol. 2014;12:87.
Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7(5):335–6.
Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011;17(1):10–2.
Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41(D1):D590–6.
McMurdie PJ, Holmes S. Phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One. 2013;8(4):e61217.
Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33:1870–4.
Revell LJ. Phytools: an R package for phylogenetic comparative biology (and other things). Methods Ecol Evol. 2012;3:217–23.
Silverman JD, Washburne AD, Mukherjee S, David LA. A phylogenetic transform enhances analysis of compositional microbiota data. eLife. 2017;6. https://doi.org/10.7554/eLife.21887.
Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, McGlinn D, et al. Vegan: community ecology package. R package version 2.5–2. https://CRAN.R-project.org/package=vegan. Accessed 16 June 2019.
Thevenot EA, Roux A, Xu Y, Ezan E, Junot C. Analysis of the human adult urinary metabolome variations with age, body mass index and gender by implementing a comprehensive workflow for univariate and OPLS statistical analyses. J Proteome Res. 2015;14:3322–35.
Mandal S, Van Treuren W, White RA, Eggesbø M, Knight R, Peddada SD. Analysis of composition of microbiomes: a novel method for studying microbial composition. Microb Ecol Health Dis. 2015;26.
Analysis of sequencing data used the computing resources of the UK MEDical BIOinformatics partnership—aggregation, integration, visualisation and analysis of large, complex data (UK MED-BIO), which is supported by the Medical Research Council (grant number MR/L01632X/1).
The direct experimental work reported here (including the animal work, 1H-NMR spectroscopy and genetic work) was supported by the Barts and the London Renal Research Fund. Subsequent analysis of sequencing data used the computing resources of the UK MEDical BIOinformatics partnership—aggregation, integration, visualisation and analysis of large, complex data (UK MED-BIO), which is supported by the Medical Research Council (grant number MR/L01632X/1). Neither funding body had any direct involvement in the experimental design or the interpretation and writing up of results.
Animal experiments were conducted in accordance with UK Home Office Animals (Scientific Procedures) Act 1986, with local ethical committee approval. All animal works were carried out at the Biological Services Unit of Queen Mary University of London at Charterhouse Square, and complied fully with all relevant animal welfare guidance and legislation.
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.
Animal data (XLSX 11 kb)
Raw NMR data (CSV 3674 kb)
MATLAB code for NMR analysis (DOCX 17 kb)
R code for NMR analysis (DOCX 17 kb)
Non-rarefied 16S abundance data (CSV 77 kb)
R code for microbiome analysis (DOCX 25 kb)
About this article
Cite this article
Randall, D.W., Kieswich, J., Swann, J. et al. Batch effect exerts a bigger influence on the rat urinary metabolome and gut microbiota than uraemia: a cautionary tale. Microbiome 7, 127 (2019). https://doi.org/10.1186/s40168-019-0738-y