Next Article in Journal
Correlating Genotyping Data of Coxiella burnetii with Genomic Groups
Next Article in Special Issue
Sow Contact Is a Major Driver in the Development of the Nasal Microbiota of Piglets
Previous Article in Journal
The Novel Zoonotic Pathogen, Anaplasma capra, Infects Human Erythrocytes, HL-60, and TF-1 Cells In Vitro
Previous Article in Special Issue
Antiviral Mechanism of Tea Polyphenols against Porcine Reproductive and Respiratory Syndrome Virus
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Altered Nasal Microbiota Composition Associated with Development of Polyserositis by Mycoplasma hyorhinis

1
IRTA, Centre de Recerca en Sanitat Animal (CReSA, IRTA-UAB), Campus de la Universitat Autònoma de Barcelona, 08193 Bellaterra, Spain
2
OIE Collaborating Centre for the Research and Control of Emerging and Re-Emerging Swine Diseases in Europe (IRTA-CReSA), 08193 Bellaterra, Spain
3
Departamento de Ciencia Animal, Escuela Técnica Superior de Ingeniería Agraria (ETSEA), Universidad de Lleida, 25198 Lleida, Spain
*
Author to whom correspondence should be addressed.
Pathogens 2021, 10(5), 603; https://doi.org/10.3390/pathogens10050603
Submission received: 1 April 2021 / Revised: 10 May 2021 / Accepted: 12 May 2021 / Published: 14 May 2021
(This article belongs to the Special Issue Microbiota and Respiratory Diseases in Pigs)

Abstract

:
Fibrinous polyserositis in swine farming is a common pathological finding in nursery animals. The differential diagnosis of this finding should include Glaesserella parasuis (aetiological agent of Glässer’s disease) and Mycoplasma hyorhinis, among others. These microorganisms are early colonizers of the upper respiratory tract of piglets. The composition of the nasal microbiota at weaning was shown to constitute a predisposing factor for the development of Glässer’s disease. Here, we unravel the role of the nasal microbiota in the subsequent systemic infection by M. hyorhinis, and the similarities and differences with Glässer’s disease. Nasal samples from farms with recurrent problems with polyserositis associated with M. hyorhinis (MH) or Glässer’s disease (GD) were included in this study, together with healthy control farms (HC). Nasal swabs were taken from piglets in MH farms at weaning, before the onset of the clinical outbreaks, and were submitted to 16S rRNA gene amplicon sequencing (V3–V4 region). These sequences were analyzed together with sequences from similar samples previously obtained in GD and HC farms. Animals from farms with disease (MH and GD) had a nasal microbiota with lower diversity than those from the HC farms. However, the composition of the nasal microbiota of the piglets from these disease farms was different, suggesting that divergent microbiota imbalances may predispose the animals to the two systemic infections. We also found variants of the pathogens that were associated with the farms with the corresponding disease, highlighting the importance of studying the microbiome at strain-level resolution.

1. Introduction

The microbiota is defined as the community of microorganisms found in the surface of a tissue, which represents their usual ecological niche [1]. Several studies have found an association between the microbiota composition and phenotypic features from the host, such as body weight, health status and disease onset, among others [2]. In animals, and in particular swine, the gut microbiota has been deeply studied in contrast to other body sites. The porcine intestinal microbiota influences several production parameters like body weight [3] and feed efficiency [4]. The alterations in its composition can affect the translocation of metabolites across a disrupted intestinal barrier, which in turn could promote metabolic pathologies in various organs, such as liver and adipose tissue [4]. Although to a lesser extent, the nasal microbiota has also been associated with the development of diseases [5]. In swine, the nasal microbiota has been suggested as a predisposing factor for the development of Glässer’s disease, a systemic infection caused by Glaesserella (Haemophilus) parasuis and characterized by fibrinous polyserositis [6].
Polyserositis is a common finding in necropsies during the post-weaning period. This pathology is not only caused by G. parasuis, but also by other bacteria, such as Mycoplasma hyorhinis [7]. These microorganisms are considered early colonizers of the respiratory tract, as they are part of the bacterial communities that conform the normal microbiota of the nasal cavity of piglets early in life [6,8,9,10]. How these etiological agents switch from members of the normal microbiota to disseminate and induce systemic infection is the key to understand the pathogenesis of these diseases. A balance between colonization by these potential pathogens and the host is established during the early stages of the piglets’ life [11] where the immune system plays a crucial role to establish this balance and maintain health. Recently, the role of particular taxa from the piglet microbiota has been shown in the polyserositis caused by G. parasuis where Lactobacillus and Prevotella were associated with health, while Moraxella, Haemophilus (Glaesserella) and Streptococcus were associated with Glässer’s disease [6,12]. Similarly, differences in the oropharyngeal microbiota have been associated with the development of respiratory diseases in pigs [13]. Higher relative abundance of the Moraxella genus was associated with respiratory pathology and Lactobacillus was associated with healthy animals [13,14].
Here, we compare the piglet nasal microbiota composition at weaning, in farms with and without polyserositis cases caused by M. hyorhinis. Furthermore, we compare the nasal microbiota composition from a M. hyorhinis-affected farm, with data from previously characterized farms with recurrent outbreaks of Glässer’s disease, and divergent changes from health between both pathological scenarios.

2. Results

2.1. Alpha Diversity Is Significantly Reduced in the Nasal Microbiota of Piglets from Farms with Polyserositis

Farms with post-weaning polyserositis caused by M. hyorhinis (MH group) or G. parasuis (GD group), together with healthy control (HC group) farms, were included in this study (Table 1). Nasal samples were obtained from piglets at weaning and the nasal microbiota was determined by 16S rRNA gene sequencing. We obtained a total of 33,012,373 reads from the sequencing data of 74 nasal samples. The minimum number of reads in a sample was 70,970 and the maximum was 121,478. After quality-control, denoising and chimera removal, a total of 5,061,172 sequences were included in the posterior analysis.
The feature table was rarefied to 6,545 sequences per sample and alpha diversity was estimated using different metrics. Microbial richness was assessed by calculating the observed features, which gives the count of the different amplicon sequence variants (ASVs) in rarefied samples. We found a total of 25,993 ASVs, with significant differences in richness among the three farm groups (Kruskal-Wallis; p = 0.027). Richness was significantly higher in the HC farms when compared with the farms with polyserositis MH and GD individually (Figure 1A). Chao index showed the same trend and detected the same differences between the groups (not shown). Alpha diversity was additionally assessed by Shannon index (Figure 1B), which considers both the richness and the evenness of the samples, and also detected differences among the three groups (Kruskal-Wallis; p = 0.017). Farms with disease showed lower diversity than HC group when compared individually (HC vs MH and HC vs GD, p = 0.01 in both cases; Figure 1B).

2.2. Piglets from Farms with Polyserositis by M. hyorhinis Showed Significatively Different Nasal Microbiota Composition Compared to Piglets from HC and GD Farms

The differences in composition (beta diversity) of the nasal microbial communities among the different study groups were explored using Principal Coordinates Analysis (PCoA). The beta diversity analysis was done using Bray Curtis distance to measure the compositional dissimilarity quantitatively between different groups. The PCoA plot showed a more similar nasal microbiota composition between animals from GD and HC farms compared with the MH farms (Figure 2A). Although strong clustering was observed when the samples were compared considering the farm of origin (R2 = 47% estimated by Adonis function), the differences among the three groups of farms were explained by 17.3% (R2, p = 0.001, 999 permutations). Moreover, to understand the qualitative dissimilarities between groups the Jaccard distances were calculated. The Adonis function on this qualitative approach retrieved values for R2 of 26.12% and 9.52% for farm or group clustering, respectively (Supplementary Figure S1A).
We also measured the qualitative (unweighted) and quantitative (weighted) measures of community dissimilarity incorporating the phylogenetic relationships between features (Unifrac). The percentage explained by the clustering of these groups was higher in the unweighted analysis (R2 = 12.07% Adonis, p = 0.01; Figure 2B) than in the weighted analysis (R2 = 10.04% Adonis, p < 0.01; Supplementary Figure S1B). The farm of origin had also a strong effect in the clustering (R2 = 36.6% Adonis, p = 0.01).

2.3. Taxonomic Assignment of the ASVs

Taxa assignment was done by Naive Bayes classification algorithm using the GreenGenes 13.8 16S gene database. Taxonomy was assigned to a total of 25,993 ASVs. The percentage of unassigned ASVs was increasing from phylum, with 0.96% of unassigned taxa, to lower taxonomic levels until species where 89.44% were found.
At phylum level, more than the 80% of the relative abundance were shared between three phyla, Proteobacteria, Bacteroidetes and Firmicutes. Proteobacteria was the most relative abundant phylum in the disease groups, while Bacteroidetes and Firmicutes were the most relatively abundant in the healthy group. Gammaproteobacteria was the most relative abundant bacterial class of the microbiota composition in the disease groups (MH and GD), while Clostridia was the most abundant in the healthy one. Indeed, Clostridiales and Pseudomonales were the most represented orders in the HC and disease groups, respectively.
At the family level, the most abundant were Moraxellaceae, Weeksellaceae, Pasteurellaceae, Ruminococcaceae, Lachnospiraceae, Prevotellaceae, Muribaculaceae (S24-7), Streptococcaceae and Mycoplasmataceae. Moraxellaceae was the most relative abundant family in all groups, with 17.69% in the HC group, 26.09% in GD and 32.30% in the MH group.
The most abundant genus in all groups was Bergeyella, with 14.45% in the HC group, 17.26% in GD and 20.85% in the MH group. Enhydrobacter was the second most abundant genus in HC and MH group with 6.30% and 18.19%, respectively. Moraxella was the second most abundant genus in GD (9.73%) and the third one in the MH group with the same percentage of 9.73. Finally, an unclassified genus from Moraxellaceae was the third most abundant in HC (5.76%) and GD (9.39%).
The 30 most abundant ASVs across all the samples belonged to Moraxellaceae (n = 4), Moraxella (n = 2), Enhydrobacter (n = 7), G. (Haemophilus) parasuis (n = 1) and Bergeyella zoohelcum (n = 17), which were distinctly distributed depending on the group (Figure 3). If we focus on the ASVs from G. parasuis and M. hyorhinis, the aethiological agents of the polyserositis observed in the farms, 45 ASVs were classified as M. hyorhinis (Figure 4A) and 42 as G. parasuis (Figure 4B). Curiously, the ASVs from these pathogens were not shared among the groups, but they were shared among different farms from the same group (Figure 4C,D). Interestingly, the amplicon 16S sequences amplified from the isolated clinical samples were highly similar (99–100%) to the two most abundant M. hyorhinis ASVs in the MH group clustering together into the phylogenetic tree (Figure 4A, strong yellow).

Differential Abundance Analysis

To find out taxa differentially abundant among different groups we used the Analysis of composition of Microbiome (ANCOM) method. When the nasal composition from the three groups was compared, 144 differential ASVs were detected (Supplementary Table S7). The most significant ASVs belonged to Moraxellaceae, Bergeyella zoohelcum, Enhydrobacter, Moraxella and Rothia nasimurium. Remarkably, most of these ASVs detected as differential were present with a high relative abundance in the MH group when compared with the other two groups. Two of the high-abundant ASVs from MH group were classified as M. hyorhinis, depicted in Figure 4A in bright yellow. On the other hand, one ASV was classified as G. parasuis but was only present in the MH group. Moreover, nine of the differentially abundant ASVs were found among the top 30 most abundant among all the groups globally, which are indicated with an asterisk in Figure 3.
To explore if there were ASVs associated to the health status, we analyzed groups MH and GD joined as a “disease group”, against control (HC). Only eight features appeared as differentially abundant between these two groups, including three from Prevotella genus and one from Streptococcus alactolyticus, all in higher abundance in healthy farms (Supplementary Table S1). Additionally, when GD group was compared to the HC, only one ASVs from the Lachnospiraceae family was detected as differentially abundant. However, when MH was compared to HC, 67 ASVs were detected (Supplementary Table S1), where the ten most significant ASVs included seven classified as B. zoohelcum, two as Enhydrobacter and one as Moraxellaceae family. These latter ASV was the only one (from these ten) associated with the control farms (absent in MH farms), since the seven ASVs from B. zoohelcum were only present in the MH group, (with relative abundances ranging from 1.73% to 3.51%) and the two ASVs classified as Enhydrobacter were absent in HC farms (Supplementary Table S1).
The differences observed at the ASVs level were not reflected in the differential analysis (ANCOM) derived from the collapsed ASVs table at the different taxonomic levels (Supplementary Tables S2–S7).

3. Discussion

Culture-independent studies for interpreting the swine microbiome in disease infections can be challenging, especially if the study aims to elucidate the differences in the microbiome before the onset of the clinical signs as a predisposing factor. Here, we studied the relationship between the nasal microbiota and the later development of the systemic infection by the early colonizers of the respiratory tract, M. hyorhinis and G. parasuis. We found lower alpha diversity in disease farms, with the nasal microbiota from piglets from farms with disease by M. hyorhinis showing higher divergence from the healthy farms than those with Glässer’s disease. As it has been previously reported, imbalances in richness and/or diversity in different microbiome landscapes are associated with changes that generally may lead to disease in the host organism [15]. In swine production, changes in the nasal microbiota of the weaning piglets may influence the subsequent development of Glässer’s disease [6].
Antimicrobial usage is one of the plethora of factors that can affect the pig microbiome in commercial farms. The animals included in this study underwent different antimicrobial treatments, which may impact the nasal microbiota composition. However, we found that different farms under the same antimicrobial treatment but with different health status (GE, RM, GM) did not cluster together, suggesting that these treatments did not seem to be determinant in the subsequent health status of the animals. Although the present study presents some limitations due to the different factors affecting the microbiome composition of the animals in commercial farms, our results reinforce the general idea that lower alpha diversity is related with disease development [6,16]. Lower bacterial richness and evenness was found in the nasal microbiota in the farms with disease, when they were compared with the animals from the healthy control farms. This finding is supported by many other studies that established that lower alpha diversity values are linked to poor health status in pigs, but also in other species [13,16,17]. Here, a lower microbial richness of the nasal cavity in the disease farms may be associated with the proliferation of the undesirable taxonomic groups that would lead to systemic infections. However, the nasal microbiota composition from the two different disease groups differed, which indicates that different microbiota imbalances may predispose the animals to different infections. In general, the composition from the farms with Glässer’s disease were more similar to the healthy ones but different from the M. hyorhinis farms. At higher taxonomic levels (phylum to order), highly similar composition between the disease groups MH and GD was observed, while at lower levels (family to species) differences between these two disease groups were evident. In fact, the distribution of the most abundant ASVs showed a clear different pattern in MH and GD farms. Nevertheless, we were able to find some taxa associated with health, such as the Prevotellaceae family, which deserves further study. This family was also represented in the healthy farms at the ASVs level by Prevotella copri and Prevotella stercorea, two species previously found in higher abundance in suckling and weaning healthy piglets [18]. Two ASVs from the Lachnospiraceae family were also associated with healthy farms, in agreement with previous studies, where ASVs from this family were correlated with higher feed conversion ratio [19].
Animals from control farms showed a lower relative abundance of M. hyorhinis and G. parasuis. However, the relative abundance of the pathogens at the species level was not related to the posterior development of the corresponding systemic disease. Surprisingly, G. parasuis abundance was higher in MH farms, while M. hyorhinis was higher in GD farms. Although we cannot rule out a synergistic association between these two pathogens to develop disease, as it has been described before [20], our data seems to support the role of specific ASVs of each pathogen in disease development. It is well known that G. parasuis comprises strains of different pathogenic potential with different consequences for the piglet health [21,22]. In fact, colonization by virulent G. parasuis strains increases the risk of developing Glässer’s disease, while the non-virulent strains can provide some protection against the disease [22]. Although there are some suggestions in the literature of the existence of M. hyorhinis strains with different degree of virulence [23], and our data also support it, this has not been demonstrated. We detected specific ASVs of the pathogens that were found only in the disease farms and were not shared by the farms with different disease. The role of these specific M. hyorhinis ASVs in disease was supported by their detection in the clinical samples from those farms. Due to the limitations of the 16S rRNA gene amplicon sequencing, more studies with whole genome sequencing are needed. This result also reinforces the importance of study the microbiota at strain level, to better understand the intrinsic characteristics of the different strains and their role in the predisposition to the systemic infection.
In summary, different changes in the nasal microbiota composition were observed in weaning piglets from farms with polyserositis caused by either M. hyorhinis or G. parasuis. We hypothesize that these changes might facilitate dissemination and the subsequent development of the systemic infections by M. hyorhinis and G. parasuis. The strain level resolution of the microbiota should be studied for virulent-strain detection, especially in the case of M. hyorhinis where there is a lack of studies addressing this issue.

4. Materials and Methods

4.1. Farm Selection and Sampling

Farms were selected based on the presence/absence of post-weaning polyserositis cases to study the nasal microbiota. Three groups of farms were included in the study, (Table 1). The first group, MH, was formed by two farms where polyserositis cases in nursery pigs were ascribed to M. hyorhinis. For such a purpose, samples from polyarthritis and/or polyserositis lesions were taken from animals at necropsy and M. hyorhinis was detected by qPCR [24] and/or isolation. In addition, four of the positive clinical samples (2 from RM and 2 from GE farms) were submitted at Servei de Genòmica, Universitat Autònoma de Barcelona to sequence the 16S rRNA gene using Sanger technology. All the samples from necropsies from MH farms were negative to isolation and/or PCR of G. parasuis [25] and S. suis [26]. In each MH farm, nasal swabs were taken from ten 3–4 week-old piglets selected from different litters (two piglets per sow) in order to avoid a sow effect bias. The second group was formed by three control farms without polyserositis or respiratory cases in the last two years previous to the study (Group HC). The third group of farms, group GD, was composed of four farms in which the polyserositis problems were attributed to G. parasuis. Raw sequencing data from groups HC and GD belonged to a previous study (SRA Accession number SRP068182 [6]).
Sampling of piglets was done under institutional authorization and followed good veterinary practices. According to European (Directive 2010/63/EU) and Spanish (Real Decreto 53/2013) normative, this procedure did not require specific approval by an Ethical Committee (Chapter I, Article 1, 5 (f) of 2010/63/EU).

4.2. DNA Extraction and 16S rRNA Gene Amplicon Sequencing

Nasal swabs were placed in 500 μL de PBS, vortex for 30 s and stored at −80 °C until used for DNA extraction. DNA was extracted from 200 μL of the initial 500 μL PBS where the swabs were resuspended and eluted in 100 μL of PBS using the Nucleospin Blood (Macherey Nagel, Bethlehem, PA, USA) kit. Quantity and quality assessment of the DNA was performed using BioDrop DUO (BioDrop Ltd., Cambridge, UK). Samples were submitted for 16S rRNA gene amplicon sequencing using the Illumina paired-end 2 × 250 bp kit (MS-102-2003 MiSeq® Reagent Kit v2, 500 cycle) following the manufacturer’s instructions. The library preparation for sequencing was performed within 24 h after the DNA extraction at Servei de Genòmica, Universitat Autònoma de Barcelona. The region amplified was V3–V4 that covers two hypervariable regions of the conserved gene, and it was sequenced according to the Illumina protocol [27] as previously described [6]. The entire sequence dataset is available at the NCBI database, SRA accession number PRJNA717778.

4.3. Microbiota In-Silico Analysis

The in-silico analysis was done using the plugin based software Quantitative Insights into Microbial Ecology (QIIME) vs 2020.11 [28]. First, after importing the raw sequencing data into Qiime2 (q2-import), we performed a quality check step and decided to trim the reads to a length of 240 bp based on the quality drop at the end of the reads. Denoising, trimming and quality-based filtering was performed with the q2-dada2 plugin [29], which also includes a chimera detection and consequent removal. Moreover, to improve the quality of our dataset, unassigned taxa was filtered out after aligning the ASVs against the Greengenes reference database (vs. 13.8) [30], clustered at 88% identity with 65% of identity and over the 50% of query [31].
With the rarefied curves of the richness, we extracted the proper sample-depth as a key parameter for the following diversity calculations. Diversity within each sample (alpha diversity) was estimated with Shannon’s entropy index, Observed Features and Chao index using q2-core-metrics plugin [32,33]. To compare the microbiota composition among the study groups (beta diversity), we calculated Jaccard, Bray Curtis and Unifrac [34,35,36] (weighted and unweighted) metrics and represented in the spatial coordinates with the Principal Components Analysis (PCoA). PERMANOVA [37] tests were performed to analyze beta diversity clusters among study groups. Also, the Adonis function [38] from Vegan package [39] was used across the matrix distances in order to determine the percentage of explanation from the grouping variables analyzed.
Taxonomic assignment of the ASVs was performed by a pre-trained Naive Bayes classifier within the q2-feature-classifier classify-sklearn plugin [40] using the Greengenes database (Vs. 13.8) [30]. Classification of ASVs from Glaesserella and Mycoplasma genera was confirmed to species level using Basic Local Alignment Search Tool (Blast) [41] algorithm and NCBI RefSeq database [42]. To find differentially abundant features across groups, analysis of compositional microbiomes (ANCOM) [43] was performed. To process and analyze the data from the taxa assignment of the ASVs and the counts of the features we used qiime2r (18), phyloseq [44], tidyverse [45] R packages. Heatmap, PCoA and Venn diagrams plots were built in Rstudio [46] using ggplot2 [47], ComplexHeatmap [48] and ggVennDiagram [49] packages. The phylogenetic tree of the ASVs from M. hyorhinis and G. parasuis were built with Maximum likelihood algorithm using Qiime2 [28].

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/pathogens10050603/s1, Figure S1: Beta diversity analysis on Jaccard (A) and weighted Unifrac distances (B) of the nasal microbiota of weaning piglets. Each circle represents the microbial composition of each sample and the size of the circles is proportional to the number of Observed Features as indicated in the legend. Healthy control farms (HC) are depicted in blue, while farms with polyserositis caused by M. hyorhinis (MH) are in orange and G. parasuis (GD) in green palettes. Ellipses are calculated with the euclidean distances of the grouped samples. Table S1: Differential abundance analysis (ANCOM) comparing the nasal microbiota of weaning piglets from the three different group of farms (GD, farms with Glässer’s disease; MH, farms with polyserositis caused by M. hyorhinis and HC, healthy control farms) at ASVs level. Different comparisons are shown in the different excel sheets as indicated in each sheet label. Table S2: Differential abundance analysis (ANCOM) comparing the nasal microbiota of weaning piglets from the three different group of farms (GD, farms with Glässer’s disease; MH, farms with polyserositis caused by M. hyorhinis and HC, healthy control farms) at phylum level. Different comparisons are shown in the different excel sheets as indicated in each sheet label. Table S3: Differential abundance analysis (ANCOM) comparing the nasal microbiota of weaning piglets from the three different group of farms (GD, farms with Glässer’s disease; MH, farms with polyserositis caused by M. hyorhinis and HC, healthy control farms) at class level. Different comparisons are shown in the different excel sheets as indicated in each sheet label. Table S4: Differential abundance analysis (ANCOM) comparing the nasal microbiota of weaning piglets from the three different group of farms (GD, farms with Glässer’s disease; MH, farms with polyserositis caused by M. hyorhinis and HC, healthy control farms) at order level. Different comparisons are shown in the different excel sheets as indicated in each sheet label. Table S5: Differential abundance analysis (ANCOM) comparing the nasal microbiota of weaning piglets from the three different group of farms (GD, farms with Glässer’s disease; MH, farms with polyserositis caused by M. hyorhinis and HC, healthy control farms) at family level. Different comparisons are shown in the different excel sheets as indicated in each sheet label. Table S6: Differential abundance analysis (ANCOM) comparing the nasal microbiota of weaning piglets from the three different group of farms (GD, farms with Glässer’s disease; MH, farms with polyserositis caused by M. hyorhinis and HC, healthy control farms) at genus level. Different comparisons are shown in the different excel sheets as indicated in each sheet label. Table S7: Differential abundance analysis (ANCOM) comparing the nasal microbiota of weaning piglets from the three different group of farms (GD, farms with Glässer’s disease; MH, farms with polyserositis caused by M. hyorhinis and HC, healthy control farms) at species level. Different comparisons are shown in the different excel sheets as indicated in each sheet label.

Author Contributions

Conceptualization, F.C.-F., M.S. and V.A.; methodology, M.B.-F. and L.F.; software, M.B.-F.; validation, F.C.-F.; formal analysis, M.B.-F. and F.C.-F.; writing—original draft preparation, M.B.-F.; writing—review and editing, F.C.-F., L.F., M.S., V.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Economy and Competitiveness (MINECO) of the Spanish Government (grant number AGL2016-77361). MBF is a recipient of an FPI fellowship, also from MINECO (Spain). The authors are also grateful to the Centres de Recerca de Catalunya (CERCA) Programme.

Institutional Review Board Statement

Ethical review and approval were waived for this study, following the European (Directive 2010/63/EU Chapter I, Article 1, 5 (f) of 2010/63/EU) and Spanish (Real Decreto 53/2013) normative.

Informed Consent Statement

Not applicable.

Data Availability Statement

The entire sequence dataset is available at the NCBI database, SRA accession number PRJNA717778.

Acknowledgments

The authors want to acknowledge the technical laboratory work of Eva Huerta and Núria Galofré-Milà, both from IRTA-CReSA.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. NIH HMP Working Group; Peterson, J.; Garges, S.; Giovanni, M.; McInnes, P.; Wang, L.; Schloss, J.A.; Bonazzi, V.; McEwen, J.E.; Wetterstrand, K.A.; et al. The NIH Human Microbiome Project. Genome Res. 2009, 19, 2317–2323. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Zheng, D.; Liwinski, T.; Elinav, E. Interaction between microbiota and immunity in health and disease. Cell Res. 2020, 30, 492–506. [Google Scholar] [CrossRef] [PubMed]
  3. Han, G.G.; Lee, J.-Y.; Jin, G.-D.; Park, J.; Choi, Y.H.; Chae, B.J.; Kim, E.B.; Choi, Y.-J. Evaluating the association between body weight and the intestinal microbiota of weaned piglets via 16S rRNA Sequencing. Appl. Microbiol. Biotechnol. 2017, 101, 5903–5911. [Google Scholar] [CrossRef] [PubMed]
  4. McCormack, U.M.; Curião, T.; Buzoianu, S.G.; Prieto, M.L.; Ryan, T.; Varley, P.; Crispie, F.; Magowan, E.; Metzler-Zebeli, B.U.; Berry, D.; et al. Exploring a Possible Link between the Intestinal Microbiota and Feed Efficiency in Pigs. Appl. Environ. Microbiol. 2017, 83. [Google Scholar] [CrossRef] [Green Version]
  5. Man, W.H.; de Steenhuijsen Piters, W.A.A.; Bogaert, D. The microbiota of the respiratory tract: Gatekeeper to respiratory health. Nat. Rev. Microbiol. 2017, 15, 259–270. [Google Scholar] [CrossRef] [PubMed]
  6. Correa-Fiz, F.; Fraile, L.; Aragon, V. Piglet nasal microbiota at weaning may influence the development of Glässer’s disease during the rearing period. BMC Genom. 2016, 17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Mycoplasma Hyorhinis—An Overview | ScienceDirect Topics. Available online: https://0-www-sciencedirect-com.brum.beds.ac.uk/topics/agricultural-and-biological-sciences/mycoplasma-hyorhinis (accessed on 15 March 2021).
  8. Clavijo, M.J.; Davies, P.; Morrison, R.; Bruner, L.; Olson, S.; Rosey, E.; Rovira, A. Temporal patterns of colonization and infection with Mycoplasma hyorhinis in two swine production systems in the USA. Vet. Microbiol. 2019, 234, 110–118. [Google Scholar] [CrossRef]
  9. Roos, L.R.; Surendran Nair, M.; Rendahl, A.K.; Pieters, M. Mycoplasma hyorhinis and Mycoplasma hyosynoviae Dual Detection Patterns in Dams and Piglets. PLoS ONE 2019, 14, e0209975. [Google Scholar] [CrossRef] [Green Version]
  10. Cerdà-Cuéllar, M.; Naranjo, J.F.; Verge, A.; Nofrarías, M.; Cortey, M.; Olvera, A.; Segalés, J.; Aragon, V. Sow Vaccination Modulates the Colonization of Piglets by Haemophilus parasuis. Vet. Microbiol. 2010, 145, 315–320. [Google Scholar] [CrossRef]
  11. Costa-Hurtado, M.; Barba-Vidal, E.; Maldonado, J.; Aragon, V. Update on Glässer’s Disease: How to Control the Disease under Restrictive Use of Antimicrobials. Vet. Microbiol. 2020, 242, 108595. [Google Scholar] [CrossRef]
  12. Mahmmod, Y.S.; Correa-Fiz, F.; Aragon, V. Variations in association of nasal microbiota with virulent and non-virulent strains of Glaesserella (Haemophilus) parasuis in weaning piglets. Vet. Res. 2020, 51, 7. [Google Scholar] [CrossRef] [Green Version]
  13. Wang, Q.; Cai, R.; Huang, A.; Wang, X.; Qu, W.; Shi, L.; Li, C.; Yan, H. Comparison of oropharyngeal microbiota in healthy piglets and piglets with respiratory disease. Front. Microbiol. 2018, 9. [Google Scholar] [CrossRef] [Green Version]
  14. Valeris-Chacin, R.; Sponheim, A.; Fano, E.; Isaacson, R.; Singer, R.S.; Nerem, J.; Leite, F.L.; Pieters, M. Relationships among fecal, air, oral, and tracheal microbial communities in pigs in a respiratory infection disease model. Microorganisms 2021, 9, 252. [Google Scholar] [CrossRef]
  15. Ojetti, V.; Gigante, G.; Ainora, M.E.; Fiore, F.; Barbaro, F.; Gasbarrini, A. Microflora imbalance and gastrointestinal diseases. Dig. Liver Dis. Suppl. 2009, 3, 35–39. [Google Scholar] [CrossRef]
  16. Prehn-Kristensen, A.; Zimmermann, A.; Tittmann, L.; Lieb, W.; Schreiber, S.; Baving, L.; Fischer, A. Reduced microbiome alpha diversity in young patients with ADHD. PLoS ONE 2018, 13. [Google Scholar] [CrossRef]
  17. Pirolo, M.; Espinosa-Gongora, C.; Bogaert, D.; Guardabassi, L. The porcine respiratory microbiome: Recent insights and future challenges. Anim. Microbiome 2021, 3, 9. [Google Scholar] [CrossRef]
  18. Amat, S.; Lantz, H.; Munyaka, P.M.; Willing, B.P. Prevotella in Pigs: The positive and negative associations with production and health. Microorganisms 2020, 8, 1584. [Google Scholar] [CrossRef]
  19. Quan, J.; Cai, G.; Ye, J.; Yang, M.; Ding, R.; Wang, X.; Zheng, E.; Fu, D.; Li, S.; Zhou, S.; et al. A global comparison of the microbiome compositions of three gut Locations in commercial pigs with extreme feed conversion ratios. Sci. Rep. 2018, 8. [Google Scholar] [CrossRef] [Green Version]
  20. Palzer, A.; Haedke, K.; Heinritzi, K.; Zoels, S.; Ladinig, A.; Ritzmann, M. Associations among Haemophilus parasuis, Mycoplasma hyorhinis, and porcine reproductive and respiratory syndrome virus infections in pigs with polyserositis. Can. Vet. J. 2015, 56, 285–287. [Google Scholar]
  21. Galofré-Milà, N.; Correa-Fiz, F.; Lacouture, S.; Gottschalk, M.; Strutzberg-Minder, K.; Bensaid, A.; Pina-Pedrero, S.; Aragon, V. A robust PCR for the differentiation of potential virulent strains of Haemophilus parasuis. BMC Vet. Res. 2017, 13, 124. [Google Scholar] [CrossRef] [Green Version]
  22. Brockmeier, S.L.; Loving, C.L.; Mullins, M.A.; Register, K.B.; Nicholson, T.L.; Wiseman, B.S.; Baker, R.B.; Kehrli, M.E. Virulence, transmission, and heterologous protection of four isolates of Haemophilus parasuis. Clin. Vaccine Immunol. 2013, 20, 1466–1472. [Google Scholar] [CrossRef] [Green Version]
  23. Lin, J.H.; Chen, S.P.; Yeh, K.S.; Weng, C.N. Mycoplasma hyorhinis in Taiwan: Diagnosis and isolation of swine pneumonia pathogen. Vet. Microbiol. 2006, 115, 111–116. [Google Scholar] [CrossRef]
  24. Clavijo, M.J.; Oliveira, S.; Zimmerman, J.; Rendahl, A.; Rovira, A. Field evaluation of a quantitative polymerase chain reaction assay for Mycoplasma hyorhinis. J. Vet. Diagn Investig. 2014, 26, 755–760. [Google Scholar] [CrossRef] [Green Version]
  25. 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–218. [Google Scholar] [CrossRef]
  26. Ishida, S.; Tien, L.H.T.; Osawa, R.; Tohya, M.; Nomoto, R.; Kawamura, Y.; Takahashi, T.; Kikuchi, N.; Kikuchi, K.; Sekizaki, T. Development of an appropriate PCR system for the reclassification of Streptococcus suis. J. Microbiol. Methods 2014, 107, 66–70. [Google Scholar] [CrossRef]
  27. Klindworth, A.; Pruesse, E.; Schweer, T.; Peplies, J.; Quast, C.; Horn, M.; Glöckner, F.O. Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. 2013, 41, e1. [Google Scholar] [CrossRef]
  28. Bolyen, E.; Rideout, J.R.; Dillon, M.R.; Bokulich, N.A.; Abnet, C.C.; Al-Ghalith, G.A.; Alexander, H.; Alm, E.J.; Arumugam, M.; Asnicar, F.; et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 2019, 37, 852–857. [Google Scholar] [CrossRef]
  29. Callahan, B.J.; McMurdie, P.J.; Rosen, M.J.; Han, A.W.; Johnson, A.J.A.; Holmes, S.P. DADA2: High-resolution rample inference from Illumina amplicon data. Nat. Methods 2016, 13, 581–583. [Google Scholar] [CrossRef] [Green Version]
  30. McDonald, D.; Price, M.N.; Goodrich, J.; Nawrocki, E.P.; DeSantis, T.Z.; Probst, A.; Andersen, G.L.; Knight, R.; Hugenholtz, P. An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J. 2012, 6, 610–618. [Google Scholar] [CrossRef]
  31. Diaz, S.; Escobar, J.S.; William-Avila, F. Identification and removal of potential contaminants in 16S rRNA gene sequence datasets from low microbial biomass samples: An example from the mosquito gut. Res. Sq. 2020. [Google Scholar] [CrossRef]
  32. Shannon, C.; Weaver, W. The Mathematical Theory of Communication. Bell Syst. Tech. J. 1948, 27, 379–423. [Google Scholar] [CrossRef] [Green Version]
  33. Eren, M.I.; Chao, A.; Hwang, W.-H.; Colwell, R.K. Estimating the Richness of a Population When the Maximum Number of Classes Is Fixed: A Nonparametric Solution to an Archaeological Problem. PLoS ONE 2012, 7, e34179. [Google Scholar] [CrossRef] [Green Version]
  34. Jaccard, P. Nouvelles Recherches sur la Distribution Florale; Rouge: Lausanne, Switzerland, 1908. [Google Scholar]
  35. Sørensen, T.J. A Method of Establishing Groups of Equal Amplitude in Plant Sociology Based on Similarity of Species Content and Its Application to Analyses of the Vegetation on Danish Commons; I kommission hos E. Munksgaard: København, Demark, 1948. [Google Scholar]
  36. Lozupone, C.; Knight, R. UniFrac: A new phylogenetic method for comparing microbial communities. Appl. Environ. Microbiol. 2005, 71, 8228–8235. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Tang, Z.-Z.; Chen, G.; Alekseyenko, A.V. PERMANOVA-S: Association test for microbial community composition that accommodates confounders and multiple distances. Bioinformatics 2016, 32, 2618–2625. [Google Scholar] [CrossRef]
  38. Anderson, M.J. A new Method for non-parametric multivariate analysis of variance. Austral. Ecol. 2001, 26, 32–46. [Google Scholar] [CrossRef]
  39. Jari Oksanen, F.; Blanchet, G.; Friendly, M.; Kindt, R.; Legendre, P.; McGlinn, D.; Peter, R.; Minchin, R.B.; O’Hara, G.; Simpson, L.; et al. Stevens, Eduard Szoecs and Helene Wagner. Vegan: Community Ecology Package. R Package Version 2.5-7. 2020. Available online: https://CRAN.R-project.org/package=vegan (accessed on 15 February 2021).
  40. Bokulich, N.A.; Kaehler, B.D.; Rideout, J.R.; Dillon, M.; Bolyen, E.; Knight, R.; Huttley, G.A.; Gregory Caporaso, J. Optimizing taxonomic classification of marker-gene amplicon sequences with QIIME 2’s q2-feature-classifier plugin. Microbiome 2018, 6, 90. [Google Scholar] [CrossRef]
  41. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic Local Alignment Search Tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  42. O’Leary, N.A.; Wright, M.W.; Brister, J.R.; Ciufo, S.; Haddad, D.; McVeigh, R.; Rajput, B.; Robbertse, B.; Smith-White, B.; Ako-Adjei, D.; et al. Reference Sequence (RefSeq) Database at NCBI: Current status, taxonomic expansion, and functional annotation. Nucleic Acids Res. 2016, 44, D733–D745. [Google Scholar] [CrossRef] [Green Version]
  43. Mandal, S.; Van Treuren, W.; White, R.A.; Eggesbø, M.; Knight, R.; Peddada, S.D. Analysis of composition of microbiomes: A novel method for studying microbial composition. Microb. Ecol. Health Dis. 2015, 26, 27663. [Google Scholar] [CrossRef] [Green Version]
  44. McMurdie, P.J.; Holmes, S. Phyloseq: An R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE 2013, 8, e61217. [Google Scholar] [CrossRef] [Green Version]
  45. Wickham, H.; Averick, M.; Bryan, J.; Chang, W.; McGowan, L.D.; François, R.; Grolemund, G.; Hayes, A.; Henry, L.; Hester, J.; et al. Welcome to the Tidyverse. J. Open Source Softw. 2019, 4, 1686. [Google Scholar] [CrossRef]
  46. RStudio|Open Source & Professional Software for Data Science Teams. Available online: https://rstudio.com/ (accessed on 15 February 2021).
  47. Wickham, H. Ggplot2: Elegant Graphics for Data Analysis; Springer: New York, NY, USA, 2009; ISBN 978-0-387-98141-3. [Google Scholar]
  48. Gu, Z.; Eils, R.; Schlesner, M. Complex Heatmaps Reveal Patterns and Correlations in Multidimensional Genomic Data. Bioinformatics 2016, 32, 2847–2849. [Google Scholar] [CrossRef] [Green Version]
  49. Gao, C.-H. Gaospecial/GgVennDiagram. Available online: https://github.com/gaospecial/ggVennDiagram (accessed on 15 February 2021).
Figure 1. Boxplots representing median and interquartile ranges of alpha diversity estimated measuring the observed features (A) or Shannon’s index (B) in nasal microbiota from piglets at weaning from healthy control farms (HC) and farms with polyserositis in the nursery caused by Mycoplasma hyorhinis (MH) or caused by Glaesserella parasuis (Glässer’s disease, GD). Outlier is indicated with a grey circle on the plot. Error bars are standard deviation. * means p < 0.05; ** means p ≤ 0.01.
Figure 1. Boxplots representing median and interquartile ranges of alpha diversity estimated measuring the observed features (A) or Shannon’s index (B) in nasal microbiota from piglets at weaning from healthy control farms (HC) and farms with polyserositis in the nursery caused by Mycoplasma hyorhinis (MH) or caused by Glaesserella parasuis (Glässer’s disease, GD). Outlier is indicated with a grey circle on the plot. Error bars are standard deviation. * means p < 0.05; ** means p ≤ 0.01.
Pathogens 10 00603 g001
Figure 2. Beta diversity analysis on Bray Curtis (A) and unweighted Unifrac distances (B) of the nasal microbiota of weaning piglets. Each circle represents the microbial composition of each sample and the size of the circles is proportional to the number of Observed Features as indicated in the legend. Healthy control farms (HC) are depicted in blue, while farms with polyserositis caused by M. hyorhinis (MH) are in orange and G. parasuis (GD) in green palettes. Ellipses are calculated with the Euclidean distances of the grouped samples.
Figure 2. Beta diversity analysis on Bray Curtis (A) and unweighted Unifrac distances (B) of the nasal microbiota of weaning piglets. Each circle represents the microbial composition of each sample and the size of the circles is proportional to the number of Observed Features as indicated in the legend. Healthy control farms (HC) are depicted in blue, while farms with polyserositis caused by M. hyorhinis (MH) are in orange and G. parasuis (GD) in green palettes. Ellipses are calculated with the Euclidean distances of the grouped samples.
Pathogens 10 00603 g002
Figure 3. Relative abundance in the nasal microbiota of weaning piglets of the top 30 most abundant ASVs among all the groups in log10 scale. The X axis (bottom) shows the distribution of each sample by farm and the Y axis represents the top 30 most abundant ASVs grouped by the taxa assignment. Top labels correspond to the study groups: Glässer’s disease (GD), healthy control (HC) and M. hyorhinis (MH) farms. * Asterisks mark the ASVs detected as statistically different in the differential abundance analysis (ANCOM) among the three study groups (GD, MH, HC).
Figure 3. Relative abundance in the nasal microbiota of weaning piglets of the top 30 most abundant ASVs among all the groups in log10 scale. The X axis (bottom) shows the distribution of each sample by farm and the Y axis represents the top 30 most abundant ASVs grouped by the taxa assignment. Top labels correspond to the study groups: Glässer’s disease (GD), healthy control (HC) and M. hyorhinis (MH) farms. * Asterisks mark the ASVs detected as statistically different in the differential abundance analysis (ANCOM) among the three study groups (GD, MH, HC).
Pathogens 10 00603 g003
Figure 4. Relative abundance of M. hyorhinis (A) and G. parasuis (B) amplicon sequence variants (ASVs) in the nasal microbiota of weaning piglets from farms with Glässer’s disease (GD), healthy controls (HC) or farms with polyserositis caused by M. hyorhinis (MH). The phylogenetic relationship of ASVs from each pathogen is represented in Maximum likelihood trees. Branches from ASVs found exclusively in GD farms are colored in green while the ones found only the MH group are colored in red. Venn diagrams representing the distribution of M. hyorhinis (C) and G. parasuis (D) ASVs in the different groups were plotted following the same color pattern.
Figure 4. Relative abundance of M. hyorhinis (A) and G. parasuis (B) amplicon sequence variants (ASVs) in the nasal microbiota of weaning piglets from farms with Glässer’s disease (GD), healthy controls (HC) or farms with polyserositis caused by M. hyorhinis (MH). The phylogenetic relationship of ASVs from each pathogen is represented in Maximum likelihood trees. Branches from ASVs found exclusively in GD farms are colored in green while the ones found only the MH group are colored in red. Venn diagrams representing the distribution of M. hyorhinis (C) and G. parasuis (D) ASVs in the different groups were plotted following the same color pattern.
Pathogens 10 00603 g004
Table 1. Main characteristics of the farms and number of samples included in the study.
Table 1. Main characteristics of the farms and number of samples included in the study.
Farm (n *)GroupHealth StatusProduction
System
Size **Treatments ***Reference
RM (10)MHPolyserositis by M. hyorhinisMulti-site650AmxThis study
GE (10)MHPolyserositis by M. hyorhinisMulti-site800AmxThis study
VL (5)HCHealthyFarrow to finish700Tlt-Ceft[6]
PT (5)HCHealthyMulti-site1000NA[6]
GM (10)HCHealthyMulti-site1200Amx[6]
MT (8)GDGlässer’s diseaseMulti-site3300Pen-Strep[6]
MC (10)GDGlässer’s diseaseFarrow to finish480Ceft[6]
RC (6)GDGlässer’s diseaseMulti-site1400Ceft[6]
EJ (10)GDGlässer’s diseaseMulti-site2000Enro[6]
* number of samples included from each farm; ** n. of sows; *** Perinatal antimicrobials: Amx, amoxicillin; Tlt, tulathromycin; Ceft, ceftiofur; Pen, Penicillin; Strep, streptomycin; Enro, enrofloxacin; NA, not available.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Blanco-Fuertes, M.; Correa-Fiz, F.; Fraile, L.; Sibila, M.; Aragon, V. Altered Nasal Microbiota Composition Associated with Development of Polyserositis by Mycoplasma hyorhinis. Pathogens 2021, 10, 603. https://0-doi-org.brum.beds.ac.uk/10.3390/pathogens10050603

AMA Style

Blanco-Fuertes M, Correa-Fiz F, Fraile L, Sibila M, Aragon V. Altered Nasal Microbiota Composition Associated with Development of Polyserositis by Mycoplasma hyorhinis. Pathogens. 2021; 10(5):603. https://0-doi-org.brum.beds.ac.uk/10.3390/pathogens10050603

Chicago/Turabian Style

Blanco-Fuertes, Miguel, Florencia Correa-Fiz, Lorenzo Fraile, Marina Sibila, and Virginia Aragon. 2021. "Altered Nasal Microbiota Composition Associated with Development of Polyserositis by Mycoplasma hyorhinis" Pathogens 10, no. 5: 603. https://0-doi-org.brum.beds.ac.uk/10.3390/pathogens10050603

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop