- Open Access
- Open Peer Review
Human Endogenous Retrovirus K(HML-2) Gag and Env specific T-cell responses are not detected in HTLV-I-infected subjects using standard peptide screening methods
Journal of Negative Results in BioMedicinevolume 12, Article number: 3 (2013)
An estimated 10–20 million individuals are infected with the retrovirus human T-cell leukemia virus type 1 (HTLV-1). While the majority of these individuals remain asymptomatic, 0.3-4% develop a neurodegenerative inflammatory disease, termed HTLV-1-associated myelopathy/tropical spastic paraparesis (HAM/TSP). HAM/TSP results in the progressive demyelination of the central nervous system and is a differential diagnosis of multiple sclerosis (MS). The etiology of HAM/TSP is unclear, but evidence points to a role for CNS-inflitrating T-cells in pathogenesis. Recently, the HTLV-1-Tax protein has been shown to induce transcription of the human endogenous retrovirus (HERV) families W, H and K. Intriguingly, numerous studies have implicated these same HERV families in MS, though this association remains controversial.
Here, we explore the hypothesis that HTLV-1-infection results in the induction of HERV antigen expression and the elicitation of HERV-specific T-cells responses which, in turn, may be reactive against neurons and other tissues. PBMC from 15 HTLV-1-infected subjects, 5 of whom presented with HAM/TSP, were comprehensively screened for T-cell responses to overlapping peptides spanning HERV-K(HML-2) Gag and Env. In addition, we screened for responses to peptides derived from diverse HERV families, selected based on predicted binding to predicted optimal epitopes. We observed a lack of responses to each of these peptide sets.
Thus, although the limited scope of our screening prevents us from conclusively disproving our hypothesis, the current study does not provide data supporting a role for HERV-specific T-cell responses in HTLV-1 associated immunopathology.
Human T-cell leukemia virus type 1 (HTLV-1) is a deltaretrovirus that has infected humans for thousands of years. An estimated 10–20 million individuals worldwide are infected with HTLV-1, the majority of whom will remain asymptomatic for life. For a minority of infected individuals, however, HTLV-1 infection is associated with severe diseases, with 1-5% developing adult T-cell leukemia/lymphoma (ALT) and an additional 0.3-4% developing immune-mediated inflammatory disease involving the central nervous system, the eyes, the lungs and/or the skeletal muscles. Inflammation of the central nervous system (CNS) in HTLV-1-infected subjects results in progressive spasticity and muscle weakness in lower extremities and is termed HTLV-1-associated myelopathy/tropical spastic paraparesis (HAM/TSP). The etiology of HAM/TSP is poorly understood, however it is associated with higher HTLV-1 proviral loads, a higher frequency of HTLV-1-specific CD8+ T-cells and greater production of IFN-γ and TNF-α than that observed in asymptomatic HTLV-1-infected subjects, suggesting the involvement of cellular immune responses in the pathogenesis[5–8]. As there is no evidence that HTLV-1 infects astrocytes, microglia or neuronal cells, it is unlikely that the damage incurred by these cell populations is mediated by direct recognition of HTLV-1 antigens presented on these cells by HTLV-1-specific CD8+ T-cells. Two potential mechanisms have been proposed to explain how HTLV-1-specific CD8+ T-cells could mediate damage to non-HTLV-1-infected cells. Firstly, this damage could occur simply as a bystander effect of HTLV-1-specific CD8+ T-cells recognizing the HTLV-1-infected T-cells that have been shown to accumulate in the CNS. Secondly, it has been suggested that there may be cross-reactivity of HTLV-1-specific CD8+ T-cells with self antigens expressed in neurons. Here, we present and explore a third possibility, that HTLV-1-infection may elicit human endogenous retrovirus (HERV) specific CD8+ T-cell responses which, in turn, may be reactive against neurons and other tissues where HERV antigens may be expressed.
Human endogenous retroviruses, the DNA remnants of ancient infectious retroviruses, comprise approximately 8% of the human genome. This complement of HERVs is diverse, with all three major branches of the retroviral tree represented: gamma-epsilon, spuma and delta-lenti-alpha-beta retroviruses. While the vast majority of HERV insertions do not contain intact open reading frames (ORFs), a recent study has illustrated the potential for even short, disrupted ORFs to serve as a source of immunologically relevant antigens, by mapping the epitope specificity of a renal cell carcinoma reactive CD8+ T-cell to a HERV-E-derived peptide. As a gross generalization, HERV antigens are not thought to be expressed in healthy tissues, but rather are associated with transformed cell lines and disease states. We have previously presented the hypothesis that HIV-1 infection would result in the induction of HERV expression, either through rescue of HERV function by HIV-1, for example, nuclear export of unspliced HERV RNA by HIV-1 Rev, or through the degradation of host restriction factors which may normally keep HERV expression in check. Indirect support for this hypothesis is provided by our observation that T-cell responses to peptides derived from diverse HERV families can be detected in HIV-1-infected, but not uninfected, subjects[12–14]. We have since focused in on the HERV-K(HML-2) lineage of young and relatively intact endogenous retroviruses, and determined that HIV-1 infection results in the induction of Gag and Env protein expression. A HERV-K(HML-2)-Env-specific CD8+ T-cell clone, isolated from an HIV-1-infected subject, specifically responded to and eliminated cells infected with diverse isolates of HIV-1, HIV-2 and SIV. The establishment of this precedent for induction of HERV expression by an exogenous retroviruses led us to consider whether HTLV-1, as a second human retrovirus, may also drive the expression of HERVs. This idea draws support from a very recent study which demonstrated that the HTLV-1-Tax protein drives transcription from the long-terminal repeats (LTR) of HERV-W8, HERV-H and HERV-K. In the current study we sought to obtain indirect evidence for the in vivo expression of HERV antigens, by screening PBMC from HTLV-1-infected subjects for T-cell responses to HERV-derived peptides, with a particular focus on the HERV-K(HML-2) lineage.
In addition to providing evidence for the HTLV- 1-induced expression of HERVs, the detection of HERV-specific T-cell responses in HTLV-1-infected subjects would open a line of inquiry into determining whether such responses may contribute to the pathogenesis of HAM/TSP. It is intriguing to note, given the similarities in clinical presentation and pathogenesis of HAM/TSP and multiple sclerosis (MS), that numerous studies have associated HERV polymorphisms or expression levels with MS (reviewed in). Specifically, HERV-H, HERV-W, and HERV-K have been associated with MS, the same families shown to be induced at the transcriptional level by HTLV- 1-Tax. We posited that HTLV-1 infection would result in the expression of HERV antigens and the induction of HERV-specific T-cell responses, which would, in turn, be reactive against neurons and other tissues which may express low levels of HERV antigens, resulting in the pathologies of HAM/TSP. Consistent with this, we have recently performed a comprehensive screening of human tissues for HERV-K(HML-2)-Gag expression and detected this protein in the central nervous system (neurons, purkinje cells and occasionally ependymal cells), as well as in the endocrine pancreas (few cells in the β-islets of Langerhans) and in the exocrine pancreas (ductal epithelial cells) (Sacha et al. Submitted).
Results and discussion
To test this hypothesis, we obtained peripheral blood mononuclear cells (PBMC) from 15 HTLV-1-infected subjects from Sao Paulo, Brazil, and from 8 age/sex matched uninfected subjects from the same geographical area. The HTLV-1-infected subjects comprised 9 individuals who were asymptomatic, 5 with HAM/TSP and 1 subject who we define as pre-ATL. In the latter subject, significant immunophenotype alterations were present at the time of sample collection, but the criteria for an ATL diagnosis were not present. Six months later, this patient was diagnosed with the lymphomatous variant of ATL (Table1). These PBMC were screened for T-cell responses by IFN-γ ELISPOT. The primary objective of these experiments was to comprehensively screen for responses to HERV-K(HML-2) Gag and Env using overlapping peptide pools spanning these antigens. We did also screen these subjects for responses to peptides derived from diverse HERV families selected based on predicted binding to common MHC-I alleles, many of which have previously been demonstrated to be immunoreactive HIV-1-infected subjects[12–14]. However, as we had no prior knowledge of the HLA-types of subjects in our cohort, the latter screening was performed without consideration of whether subjects possessed the necessary restricting alleles, imposing limitations on our ability to draw conclusions from these data. Peptides tested were thus: i) 15mer peptides, overlapping by 11 amino acids, spanning HERV-K(HML-2) Gag and Env – tested in pools of 10 peptides each at 25 μg/ml/peptide ii) 29 shorter peptides (9–13 amino acids) derived from HERV-K(HML-2), HERV-L, HERV-H and HERV-W (Table2) – tested in pools of 7 peptides each at 10 μg/ml/peptide. These peptides have been previously described[12–14], and were selected based on predicted binding to the MHC-I alleles A02, B58, B07 or B35 iii) the immunodominant HLA-A02-restricted Tax 11–19 peptide tested at 10 μg/ml iv) a pool of CMV, EBV and influenza derived peptides (CEF pool) as a positive control tested at 2.5 μg/ml/peptide v) Staphylococcus enterotoxin B (SEB) tested at 2.5 μg/ml. ELISPOTS were performed as previously described, except that stimulation periods were limited to 6 hours. Longer periods of in vitro culture of PBMC derived from HTLV-1-infected subjects resulted in very high levels of spontaneous IFN-γ production, as has been previously reported.
We observed a total of 12 HERV-specific responses in HTLV-1-infected individuals which passed our pre-established criteria of > 50 sfu/106 PBMC after background subtraction and > 3x background versus 3 HERV-specific responses in uninfected subjects. These responses were, however, of borderline magnitude and replicates between duplicate wells were poor. We sought to repeat these responses with new aliquots of cells from the same patient time-points using the same peptide pools and, in parallel, to deconvolute responses to pools by testing individual peptides. In these follow-up experiments we observed that none of these responses were reproducible. The incidence of these false positive responses was related to ‘noisier’ ELISPOTs than those typically observed using PBMC from HTLV-1-uninfected subjects. Although we sought to minimize the impact of spontaneous IFN-γ production on the assay by reducing stimulation periods to 6 hours, we observed higher levels of background in ELISPOTs performed on HTLV- 1-infected subjects than in uninfected controls. Despite this background, we observed clear and reproducible responses to the HLA-A02-restricted Tax peptide in 6 HTLV-1 subjects, and an absence of this response in uninfected subjects. Responses to CEF pool were detected in the majority of both HTLV-1-infected and uninfected subjects.
We next summed the background-subtracted responses to: HERV-K(HML-2)-Env, HERV-K(HML-2)-Gag and predicted optimal HERV peptides (HERV-Opt). In screens performed on non-HTLV-1-infected subjects ELISPOT results with negative values after background subtraction are atypical and generally reported as zero. In the present case, however, the variability caused by spontaneous cytokine production, and summed over multiple peptide pools in the cases of HERV-K(HML-2)-Gag and Env (17 and 18 pools respectively), manifested as cumulative totals which spanned a range from substantially more than to less than zero. As zeroing these negative values would give the false impression of positive responses in HTLV-1-infected subjects (instead of just cumulative noise) these data are presented as is. We observed a lack of statistically significant differences between the magnitudes of responses detected in HTLV-1-infected versus uninfected subjects for each of the peptide sets tested (Figure1). Notably, however, we did observe a trend towards the detection of responses to the HLA-A02-restricted Tax peptide in HTLV-1-infected subjects (p=0.135). The failure of this response to reach significance is likely tied to the fact that subjects were tested without a prior knowledge of their HLA types. Thus, only a minority (approximately 40%) of subjects would be expected to have the capacity to present this particular peptide. A lack of statistically significant differences in magnitudes of responses to HERV, Tax or CEF peptides were also observed when we broke the HTLV-1-infected subjects into asymptomatic versus HAM/TSP groups (data not shown). We did, however, observe a positive correlation between the magnitude of responses to the HLA-A02-restricted Tax peptide and HTLV-1 proviral load (considering all subjects, not just responders) (R=0.785, p=0.002, Table3). This is in agreement with observations from two previous reports[19, 20]. We did not observe any such correlation with magnitudes of responses to HERV-K(HML-2) Env or Gag, to the optimized HERV peptide pool or to CEF (Table3). Thus the data are indicative of a lack of detectable T-cell responses to the HERV-derived peptides considered in the present study.
The current study does not provide any data to support the hypothesis that the induction of HERV expression and HERV-specific T-cell responses in HTLV-1-infected subjects may play a role in the pathogenesis of HAM/TSP. This may indicate that, in contrast to HIV-1, HTLV-1 does not induce HERV expression. We have recently demonstrated that the HIV-1-Vif protein is necessary for the induction of HERV-K(HML-2) expression, thus the lack of a Vif equivalent in HTLV-1 may at least partially underlie this difference. There are a number of important caveats, however, which prevent the current study from comprehensively ruling out a role for HERVs in HTLV-1 pathogenesis. Firstly, the HERV-derived peptides tested in the current study only represent a small fraction of the total HERV antigens encoded in the human genome. In particular, HERV-H and HERV-W, which were shown to be transcriptionally induced to the greatest level by Tax in a recent study, were only represented by one peptide each (this study was published after our experiments were completed). Secondly, the HERV peptides, predicted to be optimal binders for either HLA-A02, B07, B58 or B35, were tested without a prior knowledge of the HLA-types of the test subjects in this cohort. HLA-A02 is common enough (~40% prevalence in this cohort) that responses to immunodominant HLA-A02-restricted HERV-derived peptides would likely have been detected in this cohort, as is evidenced by our detection of responses to the HLA-A02-restricted Tax peptide. However, the detection of responses restricted by the less common B07, B58, or B35 alleles was unlikely. Indeed, these were tested as an adjunct to the primary objective of the current study, which was to comprehensively screen for responses to HERV-K(HML-2) Gag and Env. Importantly, in our previous studies measuring responses to these predicted optimal epitopes in HIV-1-infected subjects, these peptides were tested against a cohort of subjects selected on the basis of having one or more of these restricting alleles. In the current manuscript, we do comprehensively screen for T-cell responses to HERV-K(HML-2) Gag and Env using overlapping 15mer peptides. This approach is not reliant upon prior knowledge of a subject’s HLA-type, and yet a lack of T-cell responses was still observed. When we have previously screened HIV-1-infected subjects with these same peptides we observed only infrequent responses, with 3 responses detected in the 27 subjects who were screened. Thirdly, the peptides used in the current study represent only a small fraction of the potential HERV antigens encoded by 8% of our genome. The current study is best interpreted as a starting point, effectively ruling out the presence of responses to HERV-K(HML-2) Gag and Env that are detectable using 15mer peptide screens. The recent demonstration that HERV-H and HERV-W are most potently induced at the transcriptional level by HTLV-1-Tax, with a lesser induction of HERV-K, suggests that future studies should focus on screening for T-cell responses to these HERV families. Finally, while we have focused exclusively on T-cell responses as potential mediators, there are also other potential mechanisms by which putative HTLV-1-mediated HERV induction could result in pathogenesis – for example, through antibody mediated immunopathology or through direct effects of HERV gene products on the affected tissue (toxicities of Env proteins for example). Although limited by the caveats given above, the current study argues against a role for an involvement of HERVs in the immunopathogenesis of HTLV-I infection.
Samples were collected under an IRB approved protocol. PBMC were isolated through Ficoll gradient enrichment and cryopreserved. ELISPOT assays were performed as previously described, but with a 6 hour, rather than overnight, stimulation period. Peptides were tested as follows: 15mer peptides, overlapping by 11 amino acids, spanning HERV-K(HML-2) Gag and Env – tested in pools of 10 peptides each at 25 μg/ml/peptide ii) 29 shorter peptides (9–13 amino acids) derived from HERV-K(HML-2), HERV-L, HERV-H and HERV-W (Table2) – tested in pools of 7 peptides each at 10 μg/ml/peptide iii) the immunodominant HLA-A02-restricted Tax 11–19 peptide tested at 10 μg/ml iv) a pool of CMV, EBV and influenza derived peptides (CEF pool) tested at 2.5 μg/ml/peptide. Staphylococcus enterotoxin B (SEB) was used at 2.5 μg/ml. A response was considered to be positive if it passed both the criteria of: i) > 50 sfu/million PBMC after background subtraction and ii) > 3x background. Statistics were performed using Prism Graphpad 2.0 software.
Human T-cell leukemia virus type 1
Human endogenous retrovirus
HTLV-1-associated myelopathy/tropical spastic paraparesis
Human immunodeficiency virus type 1
Peripheral blood mononuclear cells
Cytomegalovirus, epstein barr virus, influenza (peptide pool)
Adult T-cell leukemia
Human leukocyte antigen.
de The G, Bomford R: An HTLV-I vaccine: why, how, for whom. AIDS Res Hum Retroviruses. 1993, 9: 381-386. 10.1089/aid.1993.9.381.
Murphy EL, Hanchard B, Figueroa JP, Gibbs WN, Lofters WS, Campbell M, Goedert JJ, Blattner WA: Modelling the risk of adult T-cell leukemia/lymphoma in persons infected with human T-lymphotropic virus type I. Int J Cancer. 1989, 43: 250-253. 10.1002/ijc.2910430214.
Kaplan JE, Osame M, Kubota H, Igata A, Nishitani H, Maeda Y, Khabbaz RF, Janssen RS: The risk of development of HTLV-I-associated myelopathy/tropical spastic paraparesis among persons infected with HTLV-I. J Acquir Immune Defic Syndr. 1990, 3: 1096-1101.
Osame M, Usuku K, Izumo S, Ijichi N, Amitani H, Igata A, Matsumoto M, Tara M: HTLV-I associated myelopathy, a new clinical entity. Lancet. 1986, 1: 1031-1032.
Goon PK, Igakura T, Hanon E, Mosley AJ, Asquith B, Gould KG, Taylor GP, Weber JN, Bangham CR: High circulating frequencies of tumor necrosis factor alpha- and interleukin-2-secreting human T-lymphotropic virus type 1 (HTLV-1)-specific CD4+ T cells in patients with HTLV-1-associated neurological disease. J Virol. 2003, 77: 9716-9722. 10.1128/JVI.77.17.9716-9722.2003.
Montanheiro PA, Oliveira AC, Posada-Vergara MP, Milagres AC, Tauil C, Marchiori PE, Duarte AJ, Casseb J: Human T-cell lymphotropic virus type I (HTLV-I) proviral DNA viral load among asymptomatic patients and patients with HTLV-I-associated myelopathy/tropical spastic paraparesis. Brazilian journal of medical and biological research = Revista brasileira de pesquisas medicas e biologicas / Sociedade Brasileira de Biofisica. 2005, 38: 1643-1647.
Olindo S, Lezin A, Cabre P, Merle H, Saint-Vil M, Edimonana Kaptue M, Signate A, Cesaire R, Smadja D: HTLV-1 proviral load in peripheral blood mononuclear cells quantified in 100 HAM/TSP patients: a marker of disease progression. J Neurol Sci. 2005, 237: 53-59. 10.1016/j.jns.2005.05.010.
Sakai JA, Nagai M, Brennan MB, Mora CA, Jacobson S: In vitro spontaneous lymphoproliferation in patients with human T-cell lymphotropic virus type I-associated neurologic disease: predominant expansion of CD8+ T cells. Blood. 2001, 98: 1506-1511. 10.1182/blood.V98.5.1506.
Levin MC, Lee SM, Morcos Y, Brady J, Stuart J: Cross-reactivity between immunodominant human T lymphotropic virus type I tax and neurons: implications for molecular mimicry. J Infect Dis. 2002, 186: 1514-1517. 10.1086/344734.
Tristem M: Identification and characterization of novel human endogenous retrovirus families by phylogenetic screening of the human genome mapping project database. J Virol. 2000, 74: 3715-3730. 10.1128/JVI.74.8.3715-3730.2000.
Takahashi Y, Harashima N, Kajigaya S, Yokoyama H, Cherkasova E, McCoy JP, Hanada K, Mena O, Kurlander R, Tawab A: Regression of human kidney cancer following allogeneic stem cell transplantation is associated with recognition of an HERV-E antigen by T cells. J Clin Invest. 2008, 118: 1099-1109.
Garrison KE, Jones RB, Meiklejohn DA, Anwar N, Ndhlovu LC, Chapman JM, Erickson AL, Agrawal A, Spotts G, Hecht FM: T cell responses to human endogenous retroviruses in HIV-1 infection. PLoS Pathog. 2007, 3: e165-10.1371/journal.ppat.0030165.
SenGupta D, Tandon R, Vieira RG, Ndhlovu LC, Lown-Hecht R, Ormsby CE, Loh L, Jones RB, Garrison KE, Martin JN: Strong human endogenous retrovirus-specific T cell responses are associated with control of HIV-1 in chronic infection. J Virol. 2011, 85: 6977-6985. 10.1128/JVI.00179-11.
Tandon R, SenGupta D, Ndhlovu LC, Vieira RG, Jones RB, York VA, Vieira VA, Sharp ER, Wiznia AA, Ostrowski MA: Identification of human endogenous retrovirus-specific T cell responses in vertically HIV-1-infected subjects. J Virol. 2011, 85: 11526-11531. 10.1128/JVI.05418-11.
Jones RB, Garrison KE, Mujib S, Mihajlovic V, Aidarus N, Hunter DV, Martin E, John VM, Zhan W, Faruk NF, Gyenes G, Sheppard NC, Priumboom-Brees IM, Goodwin DA, Chen L, Rieger M, Muscat-King S, Loudon PT, Stanley C, Holditch SJ, Wong JC, Clayton K, Duan E, Song H, Xu Y, Sengupta D, Tandon R, Sacha JB, Brockman MA, Benko E, Kovacs C, Nixon DF, Ostrowski MA: HERV-K-specific T cells eliminate diverse HIV-1/2 and SIV primary isolates. J Clin Invest. 2012, 122 (12): 4473-4489. 10.1172/JCI64560. Epub 2012 Nov 12
Toufaily C, Landry S, Leib-Mosch C, Rassart E, Barbeau B: Activation of LTRs from Different Human Endogenous Retrovirus (HERV) Families by the HTLV-1 Tax Protein and T-Cell Activators. Viruses. 2011, 3: 2146-2159. 10.3390/v3112146.
Antony JM, Deslauriers AM, Bhat RK, Ellestad KK, Power C: Human endogenous retroviruses and multiple sclerosis: innocent bystanders or disease determinants. Biochim Biophys Acta. 2011, 1812: 162-176. 10.1016/j.bbadis.2010.07.016.
Hanon E, Goon P, Taylor GP, Hasegawa H, Tanaka Y, Weber JN, Bangham CR: High production of interferon gamma but not interleukin-2 by human T-lymphotropic virus type I-infected peripheral blood mononuclear cells. Blood. 2001, 98: 721-726. 10.1182/blood.V98.3.721.
Kubota R, Nagai M, Kawanishi T, Osame M, Jacobson S: Increased HTLV type 1 tax specific CD8+ cells in HTLV type 1-asociated myelopathy/tropical spastic paraparesis: correlation with HTLV type 1 proviral load. AIDS Res Hum Retroviruses. 2000, 16: 1705-1709. 10.1089/08892220050193182.
Wodarz D, Hall SE, Usuku K, Osame M, Ogg GS, McMichael AJ, Nowak MA, Bangham CR: Cytotoxic T-cell abundance and virus load in human immunodeficiency virus type 1 and human T-cell leukaemia virus type 1. Proceedings. 2001, 268: 1215-1221.
Jones RB, John VM, Hunter DV, Martin E, Mujib S, Mihajlovic V, Burgers PC, Luider TM, Gyenes G, Sheppard NC: HERV-K(HML-2) Gag and Env specific T cell responses are infrequently detected in HIV-1-infected subjects using standard peptide-matrix based screening. Clin Vaccine Immunol. 2012, 19 (2): 288-292. 10.1128/CVI.05583-11. Epub 2011 Dec 28
This work was supported in part by funds from the National Institutes of Health (AI076059 and AI084113). Additional support was provided by the Brazilian Program for STD and AIDS, Ministry of Health (914/BRA/3014 - UNESCO/Kallas), the São Paulo City Health Department (2004–0.168.922-7/Kallas), Fundação de Amparo a Pesquisa do Estado de São Paulo (04/15856-9/Kallas and 2010/05845-0/Kallas and Nixon), the John E. Fogarty International Center (D43 TW00003), and the Fundação de Amparo a Pesquisa do Estado de São Paulo (04/15856-9/Kallas and 2010/05845-0/Kallas and Nixon). RBJ gratefully acknowledges salary support from the Ontario HIV Treatment Network (OHTN).
The authors declare that they have no competing financial interests.
RBJ, FEL, DFN, MAO, EGK conceived of the project and designed experiments. FEL, ACS recruited study participants and AMH processed samples. RBJ, AMH, MAO performed the ELISPOT assays. RBJ and AMH analyzed the ELISPOT data. DFN, MAO, EGK supervised the project. Experiments were performed in the laboratory of EGK. RBJ wrote the manuscript with input from all co-authors. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.