Glaesserella parasuis is the etiological agent of Glässer’s disease, a common pathology in the pork industry with higher prevalence in the postweaning period. Vaccination is one of the strategies to control this disease. Here, we investigated the effect that sow vaccination against virulent strains of G. parasuis had in the nasal microbiota of their offspring. Nasal swabs from fifteen days-old piglets from vaccinated (vs-P, n = 11) and unvaccinated sows (cs-P, n = 11) were obtained and DNA was extracted for 16S amplicon sequencing. Microbiota composition was different, with lower diversity in vs-P, and a strong clustering of the groups in beta diversity analysis. Among the 1509 sequences associated to either study group, all the sequences classified as G. parasuis (10 ASVs) had lower relative abundance in the vs-P group. A list of 32 inferred metabolic pathways were statistically different between groups. A distinctive structure of the two microbial networks was detected, with modules in the cs-P not conserved in the vs-P network. In conclusion, vaccination of the sows had a large effect in the microbiota composition of their offspring that went beyond the effect on the targeted pathogen. The mechanisms underneath these changes may include alteration of the microbiota network due to the elimination of the targeted pathogen and/or immunological changes.
Glaesserella parasuis is the etiological agent of Glässer’s disease, an infectious disease with an important economic impact on swine production1. This microorganism is an early colonizer of the upper respiratory tract and is part of the healthy nasal microbiota during all the stages of pig’s life. In order to reduce the prevalence of Glässer’s disease, different strategies have been used. Antimicrobial treatments and vaccines are the main tools to control disease2, although the global concern about rising antimicrobial resistances makes vaccination the method of choice. The natural source of this microorganism for piglets resides in the dams and colonization is known to occur very early in life3. Vaccination of the sows increases the specific antibodies transferred to their offspring through colostrum intake4 and may delay the colonization by G. parasuis3. Thus, sow vaccination can be a good strategy to decrease the prevalence of the virulent strains of this microorganism in the offspring, when using a vaccine that specifically targets these strains5.
The microbiota, known as the set of microorganisms that inhabit a host6,7, has different essential functions for the host such as nutrition, preservation of the mucosal homeostasis, prevention of pathogen colonization, and maturation of the immune system8. Nasal microbiota, as gut microbiota, is affected by a wide range of factors in piglets, such as age, contact with the mother or environment9,10,11. Several studies have examined the impact of the pig microbiome on the immune response to vaccination12,13, but the impact of vaccination on the microbiota had not been extensively analyzed, with a few studies focusing on gut microbiota14,15,16. No information is yet available about the effect that specific stimulation of the immune system by vaccines could have on the nasal microbiota of the piglet.
Here, we report the effect of vaccination of sows against virulent G. parasuis on the nasal microbiota composition of their offspring. A reduction in the relative abundance of G. parasuis sequences was found in the samples from the piglets derived from vaccinated sows. In addition, other differences in the microbial composition were detected; especially, on the most relatively abundant bacterial taxa, as well as in inferred pathways, which were divergent between piglets born to vaccinated or non-vaccinated sows. Correlation Network Analysis (CNA) resulted in two complex networks with a different number of nodes and edges. The structure of the modules containing the most connected nodes in the co-occurrent ASVs control network was not conserved in the network from vaccinated sows.
Sequencing throughput stats
After the denoising and quality control steps, a total of 9817 features were found in the 22 samples of the study. A total of 9,202,583 sequences were analyzed with a minimum ASVs frequency per sample of 211,733 and maximum of 649,315, being the mean 418,299. The ASVs had a mean length of 414 bp (range 236–444), where 98% of the sequences were 428 bp length. After processing, 73% of the raw sequences passed the quality-control steps and were used in the downstream analyses.
Sow vaccination decreased diversity in the nasal microbiota of the offspring
Richness and alpha diversity metrics (Observed Features, Simpson, Shannon, Faith phylogenetic diversity and Chao indexes) were calculated and compared between piglets born to vaccinated sows (vs-P) and to non-vaccinated sows (cs-P). Nasal microbiota from vs-P was significantly less diverse than cs-P when using Simpson diversity index and Simpson evenness index (Fig. 1). The rest of indexes showed the same tendency although they did not detect statistically significant differences between both groups.
The composition of the nasal microbiota evaluated through several beta diversity metrics showed a divergency between the vs-P and cs-P groups and this divergency was statistically significant when evaluated either quantitatively or qualitatively (P < 0.001, PERMANOVA test). The first axis in the PCoA (PC1) using Bray Curtis dissimilarity index, captured the highest amount of variation of the input data (53.14%; Fig. 2A) and a high percentage of the differences were explained by the vaccination of the sows (R2 = 51.47%; Adonis function). Spatial representation of the Jaccard distance (a metric that only calculates the presence or absence of the different ASVs in each sample) showed a 16.08% of variation among the samples in the PC1 (Supplementary Fig. S1), which was mainly explained by sow vaccination (R2 = 15%, Adonis function). Similar results were also obtained with metrics that include phylogenetic relationships between the microorganisms. Thus, weighted (Fig. 2B) and unweighted (Supplementary Fig. S1) Unifrac showed PC1 with percentages of 44.12%, and 15.52%, respectively, which were in agreement with the percentage of explanation due to sow vaccination (R2 = 41.2% for the weighted and R2 = 11.5% for the unweighted distance). In all these analyses, sow vaccination explained a higher percentage of variation when the analysis included quantitative data.
Taxonomical assignment and differential abundance analysis
The denoised sequences corresponded to 9817 ASVs with different taxonomic assignment. Taxa assignation at the different levels showed a percentage of non-assigned taxa: 3.1% at phylum and 5.05% at genus. At phylum level, Firmicutes, Proteobacteria and Actinobacteria were the most relative abundant taxa in both cs-P and vs-P groups, while the ten most relatively abundant genera were Enhydrobacter, SMB53, Rothia, Moraxella, Neisserieceae sp., Lactobacillus, Clostridaceae sp., Catenibacterium, Prevotella and Fusobacterium. The top 30 most relatively abundant ASVs among all the samples from both groups were heterogenously distributed within each group and included 6 ASV classified as Moraxella, 3 as Rothia, 2 as Prevotella, and the rest of taxa represented by a unique ASV (Fig. 3).
To unravel the differences in bacterial communities inhabiting the nasal microbiota from each group, we used two different approaches. First, the analysis of composition of microbiomes (ANCOM) showed 21 ASVs that differed between each study group (Supplementary Table S1). Eighteen of those ASVs were more abundant in the control group, and included Moraxella, Staphylococcus, Prevotella and Actinobacillus; all with relative abundances lower than 0.33%. ANCOM was also performed at different taxonomic levels and detected ten different phyla that were differentially abundant according to the study group. Seven out of these ten were increased in the vs-P group: Tenericutes, SR1, Planctomycetes, Euryarcheota, Lentisphareae, TM7, Fibrobacteres, Fusobacteres, Verrucomicrobia, Spirochaetes, and Actinobacteria. On the other hand, the phyla that were more abundant in the cs-P group were GN02, Firmicutes, and Bacteroidetes. At genus level, Selenomonas, Megasphaera, Akkermansia, and Anaeroplasma were associated with the vs-P group.
Second, the differential ranking Songbird method was used as a complementary method to infer taxa associated to either cs-P or vs-P group, showing a list of 1509 ASVs. Once the low relatively abundant ASVs (less than 0.1% in either group) were removed, 186 ASVs showed a statistically significant differential abundance, with 83 of them associated to the cs-P group and 103 ASVs associated to the vs-P group. The differential abundant ASVs included the 30 most abundant ASVs across all the groups, whose association to either vs-P or cs-P is depicted in Supplementary Figure S2. The 30 ASVs most associated to the vs-P group included ASVs classified as Ruminococcaceae (4 ASVs), Megasphaera (4 ASVs), Paraprevotellaceae (3 ASVs) and Lachnospiraceae (2 ASVs). In contrast, among the 30 most associated to the cs-P group, we found ASVs assigned to Moraxella (5 ASVs), Moraxellaceae (3 ASVs), Ruminococcaceae (3 ASVs), Streptococcus (3 ASVs), Staphylococcus genus (1 ASV), Flavobacteriaceae (2 ASVs) and Rothia nasimurium (2 ASVs) (Fig. 4). In general, the differential ASVs classified within a given taxon seemed to be predominantly associated to one group or the other, e.g., Moraxella ASVs were associated with cs-P, while Megasphaera ASVs were associated with vs-P (Fig. 4). Among the most relatively abundant ASVs (Supplementary Figure S2), twelve were shared with the most associated to each group (with high log-fold change value) corresponding to the following taxa: Rothia nassimurium (2 ASVs), Prevotella copri, Flavobacteriaceae, Helcococcus, Lachnospiraceae, Megasphaera, Fusobacterium, Moraxellaceae and Moraxella (3 ASVs).
Since the vaccine used targeted virulent strains of G. parasuis, we focused our analysis on the potential effect of this vaccine on this species. We found a total of ten different ASVs classified as G. parasuis, which in global showed lower relative abundance in the vs-P group than in the cs-P group (G. parasuis relative abundance: 0.058% ± 0.029 for vs-P and 0.225% ± 0.112 for cs-P). From these 10 G. parasuis ASVs, 5 were statistically associated to the control group, where three ASVs were exclusively found in the cs-P group. These findings were found by the Songbird differential rank analysis, but they were not detected using ANCOM.
Correlation network analysis
Correlation Network Analysis (CNA) at ASV level, resulted in two different undirected graphs when samples from the different groups were analyzed separately. The topology of the two graphs showed different characteristics. The cs-P network had 3,552 nodes with 1,678,545 edges and a network density of 0.266, while the vs-P network had 3,048 nodes with 1,324,390 edges and a network density of 0.285. On the other hand, the diameter of both networks was three.
The top ten most-connected nodes in the cs-P network ranged between 1,640 and 1,670 edges, and they were grouped in three different modules. Module 6 included ASVs from Actinobacillus, Ruminococcus flavefaciens, Erysipelotrichaceae, Staphylococcus, Glaesserella parasuis and Filifactor out of 76 total component nodes; module 8 (67 ASVs) included WCHB1-25 and Corynebacterium, while module 13 (37 ASVs) included Oscillospira, and Ruminococcus gnavus. On the other hand, the most connected nodes in the vs-P were grouped in five different modules. Module 1 (124 ASVs) included the following most-connected taxa: Lachnospiraceae and Clostridiales; whereas Coprococcus, Bacteroidales and Sphaerochaeta taxa were highly-connected members of module 2 (91 ASVs). ASVs from Moraxellaceae, Neisseriaceae and Moraxella belonged to module 12 (37 ASVs), while the last two modules included Oscillospira in module 32 (23 ASVs) and Lachnospiraceae in module 41 (20 ASVs). None of the 10 most-connected nodes in the two networks was shared.
Among all the modules containing the ten most-connected nodes, module 6 from cs-P network had the highest number of these highly connected nodes, including an ASV classified as G. parasuis (Fig. 5). When we compared the networks, we found that the ASVs from module 6 in the cs-P network, clustered differently in the vs-P network. From the 76 ASVs composing the module 6 from the cs-P network, 72 were found in 31 different modules in the vs-P network. Noteworthy, among these 31 modules, we found the module 23 containing an ASV classified as G. parasuis, reinforcing that the rearrangements in the community structure were due to vaccination against this pathogen. Five out of these 31 modules contained the top-ten most-connected nodes in the vs-P network (Fig. 5). The remaining 4 nodes, assigned to G. parasuis, Neisseria, Staphylococcus and Actinobacillus porcitonsillarum, were absent in vs-P network.
Metagenome prediction and functional analysis inference
The metagenome was predicted from the 16S partial sequences in the two study groups to infer the potential functionality of the microbiota using the KEGG public database. PCA of the functional pathways showed a robust clustering of the samples by the group they belonged (Fig. 6), in agreement with the composition (beta diversity) analysis (Fig. 2) and suggesting that the microbial composition changes led to functional changes. Only 32 functional pathways were significantly different between the two study groups (Kruskall Wallis test, P value with Bonferroni correction < 0.05, Supplementary Table S2), where 15 of them presented a two-fold change (Fig. 7). The most relatively abundant pathways in the cs-P group were the isopropanol biosynthesis and L-tryptophan degradation to 2-amino-3-carboxymuconate semialdehyde, while the most relatively abundant in the vs-P group were L-arabinose degradation IV and the tetrahydromethanopterin biosynthesis. Pathways related with the biosynthesis with the peptidoglycan were in higher abundance in the cs-P group and with the biosynthesis of the LPS was higher in the vs-P group.
Vaccination against specific pathogens has been an extensive method to control infectious diseases for a long time now. Primary function of all vaccines is the stimulation of the immune system against specific antigens, while one of the main tasks of the microbiota is to contribute to the maturation of the immune system17,18,19,20. Thus, the influence of the microbiota on the vaccine response has been more studied than the effect that the vaccines may have in the host microbiota and, in particular, in the microbiota of the progeny16,21. Here, we have shown that sow vaccination modified the nasal microbiota composition of their piglets with an effect extended beyond the targeted pathogen.
Vaccination reduces the shedding of the targeted pathogens14. Accordingly, here we observed a reduction of the relative abundance of all the G. parasuis ASVs after vaccination. Some of these ASVs classified as G. parasuis were statistically associated with the control group, while two of them were absent in piglets born to vaccinated sows. This reduction may be due not only to the decreased shedding from the mother but also by the maternal immunity transferred through the colostrum to the piglets as specific antibodies22. As the vaccine used was designed against virulent strains of G. parasuis, it would be very interesting to determine if in fact, the virulent strains were specifically reduced in the piglets from vaccinated sows. Additional analyses, including shotgun sequencing metagenomics, may provide better resolution to study the different strains (variants) in the microbiota.
Sow vaccination also affected the general composition of the piglets’ nasal microbiota. The nasal microbiota of the offspring from vaccinated sows showed a slight decrease in alpha-diversity. This was especially observed with the Simpson index (biased towards abundant features23,24), indicating that the bacterial communities that are in relatively higher abundance in the nasal cavity of piglets drove the differences observed in our groups. In agreement, beta diversity analysis showed a higher percentage of explanation of the clustering in the indexes that include the quantity of the ASVs. Taken together, our findings indicate that vaccination of the mother did not change the main components of the nasal microbiota in the offspring but caused a clear disparity in the distribution of the taxa. In the differential abundance analysis, we observed that Songbird, a differential rank approach, discovered many ASVs that were not detected by ANCOM (1509 ASVs detected by Songbird vs 21 ASVs by ANCOM). It was reported that ANCOM fails to infer differential microbes in some composition scenarios, where these are actually changing25. The top ranked ASVs (most positive or negatively correlated to sow vaccination) within a given taxa tended to be associated to only one of the groups. Thus, ASVs from Moraxella, Streptococcus, Staphylococcus and Rothia were associated with cs-P, while ASVs classified as Megasphaera, YRF22 and RF16 were associated with the vs-P group. However, this specific association of the ASVs from one taxon to one group was not observed in other ASVs, reinforcing the importance of detailed analyses at the ASV level. Strains within a given taxon can be heterogeneous, as demonstrated for Glaesserella26,27 or Moraxella spp. nasal isolates28, and therefore they can play a different role in the microbiota.
Correlation Network Analysis has been a useful tool to identify global relationships between the microorganisms that compose the microbiota. Here, we observed that the network of piglets from vaccinated sows had lower complexity, demonstrated by the lower number of nodes, modules and edges, when compared to the control network. The structure of the modules containing the most connected nodes from the control network were not conserved in the structure shaped by the vaccination of the sows. Both the lower shedding of the pathogen by the mother, which in turn may produce a reduction of the pathogen in the piglet, and the effect of the maternal immunity transferred by the colostrum, might contribute to this rearrangement of the network topology.
The strong clustering between the two groups of piglets observed in the beta diversity analyses, resulted in a clear differential predicted metagenome, i.e., differential predicted functionality. Some of the metabolic functional pathways that were in higher abundance in the control group, such as isopropanol and L-tryptophan degradation are related to an immune suppressive status in mammals, as reported in the literature. On the other hand, none of the pathways in the vs-P group was related to any function of the immune system. The changes that we observed in the functionality may affect the health status or the immune response in the future, but this is difficult to assess and deserves further exploration. In general, the ASVs assigned to taxa of potential pathogens (Staphylococcus, Neisseria, Streptococcus and Mycoplasmataceae) were negatively correlated with the vs-P group, suggesting that the reduction of virulent strains of G. parasuis may contributed to drag out other associated pathogens in the nasal cavity26,29,30. Moreover, potentially beneficial bacteria, such as Prevotella, were associated with the vs-P group. Members of the Prevotella genus have been associated to positive outcomes in pig production, including growth performance and immune response31,32. Indeed, 122 from the 172 ASVs classified as Prevotella were associated with sow vaccination.
In summary, sow vaccination with a peptide directed against virulent G. parasuis influenced the nasal microbiota of the offspring, not only reducing the relative abundance of the targeted pathogen but also modifying the balances in general composition of the microbiota during at least 15 days of life. These changes in nasal microbiota composition are expected to affect its functionality.
Materials and methods
Study design and sampling
The experimental study was described previously by López-Serrano et al.5 and was performed under the approval of the Ethics Commission in Animal Experimentation of the Generalitat de Catalunya (Protocol number 9211) and was carried out in compliance with the ARRIVE guidelines. Briefly, four pregnant sows were included in the study based on low antibody levels against G. parasuis. Two sows were vaccinated intramuscularly with 2.5 ml of 100 μg/mL F4 (an antigen that is specific of virulent strains of G. parasuis) in combination with Carbopol 5984 EP polymer adjuvant (Lubrizol, Cleveland, OH, USA), while the other two sows remained unvaccinated and were used as control. Vaccination was performed at 32 and 12 days pre-farrowing. Newborn piglets were kept with their biological mother during their first hours of life (at least 12 h). Twenty-four hours after delivery, only one sow from each group was maintained in the study, with biological piglets and cross-fostered piglets from the other sow from the same study group. Nasal swabs were taken from 11 piglets from vaccinated sows (vs-P) and 11 piglets from control sows (cs-P) at 15 days of life. Nasal swabs were resuspended in 500 μl de PBS, vortexed for 30 s and stored at -80ºC until used for microbiota analysis.
DNA extraction and 16S rRNA gene amplicon sequencing
DNA was extracted from 200 μl of resuspended swabs using the Nucleospin Blood (Macherey–Nagel) kit. Quantity and quality assessment of the DNA was performed using BioDrop DUO (BioDrop Ltd). DNA samples were submitted to the Servei de Genòmica (Universitat Autònoma de Barcelona) for 16S rRNA amplicon sequencing using Illumina MiSeq pair-end 2 × 250 bp technology following the manufacturer instructions (MS-102–2003 MiSeq Reagent Kit v2, 500 cycle). The region amplified covered the V3-V4 hypervariable regions33. The complete dataset from this study is available at the NCBI database, under SRA accession number PRJNA783272.
Bioinformatics analysis was done using the Quantitative insights into microbial ecology (Qiime2) platform34. Quality control step was done using q2-dada2 plugin, where all the reads with a Phred quality-value under 20 and a length under of 234 bp in the forward read and 229 bp in the reverse read, were removed. Chimeric-reads check was also done and were removed for further analysis. In addition to the quality-control step, an alignment against the reference Greengenes (version 13.8) 16S rRNA gene database35 was done to remove sequences not matching with an identity of 80% and a query length of 50%. Richness and alpha diversity analysis were done using Observed Features (in this case, Amplicon Sequence Variants, ASVs)36, Shannon index36, Chao index37, Simpson Index and Simpson evenness38. Different beta diversity metrics were used to assess the diversity across the samples, both qualitatively with Jaccard similarity coefficient39 and quantitatively using Bray Curtis dissimilarity index40. Moreover, phylogenetic metrics weighted and unweighted Unifrac41 were calculated as well. These distance matrixes were used to perform Principal Coordinate Analysis (PCoA) using core-metrics plugin and a PERMANOVA test was done to analyze the clustering of the sample groups. To extract the percentage of variations explained by each metadata column, the Adonis function was performed on every distance matrix using Vegan42 package. Taxonomic assignation of each amplicon sequence variant was done using the Qiime2 classifier trained with the V3-V4 region from 16S gene and the Greengenes database (13.8 version)35. Differential abundant analysis was done using the analysis of composition of microbiomes (ANCOM)43 and the Songbird differential ranking25 algorithms.
Metagenome prediction from the 16S rRNA sequencing was done using PICRUSt244. This software aligns the ASVs and place them in a tree (EPA-NG, GAPPA) with reference sequences and it infers the family genes and function pathways from KEGG45 orthologs and Enzyme Classification numbers. Together with the gene family copy number references, it predicts the gene content per ASV and with the abundance of each ASV in the study dataset it can determine the gene family abundance per sample. Mapping the gene family to pathways abundances was done using MinPath46 integrated in the PICRUSt2 pipeline. Statistical Analysis of Microbial Profiles (STAMP)47 software was used to visualize and perform statistical analysis over the predicted functional pathways outputted by PICRUSt2. Statistical analysis was based on a Kruskal Wallis test and Tukey–Kramer post-hoc test of the functional pathways of the different study groups were done using a Kruskal Wallis48 test and a Tukey–Kramer post-hoc test49. After a correction with Bonferroni P values under 0.05 were considered significant.
We built two separate networks from the ASVs of each group samples (cs-P and vs-P) Correlation Network Analysis (CNA) was performed using the Sparse Co-occurrence Network Investigation for Compositional data (SCNIC) integrated in the q2-scnic plugin from Qiime250 software toolkit. All the samples that had less than 500 features and also all the features that had an average mean less than 2 across all the samples, were filtered out. Pairwise correlations between all our ASVs in our dataset were calculated with sparCC51 algorithm. From that correlation table, we built our network with a minimum of R correlation-value cutoff of 0.35. Detection of modules was done by the summarization and clustering high connected nodes. Output files with nodes membership and edged connections were analyzed in Cytoscape52 software and extracted to tables to further data analysis to query ASVs and nodes. These tables were processed in Rstudio53 and merged with previous results using tidyverse package54.
All the plots and figures were done using R script language. Rstudio53 was used as a coding work environment. The packages used to build the plots and figures were specifically: qiime2R55, ggplot256 and tidyverse54.
The animal study was performed in the BSL3 animal facilities of IRTA-CReSA (Bellaterra, Spain) following proper veterinary practices, in accordance with European (Directive 2010/63/EU) and Spanish (Real Decreto 53/2013) regulation and with the approval of the Ethics Commission in Animal Experimentation of the Generalitat de Catalunya (Protocol number 9211).
The complete dataset from this study is available at the NCBI database, under SRA accession number PRJNA783272 available at https://www.ncbi.nlm.nih.gov/sra/PRJNA783272.
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. 242, 108595 (2020).
Zhao, Z. et al. Analysis of efficacy obtained with a trivalent inactivated Haemophilus parasuis serovars 4, 5, and 12 vaccine and commercial vaccines against Glässer’s disease in piglets. Can. J. Vet. Res. 81, 22–27 (2017).
Cerdà-Cuéllar, M. et al. Sow vaccination modulates the colonization of piglets by Haemophilus parasuis. Vet. Microbiol. 145, 315–320 (2010).
Devillers, N., Farmer, C., Le Dividich, J. & Prunier, A. Variability of colostrum yield and colostrum intake in pigs. Anim. Int. J. Anim. Biosci. 1, 1033–1041 (2007).
López-Serrano, S. et al. Sow vaccination with a protein fragment against virulent Glaesserella (Haemophilus) parasuis modulates immunity traits in their offspring. Vaccines 9, 534 (2021).
NIH HMP Working Group et al. The NIH Human Microbiome Project. Genome Res. 19, 2317–2323 (2009).
Robinson, C. M. & Pfeiffer, J. K. Viruses and the microbiota. Annu. Rev. Virol. 1, 55–69 (2014).
Brestoff, J. R. & Artis, D. Commensal bacteria at the interface of host metabolism and the immune system. Nat. Immunol. 14, 676–684 (2013).
Slifierz, M. J., Friendship, R. M. & Weese, J. S. Longitudinal study of the early-life fecal and nasal microbiotas of the domestic pig. BMC Microbiol. 15, 184 (2015).
Obregon-Gutierrez, P., Aragon, V. & Correa-Fiz, F. Sow contact is a major driver in the development of the nasal microbiota of piglets. Pathogens 10, 697 (2021).
Valeris-Chacin, R. et al. Relationships among fecal, air, oral, and tracheal microbial communities in pigs in a respiratory infection disease model. Microorganisms 9, 252 (2021).
de Jong, S. E., Olin, A. & Pulendran, B. The impact of the microbiome on immunity to vaccination in humans. Cell Host Microbe 28, 169–179 (2020).
Ruck, C. E., Odumade, O. A. & Smolen, K. K. Vaccine interactions with the infant microbiome: Do they define health and disease?. Front. Pediatr. 8, 811 (2020).
Leite, F. L. L., Singer, R. S., Ward, T., Gebhart, C. J. & Isaacson, R. E. Vaccination against lawsonia intracellularis decreases shedding of Salmonella enterica serovar typhimurium in co-infected pigs and alters the gut microbiome. Sci. Rep. 8, 2857 (2018).
Kim, A. H., Hogarty, M. P., Harris, V. C. & Baldridge, M. T. The complex interactions between rotavirus and the gut microbiota. Front. Cell. Infect. Microbiol. 10, 586–751 (2020).
Guevarra, R. B. et al. Oral vaccination against Lawsonia intracellularis changes the intestinal microbiome in weaned piglets. Animals 11, 2082 (2021).
Ciabattini, A., Olivieri, R., Lazzeri, E. & Medaglini, D. Role of the microbiota in the modulation of vaccine immune responses. Front. Microbiol. 10, 1305 (2019).
Munyaka, P. M. et al. Discovery of predictors of Mycoplasma hyopneumoniae vaccine response efficiency in pigs: 16S rRNA gene fecal microbiota analysis. Microorganisms 8, 1151 (2020).
Lynn, D. J., Benson, S. C., Lynn, M. A. & Pulendran, B. Modulation of immune responses to vaccination by the microbiota: Implications and potential mechanisms. Nat. Rev. Immunol. https://doi.org/10.1038/s41577-021-00554-7 (2021).
Borey, M. et al. Links between fecal microbiota and the response to vaccination against influenza A virus in pigs. NPJ Vaccines 6, 92 (2021).
Salgado, V. R. et al. Effects of 10-valent pneumococcal conjugate (PCV10) vaccination on the nasopharyngeal microbiome. Vaccine 38, 1436–1443 (2020).
Le Doare, K., Holder, B., Bassett, A. & Pannaraj, P. S. Mother’s milk: A purposeful contribution to the development of the infant microbiota and immunity. Front. Immunol. 9, 361 (2018).
Kim, B.-R. et al. Deciphering diversity indices for a better understanding of microbial communities. J. Microbiol. Biotechnol. 27, 2089–2093 (2017).
Daly, A. J., Baetens, J. M. & De Baets, B. Ecological diversity: Measuring the unmeasurable. Mathematics 6, 119 (2018).
Morton, J. T. et al. Establishing microbial composition measurement standards with reference frames. Nat. Commun. 10, 2719 (2019).
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. 51, 7 (2020).
Olvera, A., Calsamiglia, M. & Aragon, V. Genotypic diversity of Haemophilus parasuis field strains. Appl. Environ. Microbiol. https://doi.org/10.1128/AEM.02834-05 (2006).
López-Serrano, S., Galofré-Milà, N., Costa-Hurtado, M., Pérez-de-Rozas, A. M. & Aragon, V. Heterogeneity of Moraxella isolates found in the nasal cavities of piglets. BMC Vet. Res. 16, 28 (2020).
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 10, 603 (2021).
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 Genomics 17, 1–14 (2016).
Amat, S., Lantz, H., Munyaka, P. M. & Willing, B. P. Prevotella in pigs: The positive and negative associations with production and health. Microorganisms 8, (2020).
Correa-Fiz, F., Gonçalves dos Santos, J. M., Illas, F. & Aragon, V. Antimicrobial removal on piglets promotes health and higher bacterial diversity in the nasal microbiota. Sci. Rep. 9, 6545 (2019).
Klindworth, A. et al. Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. 41, e1 (2013).
Bolyen, E. et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 37, 852–857 (2019).
McDonald, D. et al. An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J. 6, 610–618 (2012).
Shannon, C. E. & Weaver, W. The Mathematical Theory of Communication (University of Illinois Press, 1963).
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 7, e34179 (2012).
Grabchak, M., Marcon, E., Lang, G. & Zhang, Z. The generalized Simpson’s entropy is a measure of biodiversity. PLoS ONE 12, e0173305 (2017).
Jaccard, P. Nouvelles recherches sur la distribution florale. (Rouge, 1908).
Bray, J. R. & Curtis, J. T. An ordination of upland forest communities of southern Wisconsin. Ecol Monogr. 27, 325–349 (1957).
Lozupone, C. & Knight, R. UniFrac: A new phylogenetic method for comparing microbial communities. Appl. Environ. Microbiol. 71, 8228–8235 (2005).
Dixon, P. VEGAN, a package of R functions for community ecology. J. Veg. Sci. 14, 927–930 (2003).
Mandal, S. et al. Analysis of composition of microbiomes: A novel method for studying microbial composition. Microb. Ecol. Health Dis. 26, 27663 (2015).
Douglas, G. M. et al. PICRUSt2 for prediction of metagenome functions. Nat. Biotechnol. 38, 685–688 (2020).
Kanehisa, M. & Goto, S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 28, 27–30 (2000).
Ye, Y. & Doak, T. G. A parsimony approach to biological pathway reconstruction/inference for genomes and metagenomes. PLOS Comput. Biol. 5, e1000465 (2009).
Parks, D. H., Tyson, G. W., Hugenholtz, P. & Beiko, R. G. STAMP: statistical analysis of taxonomic and functional profiles. Bioinformatics 30, 3123–3124 (2014).
Kruskal, W. H. & Wallis, W. A. Use of ranks in one-criterion variance analysis. J. Am. Stat. Assoc. 47, 583–621 (1952).
Tukey, J. W. Exploratory Data Analysis (Addison-Wesley, 1977).
Shaffer, M., Thurimella, K. & Lozupone, C. A. SCNIC: Sparse Correlation Network Investigation for Compositional Data. https://doi.org/10.1101/2020.11.13.380733v1(2020).
Friedman, J. & Alm, E. J. Inferring correlation networks from genomic survey data. PLOS Comput. Biol. 8, e1002687 (2012).
Shannon, P. et al. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498–2504 (2003).
RStudio: Integrated Development for R. RStudio, PBC, Boston, MA. (2020). http://www.rstudio.com/.
Wickham, H. et al. Welcome to the tidyverse. J. Open Source Softw. 4, 1686 (2019).
Bisanz, J. Tutorial: Integrating QIIME2 and R for data visualization and analysis using qiime2R v0.99.6 (2021).
Wickham, H. ggplot2: Elegant Graphics for Data Analysis. (Springer, 2009). https://doi.org/10.1007/978-0-387-98141-3
This study was supported by the Ministry of Economy and Competitiveness (MINECO) from the Spanish Government (Grants AGL2016-77361-R and PID2019-106233RB-I00). M.B-F. is a recipient of a FPI fellowship from the same institution (BES-2017-082224). The authors are also grateful to the Centres de Recerca de Catalunya (CERCA) Programme. The authors are grateful to the personnel from the animal facility of CReSA-IRTA for their excellent animal care.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Blanco-Fuertes, M., Correa-Fiz, F., López-Serrano, S. et al. Sow vaccination against virulent Glaesserella parasuis shapes the nasal microbiota of their offspring. Sci Rep 12, 3357 (2022). https://doi.org/10.1038/s41598-022-07382-2
This article is cited by
BMC Microbiology (2022)