- Research article
- Open Access
African swine fever virus does not express viral microRNAs in experimentally infected pigs
BMC Veterinary Researchvolume 14, Article number: 268 (2018)
African swine fever virus (ASFV) is the etiological agent of African swine fever (ASF), a re-expanding devastating and highly lethal hemorrhagic viral disease. microRNAs (miRNAs) are a new class of small non-coding RNAs that regulate gene expression post-transcriptionally. The discovery of virus specific miRNAs has increased both in number and importance in the past few years. We have recently described the differential expression of several porcine miRNAs during in vivo infection with attenuated and virulent ASFV strains. Here, we have extended these studies trying to identify the presence of viral miRNAs encoded by ASFV in an in vivo infection in pigs.
Sixteen small RNA libraries were analyzed from spleen and submandibular lymph nodes obtained from eight pigs, seven infected with either the virulent E75 ASFV strain or its attenuated counterpart E75CV1, or from pigs surviving E75CV1-infection and challenged with BA71 (heterologous challenge) and one non infected as negative control. Samples were recovered at different times post-infection. Libraries were analyzed by next-generation sequencing. Some viral miRNA candidates were initially identified, which did not correspond to porcine miRNAs. Further structural analyses were carried out in order to confirm if they met the conformational requirements to be considered a viral miRNA.
The analysis of sixteen small RNA libraries prepared from two different tissues obtained from pigs experimentally infected with E75, E75CV1 or with E75CV1 plus BA71, revealed the presence of six potential miRNA sequences but none of them met the requirements to be considered as viral miRNAs. Thus, we can conclude that ASFV does not express miRNAs in vivo, at least under the experimental conditions described here.
African swine fever (ASF) is a contagious disease that affects domestic and wild pigs caused by ASF virus (ASFV) belonging to the Asfarviridae family, genus Asfivirus. ASFV is a large icosahedral, enveloped, double stranded (ds) DNA virus ranging from 170 to 193 kb . ASFV is today considered one of the main biological threats for pigs world-wide, causing enormous economical loss in affected countries. Endemic since its origin in several areas of sub-Saharan Africa, ASFV re-entered the European continent in 2007, expanding to surrounding areas and reaching the European Union (EU) borders in 2014. Since then, the virus reached some Eastern European Union countries and in 2017, ASFV was reported in Czech Republic and Romania . Lack of vaccines or treatments complicates its control , therefore gaining knowledge about ASFV pathogenesis is valuable for the future development of anti-ASFV strategies. The viral particle, with a diameter of 200 nm, has three concentric wraps with an additional external membrane obtained when coming out by budding from the infected cell. Between 150 and 175 ORFs have been identified but the functions for half of them still remain unknown. Sixteen different genomes have been identified and fully sequenced from pigs, wild boars and ticks and 23 genotypes have been so far described . The virus infects monocytes and macrophages, key cell players for the immune system. In addition, it has been observed that virus replicates in endothelial cells, hepatocytes, renal tubular epithelial cells and neutrophils, mainly at late stages of an acute ASFV-infection. miRNAs are a new class of small non-coding RNAs with posttranscriptional regulation functions of gene expression. miRNAs play an increasing role in many biological processes . The study of miRNA-mediated host-pathogen interactions has emerged in the last decade mainly due to the important role that miRNAs play in antiviral defense . In addition, miRNAs encoded by viruses have been identified to be involved in the infection process. To date, the microRNA database, miRBase , shows 308 precursors and 502 mature viral miRNAs have been identified. Viruses encoding miRNAs has been identified in several families, as Herpesviridae, Polyomaviridae, Adenoviridae, Papillomaviridae, Baculoviridae, Ascoviridae and Iridoviridae [8,9,10]. All of them are DNA viruses, with an essential nuclear phase in their replication cycle, necessary for canonical miRNA biogenesis. In addition, Ascovirus and Iridovirus, together with Asfivirus, belong to the monophyletic superfamily of nucleocytoplasmic large DNA viruses (NCLDV). Thus, ASFV was a good candidate for a virus capable of encoding miRNAs. On the other hand, recently, our group has identified differential expression of porcine miRNAs in pigs after the infection with ASFV .
Considering this background, the aim of the present study was to explore if ASFV can encode viral miRNAs in an experimental infection in its natural host by using high-throughput sequencing. To this end, in order to characterize the potential presence of ASFV-specific miRNAs, in silico prediction was carried out to determine if the ASFV genomes of two different ASFV strains with different degrees of virulence, encode possible miRNA precursors. Subsequently, to explore whether the possible miRNAs candidates are present in the ASFV, small RNA libraries from target tissues were constructed and analyzed.
In silico prediction
The Vmir prediction algorithm  was used to predict the presence of hairpin structures in the ASFV genomes from: strain E75 (Genbank accession number FN557520)  and BA71 (Genbank accession number KP055815) . The sequence of the non-pathogenic BA71V (Genbank accession number NC_001659) , highly adapted to Vero cells, was additionally included in the study. Computational prediction was conducted by using Vmir of viral miRNAs, setting as conditions a hairpin size from 50 to 100 nt and a window count of 15.
In order to explore whether these candidates were present in the ASFV genome, samples from spleen and submandibular lymph node (SLN) from seven 8-week-old Landrace X Large White male pigs were used. The animals were acquired from Hypor (A Coruña, Spain), and were fed ad libitum with conventional compound feed for growing pigs. The animals were intramuscularly inoculated in the neck with the virulent ASFV strains: BA71 or E75 or with the attenuated E75CV1, obtained by adaptation of E75 after three passages in Vero cell lineage. The same dose, 104 hemadsorbing units (HAU50), was used for each virus. One non-infected pig was included as negative control. All samples used are summarized in Table 1, noticing that some of the samples were previously used to successfully identify the expression of cellular miRNAs during the ASFV infection in vivo . The potential presence of viral miRNAs in ASFV-infected pigs was analyzed by high-throughput sequencing of small RNA libraries prepared from two tissues recovered in three conditions: i) The animals infected with virulent strain E75 necropsied at 3 and 7 dpi, ii) The animals infected with the attenuated strain E75CV1 and necropsied at 3 and 31 dpi, iii) The animals infected with the attenuated strain E75CV1 and re-inoculated at 31 dpi with the virulent strain BA71 and necropsied at 38 dpi. As negative control, small RNA libraries were sequenced from a non-infected animal at 0 dpi.
The time points chosen for each experimental infection were selected taking into account previous results  and covering acute and chronic phases of infection with different ASFV strains. Euthanasia was performed according to European Directive 2010/63/EU, using a pentobarbital overdose of 60–100 mg/kg administered via the anterior vena cava.
All animal experiments were carried out in the BSL3 facilities at CReSA, and all animal managements were performed in accordance with the guidelines of the Good Experimental Practices. All procedures were approved by the Ethical and Animal Welfare Committee of the Universitat Autònoma de Barcelona with the permit number DMAH-5796.
SYBR Green quantitative PCR (qPCR) was used with the following oligonucleotides: pF1 primer (5′ CCTCGGCGAGCGCTTTATCAC 3′) and pR2 (5′ GGAAACTCATTCACCAAATCCTT 3′), designed in a highly conserved region within the p72 ORF . qPCR mix contained 480 nM each primer, 12.5 μl of Power SYBR®Green PCR Master Mix (Thermo Fisher Scientific) and 2.5 μl of template. Nanopure autoclaved water was added to final volume of 25 μL. Amplification were carried out with these conditions: 10 min at 95 °C, 2 min at 50 °C and 40 cycles of 15 s at 95 °C and 1 min at 60 °C. Triplicates of each sample were used for the qPCR.
Small RNA library construction
Total RNA extractions were carried out in order to construct small RNA libraries from spleen and SLN samples from the animals infected with E75CV1 necropsied at 31 dpi, the E75CV1-infected animals and re-inoculated at 31 dpi with BA71, necropsied at 38 dpi, and from one non-infected animal as described in  with some modifications. A total of 8 small RNA libraries were constructed in a two-step ligation procedure with the 3′ and 5′ adaptors from IDT technologies. Briefly, reverse transcription (RT) was carried out with primer mir5 and PCR amplification with primers A5 and B3 (which included multiplex identifiers at the 5′ end) containing sequences complementary to the 3′ and 5′ adaptors used for miRNA library construction. Ion Torrent adapters were ligated to PCR products and libraries were then amplified with Ion Torrent primers. Libraries were sequenced following the manufacturer’s protocol with Ion PGM™ sequencer (Thermo Fisher Scientific, Massachusetts, USA) at DNA sequencing facilities at CRAG (Centre de Recerca en Agrigenòmica, Universitat Autònoma de Barcelona, Spain). On the other hand, cellular miRNAs expression has been previously detected using some of these libraries thus ensuring that they are valid for the present study.
Sequence processing scheme and analysis
Strict quality parameters of the individual nucleotide sequence reads has been applied and only those with the Phread base quality scores of Q > 30 in all its bases have been used. Primer sequences were removed and only those insert sequences between 15 and 29 nucleotides were used for further analysis. Sequences aligned with Sus scrofa genome were analyzed in  since they are porcine miRNAs. Local Blast against ASFV genomes was carried out and Mfold  was used in order to evaluate the secondary structure and free energy for the potential miRNA candidates. Likewise, candidates were blasted to the output of Vmir hairpin ASFV structures. In addition, candidates were checked with the algorithm developed in miREval 2.0 for detecting novel miRNAs .
The in silico prediction allowed the identification of 415 miRNA candidates in the E75 strain with a score higher than 100 (Fig. 1), in the case of BA71, 391 candidates with a score higher than 100 were identified (Fig. 2) and in the case of BA71V, 366 candidates were identified in the same conditions (data not shown).
In the animal infection, the pigs inoculated with the virulent strains E75 and BA71 developed typical clinical signs, while the pigs inoculated with the attenuated E75CV1 strain developed mild or no clinical signs . Quantitative PCR was performed from the samples of spleen and SLN from all the animals. Those animals inoculated with E75 showed low Ct values at 7 dpi coinciding with the peak of virus replication and ASF clinical signs, while in contrast, inoculation of the same dose of the homologous attenuated E75CV1 strain showed no amplification. The animals infected with E75CV1 and re-inoculated with the BA71 virulent strain showed amplification levels similar to those found after E75 infection, correlating with the lack of cross-protection observed for these two strains . As expected, the negative control animal did not show any specific PCR amplification (Table 2).
Total RNAs from SLN and spleen were used to construct small RNA libraries from the selected animals. High-throughput sequencing provided 272,967 reads from the eight new libraries that, in addition to eight previously sequenced libraries with 301,802 reads, provided a total of 574,769 reads from the sixteen libraries (Table 3). New sequencing data were submitted to ENA database  with the accession numbers ERX1811688–95. When trimmed to filter sequences ranging from 15 to 29 nucleotides, a total of 169,511 reads were obtained. After the identification of host miRNAs by comparing with miRBase , sequences were blasted to ASFV. Reads were blasted to the E75 and BA71 isolates considering only sequences with 100% of alignment and identity rendering 9 hits with the ASFV genome (Table 3). An increased number of mismatches in the extremes due to miRNA variability (isomiRs)  were allowed. The search rendered 31 hits. This processing allowed the identification of 40 viral genome hits with 8 unique sequences when allowing 100% homology (Tables 3 and 4) and 23 unique sequences when allowing < 100% homology (Tables 3 and 5). From these sequences, 25 candidates were obtained when grouped by isomiRs and 6 candidates when copy number (CN) was ≥2.
Sequences were blasted to the output of Vmir hairpins from E75, BA71 and BA71V. In addition, candidates were evaluated with miREval 2.0 software with no positive results. Therefore, after the analysis of the complementarity of the candidates with the Vmir prediction and the checking with miREval, together with the analysis with Mfold of their secondary structure and the free minimum energy of their folding, none of the candidates could be considered as a viral miRNA as they failed to fulfill the necessary requirements.
In a previous study carried out by our group, it was described that there is in vivo upregulation and downregulation of several porcine miRNAs after infection with ASFV . Using the small RNA libraries obtained for this purpose plus new small RNA libraries constructed in this work, we have studied the capability of E75, BA71 and E75CV1 isolates of ASFV to express viral miRNAs in vivo by using high-throughput sequencing.
The genomic characteristics of ASFV fit with the requirements of a virus to be able to encode viral miRNAs because it is a DNA virus with a nuclear phase [22, 23] and the in silico prediction carried out in this work found several potential hairpins to be virus- encoded miRNAs. In addition, 16 viral miRNAs have been identified in Singapore grouper iridovirus (SGIV)  and in tiger frog virus (TFV)  both members of the Iridoviridae family. Likewise, virus encoded miRNAs have been identified in ascovirus which infects mosquitoes . In this last case, a viral miRNA reduces DNA polymerase I levels by transcriptional degradation and modulates the ascovirus replication. All these viruses belong, as the asfivirus, to the monophyletic superfamily of nucleocytoplasmic large DNA viruses (NCLDV), reinforcing that ASFV is a good candidate to explore virus encoded miRNAs.
There is no clear knowledge about when viral miRNAs are expressed at highest levels after infection according to Cullen , who showed the variability of the viral miRNAs expression with time among different viruses. Therefore, we have analyzed samples from the infected animals at different times points post-infection in addition to different viruses with different virulence and the animals infected with an attenuated strain and re-inoculated with a heterologous virulent strain in order to increase the spectrum of conditions. Besides, although E75 and BA71 are classified in the West African-European cluster, the percentage of identity at genome level is 88.5%, the lower identity among viruses included in this cluster . ASFV produces lytic infection. Nevertheless, although most virus-encoded miRNAs have been related to persistent infections, viral miRNAs have also been detected during lytic infection [9, 17].
After the identification of porcine miRNAs expressed in ASFV-infected animals , high-throughput sequencing revealed other potential miRNAs sharing identity with the viral genome. Out of a total of 40, six with a copy number ≥ 2 were selected for further analysis. These candidates were blasted against the Vmir hairpins prediction from E75 and BA71 genomes. Also, BA71V, a highly adapted strain in Vero cells, was included in the analysis in order to avoid the loss of candidates because of the possibility of genome variability in cell culture passages . However, it was observed that none of the candidates matched any of the predicted hairpins, had appropriate secondary structure, or minimum free energy, showing that the candidates did not have the necessary characteristics to be identified as viral miRNAs. The candidate ASFV-C9, which presents the higher CN (9), mapped into the early structural protein p22 . p22 is a target gene of miR-30e-5p, a differentially expressed (DE) miRNA described in our previous work . This would explain why a higher number of sequences have been found. On the other hand, all isomiRs were blasted to miRBase, and none of them were identified as host miRNA. Therefore, ASFV-C9 with all other sequences can be considered as degradation products of viral RNA. Being aware of the reduced number of animals used, different tissues were used to generate a total of 14 RNA libraries, number that exceeds the majority of similar studies so far published, normally ranging from 1 to 4 and carried out in cell culture, and not in the natural host [29,30,31,32].
Although DNA viruses have the potential to encode viral miRNAs, this characteristic does not imply the viral miRNAs expression of these viruses. For example, the herpesvirus family have the highest viral miRNAs encoding rates, and almost each member is capable to codify miRNAs. However, there are members of this family as the varicella zoster virus (VZV) and other human neurotropic herpesviruses (HSV1 and HSV2) , which do not encode viral miRNAs. The knowledge of why some viruses encode miRNAs and why some phylogenetically related viruses do not, can be of great relevance for the fields of virology and RNAi .
In this study, we have found that ASFV, under these experimental conditions, does not express viral miRNAs like other viruses such as Cowpox virus or Porcine circovirus type 2 (PCV2) [34, 35]. Interestingly, in a recent and exhaustive work carried out by Jaing et al. about gene expression in blood of ASFV-infected pigs, wherein some DE host miRNA have been identified, viral miRNAs are not described . In the literature, there are examples in which were indicated that certain viruses encode miRNAs that were initially considered as non coding miRNAs, e.g. human papillomavirus . Therefore, studies with higher number of reads, together with analysis of the presence of viral miRNAs in cell culture and other in vivo conditions should be considered before fully discounting if ASFV expresses miRNAs.
In the present work, we can conclude that African swine fever virus does not express viral miRNAs in the in vivo conditions studied here. This finding is contrary to what happens with other large DNA viruses and can be relevant to the biology of this important pathogen.
- ASF :
African swine fever
- ASFV :
African swine fever virus
- CN :
- DE :
- DPI :
- ENA :
European nucleotide database
- HTS :
- NCLDV :
Nucleocytoplasmic large DNA virus
- OIE :
World Organization for Animal Health
- SLN :
Submandibular lymph node
de Villiers EP, Gallardo C, Arias M, da Silva M, Upton C, Martin R, Bishop RP. Phylogenomic analysis of 11 complete African swine fever virus genome sequences. Virology. 2010;400(1):128–36.
World Organisation for Animal Health, OIE [http://www.oie.int].
Arias M, de la Torre A, Dixon L, Gallardo C, Jori F, Laddomada A, Martins C, Parkhouse RM, Revilla Y, Rodriguez FAJ, et al. Approaches and perspectives for development of African swine fever virus vaccines. Vaccines (Basel). 2017;5(4):e35.
Achenbach JE, Gallardo C, Nieto-Pelegrin E, Rivera-Arroyo B, Degefa-Negi T, Arias M, Jenberie S, Mulisa DD, Gizaw D, Gelaye E, et al. Identification of a new genotype of African swine fever virus in domestic pigs from Ethiopia. Transbound Emerg Dis. 2017;64(5):1393–404.
Tufekci KU, Meuwissen RL, Genc S. The role of microRNAs in biological processes. Methods Mol Biol. 2014;1107:15–31.
Umbach JL, Cullen BR. The role of RNAi and microRNAs in animal virus replication and antiviral immunity. Genes Dev. 2009;23(10):1151–64.
Griffiths-Jones S, Grocock RJ, van Dongen S, Bateman A, Enright AJ. miRBase: microRNA sequences, targets and gene nomenclature. Nucleic Acids Res. 2006;34(Database issue):D140–4.
Grundhoff A, Sullivan CS, microRNAs V-e. Virology. 2011;411(2):325–43.
Sullivan CS, Grundhoff AT, Tevethia S, Pipas JM, Ganem D. SV40-encoded microRNAs regulate viral gene expression and reduce susceptibility to cytotoxic T cells. Nature. 2005;435(7042):682–6.
Kincaid RP, Sullivan CS. Virus-encoded microRNAs: an overview and a look to the future. PLoS Pathog. 2012;8(12):e1003018.
Nunez-Hernandez F, Perez LJ, Munoz M, Vera G, Accensi F, Sanchez A, Rodriguez F, Nunez JI. Differential expression of porcine microRNAs in African swine fever virus infected pigs: a proof-of-concept study. Virol J. 2017;14(1):198.
Grundhoff A. Computational prediction of viral miRNAs. Methods Mol Biol. 2011;721:143–52.
Yanez RJ, Rodriguez JM, Nogal ML, Yuste L, Enriquez C, Rodriguez JF, Vinuela E. Analysis of the complete nucleotide sequence of African swine fever virus. Virology. 1995;208(1):249–78.
Rodríguez JM, Moreno LT, Alejo A, Lacasta A, Rodríguez F, Salas ML. Genome sequence of African swine fever virus BA71, the virulent parental strain of the nonpathogenic and tissue-culture adapted BA71V. PLoS One. 2015;10(11):e0142889.
Lacasta A, Monteagudo PL, Jimenez-Marin A, Accensi F, Ballester M, Argilaguet J, Galindo-Cardiel I, Segales J, Salas ML, Dominguez J, et al. Live attenuated African swine fever viruses as ideal tools to dissect the mechanisms involved in viral pathogenesis and immune protection. Vet Res. 2015;46:135.
Zsak L, Borca MV, Risatti GR, Zsak A, French RA, Lu Z, Kutish GF, Neilan JG, Callahan JD, Nelson WM, et al. Preclinical diagnosis of African swine fever in contact-exposed swine by a real-time PCR assay. J Clin Microbiol. 2005;43(1):112–9.
Timoneda O, Nunez-Hernandez F, Balcells I, Munoz M, Castello A, Vera G, Perez LJ, Egea R, Mir G, Cordoba S, et al. The role of viral and host microRNAs in the Aujeszky's disease virus during the infection process. PLoS One. 2014;9(1):e86965.
Zuker M. Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003;31(13):3406–15.
Gao D, Middleton R, Rasko JE, Ritchie W. miREval 2.0: a web tool for simple microRNA prediction in genome sequences. Bioinformatics. 2013;29(24):3225–6.
European Nucleotide Database [http://www.ebi.ac.uk/ena].
Neilsen CT, Goodall GJ, Bracken CP. IsomiRs--the overlooked repertoire in the dynamic microRNAome. Trends Genet. 2012;28(11):544–9.
Van Etten J. Lesser known large dsDNA viruses. Preface. Curr Top Microbiol Immunol. 2009;328:v–vii.
Ballester M, Rodríguez-Cariño C, Pérez M, Gallardo C, Rodríguez JM, Salas ML, Rodriguez F. Disruption of nuclear organization during the initial phase of African swine fever virus infection. J Virol. 2011;85(16):8263–9.
Yan Y, Cui H, Jiang S, Huang Y, Huang X, Wei S, Xu W, Qin Q. Identification of a novel marine fish virus, Singapore grouper iridovirus-encoded microRNAs expressed in grouper cells by Solexa sequencing. PLoS One. 2011;6(4):e19148.
Yuan J-M, Chen Y-S, He J, Weng S-P, Guo C-J, He J-G. Identification and differential expression analysis of MicroRNAs encoded by Tiger frog virus in cross-species infection in vitro. Virology J. 2016;13:73.
Hussain M, Taft RJ, Asgari S. An insect virus-encoded microRNA regulates viral replication. J Virol. 2008;82(18):9164–70.
Cullen BR. How do viruses avoid inhibition by endogenous cellular microRNAs? PLoS Pathog. 2013;9(11):e1003694.
Camacho A, Vinuela E. Protein p22 of African swine fever virus: an early structural protein that is incorporated into the membrane of infected cells. Virology. 1991;181(1):251–7.
Anselmo A, Flori L, Jaffrezic F, Rutigliano T, Cecere M, Cortes-Perez N, Lefevre F, Rogel-Gaillard C, Giuffra E. Co-expression of host and viral microRNAs in porcine dendritic cells infected by the pseudorabies virus. PLoS One. 2011;6(3):e17374.
Glazov EA, Horwood PF, Assavalapsakul W, Kongsuwan K, Mitchell RW, Mitter N, Mahony TJ. Characterization of microRNAs encoded by the bovine herpesvirus 1 genome. J Gen Virol. 2010;91(Pt 1):32–41.
Burnside J, Bernberg E, Anderson A, Lu C, Meyers BC, Green PJ, Jain N, Isaacs G, Morgan RW. Marek's disease virus encodes MicroRNAs that map to meq and the latency-associated transcript. J Virol. 2006;80(17):8778–86.
Levy CS, Hopkins J, Russell GC, Dalziel RG. Novel virus-encoded microRNA molecules expressed by ovine herpesvirus 2-immortalized bovine T-cells. J Gen Virol. 2012;93(Pt 1):150–4.
Umbach JL, Nagel MA, Cohrs RJ, Gilden DH, Cullen BR. Analysis of human alphaherpesvirus microRNA expression in latently infected human trigeminal ganglia. J Virol. 2009;83(20):10677–83.
Skalsky RL, Cullen BR. Viruses, microRNAs, and host interactions. Annu Rev Microbiol. 2010;64:123–41.
Nunez-Hernandez F, Perez LJ, Vera G, Cordoba S, Segales J, Sanchez A, Nunez JI. Evaluation of the capability of the PCV2 genome to encode miRNAs: lack of viral miRNA expression in an experimental infection. Vet Res. 2015;46:48.
Jaing C, Rowland RRR, Allen JE, Certoma A, Thissen JB, Bingham J, Rowe B, White JR, Wynne JW, Johnson D, et al. Gene expression analysis of whole blood RNA from pigs infected with low and high pathogenic African swine fever viruses. Sci Rep. 2017;7(1):10115.
Qian K, Pietila T, Ronty M, Michon F, Frilander MJ, Ritari J, Tarkkanen J, Paulin L, Auvinen P, Auvinen E. Identification and validation of human papillomavirus encoded microRNAs. PLoS One. 2013;8(7):e70202.
Authors thank all the personnel from the animal facilities at CReSA. IRTA is supported by CERCA Programme / Generalitat de Catalunya.
This work was supported by the projects AGL2010–22358-C02 and AGL2013–48998-C2–1-R from MINECO (Spanish government). FN-H is the recipient of a PhD fellowship FI-DGR (AGAUR) from Generalitat de Catalunya.
Availability of data and materials
Sequencing data generated in this work were deposited at European Nucleotide Archive (ENA, http://www.ebi.ac.uk/ena/) with the accession numbers ERX1811688–95.
The animal experiments were carried out in the BSL3 facilities at CReSA, and animal managements were performed in accordance with the guidelines of the Good Experimental Practices. All procedures were approved by the Ethical and Animal Welfare Committee of the Universitat Autònoma de Barcelona under the permit number DMAH-5796, according to national and European regulations.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.