BACKGROUND
neonatal piglet diarrhoea is of significant importance for the pig industry because it causes economic losses due to increased morbidity and mortality, decreased weight gain, and the need for extra medications  <cit> .
obviously, it impairs the welfare of the animals in the short term, but it may also affect their health in the longer term as a consequence of disrupting the normal bacterial succession in the gastrointestinal tract  <cit> .
bacterial colonization of the mammalian gastrointestinal tract begins at birth  <cit> .
this colonization is a dynamic event, and the bacterial succession is influenced by a number of factors including: mode of delivery, surrounding environment, gestational age, and genetics  <cit> .
the colonization of the gut has a major impact on the host’s health and disease; for example, the microbiota helps in the maturation of the gastrointestinal tract and immune system, protects against pathogen colonization through competitive exclusion, and converts otherwise indigestible substances into digestible components that benefit the host  <cit> .
in accordance with these roles, studies of germ-free animals, deficient of normally developed gut microbiota, show animals with an impaired and immature intestinal immune system, in addition to changes in intestinal morphology  <cit> .
during the last ten years, increasing attention has been focused on a new neonatal porcine diarrhoea  <cit> .
in denmark, this diarrhoea has been named ‘new neonatal porcine diarrhoea’ .
what distinguishes nnpd from other types of neonatal piglet diarrhoea are the following: 1) the aetiology is unknown, however, routine diagnostic protocols show that it is not caused by known enteric pathogens such as hemolytic escherichia coli, clostridium difficile, clostridium perfringens type a or c, coronavirus, rotavirus species a or c, cryptosporidium spp., cystoisospora suis, giardia spp., or strongyloides ransomi.
2) typical strategies, such as vaccination against enterotoxigenic e. coli and cl.
perfringens type c or treatment with antibiotics, do not seem to have a noteworthy effect on the diarrhoea.
3) no obvious connection between nnpd and pig farm health status or management has been demonstrated .
it is difficult to estimate how widespread nnpd is, mainly because of the unknown aetiology combined with a fluctuating clinical presentation  <cit> , as well as differences in routine laboratory testing  <cit> , but also because of the limited number of studies focusing on this issue.
nonetheless, a diarrhoea of much resemblance to nnpd has been described in sweden and france  <cit> .
this study is part of an interdisciplinary project investigating the aetiology of nnpd.
kongsted et al.
suggested the following case-definition of nnpd: “non-hemorrhagic diarrhoea during the first week of life, without detection of known infectious pathogens, characterized by milk-filled stomachs and flaccid intestines at necropsy”  <cit> .
the same author found dissimilarities in the course and severity of nnpd among four pig farms and estimated that affected piglets had a negative average daily weight gain with an increased risk of dying, though this risk was not significant  <cit> .
jonach et al.
examined the role of four enteric bacterial pathogens by fluorescence in situ hybridization  and found that simultaneous colonization of the intestinal mucosa with non-enterotoxigenic e. coli  and enterococcus spp.
could be involved in the pathogenesis of nnpd  <cit> .
finally, several different viral assays tested negative on samples from some of the same nnpd animals that were examined in the aforementioned studies, indicating that common known viruses likely do not contribute to nnpd .
this study investigates whether nnpd is associated with the composition of the gut microbiota obtained from piglets with and without diarrhoea.
this was examined by evaluating the overall bacterial composition and relative quantitative distribution of ileal and colonic intestinal content using the gut microbiotassay: an assembly of  <dig> primer sets targeting ribosomal rna  genes , verified to function with the high-throughput quantitative real-time pcr  chip access array™ integrated fluidic circuit  from fluidigm  <cit> .
as the name implies, the assay is designed to target major bacteria phyla and selected taxonomic sub-groups of the gut microbiota.
the gut microbiotassay provides a quick overview of the distribution, as well as the relative quantity of the gut microbiota in a large number of samples simultaneously.
subsequently, pcr amplicons from ileum and colon of four case piglets and four control piglets were sequenced using 454-technology to acquire deeper taxonomic information.
this approach revealed diverse gut microbial profiles associated with piglet diarrhoeic status.
methods
animals and sample collection
danish pig farms affected by nnpd were identified from conversations with veterinarians and farm managers.
four pig farms that fulfilled the inclusion criteria listed in were included in the study .
on each pig farm, approximately  <dig> randomly chosen sows from one farrowing batch  were followed for a seven day period following farrowing.
all newborn piglets were weighed at the beginning of the trial , and piglets weighing less than 800 g were excluded.
all animals were subject to a daily clinical examination that paid special attention to fecal appearance on rectal swabs.
diarrhoea was defined as loose or watery feces.
based on these observations, piglets were characterized as either cases or controls: a case piglet had suffered from diarrhoea for at least two consecutive days, whereas a control piglet had never experienced diarrhoea.
the inclusion criterion applied to case piglets resulted in a greater number of case piglets from gilts , as piglets born to multiparous sows did not suffer from diarrhoea to the same extent as piglets born to gilts.
control piglets were as far as possible collected from litters without diarrhoea.
in total  <dig> control piglets and  <dig> case piglets were selected  and brought to the danish pig research centre, kjellerup, for euthanization and necropsy.
for ethical reasons, farmers were allowed to treat piglets for any disorder if necessary, on the condition that any such treatments were recorded for individual animals.
during necropsy, intestinal contents were collected from the distal small intestine  and the large intestine  of each animal and stored at −80°c until further analysis.table  <dig>
piglets included in the study
age in days
3
4
5
6
7
total
gilts
sows
diarrhoea
−
+
−
+
−
+
−
+
−
+
-
+
-
+
pig farm 1
0
4
13
9
pig farm 2
2
12
10
0
pig farm 3
0
10
13
4
pig farm 4
4
8
8
5
total
14
15
9
10
17
17
6
6
4
4
6
34
44
18
the number of piglets is listed for each pig farm in accordance with age in days:  control piglets without diarrhoea,  case piglets with diarrhoea for ≥2 days.
the column “total” summarizes the number of piglets born to gilts or sows  by farm.
pig farm inclusion criteria for studying nnpd <cit> :a stock of more than  <dig> sows
routinely vaccinate <dig> against etec <dig> and cl.
perfringens type c
farrowing units proven to be prrs <dig> virus negative as demonstrated in blood samples tested by elisa4/ipt <dig> or pcr.
nnpd has been a problem for a period of minimum 3– <dig> months
management-related diarrhoea has been excluded
traditional interventions such as vaccines and antibiotics have limited effect on the diarrhoea
at least 30% of the litters are affected
the piglets suffer from diarrhoea during the first week of life
random check of five case piglets should test negative for etec, cl.
perfringens type c and rotavirus by routine diagnostic examination
dna extraction
based on development work comparing various extraction methods, the maxwell®  <dig> lev blood dna kit  was determined to be the optimal method, based on providing sufficient dna yields and acceptable purities.
according to manufacturer’s recommendations, total dna was extracted from intestinal content using this kit.
a total of 200 mg of intestinal contents were suspended in 600 μl pbs and vortexed until visually homogeneous.
the samples were centrifuged for 2 min at 200 × g, and the supernatant was transferred to new tubes.
a volume of 350 μl of lysis buffer was then added, and the bacterial cells were lysed by bead beading  for  <dig>  min at  <dig>  hertz with a 5 mm steel bead .
subsequently, samples were centrifuged for 1 min at 1000 × g at 4°c, after which the supernatant was transferred to novel 2 ml tubes and mixed manually with 30 μl proteinase k .
next, samples were incubated for 30 min at 56°c, and centrifuged at  <dig> 000 × g for 1 min.
the entire suspension was transferred to the sample inlets on the cartridge, and 50 μl of elution buffer was added to the bottom of the elution tubes.
the cartridges were prepared according to the manufacturer’s instructions  <cit>  and the settings “research mode”, “lev mode”, “dna”, and “blood/cell” were selected for dna extraction using the maxwell®  <dig> instrument .
finally, tubes were centrifuged at  <dig>  <dig> × g for 3 min to settle any magnetic bead leftovers, and the dna was moved to new tubes.
dna concentration and purity was measured using the nanodrop® nd- <dig>  spectrophotometer, and dna was stored at −20°c until further analysis.
analysis of bacterial profiles with use of the gut microbiotassay
to obtain a bacterial profile of the intestinal content, this study used the gut microbiotassay designed for the aa <dig>   <cit>  with a few modifications.
in order to improve the performance of the primer sets “domain bacteria b v4-v5” and “phylum firmicutes” degenerate nucleotides were introduced: nucleotide  <dig> was changed from c to m , and nucleotide  <dig> from t to y , , in the forward or reverse primer, respectively .
all remaining primer sets were identical to those previously published as ‘the gut microbiotassay’  <cit> .
primers were purchased from eurofins mwg synthesis gmbh  and stored at −20°c.
in brief, the aa <dig>  is a qpcr platform that is capable of running 48 × 48 =  <dig> individual reactions simultaneously, and it enables quick and easy library preparation for  <dig> sequencing  <cit> .
the aa <dig>  was processed and prepared following the ‘access array system™ user guide’  <cit> , with and without adding the access array barcode library for the  <dig> gs flx titanium sequencer , as previously described  <cit> .
primers targeting bacteria at the taxonomic level of species were not tagged or included for sequencing, as these amplicons were regarded as contributing little information due to their specificity.
however, analysis at the species-level was possible using the information obtained at higher taxonomic levels.
all samples were diluted to 50 ng/μl with nuclease-free water .
primers were diluted to 4 μm with 20× access array loading buffer and nuclease-free water .
master mix was a mixture of: 10× faststart high fidelity reaction buffer with 18 mm mgcl <dig> , 25 mm mgcl <dig> , dmso , 20× access array loading reagent , 50× rox , 20× evagreen® , 10 mm pcr grade nucleotide mix , and  <dig> u/μl faststart high fidelity enzyme blend , in final concentrations of 1×,  <dig>  mm, 5%, 1×,  <dig> ×, 1×, 200 μm, and  <dig>  u/μl, respectively.
sample mix was prepared from 3 μl master mix, 1 μl 454bl  , and 1 μl dna  as the very last step before running the array ).
when the pcr reaction had finished, the barcoded pcr amplicons were harvested and stored at −20°c.
next generation sequencing
samples from two representative animals were selected from each farm: one piglet with nnpd and one piglet without nnpd, giving a total of  <dig> samples to be sequenced .
piglets were chosen to be approximately the same age .
concentrations of the respective pcr amplicons for each animal were determined using the agilent  <dig> chip .
amplicons were subsequently pooled in equal concentrations and size-separated by running the amplicons for 86 min, 90 v, in a  <dig> % seakem® le agarose gel  followed by incubation for 30 min in  <dig> % ethidium bromide for staining.
by means of uv radiation from the bio-rad universal hood ii  gel bands were visualized and bands spanning the size range of the primer products  were excised.
finally, the qiaquick gel extraction kit  was used to purify dna from the gel.
this pool of  <dig>  ng dna  derived from ileal and colonic luminal contents of  <dig> different animals was sequenced on a half picotiterplate™ by a  <dig> gs flx titanium sequencer  via lgc genomics .
data analysis
relative quantification of quantification cycle  values
data analysis was conducted as described in a previously published methodology article  <cit> : raw cq values were exported from ‘fluidigm real-time pcr analysis’ software version  <dig> .
<dig> to microsoft excel.
to even out possible variation between the aa <dig>  runs, all cq values were normalized to an interplate calibrator.
next, cq values from each sample were normalized to their respective “domain bacteria b v4-v5” primer set, thereby calculating the relative quantification.
principal component analysis on normalized cq data
principal component analysis  was conducted in the software package r, version  <dig> .
<dig>  <cit> .
all normalized data were initially transformed with the natural logarithm .
next, primer sets with less than 50% recorded cq values were excluded from the analysis, and in the same manner, samples that resulted in less than 50% registered cq values were removed.
for the remaining primer sets, missing values were substituted with the lowest registered primer-specific value.
data were scaled by the individual primer mean and standard deviation for pca.
the number of components to be investigated further was determined by visual inspection of scree plots of eigenvalues; based on an inflection point or when the components explained > 60% variance .
all possible combinations of included principal components were visualized to ensure that no important clustering was overlooked.
the number of components used for graphical representation in the pca was set at two to ease interpretation.
samples used for  <dig> sequencing were pinpointed from the pca results.
in order for samples to be considered for sequencing the samples were represented with luminal content from both ileum and colon, and none of these were allowed to have outlying coordinates in the pca, but were somewhat in the centre of its respective category .
linear discriminant analysis  based on leave-one-out cross-validation was used to classify observations into cases or controls by their fluidigm-derived microbial profile and to estimate how well this classification corresponded with the actual status of the sample.
this classification was followed up with a multivariate analysis of variance  using wilks lambda statistic.
statistics on the gut microbiotassay
bacterial profiles obtained from the gut microbiotassay were used to study the gut microbiota of piglets with and without nnpd.
for the statistical calculations all values were initially ln-transformed.
a linear mixed-effect model was then applied to each primer set separately with randomized effect of herd origin:  <dig> lnyi=agutsectioni,statusi+agilti+atreatmenti+βstatusi⋅agei+β⋅diarrhoeai+η⋅yherdi+εi,i= <dig> …, <dig> , where status indicates the status of the i’th animal, and where yherd is a normalized gaussian stochastic variable indicating the randomized effect of herd, indexed by the herd of the i’th animal, so that animals from the same herd share the random effect.
the model is applied for different responses y, as listed in the first column in table  <dig>  the following variables were included as deterministic effects in the model: gut section , status , gilt , percentage of diarrhoea , treatment , age .
the interaction between gut section and status was also included to examine if possible diarrhoea-causing agents could be traced to a specific gut section.
additionally, the interaction between status and age was included to study whether the effect of age differed between case and control piglets.
for each primer set as response variable, the model was reduced by backwards stepwise elimination, resulting in a series of primer-specific final models.
model reduction was performed as likelihood ratio tests.
model fit for each primer set was assessed through residual analysis of the model residuals corrected for random effects and the estimated random effects themselves, using standard techniques.
these analyses were consistent with standard model behavior.
the effect of gender and birth weight in relation to status was also evaluated.
significances in the final models were reported through p-values.table  <dig>
significant findings from the gut microbiotassay and  <dig> sequencing results displaying similar trends
primer set and respective  <dig> sequencing findings
gut section estimate
gut section
p-
value
status estimate
status
p
-value
gilt estimate
gilt
p
-value
percentage of diarrhoea estimate
#
percentage of diarrhoea
p
-value
age estimate
age
p
-value
antibiotics estimate
antibiotics
p
-value
domain bacteria b
268
phylum firmicutes
183
class bacilli
125
genus
enterococcus
27
genus
streptococcus
47
genus
lactobacillus
79
family clostridium cluster i
17
species
cl.
perfringens
family clostridium cluster iv
64
family clostridium cluster xiv
81
phylum bacteroidetes
126
genus
bacteroides
65
species bacteroides pyogenes †
species bacteroides rodentium †
species bacteroides xylanisolvens †
species porphyromonadaceae bacterium c <dig> †
phylum actinobacteria
74
family bifidobacteriaceae
12
phylum fusobacteria
16
class beta- and gammaproteobacteria
109
family enterobacteriaceae
89
species
e. coli
class deltaproteobacteria
1
species desulfovibrio piger †
estimated fold changes as the estimates of the exponentiated effect parameter, corresponding 95% confidence intervals , and p-values.
superscripted numbers equal the total number of otus belonging to the respective taxonomic group identified by sequencing amplicons generated by the gut microbiotassay.
otu-ratio between diarrhoeic and healthy piglets for ileum and colon, respectively, listed in parenthesis.
†indicate a significant effect of interaction .
species mark with † are not consistent findings.
#estimated fold changes for this category are calculated as: γ
percentage points, where the parameter y from column  <dig> is the estimate of the exponential eβ of the estimated effect parameter β for diarrhoea in formula .
the diarrhoeic percentage points range from  <dig> to  <dig>
previous studies have suggested that simultaneous colonization with e. coli and enterococci species may be a cause of neonatal porcine diarrhoea  <cit> .
this co-occurrence was also implied using the aforementioned model ), thus, this was further investigated.
a logistic regression model with randomized effect of herd origin, fit by penalized quasi-likelihood, was used to determine if the odds of nnpd was associated with the following variables: gut section, gilt, age, genus enterococcus, class beta- and gammaproteobacteria, family enterobacteriaceae, species e. coli, as well as possible interaction in between bacteria, and between bacteria and gut section:table  <dig>
estimated fold changes, 95% confidence intervals , and
p-
values of significant interactions from the gut microbiotassay
primer set
gut section: status estimate
gut section: status estimate
gut section: status
p-
value
status: age estimate, day 3
status: age estimate, day 4
status: age estimate, day 5
status: age estimate, day 6
status: age estimate, day 7
status: age estimate, mean age
status: age
p
-value
family clostridium cluster i
species
cl.
perfringens
phylum bacteroidetes
genus
bacteroides
class beta- and gammaproteobacteria
family enterobacteriaceae
species
e. coli
this table is an extension of table  <dig>  and the estimated fold changes are the estimates of the exponentiated effect parameter.
<dig> logitpi=agutsectioni+agilti+βage⋅agei+∑j=14βj,gutsectionibacteriaj,i+∑k=14∑ℓ≠kβk,lbacteriak,ibacteriaℓ,i+ηyherdi , where bacteria1-bacteria <dig> are genus enterococcus, class beta- and gammaproteobacteria, family enterobacteriaceae, and species e. coli, respectively, where bacteriaj,i is the value of bacteriaj animal i, where agei is the age in days of animal i, and where yherd is the random effect component corresponding to the herd of animal i.
the models were validated through graphical procedures.
a p-value of less than  <dig>  was considered significant.
statistical analyses were performed in r  <cit> .
bioinformatics analysis of  <dig> sequencing data
sequencing data available at ncbi sequence read archive  were analyzed using the as yet unpublished open source package bion-meta  .
this program facilitates quick and easy bioinformatics analysis, and bion-meta has previously been applied to a similar dataset  <cit> .
however, in addition to matching all sequences against the ribosomal small subunit  silva dataset, suited for 16s rrna gene targeting primer sets, bion-meta was advanced to encompass searches in the ribosomal large subunit  silva dataset, thereby accommodating the primer sets targeting the 23s rrna gene used in the gut microbiotassay  <cit> .
the bion-meta workflow included the following elements: 1) separation by sample barcodes and primer tags; 2) removal of primer remnants and bases at the ends with lower than 96% quality , as well as sequence filtering by length  and filtering by minimum base quality of 96%.
; 3) removal of chimeric sequences; 4) separation of each sample by matching it with the phylogenetic primer; 5) matching all sequences against the ssu and lsu silva datasets and producing a table with the highest 1% similarities for each query; and 6) mapping the similarities to the silva ssu and lsu taxonomies, identifying consensus operational taxonomic units .
to examine the gut microbiota in further details, bion-meta data were analyzed by principal coordinates analysis .
this was conducted by applying the vegan package in r using bray-curtis distances on the untransformed sequence reads followed by k-means clustering  <cit> .
due to the hierarchical taxonomic design of the gut microbiota, where several primers in taxonomic lineage potentially target the same bacterial organism, a number of primer sets were chosen for this expanded diversity analysis : phylum: actinobacteria, bacteroidetes, fusobacteria, and spirochaetes; class: deltaproteobacteria, and epsilonproteobacteria; family: clostridium cluster i, cluster iv, and cluster xiv, and enterobacteriaceae; genus: enterococcus, lactobacillus, streptococcus.
these primer sets were selected because they provide the most comprehensive taxonomic information with the least taxonomic trespassing that would result in distortion and consequentially misinterpretation of the results.
RESULTS
from four well-managed danish conventional pig farms a total of  <dig> ileal and colonic intestinal content samples were collected from  <dig> control piglets without nnpd and  <dig> case piglets suffering from nnpd .
of the case piglets,  <dig> were treated with broad-spectrum antibiotics intramuscularly, depending on their diagnosis and according to herd treatment practices.
piglets were three to seven days old and were all reared with their biological dam and siblings  <cit> .
a bacterial taxonomic profile was obtained from both gut sections of each animal using the gut microbiotassay targeting rrna genes of major bacterial groups in the mammalian intestine.
normalized cq data were used for pca.
after excluding the primer sets with less than 50% recorded cq values , as well as the “domain bacteria a” and “domain bacteria b” primer sets ,  <dig> primer sets remained for the data analysis.
in addition,  <dig> samples were removed from the dataset as a consequence of having too few cq registrations.
from the pca scores plot , four case animals  and four control animals were randomly chosen to be representatives for their respective herds.
detailed taxonomic information was acquired for these  <dig> samples using  <dig> next generation sequencing.figure  <dig> pca score plot of ileal and colonic content of piglets with  and without  nnpd.
pca scores are generated from normalized relative values of ileal and colonic content of case and control piglets obtained from the gut microbiotassay.
asterisks symbolize samples selected for  <dig> sequencing.
circles represent k-means clustering.
for further information see additional files  <dig> and  <dig>
the gut microbiotassay
the significant findings from the gut microbiotassay are presented in table  <dig> and significant effects of interaction as investigated in formula  are presented in table  <dig>  consequently these two tables are intertwined and should be read together.
table  <dig> summarizes the estimated fold changes as the estimates of the exponentiated effect parameter, 95% confidence intervals, and p-values of these results, and table  <dig> lists the same parameters for all significant interactions.
from the relative quantitative results generated by the gut microbiotassay, it was evident that luminal content from the large intestine possessed a significantly higher bacterial load than luminal content from the small intestine , regardless of diarrhoeic status .
nnpd was associated with a diminished quantity of bacteria from the phyla actinobacteria  and firmicutes .
firmicutes comprised the following: class bacilli , genus lactobacillus , genus streptococcus , and family cl.
cluster iv .
however, genus enterococcus was estimated to be more than  <dig> times more abundant in diarrhoeic piglets .
the number of bacteria from the phylum fusobacteria  was also doubled in piglets suffering from diarrhoea .
it was investigated whether the effect of diarrhoeic status differed between gut sections through a statistical interaction term .
in both case and control piglets, there were generally more members from phylum bacteroidetes and genus bacteroides in the colon compared with the ileum .
nonetheless, case piglets had a reduced number of bacteroidetes and bacteroides compared to control piglets, and the depletion of these bacteria was located to the colon , while no significant depletion was detected in the ileum .
overall, there were greater numbers of family cl.
cluster i and species cl.
perfringens in control piglets compared with case piglets, except for cl.
perfringens, which was more numerous in the colon of case piglets, though not significantly.
the numbers of cl.
cluster i and cl.
perfringens were significantly reduced in the ileum of case piglets compared with control piglets  .
control piglets generally exhibited a disparity in the number of bacteria present in the ileum versus the colon, while the difference in bacterial abundance between gut sections was more negligible in case piglets.
it should be noted that the preceding results all relate to circumstances of ‘all other things being equal’.
independent of diarrhoeic status, piglets from first parity sows possessed significantly more bacteria from family enterobacteriaceae  and species e. coli , and fewer bacteria from phylum firmicutes , hereof genus streptococcus  compared to multiparous sows .
bacteria from family cl.
cluster xiv  and class deltaproteobacteria  were depleted, but bacteria belonging to genus streptococcus , and family cl.
cluster i  were elevated in case piglets that suffered from an increasing percentage of diarrhoea in their lifetime.
in table  <dig> the estimated fold changes may be calculated as: γpercentage points, where the parameter y from column  <dig> is the estimate of the exponential eβ of the estimated effect parameter β for diarrhoea in formula .
the diarrhoeic percentage points range from  <dig> to  <dig>
diarrhoeic piglets treated with broad-spectrum antibiotics had a reduced presence of fusobacteria , and deltaproteobacteria , .
regardless of diarrhoeic status, the quantity of bacteria in the intestine generally increased with increasing age .
the older the piglet the more bacteria from the following groups: phylum firmicutes , class bacilli , genus lactobacillus , and class cl.
cluster iv .
bacteria that diminished from the gut microbiota with increasing age were as follows: genus streptococcus , species cl.
perfringens , phylum actinobacteria , and phylum fusobacteria  .
in addition, the interaction between status and age was examined .
there was a significant difference between control piglets and case piglets over time.
as piglets aged, it was estimated that apparently increasing numbers of class beta- and gammaproteobacteria , family enterobacteriaceae  and species escherichia coli  colonized the intestine of case piglets, while the numbers of these bacteria were estimated to decrease for control piglets .
figure  <dig> illustrates the abundance of class beta- and gammaproteobacteria in ileal and colonic luminal content in relation to piglet age, segregated according to diarrhoeic status.
however, while the differences in age effects between control and case piglets were highly significant , the three estimated positive slope values for diarrhoeic piglets were not significantly greater than zero .figure  <dig> age-related changes of beta- and gammaproteobacteria in luminal content of piglets with and without nnpd.
relative amount of class beta- and gammaproteobacteria  in ileal and colonic content of all piglets included in the study, plotted as function of the age of the piglets.
case and controls are piglets with and without diarrhoea, respectively.
age-related model regression lines are depicted with corresponding 95% confidence intervals as shaded areas.
horizontal bars show means for a given age and group.
the effects of birth weight and gender were investigated, and no significant results were discovered.
table  <dig> lists the significant results calculated from the logistic regression model  investigating whether nnpd was associated with genus enterococcus, class beta- and gammaproteobacteria, family enterobacteriaceae, species e. coli, possible co-occurrence of bacteria, gut section, gilt, and age.
the model fitted the data according to standard model control measures.
nnpd was found to be significantly associated with the presence of genus enterococcus , and there was also a significant effect of the co-occurrence of this genus and species e. coli .
though the interaction apparently had a slightly diminishing effect on the probability of being recognized as an nnpd piglet, it actually contributed to the risk in the majority of cases with increasing values of e. coli.
this result is explained by the logged values of enterococci, which were negative in 89% of the data.
finally, there was a significant effect of the interaction between class beta- and gammaproteobacteria and the colon, suggesting that nnpd piglets are differentiated from control piglets by colonization of the colon with class beta- and gammaproteobacteria .
all things being equal, piglets born to gilts had  <dig> times higher odds of having nnpd compared with piglets born to multiparous sows .table  <dig>
nnpd status associated with selected bacteria and relevant variables
ln
estimate
95% confidence interval
p-
value
species e. coli
genus enterococcus
enterococcus: e. coli
of the following variables tested, genus enterococcus, class beta- and gammaproteobacteria, family enterobacteriaceae, species escherichia coli, gut section, gilt, and age, as well as possible interactions between bacteria, nnpd was found to be significantly related to the variables mentioned below.
numbers are transformed with the natural logarithm .
the odds ratio can be calculated as the exponentiated estimate.
pca
additional file  <dig> summarizes the loadings, the eigenvalues, and the variance explained by the major principal components for all pca analyses.
principal component analysis of data from the  <dig> primer sets showed no clustering in relation to pig farms .
the major loading scores in the pca were explained by the primer sets: component 1: phylum bacteroidetes  and family enterobacteriaceae , component 2: genus lactobacillus , and genus bacteroides .
generally, the control large intestine was characterized by a gut microbiota that clustered together, whereas the bacterial composition of the small intestine from control piglets was more diverse .
diarrhoeic gut sections were more scattered compared to non-diarrhoeic ones in the pca plot, demonstrating that nnpd seems to be associated with a disturbed bacterial composition and larger variation between the diarrhoeic piglets.
additional file  <dig> lists lda results.
according to lda, 70% and 83% of samples were correctly classified as cases and controls, respectively, in figure  <dig> .
when looking at the microbial composition in piglets at different ages , separation became increasingly more apparent with age, in relation to diarrhoeic status as the predictive value calculated from lda improved with age .figure  <dig> pca score plot of the gut microbiota segregated on gut section according to age in days.
pca scores from normalized relative values of ileal and colonic content of case and control piglets in relation to age in days obtained from the gut microbiotassay.
blue: control colon, purple: case colon, green: control ileum, red: case ileum.
asterisks symbolize samples selected for  <dig> sequencing.
samples cluster primarily according to diarrhoeic status with increasing age.
for further information see additional files  <dig> and  <dig>
next generation sequencing
bion-meta processing including de-multiplexing, cleaning, and chimera checking resulted in  <dig>  <dig> reads for mapping .
otus displaying similar trends to data obtained from the gut microbiotassay are included in table  <dig>  to be represented the otus have to be a somewhat consistent finding in piglets, irrespective of the variables in relation to diarrhoea, and the reads should be represented in a noteworthy number.
a shannon’s diversity index was calculated from the grand data: control ileum =  <dig> ; case ileum =  <dig> ; control colon =  <dig> ; and case colon =  <dig> , however, the group differences  are negligible and their biological relevance unclear.
additionally, a pcoa was generated from the grand data that showed separation of the sequenced samples according to nnpd status  and gut section , as shown in figure  <dig>  the two major loading scores in the pcoa were explained by the following: control: genus lactobacillus and species fusobacterium varium; case: species e. coli and genus enterococcus; ileum: genus lactobacillus and species fusobacterium ulcerans; and colon: species f. varium and species fusobacterium mortiferum .figure  <dig> pcoa score plot showing similarity of composition of the gut microbiota from sequenced samples.
grand data of ileal and colonic luminal content of a case piglet and a control piglet from each herd included in the study.
numbers refer to herd origin.
samples were randomly chosen for  <dig> sequencing based on pca scores generated from normalized cq data obtained through the gut microbiotassay for  <dig> sequencing .
circles represent k-means clustering.
for further information see additional file  <dig>
discussion
all piglets included in the study had previously been tested for hemolytic e. coli, cl.
difficile, cl.
perfringens type a and c, coronavirus, rotavirus species a and c, and microscopically inspected for parasites, with the conclusion that none of these agents were related to nnpd  <cit> .
in this study it was generally expected that bacteria associated with diarrhoea would exhibit more or less consistent tendencies in case piglets compared with control piglets.
consequently, the cq results and the sequencing results displaying similar trends will be discussed in relation to diarrhoea in the subsequent section, with the main focus being on potential nnpd-causing bacteria with high statistical estimates .
sequencing results are based on a small subset of case and control piglets.
the mapped sequences are mainly referred to as most similar species or isolate identified, while acknowledging the limits of classifying bacteria at these taxonomic levels.
nonetheless, even though mapped sequencing results at species and particularly isolate level are questionable, these are still construed as indicative and possibly important bacteria, which is why they are highlighted in the consecutive text.
genus enterococcus
enterococci were  <dig>  times more abundant in diarrhoeic piglets than in healthy piglets , and moreover genus enterococcus was one of the bacterial groups positively related with status, meaning that the presence of enterococci increases the risk of suffering from nnpd .
genus enterococcus was classified to  <dig> different otus, of which enterococcus hirae, enterococcus faecium, enterococcus durans, and enterococcus mundtii were the predominant.
enterococci are ubiquitous distributed in nature and a consistent finding in the gastrointestinal tract of several animal species, including the pig  <cit> .
an association of enterococci to neonatal porcine diarrhoea has been reported in previous studies  <cit> .
another noteworthy feature is that this genus is reported to be resistance to several antibiotics  <cit> , which could explain why the efficiency of antibiotic treatment of nnpd is variable.
larsson et al.
found neonatal porcine diarrhoea  in swedish pig farms to be associated with enteroadherent e. hirae colonizing the small intestine accompanied by mucosal lesions  <cit> .
this conforms to observations published by jonach et al.
<cit>  using fish to investigate the same piglets as the present study.
the authors demonstrated small intestinal colonization by adherent enterococcus spp.
and often found them occurring together with adherent e. coli  <cit> .
according to sequencing results, e. hirae was the most abundant species, consistent with results from a former study on enterococcal communities in pig feces  <cit> .
nonetheless, this species displayed the smallest difference in numbers of otu between diarrhoeic and healthy piglets’ gut microbiota.
however, the numerical occurrence of reads may not be related to the development of diarrhoea, especially not if diarrhoea is a result of co-occurrence of different bacteria such as e. hirae and e. coli.
e. durans could be a contributor to nnpd, as it has previously been reported to co-occur with e. coli in cases of neonatal piglet diarrhoea  <cit> .
in a study from  <dig> in which foals and gnotobiotic piglets were experimentally inoculated with e. durans, this species was found to promote the proliferation of e. coli, and to adhere to the mucosa of the small intestine in either case  <cit> .
in fact, e. durans has been hypothesized to act as a primary enteric pathogen with the ability to clear the way for other pathogens  <cit> .
although the ratio of e. durans between case and control gut sections is largest for colonic content, there were 40% more reads of this species in the ileum compared with the colon, in line with previous findings  <cit> .
the preceding information implies an interesting co-occurrence of enterococci and e. coli  <cit> .
similarly, this study also found a significant effect of the interaction between species e. coli and genus enterococcus, which contributed to the risk of piglets suffering from nnpd .
furthermore, there were generally more enterococci in the ileum compared with the colon, which supports previous findings   <cit> .
both e. hirae and e. durans have been reported to cause villous atrophy, which was frequently found in the small intestine of piglets included in this study  <cit> .
e. mundtii was almost exclusively found in the gut microbiota of case piglets, primarily in antibiotic-treated piglets, indicating that antibiotic treatment created a favorable environment for this species.
interestingly, this finding did not account for one of the antibiotic-treated animals, in which no e. mundtii was detected.
phylum fusobacteria
piglets with diarrhoea harbored more bacteria classified to phylum fusobacterium .
a similar trend has been associated with acute hemorrhagic diarrhoea in dogs and ulcerative colitis in humans  <cit> .
in this study,  <dig> otus were identified.
species f. mortiferum was a consistent finding in the digesta of all piglets, but it was less abundant in control piglets compared to case piglets.
nonetheless, portrait et al.
demonstrated that this species was able to produce bacteriocin-like substance with an inhibitory effect on a number of both gram-positive and gram-negative species  <cit> .
therefore, it is difficult to say whether the increased abundance of f. mortiferum in case piglets is due to its potential pathogenicity or if the increase is a type of defense mechanism mediated by the gut microbiota.
phylum fusobacterium was one of the two bacterial groups that were significantly reduced by antibiotic treatment .
class beta- and gammaproteobacteria, family enterobacteriaceae and species escherichia coli
e. coli are commonly found in the gastrointestinal tract.
this species is normally a harmless commensal in the host, nonetheless some e. coli possess pathogenic features  <cit> .
this study found a highly significant effect of the interaction between nnpd status and this bacterial taxonomic lineage in relation to age in days .
thus, with increasing age, these bacteria proliferated in piglets with diarrhoea , whereas members of this group were significantly reduced in control piglets.
with age, the estimated difference increased, indicating that the diarrhoeic piglets have impaired ability to control and clear this group of bacteria.
however, it also suggests that nnpd is associated with an imbalanced gut microbiota, creating a favorable environment for these bacteria.
kongsted et al.
<cit>  thoroughly examined the same animals for e. coli by aerobic culturing, serogrouping, and testing for virulence factor genes.
the most prevalent finding was the presence of non-hemolytic e. coli, independent of diarrhoeic status.
fimbrial genes were evenly distributed among the genes investigated in approximately 25% of the tested e. coli from both case and control piglets.
of all piglets, only one case of classic enterotoxigenic e. coli and six cases of non-typable hemolytic e. coli were detected, all in diarrhoeic piglets  <cit> .
nnpd piglets investigated by  <dig> sequencing from different pig farms resulted in the identification of  <dig> otus and revealed a remarkably large number of reads classified to dec8a and e. coli hq <dig> .
<dig>  both with the highest prevalence in case animals in all but one pig.
the single pig that stood out was from pig farm four, which generally differed from the others.
dec8a is the name of a diarrhoeagenic e. coli  isolate classified as ehec  <dig> , serotype o111a:nm  <cit> .
ehec  <dig> comprises a group of pathogenic e. coli reported to cause various disorders in mammals, including diarrhoea.
the o111a was not part of the serogroups investigated by kongsted et al.
<cit> .
except for one nucleotide positioned in the 5′-end of the e. coli oligonucleotide probe used for fish analysis by jonach et al.
<cit> , this probe matches the dec8a 23s rrna gene  <cit> .
hence, the possibility cannot be dismissed that the e. coli adherent to the small intestinal epithelium observed by jonach et al.
<cit>  could be dec8a.
this could also explain the frequent finding of villous atrophy previously found to characterize diarrhoeic piglets  <cit> .
e. coli hq <dig> .
<dig> isolate was detected in ileal and colonic content of all animals sequenced, though in much higher numbers in diarrhoeic piglets, indicating that it might contribute to the pathogenesis of nnpd.
piglets born to gilts had a significantly higher abundance of family enterobacteriaceae and species e. coli than piglets born to multiparous sows .
these were mainly classified to be the aforementioned isolates: dec8a and e. coli hq <dig> .
<dig>  the increased risk of having nnpd if born to a gilt is consistent with farmers’ and veterinarians’ anamnesis reports, as they describe nnpd to be most prevalent among litters from first parity sows  <cit> .
while the frequency of piglets born to gilts in the data material was 38% as a product of the randomized selection of sows, and the frequency of case piglets born to gilts was 66%, the odds ratio of nnpd for being born to a gilt from the logistic regression model  was estimated at  <dig>  , reflecting the strongly significant effect of being born to a gilt .
the logistic regression model also found that piglets of different nnpd status could be separated based on colonic colonization with species e. coli from class gammaproteobacteria according to  <dig> sequencing results.
phylum firmicutes, class bacilli, genus lactobacillus, and genus streptococcus
these bacterial groups are all part of the normal gut microbiota of pre-weaned piglets  <cit> .
the gut microbiota of case piglets was inhabited by fewer lactobacillus acidophilus, streptococcus gallolyticus subsp.
gallolyticus, streptococcus alactolyticus, and streptococcus uberis compared with the gut microbiota of control piglets.
lactobacilli have been shown to colonize the intestines of piglets soon after birth and to be a stable member of the gut microbiota throughout the intestinal tract  <cit> .
a low abundance of species l. acidophilus, which is regarded as a beneficial bacterium, could be an indicator of a troubled gastrointestinal milieu in pigs, as it was also diminished in diarrhoeic piglets included in this study.
s. gallolyticus subsp.
gallolyticus has previously been isolated from the gastrointestinal tract of numerous animal species, including pigs, and is therefore most likely included in the normal gut microbiota  <cit> .
additionally, s. alactolyticus is a member of the normal gut microbiota in pigs.
there are few reports of s. uberis being isolated from pigs  <cit> , and it does not seem to play a significant role in the development of nnpd, as the abundance of s. uberis is significantly lower in case piglets versus control piglets.
phylum firmicutes, family clostridium cluster i, species clostridium perfringens, family clostridium cluster iv, and family clostridium cluster xiv
in general, there were fewer of these bacteria in intestinal content from both the ileum and colon of case piglets compared with control piglets, except that colonic content possessed a higher number of cl.
perfringens .
it can be speculated whether the increased abundance in the colon reflects the reduction of cl.
cluster i and cl.
perfringens in the ileum .
nonetheless,  <dig> sequencing of amplicons generated by primer sets targeting cl.
cluster i and higher taxonomical levels failed to demonstrate this trend.
cl.
perfringens is a normal finding in the gastrointestinal tract but it is also a potential pathogen  <cit> .
sequencing amplicons generated by the primer set targeting family cl.
cluster i revealed reads from  <dig> otus, of which only cl.
perfringens atcc  <dig> was worthy of notice.
cl.
perfringens atcc  <dig> is a type a strain that is a potential diarrhoea-causing agent  <cit> .
however, because piglets included in this study have all been tested for cl.
perfringens type a  previously  <cit> , where it was a frequent finding with higher prevalence in control piglets versus case piglets, this species is not considered to be essential to the development of nnpd.
the fact that all remaining bacterial groups were diminished in diarrhoeic piglets is most likely due to an imbalance in the gut microbiota because all of these groups have been demonstrated in digesta from healthy piglets  <cit> .
various bacteria from family cl.
cluster iv and cl.
cluster xiv are regarded to be beneficial due to their ability to produce short chain fatty acids  such as acetate, propionate, and butyrate  <cit> .
of the  <dig> otus classified to family cl.
cluster iv, deficient bacteria mainly comprised species ruminococcus bromii, which have previously been found in the porcine gastrointestinal tract  <cit> .
a reduced number of r. bromii have also been reported in patients suffering from ulcerative colitis  <cit> .
family cl.
cluster xiv was not significantly related to the nnpd status of the pigs, but there was a significant effect of suffering from diarrhoea for a prolonged period of life  that resulted in a reduction in the abundance of this family .
thus, this bacterial group was indirectly affected by nnpd, as results indicate that family cl.
cluster xiv was vulnerable to the continuously flushing effect of diarrhoea.
a scarce population of bacteria from family cl.
cluster xiv has also been described in previous studies on gut microbial communities in intestinal disorders  <cit> .
roseburia inulinivorans was the central bacterium missing from the  <dig> otus classified to family cl.
cluster xiv.
this species is able to degrade oligofructose to free fructose, which can function as a fuel for other members of the gut microbiota, a phenomenon called cross-feeding  <cit> .
phylum bacteroidetes and genus bacteroides
members of phylum bacteroidetes are a common finding in the gastrointestinal tract of mammals, and they are also an early intestinal colonizer of the healthy piglet  <cit> .
culture studies found bacteria belonging to this phylum in piglets older than 48 hours, most frequently in the large intestine  <cit> .
this result is consistent with the findings of the present study, which found phylum bacteroidetes to be more abundant in colonic contents of both case piglets and control piglets.
however, diarrhoea resulted in a significant depletion of these bacteria in the colon .
several studies have also found a reduced presence of these bacteria in different enteric disorders, such as acute non-hemorrhagic diarrhoea in dogs, experimentally induced swine dysentery in pigs, and inflammatory bowel disease in humans  <cit> .
a total of  <dig> otus were identified.
unfortunately, sequencing results did not reveal any consistent tendencies for this taxonomic lineage.
disregarding the criterion that piglets from all herds should exhibit the same trends, case piglets harbored particularly lower numbers of species bacteroides pyogenes, bacteroides rodentium, bacteroides xylanisolvens, and the unclassified porphyromonadaceae bacterium c <dig> in their intestinal luminal content versus control piglets.
phylum actinobacteria
several studies have found phylum actinobacteria to be part of the normal gut microbiota and to be scarce in the gut microbiota of different gastrointestinal disorders, such as irritable bowel syndrome in humans and acute hemorrhagic diarrhoea in dogs  <cit> .
the phylum includes the genus bifidobacterium, which is considered to be beneficial to its host, and a number of species from this genus are recognized probiotics  <cit> .
there was a significant depletion of members from this phylum in the intestinal content of nnpd-affected piglets compared with those of control piglets .
this was supported by otu counts from all diarrhoeic piglets, except for the piglet that originated from the atypical herd mentioned earlier.
of the  <dig> otus the species scarcely represented were bifidobacterium boum, and corynebacterium kutscheri and they were mainly diminished from the colon.
bifidobacterium as probiotic has been demonstrated to protect piglets against weaning diarrhoea associated with e. coli, and b. boum to inhibit shiga toxigenic e. coli virulence gene expression experimentally  <cit> .
this fits the fact that case piglets harbored an inverse proportion of members from the taxonomical lineage class beta- and gammaproteobacteria, family enterobacteriaceae and species e. coli, and from phylum actinobacteria, having a lot more of the former bacteria and fewer of the latter.
c. kutscheri has primarily been reported to be a common bacterium on the mucous membranes of rodents including the gastrointestinal tract  <cit> .
interestingly, c. kutscheri was more plentiful in the intestinal content of healthy piglets than in that of diarrhoeic piglets, irrespective of gut section.
the bacterium does not count as part of the normal porcine gut microbiota in the existing literature.
thus, even though this cannot be ruled out for certain, it is possible that piglets may have acquired this species from foraging rodents directly or indirectly.
class deltaproteobacteria
this bacterial group was not associated with the diarrhoeic status of the piglets, but there was a significant reduction in numbers of members of this class in case piglets suffering from prolonged diarrhoea .
moreover, deltaproteobacteria were also diminished by administration of antibiotics to the piglets.
these facts are intertwined, as case piglets were mainly given antibiotics in an attempt to cure prolonged diarrhoea.
finally, both results were not strongly significant .
the only otu found by sequencing amplicons generated by the primer targeting this class was species desulfovibrio piger, and this was only identified in healthy piglets from two herds.
this species has primarily been described as an opportunistic pathogen in humans but has also been detected in the feces of wild ducks  <cit> .
while the species was primarily found in colonic content of the control piglets, it is postulated that this species constitutes part of the normal gut microbiota of piglets.
pcoa of the gut microbiota
from the bion-meta processed sequencing data a pcoa was generated .
this investigation supported the preliminary pca that was produced from the gut microbiotassay data, and additionally the major loading scores confirm the findings in the section above concerning the different bacterial groups.
it diverges from the preceding discussion by taking all bacteria into account, and the clustering strongly suggests that it is possible to associate certain bacteria with a healthy gut microbiota.
the intestinal microbiota of the aging piglet
when the data were subset according to age in days, this study found that, with age, the gut microbiota clustered in relation to status, which probably reflects the significance of e. coli in nnpd .
control piglets had significantly fewer of this species compared with case piglets, and over time the abundance of e. coli diminished in control piglets whilst it increased in case piglets .
a more defined microbiota with increasing age is in line with conclusions from a previous study  <cit> , which found that seven- to nine-day-old piglets  had the lowest individual average similarity, independent of environmental or maternal relationship, compared with piglets from older age groups.
the more distinctive clustering of nnpd status with increasing age could hint age to be a relevant factor in the diagnosis of nnpd.
antibiotics
the nnpd piglets investigated in this study included antibiotic-treated piglets.
antibiotic-treated piglets were included because they typically suffered the worst cases of nnpd and were medicated for ethical reasons.
it is evident that antibiotics have an impact on the gut microbial composition  <cit> .
however, one of the characteristics of nnpd is that it responds poorly to antibiotics.
additionally, figures  <dig>   <dig> and  <dig> illustrate how it is possible to cluster piglets according to nnpd status despite the fact that some of the nnpd cases are treated with different antibiotics, indicating that the effect of nnpd on the gut microbial distribution and composition is more significant than the effect exerted by antibiotics in piglets of this age group.
considerations on study approach
case and control piglets included in this study all originate from herds affected with nnpd, hence the study lack the inclusion of true control piglets from non-affected herds.
nonetheless, at the time of selecting the herds it was virtually impossible to find a true control herd due to the unknown aetiology and the much limited knowledge of nnpd.
this also means that despite great efforts to ensure correct inclusion of nnpd herds through predefined criteria there is still a risk that piglets actually did not suffer from nnpd.
it is acknowledged that this study only examines the bacterial composition of the luminal content, and not bacteria from the gastrointestinal wall.
however, bacteria associated with the mucosa are not restricted to this site, as they are shed into the lumen together with intestinal epithelial cells and mucus.
thus, it is assumed that this study also detects these or at least a sub-fraction of these.
additionally, results from this study are affected by the various steps involved in preparing the samples, for example dna extraction, pcr, and sequencing, which all introduce variation into the final outcome  <cit> .
furthermore, the results do not reflect bacterial numbers but relative values of 16s and 23s rrna gene abundance or otus .
conclusions drawn from sequencing data are based on otus displaying similar trends to data obtained from the gut microbiotassay.
moreover, sequences had to be represented in a noteworthy number, and with somewhat consistent findings in piglets, irrespective of the variables in relation to diarrhoea.
this means that there might be a risk that important bacteria could be overlooked if they do not meet the aforementioned criteria.
finally, to determine causality based on the associations found in this study these should be investigated further for instance through bacterial isolation and inoculation studies.
CONCLUSIONS
the results of this study on nnpd indicate that bacteria could be the aetiology of this diarrhoea that affects piglets during the first week of life.
at present, it is impossible to conclude whether diarrhoea is a consequence of an absence of beneficial bacteria or if diarrhoea is an outcome of invading pathogenic bacteria or an overgrowth of opportunistic pathogenic bacteria.
e. coli appears to be a contributing factor to nnpd as nnpd piglets differ from healthy piglets by colonic colonization of bacteria from class beta- and gammaproteobacteria, primarily being species e. coli, and it appears that there is an important effect of age, both of which might be relevant in the characterization and diagnosis of nnpd.
there is reason to believe that genus enterococcus is a participating factor, most likely species e. hirae or species e. durans.
furthermore, diarrhoeic piglets appear to suffer from an imbalanced gut microbiota, in which bacteria regarded as beneficial are diminished, which particularly accounts for l. acidophilus.
the associations found in this study are currently being investigated further.
ethical
according to danish laws no ethical approval is required for studies not including treatment groups or needle injections/blood testing.
in addition no ethical approval is required for euthanization of animals by veterinarians.
therefore this study was not subject to ethical approval, but fulfilled the regulations from the danish ministry of justice.
hence, all handling of animals was performed by trained personnel and veterinarians, and euthanization was executed exclusively by veterinarians.
procedures concerning the animals were all part of routine examinations and diagnosis of animals normally used in practice.
endnotes
<dig> gilts:  <dig> and  <dig> weeks ante partum, sows:  <dig> weeks ante partum.
<dig> etec: enterotoxigenic escherichia coli.
<dig> prrs: porcine reproductive and respiratory syndrome.
<dig> enzyme-linked immunosorbent assay.
<dig> immunoperoxidase test.
additional files
additional file 1:
primer modifications introduced to improve their performance and accompanying ribosomal database project  search results.
nucleotide explanation: m = c/a, y = t/c.
as the legends states, this table summarizes the modifications introduced in two primer sets from the gut microbiotassay: “domain bacteria b v4-v5” and “phylum firmicutes”, respectively.
ribosomal database project search results are given as percentage coverage of intended target group to demonstrate the improved performance of the modified primers tested in silico.
additional file 2:
loadings, eigenvalues, and explained variance by the major principal components related to figures 1
,
3
, and
4
, respectively.
additional file 3:
bion-meta generated output.
this file describes the bion-meta pipeline of the bioinformatics analysis of the sequencing data .
in addition, it contains a link to a website holding the bion-meta output generated in this study.
finally, it also includes a link to a website from where bion-meta can be downloaded.
additional file 4:
results from linear discriminant analysis  of data related to figure
1
and figure
3
.
linear discriminant analysis  was used to classify observations into cases or controls by their fluidigm-derived microbial profile and to estimate how well this classification corresponded with the actual status of the sample.
this classification was followed up with a multivariate analysis of variance  using wilks lambda statistic.
abbreviations
aa <dig> access array  <dig>
454blaccess array barcode library for the  <dig> gs flx titanium sequencer
cqquantification cycle
decdiarrhoeagenic e. coli
ehec 2enterohemorrhagic e. coli  <dig> clonal complex
etecenterotoxigenic e. coli
fishfluorescence in situ hybridization
lnnatural logarithm
ldalinear discriminant analysis
lsularge subunit
nnpdnew neonatal porcine diarrhoea
npdneonatal porcine diarrhoea
otuoperational taxonomic units
pcaprincipal component analysis
pcoaprincipal coordinates analysis
prrsporcine reproductive and respiratory syndrome
qpcrquantitative real-time polymerase chain reaction
rrnaribosomal rna
scfashort chain fatty acids
ssusmall subunit
competing interests
the authors declare that they have no competing interests.
authors’ contributions
mlhb, lm, mb and ks participated in conceiving and designing this study.
the inclusion of pig farms and clinical examination of animals was performed by hk.
mlhb completed all the laboratory experiments.
nl advanced bion-meta to also include search results from ribosomal lsu databases, and hci assisted in the improvement of bion-meta for this study, and critically revised the manuscript.
as conducted all statistical calculations related to the qpcr results and mls contributed with the multivariate analysis of the qpcr data and the sequencing data.
mlhb drafted the manuscript, and lm, mb and ks helped along by revising and proofreading the manuscript for appropriate content.
all authors read and approved the final manuscript.

 neonatal piglet diarrhoea is of significant importance for the pig industry because it causes economic losses due to increased morbidity and mortality, decreased weight gain, and the need for extra medications. obviously, it impairs the welfare of the animals in the short term, but it may also affect their health in the longer term as a consequence of disrupting the normal bacterial succession in the gastrointestinal tract. bacterial colonization of the mammalian gastrointestinal tract begins at birth. this colonization is a dynamic event, and the bacterial succession is influenced by a number of factors including: mode of delivery, surrounding environment, gestational age, and genetics. the colonization of the gut has a major impact on the host’s health and disease; for example, the microbiota helps in the maturation of the gastrointestinal tract and immune system, protects against pathogen colonization through competitive exclusion, and converts otherwise indigestible substances into digestible components that benefit the host. in accordance with these roles, studies of germ-free animals, deficient of normally developed gut microbiota, show animals with an impaired and immature intestinal immune system, in addition to changes in intestinal morphology. during the last ten years, increasing attention has been focused on a new neonatal porcine diarrhoea. in denmark, this diarrhoea has been named ‘new neonatal porcine diarrhoea’. what distinguishes nnpd from other types of neonatal piglet diarrhoea are the following: 1) the aetiology is unknown, however, routine diagnostic protocols show that it is not caused by known enteric pathogens such as hemolytic escherichia coli, clostridium difficile, clostridium perfringens type a or c, coronavirus, rotavirus species a or c, cryptosporidium spp., cystoisospora suis, giardia spp., or strongyloides ransomi. 2) typical strategies, such as vaccination against enterotoxigenic e. coli and cl. perfringens type c or treatment with antibiotics, do not seem to have a noteworthy effect on the diarrhoea. 3) no obvious connection between nnpd and pig farm health status or management has been demonstrated. it is difficult to estimate how widespread nnpd is, mainly because of the unknown aetiology combined with a fluctuating clinical presentation, as well as differences in routine laboratory testing, but also because of the limited number of studies focusing on this issue. nonetheless, a diarrhoea of much resemblance to nnpd has been described in sweden and france. this study is part of an interdisciplinary project investigating the aetiology of nnpd. kongsted et al. suggested the following case-definition of nnpd: “non-hemorrhagic diarrhoea during the first week of life, without detection of known infectious pathogens, characterized by milk-filled stomachs and flaccid intestines at necropsy”. the same author found dissimilarities in the course and severity of nnpd among four pig farms and estimated that affected piglets had a negative average daily weight gain with an increased risk of dying, though this risk was not significant. jonach et al. examined the role of four enteric bacterial pathogens by fluorescence in situ hybridization and found that simultaneous colonization of the intestinal mucosa with non-enterotoxigenic e. coli and enterococcus spp. could be involved in the pathogenesis of nnpd. finally, several different viral assays tested negative on samples from some of the same nnpd animals that were examined in the aforementioned studies, indicating that common known viruses likely do not contribute to nnpd. this study investigates whether nnpd is associated with the composition of the gut microbiota obtained from piglets with and without diarrhoea. this was examined by evaluating the overall bacterial composition and relative quantitative distribution of ileal and colonic intestinal content using the gut microbiotassay: an assembly of primer sets targeting ribosomal rna genes, verified to function with the high-throughput quantitative real-time pcr chip access array™ integrated fluidic circuit from fluidigm. as the name implies, the assay is designed to target major bacteria phyla and selected taxonomic sub-groups of the gut microbiota. the gut microbiotassay provides a quick overview of the distribution, as well as the relative quantity of the gut microbiota in a large number of samples simultaneously. subsequently, pcr amplicons from ileum and colon of four case piglets and four control piglets were sequenced using 454-technology to acquire deeper taxonomic information. this approach revealed diverse gut microbial profiles associated with piglet diarrhoeic status. methods animals and sample collection danish pig farms affected by nnpd were identified from conversations with veterinarians and farm managers. four pig farms that fulfilled the inclusion criteria listed in were included in the study. on each pig farm, approximately randomly chosen sows from one farrowing batch were followed for a seven day period following farrowing. all newborn piglets were weighed at the beginning of the trial, and piglets weighing less than 800 g were excluded. all animals were subject to a daily clinical examination that paid special attention to fecal appearance on rectal swabs. diarrhoea was defined as loose or watery feces. based on these observations, piglets were characterized as either cases or controls: a case piglet had suffered from diarrhoea for at least two consecutive days, whereas a control piglet had never experienced diarrhoea. the inclusion criterion applied to case piglets resulted in a greater number of case piglets from gilts, as piglets born to multiparous sows did not suffer from diarrhoea to the same extent as piglets born to gilts. control piglets were as far as possible collected from litters without diarrhoea. in total control piglets and case piglets were selected and brought to the danish pig research centre, kjellerup, for euthanization and necropsy. for ethical reasons, farmers were allowed to treat piglets for any disorder if necessary, on the condition that any such treatments were recorded for individual animals. during necropsy, intestinal contents were collected from the distal small intestine and the large intestine of each animal and stored at −80°c until further analysis.table piglets included in the study age in days 3 4 5 6 7 total gilts sows diarrhoea − + − + − + − + − + - + - + pig farm 1 0 4 13 9 pig farm 2 2 12 10 0 pig farm 3 0 10 13 4 pig farm 4 4 8 8 5 total 14 15 9 10 17 17 6 6 4 4 6 34 44 18 the number of piglets is listed for each pig farm in accordance with age in days: control piglets without diarrhoea, case piglets with diarrhoea for ≥2 days. the column “total” summarizes the number of piglets born to gilts or sows by farm. pig farm inclusion criteria for studying nnpd:a stock of more than sows routinely vaccinate against etec and cl. perfringens type c farrowing units proven to be prrs virus negative as demonstrated in blood samples tested by elisa4/ipt or pcr. nnpd has been a problem for a period of minimum 3– months management-related diarrhoea has been excluded traditional interventions such as vaccines and antibiotics have limited effect on the diarrhoea at least 30% of the litters are affected the piglets suffer from diarrhoea during the first week of life random check of five case piglets should test negative for etec, cl. perfringens type c and rotavirus by routine diagnostic examination dna extraction based on development work comparing various extraction methods, the maxwell® lev blood dna kit was determined to be the optimal method, based on providing sufficient dna yields and acceptable purities. according to manufacturer’s recommendations, total dna was extracted from intestinal content using this kit. a total of 200 mg of intestinal contents were suspended in 600 μl pbs and vortexed until visually homogeneous. the samples were centrifuged for 2 min at 200 × g, and the supernatant was transferred to new tubes. a volume of 350 μl of lysis buffer was then added, and the bacterial cells were lysed by bead beading for  min at hertz with a 5 mm steel bead. subsequently, samples were centrifuged for 1 min at 1000 × g at 4°c, after which the supernatant was transferred to novel 2 ml tubes and mixed manually with 30 μl proteinase k. next, samples were incubated for 30 min at 56°c, and centrifuged at 000 × g for 1 min. the entire suspension was transferred to the sample inlets on the cartridge, and 50 μl of elution buffer was added to the bottom of the elution tubes. the cartridges were prepared according to the manufacturer’s instructions and the settings “research mode”, “lev mode”, “dna”, and “blood/cell” were selected for dna extraction using the maxwell® instrument. finally, tubes were centrifuged at × g for 3 min to settle any magnetic bead leftovers, and the dna was moved to new tubes. dna concentration and purity was measured using the nanodrop® nd- spectrophotometer, and dna was stored at −20°c until further analysis. analysis of bacterial profiles with use of the gut microbiotassay to obtain a bacterial profile of the intestinal content, this study used the gut microbiotassay designed for the aa with a few modifications. in order to improve the performance of the primer sets “domain bacteria b v4-v5” and “phylum firmicutes” degenerate nucleotides were introduced: nucleotide was changed from c to m, and nucleotide from t to y,, in the forward or reverse primer, respectively. all remaining primer sets were identical to those previously published as ‘the gut microbiotassay’. primers were purchased from eurofins mwg synthesis gmbh and stored at −20°c. in brief, the aa is a qpcr platform that is capable of running 48 × 48 =  individual reactions simultaneously, and it enables quick and easy library preparation for sequencing. the aa was processed and prepared following the ‘access array system™ user guide’, with and without adding the access array barcode library for the gs flx titanium sequencer, as previously described. primers targeting bacteria at the taxonomic level of species were not tagged or included for sequencing, as these amplicons were regarded as contributing little information due to their specificity. however, analysis at the species-level was possible using the information obtained at higher taxonomic levels. all samples were diluted to 50 ng/μl with nuclease-free water. primers were diluted to 4 μm with 20× access array loading buffer and nuclease-free water. master mix was a mixture of: 10× faststart high fidelity reaction buffer with 18 mm mgcl, 25 mm mgcl, dmso, 20× access array loading reagent, 50× rox, 20× evagreen®, 10 mm pcr grade nucleotide mix, and u/μl faststart high fidelity enzyme blend, in final concentrations of 1×,  mm, 5%, 1×, ×, 1×, 200 μm, and u/μl, respectively. sample mix was prepared from 3 μl master mix, 1 μl 454bl, and 1 μl dna as the very last step before running the array ). when the pcr reaction had finished, the barcoded pcr amplicons were harvested and stored at −20°c. next generation sequencing samples from two representative animals were selected from each farm: one piglet with nnpd and one piglet without nnpd, giving a total of samples to be sequenced. piglets were chosen to be approximately the same age. concentrations of the respective pcr amplicons for each animal were determined using the agilent chip. amplicons were subsequently pooled in equal concentrations and size-separated by running the amplicons for 86 min, 90 v, in a % seakem® le agarose gel followed by incubation for 30 min in % ethidium bromide for staining. by means of uv radiation from the bio-rad universal hood ii gel bands were visualized and bands spanning the size range of the primer products were excised. finally, the qiaquick gel extraction kit was used to purify dna from the gel. this pool of  ng dna derived from ileal and colonic luminal contents of different animals was sequenced on a half picotiterplate™ by a gs flx titanium sequencer via lgc genomics. data analysis relative quantification of quantification cycle values data analysis was conducted as described in a previously published methodology article: raw cq values were exported from ‘fluidigm real-time pcr analysis’ software version. to microsoft excel. to even out possible variation between the aa runs, all cq values were normalized to an interplate calibrator. next, cq values from each sample were normalized to their respective “domain bacteria b v4-v5” primer set, thereby calculating the relative quantification. principal component analysis on normalized cq data principal component analysis was conducted in the software package r, version.. all normalized data were initially transformed with the natural logarithm. next, primer sets with less than 50% recorded cq values were excluded from the analysis, and in the same manner, samples that resulted in less than 50% registered cq values were removed. for the remaining primer sets, missing values were substituted with the lowest registered primer-specific value. data were scaled by the individual primer mean and standard deviation for pca. the number of components to be investigated further was determined by visual inspection of scree plots of eigenvalues; based on an inflection point or when the components explained > 60% variance. all possible combinations of included principal components were visualized to ensure that no important clustering was overlooked. the number of components used for graphical representation in the pca was set at two to ease interpretation. samples used for sequencing were pinpointed from the pca results. in order for samples to be considered for sequencing the samples were represented with luminal content from both ileum and colon, and none of these were allowed to have outlying coordinates in the pca, but were somewhat in the centre of its respective category. linear discriminant analysis based on leave-one-out cross-validation was used to classify observations into cases or controls by their fluidigm-derived microbial profile and to estimate how well this classification corresponded with the actual status of the sample. this classification was followed up with a multivariate analysis of variance using wilks lambda statistic. statistics on the gut microbiotassay bacterial profiles obtained from the gut microbiotassay were used to study the gut microbiota of piglets with and without nnpd. for the statistical calculations all values were initially ln-transformed. a linear mixed-effect model was then applied to each primer set separately with randomized effect of herd origin: lnyi=agutsectioni,statusi+agilti+atreatmenti+βstatusi⋅agei+β⋅diarrhoeai+η⋅yherdi+εi,i= …,, where status indicates the status of the i’th animal, and where yherd is a normalized gaussian stochastic variable indicating the randomized effect of herd, indexed by the herd of the i’th animal, so that animals from the same herd share the random effect. the model is applied for different responses y, as listed in the first column in table  the following variables were included as deterministic effects in the model: gut section, status, gilt, percentage of diarrhoea, treatment, age. the interaction between gut section and status was also included to examine if possible diarrhoea-causing agents could be traced to a specific gut section. additionally, the interaction between status and age was included to study whether the effect of age differed between case and control piglets. for each primer set as response variable, the model was reduced by backwards stepwise elimination, resulting in a series of primer-specific final models. model reduction was performed as likelihood ratio tests. model fit for each primer set was assessed through residual analysis of the model residuals corrected for random effects and the estimated random effects themselves, using standard techniques. these analyses were consistent with standard model behavior. the effect of gender and birth weight in relation to status was also evaluated. significances in the final models were reported through p-values.table significant findings from the gut microbiotassay and sequencing results displaying similar trends primer set and respective sequencing findings gut section estimate gut section p- value status estimate status p -value gilt estimate gilt p -value percentage of diarrhoea estimate # percentage of diarrhoea p -value age estimate age p -value antibiotics estimate antibiotics p -value domain bacteria b 268 phylum firmicutes 183 class bacilli 125 genus enterococcus 27 genus streptococcus 47 genus lactobacillus 79 family clostridium cluster i 17 species cl. perfringens family clostridium cluster iv 64 family clostridium cluster xiv 81 phylum bacteroidetes 126 genus bacteroides 65 species bacteroides pyogenes † species bacteroides rodentium † species bacteroides xylanisolvens † species porphyromonadaceae bacterium c † phylum actinobacteria 74 family bifidobacteriaceae 12 phylum fusobacteria 16 class beta- and gammaproteobacteria 109 family enterobacteriaceae 89 species e. coli class deltaproteobacteria 1 species desulfovibrio piger † estimated fold changes as the estimates of the exponentiated effect parameter, corresponding 95% confidence intervals, and p-values. superscripted numbers equal the total number of otus belonging to the respective taxonomic group identified by sequencing amplicons generated by the gut microbiotassay. otu-ratio between diarrhoeic and healthy piglets for ileum and colon, respectively, listed in parenthesis. †indicate a significant effect of interaction. species mark with † are not consistent findings. #estimated fold changes for this category are calculated as: γ percentage points, where the parameter y from column is the estimate of the exponential eβ of the estimated effect parameter β for diarrhoea in formula. the diarrhoeic percentage points range from to previous studies have suggested that simultaneous colonization with e. coli and enterococci species may be a cause of neonatal porcine diarrhoea. this co-occurrence was also implied using the aforementioned model ), thus, this was further investigated. a logistic regression model with randomized effect of herd origin, fit by penalized quasi-likelihood, was used to determine if the odds of nnpd was associated with the following variables: gut section, gilt, age, genus enterococcus, class beta- and gammaproteobacteria, family enterobacteriaceae, species e. coli, as well as possible interaction in between bacteria, and between bacteria and gut section:table estimated fold changes, 95% confidence intervals, and p- values of significant interactions from the gut microbiotassay primer set gut section: status estimate gut section: status estimate gut section: status p- value status: age estimate, day 3 status: age estimate, day 4 status: age estimate, day 5 status: age estimate, day 6 status: age estimate, day 7 status: age estimate, mean age status: age p -value family clostridium cluster i species cl. perfringens phylum bacteroidetes genus bacteroides class beta- and gammaproteobacteria family enterobacteriaceae species e. coli this table is an extension of table  and the estimated fold changes are the estimates of the exponentiated effect parameter. logitpi=agutsectioni+agilti+βage⋅agei+∑j=14βj,gutsectionibacteriaj,i+∑k=14∑ℓ≠kβk,lbacteriak,ibacteriaℓ,i+ηyherdi, where bacteria1-bacteria are genus enterococcus, class beta- and gammaproteobacteria, family enterobacteriaceae, and species e. coli, respectively, where bacteriaj,i is the value of bacteriaj animal i, where agei is the age in days of animal i, and where yherd is the random effect component corresponding to the herd of animal i. the models were validated through graphical procedures. a p-value of less than was considered significant. statistical analyses were performed in r. bioinformatics analysis of sequencing data sequencing data available at ncbi sequence read archive were analyzed using the as yet unpublished open source package bion-meta. this program facilitates quick and easy bioinformatics analysis, and bion-meta has previously been applied to a similar dataset. however, in addition to matching all sequences against the ribosomal small subunit silva dataset, suited for 16s rrna gene targeting primer sets, bion-meta was advanced to encompass searches in the ribosomal large subunit silva dataset, thereby accommodating the primer sets targeting the 23s rrna gene used in the gut microbiotassay. the bion-meta workflow included the following elements: 1) separation by sample barcodes and primer tags; 2) removal of primer remnants and bases at the ends with lower than 96% quality, as well as sequence filtering by length and filtering by minimum base quality of 96%.; 3) removal of chimeric sequences; 4) separation of each sample by matching it with the phylogenetic primer; 5) matching all sequences against the ssu and lsu silva datasets and producing a table with the highest 1% similarities for each query; and 6) mapping the similarities to the silva ssu and lsu taxonomies, identifying consensus operational taxonomic units. to examine the gut microbiota in further details, bion-meta data were analyzed by principal coordinates analysis. this was conducted by applying the vegan package in r using bray-curtis distances on the untransformed sequence reads followed by k-means clustering. due to the hierarchical taxonomic design of the gut microbiota, where several primers in taxonomic lineage potentially target the same bacterial organism, a number of primer sets were chosen for this expanded diversity analysis: phylum: actinobacteria, bacteroidetes, fusobacteria, and spirochaetes; class: deltaproteobacteria, and epsilonproteobacteria; family: clostridium cluster i, cluster iv, and cluster xiv, and enterobacteriaceae; genus: enterococcus, lactobacillus, streptococcus. these primer sets were selected because they provide the most comprehensive taxonomic information with the least taxonomic trespassing that would result in distortion and consequentially misinterpretation of the results.  from four well-managed danish conventional pig farms a total of ileal and colonic intestinal content samples were collected from control piglets without nnpd and case piglets suffering from nnpd. of the case piglets, were treated with broad-spectrum antibiotics intramuscularly, depending on their diagnosis and according to herd treatment practices. piglets were three to seven days old and were all reared with their biological dam and siblings. a bacterial taxonomic profile was obtained from both gut sections of each animal using the gut microbiotassay targeting rrna genes of major bacterial groups in the mammalian intestine. normalized cq data were used for pca. after excluding the primer sets with less than 50% recorded cq values, as well as the “domain bacteria a” and “domain bacteria b” primer sets, primer sets remained for the data analysis. in addition, samples were removed from the dataset as a consequence of having too few cq registrations. from the pca scores plot, four case animals and four control animals were randomly chosen to be representatives for their respective herds. detailed taxonomic information was acquired for these samples using next generation sequencing.figure pca score plot of ileal and colonic content of piglets with and without nnpd. pca scores are generated from normalized relative values of ileal and colonic content of case and control piglets obtained from the gut microbiotassay. asterisks symbolize samples selected for sequencing. circles represent k-means clustering. for further information see additional files and the gut microbiotassay the significant findings from the gut microbiotassay are presented in table  and significant effects of interaction as investigated in formula are presented in table  consequently these two tables are intertwined and should be read together. table  summarizes the estimated fold changes as the estimates of the exponentiated effect parameter, 95% confidence intervals, and p-values of these results, and table  lists the same parameters for all significant interactions. from the relative quantitative results generated by the gut microbiotassay, it was evident that luminal content from the large intestine possessed a significantly higher bacterial load than luminal content from the small intestine, regardless of diarrhoeic status. nnpd was associated with a diminished quantity of bacteria from the phyla actinobacteria and firmicutes. firmicutes comprised the following: class bacilli, genus lactobacillus, genus streptococcus, and family cl. cluster iv. however, genus enterococcus was estimated to be more than times more abundant in diarrhoeic piglets. the number of bacteria from the phylum fusobacteria was also doubled in piglets suffering from diarrhoea. it was investigated whether the effect of diarrhoeic status differed between gut sections through a statistical interaction term. in both case and control piglets, there were generally more members from phylum bacteroidetes and genus bacteroides in the colon compared with the ileum. nonetheless, case piglets had a reduced number of bacteroidetes and bacteroides compared to control piglets, and the depletion of these bacteria was located to the colon, while no significant depletion was detected in the ileum. overall, there were greater numbers of family cl. cluster i and species cl. perfringens in control piglets compared with case piglets, except for cl. perfringens, which was more numerous in the colon of case piglets, though not significantly. the numbers of cl. cluster i and cl. perfringens were significantly reduced in the ileum of case piglets compared with control piglets. control piglets generally exhibited a disparity in the number of bacteria present in the ileum versus the colon, while the difference in bacterial abundance between gut sections was more negligible in case piglets. it should be noted that the preceding results all relate to circumstances of ‘all other things being equal’. independent of diarrhoeic status, piglets from first parity sows possessed significantly more bacteria from family enterobacteriaceae and species e. coli, and fewer bacteria from phylum firmicutes, hereof genus streptococcus compared to multiparous sows. bacteria from family cl. cluster xiv and class deltaproteobacteria were depleted, but bacteria belonging to genus streptococcus, and family cl. cluster i were elevated in case piglets that suffered from an increasing percentage of diarrhoea in their lifetime. in table  the estimated fold changes may be calculated as: γpercentage points, where the parameter y from column is the estimate of the exponential eβ of the estimated effect parameter β for diarrhoea in formula. the diarrhoeic percentage points range from to diarrhoeic piglets treated with broad-spectrum antibiotics had a reduced presence of fusobacteria, and deltaproteobacteria,. regardless of diarrhoeic status, the quantity of bacteria in the intestine generally increased with increasing age. the older the piglet the more bacteria from the following groups: phylum firmicutes, class bacilli, genus lactobacillus, and class cl. cluster iv. bacteria that diminished from the gut microbiota with increasing age were as follows: genus streptococcus, species cl. perfringens, phylum actinobacteria, and phylum fusobacteria. in addition, the interaction between status and age was examined. there was a significant difference between control piglets and case piglets over time. as piglets aged, it was estimated that apparently increasing numbers of class beta- and gammaproteobacteria, family enterobacteriaceae and species escherichia coli colonized the intestine of case piglets, while the numbers of these bacteria were estimated to decrease for control piglets. figure  illustrates the abundance of class beta- and gammaproteobacteria in ileal and colonic luminal content in relation to piglet age, segregated according to diarrhoeic status. however, while the differences in age effects between control and case piglets were highly significant, the three estimated positive slope values for diarrhoeic piglets were not significantly greater than zero.figure age-related changes of beta- and gammaproteobacteria in luminal content of piglets with and without nnpd. relative amount of class beta- and gammaproteobacteria in ileal and colonic content of all piglets included in the study, plotted as function of the age of the piglets. case and controls are piglets with and without diarrhoea, respectively. age-related model regression lines are depicted with corresponding 95% confidence intervals as shaded areas. horizontal bars show means for a given age and group. the effects of birth weight and gender were investigated, and no significant results were discovered. table  lists the significant results calculated from the logistic regression model investigating whether nnpd was associated with genus enterococcus, class beta- and gammaproteobacteria, family enterobacteriaceae, species e. coli, possible co-occurrence of bacteria, gut section, gilt, and age. the model fitted the data according to standard model control measures. nnpd was found to be significantly associated with the presence of genus enterococcus, and there was also a significant effect of the co-occurrence of this genus and species e. coli. though the interaction apparently had a slightly diminishing effect on the probability of being recognized as an nnpd piglet, it actually contributed to the risk in the majority of cases with increasing values of e. coli. this result is explained by the logged values of enterococci, which were negative in 89% of the data. finally, there was a significant effect of the interaction between class beta- and gammaproteobacteria and the colon, suggesting that nnpd piglets are differentiated from control piglets by colonization of the colon with class beta- and gammaproteobacteria. all things being equal, piglets born to gilts had times higher odds of having nnpd compared with piglets born to multiparous sows.table nnpd status associated with selected bacteria and relevant variables ln estimate 95% confidence interval p- value species e. coli genus enterococcus enterococcus: e. coli of the following variables tested, genus enterococcus, class beta- and gammaproteobacteria, family enterobacteriaceae, species escherichia coli, gut section, gilt, and age, as well as possible interactions between bacteria, nnpd was found to be significantly related to the variables mentioned below. numbers are transformed with the natural logarithm. the odds ratio can be calculated as the exponentiated estimate. pca additional file summarizes the loadings, the eigenvalues, and the variance explained by the major principal components for all pca analyses. principal component analysis of data from the primer sets showed no clustering in relation to pig farms. the major loading scores in the pca were explained by the primer sets: component 1: phylum bacteroidetes and family enterobacteriaceae, component 2: genus lactobacillus, and genus bacteroides. generally, the control large intestine was characterized by a gut microbiota that clustered together, whereas the bacterial composition of the small intestine from control piglets was more diverse. diarrhoeic gut sections were more scattered compared to non-diarrhoeic ones in the pca plot, demonstrating that nnpd seems to be associated with a disturbed bacterial composition and larger variation between the diarrhoeic piglets. additional file lists lda results. according to lda, 70% and 83% of samples were correctly classified as cases and controls, respectively, in figure . when looking at the microbial composition in piglets at different ages, separation became increasingly more apparent with age, in relation to diarrhoeic status as the predictive value calculated from lda improved with age.figure pca score plot of the gut microbiota segregated on gut section according to age in days. pca scores from normalized relative values of ileal and colonic content of case and control piglets in relation to age in days obtained from the gut microbiotassay. blue: control colon, purple: case colon, green: control ileum, red: case ileum. asterisks symbolize samples selected for sequencing. samples cluster primarily according to diarrhoeic status with increasing age. for further information see additional files and next generation sequencing bion-meta processing including de-multiplexing, cleaning, and chimera checking resulted in reads for mapping. otus displaying similar trends to data obtained from the gut microbiotassay are included in table  to be represented the otus have to be a somewhat consistent finding in piglets, irrespective of the variables in relation to diarrhoea, and the reads should be represented in a noteworthy number. a shannon’s diversity index was calculated from the grand data: control ileum = ; case ileum = ; control colon = ; and case colon = , however, the group differences are negligible and their biological relevance unclear. additionally, a pcoa was generated from the grand data that showed separation of the sequenced samples according to nnpd status and gut section, as shown in figure  the two major loading scores in the pcoa were explained by the following: control: genus lactobacillus and species fusobacterium varium; case: species e. coli and genus enterococcus; ileum: genus lactobacillus and species fusobacterium ulcerans; and colon: species f. varium and species fusobacterium mortiferum.figure pcoa score plot showing similarity of composition of the gut microbiota from sequenced samples. grand data of ileal and colonic luminal content of a case piglet and a control piglet from each herd included in the study. numbers refer to herd origin. samples were randomly chosen for sequencing based on pca scores generated from normalized cq data obtained through the gut microbiotassay for sequencing. circles represent k-means clustering. for further information see additional file discussion all piglets included in the study had previously been tested for hemolytic e. coli, cl. difficile, cl. perfringens type a and c, coronavirus, rotavirus species a and c, and microscopically inspected for parasites, with the conclusion that none of these agents were related to nnpd. in this study it was generally expected that bacteria associated with diarrhoea would exhibit more or less consistent tendencies in case piglets compared with control piglets. consequently, the cq results and the sequencing results displaying similar trends will be discussed in relation to diarrhoea in the subsequent section, with the main focus being on potential nnpd-causing bacteria with high statistical estimates. sequencing results are based on a small subset of case and control piglets. the mapped sequences are mainly referred to as most similar species or isolate identified, while acknowledging the limits of classifying bacteria at these taxonomic levels. nonetheless, even though mapped sequencing results at species and particularly isolate level are questionable, these are still construed as indicative and possibly important bacteria, which is why they are highlighted in the consecutive text. genus enterococcus enterococci were times more abundant in diarrhoeic piglets than in healthy piglets, and moreover genus enterococcus was one of the bacterial groups positively related with status, meaning that the presence of enterococci increases the risk of suffering from nnpd. genus enterococcus was classified to different otus, of which enterococcus hirae, enterococcus faecium, enterococcus durans, and enterococcus mundtii were the predominant. enterococci are ubiquitous distributed in nature and a consistent finding in the gastrointestinal tract of several animal species, including the pig. an association of enterococci to neonatal porcine diarrhoea has been reported in previous studies. another noteworthy feature is that this genus is reported to be resistance to several antibiotics, which could explain why the efficiency of antibiotic treatment of nnpd is variable. larsson et al. found neonatal porcine diarrhoea in swedish pig farms to be associated with enteroadherent e. hirae colonizing the small intestine accompanied by mucosal lesions. this conforms to observations published by jonach et al. using fish to investigate the same piglets as the present study. the authors demonstrated small intestinal colonization by adherent enterococcus spp. and often found them occurring together with adherent e. coli. according to sequencing results, e. hirae was the most abundant species, consistent with results from a former study on enterococcal communities in pig feces. nonetheless, this species displayed the smallest difference in numbers of otu between diarrhoeic and healthy piglets’ gut microbiota. however, the numerical occurrence of reads may not be related to the development of diarrhoea, especially not if diarrhoea is a result of co-occurrence of different bacteria such as e. hirae and e. coli. e. durans could be a contributor to nnpd, as it has previously been reported to co-occur with e. coli in cases of neonatal piglet diarrhoea. in a study from in which foals and gnotobiotic piglets were experimentally inoculated with e. durans, this species was found to promote the proliferation of e. coli, and to adhere to the mucosa of the small intestine in either case. in fact, e. durans has been hypothesized to act as a primary enteric pathogen with the ability to clear the way for other pathogens. although the ratio of e. durans between case and control gut sections is largest for colonic content, there were 40% more reads of this species in the ileum compared with the colon, in line with previous findings. the preceding information implies an interesting co-occurrence of enterococci and e. coli. similarly, this study also found a significant effect of the interaction between species e. coli and genus enterococcus, which contributed to the risk of piglets suffering from nnpd. furthermore, there were generally more enterococci in the ileum compared with the colon, which supports previous findings. both e. hirae and e. durans have been reported to cause villous atrophy, which was frequently found in the small intestine of piglets included in this study. e. mundtii was almost exclusively found in the gut microbiota of case piglets, primarily in antibiotic-treated piglets, indicating that antibiotic treatment created a favorable environment for this species. interestingly, this finding did not account for one of the antibiotic-treated animals, in which no e. mundtii was detected. phylum fusobacteria piglets with diarrhoea harbored more bacteria classified to phylum fusobacterium. a similar trend has been associated with acute hemorrhagic diarrhoea in dogs and ulcerative colitis in humans. in this study, otus were identified. species f. mortiferum was a consistent finding in the digesta of all piglets, but it was less abundant in control piglets compared to case piglets. nonetheless, portrait et al. demonstrated that this species was able to produce bacteriocin-like substance with an inhibitory effect on a number of both gram-positive and gram-negative species. therefore, it is difficult to say whether the increased abundance of f. mortiferum in case piglets is due to its potential pathogenicity or if the increase is a type of defense mechanism mediated by the gut microbiota. phylum fusobacterium was one of the two bacterial groups that were significantly reduced by antibiotic treatment. class beta- and gammaproteobacteria, family enterobacteriaceae and species escherichia coli e. coli are commonly found in the gastrointestinal tract. this species is normally a harmless commensal in the host, nonetheless some e. coli possess pathogenic features. this study found a highly significant effect of the interaction between nnpd status and this bacterial taxonomic lineage in relation to age in days. thus, with increasing age, these bacteria proliferated in piglets with diarrhoea, whereas members of this group were significantly reduced in control piglets. with age, the estimated difference increased, indicating that the diarrhoeic piglets have impaired ability to control and clear this group of bacteria. however, it also suggests that nnpd is associated with an imbalanced gut microbiota, creating a favorable environment for these bacteria. kongsted et al. thoroughly examined the same animals for e. coli by aerobic culturing, serogrouping, and testing for virulence factor genes. the most prevalent finding was the presence of non-hemolytic e. coli, independent of diarrhoeic status. fimbrial genes were evenly distributed among the genes investigated in approximately 25% of the tested e. coli from both case and control piglets. of all piglets, only one case of classic enterotoxigenic e. coli and six cases of non-typable hemolytic e. coli were detected, all in diarrhoeic piglets. nnpd piglets investigated by sequencing from different pig farms resulted in the identification of otus and revealed a remarkably large number of reads classified to dec8a and e. coli hq. both with the highest prevalence in case animals in all but one pig. the single pig that stood out was from pig farm four, which generally differed from the others. dec8a is the name of a diarrhoeagenic e. coli isolate classified as ehec, serotype o111a:nm. ehec comprises a group of pathogenic e. coli reported to cause various disorders in mammals, including diarrhoea. the o111a was not part of the serogroups investigated by kongsted et al.. except for one nucleotide positioned in the 5′-end of the e. coli oligonucleotide probe used for fish analysis by jonach et al., this probe matches the dec8a 23s rrna gene. hence, the possibility cannot be dismissed that the e. coli adherent to the small intestinal epithelium observed by jonach et al. could be dec8a. this could also explain the frequent finding of villous atrophy previously found to characterize diarrhoeic piglets. e. coli hq. isolate was detected in ileal and colonic content of all animals sequenced, though in much higher numbers in diarrhoeic piglets, indicating that it might contribute to the pathogenesis of nnpd. piglets born to gilts had a significantly higher abundance of family enterobacteriaceae and species e. coli than piglets born to multiparous sows. these were mainly classified to be the aforementioned isolates: dec8a and e. coli hq. the increased risk of having nnpd if born to a gilt is consistent with farmers’ and veterinarians’ anamnesis reports, as they describe nnpd to be most prevalent among litters from first parity sows. while the frequency of piglets born to gilts in the data material was 38% as a product of the randomized selection of sows, and the frequency of case piglets born to gilts was 66%, the odds ratio of nnpd for being born to a gilt from the logistic regression model was estimated at, reflecting the strongly significant effect of being born to a gilt. the logistic regression model also found that piglets of different nnpd status could be separated based on colonic colonization with species e. coli from class gammaproteobacteria according to sequencing results. phylum firmicutes, class bacilli, genus lactobacillus, and genus streptococcus these bacterial groups are all part of the normal gut microbiota of pre-weaned piglets. the gut microbiota of case piglets was inhabited by fewer lactobacillus acidophilus, streptococcus gallolyticus subsp. gallolyticus, streptococcus alactolyticus, and streptococcus uberis compared with the gut microbiota of control piglets. lactobacilli have been shown to colonize the intestines of piglets soon after birth and to be a stable member of the gut microbiota throughout the intestinal tract. a low abundance of species l. acidophilus, which is regarded as a beneficial bacterium, could be an indicator of a troubled gastrointestinal milieu in pigs, as it was also diminished in diarrhoeic piglets included in this study. s. gallolyticus subsp. gallolyticus has previously been isolated from the gastrointestinal tract of numerous animal species, including pigs, and is therefore most likely included in the normal gut microbiota. additionally, s. alactolyticus is a member of the normal gut microbiota in pigs. there are few reports of s. uberis being isolated from pigs, and it does not seem to play a significant role in the development of nnpd, as the abundance of s. uberis is significantly lower in case piglets versus control piglets. phylum firmicutes, family clostridium cluster i, species clostridium perfringens, family clostridium cluster iv, and family clostridium cluster xiv in general, there were fewer of these bacteria in intestinal content from both the ileum and colon of case piglets compared with control piglets, except that colonic content possessed a higher number of cl. perfringens. it can be speculated whether the increased abundance in the colon reflects the reduction of cl. cluster i and cl. perfringens in the ileum. nonetheless, sequencing of amplicons generated by primer sets targeting cl. cluster i and higher taxonomical levels failed to demonstrate this trend. cl. perfringens is a normal finding in the gastrointestinal tract but it is also a potential pathogen. sequencing amplicons generated by the primer set targeting family cl. cluster i revealed reads from otus, of which only cl. perfringens atcc was worthy of notice. cl. perfringens atcc is a type a strain that is a potential diarrhoea-causing agent. however, because piglets included in this study have all been tested for cl. perfringens type a previously, where it was a frequent finding with higher prevalence in control piglets versus case piglets, this species is not considered to be essential to the development of nnpd. the fact that all remaining bacterial groups were diminished in diarrhoeic piglets is most likely due to an imbalance in the gut microbiota because all of these groups have been demonstrated in digesta from healthy piglets. various bacteria from family cl. cluster iv and cl. cluster xiv are regarded to be beneficial due to their ability to produce short chain fatty acids such as acetate, propionate, and butyrate. of the otus classified to family cl. cluster iv, deficient bacteria mainly comprised species ruminococcus bromii, which have previously been found in the porcine gastrointestinal tract. a reduced number of r. bromii have also been reported in patients suffering from ulcerative colitis. family cl. cluster xiv was not significantly related to the nnpd status of the pigs, but there was a significant effect of suffering from diarrhoea for a prolonged period of life that resulted in a reduction in the abundance of this family. thus, this bacterial group was indirectly affected by nnpd, as results indicate that family cl. cluster xiv was vulnerable to the continuously flushing effect of diarrhoea. a scarce population of bacteria from family cl. cluster xiv has also been described in previous studies on gut microbial communities in intestinal disorders. roseburia inulinivorans was the central bacterium missing from the otus classified to family cl. cluster xiv. this species is able to degrade oligofructose to free fructose, which can function as a fuel for other members of the gut microbiota, a phenomenon called cross-feeding. phylum bacteroidetes and genus bacteroides members of phylum bacteroidetes are a common finding in the gastrointestinal tract of mammals, and they are also an early intestinal colonizer of the healthy piglet. culture studies found bacteria belonging to this phylum in piglets older than 48 hours, most frequently in the large intestine. this result is consistent with the findings of the present study, which found phylum bacteroidetes to be more abundant in colonic contents of both case piglets and control piglets. however, diarrhoea resulted in a significant depletion of these bacteria in the colon. several studies have also found a reduced presence of these bacteria in different enteric disorders, such as acute non-hemorrhagic diarrhoea in dogs, experimentally induced swine dysentery in pigs, and inflammatory bowel disease in humans. a total of otus were identified. unfortunately, sequencing results did not reveal any consistent tendencies for this taxonomic lineage. disregarding the criterion that piglets from all herds should exhibit the same trends, case piglets harbored particularly lower numbers of species bacteroides pyogenes, bacteroides rodentium, bacteroides xylanisolvens, and the unclassified porphyromonadaceae bacterium c in their intestinal luminal content versus control piglets. phylum actinobacteria several studies have found phylum actinobacteria to be part of the normal gut microbiota and to be scarce in the gut microbiota of different gastrointestinal disorders, such as irritable bowel syndrome in humans and acute hemorrhagic diarrhoea in dogs. the phylum includes the genus bifidobacterium, which is considered to be beneficial to its host, and a number of species from this genus are recognized probiotics. there was a significant depletion of members from this phylum in the intestinal content of nnpd-affected piglets compared with those of control piglets. this was supported by otu counts from all diarrhoeic piglets, except for the piglet that originated from the atypical herd mentioned earlier. of the otus the species scarcely represented were bifidobacterium boum, and corynebacterium kutscheri and they were mainly diminished from the colon. bifidobacterium as probiotic has been demonstrated to protect piglets against weaning diarrhoea associated with e. coli, and b. boum to inhibit shiga toxigenic e. coli virulence gene expression experimentally. this fits the fact that case piglets harbored an inverse proportion of members from the taxonomical lineage class beta- and gammaproteobacteria, family enterobacteriaceae and species e. coli, and from phylum actinobacteria, having a lot more of the former bacteria and fewer of the latter. c. kutscheri has primarily been reported to be a common bacterium on the mucous membranes of rodents including the gastrointestinal tract. interestingly, c. kutscheri was more plentiful in the intestinal content of healthy piglets than in that of diarrhoeic piglets, irrespective of gut section. the bacterium does not count as part of the normal porcine gut microbiota in the existing literature. thus, even though this cannot be ruled out for certain, it is possible that piglets may have acquired this species from foraging rodents directly or indirectly. class deltaproteobacteria this bacterial group was not associated with the diarrhoeic status of the piglets, but there was a significant reduction in numbers of members of this class in case piglets suffering from prolonged diarrhoea. moreover, deltaproteobacteria were also diminished by administration of antibiotics to the piglets. these facts are intertwined, as case piglets were mainly given antibiotics in an attempt to cure prolonged diarrhoea. finally, both results were not strongly significant. the only otu found by sequencing amplicons generated by the primer targeting this class was species desulfovibrio piger, and this was only identified in healthy piglets from two herds. this species has primarily been described as an opportunistic pathogen in humans but has also been detected in the feces of wild ducks. while the species was primarily found in colonic content of the control piglets, it is postulated that this species constitutes part of the normal gut microbiota of piglets. pcoa of the gut microbiota from the bion-meta processed sequencing data a pcoa was generated. this investigation supported the preliminary pca that was produced from the gut microbiotassay data, and additionally the major loading scores confirm the findings in the section above concerning the different bacterial groups. it diverges from the preceding discussion by taking all bacteria into account, and the clustering strongly suggests that it is possible to associate certain bacteria with a healthy gut microbiota. the intestinal microbiota of the aging piglet when the data were subset according to age in days, this study found that, with age, the gut microbiota clustered in relation to status, which probably reflects the significance of e. coli in nnpd. control piglets had significantly fewer of this species compared with case piglets, and over time the abundance of e. coli diminished in control piglets whilst it increased in case piglets. a more defined microbiota with increasing age is in line with conclusions from a previous study, which found that seven- to nine-day-old piglets had the lowest individual average similarity, independent of environmental or maternal relationship, compared with piglets from older age groups. the more distinctive clustering of nnpd status with increasing age could hint age to be a relevant factor in the diagnosis of nnpd. antibiotics the nnpd piglets investigated in this study included antibiotic-treated piglets. antibiotic-treated piglets were included because they typically suffered the worst cases of nnpd and were medicated for ethical reasons. it is evident that antibiotics have an impact on the gut microbial composition. however, one of the characteristics of nnpd is that it responds poorly to antibiotics. additionally, figures  and illustrate how it is possible to cluster piglets according to nnpd status despite the fact that some of the nnpd cases are treated with different antibiotics, indicating that the effect of nnpd on the gut microbial distribution and composition is more significant than the effect exerted by antibiotics in piglets of this age group. considerations on study approach case and control piglets included in this study all originate from herds affected with nnpd, hence the study lack the inclusion of true control piglets from non-affected herds. nonetheless, at the time of selecting the herds it was virtually impossible to find a true control herd due to the unknown aetiology and the much limited knowledge of nnpd. this also means that despite great efforts to ensure correct inclusion of nnpd herds through predefined criteria there is still a risk that piglets actually did not suffer from nnpd. it is acknowledged that this study only examines the bacterial composition of the luminal content, and not bacteria from the gastrointestinal wall. however, bacteria associated with the mucosa are not restricted to this site, as they are shed into the lumen together with intestinal epithelial cells and mucus. thus, it is assumed that this study also detects these or at least a sub-fraction of these. additionally, results from this study are affected by the various steps involved in preparing the samples, for example dna extraction, pcr, and sequencing, which all introduce variation into the final outcome. furthermore, the results do not reflect bacterial numbers but relative values of 16s and 23s rrna gene abundance or otus. conclusions drawn from sequencing data are based on otus displaying similar trends to data obtained from the gut microbiotassay. moreover, sequences had to be represented in a noteworthy number, and with somewhat consistent findings in piglets, irrespective of the variables in relation to diarrhoea. this means that there might be a risk that important bacteria could be overlooked if they do not meet the aforementioned criteria. finally, to determine causality based on the associations found in this study these should be investigated further for instance through bacterial isolation and inoculation studies.  the results of this study on nnpd indicate that bacteria could be the aetiology of this diarrhoea that affects piglets during the first week of life. at present, it is impossible to conclude whether diarrhoea is a consequence of an absence of beneficial bacteria or if diarrhoea is an outcome of invading pathogenic bacteria or an overgrowth of opportunistic pathogenic bacteria. e. coli appears to be a contributing factor to nnpd as nnpd piglets differ from healthy piglets by colonic colonization of bacteria from class beta- and gammaproteobacteria, primarily being species e. coli, and it appears that there is an important effect of age, both of which might be relevant in the characterization and diagnosis of nnpd. there is reason to believe that genus enterococcus is a participating factor, most likely species e. hirae or species e. durans. furthermore, diarrhoeic piglets appear to suffer from an imbalanced gut microbiota, in which bacteria regarded as beneficial are diminished, which particularly accounts for l. acidophilus. the associations found in this study are currently being investigated further. ethical according to danish laws no ethical approval is required for studies not including treatment groups or needle injections/blood testing. in addition no ethical approval is required for euthanization of animals by veterinarians. therefore this study was not subject to ethical approval, but fulfilled the regulations from the danish ministry of justice. hence, all handling of animals was performed by trained personnel and veterinarians, and euthanization was executed exclusively by veterinarians. procedures concerning the animals were all part of routine examinations and diagnosis of animals normally used in practice. endnotes gilts: and weeks ante partum, sows: weeks ante partum. etec: enterotoxigenic escherichia coli. prrs: porcine reproductive and respiratory syndrome. enzyme-linked immunosorbent assay. immunoperoxidase test. additional files additional file 1: primer modifications introduced to improve their performance and accompanying ribosomal database project search results. nucleotide explanation: m = c/a, y = t/c. as the legends states, this table summarizes the modifications introduced in two primer sets from the gut microbiotassay: “domain bacteria b v4-v5” and “phylum firmicutes”, respectively. ribosomal database project search results are given as percentage coverage of intended target group to demonstrate the improved performance of the modified primers tested in silico. additional file 2: loadings, eigenvalues, and explained variance by the major principal components related to figures 1, 3, and 4, respectively. additional file 3: bion-meta generated output. this file describes the bion-meta pipeline of the bioinformatics analysis of the sequencing data. in addition, it contains a link to a website holding the bion-meta output generated in this study. finally, it also includes a link to a website from where bion-meta can be downloaded. additional file 4: results from linear discriminant analysis of data related to figure  1 and figure  3. linear discriminant analysis was used to classify observations into cases or controls by their fluidigm-derived microbial profile and to estimate how well this classification corresponded with the actual status of the sample. this classification was followed up with a multivariate analysis of variance using wilks lambda statistic. abbreviations aa access array 454blaccess array barcode library for the gs flx titanium sequencer cqquantification cycle decdiarrhoeagenic e. coli ehec 2enterohemorrhagic e. coli clonal complex etecenterotoxigenic e. coli fishfluorescence in situ hybridization lnnatural logarithm ldalinear discriminant analysis lsularge subunit nnpdnew neonatal porcine diarrhoea npdneonatal porcine diarrhoea otuoperational taxonomic units pcaprincipal component analysis pcoaprincipal coordinates analysis prrsporcine reproductive and respiratory syndrome qpcrquantitative real-time polymerase chain reaction rrnaribosomal rna scfashort chain fatty acids ssusmall subunit competing interests the authors declare that they have no competing interests. authors’ contributions mlhb, lm, mb and ks participated in conceiving and designing this study. the inclusion of pig farms and clinical examination of animals was performed by hk. mlhb completed all the laboratory experiments. nl advanced bion-meta to also include search results from ribosomal lsu databases, and hci assisted in the improvement of bion-meta for this study, and critically revised the manuscript. as conducted all statistical calculations related to the qpcr results and mls contributed with the multivariate analysis of the qpcr data and the sequencing data. mlhb drafted the manuscript, and lm, mb and ks helped along by revising and proofreading the manuscript for appropriate content. all authors read and approved the final manuscript. 
nnpd was associated with a diminished quantity of bacteria from the phyla actinobacteria and firmicutes while genus enterococcus was more than times more abundant in diarrhoeic piglets. nnpd affects piglets during the first week of life and results in impaired welfare, decreased weight gain, and in the worst-case scenario death. the results of this study support previous findings that points towards genus enterococcus and species e. coli to be involved in the pathogenesis of nnpd. moreover, the results indicate that nnpd is associated with a disturbed bacterial composition and larger variation between the diarrhoeic piglets. this study examined ileal and colonic luminal contents of control piglets and nnpd piglets by means of the qpcr-based gut microbiotassay and samples by sequencing to study the composition of the bacterial gut microbiota in relation to nnpd. with increasing age, the gut microbiota of nnpd affected piglet and control piglets became more diverse. nnpdneonatalpigletdiarrhoeaqpcrmicrobiotagut microbiotassay sequencingissue-copyright-statement© the author 2015 the number of bacteria from the phylum fusobacteria was also doubled in piglets suffering from diarrhoea. independent of diarrhoeic status, piglets from first parity sows possessed significantly more bacteria from family enterobacteriaceae and species e. coli, and fewer bacteria from phylum firmicutes. previous studies have investigated the clinical manifestations, histopathology, and to some extent, microbiological findings; however, these studies were either inconclusive or suggested that enterococci, possibly in interaction with escherichia coli, contribute to the aetiology of nnpd. piglets born to gilts had times higher odds of having nnpd compared with piglets born to multiparous sows. finally, the co-occurrence of genus enterococcus and species e. coli contributed to the risk of having nnpd. in recent years, new neonatal porcine diarrhoea of unknown aetiology has emerged in denmark.
 in recent years, new neonatal porcine diarrhoea of unknown aetiology has emerged in denmark. nnpd affects piglets during the first week of life and results in impaired welfare, decreased weight gain, and in the worst-case scenario death. commonly used preventative interventions such as vaccination or treatment with antibiotics, have a limited effect on nnpd. previous studies have investigated the clinical manifestations, histopathology, and to some extent, microbiological findings; however, these studies were either inconclusive or suggested that enterococci, possibly in interaction with escherichia coli, contribute to the aetiology of nnpd. this study examined ileal and colonic luminal contents of control piglets and nnpd piglets by means of the qpcr-based gut microbiotassay and samples by sequencing to study the composition of the bacterial gut microbiota in relation to nnpd.  nnpd was associated with a diminished quantity of bacteria from the phyla actinobacteria and firmicutes while genus enterococcus was more than times more abundant in diarrhoeic piglets. the number of bacteria from the phylum fusobacteria was also doubled in piglets suffering from diarrhoea. with increasing age, the gut microbiota of nnpd affected piglet and control piglets became more diverse. independent of diarrhoeic status, piglets from first parity sows possessed significantly more bacteria from family enterobacteriaceae and species e. coli, and fewer bacteria from phylum firmicutes. piglets born to gilts had times higher odds of having nnpd compared with piglets born to multiparous sows. finally, the co-occurrence of genus enterococcus and species e. coli contributed to the risk of having nnpd.  the results of this study support previous findings that points towards genus enterococcus and species e. coli to be involved in the pathogenesis of nnpd. moreover, the results indicate that nnpd is associated with a disturbed bacterial composition and larger variation between the diarrhoeic piglets. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords nnpdneonatalpigletdiarrhoeaqpcrmicrobiotagut microbiotassay sequencingissue-copyright-statement© the author 2015 