1. Introduction
Systemic lupus erythematosus (SLE) is a chronic, multisystem, autoimmune, and inflammatory disease that is most often of unknown cause. It is characterized by the production of multiple autoantibodies, especially against nuclear components (e.g., double-stranded DNA, dsDNA); this can generate inflammatory-mediated effects in any organ and/or system. Childhood-onset forms (cSLE, onset before age 18) represent 20% of all cases [1,2]. Lupus nephritis (LN) is generally estimated to occur 10–30% more frequently in cSLE (40–70%) than in adult forms of SLE [2,3,4]. This complication is among those with the highest impacts on the life quality and survival of SLE patients since it increases the risk of developing end-stage kidney disease (ESKD) [5]. LN shows an earlier onset and more severe clinical course in cSLE compared to adult SLE forms; additionally, it has been reported as a negative predictor of survival [4].
Identifying biomarkers for the diagnosis, prognosis, and non-invasive evaluation of renal disease activity in SLE is a rapidly evolving field [6]. However, in the clinical setting, LN diagnosis and the evaluation of renal flares still rely on the assessment of traditional serum (e.g., serum creatinine, erythrocyte sedimentation rate, complement components C3 and C4, glomerular filtration rate, anti-C1q or anti-dsDNA autoantibody titers, etc.), and urinary (urine sediment, proteinuria, albumin–creatinine ratio, etc.) biomarkers, or on the histological evaluation of renal biopsies [6,7,8].
To date, genome-wide association studies (GWASs) have identified around 100 susceptibility genes for SLE development [6,9]. The relevant genetic changes have mainly been related to single-nucleotide variations in immune-response genes (e.g., HLA-DR3 in patients of European ancestry) [9,10], along with some gene copy number variations (CNVs) in non-immune-related genes, such as VANGL1 (1p13.1, MIM*610132) [11]. VANGL1 is an essential gene in the establishment of planar cell polarity (PCP), and heterozygous variants have been associated with neural tube defects (MIM#182940) and caudal regression syndrome (MIM#600145), but not with human kidney disease. A GWAS performed in 55 patients with SLE, 11 with Sjogren’s syndrome, and 11 healthy controls identified three SLE patients homozygous for a VANGL1 deletion-type CNV of ~3.17 kb (esv3587290) located at intron 7. Of these patients, two had LN. The authors furthermore evaluated the esv3587290 using the TaqMan® assay in 177 SLE patients of mainly European descent. The results revealed that the deletion was significantly associated with LN (χ2 = 27.06, 2 d.f., p < 0.0001) and demonstrated a gene–dosage effect. Moreover, the esv3587290 CNV seems to be a highly prevalent allele in the Australian Aboriginal Tiwi Islander population, who present high rates of kidney disease [11]. However, in a third SLE cohort of mainly Spanish-descent patients (N = 281, χ2 = 2.1, 1 d.f., p = 0.14), this association was not replicated [11]. Indeed, a Vangl1−/+ model mice showed spontaneous deposition of IgA and IgG, but not of IgM or complement, in the mesangium [11]. This led researchers to hypothesize that a deficiency of Vangl1 protein in heterozygous mice could alter the permeability of the glomerular endothelium for monomeric immunoglobulins. Despite the association of esv3587290 CNV with LN, no convincing evidence for a deleterious effect of this CNV on VANGL1 function was achieved in this study; however, the experimental murine models also supported that Vangl1 deficiency was only associated with the development of nephritis in the Vangl1−/+ mice injected with autoreactive serum, which further supports an altered glomerular endothelial permeability to autoreactive immunoglobulins. Whether this mechanism underlies an association of the esv3587290 CNV with LN development in humans remains unknown [11].
African American, Hispanic, and Asian SLE patients are at greater risk for developing and presenting more severe forms of LN compared to those from European-descent populations [10]. To our knowledge, no GWAS has yet identified any LN susceptibility locus linked to VANGL1 or the 1p13.1 region, even in populations of European ancestry [9,10,12,13]. Furthermore, Jiang et al. failed to replicate the association of LN with the esv3587290 CNV in a cohort of predominantly Spanish-descent patients [11]. Thus, additional work is needed to test whether this genetic variation could be a risk factor in other ethnicities or clinical presentations of SLE, such as the childhood-onset form. Finally, Jiang et al. employed whole-genome sequencing (WGS) to reveal that the esv3587290 deletion-type CNV varied in size among SLE patients [11], although the authors did not report precise nucleotide-resolution breakpoints using the Human Genome Variation Society (HGVS) nomenclature [14]. This size variability could suggest that distinct mutational events may generate the CNV. Meanwhile, given its high minor allelic frequencies in different populations (0.17–0.43) [11], we cannot discard the possibility of a one-time emergence of a single common allele.
Here, we sought to determine if the VANGL1 esv3587290 deletion-type CNV is associated with LN in a sample of Mexican cSLE patients, and to characterize the precise breakpoints of this rearrangement in our cSLE patients as well as a reference group of ethnically matched individuals.
2. Materials and Methods
2.1. Patient Selection
We enrolled a total of 66 unrelated Mexican children (11 males, 55 females; aged 7.8 to 18.6 years) who were born from Mexican parents in the central region of Mexico. These patients were diagnosed with cSLE between the years 2008 to 2022, with an average age of 11.19 ± 3.31 years, as established from the pediatric immunology service at Instituto Nacional de Pediatría, México. For inclusion, subjects were required to fulfill the criteria of the 2012 Systemic Lupus International Collaborating Clinics [15] or the 2019 European League Against Rheumatism/American College of Rheumatology Classification Criteria for Systemic Lupus Erythematosus [16]. Patients with suspected (disease onset < 5 years, parental consanguinity, or familial history of autoimmune disease in a first-degree relative) or a confirmed diagnosis of a monogenic form of cSLE were not included. Diagnosis of LN was based on the presence of proteinuria (>0.5 g/dL) or a urinary protein–creatinine (UPr–Cr) ratio > 0.2, erythrocyturia, the presence of granular casts in urinary sediment, and/or histopathological confirmation through renal biopsy.
To estimate the allelic esv3587290 CNV frequencies in our population and determine the necessary sample size, we included 181 deidentified genomic DNA samples from unrelated Mexican newborns. The DNA was obtained from residual dried blood spot (DBS) samples analyzed for newborn screening (reference group).
2.2. Genotyping of the esv3587290 CNV
Buccal cell swabs obtained from the cSLE patients, and DBS samples from the reference group, were processed by using the standard salting-out method to obtain genomic DNA. We initially performed a PCR assay using forward (5′-AGGGGAGGTGATGGACCCTA-3′) and reverse (5′-CTCAGACTGTAAGCGAAGGACA-3′) primers located inside exons 7 and 8 of VANGL1, respectively, to identify the esv3587290 CNV. This assay generated ~6 kb and ~3 kb PCR fragments from the wild-type and esv3587290 CNV alleles, respectively. This strategy was initially applied to 16 genomic DNA samples from the reference group, and the ~3 kb PCR fragments were identified in three individuals. These ~3 kb PCR fragments were gel excised, purified, and subjected to Sanger sequencing using a “primer-walking” strategy (primers are available upon request). As the entire ~3 kb sequence revealed a single identical breakpoint in these amplicons, we designed a new set of primers (VANGL1-INT7-CNVdelB-F: 5′-TGGCTGTTTCTTGTAATATCCC-3′ and VANGL1-INT7-CNVdelB-R: 5′-CCGACATGGTAAGCAAGC-3′) to amplify a shorter fragment (521 bp) encompassing the breakpoint boundaries of the esv3587290 CNV. To detect the non-deleted VANGL1 allele, we designed a set of primers (VANGL1-INT7A-F: 5′-ACTGATTGTCTGTTGATGCACATTT-3′ and VANGL1-INT7A-R: 5′-CACCCCCTAGGAGGGCAAT-3′) to amplify an internal region of intron 7 (357 bp) that is absent from the esv3587290 CNV sequence. These two mutually exclusive amplicons were generated by two separate monoplex end-point PCR assays (PCR conditions are available upon request) and resolved by agarose gel electrophoresis. Allelic and genotypic frequencies were obtained by direct counting in the cSLE patient and reference group samples as follows: wild-type homozygous, 357 bp fragment; heterozygous, 521 and 357 bp fragments; and esv3587290 CNV homozygous, 521 bp fragment. All the 521 bp esv3587290 CNV-derived amplicons from the cSLE patients and reference group individuals were subjected to direct automated Sanger sequencing and further alignment (Program Chromas Pro Version 2.1.10, Technelysium Pty Ltd., South Brisbane, QLD, Australia) with the gene (NG_016548.1) and Vang-like protein 1 isoform 1 (NM_138959.3) reference sequences to determine the precise breakpoints. The rearrangement was reported according to HGVS nomenclature (
2.3. Statistical Analysis
The sample size was calculated by applying a formula to find differences between two proportions, using reference-group allelic frequencies for the wild-type (0.765) and esv3587290 CNV (0.235) VANGL1 (Table 1), employing a confidence level of 99% and a power of 95%, and assuming a 10% loss to follow-up. From this, 33 individuals had to be included for each cSLE group (with and without LN).
The VANGL1 genotypes obtained from the reference samples and cSLE patients were tested for the Hardy–Weinberg equilibrium (HWE) using the χ2 test. Associations between the presence of LN and the presence of the esv3587290 CNV were examined by using odds ratio (OR) calculations and Pearson’s χ2 test, employing a confidence interval of 95% and a significance threshold of p ≤ 0.05. These calculations were performed using the IBM® SPSS® Statistics software Version 25.0.
This study protocol was revised and approved by the Institutional Review Research, Biosecurity, and Ethics Committees of the National Institute of Pediatrics, Mexico (Registry 2022/030; approval date 20 June 2022). We conducted this study according to the guidelines of the Declaration of Helsinki.
3. Results
LN was present in 39 of 66 (59.1%) patients at the moment of their inclusion in this study. Seventeen patients had undergone renal biopsy, and histology was available in all of them (five with class II, two with class III, eight with class IV, and two with class V of LN stratification), while thirty-three patients had 24 h proteinuria > 0.5 g/dL, thirty-four had a UPr–Cr ratio > 0.2, and sixteen had erythrocyturia and/or granular casts in the urinary sediment. Most of the patients without LN (N = 20/27) had a medical follow-up of more than 2 years, while the remaining seven underwent their last immunologist-conducted evaluation between 19 and 23 months after cSLE diagnosis.
The allelic and genotypic frequencies of the esv3587290 CNV in cSLE patients and reference samples are shown in Table 1. The proportions of genotypes in both groups were consistent with HWE, as assessed by using the χ2 test. The esv3587290 CNV alleles and cSLE patient genotypes showed significant associations with the presence of LN, as assessed using the χ2 test. More specifically, the esv3587290 CNV (OR 0.108, 95% CI 0.034–0.33, p = 0.0003) and the heterozygous genotype (OR 0.04, 95% CI 0.119–0.9811, p = 0.002) were observed to show a protective effect on the development of LN. The relatively low number of cSLE patients carrying a homozygous esv3587290 CNV genotype in the groups with and without LN precluded χ2 testing for this comparison (Table 2).
The Sanger sequencing results of all the 521 bp PCR fragments derived from the esv3587290 CNV identified in the reference and cSLE samples showed an identical breakpoint, wherein deletion of 3357 bp eliminates nearly 60% of the intron 7 sequence (5591 bp) of VANGL1 (Figure 1). We therefore designated this common deletion as NG_016548.1(NM_138959.3):c.1314+1339_1315-897del, NC_000001.10:g.116229487_116232843del (GRCh37), or NC_000001.11:g.115686866_115690222del (GRCh38) according to the HGVS nomenclature. This variant was submitted to the Leiden Open Variation Database (LOVD) of VANGL1 gene (
4. Discussion
Replication studies are mainly intended to confirm genetic associations discovered through GWAS, such as that between esv3587290 CNV and SLE-related LN [11]. These studies are needed to accumulate convincing statistical evidence that supports the association and rules out spurious findings due to uncontrolled biases [17]. To the best of our knowledge, this is the first work to explore the possible association of the esv3587290 CNV with LN in a non-European population, as previously recommended [11]. Here, we assessed an SLE population of Mexican descent, who are among the ethnicities considered to have a high risk of developing LN [18]. We further focused on a clinical form of SLE different than that previously studied in this context, namely childhood-SLE, for which LN is considered to be more prevalent and severe than in the adult form of SLE [4].
A few genetic markers have been associated with LN in Mexican populations, including SPP1 (MIM*166490) in adult SLE [allele T of rs1126616 OR 2.0 (95% CI 1.26–3.16), p = 0.003 and TT genotype under the recessive model OR 2.76 (95% CI 1.31–5.82), p = 0.011] [19] and NFE2L2 (formerly NRF2; MIM*600492) in cSLE [heterozygous A/G rs35652124 genotype OR = 1.81 (95% CI 1.04–3.12), p = 0.032) [20]. Meanwhile, other markers of LN susceptibility previously identified in European-descent populations (e.g., PDCD1, MIM*600244) [21] or murine models (CNVs of TLR7, MIM*300365) [22] have failed to show any significant association with LN in cSLE patients from Mexican populations [23,24]. These previous observations, together with our present finding that esv3587290 CNV appears to protect against the development of LN in this group of Mexican patients with cSLE (Table 2), could support the idea that LN exhibits broad genetic and phenotypic heterogeneity. Our results may also agree with the lack of an evident association between LN and VANGL1, the 1p13.1 region, or PCP pathways in a GWAS performed in Hispanic, European, African American, and Asian patients [13]. We further believe that our findings and the inability of Jiang et al. to replicate their association in the third cohort suggest that, in contrast to the previous proposal [11], esv3587290 CNV genotyping is not a viable strategy for LN risk stratification, at least for non-European SLE populations and cSLE patients.
We further characterized the precise breakpoint of the esv3587290 CNV in our study population. Whilst Jiang et al. found that this deletion varied in size when assessed using a WGS strategy [11], our Sanger sequencing approach revealed an identical deletion event in all individuals carrying one or two esv3587290 CNV copies (Figure 1). We speculate that these discrepancies may reflect different origins of the esv3587290 CNV among diverse populations, or they may be related to methodological issues. The former could imply that there is a “hot-spot” for mutational events leading to distinct rearrangements, as occur for some monogenic traits (i.e., gross deletions in the DMD gene). Alternatively, some single mutational events may occur and may even be associated with a founder effect [25], as appears to be the case for the esv3587290 CNV in our Mexican population. Further haplotypic analysis is warranted to support the latter hypothesis. Regarding the potential impact of methodologic issues, we note that the short-read WGS strategy is generally intended to approximately localize the breakpoints of a gross genomic structural rearrangement (e.g., a deletion-type CNV); to reach a nucleotide-level resolution, it would be necessary to apply long-range PCR and Sanger sequencing [26,27], as performed herein. Given this, we propose that it would be desirable to determine the precise nucleotide-level rearrangements of the esv3587290 CNV in other populations. In the PCR-based Sanger sequencing strategy, we used the forward and reverse primers designed to anneal ~250 bp away from the breakpoints estimated by WGS to amplify the esv3587290 CNV [11]. There remains some possibility that allelic drop-out may have occurred due to the non-amplification of alleles bearing different breakpoints. However, we believe that this is unlikely due to the lack of departure from HWE in both study groups, along with the similarity between our allelic frequencies (Table 1) and those previously reported for esv3587290 CNV in Latino populations (~0.3) [11].
Although mis-spliced mRNA species of VANGL1 lacking exon 2 were identified in peripheral blood mononuclear cells from two of the homozygous esv3587290 CNV SLE patients, whether the esv3587290 CNV has any effect on VANGL1 function at the kidney remains to be determined [11]. Such a finding would support the idea that this CNV plays a role in genetic susceptibility to the development of LN in humans. Interestingly, a predicted enhancer element is located inside exon 8 (encoding the 3′UTR) of VANGL1 (
Limitations of the Study
Nonetheless, our study had several limitations: It is possible that our study sample was under-representative or phenotypically heterogeneous, given the lack of renal biopsy-proven LN in most of our patients (N = 22/39 LN cSLE patients, 56.4%), which is still considered the diagnostic gold standard [7]. Obtaining such information could allow us to perform an association analysis stratified by histopathological classes, as previously recommended [11]. Also, we did not perform any ancestry analysis to determine if the reference group and cSLE patients were admixed; thus, a possibility of bias by population stratification cannot be excluded. The low number of male cSLE patients (nine with LN and two without LN) did not allow us to perform stratification analysis by gender, so this aspect should be addressed in the future. It is important to note that the above-described aspects were also not considered by Jiang et al. [11]. Genotyping errors could have biased our analysis, as we did not use a second molecular technique to validate our PCR-based genotyping assay. However, this seems unlikely given the lack of esv3587290 CNV departure from HWE in the reference group and cSLE patients, along with the similarity between the observed allelic frequencies (Table 1) and those previously reported in Latino populations [11]. Finally, the cross-sectional nature of this study precluded us from determining whether the included cSLE patients without LN could later develop nephritis, particularly in those seven cSLE patients with normal kidney function who did not complete the two years of follow-up, adding another potential confounding factor. The estimated risk of developing LN in these patients is expected to be 7% [28].
5. Conclusions
The esv3587290 CNV of the VANGL1 gene was not associated with the development of LN in a sample of Mexican cSLE patients; rather, this CNV seems to be a genetic protective factor. Further replication studies on the esv3587290 CNV in other ethnicities and clinical forms of SLE are warranted to define its role as a genetic factor in the development of LN. The esv3587290 CNV seems to be a unique 3357 bp deletion that may have originated from a single mutational event in our Mexican population.
Conceptualization, M.A.A.-O.; methodology, M.A.A.-O., L.D.-G. and B.E.-O.; software, M.A.A.-O. and L.D.-G.; validation, M.A.A.-O., A.L.R.-L., L.D.-G. and A.G.-d.A.; formal analysis, M.A.A.-O., A.L.R.-L., L.D.-G., R.G.N.-V. and A.G.-d.A.; investigation, M.A.A.-O., B.E.-O., A.L.R.-L., F.E.R.-L., R.G.N.-V. and A.G.-d.A.; resources, M.A.A.-O., A.L.R.-L., F.E.R.-L. and A.G.-d.A.; data curation, M.A.A.-O., A.L.R.-L., B.E.-O. and A.G.-d.A.; writing—original draft preparation, M.A.A.-O.; writing—review and editing, M.A.A.-O., B.E.-O., A.L.R.-L., L.D.-G., F.E.R.-L., R.G.N.-V. and A.G.-d.A.; visualization, M.A.A.-O.; supervision, M.A.A.-O.; project administration, M.A.A.-O.; funding acquisition, M.A.A.-O. All authors have read and agreed to the published version of the manuscript.
This study was conducted according to the guidelines of the Declaration of Helsinki, and was approved by the Institutional Review Research, Biosecurity, and Ethics Committees of the National Institute of Pediatrics, Mexico (Registry 2022/030). The approval date was 20 June 2022.
Written informed consent to participate in this study and to publish this paper was obtained either from the patients involved in the study or their parents.
Publicly available datasets were analyzed in this study. These data can be found in the LOVD v.3.0—Leiden Open Variation Database:
The authors gratefully acknowledge the patients and parents for their commitment. They also thank Roberto Edher Demetrio-Ríos and Karen Jimena Guerrero González for their assistance with the recruitment and patient sampling processes, as well as data collection and technical support.
The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analysis, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.
Footnotes
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.
Allelic and genotypic frequencies of esv3587290 CNV in cSLE patients and reference samples.
VANGL1 | Total cSLE Patients | Reference Group |
---|---|---|
Wild-type allele | 0.712 | 0.765 |
esv3587290 CNV- | 0.288 | 0.235 |
VANGL1 | Total cSLE Patients | Reference Group |
Wild-type homozygous | 48.5% (N = 32) | 58% (N = 105) |
Heterozygous | 45.5% (N = 30) | 37% (N = 67) |
Homozygous esv3587290 | 6% (N = 4) | 5% (N = 9) |
Abbreviations: cSLE: childhood-onset systemic lupus erythematosus; CNV: copy number variation; DBS: residual anonymized dried blood spots from Mexican unrelated neonates subjected to newborn screening program.
Association results and VANGL1 allelic and genotypic frequencies among cSLE patients with or without LN.
VANGL1 | Allelic Frequencies | Allelic Frequencies | χ2/ | p-Value |
---|---|---|---|---|
Wild-type allele | 0.769 | 0.630 | χ2 = 16.41 | p = 0.000025 |
esv3587290 CNV | 0.231 | 0.370 | ||
VANGL1 | Genotypic Frequencies | Genotypic Frequencies | χ2/ | p-Value |
Wild-type | 59% (N = 23) | 33.3% (N = 9) | ||
Heterozygous esv3587290 | 35.9% (N = 14) | 59.3% (N = 16) | χ2 = 4.089 2 | p = 0.02 |
Homozygous esv3587290 | 5.1% (N = 2) | 7.4% (N = 2) | χ2 = not calculable 3 | p = not calculable |
1 The presence of LN, esv3587290 alleles, and heterozygous genotypes were taken as the reference categories for the OR calculations. 2 Wild-type homozygous vs. heterozygous for esv3587290; 3 Wild-type homozygous vs. homozygous for esv3587290. Abbreviations: CI: confidence interval; cSLE: childhood-onset systemic lupus erythematosus; CNV: copy number variation; OR: odds ratio.
References
1. Barsalou, J.; Levy, D.M.; Silverman, E.D. An update on childhood-onset systemic lupus erythematosus. Curr. Opin. Rheumatol.; 2013; 25, pp. 616-622. [DOI: https://dx.doi.org/10.1097/BOR.0b013e328363e868] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23836073]
2. Brunner, H.I.; Holland, M.J.; Beresford, M.W.; Ardoin, S.P.; Appenzeller, S.; Silva, C.A.; Flores, F.; Goilav, B.; Aydin, P.O.A.; Wenderfer, S.E. et al. American College of Rheumatology Provisional Criteria for Clinically Relevant Improvement in Children and Adolescents with Childhood-Onset Systemic Lupus Erythematosus. Arthritis Care Res.; 2019; 71, pp. 579-590. [DOI: https://dx.doi.org/10.1002/acr.23834] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30680946]
3. Ambrose, N.; Morgan, T.A.; Galloway, J.; Ionnoau, Y.; Beresford, M.W.; Isenberg, D.A. Differences in disease phenotype and severity in SLE across age groups. Lupus; 2016; 25, pp. 1542-1550. [DOI: https://dx.doi.org/10.1177/0961203316644333] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27147622]
4. Fatemi, A.; Matinfar, M.; Smiley, A. Childhood versus adult-onset systemic lupus erythematosus: Long-term outcome and predictors of mortality. Clin. Rheumatol.; 2017; 36, pp. 343-350. [DOI: https://dx.doi.org/10.1007/s10067-016-3509-1] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28012055]
5. Borchers, A.T.; Leibushor, N.; Naguwa, S.M.; Cheema, G.S.; Shoenfeld, Y.; Gershwin, M.E. Lupus nephritis: A critical review. Autoimmun. Rev.; 2012; 12, pp. 174-194. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22982174][DOI: https://dx.doi.org/10.1016/j.autrev.2012.08.018]
6. Alduraibi, F.K.; Tsokos, G.C. Lupus Nephritis Biomarkers: A Critical Review. Int. J. Mol. Sci.; 2024; 25, 805. [DOI: https://dx.doi.org/10.3390/ijms25020805] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/38255879]
7. Abulaban, K.M.; Brunner, H.I. Biomarkers for childhood-onset systemic lupus erythematosus. Curr. Rheumatol. Rep.; 2015; 17, 471. [DOI: https://dx.doi.org/10.1007/s11926-014-0471-2] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25475594]
8. Picard, C.; Lega, J.C.; Ranchin, B.; Cochat, P.; Cabrera, N.; Fabien, N.; Belot, A. Anti-C1q autoantibodies as markers of renal involvement in childhood-onset systemic lupus erythematosus. Pediatr. Nephrol.; 2017; 32, pp. 1537-1545. [DOI: https://dx.doi.org/10.1007/s00467-017-3646-z] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28343355]
9. Vinuesa, C.G.; Shen, N.; Ware, T. Genetics of SLE: Mechanistic insights from monogenic disease and disease-associated variants. Nat. Rev. Nephrol.; 2023; 19, pp. 558-572. [DOI: https://dx.doi.org/10.1038/s41581-023-00732-x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/37438615]
10. Iwamoto, T.; Niewold, T.B. Genetics of human lupus nephritis. Clin. Immunol.; 2017; 185, pp. 32-39. [DOI: https://dx.doi.org/10.1016/j.clim.2016.09.012] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27693588]
11. Jiang, S.H.; Mercan, S.; Papa, I.; Moldovan, M.; Walters, G.D.; Koina, M.; Fadia, M.; Stanley, M.; Lea-Henry, T.; Cook, A. et al. Deletions in VANGL1 are a risk factor for antibody-mediated kidney disease. Cell Rep. Med.; 2021; 2, 100475. [DOI: https://dx.doi.org/10.1016/j.xcrm.2021.100475] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35028616]
12. Chung, S.A.; Brown, E.E.; Williams, A.H.; Ramos, P.S.; Berthier, C.C.; Bhangale, T.; Alarcon-Riquelme, M.E.; Behrens, T.W.; Criswell, L.A.; Graham, D.C. et al. Lupus nephritis susceptibility loci in women with systemic lupus erythematosus. J. Am. Soc. Nephrol.; 2014; 25, pp. 2859-2870. [DOI: https://dx.doi.org/10.1681/ASN.2013050446] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24925725]
13. Lanata, C.M.; Nititham, J.; Taylor, K.E.; Chung, S.A.; Torgerson, D.G.; Seldin, M.F.; Pons-Estel, B.A.; Tusié-Luna, T.; Tsao, B.P.; Morand, E.F. et al. Genetic contributions to lupus nephritis in a multi-ethnic cohort of systemic lupus erythematous patients. PLoS ONE; 2018; 13, e0199003. [DOI: https://dx.doi.org/10.1371/journal.pone.0199003] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29953444]
14. den Dunnen, J.T.; Dalgleish, R.; Maglott, D.R.; Hart, R.K.; Greenblatt, M.S.; McGowan-Jordan, J.; Roux, A.F.; Smith, T.; Antonarakis, S.E.; Taschner, P.E. HGVS Recommendations for the Description of Sequence Variants: 2016 Update. Hum. Mutat.; 2016; 37, pp. 564-569. [DOI: https://dx.doi.org/10.1002/humu.22981] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26931183]
15. Petri, M.; Orbai, A.M.; Alarcón, G.S.; Gordon, C.; Merrill, J.T.; Fortin, P.R.; Bruce, I.N.; Isenberg, D.; Wallace, D.J.; Nived, O. et al. Derivation and validation of the Systemic Lupus International Collaborating Clinics classification criteria for systemic lupus erythematosus. Arthritis Rheum.; 2012; 64, pp. 2677-2686. [DOI: https://dx.doi.org/10.1002/art.34473] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22553077]
16. Aringer, M.; Costenbader, K.; Daikh, D.; Brinks, R.; Mosca, M.; Ramsey-Goldman, R.; Smolen, J.S.; Wofsy, D.; Boumpas, D.T.; Kamen, D.L. et al. 2019 European League Against Rheumatism/American College of Rheumatology Classification Criteria for Systemic Lupus Erythematosus. Arthritis Rheumatol.; 2019; 71, pp. 1400-1412. [DOI: https://dx.doi.org/10.1002/art.40930] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31385462]
17. Kraft, P.; Zeggini, E.; Ioannidis, J.P. Replication in genome-wide association studies. Stat. Sci.; 2009; 24, pp. 561-573. [DOI: https://dx.doi.org/10.1214/09-STS290] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20454541]
18. Sánchez, E.; Rasmussen, A.; Riba, L.; Acevedo-Vasquez, E.; Kelly, J.A.; Langefeld, C.D.; Williams, A.H.; Ziegler, J.T.; Comeau, M.E.; Marion, M.C. et al. Impact of genetic ancestry and sociodemographic status on the clinical expression of systemic lupus erythematosus in American Indian-European populations. Arthritis Rheum.; 2012; 64, pp. 3687-3694. [DOI: https://dx.doi.org/10.1002/art.34650] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22886787]
19. Rivera-Cameras, A.; Gallegos-Arreola, M.P.; Morán-Moguel, M.C.; Salazar-Páramo, M.; Alcaraz-López, M.F.; Echeverría-González, G.; Topete-Reyes, J.F.; Franco-Chávez, S.A.; Dávalos-Rodríguez, I.P. Association of the rs1126616 and rs9138 Variants in the SPP1 Gene among Mexican Patients with Systemic Lupus Erythematosus and Lupus Nephritis. Int. J. Mol. Sci.; 2024; 25, 1000. [DOI: https://dx.doi.org/10.3390/ijms25021000] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/38256074]
20. Córdova, E.J.; Velázquez-Cruz, R.; Centeno, F.; Baca, V.; Orozco, L. The NRF2 gene variant, -653G/A, is associated with nephritis in childhood-onset systemic lupus erythematosus. Lupus; 2010; 19, pp. 1237-1242. [DOI: https://dx.doi.org/10.1177/0961203310367917] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20507872]
21. Johansson, M.; Arlestig, L.; Möller, B.; Rantapää-Dahlqvist, S. Association of a PDCD1 polymorphism with renal manifestations in systemic lupus erythematosus. Arthritis Rheum.; 2005; 52, pp. 1665-1669. [DOI: https://dx.doi.org/10.1002/art.21058] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15934088]
22. Subramanian, S.; Tus, K.; Li, Q.Z.; Wang, A.; Tian, X.H.; Zhou, J.; Liang, C.; Bartov, G.; McDaniel, L.D.; Zhou, X.J. et al. A Tlr7 translocation accelerates systemic autoimmunity in murine lupus. Proc. Natl. Acad. Sci. USA; 2006; 103, pp. 9970-9975. [DOI: https://dx.doi.org/10.1073/pnas.0603912103] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16777955]
23. Velázquez-Cruz, R.; Orozco, L.; Espinosa-Rosales, F.; Carreño-Manjarrez, R.; Solís-Vallejo, E.; López-Lara, N.D.; Ruiz-López, I.K.; Rodríguez-Lozano, A.L.; Estrada-Gil, J.K.; Jiménez-Sánchez, G. et al. Association of PDCD1 polymorphisms with childhood-onset systemic lupus erythematosus. Eur. J. Hum. Genet.; 2007; 15, pp. 336-341. [DOI: https://dx.doi.org/10.1038/sj.ejhg.5201767] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/17228327]
24. García-Ortiz, H.; Velázquez-Cruz, R.; Espinosa-Rosales, F.; Jiménez-Morales, S.; Baca, V.; Orozco, L. Association of TLR7 copy number variation with susceptibility to childhood-onset systemic lupus erythematosus in Mexican population. Ann. Rheum. Dis.; 2010; 69, pp. 1861-1865. [DOI: https://dx.doi.org/10.1136/ard.2009.124313] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20525845]
25. Shotelersuk, V.; Larson, D.; Anikster, Y.; McDowell, G.; Lemons, R.; Bernardini, I.; Guo, J.; Thoene, J.; Gahl, W.A. CTNS mutations in an American-based population of cystinosis patients. Am. J. Hum. Genet.; 1998; 63, pp. 1352-1362. [DOI: https://dx.doi.org/10.1086/302118] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/9792862]
26. Ma, S.; Zhang, Z.; Fu, Y.; Zhang, M.; Niu, Y.; Li, R.; Guo, Q.; He, Z.; Zhao, Q.; Song, Z. et al. Identification of the first Alu-mediated gross deletion involving the BCKDHA gene in a compound heterozygous patient with maple syrup urine disease. Clin. Chim. Acta.; 2021; 517, pp. 23-30. [DOI: https://dx.doi.org/10.1016/j.cca.2021.01.023] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33607070]
27. Cuenca-Guardiola, J.; de la Morena-Barrio, B.; García, J.L.; Sanchis-Juan, A.; Corral, J.; Fernández-Breis, J.T. Improvement of large copy number variant detection by whole genome nanopore sequencing. J. Adv. Res.; 2023; 50, pp. 145-158. [DOI: https://dx.doi.org/10.1016/j.jare.2022.10.012] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/36323370]
28. Pennesi, M.; Benvenuto, S. Lupus Nephritis in Children: Novel Perspectives. Medicina; 2023; 59, 1841. [DOI: https://dx.doi.org/10.3390/medicina59101841] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/37893559]
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
A ~3-kb deletion-type DNA copy number variation (CNV, esv3587290) located at intron 7 of the VANGL1 gene (1p13.1, MIM*610132) has been proposed as a genetic factor in lupus nephritis (LN) development in adult systemic lupus erythematosus (SLE) patients across European-descent populations, but its replication in other ethnicities has been inconsistent and its association with LN in childhood-onset SLE (cSLE) remains unknown. Here, we performed an exploratory association study in a sample of 66 unrelated cSLE Mexican patients (11 males, 55 females; ages 7.8 to 18.6 years). Two stratified groups were compared: cSLE patients with (N = 39) or without (N = 27) LN, as diagnosed by renal biopsy (N = 17), proteinuria (N = 33), urinary protein–creatinine ratio > 0.2 (N = 34), and erythrocyturia and/or granular casts in urinary sediment (N = 16). For esv3587290 CNV genotyping, we performed an end-point PCR assay with breakpoint confirmation using Sanger sequencing. We also determined the allelic frequencies of the esv3587290 CNV in 181 deidentified ethnically matched individuals (reference group). The obtained genotypes were tested for Hardy–Weinberg equilibrium using the χ2 test. Associations between LN and esv3587290 CNV were tested by calculating the odds ratio (OR) and using Pearson’s χ2 tests, with a 95% confidence interval and p ≤ 0.05. The esv3587290 CNV allele (OR 0.108, 95% CI 0.034–0.33, p = 0.0003) and the heterozygous genotype (OR 0.04, 95% CI 0.119–0.9811, p = 0.002) showed a significant protective effect against LN development. Finally, we characterized the precise breakpoint of the esv3587290 CNV to be NG_016548.1(NM_138959.3):c.1314+1339_1315-897del in our population. This report supports the notion that a broad genetic heterogeneity underlies the susceptibility for developing LN.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details






1 Laboratorio de Biología Molecular, Subdirección de Investigación Médica, Instituto Nacional de Pediatría, Ciudad de México 04530, Mexico;
2 Servicio de Inmunología, Instituto Nacional de Pediatría, Ciudad de México 04530, Mexico;
3 Departamento de Metodología de la Investigación, Subdirección de Investigación Médica, Instituto Nacional de Pediatría, Ciudad de México 04530, Mexico;