- 1Department of Immunology, Ophthalmology and ENT, School of Medicine, Complutense University, 12 de Octubre Health Research Institute (imas12), Madrid, Spain
- 2Hospital 12 de Octubre Health Research Institute (imas12), Madrid, Spain
- 3Pediatric Hematology and Oncology Unit, University Hospital Virgen de la Salud, Toledo, Spain
- 4Pediatric Hematology and Oncology Unit, University Hospital Miguel Servet, Zaragoza, Spain
- 5Immunodeficiencies Unit, Pediatrics, University Hospital 12 octubre, Madrid, Spain
- 6Complutense University School of Medicine, Madrid, Spain
- 7Department of Immunology, University Hospital 12 Octubre, Madrid, Spain
Cernunnos/XLF deficiency is a rare primary immunodeficiency classified within the DNA repair defects. Patients present with severe growth retardation, microcephaly, lymphopenia and increased cellular sensitivity to ionizing radiation. Here, we describe two unrelated cases with the same non-sense mutation in the NHEJ1 gene showing significant differences in clinical presentation and immunological profile but a similar DNA repair defect.
Introduction
Clinical Presentation and Laboratory Test Results
We report two patients harboring the same homozygous mutation in NHEJ1 gene. Strikingly, their clinical phenotypes differed markedly. One presented with severe combined immunodeficiency whereas the other only had isolated thrombocytopenia and macrocytosis. Clinical characteristics are summarized in Table 1. Both patients presented severe lymphopenia: immunophenotyping of Patient 1 (P1) showed severe T and B-cell lymphopenia and normal to elevated NK cells (T−B−NK+ phenotype) in contrast to Patient 2 (P2) who presented a milder immunophenotype. T-cell compartment also showed differences between the two patients: P1 had a more senescent T-cell phenotype with increase in CD4+ and CD8+ effector memory (CCR7−CD45RA−) and decrease naïve and recent thymic emigrants (RTE) (CD4+CD45RA+CD31+) T cells and a higher proportion of activated T-cells (CD3+HLA-DR+) (Table 2). Of note, P1 had a severe decrease in CD8+ population mimicking the one reported in ZAP70 or HLA class I deficiencies among patients with combined immunodeficiency/severe combined immunodeficiency.
B cells were severely reduced in both patients. The B-cell profile in P2 showed a normal proportion of naive, unswitched memory, switched memory, transitional, and plasmoblast B-cells. Serum immunoglobulins were reduced in P1. Due to age limitation, specific antibody response was only tested in P2, showing a specific antibody deficiency: IgG and IgG2 levels against pneumococcus antigen were reduced while IgG levels against tetanus toxoid antigen showed normal values after vaccination. Regarding to the immunophenotype in T and B-cells, P1 had very low TRECs and KRECs copies, whereas P2 had preserved levels (Table 2).
The clinical presentation and immunologic features of both patients lead to suspicion of a primary immunodeficiency (PID). Accordingly, an in house targeted NGS sequencing panel for 192 PID related genes (Table 3) was performed and revealed a homozygous nucleotide substitution in exon 2 of NHEJ1 gene (NM_024782, c.169C>T) that affects codon 57 and is predicted to result in a severely truncated protein (p.R57X), this mutational change has been previously described in two Cernunnos /XLF defective patients similar to P2 patient displaying microcephaly and slight lymphopenia (1–6). Sanger sequencing confirmed the variants in the patients which were inherited in autosomal recessive fashion from their healthy parents (Figure 1).
Figure 1. Pedigree and genomic sequence analysis of NHEJ1 gene. Genomic NHEJ1 DNA sequence showing c.169C>T mutation and Family tree in P1 (A) and P2 (B).
Increased sensitivity and decreased double strand breaks (DSB) rejoining to ionizing irradiation (IR) is a feature of Cernunnos /XLF deficient cells. Accordingly, a cell-survival assay revealed an increased IR sensitivity in fibroblasts from P2 similar to that observed in LIG4-deficient fibroblasts (Figure 2A). We also analyze DSB rejoining in primary fibroblasts from both patients (P1 and P2) by enumerating the rate of loss of γ-H2AX foci following exposure to γ-IR. The results showed impaired DSB rejoining after treatment with γ-irradiation and etoposide in both patients (Figure 2B). These findings are in agreement with previously published results in Cernunnos and LIG4-defective patients (7).
Figure 2. Cellular Response to DNA Damage. (A) Clonogenic survival assay in XLF-deficient fibroblasts. Cell survival following IR with γ-rays (1 and 3Gy) was assessed in primary fibroblasts from normal control (C) and patient (P2). LIG4-deficient fibroblasts (LIG4−/−) were used as a radiosensitive control. The results represent the mean and standard deviation of two separate experiments and are expressed as percentages of survival cells relative to unirradiated primary fibroblasts. (B) Top panel: Primary fibroblasts from control (C) and patients (P1 and P2) were irradiated with 3 Gy and fixed at given time points post-irradiation before staining with anti- γ -H2AX. Nuclei were stained with DAPI. Numbers of γ -H2AX foci per nucleus were determined at indicated time points after irradiation (average number of γ -H2AX foci per nucleus in 30 cells). Error bars represent the SD from 3 independent experiments. Bottom panel: γ-H2AX detection by flow cytometry was performed in PBMCs from P2 and control. Mean fluorescence intensities (MFI) are shown as histograms (unirradiated, 2 and 24 h) compared to isotype (IC). Persistence of γ H2AX signal at 24 h post-treatment in P2 is indicative of a general DNA repair defect. *p < 0.05, **p < 0.01.
Comparative testing of DNA repair was also performed on T cells from the second patient (P2). Upon exposure of PBMCs to 10 Gy, induction and resolution of DNA damage was measured at various time points and a delayed kinetics of DNA repair were observed in P2 compared to a healthy control (Figure 2B).
Methods
Cell Culture
Primary fibroblasts were grown in minimal essential medium (MEM) supplemented with 10% fetal calf serum (FCS), penicillin and streptomycin.
Immunofluorescence and Antibodies
To characterize the repair capacity of the patient's cells we scored the in situ modification of the histone variant H2AX, which is phosphorylated proximal to sites of DNA double-strand breaks. The number of phosphorylated H2AX (γH2AX) foci in a nucleus is reported to be directly proportional to the number of DSBs, and de-phosphorylation coincides with DSB repair. Primary skin fibroblasts were irradiated with ionizing irradiation (137Cs) or treated with 20 mM Etoposide for 1 h. After the indicated treatments, the slides were washed with PBS, fixed using 4% formaldehyde for 10 min at room temperature and permeabilized with 0.5% Triton X-100 for 5 min. Cells were incubated with primary antibodies for 1 h, the slides were washed with PBS and the bound antibodies were revealed by IgG Alexa fluor antibodies (Invitrogen). Nuclei were counterstained with DAPI, and slides were mounted for immunofluorescence. Images were taken with a fluorescent microscopy (Zeiss AxioImages.A1, Carl Zeiss). In a single experiment at least 30 cells per sample were counted.
Survival Assay
Primary skin fibroblasts were irradiated with ionizing radiation (137Cs). After irradiation, the cells were seeded at a density of 1 × 104 cells/mL in T75 flasks in triplicate. To evaluate cell sensitivity to γ-IR (1 and 3 Gy), adherent cells were trypsinized and counted 11 days later.
Flow Cytometry
Proportions and lymphocyte count of T-, B-, and NK-cells were determined in blood samples using conjugated mouse anti-human monoclonal antibodies and data were collected by flow cytometry using a Navios Cytometer (Beckman Coulter, Madrid, Spain) and analyzed with Kaluza 1.5a software (Beckman Coulter, Indianapolis IN, US).
PBMCs from patient and healthy controls were irradiated with 10Gy, fixed and stained for CD3, CD19, and phospho-histone H2AX. Mean fluorescence intensities (MFI) of γH2AX were evaluated on gated CD3+ lymphocytes.
Immunoglobulins
Total serum immunoglobulins (IgG, IgA, IgM, and IgE) were measured by nephelometry (Beckman Coulter, Madrid, Spain).
NGS and Sanger Sequencing
Genomic DNA was extracted from EDTA whole blood using a MagNa Pure Compact Nucleic Acid Isolation Kit (Roche, Madrid, Spain). Missense NHEJ1 mutation was identified by targeted next-generation sequencing with an in-house designed panel of 192 genes involved in primary immunodeficiencies (PID) (Ampliseq, Thermo Fisher, Madrid, Spain) and confirmed by PCR and Sanger sequencing using an ABI PRISM 3130 genetic analyzer.
TRECs and KRECs
TRECs, KRECs, and beta-actin (ACTB) copy numbers were determined from dried blood spots (DBS, punches of 3.2 mm) using triplex real-time quantitative polymerase chain reaction (RT-PCR) (TIB MOLBIOL) and run in a Light Cycler 480 II from Roche Diagnostics.
Ethics Statement
The protocols of this study were approved by the Institutional Review Board of Hospital Universitario 12 de Octubre (Madrid, Spain) and written informed consent was obtained from all subjects/caregivers in accordance with the Declaration of Helsinki.
Background
DNA non-homologous end-joining (NHEJ) is the major DNA double strand break (DSB) repair pathway in mammalian cells. NHEJ also functions during immune development rejoining the programmed DSBs introduced during V(D)J recombination (7). Most of the patients deficient in NHEJ components display radiosensitivity and severe combined immunodeficiency (SCID), a phenotype which has been called radiosensitive SCID (RS-SCID). To date, mutations in five genes encoding components of the NHEJ pathway, LIG4 (encoding DNA ligase IV), NHEJ1 (encoding Cernunnos), PRKDC (encoding DNA-PKcs), DCLRE1C (encoding Artemis), and XRCC4 (encoding XRCC4) have been identified in patients. Although RS-SCID represents the extreme phenotype, radiosensitivity coupled with variable immunodeficiency ranging from Omenn's Syndrome to CID has been observed. Additional features are also observed in some patients, including microcephaly and severe growth delay. However, all of these patients are at risk of malignancies, particularly lymphoma (6).
History
Cernunnos, was identified through both cDNA complementation of cells derived from a IR-sensitive immunodeficient patient (1) and through a yeast two-hybrid screen for XRCC4-interacting partners (8). Cernunnos is a homolog of Nej1, one of the NHEJ factors identified in yeast (9, 10). Data from several studies suggest that although Cernunnos may not be strictly required for NHEJ, the loss of Cernunnos does affect NHEJ and interestingly DNA repair defect in these patients is quite similar to that found in LIG4 syndrome patients (11).
Review of Similar Cases
Cernunnos/XLF deficiency in human results in extreme sensitivity to IR, microcephaly, and growth retardation, but the effect on the immune system is variable (1–4, 8). Mutations in NHEJ1 have been described previously in 27 patients with clinical features comparable to LIG4 deficiencies. Of contrast, XRCC4 deficiency shares the radiosensitivity and neurological impairment but not overt immunodeficiency (5, 12, 13).
Most of the Cernunnos/XLF deficient patients are hypersensitive to IR and have a significant NHEJ defect (1, 14) similar to LIG4 deficiency, which belong to NHEJ ligation complex; immunodeficiency, however, is milder compared with defects in NHEJ factors involved in hairpin opening as Artemis and DNA-PKcs.
In this study, we reported two patients who present the same c.169C>T mutation in NHEJ1 gene but different immunologic features (Table 2). P2 presented with mild T lymphopenia, hypersensitivity and NHEJ repair defect, typical for patients with Cernunnos/XLF defects (1, 5). On the other hand, P1 presented a more severe phenotype (T-B-); hypersensitivity and NHEJ repair defect, however, was similar to P2. These findings indicate that patients with same NHEJ1 mutation and DNA repair defect may show great variability in the clinical phenotype.
Patients with homozygous mutations (p.R178X) in NHEJ1 gene have been previously reported. Two patients died at 1.5 and 4 years (1, 2), while another of the patients remains alive at the age of 8 years (without HSCT) (5). However, none of these patients presented the severe T lymphopenia observed in our first patient.
One explanation for the different degree of lymphopenia found in Cernunnos-deficient patients (P1 and P2) might be the existence of alternative DNA repair proteins that only work to repair DSBs generated during lymphocyte development (3). It has been described that RAG complex and Cernunnos functionally overlap in the repair of DNA breaks during antigen receptor assembly ensuring stabilization of DNA ends after DNA cleavage by RAG (15). ATM and/or ATM-dependent DRR factors, as 53BP1, would contribute to the RAG-DSB stabilization and the recruitment of NHEJ pathway proteins. Severely impaired joining of RAG-generated DSBs in cells that are deficient for Cernunnos and either ATM, 53PBP1or H2AX has been observed (16). Thus, genetic factors such as mutations or polymorphisms in some of these proteins could affect the VDJ recombination process and explain the severe T lymphopenia observed in P1. As Cernunnos is not usually found in the context of a patient with a SCID phenotype we cannot rule out a digenic cause (particularly 53BP1, as other VDJ recombination genes were wild-type in our panel).
Discussion
Diagnosis and Treatment
The assignment of a timely and accurate diagnosis is of paramount importance in the management of patients with defects in DNA repair, as HSCT is the only curative therapy available. Usually the repair defect in these disorders is assessed by immunofluorescence assays of irradiation-induced γ-H2AX foci using skin fibroblasts. Flow cytometry (FC) can be applied as a rapid diagnostic tool for DNA repair disorders (17, 18). Therefore, we have used flow cytometry to analyze PBMCs from P2 and the results showed a DNA repair defect similar to that obtained in skin fibroblasts (data not shown) (19).
Hence, a high throughput, sensitive and reliable assay to quantify γ-H2AX foci in PBMCs isolated from blood samples would be a valuable tool to diagnose these patients and thus allow HSCT without delay.
In addition, it would also be helpful in cancer patients to individualize and to guide the dosing of ionizing radiation (IR) and/ or genotoxic agents to avoid accumulation of cells with genomic instability that could accelerate cancer development.
In the era of newborn screening an abnormal TREC assay should be followed by NGS approach as Cernunnos/XLF deficiency may present early in life as SCID, as other RS-SCID defects (20). Since genetic diagnosis takes time, functional radiosensitivity assays in peripheral blood may lead to the correct diagnosis and avoid exposure to alkylating agents during the conditioning regimen even in the absence of a genetic diagnosis (21).
The patients presented in this work are alive and well and both patients after undergoing HSCT. Of note, P2 has survived to age 7 years. It has been reported that the Cernunnos/XLF deficient patients may survive the first years of life, or even up to 18 years, without HSCT. However, due to the comorbilities that these patients face, it is highly recommended HSCT pre-emptively, rather than expect the appearance of a malignant refractory disease disregarding the age at diagnosis (21).
Given the extreme rarity of the disorder, the appearance of two unrelated homozygous cases in non-consanguineous family with the same mutation seems unlikely. However, we can offer three reasons to discard consanguinity in both kindreds: First. Geographical: Both families are more than 450 km apart each other. Both kindreds deny any consanguinity (even far) ties. Second, Molecular: HLA typing were done in both families. It is shown HLA haplotypes generated by segregation analysis in both families: Family 1, Father: A*30, B*18, DRB1*03 and A*02, B*40, DRB1*08; Mother: A*01, B*37, DRB1*13 and A*-, B*57, DRB1*07. Family 2, Father: A*02, B*39, DRB1*11 and A*29, B*51, DRB1*04; Mother: A*11, B*40, DRB1*04 and A*11, B*07, DRB1*15. In conclusion, both families did not share any HLA haplotype (Table 4). Third, No other mutations in PID genes that associate radiosensibility (LIG4, PRKDC, DCLRE1C, ATM, RNF168) were found.
Concluding Remarks
We report two unrelated cases diagnosed with Cernunnos/XLF deficiency. Both patients showed a similar DNA repair defect and increased cellular sensitivity to ionizing radiation as per in vitro assays. However, clinical presentation and immunological profile were extremely different. The first patient showed a senescent phenotype with decreased TRECs, RTE and naive T-cells counts suggesting that the sustained self-renewal of T-cell pool was impaired, while P2 presented only a slight reduction in T-cell counts as it has been reported in Cernunnos/XLF deficient patients reported up to date with normal TRECs/KRECs and IgG levels (22).
Author Contributions
JB-E and RR-G performed the laboratory work for this study and computational predictions. MR, RR-G, ML-N, RC-G, MC-P, MP, and JB-E performed the laboratory work for this study. CR, CM, ND-P, NS-R, CR, and LG-G were responsible for the clinical management of the patients. LA, MR, JB-E, and LG-G designed the research and drafted the manuscript. All authors approved the final version of this manuscript.
Funding
This research was supported by grants from the Spanish Health Research Fund FIS PI16/2053 to LA and LG-G, PI16/0044 to LG-G, MINECO grants SAF2014-54708-R, SAF2016-81876-REDT, and RTI2018-095673-B-I00. CAM grant B2017/BMD3673. The project has been co-financed with FEDER funds.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We thank both families and Pablo Morales, Isabel Castro, María José Díaz-Madroñero, Marta Martín Cabezudo, and Maria Teresa Fernandez González for their personal and technical assistance and Jaime de Inocencio for critical review of the article.
References
1. Buck D, Malivert L, de Chasseval R, Barraud A, Fondanèche MC, Sanal O, et al. Cernunnos, a novel nonhomologous end-joining factor, is mutated in human immunodeficiency with microcephaly. Cell (2006) 124:287–99. doi: 10.1016/j.cell.2005.12.030
2. Dutrannoy V, Demuth I, Baumann U, Schindler D, Konrat K, Neitzel H, et al. Clinical variability and novel mutations in the NHEJ1 gene in patients with a Nijmegen breakage syndrome-like phenotype. Hum Mutat. (2010) 31:1059–68. doi: 10.1002/humu.21315
3. Sheikh F, Hawwari A, Alhissi S, Al Gazlan S, Al Dhekri H, Rehan Khaliq AM, et al. Loss of NHEJ1 protein due to a novel splice site mutation in a family presenting with combined immunodeficiency, microcephaly, and growth retardation and literature review. J Clin Immunol. (2017) 37:575–81. doi: 10.1007/s10875-017-0423-5
4. Cipe FE, Aydogmus C, Babayigit Hocaoglu A, Kilic M, Kaya GD, Yilmaz Gulec E. Cernunnos/XLF deficiency: a syndromic primary immunodeficiency. Case Rep Pediatr. (2014) 2014:614238. doi: 10.1155/2014/614238
5. IJspeert H, Rozmus J, Schwarz K, Warren RL, van Zessen D, Holt RA, et al. XLF deficiency results in reduced N-nucleotide addition during V(D)J recombination. Blood (2016) 128:650–9. doi: 10.1182/blood-2016-02-701029
6. Du L, Peng R, Björkman A, Filipe de Miranda N, Rosner C, Kotnis A, et al. Cernunnos influences human immunoglobulin class switch recombination and may be associated with B cell lymphomagenesis. J Exp Med. (2012) 209:291–305. doi: 10.1084/jem.20110325
7. Riballo E, Woodbine L, Stiff T, Walker SA, Goodarzi AA, Jeggo PA. XLF-Cernunnos promotes DNA ligase IV-XRCC4 re-adenylation following ligation. Nucleic Acids Res. (2009) 37:482–92. doi: 10.1093/nar/gkn957
8. Ahnesorg P, Smith P, Jackson SP. XLF interacts with the XRCC4-DNA ligase IV complex to promote DNA nonhomologous end-joining. Cell (2006) 124:301–13. doi: 10.1016/j.cell.2005.12.031
9. Kegel A, Sjöstrand JO, Aström SU. Nej1p, a cell type-specific regulator of nonhomologous end joining in yeast. Curr Biol. (2001) 11:1611–7. doi: 10.1016/S0960-9822(01)00488-2
10. Callebaut I, Malivert L, Fischer A, Mornon JP, Revy P, de Villartay JP. Cernunnos interacts with the XRCC4 x DNA-ligase IV complex and is homologous to the yeast nonhomologous end-joining factor Nej1. J Biol Chem. (2006) 281:13857–60. doi: 10.1074/jbc.C500473200
11. Girard PM, Kysela B, Härer CJ, Doherty AJ, Jeggo PA. Analysis of DNA ligase IV mutations found in LIG4 syndrome patients: the impact of two linked polymorphisms. Hum Mol Genet. (2004) 13:2369–76. doi: 10.1093/hmg/ddh274
12. Menon V, Povirk LF. XLF/Cernunnos: an important but puzzling participant in the nonhomologous end joining DNA repair pathway. DNA Repair. (2017) 58:29–37. doi: 10.1016/j.dnarep.2017.08.003
13. Guo C, Nakazawa Y, Woodbine L, Björkman A, Shimada M, Fawcett H, et al. XRCC4 deficiency in human subjects causes a marked neurological phenotype but no overt immunodeficiency. J Allergy Clin Immunol. (2015) 136:1007–17. doi: 10.1016/j.jaci.2015.06.007
14. Vera G, Rivera-Munoz P, Abramowski V, Malivert L, Lim A, Bole-Feysot C, et al. Cernunnos deficiency reduces thymocyte life span and alters the T cell repertoire in mice and humans. Mol Cell Biol. (2013) 33:701–11. doi: 10.1128/MCB.01057-12
15. Lescale C, Abramowski V, Bedora-Faure M, Murigneux V, Vera G, Roth DB, et al. RAG2 and XLF/Cernunnos interplay reveals a novel role for the RAG complex in DNA repair. Nat Commun. (2016) 7:10529. doi: 10.1038/ncomms10529
16. Zha S, Guo C, Boboila C, Oksenych V, Cheng HL, Zhang Y, et al. ATM damage response and XLF repair factor are functionally redundant in joining DNA breaks. Nature (2011) 469:250–4. doi: 10.1038/nature09604
17. Johansson P, Fasth A, Ek T, Hammarsten O. Validation of a flow cytometry-based detection of γ-H2AX, to measure DNA damage for clinical applications. Cytometry B Clin Cytom. (2017) 92:534–40. doi: 10.1002/cyto.b.21374
18. Viktorisson A, Mathew ST, Hammarsten O, Johansson P. A control for the day-to-day normalization of the flow cytometry γ-H2AX assay for clinical routine. Cytometry B Clin Cytom. (2018) 94:946–9. doi: 10.1002/cyto.b.21627
19. Kwan A, Hu D, Song M, Gomes H, Brown DR, Bourque T, et al. Successful newborn screening for SCID in the Navajo Nation. Clin Immunol. (2015) 158:29–34. doi: 10.1016/j.clim.2015.02.015
20. Yazdani R, Abolhassani H, Tafaroji J, Azizi G, Hamidieh AA, Chou J, et al. Cernunnos deficiency associated with BCG adenitis and autoimmunity: first case from the national Iranian registry and review of the literature. Clin Immunol. (2017) 183:201–6. doi: 10.1016/j.clim.2017.07.007
21. Slack J, Albert MH, Balashov D, Belohradsky BH, Bertaina A, Bleesing J, et al. Outcome of hematopoietic cell transplantation for DNA double-strand break repair disorders. J Allergy Clin Immunol. (2018) 141:322–8.e10. doi: 10.1016/j.jaci.2017.02.036
Keywords: XLF/Cernunnos, NHEJ1 mutation, DNA repair, severe combined immunodeficiency, lymphomagenesis, radiosensitive SCID (RS-SCID)
Citation: Recio MJ, Dominguez-Pinilla N, Perrig MS, Rodriguez Vigil-Iturrate C, Salmón-Rodriguez N, Martinez Faci C, Castro-Panete MJ, Blas-Espada J, López-Nevado M, Ruiz-Garcia R, Chaparro-García R, Allende LM and Gonzalez-Granado LI (2019) Extreme Phenotypes With Identical Mutations: Two Patients With Same Non-sense NHEJ1 Homozygous Mutation. Front. Immunol. 9:2959. doi: 10.3389/fimmu.2018.02959
Received: 05 October 2018; Accepted: 30 November 2018;
Published: 07 January 2019.
Edited by:
Andrew R. Gennery, Newcastle University, United KingdomReviewed by:
Silvia Clara Giliani, Università degli Studi di Brescia, ItalyJohn Bernard Ziegler, Sydney Children's Hospital, Australia
Copyright © 2019 Recio, Dominguez-Pinilla, Perrig, Rodriguez Vigil-Iturrate, Salmón-Rodriguez, Martinez Faci, Castro-Panete, Blas-Espada, López-Nevado, Ruiz-Garcia, Chaparro-García, Allende and Gonzalez-Granado. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Luis Ignacio Gonzalez-Granado, bHVpc2lnbmFjaW8uZ29uemFsZXpAc2FsdWQubWFkcmlkLm9yZw==
†These authors have contributed equally to this work and share senior authorship