- Short report
- Open Access
Hemagglutinin-neuraminidase and fusion proteins of virulent Newcastle disease virus cooperatively disturb fusion–fission homeostasis to enhance mitochondrial function by activating the unfolded protein response of endoplasmic reticulum and mitochondrial stress
Veterinary Research volume 50, Article number: 37 (2019)
The fusogenically activated F and HN proteins of virulent NDV induce complete autophagic flux in DF-1 and A549 cells. However, the effect of both glycoproteins on mitochondria remains elusive. Here, we found that F and HN cooperation increases mitochondrial biogenesis but does not cause the mitochondria damage. We observed that both glycoproteins change the morphological characteristics and spatial distribution of intracellular mitochondria. F and HN cooperate cooperatively to induce ER stress and UPRmt. Our preliminary data suggested that F and HN cooperatively disturb mitochondrial fusion–fission homeostasis to enhance mitochondrial biogenesis, and eventually meet the energy demand of syncytium formation.
Introduction, methods, and results
Autophagy is a conserved catabolic process that may be selective or nonselective. It delivers cytoplasmic contents, including damaged mitochondria and foreign pathogens engulfed into specific double-membrane autophagosome vesicles, to vacuoles or lysosomes for degradation and cycling . The selective elimination of damaged mitochondria is termed as mitophagy, which is a type of macro-autophagy . Mitophagy contributes to the maintenance of a healthy mitochondrial network pool and the prevention of programmed cell death. Some viral proteins, such as the tat of HIV , impair mitochondrial homeostasis to inhibit virus replication. Other viral proteins, such as the HBx of HBV , manipulate mitochondria to benefit viral propagation. However, the roles of the fusion (F) and hemagglutinin-neuraminidase (HN) glycoproteins of Newcastle disease virus (NDV) in mitochondrial have yet to be studied. Our previous study suggested that the F and HN of virulent NDV synergistically induced significant syncytium formation accompanied with complete autophagic flux in DF-1 and A549 cells . As describe in our previous data , A549 cell line can be used as a cell model to study the function of both NDV glycoproteins (Flag-F and HA-HN) of NDV, which are constructed on the basis of the sequence of an F48E9 strain (Genbank Accession Number: MG456905), a standard velogenic strain. Considering the major roles of mitochondria in cellular functions, including suppling energy, regulating calcium levels, and controlling apoptotic cell death etc., we here investigated the relationship between the mitochondria and syncytium formation or membrane fusion induced by F and HN cooperation.
Most mitochondrial proteins are synthesized by a nuclear genome as precursors in the cytoplasm and imported across mitochondrial membranes by the translocase of protein complexes . This translocase facilitates the importation of proteins through the translocation of the outer membrane (TOM) complex. The complementary translocase of the inner membrane (TIM) complex is responsible for protein transport across the inner membrane and into the mitochondrial matrix . The TIM complex mainly includes the membranes of the mitochondrial carrier family of proteins. The TIM23 complex facilitates the translocation of matrix-targeted proteins in the mitochondrial matrix, where they must reside to function . We initially hypothesized that the F and HN co-expression of virulent NDV may induce damage to mitochondria and further degrade them via selective mitophagy. To clarify whether intracellular mitochondria still maintain integrity with syncytium formation and membrane fusion, we first determined the distribution of intracellular TIM 23 response to both glycoproteins co-expression in A549 cells via indirect immunofluorescence assay (IFA). We then examined the colocalization between intracellular TIM 23 and mitochondria marked by two specific targeted technical plasmids, namely, DsRed-Mito (Figure 1A) and EGFP-Mito (Figure 1B). Unexpectedly, we found that their colocalization was accompanied with the extended co-transfection time of F and HN, suggesting that intracellular mitochondria had structural and functional integrity (Figure 1). Statistically, we calculated the Pearson’s correlation coefficient (Pearson’s_Rr) and the Overlap coefficient (Overlap_R) of the zoom out section by using the Image J to assess the correlation of the intensity distribution between channels and the true degree of co-location, respectively. Given the Pearson’s_Rr and Overlap_R (Figure 1), we confirmed that the co-expression of both glycoproteins did not cause intracellular mitochondria damage. To further determine whether selective mitophagy was induced by both glycoproteins, we performed Western blot assay to examine the intracellular protein level of TIM23 and TOM20, which is a TOM protein. In Figure 2A, both protein levels were slightly increased compared with those of mock control, suggesting that intracellular mitochondria were not damaged. Meanwhile, wildtype F and mutant F groups did not also impair the expression of both proteins at 36 h post-transfection (hpt), compared with those of control group. These results suggested that F and HN cooperation neither caused mitochondria damage nor induced mitophagy in A549 cells.
Mitochondria, the main energy hub of the cells, generate a constant flow of electrons through multiple parallel and alternative pathways that eventually end up in the mitochondrial electron transport chain. In mammals, the oxidative phosphorylation machinery occurs in the mitochondrial electron transport chain and is catalyzed by five membrane-bound protein complexes, namely, NADH-ubiquinol oxidoreductase (Complex I), succinate-ubiquinol oxidoreductase (Complex II), ubiquinol–cytochrome c oxidoreductase (Complex III), cytochrome c oxidase (COX) (Complex IV), and ATP synthase (Complex V) [7,8,9]. To further confirm the effect of both NDV glycoproteins on the mitochondria, we determined some subunits of the five protein complexes through Western blotting assay. As expected, all of the detected protein levels did not decrease in response to the F and HN co-expression in A549 cells (Figure 2B) and DF-1 cells (Figure 2C), suggesting that the mitochondria were not damaged, thereby maintaining their structural integrity. NDUSF1, NDUSF3, SDHA, UQCRC2, COX II, COX IV, ATP5A1 and ATP5F1 increased, especially NDUFS1, NDUFS3, COX II, COX IV, and ATP5A1, compared with those in the mock control, indicating that F and HN co-expression might influence the five membrane-bound protein complex functions. TOM20, NDUFS1, NDUFS3, UQCRC1, UQCRC2, COX II, COXIV, and ATP5F1, especially COX II (Figure 2E), in the wildtype F group increased at 36 hpt, compared with those in the control group and mutant group. Those data suggested that activated F1 and HN co-expression did not damage the integrity and function of intracellular mitochondria. Previous studies suggested that the ratio of the mitochondrial DNA (mtDNA) encoded subunit II of COX subunit II (COX II) to nuclear DNA(nDNA)-encoded COX I subunit IV (COX IV) is a surrogate of the translation efficiency of mitochondrial respiratory chain and mitochondrial protein synthesis [10,11,12]. Thus, we normalized and calculated the relative ratio of protein gray density (COX II to COX IV) by using Image J. In Figure 1D, we observed a time-dependent increase in the ratio of COX II to COX IV, suggesting the co-transfection of both glycoproteins enhanced the translation efficiency of mitochondrial respiratory chain and protein synthesis. In addition, the ratio of COX II/COX IV of the wildtype of F group was higher than that of mutant control group and mock group (Figure 2F), suggesting that the cleavage site motif of F influenced the translation efficiency of the mitochondrial electron transport chain. We determined the change in the intercellular adenosine 5’-triphosphate (ATP) level change via luciferin-luciferase bioluminescence method by using an enhanced ATP kit (Beyotime company, S0027). In Figures 2G and 2H, F and HN co-expression increase the intracellular ATP level. In short, F with HN cooperated synergistically to enhance mitochondrial biogenesis, including ATP production.
The morphological characteristics and spatial organization of mitochondria are keys to mitochondrial quality control at protein and sub-organelle levels . As shown in Figure 1A (vertical lines 1 and 5) and Figure 1B (vertical lines 7 and 10), F and HN synergistically altered cellular mitochondrial distribution in a time-dependent manner. To further identify the morphological characteristics and spatial distribution of cellular mitochondria, we conducted a time-course experiment to analyze mitochondrial movement by using DsRed-Mito plasmid (Figure 2A). We examined the mitochondrial network by obtaining original fluorescence microscopic images. The nuclear aggregation of the mitochondria was observed (Figure 3A), and this finding was consistent with the results shown in Figures 1 and 4A, suggesting that the functional interaction of both glycoproteins changed the subcellular movements and spatial distribution of intracellular mitochondria. Normal mitochondria are dynamic organelles, and them form interconnected tubular networks [13, 14]. The alteration of mitochondrial distribution is often accompanied with morphological and functional changes. Considering the aggregated mitochondria around syncytium formation site (Figure 1 vertical lines 1 and 5, Figures 3A, 4A horizontal line 1), we inferred that the spatial distribution of the mitochondria facilitated syncytium formation and membrane fusion induced by both glycoproteins. To obtain more accurate and specific results, we objectively analyzed the mitochondrial network morphology by using Image J as mentioned by Valente et al. . On the basis of the fluorescent microscopy (Figure 3A) and the morphological skeleton analyses of the mitochondria (Figure 3B), we found that F and HN co-expression induced the fragmentation and hyper-fusion of the mitochondria in time-dependent manner. This observation was consistent with fluorescence results of Figure 1 (vertical line 1 and line 7) and Figure 4A (horizontal line 1). In general, mitochondrial elongation process is associated with the dimerization and activation of the ATPase function to produce additional energy [13, 16]. Considering the positive relationship of syncytium formation and the hyper-fusion of mitochondria in a time-dependent manner (Figure 1 vertical lines 1 and 7, Figures 3A, B, 4A horizonal line 1), we speculated that F and HN co-operation induced hyper-fusion of the mitochondria to maintain the energy demand of syncytium formation and membrane fusion. Mitochondrial homeostasis, including its architecture and functions, is maintained by two interlinked but distinct processes, namely, fission and fusion . Mitochondrial fusion involves two sets of key GTPase proteins in mammals: the TOM mitofusin (Mfn) GTPases (Mfn1 and Mfn2)  and the TIM membrane optic atrophy 1(OPA1) , which is a dynamin-related guanosine triphosphatase mutated in dominant optic atrophy. The Mfn GTPases (Mfn1 and Mfn2) mediate TOM fusion, whereas the OPA1 mediates TIM fusion and crista integrity, resulting in the concomitant mixing of the mitochondrial contents and the merging of two individual mitochondria. To further assess whether F and HN change fusion–fission processes, we examined the related fusion protein changes through Western blot. The co-expression of both glycoproteins could increase Mfn1, Mfn2, and OPA1 protein levels (Figures 3C and D), especially Mfn1 and OPA1, compared with those of the control group, suggesting that mitochondrial fusion activity was stimulated. We also observed that the cleavage activity of F influenced the TIM and TOM fusion processes (Figures 3E and F). Mitochondrial fission is a complex process, including two distinct steps, namely, initial mitochondrial membranes constriction and membrane scission. Mitochondria membrane scission is mainly regulated by cytosolic dynamin-related protein 1 (DRP1), a member of the dynamin superfamily of GTPases . The recruitment of DRP1 from the cytoplasm to the mitochondrial surface is a key step in fission regulation in mammalian and yeast cells. Once recruited to the mitochondria, DRP1 oligomerizes as rings on mitochondrial tubules and causes the scission of TOM and TIM. The recruitment of DRP1 in mammalian cells requires several accessory proteins, such as the mitochondrial fission protein 1 (Fis-1) and mitochondrial fission factor (Mff) . Here, we also examined the changes in the levels of related fission proteins by using Western blot technique. F and HN proteins significantly increased Mff, Fis1, and Drap1 levels compared with those of the control group (Figure 3C), suggesting that the co-expression of both glycoproteins could increase the mitochondrial fission activity. Those data indicated that F and HN synergistically disturbed intracellular mitochondrial fusion–fission homeostasis, thereby stabilizing of membrane potential and mixing of matrix proteins, and this observation was consistent with the findings in Figure 1 (vertical line 1 and line 7). The cleavage site of the mutant of F influenced the intracellular fusion–fission protein level (Figure 3E). Therefore, the co-expression of both glycoproteins of NDV synergistically altered the fusion–fission homeostasis of intracellular mitochondria. However, the complicated regulatory mechanism should be further elucidated in future research.
Endoplasmic reticulum (ER) and mitochondria are dynamic, and they undergo continuous structural and spatial reorganization in response to cellular stress . The ER-mitochondria contract serves as a platform for inter-organellar communication, which is necessary to regulate mitochondrial bioenergetic and function . In response to various stresses, cellular ER physiology is mainly regulated by unfolded protein response (UPR), which consists of three integrated signaling branches that activate the downstream of the transmembrane ER stress sensor proteins: inositol requiring enzyme 1 (IRE1), activating transcription factor6 (ATF6), and double-stranded RNA-activated protein kinase(PKR)-like ER kinase (PERK) [23, 24]. Each pathway uses a different mechanism of signal transduction to drive the transcription of UPR target genes, such as ATF6 branch by regulated proteolysis, PERK branch by translational control, and IRE1 branch by nonconventional mRNA splicing . The molecular quality control of mitochondria is regulated through mitochondrial UPR (UPRmt) response to various stresses, an adaptive response to proteotoxic stress, specifically in mitochondria [25, 26]. HSP60 is a mitochondrial chaperone responsible for the transportation and refolding of proteins from the cytoplasm to the mitochondrial matrix . HSP60 is also generally used to evaluate the activation state of the UPRmt . Here, we speculated that activated F1 and HN induced ER stress and UPRmt accompanied with membrane fusion and syncytium formation. To test this hypothesis, we co-transfected EGFP-Mito, DsRed-ER (a specific targeted ER plasmid), and both glycoproteins into A549 cells. In Figure 4A (horizontal line 4), excessive F and HN co-expression increased the interaction and tether between ER and mitochondria (yellow denoted colocalization), indicating the re-organization movement and spatial distribution of both organelles. The values of Overlap_R and plot profile also suggested that F and HN co-expression increased the true degree of co-location between the ER and the mitochondria. Western blot analysis revealed that activated F1 and HN synergistically stimulated the phosphorylation of IRE1 (IRE1 signal pathway transducer), eukaryotic initiation factor 2α (eIF2α, a PERK signal pathway transducer), and Bip (also known as GRP78, a chaperone of the heat shock protein [HSP]70 family) (both belong to ATF6 signal pathway transducers) [24, 26] (Figure 4B). This finding suggested that the overexpression of both glycoproteins could cause ER stress by activating the three branches of the UPR signaling pathway. On the contrary, cleavage site mutation influenced ER stress activation, indicating that the membrane fusion and syncytium formation ability determined the ER stress activation (Figure 4C). To assess UPRmt, we determined the protein change of HSP60 in response to the co-transfection of both glycoproteins through Western blot analysis. The protein expression of HSP60 was significantly increased in a time-dependent manner (Figure 4B), indicating UPRmt activation. Similarly, we found that the cleavage ability of F influenced the UPRmt activation (Figure 4C). On the basis of our results, we proposed a working model of the regulatory mechanism of mitochondrial homeostasis by the HN and F of NDV. Synergistic disruption of fusion–fission homeostasis by HN and F of NDV enhance mitochondiral function by activating of ER UPR and UPRmt (Figure 4D).
Mitochondria are dynamic organelles that maintain the fusion–fission homeostasis in response to intracellular and extracellular stresses . The maintenance of mitochondrial fusion–fission homeostasis is a complex and multistep cellular process, which must be maintained and regulated at protein and sub-organelle levels. This study aimed to investigate the relationship between syncytia generated by of F and HN of virulent NDV and intercellular mitochondria. We found that both glycoproteins did not induce selective autophagy to degrade the damaged mitochondria. Instead, they increased the biogenesis of mitochondria, along with syncytium formation and membrane fusion (Figures 1, 3). On the basis of the positive correlation of microscopic and Western blot analyses (Figures 1, 3, 4A), we speculated that the co-expression of both glycoproteins changed the mitochondrial distribution and morphological characteristics by disturbing mitochondrial fusion–fission homeostasis to meet the interaction energy demand of syncytium formation and membrane fusion. In this process, host cells must precisely regulate the mitochondrial distribution and morphological characteristics to meet energy demand and maintain the interactive activities of F and HN cooperation. However, the elaborate molecular mechanism should be further investigated.
The ER stress can be indirectly transmitted to the mitochondria by altering the transfer of metabolites, such as Ca2+, or directly influencing the mitochondria via the PERK-mediated UPR . In this study, we observed that F and HN co-expression could increase eIF2α levels in A549 (Figure 4B) and DF-1 cells (data not shown) in a time-dependent manner. This finding suggested that PERK-mediated UPR phosphorylated eIF2α subunit, leading to transient translational attenuation. Therefore, we speculated that the PERK-mediated ER stress integrated transcriptional and translational signaling to protect mitochondrial function, and this finding was consistent with the PERK-dependent remodeling of mitochondrial quality control pathways to prevent the potentially pathologic accumulation of misfolded or damaged proteins during ER stress [16, 22]. Recently, Lebeau et al.  reported that stress-induced mitochondrial hyper-fusion is triggered by PERK-dependent translation attenuation via the promotion of the electron transport chain, suggesting that cellular mitochondria could autonomously integrate by remodeling their morphological characteristics and distribution response to ER stress. We inferred that hyper-fusion was induced by co-expression of both glycoproteins via the PERK-mediated stress signaling pathway to maintain energy demand. However, the underlying mechanism should be elucidated in detail.
In summary, we at first reported that both glycoproteins of NDV neither damaged mitochondria nor induced mitophagy. The co-transfection of both viral glycoproteins changed mitochondrial spatial distribution and disturbed the mitochondrial fusion–fission homeostasis, along with syncytium formation and membrane fusion. To meet high-demand energy of syncytia formation, activated virulent F1 and HN synergistically enhanced mitochondrial biological function and eventually supplied energy by activating of UPR of ER and UPRmt. We speculated that the activated UPR of ER and UPRmt might represent a cellular protective mechanism to facilitate the refolding of imported proteins and the proper assembly of unfolded polypeptides in the mitochondria and the ER. Overall, our preliminary data could help to elucidate the molecular mechanism of membrane fusion and syncytium formation induced by F and HN of NDV.
Tanida I (2011) Autophagosome formation and molecular mechanism of autophagy. Antioxid Redox Signal 14:2201–2214
Rozzi SJ, Avdoshina V, Fields JA, Mocchetti I (2018) Human immunodeficiency virus Tat impairs mitochondrial fission in neurons. Cell Death Discov 4:8
Rahmani Z, Huh KW, Lasher R, Siddiqui A (2000) Hepatitis B virus X protein colocalizes to mitochondria with a human voltage-dependent anion channel, HVDAC3, and alters its transmembrane potential. J Virol 74:2840–2846
Ren S, Rehman ZU, Shi M, Yang B, Qu Y, Yang XF, Shao Q, Meng C, Yang Z, Gao X, Sun Y, Ding C (2019) Syncytia generated by hemagglutinin-neuraminidase and fusion proteins of virulent Newcastle disease virus induce complete autophagy by activating AMPK-mTORC1-ULK1 signaling. Vet Microbiol 230:283–290
Chacinska A, Koehler CM, Milenkovic D, Lithgow T, Pfanner N (2009) Importing mitochondrial proteins: machineries and mechanisms. Cell 138:628–644
Sirrenberg C, Bauer MF, Guiard B, Neupert W, Brunner M (1996) Import of carrier proteins into the mitochondrial inner membrane mediated by Tim22. Nature 384:582–585
Enriquez JA (2016) Supramolecular organization of respiratory complexes. Annu Rev Physiol 78:533–561
Chouchani ET, Pell VR, James AM, Work LM, Saeb-Parsy K, Frezza C, Krieg T, Murphy MP (2016) A unifying mechanism for mitochondrial superoxide production during ischemia–reperfusion injury. Cell Metab 23:254–263
Guaras A, Perales-Clemente E, Calvo E, Acin-Perez R, Loureiro-Lopez M, Pujol C, Martinez-Carrascoso I, Nunez E, Garcia-Marques F, Rodriguez-Hernandez MA, Cortes A, Diaz F, Perez-Martos A, Moraes CT, Fernandez-Silva P, Trifunovic A, Navas P, Vazquez J, Enriquez JA (2016) The CoQH2/CoQ ratio serves as a sensor of respiratory chain efficiency. Cell Rep 15:197–209
Moren C, Garrabou G, Noguera-Julian A, Rovira N, Catalan M, Hernandez S, Tobias E, Cardellach F, Fortuny C, Miro O (2013) Study of oxidative, enzymatic mitochondrial respiratory chain function and apoptosis in perinatally HIV-infected pediatric patients. Drug Chem Toxicol 36:496–500
Walker UA, Setzer B, Venhoff N (2002) Increased long-term mitochondrial toxicity in combinations of nucleoside analogue reverse-transcriptase inhibitors. AIDS 16:2165–2173
Chariot P, Bonne G, Authier FJ, Marsac C, Gherardi R (1994) Expression of cytochrome c oxidase subunits encoded by mitochondrial or nuclear DNA in the muscle of patients with zidovudine myopathy. J Neurol Sci 125:190–193
Wai T, Langer T (2016) Mitochondrial dynamics and metabolic regulation. Trends Endocrinol Metab 27:105–117
Twig G, Hyde B, Shirihai OS (2008) Mitochondrial fusion, fission and autophagy as a quality control axis: the bioenergetic view. Biochim Biophys Acta 1777:1092–1097
Valente AJ, Maddalena LA, Robb EL, Moradi F, Stuart JA (2017) A simple ImageJ macro tool for analyzing mitochondrial network morphology in mammalian cell culture. Acta Histochem 119:315–326
Lebeau J, Saunders JM, Moraes VWR, Madhavan A, Madrazo N, Anthony MC, Wiseman RL (2018) The PERK arm of the unfolded protein response regulates mitochondrial morphology during acute endoplasmic reticulum stress. Cell Rep 22:2827–2836
Chan DC (2012) Fusion and fission: interlinked processes critical for mitochondrial health. Annu Rev Genet 46:265–287
Santel A, Fuller MT (2001) Control of mitochondrial morphology by a human mitofusin. J Cell Sci 114:867–874
Olichon A, Baricault L, Gas N, Guillou E, Valette A, Belenguer P, Lenaers G (2003) Loss of OPA1 perturbates the mitochondrial inner membrane structure and integrity, leading to cytochrome c release and apoptosis. J Biol Chem 278:7743–7746
Smirnova E, Griparic L, Shurland DL, van der Bliek AM (2001) Dynamin-related protein Drp1 is required for mitochondrial division in mammalian cells. Mol Biol Cell 12:2245–2256
Otera H, Wang C, Cleland MM, Setoguchi K, Yokota S, Youle RJ, Mihara K (2010) Mff is an essential factor for mitochondrial recruitment of Drp1 during mitochondrial fission in mammalian cells. J Cell Biol 191:1141–1158
Rainbolt TK, Saunders JM, Wiseman RL (2014) Stress-responsive regulation of mitochondria through the ER unfolded protein response. Trends Endocrinol Metab 25:528–537
Ron D, Walter P (2007) Signal integration in the endoplasmic reticulum unfolded protein response. Nat Rev Mol Cell Biol 8:519–829
Walter P, Ron D (2011) The unfolded protein response: from stress pathway to homeostatic regulation. Science 334:1081–1086
Baker MJ, Tatsuta T, Langer T (2011) Quality control of mitochondrial proteostasis. Cold Spring Harb Perspect Bio 3:a007559
Haynes CM, Ron D (2010) The mitochondrial UPR—protecting organelle protein homeostasis. J Cell Sci 123:3849–3855
Jindal S, Dudani AK, Singh B, Harley CB, Gupta RS (1989) Primary structure of a human mitochondrial protein homologous to the bacterial and plant chaperonins and to the 65-kilodalton mycobacterial antigen. Mol Cell Biol 9:2279–2283
Jovaisaite V, Mouchiroud L, Auwerx J (2014) The mitochondrial unfolded protein response, a conserved stress response pathway with implications in health and disease. J Exp Biol 217:137–143
Smith EC, Popa A, Chang A, Masante C, Dutch RE (2009) Viral entry mechanisms: the increasing diversity of paramyxovirus entry. FEBS J 276:7217–7227
Chang A, Dutch RE (2012) Paramyxovirus fusion and entry: multiple paths to a common end. Viruses 4:613–636
This work was funded by the National Key Research and Development Program of China (No. 2018YFD0500100) and National Natural Science Foundation of China (Nos. 31530074, 31872453). We thank Prof. Sa Xiao (Northwest A& F University) for giving viral plasmids (Flag-F and HA-HN) and Dr. Chongyang Wang and Na Huo (Northwest A& F University) for technical support. We thank Dr. Yanqing Jia (Northwest A& F University) for moral encouragement.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Ren, S., Rehman, Z.U., Shi, M. et al. Hemagglutinin-neuraminidase and fusion proteins of virulent Newcastle disease virus cooperatively disturb fusion–fission homeostasis to enhance mitochondrial function by activating the unfolded protein response of endoplasmic reticulum and mitochondrial stress. Vet Res 50, 37 (2019). https://doi.org/10.1186/s13567-019-0654-y