Peizhi Deng
Zhixin Li
Bin Yi3*- 1Clinical Research Center, The Third Xiangya Hospital, Central South University, Changsha, Hunan, China
- 2Xiangya School of Medicine, Central South University, Changsha, Hunan, China
- 3Department of Nephrology, The Third Xiangya Hospital, Central South University, Changsha, Hunan, China
- 4The Affiliated Changsha Central Hospital, Research Center for Phase I Clinical Trials, Hengyang Medical School, University of South China, Changsha, Hunan, China
Background: The prevalence of immunoglobulin A nephropathy (IgAN) seems to be higher in patients with type 1 diabetes mellitus (T1DM) than that in the general population. However, whether there exists a causal relationship between T1DM and IgAN remains unknown.
Methods: This study conducted a standard two-sample Mendelian randomization (MR) analysis to assess the causal inference by four MR methods, and the inverse variance-weighted (IVW) approach was selected as the primary method. To further test the independent causal effect of T1DM on IgAN, multivariable MR (MVMR) analysis was undertaken. Sensitivity analyses incorporating multiple complementary MR methods were applied to evaluate how strong the association was and identify potential pleiotropy.
Results: MR analyses utilized 81 single-nucleotide polymorphisms (SNPs) for T1DM. The evidence supports a significant causal relationship between T1DM and increased risk of IgAN [odds ratio (OR): 1.39, 95% confidence interval (CI): 1.10–1.74 for IVW, p < 0.05]. The association still exists after adjusting for triglyceride (TG), fasting insulin (FI), fasting blood glucose (FBG), homeostasis model assessment of beta-cell function (HOMA-B) and insulin resistance (HOMA-IR), and glycated hemoglobin (HbA1c). MVMR analysis indicated that the effect of T1DM on IgAN vanished upon accounting for low-density lipoprotein cholesterol (LDL-c; OR: 0.97, 95% CI: 0.90–1.05, p > 0.05).
Conclusions: This MR study provided evidence that T1DM may be a risk factor for the onset of IgAN, which might be driven by LDL-c. Lipid-lowering strategies targeting LDL-c should be enhanced in patients with T1DM to prevent IgAN.
Introduction
Type 1 diabetes mellitus (T1DM), also known as autoimmune diabetes, is characterized by irreversible islet beta-cell destruction. According to prior research, T1DM is a major factor contributing to diabetic nephropathy (DN), and it also carries a significant risk of developing additional renal comorbidities. These past few years, the accurate and early diagnosis of non-diabetic renal disease (NDRD) has aroused widespread concern in clinical practice (1). In addition to being the most prevalent type of primary glomerulonephritis, immunoglobulin A nephropathy (IgAN) is also regarded as a major form of NDRD, especially among young adults, and responsible for the global burden of chronic kidney disease (CKD) and renal failure (2–4). Additionally, the prevalence of IgAN is still increasing, whereas the geographic distribution differs greatly (5). IgAN presented a significant prevalence in European countries, especially in France (52.7% for primary glomerular disease (PGD)), Germany (50.7%), United Kingdom (39.0%), and Czechia (37.4%) (6–10). Hence, the findings regarding the pathogenesis of IgAN and its risk factors may offer a crucial and original concept for the prevention and treatment of NDRD, particularly in places with a high prevalence.
To the best of our knowledge, past epidemiological studies have already supported an underlying correlation between T1DM and IgAN. A meta-analysis of 48 studies, including 4,876 patients who suffered from DM, has shown that IgAN (3%–59%) is one of the most common NDRD types for patients with diabetes (11). However, the participants enrolled in these studies had either T1DM or T2DM at baseline; whether there is a direct causal relationship between T1DM and IgAN is still a mystery. Furthermore, most of the existing small-sample single-centered studies are limited in assessing the potential comorbid confounders, necessitating further investigation into the precise association between T1DM and IgAN.
A better approach for investigating this is Mendelian randomization (MR), in which the causal association between modifiable exposure and disease outcome can be well evaluated using single-nucleotide polymorphisms (SNPs) as robust instrument variables (IVs) (12, 13). Since the genetic instrument is fixed at conception, this design is less susceptible to confounding and reverse causality bias. Here, we used two-sample MR analyses to estimate the probability that T1DM and IgAN might be primarily related.
Materials and methods
Study design
We performed this study using a conventional two-sample MR design, assessing the association between T1DM and IgAN, and Figure 1 provides a summary of the design that is being used for the current investigation. Firstly, we extracted available genetic IVs from the T1DM meta-analysis. Secondly, the summary data comprising all SNPs from the large-scale genome-wide association study (GWAS) for IgAN were collected. To evaluate the causal effect, we employed univariate two-sample MR and various sensitivity analyses. In addition, we conducted a multivariable MR (MVMR) analysis controlling for confounding variables. The Strengthening the Reporting of Observational Studies in Epidemiology-Mendelian Randomization (STROBE-MR) checklist is presented in Supplementary Table S1.
Figure 1 Diagram of Mendelian randomization framework in this study. IgAN, immunoglobulin A nephropathy; IVW, inverse variance-weighted; MR, Mendelian randomization; PRESSO, Pleiotropy RESidual Sum and Outlier; SNP, single nucleotide polymorphism; T1DM, type 1 diabetes mellitus.
The data source for type 1 diabetes mellitus
All identified IVs for T1DM used the data from the latest and largest GWAS, in which 18,942 T1DM patients and 501,638 control participants of European ancestry from nine cohorts were enrolled (14). After applying uniform quality control, Taliun et al. (15) imputed genotypes into the TOPMed reference panel and tested for the association with T1DM. Through the meta-analysis, they combined the association results for 61,947,369 variants and identified 81 SNPs that reached genome-wide significance [p < 5 × 10−8, linkage disequilibrium (LD); r2 < 0.001, LD distance >10,000 kb], including 48 of 59 known loci and 33 loci that have not previously been reported. A total of 81 SNPs explained 0.04 of the variance in T1DM. An F-statistic >10 indicates a significant association between the selected IVs and T1DM (16). For non-Europeans, summary statistics of T1DM were obtained from a GWAS that consisted of 1,219 East Asian ancestry cases and 132,032 controls (17). A total of 16 SNPs utilized as IVs must be strongly linked to T1DM (p < 5 × 10-8). SNPs were filtered by employing the clumping technique to ensure independence (18). Palindromic SNPs with a minor allele frequency (MAF) of <0.42 were eliminated.
The data sources for IgA nephropathy
The Medical Research Council (MRC)/Kidney Research UK National DNA Bank for Glomerulonephritis established collections in five common glomerular diseases, of which one is IgAN. The diagnosis was confirmed in all IgAN patients by direct review of renal biopsy histopathology reports and clinical IgAN records. Individuals with evidence of liver disease or Henoch–Schönlein purpura were excluded. A GWAS meta-analysis encompassing roughly 5,957 European participants yielded summary-level information regarding the genetic associations with IgAN, including 977 cases from the MRC/Kidney Research UK National DNA Bank and 4,980 controls from the 1958 British Birth Cohort and the UK National Blood Service (19). To our knowledge, the maximum sample overlapping rate between exposure data and outcome data is <1%. Summary data of IgAN were acquired for non-Europeans from a GWAS that included 71 East Asian ancestry patients and 175,288 controls (17). All of the data were taken from previously published studies that were made available to the public. As a result, neither ethical approval nor patient consent were needed for the study.
Statistical analysis
Two-sample MR analyses
There are three important assumptions of conventional MR analysis. The ideal IVs must satisfy the following: 1) must be truly associated with T1DM (in this study, defined as the genetic association with p < 5 × 10-8); 2) not associated with confounders of the T1DM–IgAN association; 3) should only be associated with the IgAN through T1DM (20).
Four different methods, namely, inverse variance-weighted (IVW), MR-Egger, weighted median, and maximum likelihood, were used in the two-sample MR analysis to evaluate the causal effects between T1DM and IgAN. Among them, IVW was used as the primary statistical approach in our study, since it is the most accurate method for estimating causal effects (21). MR-Egger may be susceptible to significant impacts of atypical genetic factors, resulting in inaccurate estimations. However, the MR-Egger test can provide non-biased estimates even if all selected IVs are invalid (22). The weighted median estimator could provide a consistent causal estimate, even when up to 50% of the information from IVs is nonfunctional. When the IVs are weak, IVW ignores the true variance of the estimate, whereas maximum likelihood provides confidence intervals (CIs) that are accurately assessed (23). The results were shown as odds ratio (OR) and it’s 95% CI. The statistical power of the MR analysis was calculated using the mRnd tool (http://cnsgenomics.com/shiny/mRnd/) (24). The results of power calculations are shown in Supplementary Table S2. In addition, the causal relationship between IgAN and T1DM was assessed using reverse MR analysis by employing data from the European and East Asian populations described above (17, 19). In order to obtain more comprehensive results, selected SNPs with a GWAS p-value cutoff lower than 5 × 10-6 were selected as IVs.
Sensitivity univariable MR analyses
We used first-order IVWs and MR-Egger to generate Cochran’s Q test to check for heterogeneity, which represents a possible violation of modeling assumptions (25). This study used the MR-Egger regression intercept examination to estimate the potential pleiotropy between exposure and outcome (26). A p-value <0.05 represented the existence of pleiotropy. Once heterogeneity or horizontal pleiotropy was noteworthy, we used MR-Pleiotropy RESidual Sum and Outlier (MR-PRESSO) to remove outlier SNPs (24). Moreover, we conducted the “leave-one-out” test to determine whether a single SNP had a significant independent effect on MR estimates. Autoimmune diseases share common genetic loci in the HLA, and this could lead to pleiotropic bias (27). Hence, we excluded the SNPs from the HLA complex and repeated the MR analysis.
Multivariable two-sample MR analyses
Previous clinical studies have demonstrated that glycolipid metabolic traits show a certain degree of correlation with T1DM and IgAN (28–33). We added summary GWAS data on these confounders and chose MVMR to assess the direct effect of T1DM while accounting for the effect of these metabolic traits. All of the summary data utilized in the study are publicly available, and the detailed information is presented in Table 1.
Table 1 Description of genome-wide association study used for each confounder.
All analyses were carried out with the packages “TwoSampleMR” (34), “MR-PRESSO” (24), and “MVMR” (35) of R (version 4.2.0). All presented p was two-sided, and statistical significance was set at the 5% level.
Results
The characteristics of SNP and participants for analyses
Table 2 presents the characteristics of the populations included in the GWAS data on exposure and outcome. A total of 81 SNPs were chosen as IVs for T1DM in the primary analysis. The F-statistic ranges from 43 to 249, reflecting a strong instrument strength for T1DM. The summary of the identified SNPs in the MR analysis is presented in Supplementary Table S3.
Table 2 Characteristics of type 1 diabetes mellitus and IgA nephropathy.
Associations between type 1 diabetes mellitus and IgA nephropathy
Univariable and multivariable MR analyses
In the main univariable analyses, we identified a significant causal relationship between exposure and outcome (p < 0.05 across four MR methods), which referred to a causal association between T1DM and increased IgAN risk (OR: 1.39, 95% CI: 1.10–1.74 for IVW; OR: 1.53, 95% CI: 1.10–2.12 for weighted median; OR: 1.39, 95% CI: 1.13–1.71 for maximum likelihood) (Table 3). Similarly, T1DM was associated with an increased IgAN risk in the MVMR analyses with statistical power [i.e., OR: 1.37, 95% CI: 1.08–1.74, adjusted for triglyceride (TG); OR: 1.38, 95% CI: 1.13–1.71, adjusted for fasting insulin (FI); OR: 1.42, 95% CI: 1.14–1.77, adjusted for fasting blood glucose (FBG); OR: 1.42, 95% CI: 1.15–1.76, adjusted for homeostasis model assessment of beta-cell function (HOMA-B); OR: 1.39, 95% CI: 1.12–1.72, adjusted for HOMA of insulin resistance (HOMA-IR); OR: 1.39, 95% CI: 1.12–1.73, adjusted for glycated hemoglobin (HbA1c)] (Table 3). However, no similar trend was shown in low-density lipoprotein cholesterol (LDL-c; OR: 0.97, 95% CI: 0.90–1.05). In East Asian participants, the MR methods did not show that genetically determined T1DM was associated with IgAN (all p > 0.05, Supplementary Table S4). In addition, the reverse MR analysis did not demonstrate that a genetic predisposition to IgAN was associated with T1DM (Supplementary Table S5).
Table 3 Univariable and multivariable two-sample Mendelian randomization estimations showing the effect of type 1 diabetes mellitus on the risk of IgA nephropathy.
Sensitivity analysis
Additionally, we conducted several sensitivity analyses to determine potential heterogeneity and horizontal pleiotropy (Table 3). Neither the Cochran’s Q test nor the MR-Egger regression analysis detected heterogeneity and horizontal pleiotropy in the primary analysis (MR-Egger Ph = 0.19, IVW Ph = 0.20, Pintercept = 0.47). We obtained consistent estimates after removing three SNPs (rs2395471, rs12665124, and rs2523679) from the human leukocyte antigen (HLA) complex (Supplementary Table S6). Moreover, the MR-PRESSO data did not reveal any outlier SNP. Four methods were used to evaluate the results of the MR analysis, and the scatter plot was generated (Figure 2A). In the leave-one-out analysis, we discovered that no single SNP drove the overall effect of T1DM on IgAN (Figure 2B). Figure 2C illustrated a relatively symmetrical distribution of variant effects for IgAN, indicating an absence of directional pleiotropy.
Figure 2 Scatter plot (A), leave-one-out test (B), and funnel plot (C) for genetically determined type 1 diabetes mellitus on IgA nephropathy risk. MR, Mendelian randomization.
Discussion
Our work utilized the large-scale GWAS data to investigate the effect of genetically predicted T1DM on IgAN risk within the MR framework and provided evidence supporting the causal effects of T1DM on IgAN, independently of a wide range of potential confounders, including TG, FI, FBG, HOMA-B, HOMA-IR, and HbA1c. However, MVMR analysis indicated the effect of T1DM on IgAN that vanished upon accounting for LDL-c.
DN, a devastating complication in diabetic patients, is regarded as the leading cause of end-stage renal disease worldwide for a long time (36). However, previous studies have shown that long-term administration of insulin to people with T1DM can cause kidney damage via a number of signal pathways, including the protein kinase B-mechanistic target of rapamycin kinase C1 (Akt-mTORC1) pathway, which makes the spectrum of renal disease among these patients complex (37). Many diabetic individuals have been found to develop IgAN, which is primarily characterized by the mesangial deposition of IgA immune complexes with the possible pathogenesis of galactose-deficient IgA1 overproduction (3, 38). And the most prevalent clinical presentations of IgAN, such as asymptomatic hematuria and progressive kidney disease, are thought to significantly increase the burden of CKD and renal failure (39). However, due to the unsystematic kidney biopsy, there is a relatively high proportion of NDRD represented by IgAN that has been misdiagnosed as DN, which hinders our ability to fully understand the association between diabetes and other renal diseases.
Interestingly, according to a case report in the 1990s, Gans et al. (40) suggested that the occurrence of IgAN in diabetic patients had to be more than coincidental because individuals with T1DM also exhibited a high prevalence of dermatitis herpetiformis and celiac disease, two autoimmune diseases involving pathogenic IgA antibodies (41). In addition, current research indicates that both T1DM and IgAN have a substantial relationship with HLA (19, 42). Some IgAN susceptibility loci are also associated with the risk of other autoimmune disorders such as T1DM (43). Given the foregoing and the obvious defects in epidemiological studies, such as unknown confounding factors, measurement error, and reverse causal effects from the environmental data, we believed that the gene may serve as the key to understanding the relationship and therefore conducted the reverse MR analysis to investigate the causal effects between T1DM and IgAN.
This study provided evidence that T1DM has a causal relationship with IgAN; however, we failed to find evidence that IgAN directly causes T1DM, suggesting that it may be mediated by other pathways. To our knowledge, there is no study available on the T1DM risk profile in individuals with IgAN, so there still remains a broad research space for an investigation into the underlying pathway, which may involve medication, inflammation, and endocrine alterations. Furthermore, given the complexity of associations in the HLA region, a more thorough sequencing may be required to depict a comprehensive gene picture of autoimmune diseases. In MVMR analysis, we adjusted for some potential comorbid confounders, among which LDL-c did not work when considering the causal role. Even though there is no precise mechanism to explain the effect of LDL-c on the association of T1DM and IgAN, numerous clinical studies indicated that the prevalence of dyslipidemia is high in patients with T1DM (44), and an unfavorable lipid profile, particularly characterized by high LDL-c concentrations, may increase the risk of renal diseases (45). Also, the lipid nephrotoxicity hypothesis has demonstrated that oxidized LDL-c could enhance macrophage infiltration, accelerate inflammatory response, and promote glomerular sclerosis (46, 47). And decreased oxidative resistance of LDL-c was found to be a hallmark of IgAN (48). Our findings might spark interest in further investigation and help reveal the significance of implementing a targeted lipid-lowering strategy in clinical practice to improve IgAN.
There are multiple strengths in this study. The primary merit is originality. We implemented a two-sample MR to explore the underlying relationship between T1DM and IgAN based on the genetic background, filling a previously unknown research gap. Additionally, the MR including multiple approaches is scientific and reasonable, which strengthened the causal inference between exposure and disease outcome. In addition, a variety of sensitivity analyses incorporating multiple complementary MR approaches were implemented to test the robustness of the association and determine potential bias from pleiotropy. Moreover, we also accounted for seven potential confounding variables by reanalyzing GWAS summary data and did not detect bias.
There are also some limitations of the analysis. The studied population only includes individuals of European ancestry and East Asian ancestry, thus, the results in this study cannot be generalized among populations of other ancestries. Also, there exists insufficient statistical power in Asians due to the small number of IgAN cases and the partially overlapping sample between T1DM and IgAN. Considering that the prevalence of IgAN exhibits epidemiological variability, more research is necessitated to identify whether there is a regional or racial difference. Additionally, there was a risk of false-positive in our study, which may be caused by the overlapping samples between exposure and outcome data. However, after calculation, we obtained a maximum overlapping rate of <1%, which may not have been sufficient to affect our results. Moreover, due to the lack of basic research regarding the effect of T1DM on the onset of IgAN, we were unable to fully adjust for confounding variables and rule out their effect of them. Thus, the molecular or cellular mechanism in this area should be further explored.
Conclusion
This MR study provides genetic evidence in support of the causal relationship between T1DM and an increased risk of IgAN, which may be driven by LDL-c. These findings shed new light on the function of T1DM in IgAN and may have implications for healthcare professionals seeking to improve lipid-lowering strategies in patients with T1DM in order to prevent IgAN.
Data availability statement
All the datasets were derived from sources in the public domain: GWAS catalog and MR-Base (https://www.mrbase.org/).
Author contributions
All authors participated in the field survey and data collection. PD and ZL drafted the manuscript. PD and ZL analyzed the data. BY and YL designed the study. BY, ZL, and YL obtained the funding. All authors participated in the field survey and data collection, critically revised the manuscript, and gave final approval to the version submitted for publication.
Funding
This work was supported by the National Science Foundation of China (NSFC) Projects (82170475, 82173905, 82070759), Hunan Provincial Natural Science Foundation of China (2021JJ31009, 2021JJ31032), Guiding Science and Technology Program Foundation of Changsha City (kzd21076), College Students’ Innovation and Entrepreneurship Project (s2021105330306, 2022105330116).
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/fendo.2022.1000627/full#supplementary-material
References
1. Bermejo S, Pascual J, Soler MJ. The current role of renal biopsy in diabetic patients. Minerva Med (2018) 109:116–25. doi: 10.23736/S0026-4806.17.05446-5
2. Donadio JV, Grande JP. IgA nephropathy. N Engl J Med (2002) 347:738–48. doi: 10.1056/NEJMra020109
3. Lai KN, Tang SC, Schena FP, Novak J, Tomino Y, Fogo AB, et al. IgA nephropathy. Nat Rev Dis Primers (2016) 2:16001. doi: 10.1038/nrdp.2016.1
4. Liu XM, Wang Q, Dong ZY, Zhang WG, Cai GY, Zhang L, et al. Clinicopathological features of nondiabetic renal diseases from different age groups: An observational cross-sectional study. Chin Med J (Engl) (2018) 131:2953–9. doi: 10.4103/0366-6999.247197
5. Schena FP, Nistor I. Epidemiology of IgA nephropathy: A global perspective. Semin Nephrol (2018) 38:435–42. doi: 10.1016/j.semnephrol.2018.05.013
6. Woo KT, Chan CM, Mooi CY, H. LC, Tan HK, Foo M, et al. The changing pattern of primary glomerulonephritis in Singapore and other countries over the past 3 decades. Clin Nephrol (2010) 74:372–83. doi: 10.5414/CNP74372
7. Moranne O, Watier L, Rossert J, Stengel B. Primary glomerulonephritis: an update on renal survival and determinants of progression. Qjm (2008) 101:215–24. doi: 10.1093/qjmed/hcm142
8. Braun N, Schweisfurth A, Lohöfener C, Lange C, Gründemann C, Kundt G, et al. Epidemiology of glomerulonephritis in northern Germany. Int Urol Nephrol (2011) 43:1117–26. doi: 10.1007/s11255-011-9955-4
9. McQuarrie EP, Mackinnon B, McNeice V, Fox JG, Geddes CC. The incidence of biopsy-proven IgA nephropathy is associated with multiple socioeconomic deprivation. Kidney Int (2014) 85:198–203. doi: 10.1038/ki.2013.329
10. Maixnerova D, Jancova E, Skibova J, Rysava R, Rychlik I, Viklicky O, et al. Nationwide biopsy survey of renal diseases in the Czech republic during the years 1994-2011. J Nephrol (2015) 28:39–49. doi: 10.1007/s40620-014-0090-z
11. Fiorentino M, Bolignano D, Tesar V, Pisano A, Biesen WV, Tripepi G, et al. Renal biopsy in patients with diabetes: a pooled meta-analysis of 48 studies. Nephrol Dial Transplant (2017) 32:97–110. doi: 10.1093/ndt/gfw070
12. Smith GD, Ebrahim S. Mendelian randomization': can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol (2003) 32:1–22. doi: 10.1093/ije/dyg070
13. Davies NM, Holmes MV, Davey Smith G. Reading mendelian randomisation studies: a guide, glossary, and checklist for clinicians. Bmj (2018) 362:k601. doi: 10.1136/bmj.k601
14. Chiou J, Geusz RJ, Okino ML, Han JY, Miller M, Melton R, et al. Interpreting type 1 diabetes risk with genetics and single-cell epigenomics. Nature (2021) 594:398–402. doi: 10.1038/s41586-021-03552-w
15. Taliun D, Harris DN, Kessler MD, Carlson J, Szpiech ZA, Torres R, et al. Sequencing of 53,831 diverse genomes from the NHLBI TOPMed program. Nature (2021) 590:290–9. doi: 10.1038/s41586-021-03205-y
16. Burgess S, Thompson SG. Avoiding bias from weak instruments in mendelian randomization studies. Int J Epidemiol (2011) 40:755–64. doi: 10.1093/ije/dyr036
17. Sakaue S, Kanai M, Tanigawa Y, Karjalainen J, Kurki M, Koshiba S, et al. A cross-population atlas of genetic associations for 220 human phenotypes. Nat Genet (2021) 53:1415–24. doi: 10.1038/s41588-021-00931-x
18. Emdin CA, Khera AV, Kathiresan S. Mendelian randomization. Jama (2017) 318:1925–6. doi: 10.1001/jama.2017.17219
19. Feehally J, Farrall M, Boland A, Gale DP, Gut I, Heath S, et al. HLA has strongest association with IgA nephropathy in genome-wide analysis. J Am Soc Nephrol (2010) 21:1791–7. doi: 10.1681/ASN.2010010076
20. Skrivankova VW, Richmond RC, Woolf BAR, Davies NM, Swanson SA, VanderWeele TJ, et al. Strengthening the reporting of observational studies in epidemiology using mendelian randomisation (STROBE-MR): explanation and elaboration. Bmj (2021) 375:n2233. doi: 10.1136/bmj.n2233
21. Burgess S, Butterworth A, Thompson SG. Mendelian randomization analysis with multiple genetic variants using summarized data. Genet Epidemiol (2013) 37:658–65. doi: 10.1002/gepi.21758
22. Bowden J, Davey Smith G, Haycock PC, Burgess S. Consistent estimation in mendelian randomization with some invalid instruments using a weighted median estimator. Genet Epidemiol (2016) 40:304–14. doi: 10.1002/gepi.21965
23. Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through egger regression. Int J Epidemiol (2015) 44:512–25. doi: 10.1093/ije/dyv080
24. Verbanck M, Chen CY, Neale B, Do R. Detection of widespread horizontal pleiotropy in causal relationships inferred from mendelian randomization between complex traits and diseases. Nat Genet (2018) 50:693–8. doi: 10.1038/s41588-018-0099-7
25. Greco MF, Minelli C, Sheehan NA, Thompson JR. Detecting pleiotropy in mendelian randomisation studies with summary data and a continuous outcome. Stat Med (2015) 34:2926–40. doi: 10.1002/sim.6522
26. Burgess S, Bowden J, Fall T, Ingelsson E, Thompson SG. Sensitivity analyses for robust causal inference from mendelian randomization analyses with multiple genetic variants. Epidemiology (2017) 28:30–42. doi: 10.1097/EDE.0000000000000559
27. Gorodezky C, Alaez C, Murguía A, Rodríguez A, Balladares S, Vazquez M, et al. HLA and autoimmune diseases: Type 1 diabetes (T1D) as an example. Autoimmun Rev (2006) 5:187–94. doi: 10.1016/j.autrev.2005.06.002
28. Vergès B. Dyslipidemia in type 1 diabetes: AMaskedDanger. Trends Endocrinol Metab (2020) 31:422–34. doi: 10.1016/j.tem.2020.01.015
29. Vergès B. Lipid disorders in type 1 diabetes. Diabetes Metab (2009) 35:353–60. doi: 10.1016/j.diabet.2009.04.004
30. Perkins BA, Sherr JL, Mathieu C. Type 1 diabetes glycemic management: Insulin therapy, glucose monitoring, and automation. Science (2021) 373:522–7. doi: 10.1126/science.abg4502
31. Choi WJ, Hong YA, Min JW, Koh ES, Kim HD, Ban TH, et al. Hypertriglyceridemia is associated with more severe histological glomerulosclerosis in IgA nephropathy. J Clin Med (2021) 10:4236. doi: 10.3390/jcm10184236
32. Mäkinen VP, Soininen P, Kangas AJ, Forsblom C, Tolonen N, Thorn LM, et al. Triglyceride-cholesterol imbalance across lipoprotein subclasses predicts diabetic kidney disease and mortality in type 1 diabetes: the FinnDiane study. J Intern Med (2013) 273:383–95. doi: 10.1111/joim.12026
33. Jia X, Pan X, Xie J, Shen P, Wang Z, Li Y, et al. The importance of sensitive screening for abnormal glucose metabolism in patients with IgA nephropathy. Clin Nephrol (2016) 85:30–7. doi: 10.5414/CN108693
34. Yavorska OO, Burgess S. MendelianRandomization: an r package for performing mendelian randomization analyses using summarized data. Int J Epidemiol (2017) 46:1734–9. doi: 10.1093/ije/dyx034
35. Sanderson E, Davey Smith G, Windmeijer F, Bowden J. An examination of multivariable mendelian randomization in the single-sample and two-sample summary data settings. Int J Epidemiol (2019) 48:713–27. doi: 10.1093/ije/dyy262
36. National Kidney Foundation. KDOQI Clinical Practice Guideline for Diabetes and CKD: 2012 update. Am J Kidney Dis (2012) 60:850–86. doi: 10.1053/j.ajkd.2012.07.005
37. Lee HJ, Mariappan MM, Norton L, Bakewell T, Feliers D, Oh SB, et al. Proximal tubular epithelial insulin receptor mediates high-fat diet-induced kidney injury. JCI Insight (2021) 6:e143619. doi: 10.1172/jci.insight.143619
38. Novak J, Julian BA, Mestecky J, Renfrow MB. Glycosylation of IgA1 and pathogenesis of IgA nephropathy. Semin Immunopathol (2012) 34:365–82. doi: 10.1007/s00281-012-0306-z
39. Rodrigues JC, Haas M, Reich HN. IgA nephropathy. Clin J Am Soc Nephrol (2017) 12:677–86. doi: 10.2215/CJN.07420716
40. Gans RO, Ueda Y, Ito S, Kohli R, Min I, Shafi M, et al. The occurrence of IgA-nephropathy in patients with diabetes mellitus may not be coincidental: a report of five cases. Am J Kidney Dis (1992) 20:255–60. doi: 10.1016/S0272-6386(12)80698-6
41. Murray JA. Celiac disease in patients with an affected member, type 1 diabetes, iron-deficiency, or osteoporosis? Gastroenterology (2005) 128:S52–56. doi: 10.1053/j.gastro.2005.02.029
42. Taplin CE, Barker JM. Autoantibodies in type 1 diabetes. Autoimmunity (2008) 41:11–8. doi: 10.1080/08916930701619169
43. Kiryluk K, Li Y, Sanna-Cherchi S, Rohanizadegan M, Suzuki H, Eitner F, et al. Geographic differences in genetic susceptibility to IgA nephropathy: GWAS replication study and geospatial risk analysis. PloS Genet (2012) 8:e1002765. doi: 10.1371/journal.pgen.1002765
44. Abed E, LaBarbera B, Dvorak J, Zhang Y, Beck J, Talsania M. Prevalence of dyslipidemia and factors affecting dyslipidemia in young adults with type 1 diabetes: evaluation of statin prescribing. J Pediatr Endocrinol Metab (2019) 32:327–34. doi: 10.1515/jpem-2018-0383
45. Osanami A, Tanaka M, Furuhashi M, Ohnishi H, Hanawa N, Yamashita T, et al. Increased LDL-cholesterol level is associated with deterioration of renal function in males. Clin Kidney J (2022) 15:1888–95. doi: 10.1093/ckj/sfac111
46. Lee HS, Song CY. Oxidized low-density lipoprotein and oxidative stress in the development of glomerulosclerosis. Am J Nephrol (2009) 29:62–70. doi: 10.1159/000151277
47. Singh J, Jain A, Bhamra R, Rathi V, Dhingra AK. The mechanistic role of different mediators in the pathophysiology of nephropathy: A review. Curr Drug Targets (2022). doi: 10.2174/1389450124666221026152647
Keywords: Mendelian randomization, causal relationship, non-diabetic renal disease, type 1 diabetes mellitus, immunoglobulin a nephropathy
Citation: Deng P, Li Z, Yi B and Leng Y (2022) A Mendelian randomization study to assess the genetic liability of type 1 diabetes mellitus for IgA nephropathy. Front. Endocrinol. 13:1000627. doi: 10.3389/fendo.2022.1000627
Received: 22 July 2022; Accepted: 24 November 2022;
Published: 14 December 2022.
Edited by:
Lingling Xu, Peking Union Medical College Hospital (CAMS), ChinaReviewed by:
Yang Wang, The First Affiliated Hospital of Xi’an Jiaotong University, ChinaDespoina Manousaki, CHU Sainte-Justine, Canada
Copyright © 2022 Deng, Li, Yi and Leng. 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: Bin Yi, eWliaW5feWJAMTYzLmNvbQ==; Yiping Leng, bHlwMDYyNkBhbGl5dW4uY29t
†These authors have contributed equally to this work and share first authorship