System immunology-based identification of blood transcriptional modules correlating to antibody responses in sheep

Correlates of disease severity in bluetongue as a model of acute arbovirus infection

Most viral diseases display a variable clinical outcome due to differences in virus strain virulence and/or individual host susceptibility to infection. Understanding the biological mechanisms differentiating a viral infection displaying severe clinical manifestations from its milder forms can provide the intellectual framework toward therapies and early prognostic markers. This is especially true in arbovirus infections, where most clinical cases are present as mild febrile illness. Here, we used a naturally occurring vector-borne viral disease of ruminants, bluetongue, as an experimental system to uncover the fundamental mechanisms of virus-host interactions resulting in distinct clinical outcomes. As with most viral diseases, clinical symptoms in bluetongue can vary dramatically. We reproduced experimentally distinct clinical forms of bluetongue infection in sheep using three bluetongue virus (BTV) strains (BTV-1IT2006, BTV-1IT2013 and BTV-8FRA2017). Infected animals displayed clinical signs varying from clinically unapparent, to mild and severe disease. We collected and integrated clinical, haematological, virological, and histopathological data resulting in the analyses of 332 individual parameters from each infected and uninfected control animal. We subsequently used machine learning to select the key viral and host processes associated with disease pathogenesis. We identified and experimentally validated five different fundamental processes affecting the severity of bluetongue: (i) virus load and replication in target organs, (ii) modulation of the host type-I IFN response, (iii) pro-inflammatory responses, (iv) vascular damage, and (v) immunosuppression. Overall, we showed that an agnostic machine learning approach can be used to prioritise the different pathogenetic mechanisms affecting the disease outcome of an arbovirus infection.

Altered transcriptomic immune responses of maintenance hemodialysis patients to the COVID-19 mRNA vaccine

Background

End-stage renal disease (ESRD) patients experience immune compromise characterized by complex alterations of both innate and adaptive immunity, and results in higher susceptibility to infection and lower response to vaccination. This immune compromise, coupled with greater risk of exposure to infectious disease at hemodialysis (HD) centers, underscores the need for examination of the immune response to the COVID-19 mRNA-based vaccines.

Methods

The immune response to the COVID-19 BNT162b2 mRNA vaccine was assessed in 20 HD patients and cohort-matched controls. RNA sequencing of peripheral blood mononuclear cells was performed longitudinally before and after each vaccination dose for a total of six time points per subject. Anti-spike antibody levels were quantified prior to the first vaccination dose (V1D0) and 7 d after the second dose (V2D7) using anti-spike IgG titers and antibody neutralization assays. Anti-spike IgG titers were additionally quantified 6 mo after initial vaccination. Clinical history and lab values in HD patients were obtained to identify predictors of vaccination response.

Results

Transcriptomic analyses demonstrated differing time courses of immune responses, with prolonged myeloid cell activity in HD at 1 wk after the first vaccination dose. HD also demonstrated decreased metabolic activity and decreased antigen presentation compared to controls after the second vaccination dose. Anti-spike IgG titers and neutralizing function were substantially elevated in both controls and HD at V2D7, with a small but significant reduction in titers in HD groups (p<0.05). Anti-spike IgG remained elevated above baseline at 6 mo in both subject groups. Anti-spike IgG titers at V2D7 were highly predictive of 6-month titer levels. Transcriptomic biomarkers after the second vaccination dose and clinical biomarkers including ferritin levels were found to be predictive of antibody development.

Conclusions

Overall, we demonstrate differing time courses of immune responses to the BTN162b2 mRNA COVID-19 vaccination in maintenance HD subjects comparable to healthy controls and identify transcriptomic and clinical predictors of anti-spike IgG titers in HD. Analyzing vaccination as an in vivo perturbation, our results warrant further characterization of the immune dysregulation of ESRD.

Funding

F30HD102093, F30HL151182, T32HL144909, R01HL138628. This research has been funded by the University of Illinois at Chicago Center for Clinical and Translational Science (CCTS) award UL1TR002003.

Immune response to the mRNA COVID-19 vaccine in hemodialysis patients: cohort study

Background

End-stage renal disease (ESRD) patients experience immune compromise characterized by complex alterations of both innate and adaptive immunity, and results in higher susceptibility to infection and lower response to vaccination. This immune compromise, coupled with greater risk of exposure to infectious disease at hemodialysis (HD) centers, underscores the need for examination of the immune response to the COVID-19 mRNA-based vaccines.

Methods

A transcriptomic analysis of the immune response to the Covid-19 BNT162b2 mRNA vaccine was assessed in 20 HD patients and cohort-matched controls. RNA sequencing of peripheral blood mononuclear cells (PBMCs) was performed longitudinally before and after each vaccination dose for a total of six time points per subject. Anti-spike antibody levels were quantified prior to the first vaccination dose (V1D0) and seven days after the second dose (V2D7) using anti-Spike IgG titers and antibody neutralization assays. Anti-spike IgG titers were additionally quantified six months after initial vaccination. Clinical history and lab values in HD patients were obtained to identify predictors of vaccination response.

Results

Transcriptomic analyses demonstrated differing time courses of immune responses, with predominant T cell activity in controls one week after the first vaccination dose, compared to predominant myeloid cell activity in HD at this time point. HD demonstrated decreased metabolic activity and decreased antigen presentation compared to controls after the second vaccination dose. Anti-spike IgG titers and neutralizing function were substantially elevated in both controls and HD at V2D7, with a small but significant reduction in titers in HD groups (p < 0.05). Anti-spike IgG remained elevated above baseline at six months in both subject groups. Anti-spike IgG titers at V2D7 were highly predictive of 6-month titer levels. Transcriptomic biomarkers after the second vaccination dose and clinical biomarkers including ferritin levels were found to be predictive of antibody development.

Conclusion

Overall, we demonstrate differing time courses of immune responses to the BTN162b2 mRNA COVID-19 vaccination in maintenance hemodialysis subjects (HD) comparable to healthy controls (HC) and identify transcriptomic and clinical predictors of anti-Spike IgG titers in HD. Analyzing vaccination as an in vivo perturbation, our results warrant further characterization of the immune dysregulation of end stage renal disease (ESRD).

Funding

F30HD102093, F30HL151182, T32HL144909, R01HL138628This research has been funded by the University of Illinois at Chicago Center for Clinical and Translational Science (CCTS) award UL1TR002003.

Transcriptomic meta-analysis reveals unannotated long non-coding RNAs related to the immune response in sheep

Long non-coding RNAs (lncRNAs) are involved in several biological processes, including the immune system response to pathogens and vaccines. The annotation and functional characterization of lncRNAs is more advanced in humans than in livestock species. Here, we take advantage of the increasing number of high-throughput functional experiments deposited in public databases in order to uniformly analyse, profile unannotated lncRNAs and integrate 422 ovine RNA-seq samples from the ovine immune system. We identified 12302 unannotated lncRNA genes with support from independent CAGE-seq and histone modification ChIP-seq assays. Unannotated lncRNAs showed low expression levels and sequence conservation across other mammal species. There were differences in expression levels depending on the genomic location-based lncRNA classification. Differential expression analyses between unstimulated and samples stimulated with pathogen infection or vaccination resulted in hundreds of lncRNAs with changed expression. Gene co-expression analyses revealed immune gene-enriched clusters associated with immune system activation and related to interferon signalling, antiviral response or endoplasmic reticulum stress. Besides, differential co-expression networks were constructed in order to find condition-specific relationships between coding genes and lncRNAs. Overall, using a diverse set of immune system samples and bioinformatic approaches we identify several ovine lncRNAs associated with the response to an external stimulus. These findings help in the improvement of the ovine lncRNA catalogue and provide sheep-specific evidence for the implication in the general immune response for several lncRNAs.

Vaccine-induced time- and age-dependent mucosal immunity to gastrointestinal parasite infection

Individuals vary broadly in their response to vaccination and subsequent challenge infection, with poor vaccine responders causing persistence of both infection and transmission in populations. Yet despite having substantial economic and societal impact, the immune mechanisms that underlie such variability, especially in infected tissues, remain poorly understood. Here, to characterise how antihelminthic immunity at the mucosal site of infection developed in vaccinated lambs, we inserted gastric cannulae into the abomasa of three-month- and six-month-old lambs and longitudinally analysed their local immune response during subsequent challenge infection. The vaccine induced broad changes in pre-challenge abomasal immune profiles and reduced parasite burden and egg output post-challenge, regardless of age. However, age affected how vaccinated lambs responded to infection across multiple immune pathways: adaptive immune pathways were typically age-dependent. Identification of age-dependent and age-independent protective immune pathways may help refine the formulation of vaccines, and indicate specificities of pathogen-specific immunity more generally.

Systems Immunology Analyses Following Porcine Respiratory and Reproductive Syndrome Virus Infection and Vaccination

This study was initiated to better understand the nature of innate immune responses and the relatively weak and delayed immune response against porcine reproductive and respiratory syndrome virus (PRRSV). Following modified live virus (MLV) vaccination or infection with two PRRSV-2 strains, we analyzed the transcriptome of peripheral blood mononuclear cells collected before and at three and seven days after vaccination or infection. We used blood transcriptional modules (BTMs)-based gene set enrichment analyses. BTMs related to innate immune processes were upregulated by PRRSV-2 strains but downregulated by MLV. In contrast, BTMs related to adaptive immune responses, in particular T cells and cell cycle, were downregulated by PRRSV-2 but upregulated by MLV. In addition, we found differences between the PRRSV strains. Only the more virulent strain induced a strong platelet activation, dendritic cell activation, interferon type I and plasma cell responses. We also calculated the correlations of BTM with the neutralizing antibody and the T-cell responses. Early downregulation (day 0-3) of dendritic cell and B-cell BTM correlated to both CD4 and CD8 T-cell responses. Furthermore, a late (day 3-7) upregulation of interferon type I modules strongly correlated to helper and regulatory T-cell responses, while inflammatory BTM upregulation correlated more to CD8 T-cell responses. BTM related to T cells had positive correlations at three days but negative associations at seven days post-infection. Taken together, this work contributes to resolve the complexity of the innate and adaptive immune responses against PRRSV and indicates a fundamentally different immune response to the less immunogenic MLV compared to field strains which induced robust adaptive immune responses. The identified correlates of T-cell responses will facilitate a rational approach to improve the immunogenicity of MLV.

Transcriptome-wide changes in gene expression, splicing, and lncRNAs in response to a live attenuated dengue virus vaccine

The tetravalent dengue vaccine candidate, TAK-003, induces a functional antibody response, but the titers of antibodies against the four serotypes of the dengue virus (DENV) can vary. Here, through a transcriptomic analysis on whole blood collected from recipients of a two-dose schedule of TAK-003, we examine gene expression, splicing, and transcript isoform-level changes for both protein-coding and noncoding genes to broaden our understanding of the immune response. Our analysis reveals a dynamic pattern of vaccine-associated regulation of long noncoding RNAs (lncRNAs), differential splicing of interferon-stimulated gene exons, and gene expression changes related to multiple signaling pathways that detect viral infection. Co-expression networks isolate immune cell-type-related and interferon-response modules that represent specific biological processes that correlate with more robust antibody responses. These data provide insights into the early determinants of the variable immune response to the vaccine, highlighting the significance of splicing and isoform-level gene regulatory mechanisms in defining vaccine immunogenicity.

Identification of Differential Responses of Goat PBMCs to PPRV Virulence Using a Multi-Omics Approach

Peste des petits ruminants (PPR) is an acute transboundary infectious viral disease of small ruminants, mainly sheep and goats. Host susceptibility varies considerably depending on the PPR virus (PPRV) strain, the host species and breed. The effect of strains with different levels of virulence on the modulation of the immune system has not been thoroughly compared in an experimental setting so far. In this study, we used a multi-omics approach to investigate the host cellular factors involved in different infection phenotypes. Peripheral blood mononuclear cells (PBMCs) from Saanen goats were activated with a T-cell mitogen and infected with PPRV strains of different virulence: Morocco 2008 (high virulence), Ivory Coast 1989 (low virulence) and Nigeria 75/1 (live attenuated vaccine strain). Our results showed that the highly virulent strain replicated better than the other two in PBMCs and rapidly induced cell death and a stronger inhibition of lymphocyte proliferation. However, all the strains affected lymphocyte proliferation and induced upregulation of key antiviral genes and proteins, meaning a classical antiviral response is orchestrated regardless of the virulence of the PPRV strain. On the other hand, the highly virulent strain induced stronger inflammatory responses and activated more genes related to lymphocyte migration and recruitment, and inflammatory processes. Both transcriptomic and proteomic approaches were successful in detecting viral and antiviral effectors under all conditions. The present work identified key immunological factors related to PPRV virulence in vitro.

Rational Vaccine Design in Times of Emerging Diseases: The Critical Choices of Immunological Correlates of Protection, Vaccine Antigen and Immunomodulation

Vaccines are the most effective medical intervention due to their continual success in preventing infections and improving mortality worldwide. Early vaccines were developed empirically however, rational design of vaccines can allow us to optimise their efficacy, by tailoring the immune response. Establishing the immune correlates of protection greatly informs the rational design of vaccines. This facilitates the selection of the best vaccine antigens and the most appropriate vaccine adjuvant to generate optimal memory immune T cell and B cell responses. This review outlines the range of vaccine types that are currently authorised and those under development. We outline the optimal immunological correlates of protection that can be targeted. Finally we review approaches to rational antigen selection and rational vaccine adjuvant design. Harnessing current knowledge on protective immune responses in combination with critical vaccine components is imperative to the prevention of future life-threatening diseases.

Influence of genetics and the pre-vaccination blood transcriptome on the variability of antibody levels after vaccination against Mycoplasma hyopneumoniae in pigs

Background

The impact of individual genetic and genomic variations on immune responses is an emerging lever investigated in vaccination strategies. In our study, we used genetic and pre-vaccination blood transcriptomic data to study vaccine effectiveness in pigs.

Results

A cohort of 182 Large White pigs was vaccinated against Mycoplasma hyopneumoniae (M. hyo) at weaning (28 days of age), with a booster 21 days later. Vaccine response was assessed by measuring seric M. hyo antibodies (Ab) at 0 (vaccination day), 21 (booster day), 28, 35, and 118 days post-vaccination (dpv). Inter-individual variability of M. hyo Ab levels was observed at all time points and the corresponding heritabilities ranged from 0.46 to 0.57. Ab persistence was higher in females than in males. Genome-wide association studies with a 658 K SNP panel revealed two genomic regions associated with variations of M. hyo Ab levels at 21 dpv at positions where immunity-related genes have been mapped, DAB2IP on chromosome 1, and ASAP1, CYRIB and GSDMC on chromosome 4. We studied covariations of Ab responses with the pre-vaccination blood transcriptome obtained by RNA-Seq for a subset of 82 pigs. Weighted gene correlation network and differential expression analyses between pigs that differed in Ab responses highlighted biological functions that were enriched in heme biosynthesis and platelet activation for low response at 21 dpv, innate antiviral immunity and dendritic cells for high response at 28 and 35 dpv, and cell adhesion and extracellular matrix for high response at 118 dpv. Sparse partial least squares discriminant analysis identified 101 genes that efficiently predicted divergent responders at all time points. We found weak negative correlations of M. hyo Ab levels with body weight traits, which revealed a trade-off that needs to be further explored.

Conclusions

We confirmed the influence of the host genetics on vaccine effectiveness to M. hyo and provided evidence that the pre-vaccination blood transcriptome co-varies with the Ab response. Our results highlight that both genetic markers and blood biomarkers could be used as potential predictors of vaccine response levels and more studies are required to assess whether they can be exploited in breeding programs.

Contagious Bovine and Caprine Pleuropneumonia: a research community's recommendations for the development of better vaccines

Contagious bovine pleuropneumonia (CBPP) and contagious caprine pleuropneumonia (CCPP) are major infectious diseases of ruminants caused by mycoplasmas in Africa and Asia. In contrast with the limited pathology in the respiratory tract of humans infected with mycoplasmas, CBPP and CCPP are devastating diseases associated with high morbidity and mortality. Beyond their obvious impact on animal health, CBPP and CCPP negatively impact the livelihood and wellbeing of a substantial proportion of livestock-dependent people affecting their culture, economy, trade and nutrition. The causative agents of CBPP and CCPP are Mycoplasma mycoides subspecies mycoides and Mycoplasma capricolum subspecies capripneumoniae, respectively, which have been eradicated in most of the developed world. The current vaccines used for disease control consist of a live attenuated CBPP vaccine and a bacterin vaccine for CCPP, which were developed in the 1960s and 1980s, respectively. Both of these vaccines have many limitations, so better vaccines are urgently needed to improve disease control. In this article the research community prioritized biomedical research needs related to challenge models, rational vaccine design and protective immune responses. Therefore, we scrutinized the current vaccines as well as the challenge-, pathogenicity- and immunity models. We highlight research gaps and provide recommendations towards developing safer and more efficacious vaccines against CBPP and CCPP.

Defining immune correlates during latent and active chlamydial infection in sheep

Ovine enzootic abortion (OEA) caused by the obligate intracellular bacterial pathogen Chlamydia abortus (C. abortus), is an endemic disease in most sheep-rearing countries worldwide. Following infection, C. abortus establishes a complex host-pathogen interaction with a latent phase in non-pregnant sheep followed by an active disease phase in the placenta during pregnancy leading to OEA. Improved knowledge of the host-pathogen interactions at these different phases of disease will accelerate the development of new diagnostic tests and vaccines to control OEA. Current evidence indicates that cellular immunity is essential for controlling C. abortus infection. We have previously described a model of mucosal (intranasal) infection of non-pregnant sheep with C. abortus that replicates the latent and active phases of OEA. We have investigated antigen-specific recall responses of peripheral blood mononuclear cells (PBMC) in sheep infected with C. abortus via the intranasal route to determine how these change during the latent and active phases of disease. By analysing cytokines associated with the major CD4+ve Thelper (Th) cell subsets (Interferon-gamma (IFN-γ)/Th1; Interleukin (IL)-4/Th2; IL-17A/Th17; IL-10/Tregulatory), we show that there is selective activation of PBMC producing IFN-γ and/or IL-10 during the latent phase following infection. These cytokines are also elevated during the active disease phase and while they are produced by sheep that are protected from OEA, they are also produced by sheep that abort, highlighting the difficulties in finding specific cellular immunological correlates of protection for complex intracellular pathogens.

The antibody response induced FMDV vaccines in sheep correlates with early transcriptomic responses in blood

Foot and mouth disease (FMD) is a highly contagious viral disease with high economic impact, representing a major threat for cloven-hooved mammals worldwide. Vaccines based on adjuvanted inactivated virus (iFMDV) induce effective protective immunity implicating antibody (Ab) responses. To reduce the biosafety constraints of the manufacturing process, a non-replicative human adenovirus type 5 vector encoding FMDV antigens (Ad5-FMDV) has been developed. Here we compared the immunogenicity of iFMDV and Ad5-FMDV with and without the ISA206VG emulsion-type adjuvant in sheep. Contrasted Ab responses were obtained: iFMDV induced the highest Ab levels, Ad5-FMDV the lowest ones, and ISA206VG increased the Ad5-FMDV-induced Ab responses to protective levels. Each vaccine generated heterogeneous Ab responses, with high and low responders, the latter being considered as obstacles to vaccine effectiveness. A transcriptomic study on total blood responses at 24 h post-vaccination revealed several blood gene module activities correlating with long-term Ab responses. Downmodulation of T cell modules’ activities correlated with high responses to iFMDV and to Ad5-FMDV+ISA206VG vaccines as also found in other systems vaccinology studies in humans and sheep. The impact of cell cycle activity depended on the vaccine types, as it positively correlated with higher responses to iFMDV but negatively to non-adjuvanted Ad5-FMDV. Finally an elevated B cell activity at 24 h correlated with high Ab responses to the Ad5-FMDV+ISA206VG vaccine. This study provides insights into the early mechanisms driving the Ab response induced by different vaccine regimens including Ad5 vectors and points to T cell modules as early biomarker candidates of different vaccine-type efficacy across species.

Foot and mouth disease virus (FMDV) is a serious pathogen of cloven hoofed mammals and is of high economic and veterinary importance. Inactivated vaccine (iFMDV) is effective but difficult to produce because of high biosafety level requirements; non-replicating adenovirus vectors carrying key FMDV antigens (Ad5-FMDV) might therefore represent an attractive alternative. Isabelle Schwartz-Cornil and colleagues use sheep to systematically compare vaccination with adjuvanted iFMDV, adjuvanted Ad5-FMDV, or non-adjuvanted Ad5-FMDV. All vaccines produce neutralizing antibody responses which are stable to at least one year, however the iFMDV group elicits the strongest response, followed by the adjuvanted Ad5-FMDV. Ad5-FMDV alone produces weak antibody titers. Blood transcriptomic analysis performed in the first 24 h following vaccination identifies a reduced T cell gene expression module as a correlate of high neutralizing antibody titers. Blood gene expression might therefore offer insights into the mechanistic underpinnings of humoral immunity as well as provide useful biomarker correlates of protection.

Blood transcriptome profile induced by an efficacious vaccine formulated with salivary antigens from cattle ticks

Ticks cause massive damage to livestock and vaccines are one sustainable alternative for the acaricide poisons currently heavily used to control infestations. An experimental vaccine adjuvanted with alum and composed by four recombinant salivary antigens mined with reverse vaccinology from a transcriptome of salivary glands from Rhipicephalus microplus ticks was previously shown to present an overall efficacy of 73.2% and cause a significant decrease of tick loads in artificially tick-infested, immunized heifers; this decrease was accompanied by increased levels of antigen-specific IgG1 and IgG2 antibodies, which were boosted during a challenge infestation. In order to gain insights into the systemic effects induced by the vaccine and by the tick challenge we now report the gene expression profile of these hosts' whole-blood leukocytes with RNA-seq followed by functional analyses. These analyses show that vaccination induced unique responses to infestations; genes upregulated in the comparisons were enriched for processes associated with chemotaxis, cell adhesion, T-cell responses and wound repair. Blood transcriptional modules were enriched for activation of dendritic cells, cell cycle, phosphatidylinositol signaling, and platelets. Together, the results indicate that by neutralizing the tick's salivary mediators of parasitism with vaccine-induced antibodies, the bovine host is able to mount normal homeostatic responses that hinder tick attachment and haematophagy and that the tick otherwise suppresses with its saliva.

Effect of simultaneous administration of foot-and-mouth disease (FMD) and anthrax vaccines on antibody response to FMD in sheep

Purpose

Foot-and-mouth disease (FMD) and anthrax are important diseases in sheep. Vaccination is a favorable strategy against both infections. Simultaneous administration of vaccines does generally not impede the immune responses of each other, although there are some exceptions, and it may help reduce the labor and costs of vaccination as well as distress on animals. Although oil adjuvant FMD vaccine has been tried with live anthrax vaccine in cattle, there are no reports on the simultaneous use of both vaccines in sheep.

Materials and Methods

In this study, FMD seronegative sheep were used to investigate the impact of the simultaneous vaccination of FMD and anthrax on FMD antibody titers of sheep. Virus neutralization test and liquid phase blocking enzyme-linked immunosorbent assay were used to determine the antibody response to the FMD vaccine.

Results

The results demonstrated that both vaccines can be used simultaneously without any interference with the FMD response. Moreover, the simultaneous administration with anthrax vaccine had a stimulating effect on the early (day 7 post-vaccination) virus neutralization antibody response to the FMD vaccine.

Conclusion

The simultaneous use of the FMD and anthrax vaccines did not hinder the response to the FMD vaccine in sheep.

Systems Immunology Characterization of Novel Vaccine Formulations for Mycoplasma hyopneumoniae Bacterins

We characterized five different vaccine candidates and a commercial vaccine in terms of safety, immunogenicity and using a systems vaccinology approach, with the aim to select novel vaccine candidates against Mycoplasma hyopneumoniae. Seven groups of six M. hyopneumoniae-free piglets were primo- and booster vaccinated with the different experimental bacterin formulations, the commercial vaccine Hyogen® as a positive control or PBS as a negative control. The experimental bacterin was formulated with cationic liposomes + c-di-AMP (Lipo_AMP), cationic liposomes + Toll-like receptor (TLR) 2/1, TLR7, and TLR9 ligands (TLR ligands; Lipo_TLR), micro-particles + TLR ligands (PLGA_TLR), squalene-in-water emulsion + TLR ligands (SWE_TLR), or DDA:TDB liposomes (Lipo_DDA:TDB). Lipo_DDA:TDB and Lipo_AMP were the most potent in terms of serum antibody induction, and Lipo_DDA:TDB, Lipo_AMP, and SWE_TLR significantly induced Th1 cytokine-secreting T-cells. Only PLGA_TLR appeared to induce Th17 cells, but was unable to induce serum antibodies. The transcriptomic analyses demonstrated that the induction of inflammatory and myeloid cell blood transcriptional modules (BTM) in the first 24 h after vaccination correlated well with serum antibodies, while negative correlations with the same modules were found 7 days post-vaccination. Furthermore, many cell cycle and T-cell BTM upregulated at day seven correlated positively with adaptive immune responses. When comparing the delivery of the identical TLR ligands with the three formulations, we found SWE_TLR to be more potent in the induction of an early innate immune response, while the liposomal formulation more strongly promoted late cell cycle and T-cell BTM. For the PLGA formulation we found signs of a delayed and weak perturbation of these BTM. Lipo_AMP was found to be the most potent vaccine at inducing a BTM profile similar to that correlating with adaptive immune response in this and other studies. Taken together, we identified four promising vaccine candidates able to induce M. hyopneumoniae-specific antibody and T-cell responses. In addition, we have adapted a systems vaccinology approach developed for human to pigs and demonstrated its capacity in identifying early immune signatures in the blood relating to adaptive immune responses. This approach represents an important step in a more rational design of efficacious vaccines for pigs.