Piglet nasal microbiota at weaning may influence the development of Glässer’s disease during the rearing period
© Correa-Fiz et al. 2016
Received: 22 December 2015
Accepted: 6 May 2016
Published: 26 May 2016
The microbiota, the ensemble of microorganisms on a particular body site, has been extensively studied during the last few years, and demonstrated to influence the development of many diseases. However, these studies focused mainly on the human digestive system, while the populations in the respiratory tract have been poorly assessed, especially in pigs. The nasal mucosa of piglets is colonized by an array of bacteria, many of which are unknown. Among the early colonizers, Haemophilus parasuis also has clinical importance, since it is also the etiological agent of Glässer’s disease. This disease produces economical losses in all the countries with pig production, and the factors influencing its development are not totally understood. Hence, the purpose of this work was to characterize the nasal microbiota composition of piglets, and its possible role in Glässer’s disease development.
Seven farms from Spain (4 with Glässer’s disease and 3 control farms without any respiratory disease) and three farms from UK (all control farms) were studied. Ten piglets from each farm were sampled at 3–4 weeks of age before weaning. The total DNA extracted from nasal swabs was used to amplify the 16S RNA gene for sequencing in Illumina MiSeq. Sequencing data was quality filtered and analyzed using QIIME software. The diversity of the nasal microbiota was low in comparison with other body sites, showing a maximum number of operational taxonomic units (OTUs) per pig of 1,603, clustered in five phyla. Significant differences were found at various taxonomical levels, when the microbiota was compared regarding the farm health status. Healthy status was associated to higher species richness and diversity, and UK farms demonstrated the highest diversity.
The composition of the nasal microbiota of healthy piglets was uncovered and different phylotypes were shown to be significantly altered in animals depending on the clinical status of the farm of origin. Several OTUs at genus level were identified over-represented in piglets from control farms, indicating their potential as probiotics. Although we provide relevant data, fully metagenomic approaches could give light on the genes and metabolic pathways involved in the roles of the nasal microbiota to prevent respiratory diseases.
Many body sites from mammalian species provide habitats for a huge number of microbial species. Studies on microbiota, the microbial composition of a body site, have increased over the last years with the arrival of NGS technology. The influence of the microbiota composition on the development of many diseases [1–3], together with the possibility of manipulating that composition to control disease [4–7]), have conferred a central role to these studies during the last few years. Most of the microbiota studies have been done in humans, although other species are also being investigated [8–11]. Additionally, the majority of the studies have focused on the bacterial gut populations, through the analysis of intestinal tract or feces, although body sites such as skin, mouth, nose or vagina have also been investigated [12–15]. Other sites related to respiratory diseases have been also analyzed, such as lungs or tonsils; especially in animal models due to invasive nature of the sampling process [16–18].
Among the animal models, pigs are important for society, mainly because pork is the world’s most consumed meat from terrestrial animals (http://www.fao.org). Pig industry has risen globally and the emergence of specific zoonotic diseases is of increasing concern. Pigs can act as reservoirs of several pathogens becoming a risk to farmers and consumers. Moreover, because pigs have demonstrated to be one of the best models for human diseases [19, 20], the unraveling of the microbiota composition is of great interest.
Most microbiota studies have focused on species related to digestive system while the populations related to respiratory have been poorly assessed for pigs. The tonsils have historically been used to explore the upper respiratory microbial populations, but the fact that nose is the first place where the microbes find a barrier to infection, studies on nasal microbiota can be relevant. The nose can harbor important microorganisms that can be pathogenic under certain circumstances.
Among the microbes widely distributed in swine population, Haemophilus parasuis is a small, pleomorphic, Gram-negative rod that colonizes the upper respiratory tract of healthy piglets . However, under some conditions (stress, lack of immunity, presence of other pathogens,...), virulent strains can spread to the lung and invade systemic sites, causing Glässer’s disease [22, 23], which is characterized by fibrinous polyserositis, including meningitis. Glässer’s disease is considered one of the bacterial diseases with highest economic impact in nursery piglets (4–8 weeks of age) and it is present in all the major swine-raising countries. This disease is normally controlled by antimicrobial treatments; however antibiotic therapies have been shown to select resistant H. parasuis strains in farms . Moreover, antimicrobial treatments can deplete members of the resident commensal bacterial population, which may be beneficial in preventing colonization with virulent H. parasuis. Actually, many members of the microbiota are known to be selected by the host to prevent colonization by pathogens, and that could be the case for H. parasuis. The knowledge of the composition of the microbiota of animals from farms without respiratory diseases and the potential differences with animals from farms with Glässer’s disease outbreaks will be of enormous relevance to better understand this disease and its control. Moreover, the different microorganisms from healthy animals can be used to further study their potential benefits.
Here, we report the nasal microbiota composition of pigs at weaning and its potential relation to susceptibility to respiratory conditions, such as Glässer’s disease.
The nasal microbiota of piglets from control and Glässer’s diseased-farms
Samples subjected to individual massive sequencing were taken from ten pigs per farm at 3–4 weeks of age, just before weaning. We analyzed samples from animals from four farms with Glässer’s disease outbreaks (GD farms) and from six control farms without any respiratory condition reported (C farms). H. parasuis colonize piglets at the early beginning of life, but Glässer’s disease normally develops at 4–8 weeks of age, after weaning. Thus, the analysis was performed with samples taken at 3–4 weeks of age to study the possible predisposition to disease.
We obtained a total of 37,376,035 sequences from four different runs. The reads were joined and the erroneous or poor quality reads (Q < 25) were eliminated. This left us with 45 ± 18 % of the reads, corresponding to a median number of reads per sample of 232,807 (range 7,338–844,521) with a mean sequence length of 456.77 ± 18.6 nt. Singletons and samples with low coverage (less than 40,000 reads per sample) were discarded, so as to have a more representative and homogenous data.
The most prevalent genus found in the nasal microbiota was Moraxella (17.2 %), followed by Weeksella (12.9 %). Haemophilus was also abundant, representing 6.1 % of the reads classified at genus level, and 27.0 % of the Pasteurallaceae family. Almost half of the genera in the Lactobacillaceae family were Lactobacillus, with 49.0 % of relative abundance within this family. Unfortunately, not all sequences could be confidently assigned to one genus (21.7 % assigned as ‘other’) and may deserve further studies to be characterized.
To unravel the differences in nasal microbiota in piglets from farms with a good healthy status (C farms) in contrast to those from farms with Glässer’s disease (GD farms), we analyzed the relative abundance of OTUs at three main levels (phylum, family and genus) by grouping samples according to the health status of the farm of origin (Fig. 1). Figure 1 shows the averaged relative abundance per health status in two groups: GD and C farms. At phylum level, there is a relative increase in Proteobacteria (32.5 % in C to 44.0 % in GD) in conjunction to a decrease in Firmicutes (21.1 to 18.5 %) in animals from farms with Glässer’s disease. The increase in Proteobacteria corresponds to a higher abundance of Pasteurellaceae (5.1 % vs 9.5 % in C and GD, respectively) and Moraxellacea (19.3 % vs 27.8 % in C and GD, respectively) families. Within the Pasteurellaceae family, the genus Haemophilus was found increased in pigs from farms with Glässer’s disease (2.8 % vs 9.4 % in C and GD, respectively). Here, Haemophilus corresponds to the unique porcine species of this genus, Haemophilus parasuis, the etiological agent of Glässer’s disease. The Moraxellaceae family showed a raised amount of the genera Moraxella (13.6 to 22.2 % in C and GD, respectively) and Enhydrobacter (3.1 to 4.3 % in C and GD, respectively). Other bacterial components of the microbiota, although present in low abundance, deserve further investigation due to their absence in animals from diseased farms. This is the case for Psychrobacter, from the Moraxellaceae family, which was only present in piglets from C farms (at 1.8 %). In addition, within Firmicutes, the principal changes were seen in the Ruminococcaceae family, which decreased from 8.0 % in C farms to 3.4 % in GD ones. These reductions are the reflection, at least in part, of the diminished amount of Oscillospira (4.5 to 1.8 % in C and GD, respectively) and Ruminococcus, with 1.0 % in animals from C farms but absent in GD farms. Lachnospiraceae family showed a significant decrease in GD farms (2.9 to 1.1 %), which was due to many low-abundant genera that diminished simultaneously. Many of the changes observed in Firmicutes could not be associated to any family, because the sequence assignment was done confidently only until order level, and corresponded to Clostridiales. In contrast to this general decreased behavior in Firmicutes, the Streptococcaceae family showed a significant raise, revealed by the increased tendency observed in Streptococcus (from 2.1 to 2.9 % in C and GD, respectively).
Effect of health status, country and farming company on beta diversity
To understand the differences across the microbiota on nasal samples, the beta diversity was estimated. Weighted and unweighted UniFrac phylogenetic distances were used to generate the beta diversity distance matrices and calculate the degree of differentiation among the samples. Samples were grouped according to different characteristics to test them as possible factors leading to clustering.
When we focus on the 3 farms with the same sow origin, the qualitative differences observed between farm health status was significant, with R 2 = 0.1821 (Fig. 3c), explaining 18 % of the variation between groups, the highest observed in all cases. Based on these results, the largest source of the variation in bacterial community composition between samples may be attributed to differences in sow origin, which could be associated to genetic background.
The weighted analysis did not show differences when comparing country of origin or health status, indicating that the differences observed between groups (those demonstrated through the unweighted analysis) may correspond to species in lower abundance, as the weighted UniFrac emphasizes common OTUs and deemphasizes rare OTUs [36, 37]. However, when analyzing only samples from the same sow origin separately, both weighted (Fig. 3d) and unweighted (Fig. 3c) UniFrac distances showed good clustering according to health condition (weighted Unifrac analysis, R 2 = 0.1917 p < 0.05; unweighted UniFrac analysis, R 2 = 0.2222 p < 0.05). Thus, sow origin may be a relevant factor to determine the nasal microbiota composition in pigs, whose modifications would affect the health status of the piglets.
The core microbiota in Spanish and British farms
Differences at several taxonomy levels were found in the nasal microbiota composition of piglets depending on the country of origin. However, some of these were due to low abundant OTUs, which may be a result of the specific moment of sampling rather than colonization of the particular taxon. The analysis of the core microbiome, defined as those OTUs present in all the samples from C farms, was made with the aim of finding permanent/essential inhabitants of the nasal cavity of pigs (colonizers).
Relative abundance of significantly different OTUs of the healthy core from Spanish (Sp) and British (UK) farms
Relative abundance (%)
The increased number of members within Firmicutes in the UK core was demonstrated to be related mainly to the order Clostridiales (Table 1). Although many of these differences can be associated to several genera, more studies should be made to make a step ahead to understand the species involved.
The apparent lack of changes in Bacteroidetes between the UK and Sp cores are the result of several compensatory modifications at family or genus levels. Bacteroidales order was increased in the UK core, with an increment in the Bacteroidaceae family and Bacteroides genus. Also, Rikenallaceae and S24-7 families showed increased levels in UK core in contrast to what is observed in Sp core. Concomitantly, an opposite tendency was observed for the Flavobacteriales order, which was increased in the Sp core, due mainly to an increase in the Weeksellaceae family, mostly related to an increase in Weeksella.
The genus Myroides (from the Flavobacteriaceae family) was especially abundant in most UK farms while absent in Spain, probably being specifically related to location (data not shown). However, this result is not found in the core analysis, because Myroides was absent in some animals from a particular farm. These particularly abundant taxa deserve further investigation on the potentially protective role they may have against diseases.
Comparison of relative abundance of core nasal microbiota composition of animals from farms with different health status
With the aim of finding microbes that could be associated to health on piglets, the core nasal microbiota composition of C pigs was compared with GD pigs, only in farms from Spain. The presence of the particular OTUs found was also contrasted in UK farms. Moreover, we decided to analyze the nasal microbiota composition of wild boars (Sus scrofa) from Catalonia (Spain). There is little information available regarding H. parasuis infection in these animals with only one report showing evidence on Glässer’s disease . So, the presence of a particular OTU in nasal microbiota of healthy pigs (Sus scrofa domesticus), absent in those pigs from farms with Glässer’s outbreaks, would become more relevant if found also in UK farms and/or in wild boars, as they seem to be more resistant to the disease .
Differential OTUs in diseased (GD) and control (C) Spanish farms; comparison with wild boars and UK farms
Relative abundance (%)
Wild Boars (Spain)
Differential OTUs in healthy or diseased animals from a Spanish farm (EJ) with Glässer’s disease (GD)
Reative abundance (%)
Elucidation of the normal bacterial communities on body sites is critical for establishing differences associated with disease. In this study, the nasal microbiota composition of piglets was explored through direct sequencing of the 16S rRNA gene, which provided an overview of the nasal microbiota of healthy piglets, and the changes in this dynamic population associated to Glässer’s disease.
A few studies exist on the microbiota of the respiratory tract of pigs, although they mainly focused on tonsils or lungs [18, 40, 41], with only a recent study analyzing the nasal microbiota of healthy pigs . Differential microbiota composition has been associated to several intestinal diseases , but this is the first report to study its relationship with Glässer’s disease, one of the main respiratory diseases of swine . The nasal microbiome of healthy piglets contained a low diverse bacterial community in comparison with other body sites, such as gut . In general, alpha diversity was found to be higher in control farms with no respiratory disease, C farms, when compared to GD farms. High microbial diversity is thought to be beneficial to the mucosal surfaces by reducing the opportunity for colonization by pathogens . To our knowledge, this is the first report to relate microbiota diversity in the pig nares with Glässer’s disease, and is in agreement with many other studies where lower diversity correlated with disease [46–49].
Regarding the phylotypes found in pigs from C and GD farms, remarkable differences were found. A notably gain in Proteobacteria and loss in Firmicutes when the animals came from farms with Glässer’s disease was seen at phylum level. Interestingly, the Pasteurallaceae family was, at least in part, responsible for the higher amount of Proteobacteria, with Haemophilus, and particularly H. parasuis, being clearly increased. The fact that H. parasuis comprises strains with different degree of virulence makes necessary a deeper characterization of the virulence of the specific strains in the pigs to assess their role on Glässer’s disease outbreaks on each farm. As an example, the healthy status of farm VL could be explained by the lack of virulent H. parasuis strains in these animals (not detected by specific PCR , data not shown), since the global diversity of species in this farm was low and close to GD farms. Importantly, characterization of H. parasuis strains, in conjunction to the fact that they are statistically more abundant on diseased farms, could become useful as a non-invasive prognostic tool.
The prevalence of pathogens in the nares has been linked to increased risk of infection by other pathogens . Relevant findings of the present study include the fact that the Mycoplasmataceae family was increased in GD farms, together with an increase in Streptococcus and Haemophilus. We hypothesized that the increased abundance of these pathogens could confer a higher risk to develop Glässer’s disease, but this should be experimentally demonstrated. On the other hand, Moraxella was shown to be a relevant component of the nasal microbiota, but it was found at higher abundance in diseased farms, declining the possibility of this genus to have a protective role in the nasal cavities, as it was hypothesized before . However, we cannot rule out a role in protection with our results, since a deeper characterization of the Moraxella species and the precise strains present in the nasal cavities would be needed. In addition, the significant gain in Proteobacteria cannot be exclusively explained by these two genera mentioned (Haemophilus and Moraxella), and needs deeper research in order to understand the potential detrimental consequences, as many pathogens belong to this phylum.
More importantly, other relevant decreases were observed in some bacterial taxa from the nasal microbiota from GD farms, which might be predictive of beneficial commensals. This could be the case of both the Ruminococcaceae and the Lachnospiraceae families, which are also commensal members of the pig gut microbiota and have been associated to health in previous reports . The same rational would apply to the order Clostridiales, which was reduced in GD farms and deserves also further attention. Moreover, the interaction between members of the gut microbiota and the respiratory system of the host have been reported [53, 54], conferring these data a more relevant interest as the gut microbiota can be manipulated with specific diet . Additionally, the importance of pigs does not lay only in its relevance in food industry, but also in the physiological similarities to human, in terms of nutritional requirements and regarding the intestinal microbial ecosystem .
A phylogeny-based metric, UniFrac , was used to compare alpha diversity and assess the differences in overall bacterial community composition (beta diversity). This UniFrac-based principal coordinate analysis revealed clustering when the health status of the farm was evaluated. Clustering was observed when samples from the same country of origin were compared (Fig. 3a), although only around 10 % of the variability observed was explained, leaving a high amount of unexplained variation. This is typical for ecological studies  where it is not possible to control all environmental variables (e.g. external environmental factors, biological interactions, etc.). The clustering of samples from the same sow origin showed the highest percentage of variation. It is worth to mention that this group is composed of farms from a single company and with similar locations. We hypothesize that since the sows have the same origin, the predisposition to Glässer’s disease could not be associated to the genetic background, as it has been stated before for other diseases [31, 32], and the differences in microbiota would be due to other management factors. However, this hypothesis need further research, mainly due to the low number of the cohort study when they were subgrouped.
When clinical symptoms at about 8 weeks of age were considered for a GD farm (EJ farm), the beta diversity analysis evidenced individual differences in animals sharing the same environment, location, sow origin and other uncontrolled factors. Thus, the microbiota inhabitants of the nares of piglets may have an effect on the predisposition to develop Glässer’s disease, demonstrated by the fact that animals that remained healthy had similar microbiota composition to animals from C farms (Fig. 3a, red spheres with white dots).
Most of the differences observed in this study were detected by unweighted UniFrac metrics. It is reported that the abundance information in weighted UniFrac matrices do not consider the variation of the weights under random sampling, thus resulting in less power for detecting differences in communities [36, 37]. So the fact that almost no differences were detected in weighted UniFrac analysis can be related to those limitations, although it has to be considered that the qualitative information obtained through unweighted analysis can be emphasized by shallow OTUs.
Pigs are continuously rooting for food in farms, and this behavior could favor the transitory presence of microorganisms that will not finally colonize the nose. The analysis of the core microbiota was made with the aim of finding permanent inhabitants of the nasal cavities of pig (colonizers). The microbiota composition showed differences depending on the localization of the farms, as it was expected (Fig. 3b). However, those differences cannot be exclusively associated to country origin, as management practices, such as antibiotic treatments, were distinct too. Those OTUs in higher abundance in UK farms should be studied carefully, as these farms demonstrated higher diversity (Additional file 2), probably linked to their healthy status. In particular, farm KD was reported to eliminate respiratory problems after eliminating early antibiotic treatments , which may have an impact in increasing bacterial diversity.
The nasal microbiota composition of wild boars, which seem to be more resistant to Glässer’s disease, was explored for the first time and contrasted to the Spanish and British cores (Table 2). Interestingly, most of the OTUs found associated to health in Spanish farms were also found in similar abundance both in UK farms and wild boars, which are not exposed to massive use of antibiotics that can disturb the normal microbiota. These healthy animals showed to have more Bacteroides, pointing them as the main responsible of the differences observed. The exception of this was observed for Haemophilus, which was found in a high relative amount in wild boars, although being healthy animals. However, as stated before, the virulence of the specific strains needs to be investigated. This observation could imply the importance of the genetic background on predisposition to both acquisition of differential nasal microbiota and developing Glässer’s disease.
The information that provides this work, unraveling the differences among communities in the nasal mucosa, must be taken into consideration for future research related to Glässer’s disease control.
Comparison of the nasal microbiota of piglets from farms experiencing Glässer’s disease and control farms revealed 8 OTUs (genus level) that were over-represented in the healthy animals, indicating a potential beneficial effect by these microorganisms. Besides, in a farm with disease, we detected differences in the nasal microbiota in animals that finally developed disease in comparison to those that remained healthy.
In conclusion, we were able to show through a confident non-invasive approach that the nasal cavity is colonized by a relatively diverse microbiota in healthy piglets. The nasal inhabitants showed to be in relation to clinical status of the farm of origin of the piglets, leading to different susceptibilities to invasive infection by H. parasuis. Several members have been identified as possible candidates to have a protective role against respiratory diseases and should be further studied. Although a 16S rRNA gene-based technique was used to reveal the microbial diversity in the nasal cavity, this approach offers only limited information on the physiological role of the microbial consortia. Fully metagenomic studies could give impressive information on the genes, metabolic pathways and resistances, among others characteristics related to the major roles that the microbiota plays in health and disease, giving more light on the every-day-less “forgotten organ” .
Characteristics of the farms included in the study
Sow farm size
Glässer’s disease (GD)
Pen + Strep
Glässer’s disease (GD)
Farrow to finish
Glässer’s disease (GD)
Glässer’s disease (GD)
Farrow to finish
Tul + Ceft
Farrow to finish
Farrow to finish
Nasal swabs were taken from the nares of 10 animals and placed into a sterile tube. Swabs were transported to the laboratory on ice where they were resuspended in 500 μl of PBS and stored at −20 °C for further analysis.
Genomic DNA was extracted using the Nucleospin Blood (Machinery Nagel) kit. Briefly, frozen swabs were vortex vigorously and 200 μl were used for DNA extraction. Purified DNA was resuspended in a final volume of 50 μl of elution buffer (5 mM Tris, pH 8.5). The quality and quantity of genomic DNA was evaluated on a BioDrop DUO (BioDrop Ltd). All the extractions were done the same day the preparation of the library was done for sequencing purposes.
16S rRNA gene sequencing
Extracted DNA was sent to IBB facilities for Illumina pair-end 2X250 bp sequencing with MiSeq following the manufacturer instructions (MS-102-2003 MiSeq® Reagent Kit v2, 500 cycle).
16S Forward Primer 5′ TCGTCGGCAGCGTCAGATGTGTATAAGAGACAGCCTACGGGNGGCWGCAG
16S Reverse Primer 5′ GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAGGACTACHVGGGTATCTAATCC
Each of the primers contained a unique six-nucleotide barcode, so that the derived sequences can be sorted into the respective sample bioinformatically in downstream analysis. Cycling conditions included initial denaturation at 95 °C for 3 min followed by 25 cycles of 95 °C for 30 s, 55 °C for 30 s, and 72 °C for 30 s and a final extension at 72 °C for 5 min. The PCR product were purified and checked to verify its size on a Bioanalyzer DNA 1000 chip (Agilent).
Each sequence was assigned to its original sample according to its barcode. The barcode was removed from sequences before further processing, so that sequences that did not contain the barcode were eliminated.
Taxonomic microbiota analysis
The QIIME (quantitative insights into microbial ecology ) software package (version 1.8) was used to process the reads following the protocols described by Kuczynski et al.  with some modifications. Firstly, reads were filtered in order to keep only high quality reads (Q > 25) as recommended in Bokulich et al.  for improving diversity estimates (using split_libraries_fastq.py script from QIIME). The high-quality-reads were paired-end joined by using fastq join [61, 62] using default values. OTUs were assigned at 97 % similarity using UCLUST algorithm  and the Greengenes database . Sequences with dissimilarity higher than 3 % or lengths shorter than 350 bp were discarded from the analysis. A representative sequence from each of the OTUs, corresponding to the cluster seed, was subjected to taxonomic assignment with RDP database .
Samples were grouped according to health status (GD or C) initially, and further analyzed based on location and sow origin. Analysis was performed at different taxonomical levels (phylum, family and genus) separately. For each taxon g-test was performed between different groups and p values were FDR-corrected for multiple hypotheses testing (with group_significance.py QIIME script). Diversity indices (including Shannon-Wiener index, Faith’s PD Whole Tree and rarefaction curves) were calculated on rarefied 16S rRNA gene sequence data for all samples at 97 % similarity using QIIME (alpha_diversity.py script). Alpha diversity between groups was compared through non-parametric tests (Monte Carlo method) at maximum depth in rarefied samples (with 999 permutations) using compare_alpha_diversity.py QIIME script.
In addition, equal number of samples was subsampled to assess the significant differences between sample types using UniFrac weighted and unweigthed distances (jackknifed_beta_divesity.py QIIME workflow script) [37, 66]. The percentage of variation between grouped samples was measured by R 2, using Adonis function of the vegan package in R software . Estimation of p-values was obtained through Monte Carlo test with 999 random permutations of the data set (compare_categories.py, QIIME script). Samples were considered to be significantly different when the accompanying P-value was < 0.05.
C: control; GD: Glässer’s disease; OTU: operational taxonomic unit.
This work was supported by grant AGL2013-45662 from the Ministerio de Economía y Competitividad of Spain. The authors would like to acknowledge Mark White for providing the samples from the UK, Oscar Cabezón (UAB) for the wild boar samples and Núria Galofré-Milà (IRTA-CReSA) for technical help.
Availability of supporting data
The entire sequence dataset is available in the NCBI database. SRA Accession number SRP068182, BioProject 305822, BioSample 04334499.
VA conceived the study and designed the experiments, LF selected farms and samples, FCF analysed the data and wrote the manuscript, with contributions from VA and LF. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Ghaisas S, Maher J, Kanthasamy A. Gut microbiome in health and disease: linking the microbiome-gut-brain axis and environmental factors in the pathogenesis of systemic and neurodegenerative diseases. Pharmacol Ther. 2015. doi:10.1016/j.pharmthera.2015.11.012.PubMedGoogle Scholar
- Mejia-Leon ME, Barca AM. Diet, microbiota and immune system in type 1 diabetes development and evolution. Nutrients. 2015;7:9171–84.View ArticlePubMedPubMed CentralGoogle Scholar
- Wu H, Tremaroli V, Backhed F. Linking microbiota to human diseases: a systems biology perspective. Trends Endocrinol Metab. 2015;26:758–70.View ArticlePubMedGoogle Scholar
- Li YT, Cai HF, Wang ZH, Xu J, Fang JY. Systematic review with meta-analysis: long-term outcomes of faecal microbiota transplantation for Clostridium difficile infection. Aliment Pharmacol Ther. 2015. doi:10.1111/apt.13492.Google Scholar
- Kelly CR, Kahn S, Kashyap P, Laine L, Rubin D, Atreja A, Moore T, Wu G. Update on fecal microbiota transplantation 2015: indications, methodologies, mechanisms, and outlook. Gastroenterology. 2015;149:223–37.View ArticlePubMedGoogle Scholar
- Xu Z, Knight R. Dietary effects on human gut microbiome diversity. Br J Nutr. 2015;113(Suppl):S1–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Lindberg JE. Fiber effects in nutrition and gut health in pigs. J Anim Sci Biotechnol. 2014;5:15–21.View ArticlePubMedPubMed CentralGoogle Scholar
- Malmuthuge N, Griebel PJ, le Guan L. The Gut Microbiome and its potential role in the development and function of newborn calf gastrointestinal tract. Front Vet Sci. 2015;2:36–45.View ArticlePubMedPubMed CentralGoogle Scholar
- Roto SM, Rubinelli PM, Ricke SC. An introduction to the avian gut microbiota and the effects of yeast-based prebiotic-type compounds as potential feed additives. Front Vet Sci. 2015;2:28–45.View ArticlePubMedPubMed CentralGoogle Scholar
- Bendtsen KM, Fisker L, Hansen AK, Hansen CH, Nielsen DS. The influence of the young microbiome on inflammatory diseases-Lessons from animal studies. Birth Defects Res C Embryo Today. 2015. doi:10.1002/bdrc.2116.PubMedGoogle Scholar
- Panagi M, Georgila K, Eliopoulos AG, Apidianakis Y. Constructing personalized longitudinal holo’omes of colon cancer-prone humans and their modeling in flies and mice. Oncotarget. 2015. doi:10.18632/oncotarget.6463.PubMedGoogle Scholar
- Avila M, Ojcius DM, Yilmaz O. The oral microbiota: living with a permanent guest. DNA Cell Biol. 2009;28:405–11.View ArticlePubMedPubMed CentralGoogle Scholar
- Frank DN, Feazel LM, Bessesen MT, Price CS, Janoff EN, Pace NR. The human nasal microbiota and Staphylococcus aureus carriage. Plos One. 2010;5:1–15.Google Scholar
- Gajer P, Brotman RM, Bai G, Sakamoto J, Schutte UM, Zhong X, et al. Temporal dynamics of the human vaginal microbiota. Sci Transl Med. 2012. doi:10.1126/scitranslmed.3003605.PubMedPubMed CentralGoogle Scholar
- Grice EA, Kong HH, Renaud G, Young AC, Bouffard GG, Blakesley RW, et al. A diversity profile of the human skin microbiota. Genome Res. 2008;18:1043–50.View ArticlePubMedPubMed CentralGoogle Scholar
- Dickson RP, Huang YJ, Martinez FJ, Huffnagle GB. The lung microbiome and viral-induced exacerbations of chronic obstructive pulmonary disease: new observations, novel approaches. Am J Respir Crit Care Med. 2013;188:1185–6.View ArticlePubMedGoogle Scholar
- Poroyko V, Meng F, Meliton A, Afonyushkin T, Ulanov A, Semenyuk E, et al. Alterations of lung microbiota in a mouse model of LPS-induced lung injury. Am J Physiol Lung Cell Mol Physiol. 2015;309:L76–83.View ArticlePubMedGoogle Scholar
- Lowe BA, Marsh TL, Isaacs-Cosgrove N, Kirkwood RN, Kiupel M, Mulks MH. Defining the core microbiome of the microbial communities in the tonsils of healthy pigs. BMC Microbiol. 2012;12:1–14.View ArticleGoogle Scholar
- Bassols A, Costa C, Eckersall PD, Osada J, Sabria J, Tibau J. The pig as an animal model for human pathologies: a proteomics perspective. Proteomics Clin Appl. 2014;8:715–31.View ArticlePubMedGoogle Scholar
- Heinritz SN, Mosenthin R, Weiss E. Use of pigs as a potential model for research into dietary modulation of the human gut microbiota. Nutr Res Rev. 2013;26:191–209.View ArticlePubMedGoogle Scholar
- Olvera A, Cerda-Cuellar M, Nofrarías M, Revilla E, Segalés J, Aragon V. Dynamics of Haemophilus parasuis genotypes in a farm recovered from an outbreak of Glässer’s disease. Vet Microbiol. 2007;123:230–7.View ArticlePubMedGoogle Scholar
- Aragon V. Exposing serum susceptibility in Haemophilus parasuis. Vet J. 2013;196:10–1.View ArticlePubMedGoogle Scholar
- Costa-Hurtado M, Aragon V. Advances in the quest for virulence factors of Haemophilus parasuis. Vet J. 2013;198:571–6.View ArticlePubMedGoogle Scholar
- Biswas K, Hoggard M, Jain R, Taylor MW, Douglas RG. The nasal microbiota in health and disease: variation within and between subjects. Front Microbiol. 2015;9:134.View ArticlePubMedGoogle Scholar
- Lemon KP, Klepac-Ceraj V, Schiffer HK, Brodie EL, Lynch SV, Kolter R. Comparative analyses of the bacterial microbiota of the human nostril and oropharynx. mBio. 2010;1:e00129–10.PubMedPubMed CentralGoogle Scholar
- Slifierz MJ, Friendship RM, Weese JS. Longitudinal study of the early-life fecal and nasal microbiotas of the domestic pig. BMC Microbiol. 2015;15:184–95.View ArticlePubMedPubMed CentralGoogle Scholar
- Bassis CM, Tang AL, Young VB, Pynnonen MA. The nasal cavity microbiota of healthy adults. Microbiome. 2014;2:27–31.View ArticlePubMedPubMed CentralGoogle Scholar
- Henao-Mejia J, Elinav E, Jin C, Hao L, Mehal WZ, Strowig T, et al. Inflammasome-mediated dysbiosis regulates progression of NAFLD and obesity. Nature. 2012;482:179–85.PubMedPubMed CentralGoogle Scholar
- Peng J, Narasimhan S, Marchesi JR, Benson A, Wong FS, Wen L. Long term effect of gut microbiota transfer on diabetes development. J Autoimmun. 2014;53:85–94.View ArticlePubMedPubMed CentralGoogle Scholar
- Srinivas G, Moller S, Wang J, Kunzel S, Zillikens D, Baines JF, Ibrahim SM. Genome-wide mapping of gene-microbiota interactions in susceptibility to autoimmune skin blistering. Nat Commun. 2013;4:2462.View ArticlePubMedPubMed CentralGoogle Scholar
- Esworthy RS, Smith DD, Chu FF. A strong impact of genetic background on Gut Microflora in mice. Int J Inflam. 2010;2010:986046.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhao L. The gut microbiota and obesity: from correlation to causality. Nat Rev Microbiol. 2013;11:639–47.View ArticlePubMedGoogle Scholar
- Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Meth. 2010;10:335–6.View ArticleGoogle Scholar
- Kuczynski J, Stombaugh J, Walters WA, González A, Caporaso JG, Knight R. Using QIIME to analyze 16S rRNA gene sequences from microbial communities. Curr Protoc Bioinformatics. 2011. doi:10.1002/0471250953.bi1007s36.PubMedPubMed CentralGoogle Scholar
- Okasnen J, Blanchet F, Kindt R, Legendre P, Minchin PR, O’Hara B, et al. Vegan: Community Ecology Package. R package version 2.3-0. 2015.Google Scholar
- Fukuyama J, McMurdie PJ, Dethlefsen L, Relman DA, Holmes S. Comparisons of distance methods for combining covariates and abundances in microbiome studies. Pac Symp Biocomput. 2012:213–24.Google Scholar
- Lozupone CA, Hamady M, Kelley ST, Knight R. Quantitative and qualitative diversity measures lead to different insights into factors that structure microbial communities. Appl Environ Microbiol. 2007;73:1576–85.View ArticlePubMedPubMed CentralGoogle Scholar
- Cuesta-Gerveno JM, Risco-Perez D, Goncalves-Blanco P, Garcia-Jimenez WL, Gil-Molino M, Fernandez-Llario P, et al. Fatal infection due to Haemophilus parasuis in a young wild boar (Sus scrofa). J Vet Diagn Invest. 2013;25:297–300.View ArticlePubMedGoogle Scholar
- Aragon V, Cerda-Cuellar M, Fraile L, Mombarg M, Nofrarias M, Olvera A, Sibila M, Solanes D, Segales J. Correlation between clinico-pathological outcome and typing of Haemophilus parasuis field strains. Vet Microbiol. 2010;142:387–93.View ArticlePubMedGoogle Scholar
- Lowe BA, Marsh TL, Isaacs-Cosgrove N, Kirkwood RN, Kiupel M, Mulks MH. Microbial communities in the tonsils of healthy pigs. Vet Microbiol. 2011;147:346–57.View ArticlePubMedGoogle Scholar
- Weese JS, Slifierz M, Jalali M, Friendship R. Evaluation of the nasal microbiota in slaughter-age pigs and the impact on nasal methicillin-resistant Staphylococcus aureus (MRSA) carriage. BMC Vet Res. 2014;10:1–10.View ArticleGoogle Scholar
- Carding S, Verbeke K, Vipond DT, Corfe BM, Owen LJ. Dysbiosis of the gut microbiota in disease. Microb Ecol Health Dis. 2015;26:26191.PubMedGoogle Scholar
- Aragon V, Segales J, Oliveira S. Glässer’s disease. In: Diseases of swine. 10th ed. Iowa: Wiley-Blackwell; 2012. p. 760–69.Google Scholar
- Kim HB, Isaacson RE. The pig gut microbial diversity: understanding the pig gut microbial ecology through the next generation high throughput sequencing. Vet Microbiol. 2015;177:242–51.View ArticlePubMedGoogle Scholar
- Keesing F, Ostfeld RS. Ecology. Is biodiversity good for your health? Science. 2015;349:235–6.View ArticlePubMedGoogle Scholar
- Lozupone CA, Stombaugh JI, Gordon JI, Jansson JK, Knight R. Diversity, stability and resilience of the human gut microbiota. Nature. 2012;489:220–30.View ArticlePubMedPubMed CentralGoogle Scholar
- Nylund L, Nermes M, Isolauri E, Salminen S, de Vos WM, Satokari R. Severity of atopic disease inversely correlates with intestinal microbiota diversity and butyrate-producing bacteria. Allergy. 2015;70:241–4.View ArticlePubMedGoogle Scholar
- Coburn B, Wang PW, Diaz Caballero J, Clark ST, Brahma V, Donaldson S, et al. Lung microbiota across age and disease stage in cystic fibrosis. Sci Rep. 2015;5:10241.View ArticlePubMedPubMed CentralGoogle Scholar
- McCoy KD, Koller Y. New developments providing mechanistic insight into the impact of the microbiota on allergic disease. Clin Immunol. 2015;159:170–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Olvera A, Pina S, Macedo N, Oliveira S, Aragon V, Bensaid A. Identification of potentially virulent strains of Haemophilus parasuis using a multiplex PCR for virulence-associated autotransporters (vtaA). Vet J. 2012;191:213–8.View ArticlePubMedGoogle Scholar
- Wertheim HF, Vos MC, Ott A, van Belkum A, Voss A, Kluytmans JA, et al. Risk and outcome of nosocomial Staphylococcus aureus bacteraemia in nasal carriers versus non-carriers. Lancet. 2004;364:703–5.View ArticlePubMedGoogle Scholar
- Schmidt B, Mulder IE, Musk CC, Aminov RI, Lewis M, Stokes CR, et al. Establishment of normal gut microbiota is compromised under excessive hygiene conditions. Plos One. 2011. doi:10.1371/journal.pone.0028284.Google Scholar
- Roussos A, Philippou N, Mantzaris GJ, Gourgoulianis KI. Respiratory diseases and Helicobacter pylori infection: is there a link? Respiration. 2006;73:708–14.View ArticlePubMedGoogle Scholar
- Ichinohe T, Pang IK, Kumamoto Y, Peaper DR, Ho JH, Murray TS, Iwasaki A. Microbiota regulates immune defense against respiratory tract influenza A virus infection. Proc Natl Acad Sci. 2011;108:5354–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Scott KP, Gratz SW, Sheridan PO, Flint HJ, Duncan SH. The influence of diet on the gut microbiota. Pharmacol Res. 2013;69:52–60.View ArticlePubMedGoogle Scholar
- Borcard D, Legendre P, Drapeau P. Partialling out the spatial component of ecological variation. Ecology. 1992;73:1045–55.View ArticleGoogle Scholar
- White M. Early medications and respiratory disease in growing pigs. 2014. Available at https://www.pig333.com/clinical-case-of-the-world/early-medications-and-respiratory-disease-in-growing-pigs_9074/.Google Scholar
- O’Hara AM, Shanahan F. The gut flora as a forgotten organ. EMBO Rep. 2006;7:688–93.View ArticlePubMedPubMed CentralGoogle Scholar
- Klindworth A, Pruesse E, Schweer T, Peplies J, Quast C, Horn M, Glockner FO. Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. 2012;41:e1.View ArticlePubMedPubMed CentralGoogle Scholar
- Bokulich NA, Subramanian S, Faith JJ, Gevers D, Gordon JI, Knight R, et al. Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing. Nat Meth. 2013;10:57–9.View ArticleGoogle Scholar
- Aronesty E. Ea utils: Command-line tools for precessing biological sequencing data. 2011. Available at http://code.google.com/p/ea-utils.Google Scholar
- Aronesty E. TOBioJ: Comparison of sequencing Utility Programs. 2013; doi:10.2174/18750362013070100001.
- Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26:2460–1.View ArticlePubMedGoogle Scholar
- DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, et al. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006;72:5069–72.View ArticlePubMedPubMed CentralGoogle Scholar
- Cole JR, Wang Q, Fish JA, Chai B, McGarrell DM, Sun Y, et al. Ribosomal DATABASE PROject: data and tools for high throughput rRNA analysis. Nucleic Acids Res. 2013;42:D633–42.View ArticlePubMedPubMed CentralGoogle Scholar
- Lozupone C, Knight R. UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol. 2005;71:8228–35.View ArticlePubMedPubMed CentralGoogle Scholar