Molecular Characterization of Interferon Regulatory Factor 10 in Red Sea Bream (Pagrus major) and Expression Analysis Following RSIV Infection
Abstract
Interferon regulatory factors (IRFs) exhibit antiviral and antibacterial capabilities. In this research, we isolated the IRF10 genes from red sea bream (Pagrus major) and examined their expression following red sea bream iridovirus (RSIV) infection. Analysis of multiple sequences identified preserved components, such as the IRF-association domain, nuclear localization signal, and DNA-binding domain. Phylogenetic examination showed that PmIRF10 grouped with IRF10 equivalents in fish. We utilized quantitative real-time PCR to explore IRF10 expression in normal tissues, revealing increased expression levels in the gill and spleen. Our results indicate that PmIRF10 plays a role in the immune response, as evidenced by the controlled expression levels of IRFs after RSIV infection in red sea bream. This finding aligns with prior research on different fish species. Overall, this investigation characterizes IRF genes in red sea bream and provides a basis for comparative analysis of IRF genes.
Keywords:
Pagrus major, Red sea bream iridovirus, Interferon regulatory factors, IRF10Ⅰ. Introduction
The development of innate and adaptive immune systems in vertebrates has enabled them to protect themselves against the infiltration of foreign pathogens. These systems encompass endogenous binding sites and the activation of various expression factors, such as interferon regulatory factors (IRFs). Initially identified in humans, IRFs constitute a protein family responsible for regulating the transcriptional and inducible gene expression of interferons (IFNs) (Antonczyk et al., 2019). Since the first discovery of IRF1 in 1988, the scope of IRFs has expanded to comprise IRF1-11. Studies have demonstrated that IRF1-9 can be found in mammals, IRF10 in birds and fish, and IRF11 exclusively in fish (Suzuki et al., 2011; Shu et al., 2015; Guan et al., 2020).
The IRF family is characterized by an N-terminal DNA-binding domain (DBD), a peptide linker (LK), and an IRF-association domain (IAD) that includes five tryptophan repeats essential for identifying similar DNA motifs within the highly conserved DBD (Antonczyk et al., 2019). Phylogenetic analysis categorizes the IRF family into four subfamilies: IRF1-G (IRF1, IRF2); IRF3-G (IRF3, IRF7); IRF4-G (IRF4, IRF8, IRF9, IRF10); and IRF5-G (IRF5, IRF6) (Nehyba et al., 2002). The primary function of the IRF family is to regulate the antiviral response in vertebrates' innate immunity by participating in the transcription of IFN and interferon-stimulated genes (ISG) (Paun & Pitha, 2007). In vertebrates, IRFs identify viral infections by recognizing pathogen-associated molecular patterns (PAMPs) through pattern recognition receptors (PRRs), such as Toll-like receptors (TLRs) (Medzhitov & Janeway Jr, 2000). Furthermore, IRFs differentiate and regulate dendritic cells (Gabriele & Ozato, 2007), and modulate crucial cytokines and transcription factors involved in T helper cell differentiation (Zhang et al., 2012).
Phylogenetic analysis reveals that IRF10 belongs to the IRF-4G subfamily, exhibiting a strong binding affinity to the interferon-stimulated response element (ISRE) domain, particularly when heterodimerized with IRF1, IRF2, or IRF4 (Alter-Koltunoff et al., 2008). Since its initial identification in birds, IRF10 has been detected in various vertebrates, with the exception of humans and mice (Huang et al., 2010). Structurally, IRF10 shares a compact configuration with IRF4 and is reported to play a role in the late stage of antiviral response in birds (Nehyba et al., 2002). IRF10 has been identified in species such as olive flounder, zebrafish, orange-spotted grouper (Epinephelus coioides), Atlantic cod (Gadus morhua), rainbow trout, common carp (Cyprinus carpio), and swamp eel (Monopterus albus) (Suzuki et al., 2011; Li et al., 2014; Qiao-Qing et al., 2014; Huang et al., 2015; Inkpen et al., 2015; Xu et al., 2016; Zhu et al., 2020).
Initially discovered in Japan in 1990 in red sea bream, red sea bream iridovirus (RSIV) is a linear double-stranded DNA virus that belongs to the Iridoviridae family and the Megalocytivirus genus (Inouye et al., 1992). RSIV inflicts significant damage on various fish species annually, as susceptible populations consistently exhibit the presence of the virus (Kim et al., 2018; Lopez-Porras et al., 2018; Girisha et al., 2020; Puneeth et al., 2021).
Red sea bream is a rapidly growing fish known for its exceptional meat quality, making it highly sought after. As a result, it has become a significant aquaculture species in Korea, Japan, and China. Numerous studies have been conducted on red sea bream, focusing on growth and stocking density (Murashita et al., 2018; Ido et al., 2019; Seo et al., 2020; Mzengereza et al., 2021). However, among the existing research on immunity, only a few reports have specifically addressed IRF (Haque et al., 2021; Joo et al., 2022; Kim et al., 2023).
Since IRFs have been shown to play a role in viral pathogen infections, gaining insight into the features and expression of the IRF gene could substantially contribute to research on innate immunity and viral infections. Consequently, this study aimed to explore the molecular properties of IRF10 (belonging to the IRF-4G subfamily) in red sea bream and to characterize their expression during RSIV infection. This report provides the first set of data on IRF10 in relation to innate immunity and viral infection in sea bream.
Ⅱ. Materials and methods
1. Fish and virus
In this study, red sea bream were sourced from a net cage culture in Tongyeong, Gyeongsangnam-do, South Korea. The fish had an average weight of 52.1 ± 4.6 g and a total length of 12.5 ± 1.6 cm. Upon arrival at the laboratory, the fish were acclimatized with aeration in a flowing water tank at 22 ± 1 °C. To ensure the absence of pathogen infection, internal and external clinical symptoms, as well as molecular diagnostics, were employed on five randomly chosen fish. RSIV was isolated from rock bream and cultured in Tongyeong. Subsequent sequencing research confirmed that the virus belonged to genotype II, which is prevalent in Korea.
2. Molecular characteristics and phylogenetic analysis
To verify the integrity of the gene sequence obtained through genome sequencing, subcloning was conducted using Escherichia coli (E. coli) JM109 and the pGEM T-easy vector (Promega, USA). Sanger sequencing validated the sequence, which was then translated into amino acid sequences using GENETYX software version 8.0 (SDC Software Development, Japan). The National Center for Biotechnology Information (http://www.ncbi.nlm.nih.gov/blast) search method was employed to identify relevant sequences, while ClustalX 2.1 and MEGA 6.0 were utilized for multiple sequence alignment and phylogenetic analysis, respectively. The Expert Protein Analysis System PROSITE Scan program (http://prosite.expasy.org) was used to detect specific domains and motifs.
3. Tissue collection and cDNA synthesis
Tissue samples from five healthy red sea breams were aseptically collected and preserved at -80 °C to investigate the mRNA distribution of PmIRF10. The samples included head kidney, trunk kidney, skin, stomach, gill, heart, liver, spleen, eye, brain, and intestine. Fish were experimentally infected intraperitoneally with 1 × 107 copies/fish of RSIV suspended in PBS. These RSIV-infected red sea breams were employed to assess the temporal expression changes of PmIRF10 mRNA. Gill, spleen, liver, and trunk kidney samples were aseptically collected from five randomly selected red sea breams at 0 (control), 1, 3, 6, 12, 24, and 36 hours, and 3, 5, and 7 days post-RSIV infection. The easy-spin™ Total RNA Extraction Kit (iNtRON Biotechnology, South Korea) was used to extract total RNA from each tissue, while the PrimeScript™ 1st strand cDNA Synthesis Kit (Takara, Japan) was employed for cDNA synthesis.
4. Reverse transcription-quantitative real-time PCR (RT–qPCR) analysis
To determine the mRNA expression levels for each gene, RT-qPCR was performed using TB Green™ Premix Ex Taq™ (Takara). The sequences and details of the primer sets employed can be found in <Table 1>. The PCR reactions were carried out using a Thermal Cycler Dice® Real-Time System III (Takara). The 2−ΔΔCt method was applied to quantify the relative mRNA expression, utilizing EF-1α (AY190693.1) as the reference gene (Livak and Schmittgen, 2001). These experiments were conducted in triplicate.
5. Statistical analysis
The expression of PmIRF10 in RSIV-infected tissues was analyzed using one-way analysis of variance (ANOVA) followed by Tukey's HSD test, employing IBM SPSS software version 19 (IBM, USA). The standard deviation (SD) of the experimental results is denoted, and all samples were examined in triplicate.
Ⅲ. Results
1. PmIRF10 Gene Structure and Amino Acid Sequence Characterization
Genome sequencing of PmIRF10 revealed molecular properties, such as mRNA structure and amino acid (aa) sequences ([Fig. 1]). The predicted molecular weight of PmIRF10 (OK340063.1) was 46 kDa, with a theoretical pI of 8.79. The PmIRF10 sequence extended across 1218 bp (406 aa) in total, encompassing a DNA-binding domain (2-115 aa) and an IRF-association domain (203-385 aa). Moreover, PmIRF10 both possessed five tryptophan residues, and the nuclear localization signal (NLS) was highly conserved.
2. Multiple alignments and phylogenetic analysis of PmIRF10
Multiple alignment analyses of PmIRF10 revealed that the five tryptophan residues and the DBD and IAD domains were highly conserved across all species examined ([Fig. 2]). PmIRF10 exhibited the highest homology with large yellow croaker (86.67%) and orange-spotted grouper (82.34%) ([Fig. 2]). Phylogenetic analysis confirmed that PmIRF10 formed distinct clusters ([Fig. 3]).
3. Analysis of PmIRF10 mRNA expression in healthy red sea bream
Expression analysis of PmIRF10 in healthy red sea bream revealed that PmIRF10 demonstrated the highest expression in the gills as well (19.7-fold), followed by the spleen (10.2-fold), liver (5.2-fold), heart (4.7-fold), and skin (4.3-fold) ([Fig. 4]).
4. PmIRF10 mRNA expression in red sea bream following RSIV infection
Hourly expression analysis of PmIRF10 following RSIV infection revealed that PmIRF10 expression was observed to be downregulated in the gills until 3 hours post-infection (hpi), after initial expression, and then strongly upregulated at 6 hpi ([Fig. 5]).
Beyond this time point, expression persisted but decreased and exhibited low regulation at 7 days post-infection (dpi). In the spleen, upregulation was observed up to 6 hpi after initial expression, followed by downregulation until 36 hpi. The highest upregulation was observed at 5 dpi, succeeded by lower regulation at 7 dpi. In the liver, initial expression was followed by rapid upregulation at 3 hpi and 6 hpi, and then downregulation at 12 hpi. Up to 36 hpi, low regulation was observed, followed by upregulation at 3 dpi, and finally, low regulation at 7 dpi. After initial expression in the kidney, the strongest regulation occurred at 6 hpi, followed by downregulation, high upregulation at 5 dpi, and low regulation at 7 dpi.
Ⅳ. Discussion
IRFs are acknowledged for their antiviral and antibacterial properties in animals. These extensively investigated transcription factors play a role in responding to viral infections and regulating apoptosis or the development of immune responses (Savitsky et al., 2010). Via the DBD domain, IRFs bind to the ISRE as homodimers and heterodimers, interacting with the ISRE half-site on the opposite side of the DNA in a proximal direction, thus activating TLR and IFN signaling pathways (Andrilenas et al., 2018). Five tryptophan residues within the DBD domain are associated with the promoters of target genes, such as IFN regulatory elements and ISRE. Although there is an abundance of research on IRFs in mammals, studies in fish are relatively limited, and the variation in gene expression among different fish species is more distinct than in mammals. Consequently, it is crucial to investigate IRFs in individual fish species.
In this study, the nucleotide sequences of IRF10 from red sea bream were identified and displayed considerable similarity to previously reported fish sequences, suggesting that they are part of the IRF family. Furthermore, the DBD domain, an essential domain of IRF, was found to contain five tryptophan residues, and both the NLS sequence and the IAD domain were verified to be highly conserved.
Following the expression analysis of PmIRF10 after RSIV infection, it was observed that all tissues exhibited upregulation at 6 hpi, subsequently downregulation, and then upregulation again at 5 dpi in the spleen, liver, and kidney. In the case of orange-spotted grouper, upon examination of immune tissue expression of IRF10 after poly (I:C) treatment, the highest expression was detected at 6 h post-treatment. Similarly, after poly (I:C) treatment in common carp, spleen, head kidney, foregut, and hindgut displayed notable upregulation at 6 h post-treatment (Huang et al., 2015; Zhu et al., 2020). Moreover, IRF10 of Atlantic cod demonstrated the highest expression at 6 hpi following poly (I:C) treatment. In Japanese flounder, the highest expression was observed 6 h after poly (I:C) administration (Suzuki et al., 2011; Inkpen et al., 2015). The consistent early expression of IRF10 across multiple fish species implies that an antiviral immune response elicits prompt feedback from IRF10. The initial upregulation in PmIRF10, succeeded by a decrease and then an increase in expression pattern, may correlate with the finding that IRF10 in zebrafish acts as a negative regulator to maintain equilibrium in the antiviral immune response (Li et al., 2014). By integrating the outcomes of this study with those from prior research, it can be concluded that PmIRF10 assumes a unique role in the immune response of red sea bream, akin to the role of IRF10 in previously investigated fish, and it has been demonstrated to play a crucial function in antiviral immunity specifically.
In this study, high expression of PmIRF10 was confirmed in the immune tissues of the gills and spleen. In previous studies, normal tissue expression in rainbow trout, Asian swamp eel, and grass carp showed high levels of immune tissue-related expression, and in Japanese flounder IRF10 was strongly expressed in the gills, intestines, head kidney, and trunk kidneys (Suzuki et al., 2011; Xu et al., 2016). However, IRF10 in common carp showed different results, with strong expression in the gonad and the lowest expression in the head kidneys (Zhu et al., 2020). Although there are differences between species in normal tissue expression analysis, it is suggested that IRF10 in red sea bream is related to the immune response, similar to the previously reported IRF10 of other fish species.
We have verified the molecular characteristics and expression analysis of IRF10 in red sea bream. Although further investigation is required to ascertain their precise functions, PmIRF10 appear to exhibit distinct roles in the antiviral immune response, consistent with observations of IRF10 in previous studies. Consequently, this research provides a foundational dataset for understanding IRF characteristics in red sea bream and expands upon existing reports in teleost fish. Moreover, PmIRF10 are crucial immune genes in red sea bream, and additional immune-related studies are warranted to thoroughly elucidate their gene functions.
References
- Alter-Koltunoff M, Goren S, Nousbeck J, Feng CG, Sher A, Ozato K, Azriel A and Levi B-Z(2008). Innate immunity to intraphagosomal pathogens is mediated by interferon regulatory factor 8 (IRF-8) that stimulates the expression of macrophage-specific Nramp1 through antagonizing repression by c-Myc. Journal of Biological Chemistry, 283(5), 2724~2733. [https://doi.org/10.1074/jbc.M707704200]
- Andrilenas KK, Ramlall V, Kurland J, Leung B, Harbaugh AG and Siggers T(2018). DNA-binding landscape of IRF3, IRF5 and IRF7 dimers: implications for dimer-specific gene regulation. Nucleic acids research, 46(5), 2509~2520. [https://doi.org/10.1093/nar/gky002]
- Antonczyk A, Krist B, Sajek M, Michalska A, Piaszyk-Borychowska A, Plens-Galaska M, Wesoly J and Bluyssen HAR(2019). Direct Inhibition of IRF-Dependent Transcriptional Regulatory Mechanisms Associated With Disease. Front Immunol, 10, 1176. [https://doi.org/10.3389/fimmu.2019.01176]
- Gabriele L and Ozato K(2007). The role of the interferon regulatory factor (IRF) family in dendritic cell development and function. Cytokine & growth factor reviews, 18(5-6), 503~510. [https://doi.org/10.1016/j.cytogfr.2007.06.008]
- Girisha S, Puneeth T, Nithin M, Kumar BN, Ajay S, Vinay T, Suresh T, Venugopal M and Ramesh K(2020). Red sea bream iridovirus disease (RSIVD) outbreak in Asian seabass (Lates calcarifer) cultured in open estuarine cages along the west coast of India: First report. Aquaculture, 520, 734712. [https://doi.org/10.1016/j.aquaculture.2019.734712]
- Guan Y, Chen X, Luo T, Ao J, Ai C and Chen X (2020). Molecular characterization of the interferon regulatory factor (IRF) family and functional analysis of IRF11 in the large yellow croaker (Larimichthys crocea). Fish & Shellfish Immunology, 107, 218~229. [https://doi.org/10.1016/j.fsi.2020.10.001]
- Haque MN, Nam S-E, Shin YK and Rhee J-S(2021). The dinoflagellate Alexandrium affine acutely induces significant modulations on innate immunity, hepatic function, and antioxidant defense system in the gill and liver tissues of red seabream. Aquatic Toxicology, 240, 105985. [https://doi.org/10.1016/j.aquatox.2021.105985]
- Huang B, Jia QQ, Liang Y, Huang WS and Nie P(2015). Interferon regulatory factor 10 (IRF10): cloning in orange spotted grouper, Epinephelus coioides, and evolutionary analysis in vertebrates. Fish & Shellfish Immunology, 46(2), 669~677. [https://doi.org/10.1016/j.fsi.2015.08.003]
- Huang B, Qi ZT, Xu Z and Nie, P(2010). Global characterization of interferon regulatory factor (IRF) genes in vertebrates: glimpse of the diversification in evolution. BMC immunology, 11(1), 1~28. [https://doi.org/10.1186/1471-2172-11-22]
- Ido A, Hashizume A, Ohta T, Takahashi T, Miura C and Miura T(2019). Replacement of fish meal by defatted yellow mealworm (Tenebrio molitor) larvae in diet improves growth performance and disease resistance in red seabream (Pargus major). Animals, 9(3), 100. [https://doi.org/10.3390/ani9030100]
- Inkpen SM, Hori TS, Gamperl AK, Nash GW and Rise ML(2015). Characterization and expression analyses of five interferon regulatory factor transcripts (Irf4a, Irf4b, Irf7, Irf8, Irf10) in Atlantic cod (Gadus morhua). Fish & Shellfish Immunology, 44(1), 365~381. [https://doi.org/10.1016/j.fsi.2015.02.032]
- Inouye K, Yamano K, Maeno Y, Nakajima K, Matsuoka M, Wada Y and Sorimachi M(1992). Iridovirus infection of cultured red sea bream, Pagrus major. Fish Pathology, 27(1), 19~27. [https://doi.org/10.3147/jsfp.27.19]
- Joo M-S, Choi K-M, Kang G, Woo W-S, Kim K-H, Sohn M-Y, Son H-J, Lee J-H and Park C-I (2022). Molecular characteristics, expression analysis and antiviral activity against RSIV of red sea bream IRF3 and IRF7. Aquaculture Reports, 27, 101415. [https://doi.org/10.1016/j.aqrep.2022.101415]
- Kim K-H, Joo M-S, Kang G, Woo W-S, Sohn M-Y, Son H-J, and Park C-I(2023). Characterization of Red Sea Bream (Pagrus major) Interferon Regulatory Factor 5 and 6 Genes and Their Expression in Response to RSIV Infection. Fishes, 8(2), 114. [https://doi.org/10.3390/fishes8020114]
- Kim K, Hwang S, Cho M, Jung S, Kim Y and Jeong H(2018). A natural infection by the red sea bream iridovirus‐type Megalocytivirus in the golden mandarin fish Siniperca scherzeri. Journal of fish diseases, 41(8), 1229~1233. [https://doi.org/10.1111/jfd.12815]
- Li S, Feng H, Wu N, Zou J, Zhang Y and Zhang Y-A(2013). IRF10 is a negative regulator of IFN transcription in zebrafish. 中國水產學會魚病專業委員會 2013 年學術研討會論文摘要匯編. [https://doi.org/10.4049/jimmunol.1400253]
- Li S, Lu L-F, Feng H, Wu N, Chen D-D, Zhang Y-B, Gui J-F, Nie P and Zhang Y-A(2014). IFN regulatory factor 10 is a negative regulator of the IFN responses in fish. The Journal of Immunology, 193(3), 1100~1109. [https://doi.org/10.4049/jimmunol.1400253]
- Livak KJ and Schmittgen TD(2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2− ΔΔCT method. Methods, 25(4), 402~408. [https://doi.org/10.1006/meth.2001.1262]
- Lopez-Porras A, Morales JA, Alvarado G, Koda SA, Camus A, Subramaniam K, Waltzek TB and Soto E(2018). Red seabream iridovirus associated with cultured Florida pompano Trachinotus carolinus mortality in Central America. Diseases of aquatic organisms, 130(2), 109~115. [https://doi.org/10.3354/dao03267]
- Medzhitov R, and Janeway Jr, C(2000). The Toll receptor family and microbial recognition. Trends in microbiology, 8(10), 452~456. [https://doi.org/10.1016/S0966-842X(00)01845-X]
- Murashita K, Matsunari H, Furuita H, Rønnestad I, Oku H and Yamamoto T(2018). Effects of dietary soybean meal on the digestive physiology of red seabream Pagrus major. Aquaculture, 493, 219~228. [https://doi.org/10.1016/j.aquaculture.2018.05.005]
- Mzengereza K, Ishikawa M, Koshio S, Yokoyama S, Yukun Z, Shadrack RS, Seo S, Duy Khoa TN, Moss A and Dossou S(2021). Effect of substituting fish oil with camelina oil on growth performance, fatty acid profile, digestibility, liver histology, and antioxidative status of red seabream (Pagrus major). Animals, 11(7), 1990. [https://doi.org/10.3390/ani11071990]
- Nehyba J, Hrdličkova, R, Burnside J and Bose Jr, HR(2002). A novel interferon regulatory factor (IRF), IRF-10, has a unique role in immune defense and is induced by the v-Rel oncoprotein. Molecular and cellular biology, 22(11), 3942~3957. [https://doi.org/10.1128/MCB.22.11.3942-3957.2002]
- Paun A and Pitha P(2007). The IRF family, revisited. Biochimie, 89(6-7), 744~753. [https://doi.org/10.1016/j.biochi.2007.01.014]
- Puneeth T, Baliga P, Girisha S, Shekar M, Nithin M, Suresh T and Kumar BN(2021). Complete genome analysis of a red seabream iridovirus (RSIV) isolated from Asian seabass (Lates calcarifer) in India. Virus research, 291, 198199. [https://doi.org/10.1016/j.virusres.2020.198199]
- Qiao-Qing X, Dai-Qin Y, Rui T, Jing W and Ming-Xian CHANG PN(2014). Gene cloning and induced expression pattern of IRF4 and IRF10 in the Asian swamp eel (Monopterus albus). Zoological Research, 35(5), 380. [https://doi.org/10.13918/j.issn.2095-8137.2014.5.380]
- Savitsky D, Tamura T, Yanai H and Taniguchi T(2010). Regulation of immunity and oncogenesis by the IRF transcription factor family. Cancer immunology, immunotherapy, 59(4), 489~510. [https://doi.org/10.1007/s00262-009-0804-6]
- Seo JS, Haque MN, Nam S-E, Kim B-M, and Rhee J-S(2020). Inorganic nitrogen compounds reduce immunity and induce oxidative stress in red seabream. Fish & Shellfish Immunology, 104, 237~244. [https://doi.org/10.1016/j.fsi.2020.05.072]
- Shu C, Sun Y and Xu T(2015). Molecular characterization of three IRF1 subfamily members reveals evolutionary significance of IRF11 in miiuy croaker. Developmental & Comparative Immunology, 53(2), 385~391. [https://doi.org/10.1016/j.dci.2015.07.009]
- Suzuki Y, Yasuike M, Kondo H, Aoki T and Hirono I(2011). Molecular cloning and expression analysis of interferon regulatory factor 10 (IRF10) in Japanese flounder, Paralichthys olivaceus. Fish & Shellfish Immunology, 30(1), 67~76. [https://doi.org/10.1016/j.fsi.2010.09.010]
- Xu Q, Jiang Y, Wangkahart E, Zou J, Chang M, Yang D, Secombes CJ, Nie P and Wang T(2016). Sequence and expression analysis of interferon regulatory factor 10 (IRF10) in three diverse teleost fish reveals its role in antiviral defense. PloS one, 11(1), e0147181. [https://doi.org/10.1371/journal.pone.0147181]
- Zhang R, Chen K, Peng L and Xiong H(2012). Regulation of T helper cell differentiation by interferon regulatory factor family members. Immunologic research, 54(1), 169~176. [https://doi.org/10.1007/s12026-012-8328-0]
- Zhu Y, Shan S, Zhao H, Liu R, Wang H, Chen X, Yang G and Li H(2020). Identification and Functional Analysis of A Non-Mammalian IRF10 in Common Carp (Cyprinus Carpio. L) in The Antiviral and Antibacterial Immune Response. [https://doi.org/10.21203/rs.3.rs-63536/v1]