Next Article in Journal
Relationships between Circulating Biomarkers and Body Composition Parameters in Patients with Metabolic Syndrome: A Community-Based Study
Next Article in Special Issue
Erythrocyte Glutathione S-Transferase Activity as a Sensitive Marker of Kidney Function Impairment in Children with IgA Vasculitis
Previous Article in Journal
Substrate Specificity Diversity of Human Terminal Deoxynucleotidyltransferase May Be a Naturally Programmed Feature Facilitating Its Biological Function
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

HLA Polymorphisms and Clinical Manifestations in IgA Vasculitis

1
Department of Pediatrics, University Hospital Centre Zagreb, University of Zagreb School of Medicine, 10000 Zagreb, Croatia
2
Tissue Typing Centre, Clinical Department for Transfusion Medicine and Transplantation Biology, University Hospital Centre Zagreb, 10000 Zagreb, Croatia
3
Department of Pediatrics, University Hospital Centre Osijek, Josip Juraj Strossmayer University of Osijek, Medical Faculty Osijek, 31000 Osijek, Croatia
4
Department of Pediatrics, University of Split School of Medicine, University Hospital Centre Split, 21000 Split, Croatia
5
Children’s Hospital Zagreb, Medical Faculty Osijek, Josip Juraj Strossmayer University of Osijek, 10000 Zagreb, Croatia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2024, 25(2), 882; https://doi.org/10.3390/ijms25020882
Submission received: 25 November 2023 / Revised: 22 December 2023 / Accepted: 9 January 2024 / Published: 10 January 2024
(This article belongs to the Special Issue Forward in Vasculitis: Genetics and Beyond)

Abstract

:
Studies concerning the genetic background of IgA vasculitis (IgAV), a small-vessel vasculitis occurring predominantly in childhood, have confirmed that the HLA-DRB1 gene showed a strong association with disease susceptibility. The objective was to investigate human leukocyte antigen (HLA) polymorphisms among Croatian patients with IgAV and their influence on disease susceptibility and clinical heterogeneity. Thus, 130 children with IgAV and 202 unrelated healthy individuals were enrolled in the study. Genomic DNA was extracted from whole peripheral blood, and HLA-A, -B, -DRB1 and -DQB1 gene polymorphism analysis was performed. HLA-A*03 (21.4% vs. 12.38%, p = 0.0092), HLA-B*37 (2.9% vs. 0.2%, p = 0.0054) and HLA-DRB1*12 (3.1% vs. 0.7%, p = 0.0216) alleles were significantly more frequent in IgAV patients than in controls. High-resolution typing revealed significantly higher frequency of HLA-DRB1*10:01 and -DRB1*11:03 among IgAV patients with gastrointestinal manifestations of the disease in comparison to controls (p = 0.0021 and p = 0.0301, respectively), while HLA-DRB1*14:01P occurred significantly more often in the group of patients who developed nephritis during the course of the disease (17.5% vs. 4.5%, p = 0.0006). Our results demonstrated that there is an association of HLA-A*03, HLA-B*37 and HLA-DRB1*12 alleles with susceptibility to IgAV in the examined Croatian pediatric population. Studies which aim to determine the HLA profile may contribute to the elucidation of the genetic background of autoimmune diseases, including IgAV.

1. Introduction

IgA vasculitis (IgAV), formerly known as Henoch-Schönlein purpura, is the most common systemic childhood vasculitis, with an estimated worldwide annual incidence of 3–55.9 cases per 100,000 children [1,2,3]. Using precise geostatistical methods, the average annual incidence of IgAV in the Republic of Croatia is estimated to be 6.79 per 100,000 children, which is similar to other European countries [2]. The main feature of IgAV is a palpable purpura or petechiae, predominantly affecting lower extremities accompanied by arthritis and/or arthralgias, abdominal pain and IgA vasculitis nephritis (IgAVN) [4]. IgAV is usually a self-limiting disease that lasts up to four weeks, with a prognosis that almost exclusively depends on the developed acute or chronic complications. Severe gastrointestinal and systemic manifestations refer to acute complications of the disease, while renal involvement is the most significant chronic complication and the main cause of morbidity and mortality among affected children [2,5,6]. According to recently published studies, children with severe cutaneous manifestations such as ulcerations and necrosis and purpuric rash that persists for more than a month, older age and with severe gastrointestinal symptoms (severe abdominal pain, intussusception, hematochezia and/or massive gastrointestinal bleeding) of IgAV are a particularly high-risk group for developing IgAVN [5,7].
A multifactorial ethiopathogenesis of IgAV has been proposed, in which various genetic and environmental factors trigger the disease. When comparing genetic variants between healthy individuals and IgAV patients in susceptibility to the disease, significant differences are found in human leukocyte antigen (HLA) genes. These genes encode molecules, which have a crucial role in presenting antigenic peptides to T cells and, therefore, facilitating the immune system to differentiate between self and non-self [8]. The HLA complex has been extensively studied for association with numerous diseases, and in the case of IgAV, previous studies have mostly pointed out the association with HLA class II genes, in particular the HLA-DRB1 locus [9,10,11,12,13].
In three genome-wide association studies (GWAS) so far, the first of them revealed the impact of the polymorphisms in the HLA-DQA1 and HLA-DQB1 intergenic zone and at the HLA-DRB1 locus (HLA-DRB1*11 and DRB1*13), and, more importantly, the second determined that the haplotype HLA-DRB1*01:01~DQB1*05:01~DQA1*01:01 was only associated with susceptibility to IgAV but not with other autoimmune diseases, highlighting that IgAV is a prototype of HLA class II disease. The third confirmed the strong association of HLA-DRB1 with susceptibility to IgAV [14,15,16]. The possible influence of HLA class I genes has also been investigated, reporting mixed results, probably due to the smaller sample size of the investigated cohort [15,17,18,19,20]. Apart from susceptibility to the disease itself, HLA class I gene polymorphisms may also affect clinical heterogeneity and disease severity, confirming the complexity of the genetic contribution in understanding the pathogenesis, course and outcome of IgAV.
For this purpose, in this study, we examined whether HLA-A, -B, -DRB1 and -DQB1 polymorphisms have an influence on the susceptibility and different phenotypes within the clinical course in a population of Croatian children with IgAV.

2. Results

The present study included 130 patients with IgAV and 202 control subjects for whom the HLA-A, -B, -DRB1 and -DQB1 polymorphism was analyzed. Information about the patients’ gender, age at disease onset and clinical features is shown in Table 1. Out of 130 children with IgAV, 71 were girls and 59 were boys, with a median age of 6.3 (4.4–8.1) years at the time of diagnosis. All patients had purpuric rash, 108 patients (83.1%) had joint involvement, 50 patients (38.5%) had gastrointestinal (GI) tract involvement, while 39 (30.0%) patients developed IgA vasculitis nephritis (IgAVN) during the course of the disease.
The results of the HLA-A, -B, -DRB1 and -DQB1 loci typing at the gene level (low-resolution level) in the entire patient group (N = 130) as well as among controls (N = 202) are presented in Table 2.
The HLA-A and -B typing results were available for 70 patients, HLA-DQB1 gene frequencies were determined for 91 patients, while the HLA-DRB1 locus was analyzed for all 130 IgAV patients.
Among sixteen different genes detected at the HLA-A locus, the most frequent one was HLA-A*02, both in the patient group and among controls (24.3% and 31.2%, respectively). The only significant difference in the frequency distribution of HLA-A genes between IgAV patients and healthy controls was observed for the HLA-A*03 gene, which was significantly more present among patients (21.4% vs. 12.38%, p = 0.0092, 95% CI—2.1% to 17.1%).
The HLA-B locus is characterized with a higher polymorphism level in comparison to other HLA class I loci, and this was proven by our results, with 28 different HLA-B genes detected in total in both study cohorts (25 different HLA-B genes among patients and 27 among controls). As opposed to the HLA-A locus, the gene frequencies at the HLA-B locus are more evenly distributed, with HLA-B*35 being the most frequently detected gene among IgAV patients, followed by HLA-B*44 and -B*51 genes (15.0%, 11.4% and 10.7%, respectively). In contrast, HLA-B*51, -B35 and -B*07 were the most frequent HLA-B genes in the control group (12.1%, 11.6% and 9.4%, respectively). This difference, however, was not significant. The only notable difference was detected for the presence of a less commonly found HLA-B gene in the Croatian population, the HLA-B*37 gene, which was present in 2.9% of IgAV patients but only among 0.2% of control subjects (p = 0.0054, 95% CI—0.51% to 6.9%).
The distribution of gene frequencies for the thirteen detected HLA-DRB1 genes also differed between IgAV patients and control subjects. The most frequent gene in the patient groups was the HLA-DRB1*11 gene (15.4%), as opposed to controls who carried the HLA-DRB1*13 gene in the majority of cases (16.3%). The only significant difference was calculated for the presence of the HLA-DRB1*12 gene, which occurred with a frequency of 3.1% in the patient group but only in 0.7% of controls (p = 0.0216, 95% CI—0.3% to 5.3%).
Finally, the HLA-DQB1 typing results analysis revealed seven different genes, with HLA-DQB1*05 being the most frequent HLA-DQB1 gene in both patient and control groups (31.9% and 25.7%, respectively). No significant difference in gene frequency for any of the observed HLA-DQB1 genes was found.
In order to determine the possible influence of the HLA-DRB1 polymorphism on different clinical manifestations of IgAV, the patient group was also divided into four subgroups: patients with GI involvement (IgAV and GI, N = 31), patients with nephritis (IgAVN, N = 20), patients with both GI involvement and nephritis (IgAVN and GI, N = 19) and patients without GI involvement and nephritis (IgAV, N = 60). The comparison of HLA-DRB1 allele frequencies found in the entire patient group and subgroups to those observed among controls, as well as the comparison of HLA-DRB1 allele frequencies distribution in different patient subgroups, are presented in Table 3 and Table 4. Only those HLA-DRB1 alleles for which a significant difference in frequency was observed in one of the comparisons are listed.
The analysis of high-resolution HLA typing results in addition to the division of patients into subgroups according to the presence of GI involvement and nephritis revealed several additional HLA-DRB1 alleles demonstrating association with IgAV. More precisely, HLA-DRB1*10:01 and -DRB1*11:03 were significantly more present among IgAV and GI patients in comparison to controls (p = 0.0021 and p = 0.0301, respectively), while HLA-DRB1*14:01P occurred significantly more often in the IgAVN group of patients than among controls (17.5% vs. 4.5%, p = 0.0006; 95% CI—3.9% to 27.6%). The increase in HLA-DRB1*10:01 allele frequency among IgAV and GI patients and -DRB1*14:01P allele frequency among IgAVN patients were also significant in comparison to those observed among patients without GI involvement and nephritis (p = 0.0050 and p = 0.0242, respectively). A comparison of HLA-DRB1 allele frequencies among IgAV and GI and IgAV patients revealed a significant increase in HLA-DRB1*04:02 allele frequency among patients with GI involvement (p = 0.0485). Finally, the alleles HLA-DRB1*04:04 and -DRB1*11:01 were detected significantly more often among patients with both GI involvement and nephritis in comparison to patients without those two clinical manifestations of the disease (p = 0.0117 and p = 0.0480, respectively). The results of multivariate logistic regression showed that HLA- DRB1*10 was significantly associated with GI symptoms and IgAV (OR = 6.67; 95% CI 1.62–27.49; p = 0.009), HLA-DRB1*14, with both IgAVN and IgAV (OR = 5.09; 95% CI 1.95–13.28; p = 0.001), and HLA-DRB1*12 with overall IgAV (OR = 5.84; 95% CI 1.37–24.88; p = 0.017) compared to controls.

3. Discussion

Although the pathogenesis of IgAV remains unclear, it is considered to be an immune complex-mediated disease, in which many genetic, immunological and environmental factors intertwine [21]. It seems that genetic predisposition is almost unquestionable in IgAV development since the disease shows a higher incidence among first-degree relatives of affected patients and sometimes familial aggregation [22,23,24]. However, we did not observe such an association here or for familial clustering, although a familial history of autoimmune rheumatic diseases was noticed among some patients.
We performed the present study to investigate the distribution of HLA class I and class II genes among Croatian children with IgAV and to examine whether these genes have an effect on the susceptibility, severity and clinical heterogeneity of IgAV in our population.
The results of HLA class I gene analysis in our cohort of patients revealed only two genes demonstrating a positive association with the disease, HLA-A*03 and HLA-B*37. The HLA-B*37 has not been identified as a risk factor for IgAV in any previous investigations, while for the HLA-A*03 gene, a positive association with IgAV patients with joint involvement was reported by Peru et al. [17]. This study also implicated HLA-B44 as a risk factor, but, although this gene occurred with an increased frequency among IgAV patients, in the present study, in comparison to controls, this difference did not reach statistical significance. The same observation is valid for the HLA-B*35 gene, which was found to be associated with IgAV among Turkish children [17] as well as among patients from Northwest Spain with renal implications [18]. This frequency of the HLA-B*35 gene was also increased in our IgAV patient group but without statistical significance. The disparity of our results regarding the HLA-B*35 gene could be explained by the small sample size, both in the present study since HLA class I typing was available for 90 IgAV patients but also in the abovementioned studies, where the sample size was 110 and 48 patients, respectively. Similarly, the HLA-B*40 (serological equivalent HLA-B60), described as a protective gene in a Mongolian study [25], was detected with a decreased frequency in our IgAV patient group but insufficiently for a significant association. It is also important to mention population differences between Mongolian and our population regarding the HLA-B frequencies. Previous studies also reported an increased susceptibility to IgAV among HLA-B*41:02-positive Spanish patients, irrespective of HLA-DRB1 status [19] and for HLA-A2- and HLA-A11-positive patients from Turkey [17], but our results did not confirm these findings. Finally, there were also reports, which observed no association of HLA class I antigens and IgAV [20].
As previously reported in GWAS studies, HLA class II genes have been implicated in susceptibility to IgAV, with special emphasis on HLA-DRB1 alleles [14,15,16]. Therefore, IgAV shares some features with giant cell arteritis (GCA) and anti-neutrophil cytoplasmic antibody (ANCA)-associated vasculitis (AAV), as these entities have also shown association with HLA class II genes [26,27,28]. In addition to GWAS studies, several studies evaluated the potential association between the HLA-DRB1 genes and IgAV in different ethnic groups. In the present study, a significant positive association with IgAV was detected for the HLA-DRB1*12 gene. This gene showed a possible susceptibility to autoimmune and rheumatic diseases, such as seronegative rheumatoid arthritis [29], undifferentiated spondyloarthritis [30], dermatomyositis and polymyositis [31,32], granulomatosis with polyangiitis [33] and sarcoidosis [34,35,36], with studies predominantly focused on patients from Asian populations. The association between IgAV and the HLA-DRB1*12 gene was previously reported in a large Chinese GWAS study where a significant independent signal of the HLA-DRB1*12:02 allele was demonstrated [16].
The division of our patient group into four subgroups regarding the presence/absence of GI or renal involvement revealed that this association is limited for patients without either GI involvement or nephritis. The difference in the strength of association shown for the HLA-DRB1*12 gene and HLA-DRB1*12:01 allele lies in the fact that the allelic variants of this gene in the Croatian population are not numerous, and HLA-DRB1*12:01 is found almost exclusively. However, there was an individual carrying the HLA-DRB1*12:02 allele in our patient group, which influenced the statistical calculation. The HLA-DRB1*12 gene has not been reported as a susceptibility factor in any of the previous studies, which, thus far, have described an increased susceptibility to IgAV for HLA-DRB1*01-positive Spanish, Italian and Iranian patients (Italian, Turkish and Iranian HLA-DRB1*11 allele carriers and Turkish patients positive for the HLA-DRB1*14 gene) [9,10,11,12,13]. On the other hand, HLA-DRB1*07 was negatively associated with the disease in the Spanish and Italian cohorts [9,10], while a decreased frequency of HLA-DRB1*10 and HLA-DRB1*03 genes was detected in the Turkish cohorts [13], so these genes seem to have a protective role in these ethnic groups. Both HLA-DRB1*03:01 and -DRB1*07:01 alleles occurred with a decreased frequency in our patient group when compared to the controls, but this decrease did not result in a significant p value. The decrease in the HLA-DRB1*03:01 allele’s frequency in our IgAV patient group is an interesting finding in general since this HLA-DRB1 allele has been identified as a risk factor in numerous autoimmune diseases thus far [37,38,39].
Opposingly, the HLA-DRB1*10 gene demonstrated a positive association with IgAV patients with GI involvement in our research. However, this finding has to be proven on a larger patient cohort and compared to data from other populations in order to determine whether this association is characteristic for Croatian patients with IgAV. HLA-DRB1 typing at high resolution revealed one more allele demonstrating a positive association with this clinical manifestation of IgAV, the HLA-DRB1*11:03 allele. Previous studies did not reveal association with the clinical features of IgAV regarding gastrointestinal or renal involvement, with any specific HLA-DRB1 gene [9,10,11,12,13], except for HLA-DRB1*11 with joint manifestations in Turkish children [13]. The same study also observed an increased frequency of the HLA-DRB1*13 gene in children with nephrotic proteinuria [13]. The exception among the conducted studies so far is an Indian one that did not reveal any association of HLA-DRB1 genes with susceptibility to IgAV, nor with clinical features of the disease [40]. Our results regarding the IgAV with renal manifestation indicated a positive association with HLA-DRB1*14:01P (HLA-DRB1*14:01/*14:54), which was implicated as a susceptibility marker in a study by Soylemezoglu et al. [13].
The analysis of our results regarding the HLA class II genes did not reveal any association with HLA-DQB1 genes, as opposed to several previous investigations, which reported the association of HLA-DQB1*05:01, -DQB1*03:01 and -DQB1*02:01 alleles with IgAV [9,15]. The association with HLA-DQB1 alleles, as the authors themselves have stated, might be a result of the linkage disequilibrium (LD) with HLA-DRB1 alleles, which show a primary association in the majority of the, thus far, investigated autoimmune diseases. The HLA-DQB1*05:01 and HLA*DQB1*03:01 are in LD with HLA-DRB1*01 and HLA-DRB1*11 genes, respectively, while HLA-DQB1*02 can be found in combination with HLA-DRB1*03 and -DRB1*07 genes. The positive association of HLA-DQB1*05:01 and HLA*DQB1*03:01 alleles, and negative association of the HLA-DQB1*02 gene (DQB1*02:01 and DQB1*02:02 alleles), exactly mirrors the direction of association, demonstrated by the corresponding HLA-DRB1 alleles, further corroborating the theory of the HLA-DRB1 locus being responsible for the association. Gluten enteropathy and multiple sclerosis are currently the only exceptions where the primary association with HLA complex is due to HLA-DQ specificities [41].
HLA genes, so far, are not useful diagnostic markers in IgAV [42]. However, the identification of genetic polymorphisms at HLA loci may be useful in determining individuals at higher risk of severe forms of the disease and poor prognosis, which particularly applies to patients with GI involvement and IgAVN, thus providing important information and improving our understanding of pathogenesis and clinical heterogeneity of IgAV.

Limitations

The main limitation of the present study is the small sample size and, for more accurate results, further studies should be performed on a larger Croatian population of IgAV patients. One potential problem that arises here is the fragmentation of care for patients with IgAV between rheumatologists and nephrologists; therefore, it is important that similar research is based on future collaborations of both specialties. In addition, some of the limitations that can generally accompany HLA studies in IgAV refer to variants in various non-HLA genes associated with immune and inflammatory response, such as genes for cytokines, chemokines, adhesion molecules, T lymphocytes, neoangiogenesis, the renin–angiotensin system and metabolism, that may also contribute to the etiopathogenesis of IgAV and could be associated with different disease phenotypes and disease severity [21,42]. Also, genetic differences within the HLA system are likely influenced by the environment, demographics, racial and regional factors.

4. Materials and Methods

4.1. Patient Selection and Study Design

A total of 130 children fulfilling the EULAR/PRINTO/PRES classification criteria for IgAV in the period between 2019 and 2022 from three Croatian pediatric rheumatology centers were enrolled in the study [4]. Samples of 3 mL peripheral blood containing ethylenediamintetraacetic acid (EDTA) were taken from each patient with IgAV. Clinical data and laboratory parameters were collected from patient medical records at the time of disease diagnosis until the last follow-up visit. The control group of 202 unrelated healthy individuals from different regions of Croatia was selected to compare the gene and allele frequencies with those obtained from IgAV patients. These subjects were typed for HLA-A, -B, -DRB1 and -DQB1 genes in previous studies [43,44]. Patient(s) and/or their parents or legal guardians of minors signed an informed consent, and study was conducted with the principles of the Declaration of Helsinki.

4.2. DNA Isolation

Genomic DNA from whole blood containing EDTA was extracted using a commercial kit (MagNA Pure LC DNA; Roche Diagnostics GmbH, Manheim, Germany) and following the manufacturer’s instructions. All DNA samples were stored at −20 °C.

4.3. HLA Class I and Class II Typing

HLA-A, -B and -DQB1 low-resolution typing was performed by the Polymerase Chain Reaction–Sequence-Specific Oligo probes (PCR-SSO) method (Immucor Inc., Stanford, CT, USA) and using the Luminex technology (Luminex Corporation, Austin, TX, USA). The software used for the generation of the HLA data was MatchIt DNA, version 2.5 (Immucor Inc., Stanford, CT, USA). The HLA-DRB1 typing was performed at the high-resolution level using the next-generation sequencing (NGS) method (GenDx, Utrecht, The Netherlands) and using the iSeq 100 Sequencing System (Illumina Inc., San Diego, CA, USA). The obtained results were analyzed using the NGSengine software 2.27.2 (GenDx, Utrecht, The Netherlands).

4.4. Statistical Analysis

Continuous data were described as mean and standard deviation (mean ± SD) and categorical variables as percentages. The differences between HLA-A, -B, -DRB1 and -DQB1 gene frequencies, as well as HLA-DRB1 allele frequencies in IgAV patient group, subgroups and controls were examined using the MedCalc software proportions calculator, version 22.014 (MedCalc Software Ltd, Ostend, Belgium) [45], which uses the “N-1” Chi-squared test, as recommended by Campbell [46] and Richardson [47]. The confidence interval is calculated according to the recommended method given by Altman et al. [48]. Additionally, multiple logistic regression analysis was performed to assess the association between the most common HLD-DRB1 genes in several phenotypes of IgAV and the control group. Statistical significance was defined as p < 0.05.

5. Conclusions

HLA genes are the main genetic factor associated with IgAV pathogenesis, with HLA class II genes being the most strongly associated. Our results also confirm this finding and, additionally, we showed that HLA genes were involved in gastrointestinal manifestations and nephritis, particular disease phenotypes that are highlighted as the most important in IgAV. Therefore, we can conclude that studies, which aim to determine the HLA profile of IgAV patients in a specific population, should be performed in order to define the precise role of HLA alleles but also HLA haplotypes in etiopathogenesis of the disease.

Author Contributions

Conceptualization, M.J.; methodology, K.S.J. and Z.G.; software, K.S.J. and M.S. (Matej Sapina); validation, M.H., M.S. (Mario Sestan), K.S.J., Z.G. and M.J.; formal analysis, M.S. (Matej Sapina); investigation, M.H., M.S. (Mario Sestan), N.K., S.S., M.F. and A.G.; resources, M.J.; data curation, M.H., M.S. (Mario Sestan), N.K., S.S., M.F., A.G. and M.J.; writing—original draft preparation, M.H. and K.S.J.; writing—review and editing, M.H., K.S.J., M.S. (Mario Sestan), N.K., S.S., M.F., A.G., Z.G. and M.J.; visualization, M.H. and K.S.J.; supervision, M.J.; project administration, M.J.; funding acquisition, M.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Croatian Science Foundation under the project IP-2019-04-8822.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Ethics Committee of University of Zagreb School of Medicine, Zagreb, Croatia (date: 18 September 2019; Protocol Number—Class: 641-01/19-02/01).

Informed Consent Statement

Patient(s) or their guardian(s)/legally authorized representative(s) provided written informed consent for the publication of patient information in the present study.

Data Availability Statement

All data are available in the manuscript. The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Acknowledgments

We thank all patients, their parents and all the medical staff involved in this study. Thanks to all healthy individuals who participated as a control group.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Piram, M.; Maldini, C.; Biscardi, S.; De Suremain, N.; Orzechowski, C.; Georget, E.; Regnard, D.; Koné-Paut, I.; Mahr, A. Incidence of IgA vasculitis in children estimated by four-source capture-recapture analysis: A population-based study. Rheumatology 2017, 56, 1358–1366. [Google Scholar] [CrossRef]
  2. Sapina, M.; Frkovic, M.; Sestan, M.; Srsen, S.; Ovuka, A.; Varga, M.B.; Kramaric, K.; Brdaric, D.; Milas, K.; Gagro, A.; et al. Geospatial clustering of childhood IgA vasculitis and IgA vasculitis associated nephritis. Ann. Rheum. Dis. 2021, 80, 610–616. [Google Scholar] [CrossRef] [PubMed]
  3. Jelusic, M.; Sestan, M.; Giani, T.; Cimaz, R. New insights and challenges associated with IgA vasculitis and IgA vasculitis with nephritis-is it time to change the paradigm of the most common systemic vasculitis in childhood? Front. Pediatr. 2022, 10, 853724. [Google Scholar] [CrossRef]
  4. Ozen, S.; Pistorio, A.; Iusan, S.M.; Bakkaloglu, A.; Herlin, T.; Brik, R.; Buoncompagni, A.; Lazar, C.; Bilge, I.; Uziel, Y.; et al. EULAR/PRINTO/PRES criteria for Henoch-Schonlein purpura, childhood polyarteritis nodosa, childhood Wegener granulomatosis and childhood Takayasu arteritis: Ankara 2008. Part II: Final classification criteria. Ann. Rheum. Dis. 2010, 69, 798–806. [Google Scholar] [CrossRef] [PubMed]
  5. Sestan, M.; Kifer, N.; Frkovic, M.; Sapina, M.; Srsen, S.; Varga, M.B.; Ovuka, A.; Held, M.; Gracanin, A.G.; Kozmar, A.; et al. Gastrointestinal involvement and its association with the risk for nephritis in IgA vasculitis. Ther. Adv. Musculoskelet. Dis. 2021, 13, 1759720X211024828. [Google Scholar] [CrossRef]
  6. Kifer, N.; Bulimbasic, S.; Sestan, M.; Held, M.; Kifer, D.; Srsen, S.; Gracanin, A.G.; Heshin-Bekenstein, M.; Giani, T.; Cimaz, R.; et al. Semiquantitative classification (SQC) and Oxford classifications predict poor renal outcome better than The International Study of Kidney Disease in Children (ISKDC) and Haas in patients with IgAV nephritis: A multicenter study. J. Nephrol. 2023, 36, 441–449. [Google Scholar] [CrossRef] [PubMed]
  7. Sestan, M.; Srsen, S.; Kifer, N.; Sapina, M.; Varga, M.B.; Ovuka, A.; Held, M.; Kozmar, A.; Frkovic, M.; Laskarin, G.; et al. Persistence and severity of cutaneous manifestations in IgA vasculitis is associated with development of IgA vasculitis nephritis in children. Dermatology 2022, 238, 340–346. [Google Scholar] [CrossRef]
  8. Dendrou, C.A.; Petersen, J.; Rossjohn, J.; Fugger, L. HLA variation and disease. Nat. Rev. Immunol. 2018, 18, 325–339. [Google Scholar] [CrossRef]
  9. Amoroso, A.; Berrino, M.; Canale, L.; Coppo, R.; Cornaglia, M.; Guarrera, S.; Mazzola, G.; Scolari, F. Immunogenetics of Henoch-Schönlein disease. Eur. J. Immunogenet. 1997, 24, 323–333. [Google Scholar] [CrossRef]
  10. Amoli, M.M.; Thomson, W.; Hajeer, A.H.; Calviño, M.C.; Garcia-Porrua, C.; Ollier, W.E.; Gonzalez-Gay, M.A. HLA-DRB1*01 association with Henoch-Schönlein purpura in patients from Northwest Spain. J. Rheumatol. 2001, 28, 1266–1270. [Google Scholar]
  11. López Mejías, R.; Genre, F.; Pérez, B.S.; Castañeda, S.; Ortego Centeno, N.; Llorca Díaz, F.J.; Ubilla García, B.; Remuzgo Martínez, S.; Mijares Díaz, V.; Pina Murcia, T.; et al. HLA-DRB1 association with Henoch-Schonlein purpura. Arthritis Rheumatol. 2015, 67, 823–827. [Google Scholar] [CrossRef]
  12. Rashidi, S.; Shiari, R.; Farivar, S. HLA-DRB1 gene polymorphisms in Iranian children with Henoch-Schönlein purpura. J. Res. Med. Sci. 2018, 23, 42. [Google Scholar]
  13. Soylemezoglu, O.; Peru, H.; Gonen, S.; Cetinyurek, A.; Buyan, N. HLA-DRB1 alleles and Henoch-Schonlein purpura: Susceptibility and severity of disease. J. Rheumatol. 2008, 35, 1165–1168. [Google Scholar]
  14. López-Mejias, R.; Carmona, F.D.; Castaneda, S.; Genre, F.; Remuzgo-Martínez, S.; Sevilla-Perez, B.; Ortego-Centeno, N.; Llorca, J.; Ubilla, B.; Mijares, V.; et al. A genome-wide association study suggests the HLA class II region as the major susceptibility locus for IgA vasculitis. Sci. Rep. 2017, 7, 5088. [Google Scholar] [CrossRef] [PubMed]
  15. Koskela, M.; Nihtilä, J.; Ylinen, E.; Kolho, K.-L.; Nuutinen, M.; Ritari, J.; Jahnukainen, T. HLA-DQ and HLA-DRB1 alleles associated with Henoch-Schönlein purpura nephritis in Finnish pediatric population: A genome-wide association study. Pediatr. Nephrol. 2021, 36, 2311–2318. [Google Scholar] [CrossRef] [PubMed]
  16. Xia, L.; Chen, M.; Zhang, H.; Zheng, X.; Bao, J.; Gao, J.; Zhu, C.; Sun, L.; Xia, H.; Zhang, X. Genome-wide association study of 7661 Chinese Han individuals and fine-mapping major histocompatibility complex identifies HLA-DRB1 as associated with IgA vasculitis. J. Clin. Lab. Anal. 2022, 36, e24457. [Google Scholar] [CrossRef]
  17. Peru, H.; Soylemezoglu, O.; Gonen, S.; Cetinyurek, A.; Bakkaloğlu, S.A.; Buyan, N.; Hasanoglu, E. HLA class 1 associations in Henoch Schonlein purpura: Increased and decreased frequencies. Clin. Rheumatol. 2008, 27, 5–10. [Google Scholar] [CrossRef] [PubMed]
  18. Amoli, M.M.; Thomson, W.; Hajeer, A.H.; Calviño, M.C.; Garcia-Porrua, C.; Ollier, W.E.; Gonzalez-Gay, M.A. HLA-B35 association with nephritis in Henoch-Schönlein purpura. J. Rheumatol. 2002, 29, 948–949. [Google Scholar]
  19. López-Mejías, R.; Genre, F.; Pérez, B.S.; Castañeda, S.; Ortego-Centeno, N.; Llorca, J.; Ubilla, B.; Remuzgo-Martínez, S.; Mijares, V.; Pina, T.; et al. Association of HLA-B*41:02 with Henoch-Schönlein purpura (IgA vasculitis) in Spanish individuals irrespective of the HLA-DRB1 status. Arthritis Res. Ther. 2015, 17, 102. [Google Scholar] [CrossRef]
  20. Ostergaard, J.R.; Storm, K.; Lamm, L.U. Lack of association between HLA and Schoenlein-Henoch purpura. Tissue Antigens 1990, 35, 234–235. [Google Scholar] [CrossRef]
  21. Jelusic, M.; Sestan, M. IgA vasculitis or Henoch-Schönlein purpura: Genetics and beyond. Pediatr. Nephrol. 2021, 36, 2149–2153. [Google Scholar] [CrossRef]
  22. Ostini, A.; Simonetti, G.D.; Pellanda, G.; Bianchetti, M.G.; Ferrarini, A.; Milani, G.P. Familial Henoch-Schönlein syndrome. J. Clin. Rheumatol. 2016, 22, 80–81. [Google Scholar] [CrossRef]
  23. Zhang, Y.; Gu, W.; Mao, J. Sibling cases of Henoch-Schönlein purpura in two families and review of literature. Pediatr. Dermatol. 2008, 25, 393–395. [Google Scholar] [CrossRef]
  24. Motoyama, O.; Iitaka, K. Familial cases of Henoch-Schönlein purpura in eight families. Pediatr. Int. 2005, 47, 612–615. [Google Scholar] [CrossRef] [PubMed]
  25. Ren, S.M.; Yang, G.L.; Liu, C.Z.; Zhang, C.; Shou, Q.; Yu, S.; Li, W.; Su, X. Association between HLA-A and -B polymorphisms and susceptibility to Henoch-Schönlein purpura in Han and Mongolian children from Inner Mongolia. Genet. Mol. Res. 2012, 11, 221–228. [Google Scholar] [CrossRef]
  26. Carmona, F.D.; Mackie, S.L.; Martín, J.-E.; Taylor, J.C.; Vaglio, A.; Eyre, S.; Bossini-Castillo, L.; Castañeda, S.; Cid, M.C.; Hernández-Rodríguez, J.; et al. A large-scale genetic analysis reveals a strong contribution of the HLA class II region to giant cell arteritis susceptibility. Am. J. Hum. Genet. 2015, 96, 565–580. [Google Scholar] [CrossRef] [PubMed]
  27. Weyand, C.M.; Hicock, K.C.; Hunder, G.G.; Goronzy, J.J. The HLA-DRB1 locus as a genetic component in giant cell arteritis. Mapping of a disease-linked sequence motif to the antigen binding site of the HLA-DR molecule. J. Clin. Investig. 1992, 90, 2355–2361. [Google Scholar] [CrossRef]
  28. Spencer, S.J.; Burns, A.; Gaskin, G.; Pusey, C.D.; Rees, A.J. HLA class II specificities in vasculitis with antibodies to neutrophil cytoplasmic antigens. Kidney Int. 1992, 41, 1059–1063. [Google Scholar] [CrossRef]
  29. Terao, C.; Ohmura, K.; Kochi, Y.; Ikari, K.; Maruya, E.; Katayama, M.; Shimada, K.; Murasawa, A.; Honjo, S.; Takasugi, K.; et al. A large-scale association study identified multiple HLA-DRB1 alleles associated with ACPA-negative rheumatoid arthritis in Japanese subjects. Ann. Rheum. Dis. 2011, 70, 2134–2139. [Google Scholar] [CrossRef] [PubMed]
  30. Liao, H.T.; Lin, K.C.; Chen, C.H.; Liang, T.-H.; Lin, M.-W.; Tsai, C.-Y.; Yu, D.T.Y.; Chou, C.-T. Human leukocyte antigens in undifferentiated spondyloarthritis. Semin. Arthritis Rheum. 2007, 37, 198–201. [Google Scholar] [CrossRef]
  31. Kang, E.H.; Go, D.J.; Mimori, T.; Lee, S.J.; Kwon, H.M.; Park, J.W.; Park, M.H.; Song, E.Y.; Ha, Y.-J.; Lee, E.Y.; et al. Novel susceptibility alleles in HLA region for myositis and myositis specific autoantibodies in Korean patients. Semin. Arthritis Rheum. 2019, 49, 283–287. [Google Scholar] [CrossRef] [PubMed]
  32. Lin, J.M.; Zhang, Y.B.; Peng, Q.L.; Yang, H.; Shi, J.; Gu, M.; Zhao, W.; Wang, G. Genetic association of HLA-DRB1 multiple polymorphisms with dermatomyositis in Chinese population. HLA 2017, 90, 354–359. [Google Scholar] [CrossRef] [PubMed]
  33. Luo, H.; Chen, M.; Yang, R.; Xu, P.C.; Zhao, M.H. The association of HLA-DRB1 alleles with antineutrophil cytoplasmic antibody-associated systemic vasculitis in Chinese patients. Hum. Immunol. 2011, 72, 422–425. [Google Scholar] [CrossRef]
  34. Sikorová, K.; Moon, S.J.; Yoon, H.Y.; Strnad, A.; Song, J.W.; Petrek, M. HLA class II variants defined by next generation sequencing are associated with sarcoidosis in Korean patients. Sci. Rep. 2022, 12, 9302. [Google Scholar] [CrossRef]
  35. Ishihara, M.; Ohno, S.; Ishida, T.; Ando, H.; Naruse, T.; Nose, Y.; Inoko, H. Molecular genetic studies of HLA class II alleles in sarcoidosis. Tissue Antigens 1994, 43, 238–241. [Google Scholar] [CrossRef]
  36. Levin, A.M.; Adrianto, I.; Datta, I.; Iannuzzi, M.C.; Trudeau, S.; Li, J.; Drake, W.P.; Montgomery, C.G.; Rybicki, B.A. Association of HLA-DRB1 with sarcoidosis susceptibility and progression in African Americans. Am. J. Respir. Cell Mol. Biol. 2015, 53, 206–216. [Google Scholar] [CrossRef]
  37. Thomson, G.; Valdes, A.M.; Noble, J.A.; Kockum, I.; Grote, M.N.; Najman, J.; Erlich, H.A.; Cucca, F.; Pugliese, A.; Steenkiste, A.; et al. Relative predispositional effects of HLA class II DRB1-DQB1 haplotypes and genotypes on type 1 diabetes: A meta-analysis. Tissue Antigens 2007, 70, 110–127. [Google Scholar] [CrossRef] [PubMed]
  38. Jacobson, E.M.; Huber, A.; Tomer, Y. The HLA gene complex in thyroid autoimmunity: From epidemiology to etiology. J. Autoimmun. 2008, 30, 58–62. [Google Scholar] [CrossRef] [PubMed]
  39. Niu, Z.; Zhang, P.; Tong, Y. Value of HLA-DR genotype in systemic lupus erythematosus and lupus nephritis: A meta-analysis. Int. J. Rheum. Dis. 2015, 18, 17–28. [Google Scholar] [CrossRef]
  40. Aggarwal, R.; Gupta, A.; Naru, J.; Berka, N.; Nanda, N.; Suri, D.; Singh, S. HLA-DRB1 in Henoch-Schonlein purpura: A susceptibility study from north India. Hum. Immunol. 2016, 77, 555–558. [Google Scholar] [CrossRef]
  41. Liu, B.; Shao, Y.; Fu, R. Current research status of HLA in immune-related diseases. Immun. Inflamm. Dis. 2021, 9, 340–350. [Google Scholar] [CrossRef] [PubMed]
  42. López-Mejías, R.; Castañeda, S.; Genre, F.; Remuzgo-Martínez, S.; Carmona, F.D.; Llorca, J.; Blanco, R.; Martín, J.; González-Gay, M.A. Genetics of immunoglobulin-A vasculitis (Henoch-Schönlein purpura): An updated review. Autoimmun. Rev. 2018, 17, 301–315. [Google Scholar] [CrossRef]
  43. Grubic, Z.; Burek Kamenaric, M.; Mikulic, M.; Stingl Jankovic, K.; Maskalan, M.; Zunec, R. HLA-A, HLA-B and HLA-DRB1 allele and haplotype diversity among volunteer bone marrow donors from Croatia. Int. J. Immunogenet. 2014, 41, 211–221. [Google Scholar] [CrossRef]
  44. Grubic, Z.; Maskalan, M.; Svilicic, D.; Stingl Jankovic, K.; Zunec, R. Determination of HLA-A, -B, and -DRB1 allele and haplotype frequencies in the Croatian population based on a family study. Arch. Immunol. Ther. Exp. 2016, 64 (Suppl. S1), 83–88. [Google Scholar] [CrossRef] [PubMed]
  45. MedCalc Software Ltd. Comparison of Proportions Calculator. Version 22.014. Available online: https://www.medcalc.org/calc/comparison_of_proportions.php (accessed on 30 December 2023).
  46. Campbell, I. Chi-squared and Fisher-Irwin tests of two-by-two tables with small sample recommendations. Stat. Med. 2007, 26, 3661–3675. [Google Scholar] [CrossRef]
  47. Richardson, J.T.E. The analysis of 2 × 2 contingency tables—Yet again. Stat. Med. 2011, 30, 890. [Google Scholar] [CrossRef] [PubMed]
  48. Altman, D.G.; Machin, D.; Bryant, T.N.; Gardner, M.J. (Eds.) Statistics with confidence. In BMJ Books, 2nd ed.; BMJ Publishing Group: London, UK, 2000; p. 49. [Google Scholar]
Table 1. Demographic and clinical features of patients with IgAV.
Table 1. Demographic and clinical features of patients with IgAV.
Characteristics of the PatientsN = 130% of the Cohort
Demographics
  mean age (years)6.3 (4.4–8.1)
  female7154.6
  male5945.4
  ratio male/female1:1.2
Clinics
Cutaneous manifestations130100
  palpable purpura10379.2
  petechiae2015.4
  purpura and petechiae75.4
  rash extended above waist5743.8
  bullae, ulcerations, and necrotic lesions75.4
Joint involvement10883.1
  arthritis3224.6 *
  arthralgias3023.1 *
  arthritis and arthralgias4635.4 *
Gastrointestinal involvement5038.5
  abdominal pain3829.2 *
  positive FOBT2116.2 *
  intussusception32.3 *
Nephritis3930.0
  hematuria129.2 *
  proteinuria107.7 *
  hematuria and proteinuria1410.8 *
Orchitis1016.9 **
Legend: Data are presented as a whole number (%); FOBT—fecal occult blood test; *—expressed as a total number since a patient can have more than one clinical manifestation of affected organ system; **—applicable only for boys.
Table 2. The distribution of HLA-A, -B, -DRB1 and -DQB1 gene frequencies in IgAV patient and control groups.
Table 2. The distribution of HLA-A, -B, -DRB1 and -DQB1 gene frequencies in IgAV patient and control groups.
Gene Frequency (%)
HLA-A*IgAV Patients (N = 70)Controls (N = 202)p Value (95% CI)
0115.014.11/
0224.331.19/
0321.412.380.0092 (2.1–17.1%)
117.15.94/
231.41.73/
2410.713.12/
252.14.21/
265.03.22/
290.70.99/
301.40.74/
311.42.48/
323.64.21/
331.42.72/
660.70.25/
682.92.48/
690.70.25/
HLA-B*IgAV Patients (N = 70)Controls (N = 202)p Value (95% CI)
0710.09.4/
085.79.2/
135.03.5/
14 (B64)0.00.7/
14 (B65)0.73.5/
15 (B62)6.44.7/
15 (B63)0.71.0/
185.08.2/
273.63.5/
3515.011.6/
372.90.20.0054 (0.51–6.9%)
384.34.5/
391.42.7/
40 (B60)0.72.7/
40 (B61)3.65.0/
410.70.5/
4411.47.7/
470.70.0/
480.00.2/
492.12.0/
501.41.0/
5110.712.1/
521.41.2/
530.00.2/
552.91.5/
561.40.5/
571.41.5/
580.71.2/
HLA-DRB1*IgAV Patients (N = 130)Controls (N = 202)p Value (95% CI)
0113.510.4/
037.710.6/
049.29.9/
078.110.9/
084.62.5/
090.00.2/
101.91.0/
1115.414.9/
123.10.70.0216 (0.3%–5.3%)
1312.316.3/
147.34.5/
159.28.4/
167.79.7/
HLA-DQB1*IgAV Patients (N = 91)Controls (N = 202)p Value (95% CI)
0215.919.3/
03 (DQ7)20.918.1/
03 (DQ8)6.68.4/
03 (DQ9)1.12.5/
043.31.7/
0531.925.7/
0620.324.3/
Legend: in parentheses—serological equivalent at HLA-B locus; ‘/’—not significant p value.
Table 3. Distribution of HLA-DRB1 gene (A) and allele (B) frequencies in patient subgroups (IgAV and GI (N = 31), IgAVN (N = 20), IgAVN and GI (N = 19) and IgAV (N = 60)) and controls (N = 202).
Table 3. Distribution of HLA-DRB1 gene (A) and allele (B) frequencies in patient subgroups (IgAV and GI (N = 31), IgAVN (N = 20), IgAVN and GI (N = 19) and IgAV (N = 60)) and controls (N = 202).
Gene/Allele Frequency (%)
HLA-DRB1*IgAV and GI
(n = 31)
IgAVN
(n = 20)
IgAVN and GI (n = 19)IgAV
(n = 60)
Controls
(n = 202)
(A)
106.502.601.0
121.62.52.64.20.7
146.517.52.65.84.5
163.215.02.69.29.7
(B)
04:023.20002.7
04:0402.55.301.2
10:016.502.601.0
11:019.72.513.24.28.4
11:033.2001.70.5
12:011.62.52.63.30.7
14:01P6.517.52.65.84.5
Legend: IgAV and GI—IgA vasculitis with gastrointestinal tract involvement, IgAVN—IgA vasculitis nephritis, IgAVN and GI—IgA vasculitis nephritis with gastrointestinal tract involvement; IgAV—IgA vasculitis; HLA-DRB1*14:01P—HLA-DRB1*14:01/*14:54.
Table 4. The comparison of HLA-DRB1 gene (A) and allele (B) frequencies according to patient subgroups (IgAV and GI (N = 31), IgAVN (N = 20), IgAVN and GI (N = 19) and IgAV (N = 60)) and controls (N = 202).
Table 4. The comparison of HLA-DRB1 gene (A) and allele (B) frequencies according to patient subgroups (IgAV and GI (N = 31), IgAVN (N = 20), IgAVN and GI (N = 19) and IgAV (N = 60)) and controls (N = 202).
p Value (95% CI)
HLA-DRB1*IgAV and GI vs. ControlsIgAVN vs. ControlsIgAVN and GI vs. ControlsIgAV vs. Controls
(A)
100.0021 (1.3–14.5%)///
12///0.0073 (0.7–8.7%)
14/0.0006 (3.9–27.6%)//
16////
(B)
04:02////
04:04////
10:010.0021 (1.3–14.5%)///
11:01////
11:030.0301 (0.1–10.5%)///
12:01///0.0302 (0.1–7.5%)
14:01P/0.0006 (3.9–27.6%)//
p Value (95% CI)
HLA-DRB1*IgAV and GI
vs. IgAV
IgAVN
vs. IgAV
IgAVN and GI
vs. IgAV
IgAV and GI
vs. IgAVN
IgAV and GI
vs. IgAVN and GI
IgAVN vs.
IgAVN and GI
(A)
100.0050
(1.5–15.5%)
/////
12//////
14/0.0242
(1.2–26.4%)
///0.0311
(0.9–29.5%)
16///0.0311
(0.7–26.1%)
//
(B)
04:020.0485
(0.7–11.0%)
/////
04:04//0.0117
(0.4%–17.3%)
///
10:010.0050
(1.5–15.5%)
/////
11:01//0.0480
(0.1–23.4%)
///
11:03//////
12:01//////
14:01P/0.0242
(1.2–26.4%)
///0.0311
(0.9–29.5%)
Legend: ‘/’—not significant p value; HLA-DRB1*14:01P—HLA-DRB1*14:01/*14:54.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Held, M.; Stingl Jankovic, K.; Sestan, M.; Sapina, M.; Kifer, N.; Srsen, S.; Frkovic, M.; Gagro, A.; Grubic, Z.; Jelusic, M. HLA Polymorphisms and Clinical Manifestations in IgA Vasculitis. Int. J. Mol. Sci. 2024, 25, 882. https://doi.org/10.3390/ijms25020882

AMA Style

Held M, Stingl Jankovic K, Sestan M, Sapina M, Kifer N, Srsen S, Frkovic M, Gagro A, Grubic Z, Jelusic M. HLA Polymorphisms and Clinical Manifestations in IgA Vasculitis. International Journal of Molecular Sciences. 2024; 25(2):882. https://doi.org/10.3390/ijms25020882

Chicago/Turabian Style

Held, Martina, Katarina Stingl Jankovic, Mario Sestan, Matej Sapina, Nastasia Kifer, Sasa Srsen, Marijan Frkovic, Alenka Gagro, Zorana Grubic, and Marija Jelusic. 2024. "HLA Polymorphisms and Clinical Manifestations in IgA Vasculitis" International Journal of Molecular Sciences 25, no. 2: 882. https://doi.org/10.3390/ijms25020882

APA Style

Held, M., Stingl Jankovic, K., Sestan, M., Sapina, M., Kifer, N., Srsen, S., Frkovic, M., Gagro, A., Grubic, Z., & Jelusic, M. (2024). HLA Polymorphisms and Clinical Manifestations in IgA Vasculitis. International Journal of Molecular Sciences, 25(2), 882. https://doi.org/10.3390/ijms25020882

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop