Kissy Guevara-Hoyer1,2,3*
Jesús Fuentes-Antrás4,5
Eduardo de la Fuente-Muñoz1,2,3
Miguel Fernández-Arquero1,2,3
Fernando Solano6Pedro Pérez-Segura4Esmeralda Neves7,8
Alberto Ocaña4,5
Rebeca Pérez de Diego3,9
Silvia Sánchez-Ramón1,2,3- 1Cancer Immunomonitoring and Immuno-Mediated Pathologies Support Unit, IdSSC, Department of Clinical Immunology, San Carlos Clinical Hospital, Madrid, Spain
- 2Department of Clinical Immunology, IML and IdSSC, San Carlos Clinical Hospital, Madrid, Spain
- 3Department of Immunology, Ophthalmology and ENT, School of Medicine, Complutense University, Madrid, Spain
- 4Oncology Department, San Carlos Clinical Hospital, Madrid, Spain
- 5Experimental Therapeutics and Translational Oncology Unit, Medical Oncology Department, San Carlos University Hospital, Madrid, Spain
- 6Department of Hematology, General University Hospital Nuestra Señora del Prado, Talavera de la Reina, Spain
- 7Department of Immunology, Centro Hospitalar e Universitário do Porto, Porto, Portugal
- 8Unit for Multidisciplinary Research in Biomedicine (UMIB), Hospital and University Center of Porto, Porto, Portugal
- 9Laboratory of Immunogenetics of Human Diseases, IdiPAZ Institute for Health Research, Madrid, Spain
Common variable immunodeficiency (CVID) represents the largest group of primary immunodeficiencies that may manifest with infections, inflammation, autoimmunity, and cancer, mainly B-cell non-Hodgkin’s lymphoma (NHL). Indeed, NHL may result from chronic or recurrent infections and has, therefore, been recognized as a clinical phenotype of CVID, although rare. The more one delves into the mechanisms involved in CVID and cancer, the stronger the idea that both pathologies can be a reflection of the same primer events observed from different angles. The potential effects of germline variants on specific somatic modifications in malignancies suggest that it might be possible to anticipate critical events during tumor development. In the same way, a somatic alteration in NHL could be conditioning a similar response at the transcriptional level in the shared signaling pathways with genetic germline alterations in CVID. We aimed to explore the genomic substrate shared between these entities to better characterize the CVID phenotype immunodeficiency in NHL. By means of an in-silico approach, we interrogated the large, publicly available datasets contained in cBioPortal for the presence of genes associated with genetic pathogenic variants in a panel of 50 genes recurrently altered in CVID and previously described as causative or disease-modifying. We found that 323 (25%) of the 1,309 NHL samples available for analysis harbored variants of the CVID spectrum, with the most recurrent alteration presented in NHL occurring in PIK3CD (6%) and STAT3 (4%). Pathway analysis of common gene alterations showed enrichment in inflammatory, immune surveillance, and defective DNA repair mechanisms similar to those affected in CVID, with PIK3R1 appearing as a central node in the protein interaction network. The co-occurrence of gene alterations was a frequent phenomenon. This study represents an attempt to identify common genomic grounds between CVID and NHL. Further prospective studies are required to better know the role of genetic variants associated with CVID and their reflection on the somatic pathogenic variants responsible for cancer, as well as to characterize the CVID-like phenotype in NHL, with the potential to influence early CVID detection and therapeutic management.
1 Introduction
Common variable immunodeficiency (CVID) is the most prevalent symptomatic primary immunodeficiency (PID) and is characterized by an increased predisposition to recurrent infections resulting from the low production of antibodies against pathogens (1). CVID also entails an increased risk of autoimmune, inflammatory, and malignant diseases (2). Malignant neoplasms are a leading cause of death in CVID patients and may be the first clinical manifestation of the disease (3–5). Hematologic and gastric cancers are the most frequent, with an estimated increased risk of 10- to 47-fold with respect to the general population (6–8). The incidence of cancer diagnoses does not seem to be age-dependent, and the identification of robust clinical predictors and diagnostic biomarkers represents an unmet need for CVID patients (9).
Multiple carcinogenetic mechanisms have been proposed to operate in CVID. They can be classified as cell-intrinsic mechanisms, encompassing defects of DNA repair, T–B co-stimulation, immunoglobulin gene recombination (VDJ), class-switch recombination, and somatic hypermutation (SHM); and cell-extrinsic mechanisms, including inadequate immune responses that facilitate chronic infections, typically caused by EBV, HPV, and Helicobacter pylori (1, 10–13), defective immune surveillance against tumors, and dysbiosis and chronic inflammation, among others (13–16). To date, only a high level of clinical suspicion enables the early diagnosis of cancer in CVID patients (4, 9, 17–20). Conversely, patients with a new diagnosis of malignancy are not routinely screened for a potential underlying CVID (17, 21). Although fortunately low, the overall incidence of malignancies in CVID patients has increased in the last decades (9, 17, 22).
We have explored in this work the genetic crossovers between CVID and NHL at the somatic level, notwithstanding that a proportion of somatic variants may underlie at the germline level and may condition the accumulation of mutagenic variants in NHL. The two conditions are epidemiologically related in the direction of CVID toward NHL. Indeed, an extensive meta-analysis by Kiaee et al., deciphering the landscape of malignancy within CVID, encompassed 48 studies worldwide with a total cohort of 8,123 CIVD patients, of which 790 cases were associated with malignancy. NHL stood out with the highest prevalence (41%) in patients with CVID with associated malignancy (4). It should be noted that despite knowing this high predominance, the mechanisms that inter non-synonymous mutations in PIK3CD, twine these two pathologies remain an enigma (23–25). The association the other way around (NHL towards CVID) remains unknown. In NHL, B cells are the subject and the target of the disease, and indeed, a variable degree of immunodeficiency both in clonal and non-clonal cells is observed. We sought to determine as a working hypothesis that there must be gene networks associating CVID and lymphoma.
NHL entails a heterogeneous group of lymphocytic disorders ranging in aggressiveness from very indolent cellular proliferation to highly aggressive and rapidly proliferative processes (26). Even though it is true that several genes and pathways are involved in the development of lymphoma, not all of them are necessarily involved in the associated immunodeficiency. Indeed, there is a diverse clinical and immunological profile in NHL patients, which will determine different predispositions to recurrent infections and the need of starting replacement therapy with immunoglobulins to manage infectious diseases (26–28).
A single and specific defect is not known to cause CVID. To date, the multiple clinical and immunological features derived from the disease are not attributed, as in other PIDs, to the alteration of a single gene expression. CVID appears to be the result of several factors contributing to a defect in antibody production, where genetic, epigenetic, and environmental factors are involved (1, 21, 23, 29, 30). It has been described that around 20% of CVIDs are associated with a monogenic defect, so approximately the remaining 80% might account for a not yet identified gene or to the combined effects of digenic or oligogenic lesions (probably common and rare prevalence ones) and triggered by external factors (23, 31–33).
On the other hand, it should not be forgotten that there is a major pool of patients initially diagnosed with “CVID” who, after performing the genetic study, present a pathogenic variant of the established monogenic IEI (34–41), such as the case of NFKB1 haploinsufficiency or activated PI3K delta syndrome (described below). These patients are classified as “CVID phenotype” by the IUIS (40, 42), although they really could encompass a more complex genetic scenario (2, 43, 44). In the face of these findings, the differential diagnosis of both pathologies should always be examined.
The current challenge encompasses identifying the bases of the germline predisposition variants present in CVID and their role in cancer development. There is increasing evidence of how germline pathogenic variants can act as an oncogenic modifier, thus determining complimentary somatic variants necessary for the development of malignancy, as well as how these variants can behave as co-oncogenes through interactions with existing somatic pathogenic variants, conditioning tumorigenesis (45–49). The identification of molecular predictors of malignancy in CVID could enable the implementation of precision medicine in this population and substantially impact follow-up strategies and treatment decisions (21, 50). However, the design of adequately powered genome-wide studies in CVID is hindered by the limited casuistry and high costs of untargeted sequencing. The combination of publicly available genomic repositories and web tools for enrichment analyses may facilitate hypothesis generation.
We have undertaken an approximation between possible genetic interactions that justify a greater predisposition of this population to the development of malignancy. For this purpose, we performed an in-silico analysis to identify potential common genomic substrates of CVID and NHL that can favor the design of prospective studies with clinically and immunologically annotated cohorts.
2 Methods
2.1 Overall study design
This is a non-interventional study based on the analysis of real-world data contained in cBioPortal and using web-based analytical tools. The study was carried out in the Cancer Immunomonitoring and Immuno-Mediated Pathologies Support Unit of the Clinical Immunology Department, in close collaboration with the Oncology Department, San Carlos Clinical Hospital (HCSC), in Madrid, Spain.
2.2 Data collection and processing
2.2.1 Gene alterations associated with CVID
A systematic search of electronic databases was conducted to identify gene alterations recurrently associated with CVID following the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines (Supplementary Figure 1). English-language articles published in peer-reviewed journals and conference abstracts from 1979 to 30 August 2021 were identified in Medline, EMBASE, Cochrane Central Register of Controlled Trials (CENTRAL), and Cochrane Database of Systematic Reviews, using the terms “Common Variable” and “Gene” and excluding preclinical studies. Two investigators (KG-H, JF-A) extracted data from the included studies and a third investigator (SS-R) decided on unclear or conflicting data. The titles and abstracts were evaluated, and potentially relevant publications were retrieved in full. Overall, 3,420 articles were obtained from the database and manual searches, of which 311 articles were reviewed in full. A total of 50 of the most frequent genes associated with CVID with the main described pathogenic variants were selected for analysis (Supplementary Table 1).
2.2.2 Datasets, gene ontology, and functional analysis
We used data contained at cBioPortal (https://www.cbioportal.org, accessed in September 2021) from six studies in NHL (n = 1,309) to explore the distribution of genes with somatic variants occurring in the main CVID-associated germline genes associated with the pathogenic variant set. The pathogenic variants related to CVID and CVID phenotype found in the altered genes in NHL were verified through bibliographic sources as well as consulted in the ClinVar-NCBI tool (https://www.ncbi.nlm.nih.gov/clinvar/) (Supplementary Table 2). Genes such as KMT2C or the MSH family, despite not presenting clear CVID-associated pathological variables described in the literature, have been depicted as possibly harmful with a high frequency of heterozygotes compared with controls. These genes play an essential role in the associated CVID pathophysiology, such as somatic hypermutation or germline variation in cancer-susceptibility genes (51–54), for which they were considered in the analysis. Data from the VAF of the somatic pathogenic variants analyzed could not be fully retrieved. The biological functions of each gene were obtained using the 2018 Molecular_function Gene Ontology Terms through the publicly available EnrichR online platform (https://maayanlab.cloud/Enrichr/, accessed on 8 September 2021).
We used the online tool STRING (http://www.string-db.org) to construct interactome maps of genes of interest. The closer the local clustering coefficient is to 1, the more likely it is for the network to form clusters. The PPI enrichment p-value indicates statistical significance. Proteins are considered hubs when they have more interactions than the average. Co-occurrence analysis for gene alterations was evaluated using the cBioPortal online platform (http://www.cbioportal.org, accessed on 8 September 2021) (Supplementary Table 3).
This tool calculates the odds ratio (OR) for each pair of query genes, indicating the likelihood that the alterations for the two genes are co-occurrent in the selected cases, by the application of a Fisher’s exact test (statistical significance p < 0.05). Data analysis was performed using R (version 1.3.1093). The interaction network was mapped as default by the STRING database analysis.
3 Results
3.1 Study population
The six NHL datasets (Table 1) (45, 56–60) combined accounted for a total of 1,309 samples. The distribution according to sex was 50.8% men and 40.4% women (8.8% not specified). The mean age ranged from 65 to 70 years old, distributed from a minimum of 3 to a maximum of 93 years of age. The age of diagnosis was not available for analysis.
Table 1 NHL datasets available in cBioPortal (www.cbioportal.org; accessed in September 2021).
3.2 CVID-associated gene set in NHL
The prevalence of alterations in the main 50 CVID-associated genes and their genetic variants was investigated in the NHL mutation datasets. A total of 323 from 1,309 samples (25%) harbored at least one alteration, ranging from 21.38% to 62.50% according to the series. The main genetic alterations associated with CVID found in the different studies and the proportions and numbers of altered samples are shown in Table 2.
Table 2 Main CVID-associated genes and their alteration in the NHL datasets.
The prevalence of gene alterations in NHL was compared with that reported as genetic alterations in CVID (50, 61). Approximately up to 25% of patients with CVID carry a germline alteration in non-consanguineous populations (31). However, the prevalence of specific gene alterations remains largely unknown, particularly those categorized as rare functional variants and those which independently do not exhibit a causal relationship (Table 3) (21, 61, 62). The most prevalent gene alterations in CVID were PIK3CD (2.6%), LRBA (2.6%), and NFKB2 (0.5%), while the most frequent in NHL were PIK3CD (6%), KMT2C (5%), and STAT3 (4%). Supplementary Table 2 shows the genes with pathogenic variables associated with the CVID-like phenotype and their frequency found in NHL samples.
Table 3 Frequencies of the main CVID-associated gene alterations in the NHL and CVID populations.
3.3 Characterization of the functional networks, gene ontology, and biological processes of the genes altered in both NHL and CVID
In order to unveil the putative functional connection between the most prevalent gene alterations in both NHL and CVID, we interrogated the STRING database and generated protein–protein interaction maps adjusted by the type and strength of the interaction (Figure 1). A confidence threshold of 0.7 was set between the nodes (PPI enrichment p-value: <1.0e−16), with 14 nodes lacking predicted interactions. Notably, PIK3R1 appeared as a central node connecting the PI3K pathway and downstream mediators with different immune mediators such as STAT3, NFKB2, CTLA4, or CD19.
Figure 1 Protein–protein interaction map displaying the significant functional network integrated by the selected genes involved in NHL and CVID. Known interactions are displayed in blue and pink lines. Predictive interactions include the following: green line, gene neighborhood; red line, gene fusions; blue line, gene co-occurrence; yellow line, text mining; black line, co-expression; and gray line, protein homology. Connections are filtered by the highest confidence (0.7). PPI enrichment p-value: < 1.0e−16.
To further investigate the functional characterization of the set of CVID-associated genes, we interrogated the EnrichR database and found that inflammation, immune surveillance, and defective DNA repair appeared as significantly enriched biological processes (Table 4), with recurrent functions in cytokine and TNF signaling.
Table 4 Relationships organized in biological pathways and processes mainly related to CVID-associated genes expressed in NHL (sorted by p-value; p < 0.05).
3.4 Co-occurrence of alterations in CVID-associated genes in NHL
We next analyzed the co-occurrence of genetic alterations in the CVID-associated gene set. Sixteen pairs presented statistical significance co-occurrence (Table 5). Regarding the relationships organized in biological pathways and processes related to CVID-associated genes expressed in NHL, the proteins encoded by 6 of the 16 significant pairs were found to contribute jointly to a shared function. In contrast, 10 of the 16 pairs that represented independent processes did not associate with common pathways. No mutual exclusivity was found in the analysis.
Table 5 Co-occurrence of mutations in the analyzed population calculated by the odds ratio method in cBioPortal.
4 Discussion
Unraveling the molecular interplay underlying CVID is a difficult task due to the lack of consistent diagnostic protocols and the need of high clinical suspicion. Patients are often studied late after recurrent infectious episodes or associated complications including malignancy. In particular, NHL is the most prevalent malignancy in CVID patients, and its concurrence may lead to worse outcomes through either poor treatment responses and recurrent infectious complications (3, 19, 22, 24, 63). Therefore, understanding how CVID and NHL are interrelated is an unmet need that requires preliminary data to substantiate future prospective endeavors.
Similar to other malignancies, NHL B cells coordinate signaling pathways in the tumor microenvironment, promoting tumorigenesis and generating a state of immunosuppression that favors the tumor growth and progression (26, 64). B cells are the main players in both NHL and CVID, and unlike solid tumors, an alteration in both somatic and germline might condition a similar response at the transcriptional level in the shared signaling pathways, generating, in this setting, a somatic profile of immunodeficiency in NHL similar to that derived from the germline of CVID (65–67). This phenomenon may be a reflection of specific somatic pathogenic variants in the aberrant clonal B cells responsible for orchestrating key processes shared between both pathologies, such as defects in maturation and differentiation of the B lymphocyte, T–B co-stimulation, immunoglobulin gene recombination (VDJ), class-switch recombination, and somatic hypermutation (45, 52, 58, 68–71).
Pathway enrichment studies have determined that the molecular mechanisms involved in both Hodgkin’s and non-Hodgkin’s lymphoma were more enriched in those pathways associated with immune response. The vital role of pathogenesis and immune escape mechanisms and the development of an immunosuppressive tumor microenvironment have been highlighted by these analyses (64, 68, 72, 73). By way of illustration, the paths associated with the members of the Bcl-2 family, which are proteins that play an essential role in regulating cell apoptosis, survival, and proliferation, are highly expressed in most NHL (74). BCL2 is an antiapoptotic factor, and its deregulation in NHL is associated with the constitutive activation of NF-κB, which in turn influences the therapeutic response and prognosis of these patients (75, 76). On the other hand, overexpression of BCL-xL has been reported in 80% of all NHL, relating it to inhibition of cell death and, therefore, the development of malignancy (74).
Bcl‐2 and Bcl‐xl share functional homology, and their transcription is intimately regulated by NF‐κB activity (35, 77, 78). NF-κB1 loss-of-function (LoF) has been associated with an impaired function of the Bcl-2 protein in memory B cells of CVID patients, which would predispose to apoptosis and to defective differentiation and maturation of immunoglobulin-producing B lymphocytes (a typical feature of CVID) (36, 77). The latter mechanism would be contrary to what occurs in NHL, where apoptosis would be suppressed contributing to cell proliferation and malignancy (74, 75). Some series describe NF-κB as the most common monogenic cause of CVID in Europeans (34, 36). However, studies such as those by Li et al. denote that NFKB1 variants are common in healthy individuals and cannot necessarily be considered causal CVID (35).
Several signals that lead to activation of the NF-κB pathway to regulate B-cell survival in mature B cells highlighted the ligands BAFF (B-cell-activating factor of the tumor necrosis factor family) and APRIL (a proliferating inducing ligand), both joining to TACI (transmembrane activator and CAML interaction) (79). The overexpression of BAFF, APRIL, or their receptor TACI is noted in hematological malignancies such as B-cell NHL (80, 81). This would contribute to the survival and prognosis of tumor cells (37, 81, 82). TACI/APRIL function is also related to the class switch recombination of immunoglobulins (38, 39). Several variants have been considered CVID causative or modifier (83, 84), depending on the occurrence of similar variants with uncertain significance in the general population (84–87).
Since the pathways associated with tumorigenesis (such as cell death, cell proliferation, and response to stress) are shared by a large number of genes with very different or even opposite roles (69, 88), it is complicated, if not impossible, to draw a perfect line that separates the immunological pathways involved in the development of lymphoma from those of underlying immunodeficiency. The ways involved could be similar but up- or downregulated by the penetrance and expressivity of the differential genes (89–94). Our study provides first evidence of the common genomic landscape of CVID and NHL. By means of an extensive systematic search followed by an in-silico analysis of publicly available genomic repositories, we have identified a set of genes recurrently altered in both diseases and therefore suggestive of common functional features of immunodeficiency. Among the CVID-associated genetic variants, non-synonymous mutations in PIK3CD, KMT2C, and STAT3 were the most prevalent in NHL. Functional enrichment analysis revealed a robust centralizing role for PIK3R1 and downstream effectors of the PI3K pathway.
PIK3CD exerts a key role in the activation of signaling cascades involved in cell proliferation and survival (95). Likewise, it mediates immune responses associated with the development, migration, and function of B lymphocytes and TCR signaling at the immune synapse and participates in the activation of the NK cell and in neutrophil chemotaxis, as well as in the production of cytokines in response to TLR4 and TLR9, necessary for antigenic presentation (96).
The second candidate KMT2C represents a catalytic subunit of the MLL2/3 coactivator complex, specific nuclear receptors for epigenetic transcriptional activation (54, 97). Pathogenic variants in KMT2C have been linked to leukemogenesis and disorders of cell development (98). So far, KMT2C does not present clear CVID-associated pathological variables described in the literature; nevertheless, this gene has been expressed as possibly harmful with a high frequency of heterozygotes compared with controls, where their function plays an essential role in the associated CVID pathophysiology such as germline variation in cancer susceptibility (49, 99).
Thirdly, an interesting finding was the prevalence of the STAT-3 gene at 4% in NHL, but unknown for CVID (63, 100). STAT-3 is a signal transducer and activator of transcription that mediates cellular responses associated with interleukins and growth factors (101). Among its multiple functions, STAT-3 intervenes in the cell cycle by inducing the expression of critical genes for the progression from G1 to S phase, is involved in the response to activated FGFR1-4 and anti-inflammatory response, and modulates differentiation toward Th17 or regulatory (TReg) cells through IL-6 signaling cascade and on the promoters of several acute-phase protein genes. In a general framework, hematological malignancies associated with STAT-3 alteration are linked to gain of function (GOF) (102). STAT genes are relevant for lymphocyte survival and tumorigenesis (100–104). STAT-3 mutations are rare in B-cell NHL. However, STAT-3 GOF mutations have been associated with increased STAT-3 phosphorylation and transcriptional activity, which may foster malignant cell survival and proliferation (101). Our study detected nine pathogenic variants associated with the GOF of STAT-3 that have been described to be associated with an immunodeficiency phenotype (105, 106). In the same way, other authors have related GOF mutations in STAT-3 with several immunodeficiency phenotypes, considering it as a candidate causative gene of the CVID phenotype (63, 100).
Of note, the loss-of-function mutations in PIK3R1 and PTEN lead to a similar clinical phenotype to that of GOF mutations in the PIK3CD gene, characterized by decreased B lymphocytes and naïve T lymphocytes, lymphoproliferation, and autoimmunity (107, 108). These genes are within our top 7 CVID-expressed genetic variants in NHL mutations and have been previously associated with lymphoma predisposition (40, 109). Pathogenic variants in the PTEN gene are a known substrate of the Cowden syndrome. Although atypical, its presentation as a CVID phenotype has been reported (110, 111). On the other hand, patients with pathogenic variants that activate PI3KR1 signaling have been associated with antibody deficiency partly due to the altered function of p85α (112). Low B-lymphocyte counts can also result from severely impaired B-cell development and differentiation (107). The degree of involvement is associated with the type of inheritance, the more profound being in homozygous patients (agammaglobulinemia) and variable behavior (hypogammaglobulinemia) in patients with heterozygous mutations at the PIK3R1 splice site (107). One patient with heterozygous PIK3R1 mutation presented the CVID phenotype, and this gene is considered causative of CVID phenotype immunodeficiency (63, 113, 114).
Patients with activated PI3Kdelta syndrome (APDS) have an increased risk of both HL and NHL. The incidence of lymphoma is as high as 13%–30% (115, 116), and a vast majority is related to EBV infection (117, 118). Although the immune escape strategies in lymphoma may vary between individuals (119), the most significant known stimuli so far related to the pathogenesis of B-cell lymphoma in APDS are related to the upregulation of the mTOR signaling pathway with the downstream effector of PI3K/AKT, increased production of transcription factors involved in the process of apoptosis, affecting the regulation of the cell cycle conducting to uncontrolled cell survival and malignant transformation (120, 121). For many years, APDS was defined as a CVID-like syndrome, especially due to its clinical similarity to that of monogenic CVID forms (107, 115, 120). However, APDS patients are more predisposed to viral infections than CVID, highlighting EBV infection (117, 122). The clinical features of cellular immunodeficiency play an essential role in APDS, such as the increased senescent T lymphocytes with shortened telomeres, especially of APDS type I (APDS1) (96, 123). The differentiation of PI3K signaling defect between both pathologies is crucial.
Our aim was to investigate the plausible common genomic grounds between CVID phenotype immunodeficiency in NHL. Therefore, an exploratory analysis of germline CVID-related gene alterations was performed in NHL tumor samples. A key immune system feature is the recombination of the (VDJ) gene segment in the generation of T- and B-cell antigen receptor (TCR and BCR) repertoire, as well as SHM of B cells in the generation of immunoglobulin isotypes (124). These processes are orchestrated and tightly regulated by the expression of specific genes, as well as their resulting proteins, which also share a function within the DNA and transcription repair mechanisms (70, 125), the alteration of the latter being widely known as cancer-predisposing factors (10, 71, 126–128).
As the availability of next-generation sequencing data increases, there is more evidence of the interrelationship of somatic and germline genetic conditions associated with immunodeficiency and cancer (46). Somatic pathogenic variants can reflect or be a consequence of alterations in germinal variants (45, 47, 65). Several studies have recently explored how germline variants can act as an oncogenic modifier, thus determining complementary somatic pathogenic variants necessary for the development of malignancy. Similarly, germline variants not directly associated with oncogenesis can behave as co-oncogenes through interactions with existing somatic pathogenic variants, thus enhancing tumorigenesis (47, 65). This phenomenon suggests that the inherited germline variants would govern where and how malignancy will develop. Germline variants of somatic alteration genes in NHL may themselves predispose to NHL in patients with CVID (3, 45, 113, 129). Since the clinical phenotype caused by germline variants can vary depending on the genetic alterations and the patient, the combination of variants could be pathogenic and lead to an increased susceptibility to developing NHL (11, 19, 58, 65, 69, 129).
The pioneer PanCanAtlas Germline Working Group has investigated the pathogenic germline variants in more than 10,000 adults within 33 types of cancer, including diffuse large B-cell lymphoma. Although the results are preliminary and the data are not yet public, the results are promising to enlighten the basis of germline predisposition variants and their role in cancer (130). Seidel et al. performed the analysis of germline variants in childhood cancers, suggesting that the development of most cancers seemed to be more related to a cell-intrinsic defect of the immune system in the detection, control, and elimination of cancer cells, rather than being directly caused by extrinsic associated factors (131). In a recent study, Hauck et al. explained how the germline genetic variants associated with IEI behave as oncogenes when they arise as somatic pathogenic variants and cause specific cancerous entities, mainly CARD11, IKZF1, GATA2, PMS2, MSH6, etc. alterations (46).
In the work of Lincoln et al., 34% of the patients (207/608) presented GPV and second neoplasms and 15% (32/207) of them carried GPV associated with specific recommendations for detection or risk reduction for their subsequent cancers (132). Given that CVID patients are at theoretical high risk for developing second neoplasms (4, 6, 10), screening of CVID-associated genes at malignancy diagnosis might modify current protocols in order to improve follow-up and prognosis.
An interesting study by Liu et al. proposed 172 possible CVID candidate genes functionally similar to known CVID genes based on their interactions and biological distance. Liu et al. studied their differential expression finding upregulation and downregulation in patients compared to controls (133). This suggests that CVID could behave as a mixed model with the presence of punctual alterations due to cumulative effects of polygenic determinants, gene–gene interactions, and/or regulatory variation in non-coding regions detrimental to a particular immune pathway rather than as an effect derived from a specific point of mutation. This overexpression could be fundamental as a complementary co-oncogene factor in the development of malignancy in CVID patients (134, 135).
Our data show that up to one-fifth of NHL patients had similar CVID-related genetic pathogenic variants in the somatic line, some with a frequency up to 6%, as is the case of PIK3CD. Several of these genes were significantly associated with different important pathways of the immune response, inflammation, and cell repair (45, 50, 56, 98, 136). CVID-associated genes expressed in NHL patients are also involved in DNA repair, genome integrity, immunosurveillance, and predisposition to cancer (7, 18, 113, 137). Consequently, these genes could also be described as NHL-associated GPV, as well as with other malignancies (136, 138). Interestingly, our data help to identify and prioritize those pleiotropic genes that could simultaneously affect several relevant immunological pathways. Two of the clearest examples in our study were the NFKB1 and NFKB2 genes that were found to be involved in pathways such as PKMT methylation of histone lysines, DEx/H-box helicase activation of type I IFN and inflammatory cytokine production, RIP-mediated NF-κB activation via ZBP1, and interleukin-1 processing, among others. Other examples are MLH1, MSH2, MSH6, and PMS2, genes involved in DNA damage/repair that are highly involved in immunodeficiency, NHL, constitutional mismatch repair deficiency (CMMRD), and Lynch syndrome and also associated with colorectal and endometrial cancer (129, 139, 140).
Another intriguing result of our study is the relationships among organized biological pathways and CVID-associated genetic variant processes expressed as somatic pathogenic variants in NHL, mainly the core role that PIK3R1 seems to have as an interconnector of altered signaling pathways (Supplementary Figure 2). PIK3R1 is the seventh in the gene ranking for NHL and CVID, with a prevalence of 1.3% in NHL and 0.48% in CVID (61). However, according to the predicted protein–protein interaction network, it may centralize relevant immune signaling pathways associated with immunodeficiency and malignancy (107, 112, 141). These influential hub genes interact with other nodes and thus exert a role of checkpoint alterations at the heart of both disorders (107, 112, 141).
Patients with APDS have a greater predisposition to B-cell lymphoma, for which the PIK3R1 involvement may be considered a risk factor for carcinogenesis (107). Moreover, the spectrum of malignancy could be broader, since the involvement of p85α can also affect the function of p110α, p110β, and PTEN (141, 142), leading to hyperactivation of PI3K signaling in other cell types with a greater predisposition to malignancy (107, 141, 142).
In CVID published cohorts, the estimated prevalence of the direct mutation in PIK3R1 is 0.481%, with its associated phenotype being PIK3CD which is the most frequent pathogenic variant in both CVID and NHL with the prevalence of 2.674% and 6.0%, respectively (61). In the same way, as previously mentioned, PIK3R1 seems to be the common link among CVID-related genes in NHL, further representing a frequently associated gene alteration in many tumors (141). PIK3R1 is altered in 3.62% of solid cancer patients and 1.52% of lymphoma patients (112, 143). In our match with NHL, the prevalence of PIK3R1 was 1.3%, slightly lower than that described in previous studies (112, 143). Our results point to PIK3R1 as a hotspot in the relationship between CVID and NHL, not just directly associated with its mutation, but rather with a defect in the common signaling pathways for which it is crucial (63, 96, 107, 113, 114, 141). This hotspot would not only condition the development of malignancy but also influence the immunodeficiency phenotype linked with CVID. These findings encourage further investigation on the pathogenic status and role of PIK3R1 in CVID patients.
Our study is based on an approach to better understand the genetic link associated with signaling pathways shared between both pathologies (NHL and CVID) in a population where malignancy is a leading cause of mortality. Our preliminary data shown here align with the hypothesis that both conditions are indeed genetically bidirectionally related, and even if all the variants found were exclusively somatic variants within NHL cells, they could induce the associated immunodeficiency. It is a shame not to be able to retrieve the VAF from the somatic pathogenic variants analyzed. These data would have implications and relevance in understanding the behavior of somatic variants in malignancy as a reflection of immunodeficiency and vice versa.
We look forward to validating our results in independent large cohorts of immunodeficiency and cancer, as well as the detailed description of pathogenic germline variants in tumors. More evidence is required such as studies of long cohorts of immunodeficiency and cancer, as well as the detailed description of pathogenic germline variants in tumors.
Nevertheless, the continuous progress in genetic sequencing and in the availability of cancer genomic datasets opens new avenues to deepen our understanding of the human immune system and its association with cancer. Our results suggest the putative gene alterations potentially connecting CVID and NHL. Whether these represent a common substrate of CVID phenotype immunodeficiency or rather a risk predisposition for NHL of CVID patients remains unknown and requires prospective studies. These studies, however, may not need to rely on extensive whole-exome sequencing but on customized targeted panels or dPCR approaches leveraging our data, thus facilitating the generation of new evidence in cohorts with a limited number of patients. Beyond conventional cohort studies, the creation of large-scale repositories with comprehensive clinical and immunological annotation (i.e., flow cytometry, immunohistochemistry, RNA-seq immunological signatures) would be a critical step in bridging the elusive gap between immunodeficiency and cancer.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Ethics statement
The studies involving human participants were reviewed and approved by Hospital Clínico San Carlos institutional research Ethics Committee (21/504-E). Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.
Author contributions
KG-H, JF-A, RPD, and SS-R contributed to the design of the study and analyzed the data. AO, RPD, PP-S, FS, EN, and MF-A contributed to the analysis of the results and critical review. EF-M contributed to the data collection and designed the tables. KGH and SSR wrote the manuscript. KG-H designed the figures. All authors have read and agreed to the published version of the manuscript.
Funding
KG-H was supported by The European Social Fund (ESF) through a Río Hortega Grant for Health Research Projects by the Carlos III Health Institute (ISCIII) (CM20/00098). JF-A was supported by a Fundación Sociedad Española de Oncología Médica (SEOM)/EUSA Pharma grant.
Acknowledgments
We appreciate the insightful and constructive comments made by the reviewers that allowed us to improve the quality of our manuscript.
Conflict of interest
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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2022.937872/full#supplementary-material
Glossary
APDS: activated PI3Kdelta syndrome
APRIL: a proliferating inducing ligand
ATP6AP1: ATPase H+ Transporting Accessory Protein 1
BAFF: B-cell-activating factor of the tumor necrosis factor family
CD19: cluster of differentiation 19
CD20: cluster of differentiation 20
CD21: cluster of differentiation 21
CD81: cluster of differentiation 81
CENTRAL: Cochrane Central Register of Controlled Trials
CMMRD: constitutional mismatch repair deficiency
CVID: common variable immunodeficiency
DOCK8: dedicator of cytokinesis 8
DLBCL: diffuse large B-cell lymphoma
EBV: Epstein–Barr virus
GSEA: gene set enrichment analysis
GPV: germline pathogenic variants
HPV: human papillomavirus
ICOS: inducible T-cell costimulatory
IEI: inborn errors of immunity
IKZF1: IKAROS family zinc finger 1
IRF2BP2: interferon regulatory factor 2-binding protein 2
LRBA: lipopolysaccharide (LPS)-responsive and beige-like anchor protein
MOGS: mannosyl-oligosaccharide glucosidase
NFKB1: nuclear factor kappa B subunit 1
NFKB2: nuclear factor kappa B subunit 2
NGS: next-generation sequencing
NHL: non-Hodgkin’s lymphoma
OR: odds ratio
PID: primary immunodeficiency
PIK3CD: phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit delta
PIK3R1: phosphoinositide-3-kinase regulatory subunit 1
PRISMA: Preferred Reporting Items for Systematic Reviews and Meta-Analyses
PRKCD: protein kinase C delta
PTEN: phosphatase and tensin homolog
SHM: somatic hypermutation
STAT3: signal transducer and activator of transcription 3
STXBP2: syntaxin-binding protein 2
TACI: transmembrane activator and CAML interaction
TCGA: The Cancer Genome Atlas
TRNT1: TRNA nucleotidyl transferase 1
TNFSF12: TNF-related weak inducer of apoptosis (TWEAK), TNF superfamily member 12
References
1. Yazdani R, Habibi S, Sharifi L, Azizi G, Abolhassani H, Olbrich P, et al. Common variable immunodeficiency: Epidemiology, pathogenesis, clinical manifestations, diagnosis, classification, and management. J Investig Allergol Clin Immunol (2020) 30:14–34. doi: 10.18176/jiaci.0388
2. Chapel H, Lucas M, Patel S, Lee M, Cunningham-Rundles C, Resnick E, et al. Confirmation and improvement of criteria for clinical phenotyping in common variable immunodeficiency disorders in replicate cohorts. J Allergy Clin Immunol (2012) 130:1197–98.e9. doi: 10.1016/j.jaci.2012.05.046
3. Kralickova P, Milota T, Litzman J, Malkusova I, Jilek D, Petanova J, et al. CVID-associated tumors: Czech nationwide study focused on epidemiology, immunology, and genetic background in a cohort of patients with CVID. Front Immunol (2019) 9:3135. doi: 10.3389/fimmu.2018.03135
4. Kiaee F, Azizi G, Rafiemanesh H, Zainaldain H, Sadaat Rizvi F, Alizadeh M, et al. Malignancy in common variable immunodeficiency: A systematic review and meta-analysis. Expert Rev Clin Immunol (2019) 15:1105–13. doi: 10.1080/1744666X.2019.1658523
5. What is the increased risk for malignancy in common variable immunodeficiency (CVID)?. Available at: https://www.medscape.com/answers/885935-190810/what-is-the-increased-risk-for-malignancy-in-common-variable-immunodeficiency-cvid (Accessed June 12, 2020).
6. Kersey JH, Shapiro RS, Filipovich AH. Relationship of immunodeficiency to lymphoid malignancy. Pediatr Infect Dis J (1988) 7:S10–12.
7. Verhoeven D, Stoppelenburg AJ, Meyer-Wentrup F, Boes M. Increased risk of hematologic malignancies in primary immunodeficiency disorders: Opportunities for immunotherapy. Clin Immunol (2018) 190:22–31. doi: 10.1016/j.clim.2018.02.007
8. Gullo I, Costa C, Silva SL, Ferreira C, Motta A, Silva SP, et al. The dysfunctional immune system in common variable immunodeficiency increases the susceptibility to gastric cancer. Cells (2020) 9:E1498. doi: 10.3390/cells9061498
9. Haas OA. Primary immunodeficiency and cancer predisposition revisited: Embedding two closely related concepts into an integrative conceptual framework. Front Immunol (2019) 9:3136. doi: 10.3389/fimmu.2018.03136
10. Bomken S, van der Werff Ten Bosch J, Attarbaschi A, Bacon CM, Borkhardt A, Boztug K, et al. Current understanding and future research priorities in malignancy associated with inborn errors of immunity and DNA repair disorders: The perspective of an interdisciplinary working group. Front Immunol (2018) 9:2912. doi: 10.3389/fimmu.2018.02912
11. Gangemi S, Allegra A, Musolino C. Lymphoproliferative disease and cancer among patients with common variable immunodeficiency. Leuk Res (2015) 39:389–96. doi: 10.1016/j.leukres.2015.02.002
12. Ballow M. Approach to the patient with recurrent infections. Clin Rev Allergy Immunol (2008) 34:129–40. doi: 10.1007/s12016-007-8041-2
13. Baldovino S, Montin D, Martino S, Sciascia S, Menegatti E, Roccatello D. Common variable immunodeficiency: Crossroads between infections, inflammation and autoimmunity. Autoimmun Rev (2013) 12:796–801. doi: 10.1016/j.autrev.2012.11.003
14. Hauck F, Gennery AR, Seidel MG. Editorial: The relationship between cancer predisposition and primary immunodeficiency. Front Immunol (2019) 10:1781. doi: 10.3389/fimmu.2019.01781
15. Arieira C, Dias de Castro F, Moreira MJ, Cotter J. Common variable immunodeficiency-associated inflammatory enteropathy: The new era of biological therapy. GE - Port J Gastroenterol (2018) 25:322–6. doi: 10.1159/000486448
16. Belkaid Y, Hand T. Role of the microbiota in immunity and inflammation. Cell (2014) 157:121–41. doi: 10.1016/j.cell.2014.03.011
17. Mayor P, Eng K, Singel K, I. Abrams S, Odunsi K, Moysich K, et al. Cancer in primary immunodeficiency diseases: Cancer incidence in the united states immune deficiency network registry. J Allergy Clin Immunol (2017) 141:1028-35. doi: 10.1016/j.jaci.2017.05.024
18. Salavoura K, Kolialexi A, Tsangaris G, Mavrou A. Development of cancer in patients with primary immunodeficiencies. Anticancer Res (2008) 28:1263–9.
19. Yakaboski E, Fuleihan RL, Sullivan KE, Cunningham-Rundles C, Feuille E. Lymphoproliferative disease in CVID: A report of types and frequencies from a US patient registry. J Clin Immunol (2020) 40:524–30. doi: 10.1007/s10875-020-00769-8
20. Wehr C, Houet L, Unger S, Kindle G, Goldacker S, Grimbacher B, et al. Altered spectrum of lymphoid neoplasms in a single-center cohort of common variable immunodeficiency with immune dysregulation. J Clin Immunol (2021) 41:1250–65. doi: 10.1007/s10875-021-01016-4
21. Abolhassani H, Hammarström L, Cunningham-Rundles C. Current genetic landscape in common variable immune deficiency. Blood (2020) 135:656–67. doi: 10.1182/blood.2019000929
22. Quinti I, Agostini C, Tabolli S, Brunetti G, Cinetto F, Pecoraro A, et al. Malignancies are the major cause of death in patients with adult onset common variable immunodeficiency. Blood (2012) 120:1953–4. doi: 10.1182/blood-2012-05-431064
23. Chapel H, Cunningham-Rundles C. Update in understanding common variable immunodeficiency disorders (CVIDs) and the management of patients with these conditions. Br J Haematol (2009) 145:709–27. doi: 10.1111/j.1365-2141.2009.07669.x
24. Resnick ES, Moshier EL, Godbold JH, Cunningham-Rundles C. Morbidity and mortality in common variable immune deficiency over 4 decades. Blood (2012) 119:1650–7. doi: 10.1182/blood-2011-09-377945
25. Mayor PC, Eng KH, Singel KL, Abrams SI, Odunsi K, Moysich KB, et al. Cancer in primary immunodeficiency diseases: Cancer incidence in the United States Immune deficiency network registry. J Allergy Clin Immuno (2018) 141(3):1028–35. doi: 10.1016/j.jaci.2017.05.024
26. Elhelbawy NG, Nassar AAH, Eltorgoman AEA, Saber SM, Badr EA. Immunological microenvironment gene expression in patients with diffuse Large b cell non Hodgkin lymphoma. Biochem Biophys Rep (2020) 21:100731. doi: 10.1016/j.bbrep.2020.100731
27. Casulo C, Maragulia J, Zelenetz AD. Incidence of hypogammaglobulinemia in patients receiving rituximab and the use of intravenous immunoglobulin for recurrent infections. Clin Lymphoma Myeloma Leuk (2013) 13:106–11. doi: 10.1016/j.clml.2012.11.011
28. Ansell SM. Harnessing the power of the immune system in non-Hodgkin lymphoma: Immunomodulators, checkpoint inhibitors, and beyond. Hematol Am Soc Hematol Educ Program (2017) 2017:618–21. doi: 10.1182/asheducation-2017.1.618
29. Rae W. Indications to epigenetic dysfunction in the pathogenesis of common variable immunodeficiency. Arch Immunol Ther Exp (Warsz) (2017) 65:101–10. doi: 10.1007/s00005-016-0414-x
30. Li J, Wei Z, Li YR, Maggadottir SM, Chang X, Desai A, et al. Understanding the genetic and epigenetic basis of common variable immunodeficiency disorder through omics approaches. Biochim Biophys Acta (2016) 1860:2656–63. doi: 10.1016/j.bbagen.2016.06.014
31. Ameratunga R, Woon S-T. Perspective: Evolving concepts in the diagnosis and understanding of common variable immunodeficiency disorders (CVID). Clin Rev Allergy Immunol (2020) 59:109–21. doi: 10.1007/s12016-019-08765-6
32. Piqueras B, Lavenu-Bombled C, Galicier L, Bergeron-van der Cruyssen F, Mouthon L, Chevret S, et al. Common variable immunodeficiency patient classification based on impaired b cell memory differentiation correlates with clinical aspects. J Clin Immunol (2003) 23:385–400. doi: 10.1023/a:1025373601374
33. Kutukculer N, Azarsiz E, Karaca NE, Ulusoy E, Koturoglu G, Aksu G. A clinical and laboratory approach to the evaluation of innate immunity in pediatric CVID patients. Front Immunol (2015) 6:145. doi: 10.3389/fimmu.2015.00145
34. Fliegauf M, Krüger R, Steiner S, Hanitsch LG, Büchel S, Wahn V, et al. A pathogenic missense variant in NFKB1 causes common variable immunodeficiency due to detrimental protein damage. Front Immunol (2021) 12:621503. doi: 10.3389/fimmu.2021.621503
35. Li J, Lei W-T, Zhang P, Rapaport F, Seeleuthner Y, Lyu B, et al. Biochemically deleterious human NFKB1 variants underlie an autosomal dominant form of common variable immunodeficiency. J Exp Med (2021) 218:e20210566. doi: 10.1084/jem.20210566
36. Tuijnenburg P, Lango Allen H, Burns SO, Greene D, Jansen MH, Staples E, et al. Loss-Of-Function nuclear factor κb subunit 1 (NFKB1) variants are the most common monogenic cause of common variable immunodeficiency in europeans. J Allergy Clin Immunol (2018) 142:1285–96. doi: 10.1016/j.jaci.2018.01.039
37. Gupta M, Ansell SM, Witzig TE, Ziesmer SC, Cerhan JR, Dillon SR, et al. APRIL-TACI interactions mediate non-Hodgkin lymphoma b cell proliferation through akt regulated cyclin D1 and P21. Blood (2007) 110:3585. doi: 10.1182/blood.V110.11.3585.3585
38. Romberg N, Chamberlain N, Saadoun D, Gentile M, Kinnunen T, Ng YS, et al. CVID-associated TACI mutations affect autoreactive b cell selection and activation. J Clin Invest (2013) 123:4283–93. doi: 10.1172/JCI69854
39. Salzer U, Grimbacher B. TACI deficiency - a complex system out of balance. Curr Opin Immunol (2021) 71:81–8. doi: 10.1016/j.coi.2021.06.004
40. Bousfiha A, Jeddane L, Picard C, Al-Herz W, Ailal F, Chatila T, et al. Human inborn errors of immunity: 2019 update of the IUIS phenotypical classification. J Clin Immunol (2020) 40:66–81. doi: 10.1007/s10875-020-00758-x
41. Kienzler A -K., Hargreaves CE, Patel SY. The role of genomics in common variable immunodeficiency disorders. Clin Exp Immunol (2017) 188:326–32. doi: 10.1111/cei.12947
42. Tangye SG, Al-Herz W, Bousfiha A, Cunningham-Rundles C, Franco JL, Holland SM, et al. Human inborn errors of immunity: 2022 update on the classification from the international union of immunological societies expert committee. J Clin Immunol (2022) 24:1–35. doi: 10.1007/s10875-022-01289-3
43. Guevara-Hoyer K, Vasconcelos J, Marques L, Fernandes AA, Ochoa-Grullón J, Marinho A, et al. Variable immunodeficiency study: Evaluation of two European cohorts within a variety of clinical phenotypes. Immunol Lett (2020) 223:78–88. doi: 10.1016/j.imlet.2020.03.006
44. de Valles-Ibáñez G, Esteve-Solé A, Piquer M, González-Navarro EA, Hernandez-Rodriguez J, Laayouni H, et al. Evaluating the genetics of common variable immunodeficiency: Monogenetic model and beyond. Front Immunol (2018) 9:636. doi: 10.3389/fimmu.2018.00636
45. Reddy A, Zhang J, Davis NS, Moffitt AB, Love CL, Waldrop A, et al. Genetic and functional drivers of diffuse Large b cell lymphoma. Cell (2017) 171:481–94.e15. doi: 10.1016/j.cell.2017.09.027
46. Barmettler S, Sharapova SO, Milota T, Greif PA, Magg T, Hauck F. Genomics driving diagnosis and treatment of inborn errors of immunity with cancer predisposition. J Allergy Clin Immunol Pract (2022) S2213-2198(22):00290–2. doi: 10.1016/j.jaip.2022.03.012
47. Agarwal D, Nowak C, Zhang NR, Pusztai L, Hatzis C. Functional germline variants as potential Co-oncogenes. NPJ Breast Cancer (2017) 3:1–4. doi: 10.1038/s41523-017-0051-5
48. Mamidi TKK, Wu J, Hicks C. Integrating germline and somatic variation information using genomic data for the discovery of biomarkers in prostate cancer. BMC Cancer (2019) 19:229. doi: 10.1186/s12885-019-5440-8
49. Qing T, Mohsen H, Marczyk M, Ye Y, O’Meara T, Zhao H, et al. Germline variant burden in cancer genes correlates with age at diagnosis and somatic mutation burden. Nat Commun (2020) 11:2438. doi: 10.1038/s41467-020-16293-7
50. Derpoorter C, Bordon V, Laureys G, Haerynck F, Lammens T. Genes at the crossroad of primary immunodeficiencies and cancer. Front Immunol (2018) 9:2544. doi: 10.3389/fimmu.2018.02544
51. Martin A, Li Z, Lin DP, Bardwell PD, Iglesias-Ussel MD, Edelmann W, et al. Msh2 ATPase activity is essential for somatic hypermutation at a-T basepairs and for efficient class switch recombination. J Exp Med (2003) 198:1171–8. doi: 10.1084/jem.20030880
52. Rada C, Ehrenstein MR, Neuberger MS, Milstein C. Hot spot focusing of somatic hypermutation in MSH2-deficient mice suggests two stages of mutational targeting. Immunity (1998) 9:135–41. doi: 10.1016/s1074-7613(00)80595-6
53. Offer SM, Pan-Hammarström Q, Hammarström L, Harris RS. Unique DNA repair gene variations and potential associations with the primary antibody deficiency syndromes IgAD and CVID. PLoS One (2010) 5:e12260. doi: 10.1371/journal.pone.0012260
54. Maurya S, Yang W, Tamai M, Zhang Q, Erdmann-Gilmore P, Bystry A, et al. Loss of KMT2C reprograms the epigenomic landscape in hPSCs resulting in NODAL overexpression and a failure of hemogenic endothelium specification. Epigenetics (2021) 17(2):220-38. doi: 10.1080/15592294.2021.1954780
55. Jiang P, Gu S, Pan D, Fu J, Sahu A, Hu X, et al. Signatures of T cell dysfunction and exclusion predict cancer immunotherapy response. Nat Med (2018) 24:1550–8. doi: 10.1038/s41591-018-0136-1
56. Lohr JG, Stojanov P, Lawrence MS, Auclair D, Chapuy B, Sougnez C, et al. Discovery and prioritization of somatic mutations in diffuse Large b-cell lymphoma (DLBCL) by whole-exome sequencing. Proc Natl Acad Sci U.S.A. (2012) 109:3879–84. doi: 10.1073/pnas.1121343109
57. Hoadley KA, Yau C, Hinoue T, Wolf DM, Lazar AJ, Drill E, et al. Cell-Of-Origin patterns dominate the molecular classification of 10,000 tumors from 33 types of cancer. Cell (2018) 173:291–304.e6. doi: 10.1016/j.cell.2018.03.022
58. Morin RD, Mungall K, Pleasance E, Mungall AJ, Goya R, Huff RD, et al. Mutational and structural analysis of diffuse Large b-cell lymphoma using whole-genome sequencing. Blood (2013) 122:1256–65. doi: 10.1182/blood-2013-02-483727
59. Morin RD, Mendez-Lago M, Mungall AJ, Goya R, Mungall KL, Corbett RD, et al. Frequent mutation of histone-modifying genes in non-Hodgkin lymphoma. Nature (2011) 476:298–303. doi: 10.1038/nature10351
60. Chapuy B, Stewart C, Dunford AJ, Kim J, Kamburov A, Redd RA, et al. Molecular subtypes of diffuse Large b cell lymphoma are associated with distinct pathogenic mechanisms and outcomes. Nat Med (2018) 24:679–90. doi: 10.1038/s41591-018-0016-8
61. Bogaert DJA, Dullaers M, Lambrecht BN, Vermaelen KY, De Baere E, Haerynck F. Genes associated with common variable immunodeficiency: One diagnosis to rule them all? J Med Genet (2016) 53:575–90. doi: 10.1136/jmedgenet-2015-103690
62. Odnoletkova I, Kindle G, Quinti I, Grimbacher B, Knerr V, Gathmann B, et al. The burden of common variable immunodeficiency disorders: A retrospective analysis of the European society for immunodeficiency (ESID) registry data. Orphanet J Rare Dis (2018) 13:201. doi: 10.1186/s13023-018-0941-0
63. Ramirez NJ, Posadas-Cantera S, Caballero-Oteyza A, Camacho-Ordonez N, Grimbacher B. There is no gene for CVID - novel monogenetic causes for primary antibody deficiency. Curr Opin Immunol (2021) 72:176–85. doi: 10.1016/j.coi.2021.05.010
64. Merdan S, Subramanian K, Ayer T, Van Weyenbergh J, Chang A, Koff JL, et al. Gene expression profiling-based risk prediction and profiles of immune infiltration in diffuse Large b-cell lymphoma. Blood Cancer J (2021) 11:1–12. doi: 10.1038/s41408-020-00404-0
65. Carter H, Marty R, Hofree M, Gross AM, Jensen J, Fisch KM, et al. Interaction landscape of inherited polymorphisms with somatic events in cancer. Cancer Discovery (2017) 7:410–23. doi: 10.1158/2159-8290.CD-16-1045
66. Riaz IB, Faridi W, Patnaik MM, Abraham RS. A systematic review on predisposition to lymphoid (B and T cell) neoplasias in patients with primary immunodeficiencies and immune dysregulatory disorders (Inborn errors of immunity). Front Immunol (2019) 10:777. doi: 10.3389/fimmu.2019.00777
67. Echeverri C, Fisher S, King D, Craig FE. Immunophenotypic variability of b-cell non-Hodgkin lymphoma: A retrospective study of cases analyzed by flow cytometry. Am J Clin Pathol (2002) 117:615–20. doi: 10.1309/AAYH-1FK8-38PL-Q6DT
68. Kuang Z, Guo L, Li X. Identification of key genes and pathways associated with classical Hodgkin lymphoma by bioinformatics analysis. Mol Med Rep (2017) 16:4685–93. doi: 10.3892/mmr.2017.7158
69. Sanchez-Vega F, Mina M, Armenia J, Chatila WK, Luna A, La KC, et al. Oncogenic signaling pathways in the cancer genome atlas. Cell (2018) 173:321–37.e10. doi: 10.1016/j.cell.2018.03.035
70. Bassing CH, Swat W, Alt FW. The mechanism and regulation of chromosomal V(D)J recombination. Cell (2002) 109 Suppl:S45–55. doi: 10.1016/s0092-8674(02)00675-x
71. Liao MJ, Van Dyke T. Critical role for atm in suppressing V(D)J recombination-driven thymic lymphoma. Genes Dev (1999) 13:1246–50. doi: 10.1101/gad.13.10.1246
72. Care MA, Westhead DR, Tooze RM. Gene expression meta-analysis reveals immune response convergence on the ifnγ-STAT1-IRF1 axis and adaptive immune resistance mechanisms in lymphoma. Genome Med (2015) 7:96. doi: 10.1186/s13073-015-0218-3
73. Sheng L, Fu D, Cao Y, Huo Y, Wang S, Shen R, et al. Integrated genomic and transcriptomic analyses of diffuse Large b-cell lymphoma with multiple abnormal immunologic markers. Front Oncol (2022) 12:790720. doi: 10.3389/fonc.2022.790720
74. Morales-Martínez M, Vega MI. Roles and regulation of BCL-xL in hematological malignancies. Int J Mol Sci (2022) 23:2193. doi: 10.3390/ijms23042193
75. Iqbal J, Meyer PN, Smith LM, Johnson NA, Vose JM, Greiner TC, et al. BCL2 predicts survival in germinal center B-Cell–like diffuse Large b-cell lymphoma treated with CHOP-like therapy and rituximab. Clin Cancer Res Off J Am Assoc Cancer Res (2011) 17:7785–95. doi: 10.1158/1078-0432.CCR-11-0267
76. Tessoulin B, Papin A, Gomez-Bougie P, Bellanger C, Amiot M, Pellat-Deceunynck C, et al. BCL2-family dysregulation in b-cell malignancies: From gene expression regulation to a targeted therapy biomarker. Front Oncol (2018) 8:645. doi: 10.3389/fonc.2018.00645
77. Del Pino Molina L, Torres Canizales JM, Pernía O, Rodríguez Pena R, Ibanez de Caceres I, López Granados E. Defective bcl-2 expression in memory b cells from common variable immunodeficiency patients. Clin Exp Immunol (2021) 203:341–50. doi: 10.1111/cei.13522
78. Lee HH, Dadgostar H, Cheng Q, Shu J, Cheng G. NF-kappaB-Mediated up-regulation of bcl-X and bfl-1/A1 is required for CD40 survival signaling in b lymphocytes. Proc Natl Acad Sci U.S.A. (1999) 96:9136–41. doi: 10.1073/pnas.96.16.9136
79. Chiu A, Xu W, He B, Dillon SR, Gross JA, Sievers E, et al. Hodgkin Lymphoma cells express TACI and BCMA receptors and generate survival and proliferation signals in response to BAFF and APRIL. Blood (2007) 109:729–39. doi: 10.1182/blood-2006-04-015958
80. Briones J, Timmerman JM, Hilbert DM, Levy R. BLyS and BLyS receptor expression in non-hodgkin’s lymphoma. Exp Hematol (2002) 30:135–41. doi: 10.1016/s0301-472x(01)00774-3
81. Schwaller J, Schneider P, Mhawech-Fauceglia P, McKee T, Myit S, Matthes T, et al. Neutrophil-derived APRIL concentrated in tumor lesions by proteoglycans correlates with human b-cell lymphoma aggressiveness. Blood (2007) 109:331–8. doi: 10.1182/blood-2006-02-001800
82. Moreaux J, Veyrune J-L, De Vos J, Klein B. APRIL is overexpressed in cancer: Link with tumor progression. BMC Cancer (2009) 9:83. doi: 10.1186/1471-2407-9-83
83. Salzer U, Bacchelli C, Buckridge S, Pan-Hammarström Q, Jennings S, Lougaris V, et al. Relevance of biallelic versus monoallelic TNFRSF13B mutations in distinguishing disease-causing from risk-increasing TNFRSF13B variants in antibody deficiency syndromes. Blood (2009) 113:1967–76. doi: 10.1182/blood-2008-02-141937
84. Pan-Hammarström Q, Salzer U, Du L, Björkander J, Cunningham-Rundles C, Nelson DL, et al. Reexamining the role of TACI coding variants in common variable immunodeficiency and selective IgA deficiency. Nat Genet (2007) 39:429–30. doi: 10.1038/ng0407-429
85. Poodt AEJ, Driessen GJA, de Klein A, van Dongen JJM, van der Burg M, de Vries E. TACI mutations and disease susceptibility in patients with common variable immunodeficiency. Clin Exp Immunol (2009) 156:35–9. doi: 10.1111/j.1365-2249.2008.03863.x
86. Freiberger T, Ravčuková B, Grodecká L, Pikulová Z, Štikarovská D, Pešák S, et al. Sequence variants of the TNFRSF13B gene in Czech CVID and IgAD patients in the context of other populations. Hum Immunol (2012) 73:1147–54. doi: 10.1016/j.humimm.2012.07.342
87. . Available at: https://gnomad.broadinstitute.org/gene/ENSG00000240505?dataset=gnomad_r2_1 (Accessed June 13, 2022). TNFRSF13B | gnomAD v2.1.1 | gnomAD.
88. Sever R, Brugge JS. Signal transduction in cancer. Cold Spring Harb Perspect Med (2015) 5:a006098. doi: 10.1101/cshperspect.a006098
89. Chen S, Benbarche S, Abdel-Wahab O. Splicing factor mutations in hematologic malignancies. Blood (2021) 138:599–612. doi: 10.1182/blood.2019004260
90. Lu Y, Yang X. The pivotal roles of neddylation pathway in immunoregulation. Immun Inflammation Dis (2020) 8:782–92. doi: 10.1002/iid3.335
91. Vijg J, Dong X. Pathogenic mechanisms of somatic mutation and genome mosaicism in aging. Cell (2020) 182:12–23. doi: 10.1016/j.cell.2020.06.024
92. Garcia-Gomez A, Rodríguez-Ubreva J, Ballestar E. Epigenetic interplay between immune, stromal and cancer cells in the tumor microenvironment. Clin Immunol Orlando Fla (2018) 196:64–71. doi: 10.1016/j.clim.2018.02.013
93. Rodríguez-Ubreva J, Arutyunyan A, Bonder MJ, Del Pino-Molina L, Clark SJ, de la Calle-Fabregat C, et al. Single-cell atlas of common variable immunodeficiency shows germinal center-associated epigenetic dysregulation in b-cell responses. Nat Commun (2022) 13:1779. doi: 10.1038/s41467-022-29450-x
94. van Schouwenburg PA, Davenport EE, Kienzler A-K, Marwah I, Wright B, Lucas M, et al. Application of whole genome and RNA sequencing to investigate the genomic landscape of common variable immunodeficiency disorders. Clin Immunol Orlando Fla (2015) 160:301–14. doi: 10.1016/j.clim.2015.05.020
95. Thouenon R, Moreno-Corona N, Poggi L, Durandy A, Kracker S. Activated PI3Kinase delta syndrome-a multifaceted disease. Front Pediatr (2021) 9:652405. doi: 10.3389/fped.2021.652405
96. Wang S, Wang W, Zhang X, Gui J, Zhang J, Guo Y, et al. A somatic mutation in PIK3CD unravels a novel candidate gene for lymphatic malformation. Orphanet J Rare Dis (2021) 16:208. doi: 10.1186/s13023-021-01782-9
97. Bortolini Silveira A, Bidard F-C, Tanguy M-L, Girard E, Trédan O, Dubot C, et al. Multimodal liquid biopsy for early monitoring and outcome prediction of chemotherapy in metastatic breast cancer. NPJ Breast Cancer (2021) 7:115. doi: 10.1038/s41523-021-00319-4
98. Mendiratta G, Ke E, Aziz M, Liarakos D, Tong M, Stites EC. Cancer gene mutation frequencies for the U.S. population. Nat Commun (2021) 12:5961. doi: 10.1038/s41467-021-26213-y
99. Bodian DL, McCutcheon JN, Kothiyal P, Huddleston KC, Iyer RK, Vockley JG, et al. Germline variation in cancer-susceptibility genes in a healthy, ancestrally diverse cohort: Implications for individual genome sequencing. PLoS One (2014) 9:e94554. doi: 10.1371/journal.pone.0094554
100. Russell MA, Pigors M, Houssen ME, Manson A, Kelsell D, Longhurst H, et al. A novel De novo activating mutation in STAT3 identified in a patient with common variable immunodeficiency (CVID). Clin Immunol Orlando Fla (2018) 187:132–6. doi: 10.1016/j.clim.2017.11.007
101. Zhu F, Wang KB, Rui L. STAT3 activation and oncogenesis in lymphoma. Cancers (2020) 12:19. doi: 10.3390/cancers12010019
102. Vogel TP, Milner JD, Cooper MA. The ying and yang of STAT3 in human disease. J Clin Immunol (2015) 35:615–23. doi: 10.1007/s10875-015-0187-8
103. Aaronson DS, Horvath CM. A road map for those who don’t know JAK-STAT. Science (2002) 296:1653–5. doi: 10.1126/science.1071545
104. Adamson AS, Collins K, Laurence A, O’Shea JJ. The current STATus of lymphocyte signaling: New roles for old players (STATs in lymphocyte signaling). Curr Opin Immunol (2009) 21:161–6. doi: 10.1016/j.coi.2009.03.013
105. Gonzalez-Mancera MS, Johnson B, Mirsaeidi M. STAT3 gain-of-Function mutation in a patient with pulmonary mycobacterium abscessus infection. Respir Med Case Rep (2020) 30:101125. doi: 10.1016/j.rmcr.2020.101125
106. Faletti L, Ehl S, Heeg M. Germline STAT3 gain-of-Function mutations in primary immunodeficiency: Impact on the cellular and clinical phenotype. BioMed J (2021) 44:412–21. doi: 10.1016/j.bj.2021.03.003
107. Deau M-C, Heurtier L, Frange P, Suarez F, Bole-Feysot C, Nitschke P, et al. A human immunodeficiency caused by mutations in the PIK3R1 gene. J Clin Invest (2015) 125:1764–5. doi: 10.1172/JCI81746
108. Aggarwal V, Banday AZ, Jindal AK, Das J, Rawat A. Recent advances in elucidating the genetics of common variable immunodeficiency. Genes Dis (2020) 7:26–37. doi: 10.1016/j.gendis.2019.10.002
109. Maffucci P, Filion CA, Boisson B, Itan Y, Shang L, Casanova J-L, et al. Genetic diagnosis using whole exome sequencing in common variable immunodeficiency. Front Immunol (2016) 7:220. doi: 10.3389/fimmu.2016.00220
110. Browning MJ, Chandra A, Carbonaro V, Okkenhaug K, Barwell J. Cowden’s syndrome with immunodeficiency. J Med Genet (2015) 52:856–9. doi: 10.1136/jmedgenet-2015-103266
111. Taylor H, Laurence ADJ, Uhlig HH. The role of PTEN in innate and adaptive immunity. Cold Spring Harb Perspect Med (2019) 9:a036996. doi: 10.1101/cshperspect.a036996
112. PIK3R1 - my cancer genome. Available at: https://www.mycancergenome.org/content/gene/pik3r1/ (Accessed November 23, 2021).
113. Ochs HD. Common variable immunodeficiency (CVID): New genetic insight and unanswered questions. Clin Exp Immunol (2014) 178:5–6. doi: 10.1111/cei.12491
114. Chapel H. Common variable immunodeficiency disorders (CVID) — diagnoses of exclusion, especially combined immune defects. J Allergy Clin Immunol Pract (2016) 4:1158–9. doi: 10.1016/j.jaip.2016.09.006
115. Singh A, Joshi V, Jindal AK, Mathew B, Rawat A. An updated review on activated PI3 kinase delta syndrome (APDS). Genes Dis (2020) 7:67–74. doi: 10.1016/j.gendis.2019.09.015
116. Durandy A, Kracker S. Increased activation of PI3 kinase-δ predisposes to b-cell lymphoma. Blood (2020) 135:638–43. doi: 10.1182/blood.2019002072
117. Carpier J-M, Lucas CL. Epstein-Barr Virus susceptibility in activated Pi3kδ syndrome (APDS) immunodeficiency. Front Immunol (2017) 8:2005. doi: 10.3389/fimmu.2017.02005
118. Nunes-Santos CJ, Uzel G, Rosenzweig SD. PI3K pathway defects leading to immunodeficiency and immune dysregulation. J Allergy Clin Immunol (2019) 143:1676–87. doi: 10.1016/j.jaci.2019.03.017
119. Péricart S, Tosolini M, Gravelle P, Rossi C, Traverse-Glehen A, Amara N, et al. Profiling immune escape in hodgkin’s and diffuse Large b-cell lymphomas using the transcriptome and immunostaining. Cancers (2018) 10:E415. doi: 10.3390/cancers10110415
120. Tessarin G, Rossi S, Baronio M, Gazzurelli L, Colpani M, Benvenuto A, et al. Activated phosphoinositide 3-kinase delta syndrome 1: Clinical and immunological data from an Italian cohort of patients. J Clin Med (2020) 9:3335. doi: 10.3390/jcm9103335
121. Takeda AJ, Zhang Y, Dornan GL, Siempelkamp BD, Jenkins ML, Matthews HF, et al. Novel PIK3CD mutations affecting n-terminal residues of P110δ cause activated PI3Kδ syndrome (APDS) in humans. J Allergy Clin Immunol (2017) 140:1152–1156.e10. doi: 10.1016/j.jaci.2017.03.026
122. Cohen JI. Primary immunodeficiencies associated with EBV disease. Curr Top Microbiol Immunol (2015) 390:241–65. doi: 10.1007/978-3-319-22822-8_10
123. Singh G. Serum free light chain assay and κ/λ ratio: Performance in patients with monoclonal gammopathy-high false negative rate for κ/λ ratio. J Clin Med Res (2017) 9:46–57. doi: 10.14740/jocmr2802w
124. Guevara-Hoyer K, Ochoa-Grullón J, Fernández-Arquero M, Cárdenas M, Pérez de Diego R, Sánchez-Ramón S. Serum free immunoglobulins light chains: A common feature of common variable immunodeficiency? Front Immunol (2020) 11:2004. doi: 10.3389/fimmu.2020.02004
125. Levy Y, Gupta N, Le Deist F, Garcia C, Fischer A, Weill JC, et al. Defect in IgV gene somatic hypermutation in common variable immuno-deficiency syndrome. Proc Natl Acad Sci U.S.A. (1998) 95:13135–40. doi: 10.1073/pnas.95.22.13135
126. Cabañas Morafraile E, Pérez-Peña J, Fuentes-Antrás J, Manzano A, Pérez-Segura P, Pandiella A, et al. Genomic correlates of DNA damage in breast cancer subtypes. Cancers (2021) 13(9):2117. doi: 10.3390/cancers13092117
127. Mortaz E, Tabarsi P, Mansouri D, Khosravi A, Garssen J, Velayati A, et al. Cancers related to immunodeficiencies: Update and perspectives. Front Immunol (2016) 7:365. doi: 10.3389/fimmu.2016.00365
128. Torgovnick A, Schumacher B. DNA Repair mechanisms in cancer development and therapy. Front Genet (2015) 6:157. doi: 10.3389/fgene.2015.00157
129. Leeksma OC, de Miranda NF, Veelken H. Germline mutations predisposing to diffuse Large b-cell lymphoma. Blood Cancer J (2017) 7:e532. doi: 10.1038/bcj.2017.15
130. Huang K-L, Mashl RJ, Wu Y, Ritter DI, Wang J, Oh C, et al. Pathogenic germline variants in 10,389 adult cancers. Cell (2018) 173:355–370.e14. doi: 10.1016/j.cell.2018.03.039
131. Kindler O, Quehenberger F, Benesch M, Seidel MG. The iceberg map of germline mutations in childhood cancer: Focus on primary immunodeficiencies. Curr Opin Pediatr (2018) 30:855–63. doi: 10.1097/MOP.0000000000000680
132. Lincoln SE, Nussbaum RL, Kurian AW, Nielsen SM, Das K, Michalski S, et al. Yield and utility of germline testing following tumor sequencing in patients with cancer. JAMA Netw Open (2020) 3:e2019452. doi: 10.1001/jamanetworkopen.2020.19452
133. Liu G, Bolkov MA, Tuzankina IA, Danilova IG. Identification of candidate disease genes in patients with common variable immunodeficiency. Quant Biol (2019) 7:190–201. doi: 10.1007/s40484-019-0174-9
134. Liu M, Liu X, Suo P, Gong Y, Qu B, Peng X, et al. The contribution of hereditary cancer-related germline mutations to lung cancer susceptibility. Transl Lung Cancer Res (2020) 9:646–58. doi: 10.21037/tlcr-19-403
135. Liu Y, Zeng G. Cancer and innate immune system interactions: Translational potentials for cancer immunotherapy. J Immunother Hagerstown Md (2012) 35:299–308. doi: 10.1097/CJI.0b013e3182518e83
136. Attarbaschi A, Carraro E, Abla O, Barzilai-Birenboim S, Bomken S, Brugieres L, et al. Non-Hodgkin lymphoma and pre-existing conditions: Spectrum, clinical characteristics and outcome in 213 children and adolescents. Haematologica (2016) 101:1581–91. doi: 10.3324/haematol.2016.147116
137. Abolhassani H, Wang Y, Hammarström L, Pan-Hammarström Q. Hallmarks of cancers: Primary antibody deficiency versus other inborn errors of immunity. Front Immunol (2021) 12:720025. doi: 10.3389/fimmu.2021.720025
138. Ramroop JR, Gerber MM, Toland AE. Germline variants impact on somatic events during tumorigenesis. Trends Genet TIG (2019) 35:515–26. doi: 10.1016/j.tig.2019.04.005
139. Porter CC. Germ line mutations associated with leukemias. Hematol Am Soc Hematol Educ Program (2016) 2016:302–8. doi: 10.1182/asheducation-2016.1.302
140. Carsetti R, Di Sabatino A, Rosado MM, Cascioli S, Piano Mortari E, Milito C, et al. Lack of gut secretory immunoglobulin a in memory b-cell dysfunction-associated disorders: A possible gut-spleen axis. Front Immunol (2019) 10:2937. doi: 10.3389/fimmu.2019.02937
141. Chen L, Yang L, Yao L, Kuang X-Y, Zuo W-J, Li S, et al. Characterization of PIK3CA and PIK3R1 somatic mutations in Chinese breast cancer patients. Nat Commun (2018) 9:1357. doi: 10.1038/s41467-018-03867-9
142. Elgizouli M, Lowe DM, Speckmann C, Schubert D, Hülsdünker J, Eskandarian Z, et al. Activating Pi3kδ mutations in a cohort of 669 patients with primary immunodeficiency. Clin Exp Immunol (2016) 183:221–9. doi: 10.1111/cei.12706
143. AACR Project GENIE Consortium. AACR project GENIE: Powering precision medicine through an international consortium. Cancer Discovery (2017) 7:818–31. doi: 10.1158/2159-8290.CD-17-0151
144. The Rockefeller University. St. Giles laboratory of human genetics of infectious disease. Available at: https://hgc.rockefeller.edu/index.php (Accessed January 10, 2022).
Keywords: CVID, non-Hodgkin’s lymphoma, genomic, in silico, malignancy
Citation: Guevara-Hoyer K, Fuentes-Antrás J, de la Fuente-Muñoz E, Fernández-Arquero M, Solano F, Pérez-Segura P, Neves E, Ocaña A, Pérez de Diego R and Sánchez-Ramón S (2022) Genomic crossroads between non-Hodgkin’s lymphoma and common variable immunodeficiency. Front. Immunol. 13:937872. doi: 10.3389/fimmu.2022.937872
Received: 06 May 2022; Accepted: 07 July 2022;
Published: 05 August 2022.
Edited by:
Mikko Risto Juhana Seppänen, Helsinki University Central Hospital, FinlandReviewed by:
Saul Oswaldo Lugo Reyes, National Institute of Pediatrics, MexicoRik Schrijvers, KU Leuven, Belgium
Copyright © 2022 Guevara-Hoyer, Fuentes-Antrás, de la Fuente-Muñoz, Fernández-Arquero, Solano, Pérez-Segura, Neves, Ocaña, Pérez de Diego and Sánchez-Ramón. 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: Kissy Guevara-Hoyer, kissgh@gmail.com