Large yellow croaker (Larimichthys crocea) mitofusin 2 inhibits type I IFN responses by degrading MAVS via enhanced K48-linked ubiquitination

In mammals, mitofusin 2 (MFN2) is involved in mitochondrial fusion, and suppresses the virus-induced RIG-I-like receptor (RLR) signaling pathway. However, little is known about the function of MFN2 in non-mammalian species. In the present study, we cloned an MFN2 ortholog (LcMFN2) in large yellow croaker (Larimichthys crocea). Phylogenetic analysis showed that MFN2 emerged after the divergence of amphioxus and vertebrates. The protein sequences of MFN2 were well conserved from fish to mammals. LcMFN2 was expressed in all the tissues/organs examined at different levels, and its expression was upregulated in response to poly(I:C) stimulation. Overexpression of LcMFN2 inhibited MAVS-induced type I interferon (IFN) promoter activation and antiviral gene expression. In contrast, knockdown of endogenous LcMFN2 enhanced poly(I:C) induced production of type I IFNs. Additionally, LcMFN2 enhanced K48-linked polyubiquitination of MAVS, promoting its degradation. Also, overexpression of LcMFN2 impaired the cellular antiviral response, as evidenced by the increased expression of viral genes and more severe cytopathic effects (CPE) in cells infected with spring viremia of carp virus (SVCV). These results indicated that LcMFN2 inhibited type I IFN response by degrading MAVS, suggesting its negative regulatory role in cellular antiviral response. Therefore, our study sheds a new light on the regulatory mechanisms of the cellular antiviral response in teleosts. Supplementary Information The online version contains supplementary material available at 10.1007/s42995-023-00189-8.


Introduction
The constant threat of invading pathogens into host cells has driven the evolution of the innate immune system to quickly detect and respond to conserved pathogen-associated molecular patterns (PAMPs) (Buchmann 2014;Kumar et al. 2011;Medzhitov and Janeway 2002;Netea et al. 2004).Cytoplasmic retinoic acid-inducible gene I-like receptors (RLRs) are a type of pattern recognition receptors (PRRs) that sense virus-derived dsRNA and trigger antiviral responses in immune and other cells, such as epithelial and gonad cells (Chen and Hur 2022;Kawai and Akira 2008;Takeuchi and Akira 2009).The RLR family, which belongs to the DExD/H box RNA helicases superfamily, includes three members: RIG-I, melanoma differentiation-associated gene 5 (MDA5), and laboratory of genetics and physiology 2 (LGP2).Members of this family sense viral nucleic acids using a C-terminal helicase domain (Bruns and Horvath 2012;Takahasi et al. 2009).Upon viral recognition, RIG-I and MDA5 activate the downstream mitochondrial antiviralsignaling protein (MAVS), a mitochondrial outer membrane protein (Baril et al. 2009;Seth et al. 2005), and promote its oligomerization (Buskiewicz et al. 2016;Onomoto et al. 2021;Wu and Hur 2015).Conversely, LGP2 interacts with RIG-I or MDA5 to control signal transduction (Gong et al.Edited by Xin Yu. 2022; Satoh et al. 2010;Zhang et al. 2018).Then, MAVS recruits TNF receptor-associated factor 3 (TRAF3) and TANK-binding kinase 1 (TBK1), which phosphorylate and activate interferon regulatory factor 3 (IRF3) and IRF7 (Liu et al. 2015).After being phosphorylated, IRF3 and IRF7 homodimerize or heterodimerize and translocate into the nucleus, where they bind to the interferon-stimulated response element (ISRE) within the IFN promoter, leading to the early production of type I IFN (Dalskov et al. 2020;Ding et al. 2016;Feng et al. 2016;Honda and Taniguchi 2006).The secretion of type I IFNs then activates the Janus kinase/signal transducer and activator of transcription (JAK/ STAT) pathway, resulting in the transcription of hundreds of IFN-stimulated genes (ISGs), including many antiviral effectors, such as myxovirus resistance protein (Mx), virus inhibitor protein, endoplasmic reticulum-related factor, interferon-inducible factor (Viperin), protein kinase R (PKR), and ISG15 (Au-Yeung et al. 2013;Raftery and Stevenson 2017).
Mitochondria, which are the double-membrane organelles commonly referred to as "powerhouses" of eukaryotic cells, maintain a constantly changing morphology through fission and fusion events that require the coordination of a complex group of proteins (Berman et al. 2008;Davey and Clark 1996;Mills et al. 2017;Sabouny and Shutt 2020).Mitochondrial fusion is mediated by outer and inner mitochondrial membrane proteins, including the GTPases MFN1 and MFN2, which are localized on the outer mitochondrial membrane (OMM), and the Opa1, in the inner mitochondrial membrane (IMM) (Chen et al. 2003a;Lee et al. 2004).Conversely, mitochondrial elongation and fission are regulated by the interaction of dynamin-related protein 1 (Drp1) and other OMM proteins, such as MiD51 and Fis1 (Lee et al. 2004;Palmer et al. 2011;Yu et al. 2019).
Growing evidence suggests that mitochondrial fusion is linked to the antiviral response and modulation of RLR signaling (Castanier et al. 2010;Kim et al. 2018).For example, mitochondrial fragmentation reduces virusinduced IRF3 and NF-κB-dependent signal transmissions and impairs host antiviral activity against RNA viruses (Onoguchi et al. 2010).Furthermore, in fusion-deficient mouse embryo fibroblasts lacking both MFN1 and MFN2, virus-induced RLR-dependent production of IFN-β and cytokine interleukin-6 (IL-6) was drastically reduced, resulting in the increased viral replication (Chen et al. 2003a).There are some discrepancies between MFN1 and MFN2.Although deletion of MFN2 in mouse embryonic fibroblast (MEF) cell lines reduces vesicular stomatitis virus (VSV)-induced IFN-β production, manipulating MFN1 expression levels did not yield similar phenotypes (Yasukawa et al. 2009).This suggests that MFN2 functions as a negative regulator of RLR signaling.Also, recent study has revealed that inhibiting MFN2 expression with epigoitrin leads to upregulation of influenza virus-induced MAVS and IFN-β activation (Luo et al. 2019).Additionally, MFN2 has been shown to inhibit the RIG-I/IRF7 signaling pathway, resulting in decreased Seneca Valley virus (SVV)-induced IFN-λ3 production (Deng et al. 2022).
While the importance of MFN proteins in modulating RLR signaling pathways has been extensively studied in mammals, their function in non-mammalian species remains largely unknown.In the current study, we cloned an MFN2 ortholog (LcMFN2) from a large yellow croaker (Larimichthys crocea), which is a commercially important mariculture species in East Asia.The fish suffer from infectious diseases caused by Pseudomonas plecoglossicida, Cryptocaryon irritans, and large yellow croaker iridovirus (Chen et al. 2003b;Li et al. 2020a, b;Nelapati et al. 2012;Yin et al. 2019).The expression patterns of LcMFN2 under normal and poly(I:C)-stimulated-conditions were examined.Moreover, we demonstrated that LcMFN2 plays a negative regulatory role in MAVS-mediated type I IFN signaling.Mechanically, LcMFN2 interacts with LcMAVS and promotes its degradation through the ubiquitin-proteasome pathway.These findings contribute significantly to our understanding of the regulatory mechanisms of the cellular antiviral response in teleosts.

Sequence characterization and phylogenetic analysis
A search was performed in the genomic database (GCF_000972845.2) (Ao et al. 2015) for MFN2, named LcMFN2 in large yellow croaker, using TBLASTN, and its sequence was confirmed by Sanger sequencing.The openreading frame (ORF) of LcMFN2 was 2277 base pairs (bp) in length, encoding 759 amino acids (aa) with a calculated molecular mass of 86 kDa.Analysis of the genome data indicates that LcMFN2 is located on primary assembly XV, and spans an approximately 11.81 kb chromosome contig.In addition, the coding sequence of LcMFN2 is arranged in 17 exons and 16 introns, with all intron-exon boundaries following the typical GT/AG pattern, similar to its mammalian counterparts (Supplementary Table S3).Domain analysis revealed that the LcMFN2 protein contains highly conserved functional domains, including an N-terminal GTPase domain, a heptad repeat 1 (HR1) domain, a proline-rich (PR) domain, a transmembrane (TM) domain, and a C-terminal HR2 domain (Supplementary Fig. S1A and B).The protein structure of LcMFN2 is comparable to that of human MFN2, with average root-mean-square deviations (RMSD) of 4.701 (rigid alignment) and 2.41 (flexible alignment), respectively (Supplementary Fig. S1C).However, the amino acid sequence of LcMFN2 shows marked differences compared to that of vertebrate MFN1 and non-chordate metazoa MFN2/MFN2-like, particularly at the C-terminus, with a lowest similarity of < 60% (Fig. 1A).
Phylogenetic analysis indicates that MFN genes from non-chordate metazoans formed the basal branch, suggesting that these genes encode ancestral MFN proteins (Fig. 1B).Nevertheless, these genes were incorrectly annotated as MFN2/MFN2-like in the NCBI (National Center for Biotechnology Information, https:// www.ncbi.nlm.nih.gov/) database.MFN1 and MFN2 from urochordates (Styela clava), cephalochordates (amphioxus), and vertebrates (cartilaginous fish, bony fish, amphibians, reptiles, birds, and mammals) formed  S2.A Sequence comparison of MFN coding regions among various species, encompassing nonchordate metazoa and chordates.The color green denotes 100% similarity, yellow-green denotes 80-100% similarity, yellow denotes 60-80% similarity, and gray denotes less than 60% similarity.B Phylogenetic relationships among MFN sequences.The phylogenetic tree was inferred using the neighbor-joining method in MEGA11 with the Jones-Taylor-Thornton (JTT) model, and 10,000 bootstrap replicates.The branches were labeled with the bootstrap values the sister group of the MFN genes from non-chordate metazoans.Two sequences from B. floridae and S. clave, annotated as MFN2-like, were found cluster outside the vertebrate MFN genes.
The time-calibrated phylogenetic tree shows that the chordate MFN genes diverged from the ancestral MFN genes ~ 522 million years ago (MYA), whereas vertebrate MFN genes emerged ~ 394 MYA (Fig. 2).Vertebrate MFN1 and MFN2 diverged at ~ 356 MYA and 310 MYA, respectively (Fig. 2).Taking together, our analysis suggested that MFN2 emerged after the divergence of amphioxus and vertebrates, and its amino acid sequence has remained highly conserved from fish to mammals, indicating that its functional role has been maintained across a wide range of species over a long period of evolutionary time.

Expression analysis of the LcMFN2 gene
Gene expression analysis showed that LcMFN2 was expressed in all examined tissues and organs with different expression levels.LcMFN2 expression level was highest in the liver, followed by the gill, skin, heart, muscle, brain, stomach, blood, and spleen, with the lowest level detected in the intestine and head kidney (Fig. 3A).
Challenge experiments were carried out to assess whether LcMFN2 could be induced by poly(I:C) stimulation.As shown in Fig. 3B-D, poly(I:C) stimulation led to significantly increased levels of LcMFN2 expression in the head kidney, spleen, and intestine after 6 h postinfection (hpi), with fold increases of 1.63, 1.99, and 1.70, respectively, compared to the control group.By 12 hpi, LcMFN2 expression was consistently higher in all of the Fig. 2 Ancestral states of MFN genes were estimated using the RelTime method in MEGA11.The numbers and horizontal bars at each node indicate the presumed divergence time in MYA and 95% confidence interval (CI), respectively tested tissues/organs, with fold increases of 5.66 (head kidney), 7.69 (spleen), and 2.49 (intestine) compared to the control group.At 24 hpi, LcMFN2 expression was persistently elevated in all the examined tissues/organs, being 15.02 (head kidney), 48.09 (spleen), and 9.47 (intestine)-fold higher than in the control group, respectively.At 48 hpi, significant increases were also observed in all the tested tissues/organs, with the head kidney being 20.73-fold, the spleen being 26.34-fold, and the intestine being 6.67-fold higher than the control, respectively.

LcMFN2 inhibits LcMAVS-induced type I IFN signaling
To investigate whether MFN2 plays a role in poly I:Cinduced IFN production, we employed siRNA to silence LcMFN2.Our results showed that MFN2 expression was reduced by ~ 63% (siMFN2-2) in large yellow croaker head kidney (LYCK) cells (Wang et al. 2014) transfected with siRNA oligos compared to control cells (Fig. 4A).Subsequently, qPCR analysis demonstrated that MFN2 knockdown enhanced poly(I:C)-induced expression of IFNi (Fig. 4B) and IFNd (Fig. 4C).Previous studies have shown the regulatory role of mammalian MFN2 in the antiviral responses elicited by MAVS (Yasukawa et al. 2009).Therefore, we investigated how MFN2 regulates the activation of MAVS-induced IFN promoters in fish.Luciferase activity assay was performed, employing LcIFNi and LcIFNd as reporter genes, with LcIFNi being classified as a group II fish IFN and LcIFNd being classified as a group I fish IFN (Chen et al. 2021(Chen et al. , 2022)).As shown in Fig. 4D, E, both LcMAVS and LcTBK1 significantly activated LcIFN promoters.However, co-transfection with LcMFN2 decreased LcMAVS-induced activation of LcIFN promoters, whereas LcTBK1 had no effect.Previous studies have shown that overexpressing fish MAVS may boost the antiviral activity of epithelioma papulosum cyprini (EPC) cells (Chen et al. 2017; LcMFN2 was measured at 3, 6, 12, 24, and 48 h post-stimulation/ infection and compared to a control group treated with PBS.The expression of LcMFN2 in each tissue was normalized to that of β-actin and presented as fold changes relative to control fish injected with the same volume of PBS.The data are presented as mean ± SD (N = 3), with asterisks (*) denoting a significant difference (*P < 0.05, **P < 0.01, ***P < 0.001) above each bar Zhou et al. 2015).Thus, we investigated the impact of LcMFN2 on MAVS-induced antiviral gene expression.Our results showed that overexpression of LcMFN2 significantly inhibited LcMAVS-induced EPC antiviral gene expression (Fig. 4F).Then, we tested the interaction between LcMFN2 and LcMAVS.Human embryonic kidney 293T (HEK293T) cells were co-transfected with FLAG-tagged LcMFN2 and His-tagged LcMAVS.The Co-IP experiment revealed that LcMFN2 could interact with LcMAVS, suggesting that this inhibitory impact may result from protein-protein interaction (Fig. 4G).Taken together, these findings suggest strongly that MFN2 serves as a negative regulator of RLRs signaling molecules, particularly MAVS-induced IFN production.

LcMFN2 promotes the degradation of LcMAVS by enhancing K48-linked ubiquitination
Based on our previous finding that LcMFN2 interacts with LcMAVS and inhibits type I IFN production, we speculated that LcMFN2 could degrade LcMAVS.To investigate the role of LcMFN2 in regulating LcMAVS protein expression, we co-transfected EPC cells with LcMFN2 and LcMAVS expression plasmids and examined LcMAVS protein expression using immunoblot assays with anti-His-tag antibody.Our results showed that co-transfection with LcMFN2 promoted the degradation of LcMAVS in a dosedependent manner (Fig. 5A).To further confirm the effect of LcMFN2 on the protein stability of LcMAVS, we conducted a cycloheximide (CHX) assay to measure the degradation rate of LcMAVS.As shown in Fig. 5B, the overexpression of LcMFN2 expedited the degradation of LcMAVS in the presence of CHX, resulting in a reduction of LcMAVS's halflife.Eukaryotic cells mainly undergo protein degradation through two pathways, the autophagolysosome pathway and the ubiquitin-proteasome pathway.Therefore, we determined the degradation pathway of LcMAVS by using specific inhibitors, MG132 for the ubiquitin-proteasome system and NH 4 CI for the autophagosome-lysosome system (Guo et al. 2016;D'Eletto et al. 2009).As illustrated in Fig. 5C, LcMFN2-mediated degradation of LcMAVS could be suppressed by MG132 but not by NH 4 CI, demonstrating that overexpression of LcMFN2 results in the degradation of LcMAVS via the ubiquitin-proteasome pathway.An additional discovery supporting this conclusion is that overexpression of LcMFN2 enhanced the polyubiquitination level of LcMAVS (Fig. 5D).To investigate which type of Ub chain structure was affecting LcMAVS, we utilized two mutants: Ub-K48 and Ub-K63, in which six of the seven lysine residues in the Ub molecule were replaced with arginine, leaving only one lysine residue intact.Our results showed a significant increase in K48-linked ubiquitination of LcMAVS under the regulation of LcMFN2, with K63linked Ub chains remaining unchanged (Fig. 5E, F).These data suggest strongly that LcMFN2 specifically induces K48-ubiquitination of LcMAVS, ultimately leading to its proteasomal degradation, which is consistent with previous research that has underscored the importance of K48ubiquitin chains in directing protein substrates toward proteasomal degradation.

Overexpression of LcMFN2 impairs the cellular antiviral response
To investigate the impact of LcMFN2 on the cellular antiviral response, EPC cells were transfected with an FLAG-tagged LcMFN2 expression plasmid or a CMV-FLAG empty vector as a control and subsequently infected with SVCV.At 72 hpi, control cells showed a moderate CPE caused by SVCV, whereas cells transfected with LcMFN2 expression plasmid exhibited a severe CPE (Fig. 6A, B).The effect of LcMFN2 overexpression on the expression of viral genes was also examined (Fig. 6C), and the results showed that viral glycoprotein (G) was increased by 4.38-fold; RNAdependent RNA polymerase (L) was increased by 7.01-fold; matrix protein (M) was increased by 5.53-fold; nucleoprotein (N) was increased by 4.55-fold; and phosphoprotein (P) was increased by 4.04-fold in cells expressing LcMFN2.

Discussion
Recent evidence highlights the mitochondrial fusion protein MFN2 in regulating host immune responses, especially in the virus-induced RLR signaling pathway (Eisenächer and Krug 2012;Tan et al. 2022;Xu et al. 2020).In this study, we observed an increase in the expression level of LcMFN2 following poly(I:C) stimulation, suggesting its potential involvement in the regulation of antiviral immune responses.
Since MFN2 gene identification is scarce in nonmammalian species, we analyzed the sequence of MFN2 in the large yellow croaker.Our analysis showed that the protein structure of the large yellow croaker MFN2 is similar to that of human MFN2, suggesting that fish MFN2 may have similar physiological functions to mammalian MFN2.It may be inferred that the well-conserved protein sequences of MFN2 from fish to mammals are subject to purifying selection, which eliminates deleterious mutations over time and reduces the genetic diversity (Cvijović et al. 2018).However, the C-terminal region of MFN2 in vertebrates is highly divergent from MFN1 in vertebrates and MFN proteins in early metazoans.The accumulation of non-synonymous mutations in the C-terminus may have contributed to the functional diversification of MFNs (Hughes 1994).Interestingly, our investigation has revealed that the MFN genes may have undergone two rounds of duplication.In early chordates, such as tunicates (Styela clava) and cyclostomes (lampreys), only one copy of the MFN gene was found, whereas most basic metazoans had two copies, indicating that the first MFN gene duplication occurred before the emergence of chordates.Furthermore, according to the results of phylogenetic analysis, the ancestral MFN genes constituted the basal branch, with MFN1 and MFN2 clustering together in vertebrates.Thus, it is likely that the subsequent duplication event occurred after the amphioxus-vertebrate split, leading to the divergence of MFN1 and MFN2, which underwent subfunctionalization or neofunctionalization (Chen et al. 2020;Naon and Scorrano 2014).
The current study demonstrated that LcMFN2 participates in regulating the RLR-induced type I IFN signaling by negatively regulating the MAVS-mediated signaling transduction, suggesting that the function of LcMFN2 is similar to its mammalian counterparts (Deng et al. 2022;Yasukawa et al. 2009).Interestingly, LcMFN2 selectively inhibits the promoter activity induced by LcMAVS but has minimal effect on the promoter activity induced by LcTBK1, which may be attributed to the upstream molecular position of LcMFN2 concerning LcTBK1.This observation is consistent with a previous report demonstrating that the knockdown of MFN2 did not affect the expression of NF-κB induced by TRAF6 or the expression of IFN-β induced by TBK1 (Yasukawa et al. 2009).The current After 24 h of transfection, the cells were immunoprecipitated with an anti-FLAG antibody, followed by western blotting using the appropriate antibodies study also provides evidence that LcMFN2 could promote the degradation of LcMAVS via the ubiquitin-proteasome pathway in a dose-dependent manner, suggesting that LcMFN2 may function as a fine-tuner of MAVS activity, providing a mechanism for regulating the intensity and duration of antiviral responses.Generally, K63-linked ubiquitination enhances substrate protein stability, whereas K48-linked ubiquitination promotes substrate protein degradation (Mallette and Richard 2012;Newton et al. 2008;Tracz and Bialek 2021).The ubiquitination of MAVS is controlled by a variety of E3 ubiquitin ligases and deubiquitinases.For example, it has been demonstrated that ring finger protein 125 (RNF125) and tripartite motif-containing protein 25 (TRIM25), two E3 ubiquitin ligases, may target MAVS for K48-linked ubiquitination and degradation (Castanier et al. 2012;You et al. 2009).Recently, Hou et al. (2021) reported that deubiquitinase ubiquitin-specific peptidase 18 (USP18) could stabilize MAVS by promoting K63-linked polyubiquitin chains from MAVS, thereby preventing its proteasomal degradation.The limitation of our current study is that we did not identify the individual E3 ubiquitinates responsible for MAVS degradation.Future research should aim to determine the specific ubiquitinates involved in this process.In conclusion, the MFN2 gene was cloned and identified in L. crocea, and its expression was upregulated by poly(I:C) treatment.LcMFN2 negatively regulated the MAVS-mediated signaling pathway by degrading LcMAVS through the ubiquitin-proteasome pathway, inhibiting the expression of type I IFNs and ISGs (Fig. 7).
These findings provide new insights into the mechanisms governing the cellular antiviral response in teleost.

Bioinformatics analysis
The conserved domain structures were annotated by querying against the Conserved Domain Database (CDD) of the NCBI, and manually checked by comparing with the protein domain of human MFN2.Protein structure prediction was performed using Alphafold2, and the model with the highest confidence (pLDDT) was selected for further analysis (Mirdita et al. 2022).Rigid and flexible comparisons of protein structure were carried out with PyMOL software (version 2.5) and FATCAT 2.0 web service (https:// fatcat.godzi klab.org/), respectively (Li et al. 2020a, b).Multiple sequence alignment was performed using the Geneious alignment plugin in Geneious Prime (version 2022.1.1)(Kearse et al. 2012).The best-fit substitution model for phylogenetic analysis was determined by MEGA11 software based on the aligned sequences.Neighbor-joining trees were constructed using MEGA11 with the Jones-Taylor-Thornton (JTT) model and 10,000 bootstrap replicates (Tamura et al. 2021).The time-based phylogenetic tree was generated using the RelTime-ML method in MEGA11, which is effective in calculating divergence times for duplicated genes.Three calibration constraints were applied, including the estimated divergence time of: (1) Danio rerio and Lepisosteus oculatus between 298.8 and 342.5 MYA, retrieved from the TimeTree5 database (http:// timet ree.org); (2) Hydra vulgaris and Acropora digitifera between 543.0 and 610.0 MYA; and (3) Podarcis muralis and Gallus gallus between 276.0 and 286.8 MYA.Two MFN genes from Amphimedon queenslandica were used as outgroups.The trees were edited for visualization in Adobe Illustrator CC2019.

Animal, cell lines, and virus
Large yellow croakers weighing 100 g were obtained from a mariculture farm in Ningde, Fujian, China, and acclimatized for one week in aquaria with water at 2% salinity and 22 °C as approved by the College of Life Science at Fujian Agriculture and Forestry University.For challenge studies, 45 fish were divided into three groups of 15 and injected intraperitoneally with either poly(I:C) (Sigma-Aldrich, 0.2 mL/100 g) in phosphate-buffered saline (PBS, pH 7.4) or PBS alone.Tricaine methanesulfonate (MS-222, 100 mg/L, Sigma-Aldrich) was used to anesthetize the fish during sample collection.Tissues/organs, including the head kidney, spleen, and intestine, were collected at 3, 6, 12, 24, and 48 h post-injection and promptly frozen in liquid nitrogen until use.

RNA extraction, cDNA synthesis, and qRT-PCR analysis
Total RNA was extracted using Trizol regent (Invitrogen), the manufacturer's protocol.The concentration of total RNA was determined by calculating the absorbance ratio at 260/280 using a Nanodrop Spectrophotometry (Thermo Scientific, USA), and the quality of the RNA was determined by electrophoresis on a 1% agarose gel.To generate the first-strand cDNA for gene amplification and expression analysis, 2 μg of total RNA was reverse-transcribed using the HiScript III first-strand cDNA Synthesis Kit (+ gDNA wiper) (Vazyme Biotech).
Real-time PCR was performed using QuantStudio™ 5 Real-time PCR system (Applied Biosystems).Amplification reactions were carried out in triplicate, and each reaction containing 10 μL of 2 × SYBR qPCR master mix (Vazyme Biotech), 1 μL of the diluted cDNA, 0.4 μL of each primer, and 8.2 μL ddH 2 O.The C p value of each sample was computed automatically by software, and relative fold changes were calculated using the 2 −ΔΔCt method.The PCR amplification efficiency was calculated using the standard curve method, and optimized to over 95%.Primer specificity was evaluated by examining the melt curves.

Gene clone and plasmid construction
Supplementary Table S1 for the full-length ORF of LcMFN2 were designed based on the predicted 5′ upstream and 3′ downstream sequences of LcMFN2 (GenBank accession No. XM_019259949.2).The ORF of LcMFN2 was amplified from the spleen cDNA template according to the instructions of Ex-Taq DNA polymerase (Takara).After gel purification, the PCR products were cloned into the pMD19-T vector (Takara) and subsequently transformed into Escherichia coli DH5α.Putative positive clones were sequenced using the vector-specific primers M13F/M13R.
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 licence, and indicate if changes were made.The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material.If material is not included in the article's Creative Commons licence 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 licence, visit http:// creat iveco mmons.org/ licen ses/ by/4.0/.

Fig. 1
Fig. 1 Sequence comparison and phylogenetic relationships of MFNs from Supplementary TableS2.A Sequence comparison of MFN coding regions among various species, encompassing nonchordate metazoa and chordates.The color green denotes 100% similarity, yellow-green denotes 80-100% similarity, yellow denotes

Fig. 3
Fig. 3 The expression pattern of the LcMFN2 gene.A The expression levels of LcMFN2 in tissues/organs were measured and normalized to the expression of the housekeeping gene β-actin, and the results were reported in arbitrary units.The lowest level of LcMFN2 expression was observed in the head kidney, which was defined as 1.The abbreviations are L, Liver; G, gill; SK, skin; HE, heart; MU, muscle; BR, brain; ST; stomach; BL, blood; SP, spleen; IN, intestine; HK, head kidney.B-D The relative mRNA level of

Fig. 4
Fig. 4 LcMFN2 interacts with LcMAVS and inhibits LcMAVSmediated signaling.A Effects of RNA interference on the expression of endogenous MFN2.LYCK cells were transfected with three different sequences of LcMFN2 siRNA or control siRNA (siNC) for 24 h, and then collected for qPCR.B, C Effects of siRNA-treated LYCK cells in the poly(I:C)-induced IFNs transcription.Cells were transfected with 100 nM si-MFN 2-2 for 24 h, and then transfected with poly(I:C) (1 μg/mL) for 24 h.The transcription of IFNi (B) and IFNd (C) was detected by qPCR.D, E Overexpression of LcMFN2 selectively inhibits LcMAVS, but not LcTBK1, mediated IFN promoter activation.EPC cells were co-transfected with 400 ng LcIFN promoter reporters, LcIFNipro-luc (D) and LcIFNdproluc (E), and plasmids expressing His-tagged LcMFN2 (250 ng) or LcTBK1 (250 ng), in combination with Flag-tagged LcMAVS or empty vector pcDNA™3.1/myc-Hisat 250 ng each.Luciferase

Fig. 5
Fig. 5 LcMFN2 degrades LcMAVS through the ubiquitinproteasome pathway.A LcMFN2 mediated the degradation of LcMAVS.EPC cells were transfected with 2 μg LcMAVS-Flag together with 2 μg empty vector or LcMAVS-His for 24 h.The protein level of LcMAVS was detected by western blotting.B CHX chase analysis of exogenous LcMAVS in LcMFN2 transfected cells or control cells.EPC cells were treated with 100 μg/mL CHX for the specified time period, and the expression levels of LcMAVS were quantified.C The protein degradation of LcMAVS was prevented by treatment with MG132.EPC cells were transfected with 2 μg of Flag-LcMAVS and 2 μg of an empty vector or His-LcMAVS for 24 h, and

Fig. 6
Fig. 6 Overexpression of LcMFN2 impairs the cellular antiviral response.A EPC cells were transfected with LcMFN2 expression plasmid and infected with SVCV.After 72 h post-infection, the CPE was observed microscopically.B The CPE was detected by staining the cells with 1% crystal violet, and the cell viability was assessed by measuring the optical density at a wavelength of 562 nm.C The

Fig. 7
Fig. 7 MFN2 regulates the MAVS-mediated antiviral response in large yellow croaker.Upon recognizing viral PAMPs such as dsRNA, RLRs recruit MAVS to activate TBK1, which phosphorylates and activates IRF3.Then, homodimerized IRF3 translocates to the nucleus, binding to ISRE in the IFN promoter, leading to early