Skip to main content

ORIGINAL RESEARCH article

Front. Immunol., 15 March 2024
Sec. Molecular Innate Immunity

Heterozygous mutations in factor H aggravate pathological damage in a stable IgA deposition model induced by Lactobacillus casei cell wall extract

Jingyi Li,,,,&#x;Jingyi Li1,2,3,4,5†Yaping Dong,,,,&#x;Yaping Dong1,2,3,4,5†Feifei Chen,,,,&#x;Feifei Chen1,2,3,4,5†Hongyu Yang,,,,Hongyu Yang1,2,3,4,5Pei Chen,,,,Pei Chen1,2,3,4,5Hongyu Li,,,,Hongyu Li1,2,3,4,5Sufang Shi,,,,Sufang Shi1,2,3,4,5Xujie Zhou,,,,Xujie Zhou1,2,3,4,5Li Zhu,,,,Li Zhu1,2,3,4,5Yuemiao Zhang,,,,Yuemiao Zhang1,2,3,4,5Lijun Liu,,,,Lijun Liu1,2,3,4,5Xinfang XieXinfang Xie6Feng Yu*Feng Yu7*Jing JinJing Jin8Jicheng Lv,,,,*Jicheng Lv1,2,3,4,5*Hong Zhang,,,,Hong Zhang1,2,3,4,5
  • 1Renal Division, Peking University First Hospital, Beijing, China
  • 2Institute of Nephrology, Peking University, Beijing, China
  • 3Key Laboratory of Renal Disease, Ministry of Health of China, Beijing, China
  • 4Key Laboratory of Chronic Kidney Disease Prevention and Treatment (Peking University), Ministry of Education, Beijing, China
  • 5State Key Laboratory of Vascular Homeostasis and Remodeling, Peking University, Beijing, China
  • 6Department of Nephrology, The First Affiliated Hospital of Xi’an Jiaotong University, Xi'an, China
  • 7Department of Nephrology, Peking University International Hospital, Beijing, China
  • 8Northwestern University Feinberg School of Medicine, Division of Nephrology, Chicago, IL, United States

Introduction: Activation of complement through the alternative pathway (AP) has a key role in the pathogenesis of IgA nephropathy (IgAN). We previously showed, by intraperitoneal injection of Lactobacillus casei cell wall extract (LCWE), C57BL/6 mice develop mild kidney damage in association with glomerular IgA deposition. To further address complement activity in causing glomerular histological alterations as suggested in the pathogenesis of IgAN, here we used mice with factor H mutation (FHW/R) to render AP overactivation in conjunction with LCWE injection to stimulate intestinal production of IgA.

Methods: Dose response to LCWE were examined between two groups of FHW/R mice. Wild type (FHW/W) mice stimulated with LCWE were used as model control.

Results: The FHW/R mice primed with high dose LCWE showed elevated IgA and IgA-IgG complex levels in serum. In addition to 100% positive rate of IgA and C3, they display elevated biomarkers of kidney dysfunction, coincided with severe pathological lesions, resembling those of IgAN. As compared to wild type controls stimulated by the same high dose LCWE, these FHW/R mice exhibited stronger complement activation in the kidney and in circulation.

Discussion: The new mouse model shares many disease features with IgAN. The severity of glomerular lesions and the decline of kidney functions are further aggravated through complement overactivation. The model may be a useful tool for preclinical evaluation of treatment response to complement-inhibitors.

Introduction

IgA nephropathy (IgAN) is the most common form of primary glomerular disease worldwide (1). More than half of the patients will eventually progress to end-stage renal disease (2). The hallmark of IgAN is IgA deposition in the mesangial area, frequently accompanied by the co-deposition of complement C3. The pathophysiology of the disease is summarized in a four-hit model. Central to the pathogenesis is an increase in circulating galactose deficient-IgA1 [Gd-IgA1] (hit 1). These abnormally glycosylated IgA1 molecules form immune complexes with autoantibodies in the systemic circulation (hit 2). These immune complexes eventually deposit in the mesangial area of the glomerulus (hit 3) to trigger inflammatory responses, including complement activation leading to tissue damage (hit 4) (3). However, the existing animal models only reflect the disease process to a limited extent although they had made a great contribution to the research of IgA nephropathy (46). We previously reported a mouse model with intraperitoneal injection of Lactobacillus casei cell wall extract (LCWE) emulsified with Complete Freund’s Adjuvant (CFA) to cause chronic inflammation. In this model, we noted a continuous and stable deposition of IgA in the glomerular mesangial areas, with high circulating levels of IgA and IgA-IgG complexes (7). While this model showed consistent IgA deposition, complement C3 co-deposition is rare in this model and there is only mild kidney damage.

Large international genome–wide association studies (GWAS) identified variants of complement factor H-related (CFHR) genes being associated with IgAN. The genome–wide significant effect of CFHR3,1 gene deletion to reduce the risk for IgAN qualifies activators and regulators of the alternative pathway (AP) as major players in the pathogenesis of the disease (8). However, the precise role of factor H and CFHRs in IgAN progression remains unclear. As an important regulator of the AP, factor H prevents the excessive activation of complement (9). A point mutation (W1206R) of factor H results in an increased binding to C3b, impairing its interaction with host cells (10). Mice harboring this mutation have a higher propensity of localized complement activation to cause complement-mediated kidney injury. Homozygous mice carrying this mutation develop severe thrombotic microangiopathy (TMA). Although heterozygous animals stay largely normal, they are susceptible to complement activation in the kidney.

In the current study, we exploited the factor H heterozygous mice (FHW/R) by challenging them with LCWE. LCWE primes the animals for intestinal IgA production, in causing IgA deposition in the kidney. When combined with FHW/R that predispose the animals for complement overactivation, the model developed diverse histological lesions in the kidney that closely resemble IgAN’s MEST-C alterations.

Materials and methods

Mouse procedures

All mice were raised and maintained under specific pathogen-free (SPF) conditions. All experiments were performed in accordance with local guidelines for laboratory animal care and the study were approved and supervised by the Laboratory Animal Care and Use Committee of Peking University First Hospital (No. 202103).

Complement factor H heterozygous mutant mice (FHW/R) and wild type mice (FHW/W) were generated by mating FHW/R mice with C57BL/6J mice (FHW/W) (10). In this study, only male mice were included.

The extraction of Lactobacillus casei cell wall extract (LCWE) and its administration were prepared as previously described (7). The FHW/R mice were randomly allocated to three groups: PBS group (FHW/R-Buffer), low dose LCWE emulsified with CFA group (0.1ug/g, FHW/R-LCWELow), and high dose LCWE emulsified with CFA group (0.4ug/g, FHW/R-LCWEHigh). FHW/W mice were randomly assigned to one of the two groups: PBS group (FHW/W-Buffer), and high dose LCWE emulsified with CFA group (0.4ug/g, FHW/W-LCWEHigh). The FHW/W-Buffer is only used as a control for western blot.

Intraperitoneal injection (i.p.) of LCWE commenced when the mice reached 8 weeks of age. Collections of serum and urine samples started at 2 months of age. Blood samples were obtained from the orbital vein, while urine samples were collected using metabolic cages. In order to observe the sequential changes of pathological damage in the model, mice in each group were sacrificed at 18, 26, and 30 weeks of age.

Histology, immunofluorescence, immunohistochemistry and electron microscopy

Paraffin-embedded kidney specimens 1.5µm in thickness were stained with periodic acid-Schiff (PAS) using standard protocols for morphological analysis. Pathological lesions were evaluated using PAS staining. As there is no established universal standard for renal lesions in IgA nephropathy model mice, we assessed the lesions (MEST-C) according to the Oxford classification of IgA nephropathy (Supplementary Table S1). In addition, we evaluated arteriolar lesions (A) including wall thickening and onion-like lesions.

For immunofluorescence staining, 2-µm thick cryostat frozen kidney sections were incubated for 1 hour at 37°C with Alexa Fluor 488-conjugated goat antibodies specific for mouse IgA (SouthernBiotech, 1040-30, 1:200, 5ug/ml), FITC-conjugated goat anti-mouse C3 antibody (MP Biomedicals, 0855500, 1:1000, 4ug/ml), FITC-conjugated goat anti-mouse C1q antibody (HycultBiotech, HM-1096F, 1:500, 2ug/ml), rabbit anti-mouse Bb antibody (Genetex, GTX81508, 1:200, 5ug/ml) and rabbit anti-mouse MBL antibody (abcam, ab190834, 1:100, 10ug/ml). For the staining of the latter three antibodies, after washing three times with PBS, incubate the matching fluorescent secondary antibodies at 37°C for 1 hour. Nuclei were stained with DAPI (Abcam). The images were collected using a fluorescent microscope (Leica, Germany). The intensity of staining was evaluated and scored on a scale of - to +++: -, no staining and trace staining; +, mild staining; ++, moderate staining; and +++, high staining on a high-power field. For IgA and C3 colocalization staining, we used the another C3 antibody (abcam, ab200999, 1:100, 10ug/ml).

Immunohistochemistry staining was performed for CFH and C5b-9. First, 2-µm thick paraffin embedded renal tissues were deparaffinized in xylene and rehydrated in grading alcohols, then sections were treated by appropriate antigen retrieval methods for 10 minutes at 37°C, including proteinase k (ZSGB-BIO, ZLI-9016, 0.5mg/ml) and proteinase XXIV (Sigma, P0652, 0.5mg/ml). After quenching endogenous peroxidase activity for 10 minutes at room temperature, sections were incubated with primary antibody for 1 hour at 37°C. All the antibodies’ optimal dilutions were predetermined. For CFH (Merck, 341276, 1:5000), after incubation with fluorescent secondary antibodies, nuclei were stained with DAPI (Abcam). And for C5b-9 (abcam, ab55811, 1:250, 20ug/ml), after incubating with the HRP-conjugated secondary antibody, the sections were developed using a fresh hydrogen peroxide plus 3-3-diaminobenzidine tetrahydrochloride solution for an appropriate duration. Finally, the sections were counterstained with hematoxylin, anti-blue with lithium carbonate, dehydrated and cleared in alcohols and xylene. The examination of the sections was conducted using light microscopy. As negative controls, normal rabbit IgG or mouse IgG was used to replace the primary antibodies.

For transmission electron microscopy, kidneys were fixed with 2.5% glutaraldehyde in PBS, postfixed with 1% osmium tetroxide in PBS, dehydrated in a graduated series of ethanol dilutions, and embedded in Epon 812 resin. Then blocks were cut using an ultramicrotome (Ultracut; Leica). Ultrathin sections placed on 200-mesh copper were stained with 1% uranyl acetate in 50% ethanol and Reynolds lead citrate. Finally, the prepared sections were examined using a transmission electron microscope (JEM-1230; JEOL, Tokyo, Japan). The assessment of the EM images included the presence and location of glomerular electron-dense deposits(mainly), qualitative assessment of mesangial hypercellularity, endocapillary hypercellularity, podocyte foot process effacement, thickening of the glomerular basement membrane, widen of subendothelial loose layer.

Detection of immunoglobulins and immune complex levels by enzyme-linked immunosorbent assay

Serum levels of IgA levels were measured through sandwich ELISA method as described previously (7). In brief, plates were coated overnight with 2.5 μg/ml goat F(ab’)2 anti-mouse Ig (SouthernBiotech, 1012-01, 1:400) in sodium carbonate buffer (pH 9.6). After washing, the plates were blocked with 1% BSA for 2 hours. Then diluted serum samples and mouse IgA (SouthernBiotech, 0106-01, the first standard concentration is 100ng/ml), used as the standard, were added and incubated for 1 hour. Following another round of washing, the plates were incubated with horseradish peroxidase (HRP)-conjugated goat anti-mouse IgA (SouthernBiotech, 1040-05, 1:5000, 2ug/ml) antibodies for 1 hour. To determine serum levels of IgA-IgG complexes, a cross-capture ELISA was performed. The 96-well plates were coated with 2.5 μg/ml of goat anti-mouse IgA (SouthernBiotech, 1040-01, 1:400) overnight. After washing and blocking, the diluted serum samples were added and allowed to bind for 1 hour. Then the HRP-conjugated goat anti-mouse IgG antibody (SouthernBiotech, 1030-05, 1:5000, 0.2ug/ml) was applied. All reactions were carried out at room temperature. Finally, the plates were developed using 3,3’,5,5’-tetramethylbenzidine (TMB) and the reactions were stopped with 1 M sulfuric acid. The results were analyzed using an ELISA reader (Bio-Rad 550) at 450 and 570 nm wavelengths.

Biochemical analyses and measurement of blood pressure

Plates were coated with 2.5 μg/ml of goat anti-mouse albumin (BETHYL, A90-134A, 1:400) antibody to detect urine albumin levels. After sample incubating and washing, the HRP-conjugated goat anti-mouse albumin (BETHYL, A90-134P, 1:10000, 0.1ug/ml) was added at 1:10000 dilution for 1 hour. The reaction was developed by adding TMB and stopped with 1 M sulfuric acid. Urinary creatinine levels were determined using a creatinine assay kit (BioAssay Systems, DICT-500). Proteinuria was defined as the urine albumin/creatinine ratio (ACR, mg/g). Hematuria was detected by urine test strips (URIT 1V), which is then scored semi-quantitatively (Supplementary Figure S1). Serum creatinine levels were measured using a serum creatinine assay kit (Nanjing Jiancheng Bioengineering Institute, C011-2-1). Blood urea nitrogen levels were assessed using a chemistry kit (BioAssay Systems, DIUR-100). All procedures are the same as the kit instructions. Blood pressure was measured by a noninvasive tail-cuff method (Softron Biotechnology, BP98AWU). One mouse’s blood pressure was measured three times, and the average value was finally calculated.

Western blotting

Equal amounts of diluted serum were solubilized in SDS sample buffer under nonreducing conditions and subjected to electrophoresis in 4-12% gradient gel to separated monomeric IgA (mIgA) and polymeric IgA (pIgA). The proteins were then transferred to polyvinylidene fluoride membranes (PVDF, Millipore) and subjected to Western blot analysis using HRP-conjugated goat anti–mouse IgA (SouthernBiotech, 1040-05, 1:5000, 0.2ug/ml) antibody. The same method was used for serum C3 detection with HRP-conjugated goat antibodies specific for mouse C3 (MP Biomedicals, 0855557, 1:10000) under reducing conditions. The detection of serum Ba were performed under reducing conditions. And complement factor H (CFH) were detected under nonreducing conditions. Goat anti-mouse FB (Sigma) or CFH (Merck, 341276, 1:5000) antibodies were first incubated on the PVDF membranes followed by HRP-conjugated donkey anti-goat IgG antibody (Abcam, ab97110, 1:10000, 0.1ug/ml). Membranes were developed by enhanced chemical luminescence treatment (GE Healthcare). The grayscale value of Western blotting graphics was quantified by Image J software. Levels of target fragment were reported as fold change relative to control mice.

Statistics

Normality test is used to determine the distribution of sample data. Continuous variables were presented as the mean ± SD (normally distributed variables) or the median with 25th and 75th centiles (non-normally distributed variables) for data distribution. In this study, we conducted comparisons among three groups - HW/R-Buffer, HW/R-LCWELow, and HW/R-LCWEHigh- under the same genetic background. Furthermore, within the context of the same LCWE dose, we examined two groups: HW/R-LCWEHigh and HW/W-LCWEHigh. For comparisons involving two groups of continuous variables, T-tests were employed for normally distributed data, and Mann-Whitney U tests were utilized for non-normally distributed data. When dealing with ordinal categorical variables between two groups, Mann-Whitney U tests were applied. For comparisons related to positive rates in two groups, the chi-square test (specifically Fisher’s exact test) was employed. In cases where three groups were compared for continuous variables, One-way ANOVA was used for normally distributed data, while Kruskal–Wallis test were employed for non-normally distributed data. Ordinal categorical variables, such as IgA and complement C3 deposition intensity, were compared using Kruskal–Wallis test, and positive rates were compared using the chi-square test (Fisher’s exact test). Bonferroni correction was used when three groups are compared.

Results

Levels of serum IgA and IgA–IgG immune complexes increase following LCWE injection

Compared with the PBS group, the body weight of mice induced by LCWE decreased in the first week, but then gradually increased (Supplementary Figure S2). IgA levels and IgA-IgG complex levels in mouse sera were measured for up to 30 weeks. As shown in Figures 1A, B, both IgA and IgA-IgG levels were elevated in the LCWE groups at 14 weeks of age, two weeks after the completion of LCWE injections at 12 weeks. Over the next 16 weeks, IgA and IgA-IgG levels stayed significantly higher in all dose groups of LCWE as compared to those in the buffer group (P < 0.05). Among FHW/R mice, there was a clear correlation between the LCWE doses and the levels of circulating IgA or IgA-IgG after 14 weeks of age (P < 0.05). Meanwhile, there was no difference in IgA levels and IgA-IgG levels between FHW/R and FHW/W in the high LCWE dose group.

Figure 1
www.frontiersin.org

Figure 1 Serum levels of IgA and IgA–IgG complexes. (A, B) Quantification of IgA (A) and IgA–IgG complexes (B) at different time points by ELISA in FHW/R-Buffer, FHW/R-LCWELow, FHW/R-LCWEHigh, FHW/W-LCWEHigh. Data are from all mice raised to the 30 weeks old. N=8 mice in each group. (C, D) Detection of IgA by western blot in serum. Randomly selected samples of each group. N = 3-4 mice per group. *P < 0.05, **P < 0.01, ***P < 0.001. ns, not significant.

In addition, we measured the level of high molecular weight polyIgA (pIgA) (Figures 1C, D and Supplementary Figure S3). Both pIgA and mIgA levels in the FHW/R-LCWEHigh group were higher than those in the FHW/R-LCWELow group (P < 0.05). This is consistent with polymeric to monomeric IgA ratios that appear also higher in the FHW/R-LCWEHigh group at 18 weeks, albeit this has not reached statistical significance (Supplementary Figure S3A, 2C). When FHW/R- LCWEHigh and FHW/W- LCWEHigh were compared, there was no significant difference in pIgA levels at 18 weeks and 30 weeks (Figure 1D). We observed a slightly higher level of mIgA in the FHW/R-LCWEHigh group, with low poly-to-mono-IgA ratios at 18 weeks (Supplementary Figure S3B).

The FHW/R- LCWEHigh group exhibits strong mesangial IgA and C3 deposition

We detected IgA and C3 in the kidneys by immunofluorescence staining (Figures 2A–H). In order to observe the renal phenotype of mice sequentially, we sacrificed mice at 18 weeks, 26 weeks, and 30 weeks respectively (Figure 2I). All mice primed with LCWE had IgA signals in the glomerulus. By contrast, only 22.2% of the FHW/R-Buffer group (4/18) showed IgA positivity, all with relatively low signal intensity. And at 30 weeks, the positivity rate of the FHW/R-Buffer group is 25% (Figure 2J). Regarding C3, the FHW/R-LCWEHigh group had an overall positive rate of 90% (18/20), which reached 100% by 30 weeks of age (8/8). In the FHW/R-LCWELow group, C3 positivity was 70% (14/20), climbing to 75% (6/8) by 30 weeks. Importantly, the FHW/W-LCWEHigh group exhibited a relatively lower C3 positive rate (2/8) despite all having IgA deposition. These findings indicate that factor H mutation in FHW/R promotes C3 co-deposition with IgA in the glomeruli (Supplementary Figures S4A, B, C), consistent with the role for factor H that as a potent inhibitor of AP.

Figure 2
www.frontiersin.org

Figure 2 Immunofluorescence staining of IgA and C3 in kidneys. (A–D) Representative images of IgA in every groups. (E–H) Representative images of C3. (I) After LCWE induction, mice were sacrificed at different times. (J–K) Semi-quantitative scoring of IgA and C3 fluorescence intensity. Mann-Whitney U tests and Kruskal–Wallis test (Bonferroni correction) were used when comparing between two groups and between three groups respectively. Data are from all mice raised to the 30 weeks old. N = 8 mice per group. Scale bars: 50 μm. *P < 0.05, **P < 0.01, ***P < 0.001.

We also observed interesting differences in the deposition patterns of IgA and C3. In the FHW/R-LCWELow and FHW/W-LCWEHigh groups, the majority of mice displayed diffuse deposition of IgA and C3 in the mesangial area (Figures 2B, C). However, FHW/R-LCWEHigh mice displayed varying patterns of deposition, ranging from diffuse mesangial deposition, with or without additional deposits in the capillary loop (Supplementary Figure S4D). Furthermore, mice with more severe glomerular lesions tended to have granular deposits of IgA (Supplementary Figure S4G).

Factor H heterozygous mice have more severe lesions

As expected, FHW/R-Buffer mice all had normal kidney histology (Figures 3A–C), whereas the FHW/R-LCWELow group developed glomerular and vascular lesions (Figures 3E–G). The glomerular lesions predominantly manifested as mesangial hypercellularity (M1) based on Oxford classification of IgAN. And mesangial cell proliferation and electron-dense deposition can be seen under electron microscopy (Figure 3H). In the FHW/R-LCWEHigh group, in addition to M1, the majority of animals also presented with endocapillary hypercellularity (E1, 17/20, 85%) (Figure 3I), with some showing crescent formation (C1, 3/20, 15%) (Supplementary Figures S4E, H), and tubular atrophy/interstitial fibrosis (T1, 4/20, 20%) (Figure 3K). Consistent with clinical IgAN that frequently has arteriolar lesions, 35% (7/20) FHW/R-LCWEHigh mice developed arteriolar lesions, with 2 mice showing onion-like lesions (10%) (Detailed information is in Table 1) (Figure 3J). In addition, we found no mice had microthrombi in glomerular capillary lumen and renal arteriolar thrombosis (0/20). Overall, mesangial hypercellularity is the dominant type of lesion of the FHW/R-LCWELow group, whereas FHW/R-LCWEHigh animals exhibited a greater variety of lesion types, indicating severity of kidney damage. A comparison between FHW/R-LCWELow and FHW/R-LCWEHigh groups revealed pathological severity in correlation with LCWE dose.

Figure 3
www.frontiersin.org

Figure 3 Representative micrographs of PAS staining and electron micrographs from each group are shown. (A–D) Normal glomerulus, renal vascular and renal interstitium from FHW/R-Buffer group mice. (E–H) Mice in FHW/R-LCWELow group showed mesangial hypercellularity, and the arterioles from a mouse showed thickening of arteriolar walls and stenosis of lumens. (I–L) Mice from FHW/R-LCWEHigh group had diverse lesions, including mesangial hypercellularity and endothelial hypercellularity (I), arteriolar onion-like lesions (J, the red arrow), tubular atrophy/interstitial fibrosis (K). Mice exhibited mesangial cell proliferation and marked mesangial expansion. Electron-dense deposits were apparent in the mesangial area (L, the red arrow). (M-P) FHW/W-LCWEHigh group mice showed milder pathological lesions compared with FHW/R-LCWEHigh group mice. Lesions are mainly mild mesangial hypercellularity (M). In the PAS staining image, the magnification of the renal interstitium is 200X and the rest is 400X. Scale bars in PAS figures: 50 μm. Scale bars in electron micrographs: 2μm.

Table 1
www.frontiersin.org

Table 1 Pathological lesions of all mice from four groups.

Interestingly, FHW/W-LCWEHigh had only mild mesangial hypercellularity (M1,14/18), with 2 out of 18 mice having endocapillary hypercellularity (E1) and 1 mouse had arteriolar lesions (A1) (Figures 3M–P). In contrast, FHW/R-LCWEHigh animals exhibited a greater diversity and severity of glomerular lesions, suggesting complement activation contributes to glomerular damage.

Transmission electron microscopy (TEM) of FHW/R-Buffer kidneys showed normal basement membranes and podocyte foot processes. No evidence of mesangial cell proliferation or other glomerular injuries was observed (Figure 3D). In contrast, mice that had been primed with LCWE, particularly those in the FHW/R-LCWEHigh group, exhibited mesangial cell proliferation and marked mesangial expansion (Figure 3L). Electron-dense deposits were apparent in the mesangial area across all three dose groups, and leukocyte infiltration (neutrophils) was observed in the capillary loops of some LCWE-primed mice (Supplementary Figure S4F). In addition, foot process effacement was observed in FHW/R-LCWEHigh mice, indicating severity of the lesions (Supplementary Figure S4I).

FHW/R-LCWEHigh mice developed proteinuria and hematuria

We detected urinary albumin levels in mice at 18, 26, and 30 weeks of age (Figures 4A, B). Consistent with more severe pathological damage, the FHW/R-LCWEHigh group also had elevated levels of urinary albumin. At 18 weeks old of mice, two weeks after the completion of LCWE injections, urine albumin levels were higher in the FHW/R-LCWEHigh group than in the FHW/R-Buffer, in the FHW/R-LCWELow, and in the FHW/W-LCWEHigh group (23.11 ± 7.72 vs 9.39 ± 3.44 vs 14.8 ± 5.41 vs 14.01 ± 4 mg/g, P<0.05, Figure 4A). After correcting for urine creatinine, FHW/R-LCWEHigh group is still higher than other groups (P<0.05, Figure 4B).

Figure 4
www.frontiersin.org

Figure 4 Kidney function analysis of all the mice. (A–C) The albumin levels and hematuria degree from all mice at different time point. Each point represents 1 mouse. The A depicted the results of urinary albumin and the B showed the urinary albumin normalized to creatinine. T-tests and one-way ANOVA (Bonferroni correction) were used when comparing between two groups and between three groups respectively. C shows the results of hematuria. Mann-Whitney U tests and Kruskal–Wallis test (Bonferroni correction) were used when comparing between two groups and between three groups respectively.

We also measured hematuria using urine test strips. At 18 weeks and 30 weeks of age, approximately half of the mice in both FHW/R-LCWEHigh and FHW/R-LCWELow groups had hematuria (Figure 4C). Whereas none of the FHW/W-LCWEHigh mice exhibited hematuria, indicating the critical contribution of complement activation in the progression of the disease.

We tested serum creatinine and BUN levels in serum samples of mice raised to 30 weeks of age. As shown in Supplementary Figure S5, mice in the FHW/R-LCWEHigh group exhibited higher serum BUN levels compared to those in the FHW/R-Buffer and FHW/R-LCWELow group. Notably, no statistically significant difference was found when compared to mice in the FHW/W-LCWEHigh group (P=0.054). However, no statistically significant difference was observed in serum creatinine levels among the four groups.

Blood pressure was assessed at 30 weeks of age, and as illustrated in Supplementary Figure S6, no significant differences were observed in systolic blood pressure (108 ± 10 vs. 103 ± 9 vs. 94 ± 15 vs. 100 ± 15 mmHg), diastolic blood pressure (44 ± 4 vs. 46 ± 7 vs. 41 ± 13 vs. 54 ± 19 mmHg), and mean arterial pressure (65 ± 6 vs. 63 ± 8 vs. 58 ± 13 vs. 52 ± 8 mmHg) among the four groups of mice.

Glomerular complement activation in FHW/R-LCWEHigh mouse

Co-deposition of C3 with IgA was observed in the glomerulus (Supplementary Figure S4C). FHW/R-LCWEHigh mice had higher C3 positivity than their FHW/W-LCWEHigh counterparts (100% vs. 25%, Figure 2K) at 30 weeks, consistent with the expectation of the factor H heterozygous mice being more susceptible to complement activation. While factor H signals were weak in FHW/R-Buffer animals, the addition of LCWE priming significantly enhanced factor H deposition in the kidney. This appeared more prominently in the FHW/R-LCWEHigh group than in the FHW/W-LCWEHigh and the FHW/R-LCWELow groups (Figures 5A–D, I). These results indicate the dynamics of complement activation in response to IgA deposition being driven by inflammation, whereas factor H counteracts the overactivation of C3.

Figure 5
www.frontiersin.org

Figure 5 The deposition of complement C5b-9 and H in the kidney. (A–D) Representative micrographs of CFH staining from each group are shown. (E–H) Representative micrographs of C5b-9 staining from each group. Mice in FHW/R-LCWEHigh group showed more FH deposition and more C5b-9 activation (I, J) For H, Mann-Whitney U tests and Kruskal–Wallis test (Bonferroni correction) were used when comparing between two groups and between three groups respectively. For C5b-9, Fisher’s exact test was used. Scale bars: 50 μm. **P < 0.01, ***P < 0.001.

Furthermore, C5b-9 participates in the membrane attack complex (MAC) in the final step of complement activation. We performed immunohistochemistry (IHC) to detect C5b-9 in the kidney (Figures 5E–H). In the FHW/R-Buffer group, only 1 out of 18 mice had weak signals of C5b-9. Similarly, in the FHW/W-LCWEHigh group, only 2 out of 18 mice showed weak positivity. Meanwhile, C5b-9 positive rate increased to 25% (5/20) in the FHW/R-LCWELow group, illustrating the importance of FH in suppressing complement activity. Notably, C5b-9 positive rate in the FHW/R-LCWEHigh group reached 55% (11/20) (Figure 5J). FHW/R-LCWEHigh mice that exhibited S1, T1 or C1 lesions generally had a higher C5b-9 deposition rate (5/5) as compared to the mice with just M1 or M1/E1 lesions (6/15). This finding indicates that more severe lesions were associated with increased complement activation.

To better understand the complement pathways involved in the renal damage, we performed immunofluorescence staining of C1q (representing the classical pathway), Bb (representing the alternative pathway) and MBL (representing the lectin pathway) using frozen sections of kidney tissue (Supplementary Figure S7). The C1q positivity rate is low, only 2 and 1 mice were found to have weak positive deposition of C1q in HW/R-LCWEHigh group (2/20) and HW/W-LCWEHigh group (1/18) respectively. Supplementary Figure S7A shows that one mouse in HW/R-LCWEHigh group is weakly positive for C1q. Under the induction of LCWE, 2 mice in HW/R-LCWELow group (2/20), 4 mice in HW/R-LCWEHigh group (4/20) and 2 mice in HW/W-LCWEHigh group (2/18) showed MBL positive (Supplementary Figure S7B). The overall MBL positive rate was not high and there was no statistical difference between the groups (Supplementary Figure S7H). In contrast, the activation of the alternative pathway in HW/R-LCWEHigh group is stronger. 1 mouse in HW/R-Buffer group, 3 mice in HW/R-LCWELow group (3/20), 15 mice in HW/R-LCWEHigh group (15/20) and 3 mice in HW/W-LCWEHigh group (3/18) showed Bb positive (Supplementary Figures S7C, D). The Bb positive rate of HW/R-LCWEHigh group was significantly higher than that of the other three groups (Supplementary Figure S7I).

FHW/R-LCWEHigh mice exhibit complement activation in circulation

Next, we performed Western blotting to evaluate the status of complement activation in circulation by measuring the ratio between intact C3α-chain and the resulting C3α-chain fragments following complement activation. At 18 weeks (Figure 6A) and 30 weeks (Figure 6B), the FHW/R-LCWEHigh group showed higher levels of cleaved C3α fragments than the FHW/R-LCWELow group (P < 0.05, Supplementary Figures S8A, B). The difference in ratios had reached statistical differences at 30 weeks (Supplementary Figure S8F). And the FHW/R-LCWEHigh group showed higher C3α-chain fragments/intact C3α-chain ratio than FHW/W-LCWEHigh group at 18 weeks and 30 weeks (P < 0.05, Supplementary Figures S8G, H). Overall, the FHW/R-LCWEHigh group had the highest active fragments to intact C3α ratio, indicating complement activation.

Figure 6
www.frontiersin.org

Figure 6 Complement activation in circulation. (A, B): western blot of C3 at 18 weeks and 30 weeks respectively. (C, D): western blot of Ba at 18 weeks and 30 weeks respectively. (E, F): western blot of FH at 18 weeks and 30 weeks respectively.

In addition to C3, we performed Western blotting of factor H (FH) and factor B (FB) of the alternative pathway. FHW/R-LCWEHigh had higher circulating factor H levels than FHW/R-LCWELow (Figure 6 and Supplementary Figures S8I, J). However, there was no statistical difference in factor H levels between the FHW/R-LCWEHigh group and the FHW/W-LCWEHigh group at 18 weeks and 30 weeks.

With respect to factor B, we measured factor B cleavage fragment Ba. Our findings were consistent with those for the C3α-chain cleavage fragments. The levels of Ba were elevated in the FHW/R-LCWEHigh group as compared to the FHW/R-Buffer and FHW/R-LCWELow groups. Moreover, FHW/R-LCWEHigh animals exhibited higher Ba levels than their FHW/W-LCWEHigh counterparts, consistent with the notion that FH mutant mice are more predisposed to complement activation.

Discussion

We have previously developed a continuous and stable IgA deposition model induced by LCWE, which partially replicates the upstream mechanisms (hits 1, 2, and 3) in IgA nephropathy. This model exhibits increased intestinal permeability, elevated levels of IgA and IgA-IgG, as well as mesangial IgA deposition, accompanied by partial C3 deposition. However, the pathological manifestations of this model are relatively mild and limited, as these mice do not present proteinuria or hematuria. In the current study, we investigated factor H heterozygous mutant mice, which have a genetic basis that predisposes them to complement activation. With the addition of LCWE stimulation, these mice developed severe and diverse pathological damages in the glomerulus, closely resembling the characteristics of IgA nephropathy. There was no overmortality and mice at 30 weeks were morphologically not substantially different from control mice.

Complement activation through the alternative pathway plays a crucial role in the pathogenesis of IgA nephropathy (1, 11). Numerous studies conducted over the years have discovered elevated circulating levels and mesangial deposition of complement proteins of the alternative pathway (1223). Additionally, and most importantly, large genome-wide association studies (GWAS) of IgA nephropathy have consistently detected the susceptible locus 1q32, which contains the CFH and CFHRs genes (8) (24, 25). In CFH and CFHRs, studies had found that there are regulatory genetic variations affecting the expression levels of circulating complement proteins (26, 27), and there are also genetic variations in the coding region affecting protein structure (26, 28). All of these studies jointly revealed the important roles of CFH and CFHRs in regulating complement activation of the alternative pathway through multiple mechanisms and influence the severity of renal tissue damage. In this study, we developed mouse model based on factor H heterozygous mutant mice. By comparing and FHW/W-LCWEHigh mice under the same dose induction, we discovered that the FHW/R-LCWEHigh group developed more severe and diverse lesions and had higher levels of proteinuria, more pronounced complement activation in the kidney. This observation indicates alternative pathway activation contributes to the aggravation of pathological damage. Numerous studies emphasize the close association between factor H-related protein and IgA nephropathy. A mouse model based on humanized factor H-related protein may be better able to understand the role of alternative complement pathway in human IgAN.

Our model also corroborates our previous findings of both IgA levels and IgA-IgG complex levels being elevated after mucosal inflammation (7). And in this study, we conducted a comparative analysis between FHW/R-LCWEHigh group mice and FHW/R-LCWELow group mice, which revealed that high-dose LCWE had a more pronounced effect compared to low-dose LCWE. Specifically, FHW/R-LCWEHigh group mice exhibited higher IgA levels and IgA-IgG levels, as well as more severe pathological lesions in the kidneys. This was attributed to higher levels of IgA and IgA-IgG complexes on the one hand, and higher levels of mutant factor H on the other hand.

Intrarenal arterial lesions are common in patients with IgAN. Prior studies, including our own cohort study, consistently demonstrated that arteriolar microangiopathic lesion(MA) was an independent risk factor for kidney failure (29) and suggested a potential involvement of complement activation in the development of arteriole lesions, especially those with MA (15, 16, 29). In previously published studies, mice with homozygous mutations in factor H exhibited aHUS, a type of thrombotic microangiopathy (TMA) while mouse with heterozygous mutant didn’t show arteriolar lesions (10). In this study, we observed that 35% of the mice in the FHW/R-LCWEHigh group developed renal arteriolar lesions, with some exhibiting concomitant MA. In contrast, FHW/W-LCWEHigh mice did not have MA, suggesting the alternative complement activation after IgA kidney deposition contribute to the endothelial lesions.

Our study demonstrates that complement activation after IgA kidney deposition contribute to the renal phenotype in mice. Mutation in the complement system might render indolent IgA deposits pathologic. Our model is closer to the four-hit theory to a certain extent, it may help to advance our knowledge of IgAN and may be a useful tool for preclinical evaluation of treatment response to complement-inhibitors. Also it should be cautious that our present model cannot fully reflect and represent human IgAN for not all patients have complement gene mutations.

In summary, our study established a new animal model of IgAN, highlighting the significant role of complement activation in the pathogenesis of the disease. This model will serve as a valuable tool for advancing research on IgAN towards the understanding of complement-mediated renal damage. The new model also provides a useful tool for the development of therapies for IgAN.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The animal study was approved by Laboratory Animal Care and Use Committee of Peking University First Hospital. The study was conducted in accordance with the local legislation and institutional requirements.

Author contributions

JyL: Data curation, Formal analysis, Investigation, Methodology, Resources, Software, Validation, Visualization, Writing – original draft. YD: Data curation, Investigation, Methodology, Writing – original draft. FC: Data curation, Resources, Writing – original draft. HY: Funding acquisition, Supervision, Writing – review & editing. PC: Funding acquisition, Supervision, Validation, Writing – review & editing. HL: Investigation, Writing – review & editing. SS: Supervision, Writing – review & editing. XZ: Supervision, Writing – review & editing. LZ: Supervision, Writing – review & editing. YZ: Supervision, Writing – review & editing. LL: Supervision, Writing – review & editing. XX: Supervision, Writing – review & editing. FY: Resources, Supervision, Writing – review & editing. JJ: Project administration, Supervision, Writing – review & editing. JcL: Conceptualization, Funding acquisition, Project administration, Supervision, Writing – review & editing. HZ: Conceptualization, Funding acquisition, Project administration, Supervision, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by the Beijing Natural Science Foundation (7224346), National Natural Science Foundation of China (81925006, 82370712, 82100737, 82070733), and CAMS Innovation Fund for Medical Sciences (2019-12M-5-046). National High level hospital clinical research funding (high quality clinical research project of Peking University First Hospital No. 2022CR80, 2022CX15).

Acknowledgments

The authors are grateful to professor Wen-Chao Song from University of Pennsylvania for kindly providing the FH mutant mouse. And we would like to express our sincere appreciation to Jinwei Wang from Peking university First hospital for his guidance in statistical methods throughout this research.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2024.1368322/full#supplementary-material

References

1. Wyatt RJ, Julian BA. IgA nephropathy. New Engl J Med. (2013) 368:2402–14. doi: 10.1056/NEJMra1206793

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Pitcher D, Braddon F, Hendry B, Mercer A, Osmaston K, Saleem MA, et al. Long-term outcomes in igA nephropathy. Clin J Am Soc Nephrology: CJASN. (2023) 18:727–38. doi: 10.2215/CJN.0000000000000135

CrossRef Full Text | Google Scholar

3. Floege J, Barbour SJ, Cattran DC, Hogan JJ, Nachman PH, Tang SCW, et al. Management and treatment of glomerular diseases (part 1): conclusions from a Kidney Disease: Improving Global Outcomes (KDIGO) Controversies Conference. Kidney Int. (2019) 95:268–80. doi: 10.1016/j.kint.2018.10.018

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Okazaki K, Suzuki Y, Otsuji M, Suzuki H, Kihara M, Kajiyama T, et al. Development of a model of early-onset IgA nephropathy. J Am Soc Nephrology: JASN. (2012) 23:1364–74. doi: 10.1681/ASN.2011121160

CrossRef Full Text | Google Scholar

5. McCarthy DD, Kujawa J, Wilson C, Papandile A, Poreci U, Porfilio EA, et al. Mice overexpressing BAFF develop a commensal flora-dependent, IgA-associated nephropathy. J Clin Invest. (2011) 121:3991–4002. doi: 10.1172/JCI45563

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Berthelot L, Papista C, Maciel TT, Biarnes-Pelicot M, Tissandie E, Wang PH, et al. Transglutaminase is essential for IgA nephropathy development acting through IgA receptors. J Exp Med. (2012) 209:793–806. doi: 10.1084/jem.20112005

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Wan F, Wang H, Wang M, Lv J, Zhao M, Zhang H. Sustained release of Lactobacillus casei cell wall extract can induce a continuous and stable IgA deposition model. J Pathol. (2022) 257:262–73. doi: 10.1002/path.5884

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Xie J, Kiryluk K, Li Y, Mladkova N, Zhu L, Hou P, et al. Fine mapping implicates a deletion of CFHR1 and CFHR3 in protection from igA nephropathy in han Chinese. J Am Soc Nephrol. (2016) 27:3187–94. doi: 10.1681/ASN.2015111210

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Józsi M, Zipfel PF. Factor H family proteins and human diseases. Trends Immunol. (2008) 29:380–7. doi: 10.1016/j.it.2008.04.008

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Ueda Y, Mohammed I, Song D, Gullipalli D, Zhou L, Sato S, et al. Murine systemic thrombophilia and hemolytic uremic syndrome from a factor H point mutation. Blood. (2017) 129:1184–96. doi: 10.1182/blood-2016-07-728253

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Tortajada A, Gutierrez E, Pickering MC, Praga Terente M, Medjeral-Thomas N. The role of complement in IgA nephropathy. Mol Immunol. (2019) 114:123–32. doi: 10.1016/j.molimm.2019.07.017

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Miyazaki R, Kuroda M, Akiyama T, Otani I, Tofuku Y, Takeda R. Glomerular deposition and serum levels of complement control proteins in patients with IgA nephropathy. Clin Nephrol. (1984) 21:335–40.

PubMed Abstract | Google Scholar

13. Rauterberg EW, Lieberknecht HM, Wingen AM, Ritz E. Complement membrane attack (MAC) in idiopathic IgA-glomerulonephritis. Kidney Int. (1987) 31:820–9. doi: 10.1038/ki.1987.72

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Tomino Y, Sakai H, Nomoto Y, Endoh M, Arimori S, Fujita T. Deposition of C4-binding protein and beta 1H globulin in kidneys of patients with IgA nephropathy. Tokai J Exp Clin Med. (1981) 6:217–22.

PubMed Abstract | Google Scholar

15. Chua JS, Zandbergen M, Wolterbeek R, Baelde HJ, van Es LA, de Fijter JW, et al. Complement-mediated microangiopathy in IgA nephropathy and IgA vasculitis with nephritis. Modern pathology: an Off J United States Can Acad Pathology Inc. (2019) 32:1147–57. doi: 10.1038/s41379-019-0259-z

CrossRef Full Text | Google Scholar

16. Li J, Guo L, Shi S, Zhou X, Zhu L, Liu L, et al. The role of complement in microangiopathic lesions of igA nephropathy. Kidney Int Rep. (2022) 7:1219–28. doi: 10.1016/j.ekir.2022.03.028

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Medjeral-Thomas NR, Troldborg A, Constantinou N, Lomax-Browne HJ, Hansen AG, Willicombe M, et al. Progressive igA nephropathy is associated with low circulating mannan-binding lectin-associated serine protease-3 (MASP-3) and increased glomerular factor H-related protein-5 (FHR5) deposition. Kidney Int Rep. (2018) 3:426–38. doi: 10.1016/j.ekir.2017.11.015

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Wei M, Guo WY, Xu BY, Shi SF, Liu LJ, Zhou XJ, et al. Collectin11 and complement activation in igA nephropathy. Clin J Am Soc Nephrology: CJASN. (2021) 16:1840–50. doi: 10.2215/CJN.04300321

CrossRef Full Text | Google Scholar

19. Chiu YL, Lin WC, Shu KH, Fang YW, Chang FC, Chou YH, et al. Alternative complement pathway is activated and associated with galactose-deficient igA(1) antibody in igA nephropathy patients. Front Immunol. (2021) 12:638309. doi: 10.3389/fimmu.2021.638309

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Zwirner J, Burg M, Schulze M, Brunkhorst R, Gotze O, Koch KM, et al. Activated complement C3: a potentially novel predictor of progressive IgA nephropathy. Kidney Int. (1997) 51:1257–64. doi: 10.1038/ki.1997.171

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Tanaka C, Suhara Y, Kikkawa Y. [Circulating immune complexes and complement breakdown products in childhood IgA nephropathy]. Nihon Jinzo Gakkai shi. (1991) 33:709–17.

PubMed Abstract | Google Scholar

22. Zhu L, Guo WY, Shi SF, Liu LJ, Lv JC, Medjeral-Thomas NR, et al. Circulating complement factor H-related protein 5 levels contribute to development and progression of IgA nephropathy. Kidney Int. (2018) 94:150–8. doi: 10.1016/j.kint.2018.02.023

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Medjeral-Thomas NR, Lomax-Browne HJ, Beckwith H, Willicombe M, McLean AG, Brookes P, et al. Circulating complement factor H-related proteins 1 and 5 correlate with disease activity in IgA nephropathy. Kidney Int. (2017) 92:942–52. doi: 10.1016/j.kint.2017.03.043

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Gharavi AG, Kiryluk K, Choi M, Li Y, Hou P, Xie J, et al. Genome-wide association study identifies susceptibility loci for IgA nephropathy. Nat Genet. (2011) 43:321–7. doi: 10.1038/ng.787

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Kiryluk K, Li Y, Scolari F, Sanna-Cherchi S, Choi M, Verbitsky M, et al. Discovery of new risk loci for IgA nephropathy implicates genes involved in immunity against intestinal pathogens. Nat Genet. (2014) 46:1187–96. doi: 10.1038/ng.3118

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Zhu L, Zhai YL, Wang FM, Hou P, Lv JC, Xu DM, et al. Variants in complement factor H and complement factor H-related protein genes, CFHR3 and CFHR1, affect complement activation in igA nephropathy. J Am Soc Nephrology: JASN. (2015) 26:1195–204. doi: 10.1681/ASN.2014010096

CrossRef Full Text | Google Scholar

27. Ansari M, McKeigue PM, Skerka C, Hayward C, Rudan I, Vitart V, et al. Genetic influences on plasma CFH and CFHR1 concentrations and their role in susceptibility to age-related macular degeneration. Hum Mol Genet. (2013) 22:4857–69. doi: 10.1093/hmg/ddt336

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Zhai YL, Meng SJ, Zhu L, Shi SF, Wang SX, Liu LJ, et al. Rare variants in the complement factor H-related protein 5 gene contribute to genetic susceptibility to igA nephropathy. J Am Soc Nephrology: JASN. (2016) 27:2894–905. doi: 10.1681/ASN.2015010012

CrossRef Full Text | Google Scholar

29. Cai Q, Shi S, Wang S, Ren Y, Hou W, Liu L, et al. Microangiopathic lesions in igA nephropathy: A cohort study. Am J Kidney Dis. (2019) 74:629–39. doi: 10.1053/j.ajkd.2019.03.416

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: IgA nephropathy, complement factor H, mouse model, the alternative complement pathway, proteinuria

Citation: Li J, Dong Y, Chen F, Yang H, Chen P, Li H, Shi S, Zhou X, Zhu L, Zhang Y, Liu L, Xie X, Yu F, Jin J, Lv J and Zhang H (2024) Heterozygous mutations in factor H aggravate pathological damage in a stable IgA deposition model induced by Lactobacillus casei cell wall extract. Front. Immunol. 15:1368322. doi: 10.3389/fimmu.2024.1368322

Received: 10 January 2024; Accepted: 27 February 2024;
Published: 15 March 2024.

Edited by:

Christoph Daniel, University Hospital Erlangen, Germany

Reviewed by:

Kevin James Marchbank, Newcastle University, United Kingdom
Claudia Seikrit, RWTH Aachen University, Germany

Copyright © 2024 Li, Dong, Chen, Yang, Chen, Li, Shi, Zhou, Zhu, Zhang, Liu, Xie, Yu, Jin, Lv and Zhang. 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: Jicheng Lv, mailto:jichenglv75@gmail.com; Feng Yu, mailto:yufengevert1@sina.com

†These authors share first authorship

Disclaimer: 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.