Abstract
Mycoplasma hyopneumoniae is the most costly pathogen for swine production. Although several studies have focused on the host-bacterium association, little is known about the changes in gene expression of swine cells upon infection. To improve our understanding of this interaction, we infected swine epithelial NPTr cells with M. hyopneumoniae strain J to identify differentially expressed mRNAs and miRNAs. The levels of 1,268 genes and 170 miRNAs were significantly modified post-infection. Up-regulated mRNAs were enriched in genes related to redox homeostasis and antioxidant defense, known to be regulated by the transcription factor NRF2 in related species. Down-regulated mRNAs were enriched in genes associated with cytoskeleton and ciliary functions. Bioinformatic analyses suggested a correlation between changes in miRNA and mRNA levels, since we detected down-regulation of miRNAs predicted to target antioxidant genes and up-regulation of miRNAs targeting ciliary and cytoskeleton genes. Interestingly, most down-regulated miRNAs were detected in exosome-like vesicles suggesting that M. hyopneumoniae infection induced a modification of the composition of NPTr-released vesicles. Taken together, our data indicate that M. hyopneumoniae elicits an antioxidant response induced by NRF2 in infected cells. In addition, we propose that ciliostasis caused by this pathogen is partially explained by the down-regulation of ciliary genes.
Similar content being viewed by others
Introduction
Respiratory diseases are among the major health problems in the pig farming industry. Mycoplasma hyopneumoniae is the causative agent of swine enzootic pneumonia, a chronic respiratory disease that affects herds worldwide. Although M. hyopneumoniae does not cause high mortality, it is considered the most expensive pathogen for swine production1. This is mainly due to the costs of treatment and vaccination and to losses related to decreased animal performance. In addition, M. hyopneumoniae is essential for the establishment of secondary pathogens in the host, which leads to a significant increase in mortality2. M. hyopneumoniae attaches to the cilia of the tracheal epithelial cells with participation of adhesins3, resulting in ciliostasis and cell death4. Besides adhesins, virulence factors are not well understood in this bacterium. Nevertheless, a recent study from our group indicated hydrogen peroxide production from glycerol and myo-inositol metabolism as important traits that might be related with pathogenesis and with the predominance of M. hyopneumoniae in the swine respiratory tract5.
MicroRNAs (miRNAs) belong to a class of small non-coding RNAs (ncRNAs) of 18–24 nucleotides (nt) in part responsible for post-transcriptional gene regulation in eukaryotes. These evolutionarily conserved molecules influence fundamental biological processes, including cell proliferation, differentiation, apoptosis, immune response, and metabolism6, 7. The binding of miRNAs to target mRNAs changes the mRNA stability and translation efficiency8, leading to degradation, suppression or up-regulation of the target mRNAs6, 9. Interactions between miRNA and mRNA are complex; one single miRNA can target a large number of genes belonging to diverse functional groups. Alternatively, the 3’-UTR of a single mRNA can be targeted by multiple miRNAs10, 11. By modulating miRNA abundance, it is thus possible to fine-tune the expression of proteins within the cell in a very precise manner6, 11.
Recently, it was widely reported that miRNAs can be packed into exosomes and transferred to neighboring or distant cells to regulate cell function7, 12,13,14. Exosomes are small membrane vesicles (50–150 nm) released from eukaryotic cells both constitutively and upon induction, under normal and pathological conditions13, 15. These vesicles are involved in several cellular functions and have the potential to selectively interact with specific target cells16, 17. In addition to miRNAs, exosomes can transmit information among cells by transferring proteins, lipids and nucleic acids that seem to be selected non-randomly, with some specific populations of molecules being preferentially packaged into the vesicles13, 15. As an efficient cellular signaling and communication system, the release of exosomes by infected host cells has been recognized as a common phenomenon, in some cases beneficial to the host and in others beneficial to the pathogen18.
Host-pathogen interactions result in signaling and physiological modifications in the host cells that induce differential miRNA expression and miRNA-mediated post-transcriptional regulation of genes involved in immune response and several other cellular pathways19, 20. Therefore, simultaneous identification of differentially expressed miRNAs and mRNAs provides a comprehensive view on host-pathogen interactions during the infection and the disease establishment process. In recent years, efforts have been made to identify miRNAs regulated by infection in several mammalian hosts8, 20,21,22. However, the identification of miRNAs during infection of swine cells with M. hyopneumoniae has not been investigated so far. To improve our understanding on the M. hyopneumoniae-host interaction, we sequenced and analyzed both the mRNAs and miRNAs of a swine tracheal epithelial cell line infected with M. hyopneumoniae strain J. In addition, we identified miRNAs differentially expressed (DE) in the extracellular milieu and in exosome-like vesicles released by the infected cells, which play an important role in cell-cell communication and in the dissemination of host and pathogen-derived molecules during infection15. The simultaneous identification of miRNAs and mRNAs will help us draw a full picture of the changes in gene expression and the possible regulatory mechanisms of host cells during the disease establishment.
Results and discussion
M. hyopneumoniae strain J adhered to NPTr cells
To analyze the differential expression of New-born Pig Trachea (NPTr) cells during the infection with M. hyopneumoniae, we first observed the infection by immunofluorescence microscopy. These analyses were performed to detect the adherence of M. hyopneumoniae strain J to NPTr cells. Supplementary Figure S1 shows the co-localization of M. hyopneumoniae with NPTr cells, corroborating the success of the infection. As few bacteria adhered to the cells within a short period of time (1 vs. 24 h), we chose to analyze the transcriptional alterations of NPTr cells at 24 h post-infection. M. hyopneumoniae strain J was chosen because previous infection assays showed that highly virulent strains, such as 7448 or 7422, damaged the host cells and we were not able to recover RNA with good sequencing quality. Although this strain is considered attenuated and incapable of causing disease in vivo23, 24, our results show that M. hyopneumoniae is capable of adhering to the swine epithelial cells, as previously reported by Burnett et al.25.
mRNA expression profiles
A total of 6 mRNA libraries were generated with the Illumina HiSeq2500 platform. A summary of the samples is provided in Table 1 and a diagram to explain the experimental design is given in Supplementary Fig. S2. The raw reads were submitted to the NCBI Sequence Read Archive under accession number PRJNA545822. After removing adapters and filtering low quality reads, mRNA-seq yielded around 40 million paired-end reads in all 6 samples (approx. 97% of raw reads). Trimming and mapping information for each mRNA sample is available in Supplementary Table S1. Around 85% of the filtered reads were mapped against the porcine genome (Sscrofa10.2—Ensembl release 89) and 40% against annotated genes (Supplementary Table S1).
In the present study, we detected a total of 20,274 (out of 23,215) genes expressed with at least 10 counts across the 6 mRNA libraries from NPTr cells (Supplementary Table S2). Since we used an attenuated strain of M. hyopneumoniae to infect the swine cells, we did not expect extreme changes in gene expression. Indeed, most of the significant log2 of Fold Changes (LogFC) of the DE genes analyzed in our results were below 1 (for up-regulated genes) or above -1 (for down-regulated genes) (Supplementary Table S3). Nevertheless, even if the strain used (J) is considered attenuated, it generated a considerable immune response. A similar situation occurs in vaccine development, in which attenuated strains are used to generate immune response without causing lethal effects. Thus, we believe that infections with highly virulent strains of M. hyopneumoniae would be able to generate a more pronounced effect on gene expression than what was observed here. Regardless, our results provide a general overview of M. hyopneumoniae infection in the host cells.
In total, we detected 1,268 DE genes (p-adj \(< 0.05\), 517 up-regulated and 751 down-regulated), from which 502 were common to two well known methods for the detections of DE genes, 721 were exclusive to DESeq2 and 45 were exclusive to EdgeR. Information from the top 15 up-regulated and down-regulated DE genes is provided in Table 2. The complete DE gene results are found in Supplementary Table S3. The results of gene ontology (GO) enrichment analysis for up-regulated and down-regulated genes in DESeq2 are shown separately in Fig. 1 (p-adj \(< 0.05\) and absolute fold enrichment \(\ge \) 0.1) and the complete results are found in Supplementary Table S4. For the up-regulated genes, the enriched terms either in biological process (BP), molecular function (MF) or cellular component (CC) were related to protein synthesis (translation/ribosome), oxidation reduction activity and cell-cell communication (such as exosomes, anchoring and focal adhesion functions) (Fig. 1A). For the down-regulated genes, the majority of the overrepresented terms in all three GO categories were related to cell cycle, cell division, cilia and cytoskeleton (Fig. 1B).
Significantly enriched Gene Ontology functions of DE mRNAs. A. Up-regulated genes were enriched in terms related to protein synthesis, oxidation-reduction activity and cell-cell communication (such as exosomes, anchoring and focal adhesion), either in biological process (BP), molecular function (MF) or cellular component (CC). B. Down-regulated genes showed an enrichment of terms related to cytoskeleton, ciliary function, cell cycle and cell division in all the three categories of GO.
M. hyopneumoniae elicited an antioxidant response and induced the accumulation of NRF2 in the nuclei of NPTr cells
Among the up-regulated genes, we found several ones related to immune response and inflammation, such as C3 complement, SAA3, chemokines (CXCL2 and CCL20) and galectins (LGALS2 and LGALS8) (Supplementary Table S3). Interestingly, we also detected 64 up-regulated genes related to redox homeostasis and antioxidant defense (Table 3). We observed that 46 out of 64 have already been described as targets of the nuclear factor erythroid 2-related factor 2 (NRF2) in closely related species (Table 3). This transcription factor is involved in the protection of the cell against oxidative damage through transcription activation of cytoprotective genes26, 27. More specifically, several studies have shown the protective role of NRF2 in bacterial lung infections in rodents, being a critical factor for assembling the innate immune response in the host28,29,30,31.
Indeed, M. hyopneumoniae infection induced the expression of genes related to the two biggest redox systems in NPTr cells - glutathione (GGT1, GGT5, GGLC, GSR) and thioredoxin (TXNRD1, PRDX5, PRDX6) -, and also genes coding for NADPH-regenerating enzymes (used by the aforementioned redox systems), as analogously reported for activation of NRF2 targets in mice32. Moreover, a number of antioxidant genes and genes coding for detoxification enzymes (such as the AKR gene family, NQO1, HMOX and GST) which were up-regulated during M. hyopneumoniae infection, were also reported to be activated by NRF232 (Table 3). We performed a Fisher’s exact test to compare the proportion of genes within the genome expected to be targets of NRF2 with the proportion of up-regulated genes putatively regulated by this transcription factor. The results indicate an extremely significant correlation of an NRF2 target being an upregulated gene (p-value \(< 0.00001\), Supplementary Table S5). The up-regulation of several targets of NRF2 was further validated by reverse-transcription quantitative PCR (RT-qPCR) and the results are available in Supplementary Fig. S3.
One of the many compounds known to activate the NRF2 pathway is hydrogen peroxide33. The production of hydrogen peroxide was previously described to have an important cytotoxic effect in several Mycoplasma species, such as M. pneumoniae34 and M. mycoides subsp. mycoides35. More specifically, the cytotoxity of M. mycoides subsp. mycoides was correlated to the bacterium’s ability to translocate hydrogen peroxide directly into the host cell36. Previously, we have identified that M. hyopneumoniae was capable of producing this toxic product from glycerol metabolism5. Although the production of hydrogen peroxide by the M. hyopneumoniae strain J was not detected in that work, novel analyses indicate that in conditions where glucose is scarce, the attenuated strain J is able to produce this toxic metabolite (Supplementary Fig. S4). In addition, we were able to detect hydrogen peroxide in the medium of the NPTr cells infected with M. hyopneumoniae in the presence of glycerol (Supplementary Fig. S4C). Although the production was higher in the cells infected with the pathogenic strain 7448, M. hyopneumoniae strain J was also able to produce the toxic product, indicating that both strains could potentially cause oxidative damage to the host cells. These results are in accordance with the gene expression of the putative enzyme responsible for hydrogen peroxide production (GlpO) in M. hyopneumoniae strain J, which did not differ from the pathogenic strains5. During infection in general, bacteria need to compete with host cells (and other organisms in the environment) for glucose and other energy sources. For instance, in M. pneumoniae, Halbedel et al. (2004)37 showed that even though glucose was the most efficient carbon source for biomass yield (as is the case for M. hyopneumoniae), the authors propose that glycerol is not only a carbon source, but it could be used by this species as an indicator that they reached their preferred ecological niche, a lipid-rich mucosal surface. Thus, it is plausible to say that M. hyopneumoniae also uses glycerol as a carbon source in vivo, however whether this is a result of competition against other fast-growing bacteria or if this species is targeting a glycerol-rich niche, we cannot affirm at this point. What we can affirm is that glycerol is not the most efficient carbon source for both M. pneumoniae and M. hyopneumoniae, and yet these species make use of it for energy uptake and also for other advantages, such as hydrogen peroxide production. Thus, it is possible that the hydrogen peroxide produced by M. hyopneumoniae strain J as a result of glycerol metabolism might be one of the triggers that activated the transcription factor NRF2.
It has been demonstrated in other species that NRF2 is largely regulated at the post-transcriptional level by its sub-cellular distribution, which is controlled by the Kelch-like ECH-associated protein (KEAP1). Under normal conditions, a portion of NRF2 is retained in the cytoplasm via its interaction to KEAP1 and it is subsequently ubiquitinated and degraded by the proteasome. In response to oxidative stress, reactive cysteines in KEAP1 are modified, generating conformational changes in the complex and releasing NRF2, which is translocated and accumulates into the nucleus27, 32, 38. In our analysis, we did not detect a difference in the expression of the mRNAs of NRF2 (base mean expression around 8000 reads, a logFC between infected and non-infected conditions close to zero and a non-significant adjusted p-value of 0.99) and KEAP1 (base mean expression around 4500 reads, logFC between infected and non-infected conditions close to zero and non-significant adjusted p-value of 0.94). However, we were able to demonstrate by confocal immunofluorescence microscopy (Fig. 2) that both the attenuated and the virulent strains of M. hyopneumoniae induced a statistically significant accumulation of this transcription factor in the nuclei of NPTr cells (Fig. 2B). In the nucleus, NRF2 is able to recognize and bind to antioxidant response element (ARE) motifs in the promoter region of target genes, activating their transcription39. In this work, we detected at least one conserved ARE sequence upstream the start codon in the promoter regions of 44 out of the 46 NRF2 predicted targets (with stringent search to TG/TAnnnnGC) with the use of fuzznuc software from EMBOSSv6.6.0 package40 (Supplementary Table S6).
Localization and expression pattern of the NRF2 protein in cells infected with M. hyopneumoniae.A. Results of the immunofluorescence microscopy analysis indicating the accumulation of NRF2 in the nuclei of infected cells after 1 h of incubation. Both attenuated (J) and virulent (7448) strains of M. hyopneumoniae induced the accumulation of the transcription factor in the nuclei of the epithelial cells. Eukaryotic cell NRF2 was labeled with anti-NRF2 antibody (green) and nuclei were stained with DAPI (blue). Scale bars: \(50\,\mu \hbox {m}\). B. Total NRF2 green fluorescence per nucleus of control and infected cells. Both M. hyopneumoniae strains significantly increased the concentration of NRF2 in the nuclei of infected cells. Boxplots represent the measures of at least 25 nuclei in each condition. Outliers are represented as black dots. CTCF—corrected total cell fluorescence. NPTr - uninfected cells. MHP J - NPTr cells infected with M. hyopneumoniae strain J. MHP 7448 - NPTr cells infected with M. hyopneumoniae strain 7448. **** \(\hbox {p} < 0.0001.\)
In this way, it seems that M. hyopneumoniae infection has the potential to cause oxidative stress to the host cells, which in turn activate antioxidant response genes induced by NRF2 to fight the infection and maintain cellular homeostasis. We believe that this oxidative stress was in part related to the hydrogen peroxide produced by this bacterium, although more experiments are needed to prove the association between this mechanism and the up-regulation of antioxidant response genes.
M. hyopneumoniae induced down-regulation of cytoskeleton and ciliary genes as well as a decrease of actin stress fibers in NPTr cells
We identified several down-regulated genes related to ciliary function, cytoskeleton and cell cycle/cell division. The impairment of the ciliary motility is a well-known effect caused by several Mycoplasma respiratory species41, such as M. pneumoniae and M. gallisepticum42, 43. It is also well established that M. hyopneumoniae attaches to cilia of epithelial cells and promotes ciliostasis and loss of cilia, causing damage to the mucociliary apparatus4, 24. Therefore, our hypothesis is that one of the reasons for epithelial damage, besides physical adhesion, could be associated with modulations in gene expression induced by the infection, which is a running hypothesis for at least two epithelial pathogens from the genus Mycobacterium44. In agreement with this hypothesis, we were able to identify the down-regulation of genes coding for axonemal dyneins (DNAH11, DNAH12, DNAI2, DNAL1), which are essential for the ciliary motility45. In addition, genes necessary for axonemal dynein assembly (DYX1C1)46, genes related to ciliogenesis (CEP162, DCDC2, MACF1, IFT57)47,48,49,50, ciliary polarization (INTU)51, ciliary beating (MYO1D)52 and several others in which the mutation or knockout is associated with ciliopathies (LRRC6, MNS1, AK7)53,54,55 were down-regulated in the infected cells (Table 4). Interestingly, in line with our results, a recent study compared the transcriptional response of unvaccinated and vaccinated chicken infected with M. gallisepticum and the authors identified enrichment of GO terms in down-regulated genes related to cilia and cytoskeleton in unvaccinated animals56. Protein functions encoded by the top down-regulated genes were involved in microtubule assembly and stability, axonemal dynein complex assembly, and formation and motor movement of cilia, indicating that at least in one Mycoplasma species the ciliary damage caused by infection could be also explained by the down-regulation of genes involved in the ciliary function.
Besides ciliary genes, we also detected the down-regulation of cytoskeleton-related genes, both from microtubules and actin filaments, involved in the organization, rearrangement and stability of these structures (Table 4). It was previously described that the intracellular species M. penetrans is able to trigger reorganization of the host cell cystoskeleton, promoting aggregation of tubulin and \(\alpha \)-actinin and condensation of phosphorylated proteins57. To investigate if M. hyopneumoniae indeed affected the host cell cytoskeleton, we verified the organization of actin fibers in infected cells by confocal immunofluorescence microscopy. The actin stress fibers were notably less evident in the infected cells, as opposed to the control condition, in which they were abundant and more evenly distributed (Fig. 3). However, whenever present in the infected conditions, these stress fibers were either disorganized and/or at the periphery of the cells. These results corroborate studies from Raymond et al. (2018)58, which suggested that this species induces cytoskeletal rearrangements in the porcine respiratory tract. Although actin is not a major component of the ciliary axoneme, actin cytoskeleton has been implicated in every stage of ciliogenesis and many aspects of ciliary function59, directly associating these two down-regulated functions. In addition, it has recently been shown that M. hyopneumoniae expresses surface-accessible actin-binding proteins and that the host’s extracellular actin may act as a receptor for this bacterium in PK-15 epithelial cells60, indicating its importance for successful infection. Furthermore, the reduction of visible actin stress fibers caused by M. hyopneumoniae may also be related to the activation of NRF2, since the actin cytoskeleton is a scaffold necessary to maintain the transcription factor in the cytoplasm61.
Organization of actin fibers in cells infected with M. hyopneumoniae. Results of the immunofluorescence microscopy analysis indicating the reduction and change in the pattern of actin stress fibers in infected cells. Eukaryotic cell actin was labeled with phalloidin (red) and nuclei were stained with DAPI (blue). Both attenuated (J) and virulent (7448) strains of M. hyopneumoniae altered the organization and abundance of actin fibers after 24 h of infection, however this effect can already be seen after 1 h of incubation (not shown). NPTr - uninfected cells. MHP J - NPTr cells infected with M. hyopneumoniae strain J. MHP 7448 - NPTr cells infected with M. hyopneumoniae strain 7448. Scale bars: \(50\,\mu \hbox {m}.\)
We also identified the down-regulation of cytoskeleton-related genes that play a role during cell division, such as BUB1B, CENP-I, NEK2 and SPAG5 (Supplementary Table S3). Many of these genes are involved with the mitotic spindle and chromosome segregation62,63,64,65,66,67,68,69. Genes encoding microtubule dependent motor proteins that physically affect chromosome segregation, such as kinesins (usually up-regulated during mitosis)70, as well as genes related to cell cycle progression were down-regulated during infection (Supplementary Table S3). These results suggest a repression of cell division of infected cells. The manipulation of cell cycle by pathogens has been extensively reported, with different pathogens being able to arrest different points of the cell cycle71,72,73,74. Therefore, it is possible that M. hyopneumoniae infection also interferes with the host cell cycle.
miRNA expression profiles
A total of 14 small RNA (sRNA) libraries were generated with the Illumina HiSeq2500 platform (Supplementary Table S7). A complete description of the samples is provided in Table 1 and Fig. S2. The raw reads were submitted to the NCBI Sequence Read Archive under accession number PRJNA545822. After removing adapters and filtering low quality reads, sRNA-seq yielded from 1 to 21 million single-end clean reads. Trimming and mapping information for each sRNA sample is given in Supplementary Table S7.
For intracellular sRNAs (INTRA: samples S1 to S6), we were able to map around 93% of the reads against the porcine genome. We also mapped all sRNA reads against the M. hyopneumoniae strain 7448 genome and samples had no more than 0.05% of unique mycoplasmal reads (Supplementary Table S7). All intracellular sRNA samples had similar sequencing depth and showed a typical miRNA distribution length curve ranging from 18 to 26 nt, with a peak at 22 nt (Supplementary Fig. S5A). We also performed a homology analysis of raw reads and we provide the percentage of intracellular clean reads clustered by classical types of sRNAs (Supplementary Fig. S5B), showing an enrichment of miRNA-homolog reads (around 40%).
As expected, total extracellular sRNA samples (EXTRA: S7 to S10) had more RNA degradation due to the presence of RNAses and degraded mRNAs in the extracellular environment. Extracellular sRNA samples of cells infected with the bacterium (S9 and S10) had up to 20% of reads mapped to M. hyopneumoniae, as they were presumed to have more remnants of mycoplasmal cells (Supplementary Table S7). However, the great majority of reads mapping to the bacteria were mainly product of mRNA degradation, as they did not have any specific signature for sRNAs when compared to previous results from Siqueira et al. (2016)75. In this way, we filtered out (i) sequences \(< 18\) nt and (ii) sequences mapped to the M. hyopneumoniae genome. Sample S8 (EXTRA) had a sequencing depth smaller than its replicate (S7); however, the distribution of counts was overall similar between replicates and we included it in the further analyses.
Extracellular exosomal sRNA samples (EXO: S11 and S12) had some RNA degradation (\(<18\) nt), but maintained a pronounced peak at 22 nt. Extracellular sRNAs from vesicle-free supernatant (SN: samples S13 and S14) were more problematic, probably due to (i) too many pre-processing steps and (ii) the presence of RNAses in the extracellular environment. Sample S14 did not yield a minimum amount of reads necessary for the subsequent steps and this condition (SN) was not used for further analyses.
Annotation of known and novel miRNAs
In order to perform miRNA prediction with mirDeep276, we only took into account reads from intracellular samples. We predicted a total of 1,041 miRNAs, which were further clustered into 773 groups. From these, 478 were completely novel (Supplementary Table S8). We created a Sus scrofa miRNA database (ssc-miRNA-DB) with 1,906 different entries from three different sources: 411 known miRNAs in Sus scrofa from miRBase (release 21)77, 722 annotated by Martini et al., (2014)78 and the 773 clusters of mature miRNAs predicted by our analysis. More than 50% of intracellular sRNA clean reads were aligned against the ssc-miRNA-DB (Supplementary Table S7). All other samples were also mapped against this database. The complete pipeline used in this study for miRNA prediction is described in Supplementary Fig. S6.
miRNA differential expression
In all sRNA samples, we detected (with at least 10 counts across all libraries) 290 from the 411 miRNAs from miRBase and 214 of the miRNAs described by Martini et al., (2014)78. Also, 263 predicted miRNAs had more than 10 counts across all libraries. We only took into account the 491 miRNAs that had more than 50 counts across all libraries for differential expression analysis (Supplementary Table S9).
In total, we identified 170 DE miRNAs (121 up-regulated, 48 down-regulated and 1 ambiguous). Table 5 shows 10 selected up-regulated and down-regulated miRNAs detected in this study, and the complete list of DE miRNAs is summarized in Supplementary Table S10. Several homologs of these miRNAs have already been linked to bacterial infection response or immune system response in the literature and these references are listed in Table 5. With the exception of one ambiguous case (ssc-miR-9842-5p), whenever a miRNA was detected as DE in more than one condition (INTRA, EXTRA or EXO), the change in expression (either up- or down-regulated) was in accordance between them (Supplementary Table S10).
Targets of DE miRNAs were enriched in genes related to redox homeostasis, translation and cytoskeleton
In order to better understand the biological functions that could be involved with the 170 DE miRNAs, we performed different analyses to predict their potential targets. In this sense, all DE miRNAs had at least one DE gene as a predicted target. The complete pipeline used for target prediction in this study is provided in Supplementary Fig. S7.
We predicted a total of 79,276 interaction pairs between miRNAs and mRNAs. In this way, based only on the predictions of the software used here, a miRNA could potentially target, on average, 465 genes in the entire porcine genome. However, we only considered as targets the mRNAs that were detected as DE in this study, which significantly decreased our list to a total of 4,287 interaction pairs. We chose to focus on anticorrelations between miRNA and mRNA expression (and kept 1,939 interaction pairs), as the main mode of action of miRNAs is a destabilization of mRNAs79 and the fact that most of experimental validations in the literature are related to interaction pairs with inversed regulation. In this context, a permissive interaction is generally described as one that occurs between a down-regulated miRNA and an up-regulated mRNA, while a repressive interaction is one in which the miRNA is up-regulated with consequent down-regulation of the target mRNA78 (Fig. 4A). In our results, permissive interactions represented 267 genes and 50 miRNAs, while repressive interactions occurred between 425 genes and 121 miRNAs, accounting for 598 permissive and 1341 repressive interactions between DE mRNAs and miRNAs. Supplementary Table S11 summarizes the permissive and repressive target pairs predicted in this study.
Correlation of GO terms between DE mRNAs and targets of DE miRNAs. A A permissive interaction occurs between a down-regulated miRNA (depicted in blue) and an up-regulated mRNA (depicted in red), whereas a repressive interaction is one in which the miRNA is up-regulated (red) and its target mRNA is down-regulated (blue). B. We performed GO enrichment analyses with the complete up- and down-regulated lists of mRNAs and also with the subset of miRNA targets among each of these lists. C. A correlation of some of the GO terms from DE mRNAs and targets of DE miRNAs was detected, indicating that these functions might also be regulated by miRNAs.
We performed GO analyses to compare the DE genes and targets of DE miRNAs in order to investigate whether their functions could be correlated (Fig. 4B and Supplementary Fig. S8). Interestingly, we were able to identify a correlation between the enriched terms in miRNA targets with some of the GO terms detected in the mRNA up-regulated or down-regulated GO results (Fig. 4C). Target genes from permissive interactions were associated with ribosome/translation and oxidation-reduction activity, whereas target genes from the repressive interactions had enriched terms related to cytoskeleton and ciliary function (Fig. 4C). It is important to highlight the relevance of the miRNAs found in exosome-like vesicles and in extracellular samples in the identification of GO enriched terms of the miRNA targets, since only a small part of the permissive and repressive interactions involved intracellular miRNAs. Complete GO enrichemnt of the miRNA targets is found in Supplementary Table S12.
Regulatory network reconstruction and analysis
In order to gain a broader view of the host’s response to the presence of M. hyopneumoniae, we built a general regulatory network by integrating mRNA gene expression with predicted miRNA-mRNA interactions collected and analyzed in this work (Supplementary Fig. S9). We also included information about physical and genetic validated interactions from the BioGRID v3.4 database80, 81. As previously mentioned, at this point we only took into account interaction pairs that were either repressive or permissive. The global interaction network was composed by 774 nodes and 1965 arcs and our objective was to identify miRNA-mRNA expression patterns. We performed a functional analysis with the use of CLueGO82 and we also detected within the permissive interactions the enrichment of several processes related to immune response and inflammation (Supplementary Fig. S10).
Interaction networks of DE miRNAs and genes involved with cytoskeleton/cilia and redox homeostasis. Rectangles depict miRNAs and ellipses represent genes. Blue ellipses/rectangles indicate that the gene/miRNA was down-regulated and red indicate gene/miRNAs detected as up-regulated. BioGrid interactions are seen as green dashed arcs. These represent both physical and genetic interactions between either genes or proteins from such genes. A. Cytoskeleton/cilia interaction network. We detected many repressive interactions involving genes related to cytoskeleton and ciliary function (blue arcs). B. Redox homeostasis interaction network. Permissive interactions of genes related to redox homeostasis are shown as red arcs. The great majority of the miRNAs targeting these genes are DE only in exosomes (rectangles partially colored in blue, far right). NRF2 targets: genes described in other species to be activated by NRF2 are highlighted in grey.
Furthermore, we created two separate regulatory networks, one related to cytoskeleton and cilia (repressive interactions, Fig. 5A) and one related to redox homeostasis (permissive interactions, Fig. 5B) to better refine and understand the possible co-regulation of several of these targets. The networks show a high level of connectivity and we were able to detect miRNAs that interacted only with NRF2 activated targets (such as ssc-miR-31) (Table 6). Furthermore, we were able to detect miRNAs that seemed to regulate more generally redox homeostasis genes and also miRNAs whose targets came from sets of genes with distinct functions (such as glycolysis, immune system defense, ribosomes, among others), indicating that this response might be related to several other important functions within the cell. We believe that this network can be a powerful tool for analyzing the influence of M. hyopneumoniae on the host gene expression in future studies.
Permissive interactions were enriched in genes regulated by NRF2
Analyzing more in detail the permissive interactions between mRNAs and miRNAs, we identified several miRNAs targeting NRF2 regulated genes (Table 6). GGT1, for instance, is predicted to be targeted by 7 different miRNAs, while ssc-mir-31 is predicted to target 4 different genes induced by NRF2 (AKR1C4, AKR1CL1, GGT1, and TXNRD1).
It was previously reported that NRF2 can be regulated independently of KEAP1 by miRNAs in human breast cancer cells83. Our miRNA target analysis predicted NRF2 as target of three miRNAs down-regulated in exosome-like vesicles: ssc-mir-340, ssc-mir-19a, and ssc-mir-19b. The first was also predicted to target the transporter gene SLC6A6, while ssc-mir-19a and ssc-mir-19b were predicted to target GCLC, necessary for glutathione synthesis. In addition, miR-340 has been previously identified negatively regulating NRF2 expression in human84. Conversely, miR-19a and miR-19b have been reported to be down-regulated in the presence of hydrogen peroxide in rats85, 86, linking these miRNAs to a response to oxidative stress.
In the literature, other miRNAs were also shown to negatively regulate the expression of NRF2: miR-2883, miR-10187, miR-92a88, miR-27b89, and miR-34a90. All of them were down-regulated in our data and were predicted by this study to regulate NRF2 activated genes (Table 6). In this way, it seems that during M. hyopneumoniae infection, there is a global change in gene expression in an attempt to activate antioxidant genes, in association with the down-regulation of miRNAs that negatively regulate these genes.
Repressive interactions were related to cytoskeleton and cilia
As concerns the repressive interactions of miRNAs and mRNAs, we identified within the targets of up-regulated miRNAs an enrichment of genes related to cytoskeleton and cilia (Supplementary Table S11). For instance, the gene TNS3, related to cytoskeleton, is predicted as target of 16 up-regulated miRNAs (Supplementary Table S11). In addition, several genes encoding for dyneins were predicted as targets of up-regulated miRNAs, with DNAL1 being target of 7 different miRNAs. Mutations and down-regulation of genes coding ciliary proteins, especially dyneins, were related to primary ciliary dyskinesia46, 48, 53, 55, 91, but had not yet been reported in M. hyopneumoniae infection. More investigations on this matter should be carried out, since these results seem to be related to one of the main adverse effects caused by this bacterium: ciliostasis.
Most DE miRNAs were detected in exosome-like vesicles
Analyzing the up-regulated genes, we identified that besides antioxidant response and protein synthesis, there is an enrichment in GO terms related to cell-cell communication (Fig. 1). The up-regulation of cell-cell communication is noteworthy considering that most DE miRNAs were detected only in the extracellular and exosome-like samples (Supplementary Table S10). Moreover, most of the down-regulated miRNAs of exosome-like vesicles were DE only in these vesicles, indicating that the cells might be selecting specific populations of miRNAs to be packaged into these structures, which could be seen as a specific message they are trying to communicate to other cells. In the same direction, the increase in the number of proteins secreted in vesicles was reported in NPTr cells infected with M. hyopneumoniae92, indicating that the communication via exosome-like vesicles is important for swine cells during infection. This corroborates our data, suggesting that M. hyopneumoniae induces modification on the protein and RNA composition of NPTr-released vesicles and is likely important in the cross-talks between epithelial cells.
The relationship between cell–cell communication and DE miRNAs released by infected cells becomes especially interesting when considering genes related to the NRF2 pathway: we identified DE miRNAs predicted to regulate such genes only in extracellular and exosome-like vesicle samples, with no difference of expression in intracellular samples (Table 6). This might suggest the existence of a mechanism in which infected cells send signals to neighboring cells in order to prevent the repression of these genes or degradation of their RNA products. This becomes even more important due to the fact that exosomal miRNAs were already shown to regulate the inflammatory response in receptor cells from mice93. While miRNAs identified in the extracellular total samples might contain sRNAs that reflect degradation remnants due to the presence of RNAses and mRNAs in the extracellular environment, we believe that this is less likely to happen in exosomes, since these vesicles have a membrane that protects their content from extracellular degradation. As exosomes have an important role in cell communication, we should consider the difference in miRNA expression in these vesicles as relevant and speculate on how they might interfere in the gene regulation of neighboring cells. In addition, as exosomes may also affect other cell types, they may contain miRNAs that target genes that were not detected as DE in our samples.
Concluding remarks
M. hyopneumoniae is considered a pathogen with huge negative impact on swine production. However, apart from studies related to adherence to the host cells, little is known about the relationship between this pathogen and the swine host. In this work, we analyzed the changes that M. hyopneumoniae induced in gene and miRNA expression in tracheal epithelial cells. As far as we know, this is the first study to establish a link between gene expression of the swine cells and the most deleterious pathogenic effects of M. hyopneumoniae, namely its cytotoxic epithelial damage (possibly via hydrogen peroxide production) and induced ciliostasis. However, we are aware that infection with this species involves a large component of the immune system, which is even said to be responsible for the major tissue damage related to the infection. Thus, the results found here reflect the effects of M. hyopneumoniae on epithelial cells, while the general picture of the respiratory tract might present different responses to the infection. In this way, it is important to observe that besides the hydrogen peroxide production by this bacterium, other factors can trigger the activation of the cell antioxidant defense during infection. M. hyopneumoniae infection is characterized by the infiltration of a large number of leukocytes in the lung tissue, which produce reactive oxygen and nitrogen species, causing tissue damage38, 94. Accordingly, systemic infections with M. hyopneumoniae have a great potential of causing oxidative stress, so that the transcription of genes activated by NRF2 seems to be important to fight the infection and maintain cellular homeostasis.
In conclusion, we conducted the first study that analyzes mRNA and miRNA differential expression by NGS in epithelial tracheal cells infected with M. hyopneumoniae. Our results bring new insights into the interaction between this bacterium and swine epithelial cells, notably the host cellular response through the activation of genes related to antioxidant response and the repression of cytoskeleton and ciliary genes (possibly related to ciliostasis), and open several perspectives related to the understanding of the pathogenicity of this bacterial species.
Methods
Cultivation of NPTr cells
NPTr cells95 negatively tested for mycoplasma, purchased from Instituto Zooprofilattico Sperimentale della Lombardia e dell’Emilia Romagna, Brescia, Italy, were grown in Minimum Essential Medium with Earle’s Balanced Salts (MEM/EBSS) and 2 mM L-Glutamine (GE Healthcare) supplemented with 10% (v/v) heat-inactivated fetal bovine serum (FBS) (Gibco) and 1% penicillin streptomycin (GE Healthcare). Cells were cultured at \(37\,{^{\circ }}\hbox {C} \) with 5% CO2 in a humid atmosphere in T75 \(cm^{2}\) flasks and 6-well tissue culture plates. Sub-passages were made when cells reached 70-80% confluence. To isolate exosome-like vesicles, the culture medium was removed, confluent cells were washed with phosphate-buffered saline (PBS) and kept in medium without FBS and antibiotics (exosome depleted) for 24 h.
Infection of NPTr cells with M. hyopneumoniae
Mycoplasma hyopneumoniae strain J (ATCC25934) was cultivated in Friis medium96 at \(37\,{^{\circ }}\hbox {C}\) for 48 h with gentle agitation in a roller drum. Prior to NPTr infection, Mycoplasma cells were pelleted, washed with PBS and resuspended in MEM/EBSS medium. Based on the methodology described by Assunção et al. (2005)97, we calculated an MOI of 5 M. hyopneumoniae for the immunofluorescence experiments and an MOI of 2 for the RNA extraction experiments. NPTr cells were challenged with M. hyopneumoniae for 24 h.
Immunofluorescence microscopy
NPTr cells were grown in cover slips coated with poly-L-lysine (Sigma-Aldrich) for 48 h and infected (or not, in the control conditions) with M. hyopneumoniae for 24 h or 1 h in the case of the NFR2 localization experiment. The cells were washed with 1x PBS three times and fixed with 4% paraformaldehyde for 15 min. After permeabilization with 0.2% Triton X-100 for 10 min and blocking with 4% BSA in PBS for 20 min, the cells were incubated with the primary antibodies (rabbit polyclonal anti-SPAse I of M. hyopneumoniae98, 1:10,000, and mouse monoclonal anti-Sodium/Potassium ATPase alpha (Invitrogen), 1:20, or rabbit polyclonal anti-NRF2 (Thermo Fisher Scientific), 1:200,) overnight at \(4\,{^{\circ }}\hbox {C}\) followed by secondary antibodies (Alexa Fluor488-conjugated anti-rabbit (Invitrogen), 1:1,000, and Alexa Fluor555-conjugated anti-mouse (Invitrogen), 1:1,000) incubation for 1 h in the dark. The cells were washed with PBS and the DNA was stained with 100 nM DAPI and actin was marked with 50 nM phalloidin for 20 min. The slides were mounted with Fluoromount (Sigma-Aldrich). Samples were imaged using an Olympus FluoView 1000 confocal microscope. NRF2 fluorescence was quantified by drawing margin to each individual nucleus and measuring the green fluorescence against the cell background. At least 25 nuclei were quantified in each condition with the ImageJ software99. The Corrected Total Fluorescence (CTCF) of each nucleus was calculated with the following formula: CTCF = Integrated density – (Area of the selected nucleus x fluorescence of the background reading). Statistical significance was analyzed by One-way ANOVA followed by Dunnett’s multiple comparison test (\(\hbox {p} < 0.05\)) in GraphPad Prism 7.0 software.
Exosome extraction
Exosome-like vesicles were purified from 250 mL (25 F75 \(cm^{2}\) flasks) of NPTr conditioned medium (as described by Forterre et al.17). Briefly, after 24 hours of incubation, the medium was recovered and cell debris and organelles were eliminated by centrifugation at 2,000 x g for 20 min and at 10,000 x g for 30 min. The resulting supernatant was filtered through a \(0.22\,\mu \hbox {m}\) filter in order to remove large particles. Exosome-like vesicles were pelleted by ultracentrifugation at 100,000 x g for 90 min \(4\,{^{\circ }}\hbox {C}\) (Beckman-Coulter, Optima MAX-XP ultracentrifuge, MLA-55 rotor). The supernatant was recovered and saved for posterior use (SN) and the exosome sediment was washed with 25 mL of cold PBS, centrifuged at 100,000 x g for 70 min \(4 \,{^{\circ }}\hbox {C}\) and resuspended in \(50\,\mu \hbox {L}\) PBS. Vesicle extractions were performed in duplicates in both conditions of NPTr cells non-infected and infected with M. hyopneumoniae. SN portion was submitted to nucleic acid precipitation with 1/10 volume of 3 M sodium acetate and 3 volumes of 100% cold ethanol at \(-20\,{^{\circ }}\hbox {C}\) overnight. After centrifugation, the nucleic acid precipitated was allowed to dry at \(60\,{^{\circ }}\hbox {C}\) for 1 hour and was ressuspended in \(100\,\mu \hbox {L}\) of ultrapure RNAse free water.
Sample preparation, RNA extraction and sequencing
For mRNA sequencing, a total of 6 samples were prepared: 3 in a control group (CTL) and 3 in the infected group (INF). For miRNA sequencing, we prepared the samples as follows: total intracellular miRNA (INTRA: 3 CTL vs 3 INF), total extracellular small RNA (EXTRA: 2 CTL vs 2 INF), extracellular exosomal miRNA (EXO: 1 pool CTL vs 1 pool INF) and extracellular miRNA from vesicle-free supernatant (SN: 1 pool CTL vs 1 pool INF: a single library was constructed from a pool of 50 biological replicates). Total RNA for mRNA sequencing and sRNA enriched (\(< 200\)bp) for sRNA sequencing and miRNA analyses were extracted with mirVana kit (Ambion), according to the manufacturer’s instructions. Total extracellular RNA for sRNA sequencing was directly extracted from the culture NPTr cell supernatant after centrifugation of cell debris. Exosome sRNA was extracted after vesicle purification. RNA quality was assessed with Bioanalyzer 2100 (Agilent Genomics). A total of 6 mRNA libraries were prepared (one for each sample) using TruSeq Stranded Total RNA Sample Preparation kit (Illumina) and a total of 14 small RNA libraries were prepared using TruSeq Small RNA Sample Prep kit (Illumina). After quality control, the sequencing of all libraries was performed by HiSeq2500 platform (Illumina). RNA library preparation and sequencing was performed by the Duke University GCB Sequencing Platform (Durham, USA).
RT-qPCR of selected differentially expressed mRNAs
For RT-qPCR, total RNA of NPTr cells infected and non-infected with M. hyopneumoniae was isolated with TRizol (Invitrogen) according to the manufacturer’s instructions. A first-strand cDNA synthesis reaction was conducted by adding \(1\,\mu \hbox {g}\) of total RNA to 500 ng of oligo(dT)15 primer (Promega) and 10 mM deoxynucleotide triphosphates. The mixture was heated for \(65\,{^{\circ }}\hbox {C}\) for 5 min and then incubated on ice for 5 min. First-strand buffer (Invitrogen), 0.1 M dithiothreitol and 200 U M-MLV RT (Moloney Murine Leukemia Virus Reverse Transcriptase – Invitrogen) were then added to a total volume of \(20\,\mu \hbox {L}\). The reaction was incubated at \(25\,{^{\circ }}\hbox {C}\) for 10 min and at \(37\,{^{\circ }}\hbox {C}\) for 50 min followed by 15 min at \(70\,{^{\circ }}\hbox {C}\) for enzyme inactivation. Quantitative PCR (qPCR) assay was performed using 1:5 cDNA as template and Platinum SYBR Green qPCR SuperMix-UDG (Invitrogen) with specific primers (Supplementary Table S13) on 7500 Real-Time PCR Systems (Applied Biosystems). The qPCR reactions were carried out at \(90\,{^{\circ }}\hbox {C} \) for 2 min and \(95\,{^{\circ }}\hbox {C}\) for 10 min followed by 40 cycles of \(95\,{^{\circ }}\hbox {C}\) for 15 s and \(60\,{^{\circ }}\hbox {C}\) for 1 min. A melting curve analysis was done to verify the specificity of the synthesized products and the amplification efficiency for each primer was calculated with the LinRegPCR software application100. The mRNA levels were normalized against porcine DNA topoisomerase II beta (TOP2B) and the relative expression of mRNA was calculated by \(2^{-\Delta \Delta {\mathrm{CT}}}\) method. The threshold cycle (CT) of each target test represents the average of three reactions. Three independent biological replicates were done for each target gene. Statistical analyses were performed using the GraphPad Prism 7.0 software with a two-tailed unpaired t-test to compare the relative expression between infected and non-infected NPTr cells (\(\hbox {p} < 0.05\)).
Computational analysis of sequencing data
Preprocessing of raw reads
The reads from sRNA-seq and total RNA-seq were processed separately. The raw reads were filtered for low quality, and adapter sequences were trimmed with CUTADAPT101. Total mRNA-seq clean reads were mapped against the porcine reference genome from Ensembl (Sscrofa10.2) with STAR102 and the raw counts of genes were obtained with HTSeq103. Only uniquely mapped reads were used for further analyses.
For miRNA prediction, clean sRNA-seq reads were mapped against the porcine genome with Bowtie104. We used intracellular samples as input for miRDeep276 and kept predictions that had a score of at least 5. After this, we collapsed similar predictions and obtained a total of 773 clusters (773 miRNAs), of which 478 were novel miRNAs. We created a porcine miRNA DB (ssc-miRNA-DB), by including the 411 annotated porcine miRNAs in version 21 of miRBase77, 105 along with the 722 miRNAs characterized by Martini et al. (2014)78 and our novel predictions. Reads from all samples were mapped against this database and a matrix of counts was generated in order to identify DE miRNAs.
Differential expression of mRNA data
Raw counts were used as input for DE gene analysis. We detected possible DE genes in both EdgeR106 and DESeq2107 packages in R. DESeq2 was run for genes that had a total count of at least 10 in all libraries, with the method’s default normalization. EdgeR was used with TMM normalization and general linear model fit, only for genes with cpm greater than 1 in at least 2 libraries. After testing, the adjusted p-values (p-adj) for both methods were adjusted with the Benjamini-Hochberg108 correction for multitesting. Genes were considered DE when p-adj \(< 0.05\) and genes with the most pronounced logFC were selected individually for further investigation. Overall, DESeq2 detected more significant p-adj with less accentuated LogFC, while EdgeR detected less significant false discovery rates (FDR) with more extreme LogFCs. Both techniques have been widely used separately and together in several publications, and in order to select good candidates for testing, we took into account the results of both methods. For GO functional analyses, we also used the complete list of DE genes whenever LogFC was greater than 0.1 (up-regulated) or smaller than \(-0.1\) (down-regulated). The complete lists from each method were used separately, and we compared the overall outcomes to check the robustness of our results.
Differential expression of miRNA data
The same pipeline used for DE mRNA was performed for total intracellular and total extracellular miRNAs, since we had biological replicates in both cases. In the case where we had no replicates (instead, a single library was constructed from a pool of 50 biological replicates), we used GFOLD109 which provides a generalized fold change for ranking DE genes. GFOLD is said to overcome the shortcomings of p-value and fold change of the existing methods and can provide a more stable and biological meaningful gene ranking when a single biological replicate is available. In this case, we selected miRNAs with GFOLD \(> 2\) or \(< -2\) for functional analyses.
miRNA target prediction
DE miRNAs were used as input to detect putative interactions with the UTRs of Ensembl transcripts in the porcine genome. We used three methods to detect target pairs: miRanda110, TargetScan111 and PITA112, and one method to validate the hybridization of a target pair, RNAhybrid113. We kept only targets that were predicted by at least two distinct tools and used the following thresholds: score in miRanda \(> 140\), DDG from PITA \(< -5\), score in RNAhybrid \(< -15\). Since TargetScan does not provide a continuous scoring system, we only validated whenever there was a prediction of at least 6mers. Based on these, we chose from the target list only genes that were detected as DE in this study, and subsequentially we only considered target pairs of miRNA-mRNA that had inversed LogFC expression.
Functional analysis of DE mRNAs and targets of DE miRNAs
Functional analysis took into account as input either the list of DE mRNAs itself or the list of targets predicted for the DE miRNAs. We performed a GO enrichment analysis114 to find out which functions were over or underrepresented in each gene list. P-values for enriched GO terms were adjusted with the Benjamini-Hochberg108 correction for multitesting. GO terms and pathways with p-adj \(< 0.05\) were defined as significantly enriched. The GO terms were reduced to representative non redundant terms with the use of the REVIGO tool115.
Regulatory network reconstruction and analysis
We created a general regulatory network of the host response to the bacterial infection with the DE miRNAs and target mRNAs detected in this study (Supplementary Fig. S9 and Supplementary File S1). In this network, we also included information about interactions from the BioGRID v3.4 database80, 81, a general repository that includes experimentally validated physical and genetic interactions. We used human-based official gene symbols to include information from BioGRID, Cytoscape116 to draw the networks, and the ClueGO plugin82 to perform functional enrichment analysis. We further manually curated the networks for genes and miRNAs related with redox homeostasis (permissive regulatory network), as well as cytoskeleton and cilia (repressive regulatory network).
Accession codes
All sequences were submitted to the Sequence Read Archive from NCBI: SRA accession: PRJNA545822.
References
Maes, D., Verdonck, M., Deluyker, H. & de Kruif, A. Enzootic pneumonia in pigs. Vet. Q.18, 104–109 (1996).
Maes, D. et al. Update on Mycoplasma hyopneumoniae infections in pigs: knowledge gaps for improved disease control. Transbound. Emerg. Dis.65, 110–124 (2018).
Djordjevic, S. P., Cordwell, S. J., Djordjevic, M. A., Wilton, J. & Minion, F. C. Proteolytic processing of the Mycoplasma hyopneumoniae cilium adhesin. Infect. Immun.72, 2791–2802 (2004).
DeBey, M. C. & Ross, R. F. Ciliostasis and loss of cilia induced by Mycoplasma hyopneumoniae in porcine tracheal organ cultures. Infect. Immun.62, 5312–5318 (1994).
Ferrarini, M. G. et al. Hydrogen peroxide production and myo-inositol metabolism as important traits for virulence of Mycoplasma hyopneumoniae. Mol. Microbiol.108, 683–696 (2018).
Bartel, D. P. Micrornas: genomics, biogenesis, mechanism, and function. Cell116, 281–297 (2004).
Luo, Y. et al. Detection of dietetically absorbed maize-derived micrornas in pigs. Sci. Rep.7, 645 (2017).
Lawless, N., Foroushani, A. B., McCabe, M. S., O’Farrelly, C. & Lynn, D. J. Next generation sequencing reveals the expression of a unique mirna profile in response to a gram-positive bacterial infection. PLoS ONE8, e57543 (2013).
Duval, M., Cossart, P. & Lebreton, A. Mammalian micrornas and long noncoding rnas in the host-bacterial pathogen crosstalk. In Seminars in cell & developmental biology, vol. 65, 11–19 (Elsevier, 2017).
Lim, L. P. et al. Microarray analysis shows that some micrornas downregulate large numbers of target mrnas. Nature433, 769 (2005).
Kim, V. N. & Nam, J.-W. Genomics of microrna. Trends Genet.22, 165–173 (2006).
Valadi, H. et al. Exosome-mediated transfer of mrnas and micrornas is a novel mechanism of genetic exchange between cells. Nat. Cell Biol.9, 654 (2007).
Rome, S. Are extracellular micrornas involved in type 2 diabetes and related pathologies?. Clin. Biochem.46, 937–945 (2013).
Chen, L. et al. Epigenetic regulation of connective tissue growth factor by microrna-214 delivery in exosomes from mouse or human hepatic stellate cells. Hepatology59, 1118–1129 (2014).
Schorey, J. S., Cheng, Y., Singh, P. P. & Smith, V. L. Exosomes and other extracellular vesicles in host–pathogen interactions. EMBO Rep.16, 24–43 (2015).
Rana, S., Yue, S., Stadel, D. & Zöller, M. Toward tailored exosomes: the exosomal tetraspanin web contributes to target cell selection. Int. J. Biochem. Cell Biol.44, 1574–1584 (2012).
Forterre, A. et al. Myotube-derived exosomal mirnas downregulate sirtuin1 in myoblasts during muscle cell differentiation. Cell Cycle13, 78–89 (2014).
Silverman, J. M. & Reiner, N. E. Exosomes and other microvesicles in infection biology: organelles with unanticipated phenotypes. Cell. Microbiol.13, 1–9 (2011).
Muxel, S. M., Laranjeira-Silva, M. F., Zampieri, R. A. & Floeter-Winter, L. M. Leishmania (Leishmania) amazonensis induces macrophage mir-294 and mir-721 expression and modulates infection by targeting nos2 and l-arginine metabolism. Sci. Rep.7, 44141 (2017).
Bao, H. et al. Genome-wide whole blood micrornaome and transcriptome analyses reveal mirna-mrna regulated host response to foodborne pathogen Salmonella infection in swine. Sci. Rep.5, 12620 (2015).
Staedel, C. & Darfeuille, F. Micro rna s and bacterial infection. Cell. Microbiol.15, 1496–1507 (2013).
Maudet, C., Mano, M. & Eulalio, A. Micrornas in the interaction between host and bacterial pathogens. FEBS Lett.588, 4140–4147 (2014).
Zielinski, G. & Ross, R. Effect of growth in cell cultures and strain on virulence of Mycoplasma hyopneumoniae for swine. Am. J. Vet. Res.51, 344–348 (1990).
Young, T. F., Thacker, E. L., Erickson, B. Z. & Ross, R. F. A tissue culture system to study respiratory ciliary epithelial adherence of selected swine mycoplasmas. Vet. Microbiol.71, 269–279 (2000).
Burnett, T. A. et al. P159 is a proteolytically processed, surface adhesin of Mycoplasma hyopneumoniae: defined domains of p159 bind heparin and promote adherence to eukaryote cells. Mol. Microbiol.60, 669–686 (2006).
Cho, H.-Y., Reddy, S. P. & Kleeberger, S. R. Nrf2 defends the lung from oxidative stress. Antioxid. Redox Signal.8, 76–87 (2006).
Ma, Q. Role of nrf2 in oxidative stress and toxicity. Annu. Rev. Pharmacol. Toxicol.53, 401–426 (2013).
Thimmulappa, R. K. et al. Nrf2 is a critical regulator of the innate immune response and survival during experimental sepsis. J. Clin. Investig.116, 984–995 (2016).
Athale, J. et al. Nrf2 promotes alveolar mitochondrial biogenesis and resolution of lung injury in staphylococcus aureus pneumonia in mice. Free Radical Biol. Med.53, 1584–1594 (2012).
Reddy, N. M. et al. Innate immunity against bacterial infection following hyperoxia exposure is impaired in nrf2-deficient mice. J. Immunol.183, 4601–4608 (2009).
Gomez, J. C., Dang, H., Martin, J. R. & Doerschuk, C. M. Nrf2 modulates host defense during Streptococcus pneumoniae pneumonia in mice. J. Immunol.197, 2864–2879 (2016).
Kensler, T. W., Wakabayashi, N. & Biswal, S. Cell survival responses to environmental stresses via the keap1-nrf2-are pathway. Annu. Rev. Pharmacol. Toxicol.47, 89–116 (2007).
Fourquet, S., Guerois, R., Biard, D. & Toledano, M. B. Activation of nrf2 by nitrosative agents and h2o2 involves keap1 disulfide formation. J. Biol. Chem.285, 8463–8471 (2010).
Hames, C., Halbedel, S., Hoppert, M., Frey, J. & Stülke, J. Glycerol metabolism is important for cytotoxicity of Mycoplasma pneumoniae. J. Bacteriol.191, 747–753 (2009).
Vilei, E. M. & Frey, J. Genetic and biochemical characterization of glycerol uptake in Mycoplasma mycoides subsp. mycoides sc: Its impact on h2o2production and virulence. Clin. Diagn. Lab. Immunol.8, 85–92 (2001).
Bischof, D. F., Janis, C., Vilei, E. M., Bertoni, G. & Frey, J. Cytotoxicity of mycoplasma mycoides subsp. mycoides small colony type to bovine epithelial cells. Infect. Immun.76, 263–269 (2008).
Halbedel, S., Hames, C. & Stülke, J. In vivo activity of enzymatic and regulatory components of the phosphoenolpyruvate: sugar phosphotransferase system in mycoplasma pneumoniae. J. Bacteriol.186, 7936–7943 (2004).
Biswal, S., Thimmulappa, R. K. & Harvey, C. J. Experimental therapeutics of nrf2 as a target for prevention of bacterial exacerbations in copd. Proc. Am. Thorac. Soc.9, 47–51 (2012).
Nioi, P., McMahon, M., Itoh, K., Yamamoto, M. & Hayes, J. D. Identification of a novel nrf2-regulated antioxidant response element (are) in the mouse nad (p) h: quinone oxidoreductase 1 gene: reassessment of the are consensus sequence. Biochem. J.374, 337–348 (2003).
Rice, P., Longden, I. & Bleasby, A. Emboss: the European molecular biology open software suite (2000).
Araake, M. Comparison of ciliostasis by mycoplasmas in mouse and chicken tracheal organ cultures. Microbiol. Immunol.26, 1–14 (1982).
Abu-Zahr, M. & Butler, M. Growth, cytopathogenicity and morphology of Mycoplasma gallisepticum and M. gallinarum in tracheal explants. J. Comp. Pathol.86, 455–463 (1976).
Chaudhry, R. et al. Pathogenesis of Mycoplasma pneumoniae: an update. Indian J. Med. Microbiol.34, 7 (2016).
Matsuyama, M. et al. Transcriptional response of respiratory epithelium to nontuberculous mycobacteria. Am. J. Respir. Cell Mol. Biol.58, 241–252 (2018).
Kobayashi, D. & Takeda, H. Ciliary motility: the components and cytoplasmic preassembly mechanisms of the axonemal dyneins. Differentiation83, S23–S29 (2012).
Tarkar, A. et al. Dyx1c1 is required for axonemal dynein assembly and ciliary motility. Nat. Genet.45, 995 (2013).
Wang, W.-J. et al. Cep162 is an axoneme-recognition protein promoting ciliary transition zone assembly at the cilia base. Nat. Cell Biol.15, 591 (2013).
Schueler, M. et al. Dcdc2 mutations cause a renal-hepatic ciliopathy by disrupting wnt signaling. Am. J. Human Genet.96, 81–92 (2015).
May-Simera, H. L. et al. Loss of macf1 abolishes ciliogenesis and disrupts apicobasal polarity establishment in the retina. Cell Rep.17, 1399–1413 (2016).
Shi, L., Koll, F., Arnaiz, O. & Cohen, J. The ciliary protein ift 57 in the macronucleus of paramecium. J. Eukaryot. Microbiol.65, 12–27 (2018).
Yasunaga, T. et al. The polarity protein inturned links nphp4 to daam1 to control the subapical actin network in multiciliated cells. J. Cell Biol.211, 963–973 (2015).
Hegan, P. S., Ostertag, E., Geurts, A. M. & Mooseker, M. S. Myosin id is required for planar cell polarity in ciliated tracheal and ependymal epithelial cells. Cytoskeleton72, 503–516 (2015).
Fernandez-Gonzalez, A., Kourembanas, S., Wyatt, T. A. & Mitsialis, S. A. Mutation of murine adenylate kinase 7 underlies a primary ciliary dyskinesia phenotype. Am. J. Respir. Cell Mol. Biol.40, 305–313 (2009).
Zhou, J., Yang, F., Leu, N. A. & Wang, P. J. Mns1 is essential for spermiogenesis and motile ciliary functions in mice. PLoS Genet.8, e1002516 (2012).
Inaba, Y. et al. Transport of the outer dynein arm complex to cilia requires a cytoplasmic protein lrrc6. Genes Cells21, 728–739 (2016).
Arachchige, S. N. K. et al. Differential response of the chicken trachea to chronic infection with virulent Mycoplasma gallisepticum strain ap3as and vaxsafe mg (strain ts-304): a transcriptional profile. Infect. Immun.88, (2020).
Girón, J. A., Lange, M. & Baseman, J. B. Adherence, fibronectin binding, and induction of cytoskeleton reorganization in cultured human cells by Mycoplasma penetrans. Infect. Immun.64, 197–208 (1996).
Raymond, B. et al. Mycoplasma hyopneumoniae resides intracellularly within porcine epithelial cells. Sci. Rep.8, 17697 (2018).
Mirvis, M., Stearns, T. & Nelson, W. J. Cilium structure, assembly, and disassembly regulated by the cytoskeleton. Biochem. J.475, 2329–2353 (2018).
Raymond, B. et al. Extracellular actin is a receptor for Mycoplasma hyopneumoniae. Front. Cell. Infect. Microbiol.8, 54 (2018).
Kang, M.-I., Kobayashi, A., Wakabayashi, N., Kim, S.-G. & Yamamoto, M. Scaffolding of keap1 to the actin cytoskeleton controls the function of nrf2 as key regulator of cytoprotective phase 2 genes. Proc. Nat. Acad. Sci.101, 2046–2051 (2004).
Mack, G. J. & Compton, D. A. Analysis of mitotic microtubule-associated proteins using mass spectrometry identifies astrin, a spindle-associated protein. Proc. Nat. Acad. Sci.98, 14434–14439 (2001).
Kim, Y. H., Choi, J. Y., Jeong, Y., Wolgemuth, D. J. & Rhee, K. Nek2 localizes to multiple sites in mitotic cells, suggesting its involvement in multiple cellular functions during the cell cycle. Biochem. Biophys. Res. Commun.290, 730–736 (2002).
Woolner, S., O’Brien, L. L., Wiese, C. & Bement, W. M. Myosin-10 and actin filaments are essential for mitotic spindle function. J. Cell Biol.182, 77–88 (2008).
Bolanos-Garcia, V. M. & Blundell, T. L. Bub1 and bubr1: multifaceted kinases of the cell cycle. Trends Biochem. Sci.36, 141–150 (2011).
Wolter, P. et al. Gas2l3, a target gene of the dream complex, is required for proper cytokinesis and genomic stability. J. Cell Sci.125, 2393–2406 (2012).
Xu, X.-L. et al. The microtubule-associated protein aspm regulates spindle assembly and meiotic progression in mouse oocytes. PLoS ONE7, e49303 (2012).
Matson, D. R. & Stukenberg, P. T. Cenp-i and aurora b act as a molecular switch that ties rzz/mad1 recruitment to kinetochore attachment status. J. Cell Biol.205, 541–554 (2014).
Wadsworth, P. Tpx2. Curr. Biol.25, R1156–R1158 (2015).
Cho, R. J. et al. Transcriptional regulation and function during the human cell cycle. Nat. Genet.27, 48 (2001).
Shirin, H. et al. Helicobacter pylori inhibits the g1 to s transition in ags gastric epithelial cells. Cancer Res.59, 2277–2281 (1999).
Marchès, O. et al. Enteropathogenic and enterohaemorrhagic Escherichia coli deliver a novel effector called cif, which blocks cell cycle g2/m transition. Mol. Microbiol.50, 1553–1567 (2003).
Jones, A., Jonsson, A.-B. & Aro, H. Neisseria gonorrhoeae infection causes a g1 arrest in human epithelial cells. FASEB J.21, 345–355 (2007).
Alekseeva, L. et al. Staphylococcus aureus-induced g2/m phase transition delay in host epithelial cells increases bacterial infective efficiency. PLoS ONE8, e63279 (2013).
Siqueira, F. M. et al. Mycoplasma non-coding rna: identification of small rnas and targets. BMC Genom.17, 743 (2016).
Friedländer, M. R., Mackowiak, S. D., Li, N., Chen, W. & Rajewsky, N. mirdeep2 accurately identifies known and hundreds of novel microrna genes in seven animal clades. Nucl. Acids Res.40, 37–52 (2011).
Kozomara, A. & Griffiths-Jones, S. mirbase: annotating high confidence micrornas using deep sequencing data. Nucl. Acids Res.42, D68–D73 (2013).
Martini, P. et al. Tissue-specific expression and regulatory networks of pig micrornaome. PLoS ONE9, e89755 (2014).
Bartel, D. P. Metazoan micrornas. Cell173, 20–51 (2018).
Stark, C. et al. Biogrid: a general repository for interaction datasets. Nucl. Acids Res.34, D535–D539 (2006).
Chatr-Aryamontri, A. et al. The biogrid interaction database: 2017 update. Nucl. Acids Res.45, D369–D379 (2017).
Bindea, G. et al. Cluego: a cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics25, 1091–1093 (2009).
Yang, M., Yao, Y., Eades, G., Zhang, Y. & Zhou, Q. Mir-28 regulates nrf2 expression through a keap1-independent mechanism. Breast Cancer Res. Treat.129, 983–991 (2011).
Shi, L. et al. mir-340 reverses cisplatin resistance of hepatocellular carcinoma cell lines by targeting nrf2-dependent antioxidant pathway. Asian Pac. J. Cancer Prev.15, 10439–10444 (2015).
Xu, J. et al. mir-19b attenuates h2o2-induced apoptosis in rat h9c2 cardiomyocytes via targeting pten. Oncotarget7, 10870 (2016).
Hong, J. et al. Transcriptional downregulation of microrna-19a by ros production and nf-\(\kappa \)b deactivation governs resistance to oxidative stress-initiated apoptosis. Oncotarget8, 70967 (2017).
Gao, A.-M., Zhang, X.-Y. & Ke, Z.-P. Apigenin sensitizes bel-7402/adm cells to doxorubicin through inhibiting mir-101/nrf2 pathway. Oncotarget8, 82085 (2017).
Liu, H. et al. Regulation of mir-92a on vascular endothelial aging via mediating nrf2-keap1-are signal pathway. Eur. Rev. Med. Pharmacol. Sci.21, 2734–2742 (2017).
Yang, H. et al. Activation of a novel c-myc-mir27-prohibitin 1 circuitry in cholestatic liver injury inhibits glutathione synthesis in mice. Antioxid. Redox Signal.22, 259–274 (2015).
Alural, B., Ozerdem, A., Allmer, J., Genc, K. & Genc, S. Lithium protects against paraquat neurotoxicity by nrf2 activation and mir-34a inhibition in sh-sy5y cells. Front. Cell. Neurosci.9, 209 (2015).
Milara, J., Armengot, M., Mata, M., Morcillo, E. J. & Cortijo, J. Role of adenylate kinase type 7 expression on cilia motility: possible link in primary ciliary dyskinesia. Am. J. Rhinol. Allergy24, 181–185 (2010).
Zimmer, FMd. A. L., Paludo, G. P., Moura, H., Barr, J. R. & Ferreira, H. B. Differential secretome profiling of a swine tracheal cell line infected with mycoplasmas of the swine respiratory tract. J. Proteom.192, 147–159 (2019).
Alexander, M. et al. Exosome-delivered micrornas modulate the inflammatory response to endotoxin. Nat. Commun.6, 7321 (2015).
Muneta, Y. et al. Immune response of gnotobiotic piglets against Mycoplasma hyopneumoniae. J. Vet. Med. Sci.70, 1065–1070 (2008).
Ferrari, M. et al. Establishment and characterization of two new pig cell lines for use in virological diagnostic laboratories. J. Virol. Methods107, 205–212 (2003).
Friis, N. Some recommendations concerning primary isolation of Mycoplasma suipneumoniae and Mycoplasma flocculare a survey. Nordisk veterinaermedicin27, 337–339 (1975).
Assunção, P. et al. Evaluation of Mycoplasma hyopneumoniae growth by flow cytometry. J. Appl. Microbiol.98, 1047–1054 (2005).
Moitinho-Silva, L. et al.Mycoplasma hyopneumoniae type i signal peptidase: expression and evaluation of its diagnostic potential. Vet. Microbiol.154, 282–291 (2012).
Schneider, C. A., Rasband, W. S. & Eliceiri, K. W. Nih image to imagej: 25 years of image analysis. Nat. Methods9, 671 (2012).
Ruijter, J. et al. Amplification efficiency: linking baseline and bias in the analysis of quantitative pcr data. Nucl. Acids Res.37, e45–e45 (2009).
Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J.17, 10–12 (2011).
Dobin, A. et al. Star: ultrafast universal rna-seq aligner. Bioinformatics29, 15–21 (2013).
Anders, S., Pyl, P. T. & Huber, W. Htseq—a python framework to work with high-throughput sequencing data. Bioinformatics31, 166–169 (2015).
Langmead, B., Trapnell, C., Pop, M. & Salzberg, S. L. Ultrafast and memory-efficient alignment of short dna sequences to the human genome. Genome Biol.10, R25 (2009).
Griffiths-Jones, S. The microrna registry. Nucl. Acids Res.32, D109–D111 (2004).
Robinson, M. D., McCarthy, D. J. & Smyth, G. K. edger: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics26, 139–140 (2010).
Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for rna-seq data with deseq2. Genome Biol.15, 550 (2014).
Benjamini, Y. & Hochberg, Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc.: Ser. B (Methodol.)57, 289–300 (1995).
Feng, J. et al. Gfold: a generalized fold change for ranking differentially expressed genes from rna-seq data. Bioinformatics28, 2782–2788 (2012).
Betel, D., Koppal, A., Agius, P., Sander, C. & Leslie, C. Comprehensive modeling of microrna targets predicts functional non-conserved and non-canonical sites. Genome Biol.11, R90 (2010).
Lewis, B. P., Shih, I.-H., Jones-Rhoades, M. W., Bartel, D. P. & Burge, C. B. Prediction of mammalian microrna targets. Cell115, 787–798 (2003).
Kertesz, M., Iovino, N., Unnerstall, U., Gaul, U. & Segal, E. The role of site accessibility in microrna target recognition. Nat. Genet.39, 1278 (2007).
Krüger, J. & Rehmsmeier, M. Rnahybrid: microrna target prediction easy, fast and flexible. Nucl. Acids Res.34, W451–W454 (2006).
Ashburner, M. et al. Gene ontology: tool for the unification of biology. Nat. Genet.25, 25 (2000).
Supek, F., Bošnjak, M., Škunca, N. & Šmuc, T. Revigo summarizes and visualizes long lists of gene ontology terms. PLoS ONE6, e21800 (2011).
Shannon, P. et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res.13, 2498–2504 (2003).
Hayes, J. D. & McMahon, M. Nrf2 and keap1 mutations: permanent activation of an adaptive response in cancer. Trends Biochem. Sci.34, 176–188 (2009).
Hayes, J. D. & Dinkova-Kostova, A. T. The nrf2 regulatory network provides an interface between redox and intermediary metabolism. Trends Biochem. Sci.39, 199–218 (2014).
Brzóska, K., Stępkowski, T. M. & Kruszewski, M. Basal pir expression in hela cells is driven by nrf2 via evolutionary conserved antioxidant response element. Mol. Cell. Biochem.389, 99–111 (2014).
Li, D., Ma, S. & Ellis, E. M. Nrf2-mediated adaptive response to methyl glyoxal in hepg2 cells involves the induction of akr7a2. Chem. Biol. Interact.234, 366–371 (2015).
Graham, D. B. et al. Nitric oxide engages an anti-inflammatory feedback loop mediated by peroxiredoxin 5 in phagocytes. Cell Rep.24, 838–850 (2018).
Reszka, E. et al. Expression of nrf2 and nrf2-modulated genes in peripheral blood leukocytes of bladder cancer males. Neoplasma60, 123–128 (2013).
Dai, L.-L., Gao, J.-X., Zou, C.-G., Ma, Y.-C. & Zhang, K.-Q. mir-233 modulates the unfolded protein response in C. elegans during Pseudomonas aeruginosa infection. PLoS Pathog.11, e1004606 (2015).
Derrick, T. et al. Conjunctival microrna expression in inflammatory trachomatous scarring. PLoS Negl. Trop. Dis.7, e2117 (2013).
Derrick, T. et al. Inverse relationship between microrna-155 and-184 expression with increasing conjunctival inflammation during ocular Chlamydia trachomatis infection. BMC Infect. Dis.16, 60 (2015).
Liang, G. et al. Altered microrna expression and pre-mrna splicing events reveal new mechanisms associated with early stage Mycobacterium avium subspecies paratuberculosis infection. Sci. Rep.6, 24964 (2016).
Siddle, K. J. et al. Bacterial infection drives the expression dynamics of micrornas and their isomirs. PLoS Genet.11, e1005064 (2015).
Naqvi, A. R., Fordham, J. B. & Nares, S. mir-24, mir-30b, and mir-142-3p regulate phagocytosis in myeloid inflammatory cells. J. Immunol.194, 1916–1927 (2015).
Liu, Z. et al. Up-regulated microrna-146a negatively modulate Helicobacter pylori-induced inflammatory response in human gastric epithelial cells. Microb. Infect.12, 854–863 (2010).
Li, N. et al.H. pylori related proinflammatory cytokines contribute to the induction of mir-146a in human gastric epithelial cells. Mol. Biol. Rep.39, 4655–4661 (2012).
Liu, Z. et al. Microrna-146a negatively regulates ptgs2 expression induced by helicobacter pylori in human gastric epithelial cells. J. Gastroenterol.48, 86–92 (2013).
Cheng, Y. et al. Downregulation of mir-27a* and mir-532-5p and upregulation of mir-146a and mir-155 in lps-induced raw264. 7 macrophage cells. Inflammation35, 1308–1313 (2012).
Derrick, T. et al. mirnas that associate with conjunctival inflammation and ocular Chlamydia trachomatis infection do not predict progressive disease. Pathog. Dis.75, (2017).
Zheng, L. et al. Differential microrna expression in human macrophages with Mycobacterium tuberculosis infection of beijing/w and non-beijing/w strain types. PLoS ONE10, e0126018 (2015).
Matsushima, K. et al. Microrna signatures in Helicobacter pylori-infected gastric mucosa. Int. J. Cancer128, 361–370 (2011).
Xue, X. et al. Downregulation of micro rna-107 in intestinal cd 11c+ myeloid cells in response to microbiota and proinflammatory cytokines increases il-23p19 expression. Eur. J. Immunol.44, 673–682 (2014).
Xue, F. et al. mir-31 regulates interleukin 2 and kinase suppressor of ras 2 during t cell activation. Genes Immun.14, 127 (2013).
Podolska, A. et al. Profiling micrornas in lung tissue from pigs infected with Actinobacillus pleuropneumoniae. BMC Genom.13, 459 (2012).
Noto, J. M. et al. Strain-specific suppression of microrna-320 by carcinogenic Helicobacter pylori promotes expression of the antiapoptotic protein mcl-1. Am. J. Physiol.-Gastrointest. Liver Physiol.305, G786–G796 (2013).
Acknowledgements
We thank Henrique B. Ferreira, Lucas Moitinho-Silva and Jéssica A. Paes for the anti-SPAse I antibody and the Centro de Microscopia e Microanálise (CMM-UFRGS) for technical support with confocal microscopy. This work was performed using the computing facilities of NLHPC (ECM-02) in Chile, the CC/PRABI/LBBE in France and the France Génomique e-infrastructure (ANR-10-INBS-09-08). This work was supported by grants from CAPES-COFECUB 782/13 and the European Research Council under the European Community’s Seventh Framework Programme (FP7 /2007-2013)/ERC grant agreement no. [247073]10. SGM was the recipient of a CAPES and a CAPES-Cofecub doctoral fellowship and a CAPES post-doctoral fellowship (CAPES-Biologia Computacional - Process Number: 611 23038.010043/2013-02). CM was the recipient of a doctoral fellowship CONICYT BECAS CHILE DOCTORADO 2016/FOLIO 72170320.
Author information
Authors and Affiliations
Contributions
S.G.M., M.G.F., A.Z., and M.F.S. conceived the experiments. S.G.M. conducted the experiments with help from M.G.F., S.R. and F.T. ADG designed and performed the processing of total RNA-seq data. C.M. designed and performed the computational prediction of miRNAs and its targets. M.G.F. performed the statistical analysis of mRNA and miRNA sequencing data. S.G.M., M.G.F. and C.M. reconstructed and refined the regulatory networks. S.G.M., M.G.F., C.M., A.D.G., and L.G. analyzed the results. S.G.M., M.G.F., A.Z., and M.F.S. wrote the manuscript with inputs from all other authors. All authors reviewed the manuscript.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Mucha, S.G., Ferrarini, M.G., Moraga, C. et al. Mycoplasma hyopneumoniae J elicits an antioxidant response and decreases the expression of ciliary genes in infected swine epithelial cells. Sci Rep 10, 13707 (2020). https://doi.org/10.1038/s41598-020-70040-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-020-70040-y
This article is cited by
-
Perspectives for improvement of Mycoplasma hyopneumoniae vaccines in pigs
Veterinary Research (2021)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
