Multiple Processes May Involve in the IgG4-RD Pathogenesis: An Integrative Study via Proteomic and Transcriptomic Analysis

Immunoglobulin G4-related disease (IgG4-RD) is a newly defined disease entity, while the exact pathogenesis is still not clear. Identifying the characters of IgG4-RD in proteomic and transcriptomic aspects will be critical to investigate the potential pathogenic mechanisms of IgG4-RD. We performed proteomic analysis realized with iTRAQ technique for serum samples from eight treatment-naive IgG4-RD patients and eight healthy volunteers, and tissue samples from two IgG4-RD patients and two non-IgG4-RD patients. Transcriptomic data (GSE40568 and GSE66465) was obtained from the GEO Dataset for validation. The weighted correlation network analysis (WGCNA) was applied to detect the gene modules correlated with IgG4-RD. KEGG pathway analysis was used to investigate pathways enriched in IgG4-RD samples. As a result, a total of 980 differentially expressed proteins (DEPs) in tissue and 94 DEPs in serum were identified between IgG4-RD and control groups. Three hundred fifty-four and two hundred forty-seven genes that most correlated with IgG4-RD were detected by WGCNA analysis in tissue and PBMC, respectively. We also found that DEPs in IgG4-RD samples were enriched in several immune-related activities including bacterial/viral infections and platelet activation as well as some immune related signaling pathways. In conclusion, we identified multiple processes/factors and several signaling pathways that may involve in the IgG4-RD pathogenesis, and found out some potential therapeutic targets for IgG4-RD.


INTRODUCTION
Immunoglobulin G4-related disease (IgG4-RD) is a newly defined immune-mediated disease with common clinical, serological, and pathological features (1). Common features of IgG4-RD include serum IgG4 level elevation, multiple organ involvements, dense infiltration of IgG4+ plasma cells, and significant tissue infiltrates (2). This disease affects men more often than women and age at diagnosis ranges from 50 to 70 years (3). Although most patients do respond to steroids well, the relapse rate can be nearly 50% (4).
IgG4-RD is also characterized by alterations in acquired immune system, in which aberrant expansion of plasmablasts, CD4+ cytotoxic T cells, and follicular T helper cells have been observed (5)(6)(7). In addition, several inflammatory factors such as TGF-β, IL-4, and IL-10 have also been identified to play a role in the pathogenesis of IgG4-RD (8). The increase of these cytokines promotes eosinophilia in serum or certain tissues, high levels of IgG4-producing plasma cells, elevated production of IgE, and fibrosis, with inflammatory cell infiltrates ultimately causing organ damage (9). Furthermore, several autoantibodies, including anti-carbonic anhydrase II and anti-lactoferrin, are often present in patients with IgG4-RD, especially those with IgG4-related autoimmune pancreatitis (AIP) (10). At present, however, the exact pathogenic mechanism remains unclear, which is an important issue in IgG4-RD studies.
Due to the lack of ideal animal models, and limited sample origin, high throughput, and bioinformatics techniques may help understand the underlying pathogenesis of IgG4-RD more deeply. Transcriptome-wide profiling, as a downstream level of genome-scale mapping, can reveal a systemic dynamics of molecular interaction (11). Recently, studies have utilized transcript profiling in labial salivary glands (LSGs) to distinguish molecular features between IgG4-RD and Sjögren's syndrome (SS), a disease with common phenotypic elements (12,13). Among other findings, active involvement of Th2-(IL-4, IL-5, and IL-21), T follicular helper cell (Tfh)-(BCL-6 and CXCR5) and Treg-related transcripts (IL-10, FOXP3, CCL18, and TGF-β1) in patients with IgG4-RD were observed. These data showed how elevated levels of such cytokines and chemokines can induce IgG4 plasma cell infiltration, high IgG4 levels in the periphery, and impact tissue fibrosis in the LSG of IgG4-RD patients (13). Further, researchers using high-throughput RNA sequencing technology revealed the molecular differences and effects from prednisone treatment among IgG4-related disease with salivary gland lesions (RD-SG), without SG lesions (RD-nonSG), and IgG4-related retroperitoneal fibrosis (RF) (14). However, the molecular mechanisms and appropriate therapeutic strategies underlying the pathogenesis of IgG4-RD are still unclear.
Proteins are effectors of biological function, and exert critical important roles in the pathogenesis of diseases. Investigation into proteins is crucial for the development of methods to realize early disease diagnosis, prognosis assessment and to monitor the disease development (15). Proteomics involves the applications of technologies for the identification and quantification of overall proteins present content of a cell, tissue or an organism. However, the proteomic study in the field of IgG4-RD is still blank. Thus, we detected both the serologic and tissue proteasome of IgG4-RD patients, and sought for the potential pathogenic information underlying the changes of expression level of proteins in IgG4-RD. Besides analysis at protein level, we also applied Abbreviations: IgG4, Immunoglobulin G4; iTRAQ, Isobaric tags for relative and absolute quantitation; CV, coefficient of variation; WGCNA, weighted gene coexpression network analysis; KEGG, Kyoto Encyclopedia of Genes and Genomes; DEP, differentially expressed protein; FDR, false discovery rate; LSG, labial salivary glands; PBMC, peripheral blood mononuclear cell. BCR, B cell receptor; PAMP, pathogen associated molecular patterns; DAMP, danger associated molecular patterns. PRR, pattern recognition receptor; MMP, matrix metalloproteinase.
Weighted correlation network analysis (WGCNA) methods, a powerful analysis tool that can be utilized for constructing a weighted correlation network and finding modules comprised of highly correlated genes (16), to analyze two public datasets at transcriptomic levels in involved tissue and PBMC, respectively, of IgG4-RD patients.
In conclusion, based on proteomic and transcriptomic analyses, we have not only identified several differently expressed proteins in serum and tissue samples from IgG4-RD patients compared with healthy people, but also illustrated some features of immuno-inflammatory reactions in IgG4-RD, which also helped provide information of its potential therapeutic targets. These results may provide clues to the elucidation of the pathogenesis of, and the development of therapeutic agents for IgG4-RD.

Patients and Treatment
We studied eight diagnostic serum samples from eight treatmentnaive IgG4-RD patients (25-70 years old) at Department of Rheumatology and Immunology, Wuhan Tongji Hospital ( Table 1). The diagnosis of IgG4-RD has been made according to diagnostic criteria for IgG4-related disease (IgG4-RD) (17). Meanwhile, eight serum samples from healthy controls (HC) were collected and stored at −80 • C for further analysis. Tissue samples of submandibular glands were obtained from two IgG4-RD patients and two non-IgG4-RD patients (adjacent normal edge of the surgical specimens), and stored in liquid nitrogen. All patients gave informed consent to the use of data records for research and to additional laboratory analysis on serum and tissue samples.

Sample Preparation
In this study, the Isobaric tags for relative and absolute quantitation (iTRAQ) technology was applied to investigate the proteasome of serum and tissue samples. First, the ProteoMiner Protein Enrichment Kit (Bio-rad laboratories, Hercules, CA, USA) was applied to deplete the high abundance proteins. Then, the protein solution (100 ug) with 8 M urea was diluted 4 times with 100 mM TEAB buffer. Trypsin Gold (Promega, Madison, WI, USA) was used to digest the proteins with the ratio of protein: trypsin = 40: 1 at 37 • C overnight. After trypsin digestion, peptides were desalted with Strata X C (Phenomenex), and vacuum-dried according to the manufacturer's protocol. Peptide labeling was performed by iTRAQ Reagent. Peptide Fractionation were realized with a HPLC Pump system (Shimadzu LC-20AB) coupled with a high pH RP column. Supernatants of fractions were loaded on UHPLC system (Thermo Scientific TM UltiMate TM 3000) equipped with a trap and an analytical column, and peptides separated from nanoHPLC were subjected into the tandem mass spectrometry QEXACTIVE HF X (Thermo Fisher Scientific, San Jose, CA) for data-dependent acquisition (DDA) detection by nano-electrospray ionization. The parameters for MS analysis are listed as following: electrospray voltage: 2.0 kV; precursor

Protein Identification and Quantification
Protein identification and quantification were realized by software IQuant (18). The propensity score matchings (PSMs) were pre-filtered with false discovery rate (FDR) ≤ 1% to assess the confidence of peptides. Then, the identified peptide sequences were assembled into proteins. After protein inference, the protein will be estimated with FDR ≤ 0.01.

Microarray Data Collection
We downloaded DNA microarray dataset GSE40568 and GSE66465 from Gene Expression Omnibus (GEO). Specifically, LSG samples in GSE40568 were obtained from Japanese patients with IgG4-RD (n = 5) as well as from Japanese patients with SS (n = 5) and HCs (n = 3) who had been followed up at the University of Tsukuba Hospital (Ibaraki, Japan), Tokyo Women's Medical University Hospital (Tokyo, Japan), and Kyushu University Hospital (Fukuoka, Japan) (13). PBMC samples from peripheral blood mononuclear cell (PBMC) of IgG4-RD were obtained from patients with IgG4-RD before (n = 2) and after steroid (n = 2) therapy who registered in the research project of the Research Program for Intractable Disease of the Ministry of Health, Labor, and Welfare (MHLW) of Japan and HCs (n = 4) (19).

WGCNA Analysis
The coefficient of variation (CV) of each gene were calculated after expression matrix were imported and normalized. Genes with CV >5% were log2 transformed, and the corresponding expression data was applied as input for WGCNA analysis. Then weighted co-expression networks were constructed by employing blockwiseModules function in the WGCNA package (https://horvath.genetics.ucla.edu/html/CoexpressionNetwork/ Rpackages/WGCNA/). In this study, we construct a scale-free network (R 2 = 0.9) based on the criteria that soft-thresholding power β were set as 20 (IgG4-RD LSG samples in GSE40568 dataset, Figure S1) and 12 (IgG4-RD PBMC sample in GSE66465 dataset, Figure S2) correspondingly. Genes, that possess edges with adjacency value of >0.2 in the module most correlated to IgG4-RD, were extracted for enrichment analysis.

Statistical Analysis
Log2 transformed data were used to calculate the difference of proteins between IgG4-RD and HC samples with "t-test" function in R package (version 3.5.1). Expression matrix from GEO datasets were extracted and normalized by using R package "GEOquery". WGCNA analysis was realized with R package "WGCNA" (version 1.68) (16). KEGG/GO analyses and network construction were realized and visualized with Cytoscape (version 3.4.0) and ClueGO plugin or R package "clusterprofiler" (version 3.12.0) (20,21). Information of targeted drugs of hub proteins were obtaind via Therapeutic Target Database (22). Without specific indication in the manuscript, all DEPs or genes in gene-module were input to make enrichment analyses.
To understand the function of these DEPs which may involve in IgG4-RD, functional enrichment analysis of these DEPs identified in tissue/serum between IgG4-RD and control samples were carried out. Results indicated that most tissue upregulated DEPs were involved in terms including immune related cells activation (e.g., immune response-activating cell surface receptor signaling pathway) and cell adhesion (e.g., leukocyte cell-cell adhesion), and infection related processes such as human immunodeficiency virus 1 infection, Epstein-Barr virus infection, and Salmonella infection etc. (Figure 2A). However, down-regulated DEPs in tissue were mainly involved in processed related to cell junction (e.g., cell junction assembly) ( Figure 2B). Given the number of downregulated DEPs in serum is small, we focus on upregulated DEPs for further analysis and terms such as protein activation cascade, platelet activation, and extracellular structure organization were outstood in both GO biological process (BP) and KEGG pathway (Figure 3).

Weighted Gene Co-expression Network Analysis of Transcriptomic Data
A total of 4,364 genes with coefficient of variation (CV) >5% as input for the construction of WGCNA using GSE40568 dataset from LSC samples. Combined with the topological overlap matrix with the hierarchical average linkage clustering method, we detected the 13 modules in IgG4-RD (n = 5), pSS (n = 5), and HC samples (n = 3) (Figure 4A). Among them, module "turquoise" with 934 genes showed strongest correlation with IgG4-RD phenotype (correlation = 0.81, p = 7e-5) ( Figure 4A,  Figures S3A,B). Then, a total of 354 genes with edge-adjacencyvalue more than 0.2 in module "turquoise" were exported for further functional enrichment analysis (Tables S3, S4).
Based on GSE66456 dataset from PBMC samples, we identified 2,306 genes with coefficient of variation (CV) >5% as input for WGCNA. In this study, 14 modules were detected by WGCNA in IgG4-RD before treatment (IgG4-RD_BT, n = 2), IgG4-RD after treatment (IgG4-RD_AT, n = 2) and HC samples (HC, n = 4) based on the criteria referred in method section ( Figure 4B). Among them, module "yellow" with 360 genes showed strongest correlation with IgG4-RD_BT phenotype ( Figure 2B, Figures S4A,B). Further, 247 genes with edge-adjacency-value more than 0.2 in module "yellow" were exported for functional enrichment (Tables S5, S6). In general, we discovered three major biological processes and several signaling pathways that may involve in IgG4-RD.

Antibody Mediated Autoimmune Responses Are the Character of IgG4-RD
Antigen-receptor mediated signaling pathway was significantly enriched in IgG4-RD tissues. Besides T cell-related signaling pathways, using proteomic data by KEGG analysis, we found several KEGG pathways related to B cell-related immune processes, such as "Fc gamma R-mediated phagocytosis, " "B cell receptor signaling pathway, " "Antigen processing and presentation, " and "Leukocyte transendothelial migration" were significantly enriched in IgG4-RD tissue (Figure 2A). Meanwhile, we also found term like "Fc gamma R-mediated phagocytosis" was enriched in genes from "turquoise" module in WGCNA analysis of the transcriptomic data of IgG4-RD LSG (Table S7). In addition, KEGG term "Systemic lupus erythematosus" was enriched in transcriptomic data of IgG4-RD derived PBMC (p < 0.0001, Table S8). These results indicated autoimmune characters of IgG4-RD, and also revealed an important role of antibodies in the autoimmune responses of IgG4-RD (23).

Potential Infection and Infection-Related Responses May Be the Trigger in IgG4-RD
Infectious agents are the main origin of pathogen associated molecular patterns (PAMPs), and can mediate the release of danger associated molecular patterns (DAMPs). Both PAMPs and DAMPs are ligands of pattern recognition receptors (PRRs), which can be found in many cellular components involved in immune reactions. Activation of PRRs can modulate the functional states of immune-related cells, which further influence the process of immune responses. Based on our proteomic data, many infection related terms such as "Human  immunodeficiency virus 1 infection, " "Epstein-Barr virus infection" and "Bacterial invasion of epithelial cells" were enriched in IgG4-RD tissue (Figure 2A, Table 2). Meanwhile, we also found genes from module "turquoise" in transcriptomic data also showed correlation with infection terms such as "Human papillomavirus infection, " "Epithelial cell signaling in Helicobacter pylori infection, " "Bacterial invasion of epithelial cells, " etc ( Table 3, Table S7). Most of these terms showed relationship with bacterial infection, which may echo the fact that most involved organs in IgG4-RD are exocrine organs, like pancreas, and submandibular glands (24). Thus, we may infer that infectious agents and infection related processes may play important roles in the development and progress of IgG4-RD.

Platelet Activation Were Observed in IgG4-RD Samples
In proteomic data, up-regulated DEPs from both serum and tissue samples of IgG4-RD patients were enriched in "Platelet activation" pathway (Figures 2A, 3). Meanwhile, Enrichment of term "Platelet activation" was also observed in genes derived from "turquoise" module from tissue transcriptomic dataset ( Figure 5). This result is consistent with the research on platelets that may serve as the immune components to possess modulate functions (25).

Multiple Signaling Pathways May Participate the Pathogenesis of IgG4-RD
Understanding the significantly enriched signaling pathways in IgG4-RD may help to search the potential therapeutic targets. In proteomic data, "Rap1 signaling pathway" and "NF-κB signaling pathway" were enriched in IgG4-RD tissue (Figure 2B), while terms such as "MAPK signaling pathway, " "PI3K-Akt signaling pathway, " "TGF-β signaling pathway, " "Ras signaling pathway, " and "Rap1 signaling pathway" were enriched in tissue transcriptomic data ( Table 4, Table S7). Term "TGF-β signaling pathway" was also enriched in transcriptomic data from IgG4-RD PBMC (Table S8). In this study, only "Rap1 signaling pathway" was observed enriched in both tissue proteomic data and IgG4-RD LSG samples in GSE40568 dataset (Figures 2B,  6). Activation of Rap1 signaling pathway can lead to the production of proinflammatory cytokines and modulate the expression level of MMPs, which are critical in the modulation of extracellular matrix, and influence the fibrogenic process (26). Beside terms illustrated above, there were also many other terms, like "Autophagy, " "Necroptosis, " etc., that were enriched in different datasets (Tables S7, S8).

Potential Therapeutic Target Identified by Tissue DEPs Related Biological Processes
In order to find out the potential therapeutic targets for IgG4-RD, we extracted all proteins (not only DEPs) involved in top 30 KEGG/GO BP terms enriched from all DEPs in our tissue proteomic data (Figure 7). One thousand six hundred seventy-four proteins in KEGG top 30 terms and, 2,291 proteins in GO BP top 30 terms were identified. Based on the involvements of these proteins in these different biological processes, we constructed networks and calculated the degree (number of connections, number of pathways in which the   (27)(28)(29). Medication targeting these proteins may help the treatment to IgG4-RD. Therefore, we also obtained the information of targeted drugs to the hub proteins from Therapeutic Target Database (Table 5) (22).

DISCUSSION
To our knowledge, the pathogenesis of IgG4-RD is still not clear and this is the first report providing new insights to help the illustration of potential pathogenic mechanisms underlying IgG4-RD based on proteomic and transcriptomic data. Existence of strong immune responses in IgG4-RD is well-known. Different T cell subsets may interact with B cell subsets in involved tissue, and further lead to the down-streaming fibrogenesis, which can also be partly reflected by our study (Figure 2) (30)(31)(32)(33)(34). In addition, this study indicated that other subclasses of immunoglobulin G could also be pathogenic in IgG4-RD (35). The membrane form of IgG4 is exactly the B cell receptor on IgG4+ plasmablasts/plasma cells, which can mediate the antigen capture and further lead to the antigen presentation from plasmablasts/plasma cells.
Our study detected significant enrichment of "Fc gamma R-mediated phagocytosis" process in IgG4-RD tissue (Figure 2A,  Table S7), which may point out the potential role in the pathogenesis of IgG4-RD. Activation of Fcγ receptors FIGURE 5 | Gene network of "Platelet activation" term derived from "turquoise" module in IgG4-RD tissue transcriptomic dataset, proteomic data from tissue, and serum.
(FcγR) on phagocytes can promote phagocytosis and the following antigen presentation process (36). Compared with high binding affinity of IgG1 and IgG3 to FcγR, the binding affinity of IgG4 to FcγR seems to be much milder. In IgG4-RD tissue, our proteomic study showed significant elevated IgG3 level, while IgG1 showed only the tendency of elevation. These clues indicate the potential role of IgG3 in IgG4-RD, at least, partly by Fc gamma receptor mediated phagocytosis, which can be another origin of antigen presentation, together with that mediated by IgG4+ plasmablasts/plama cells.
In this study, besides immunoglobulins produced by plasmablasts/plama cells, we found level of several cytokines, including CXCL13, IL-27, and IL-18, elevated significantly in IgG4-RD tissue. Cytokines exert essential effects on the immunoinflammatory process. IL-27 can act as antagonists to suppress Th1, Th2, Th9, and Th17 responses, while promote the proliferation and the expression of T-bet, EOMES, and IL-12Rβ2 associated with increased production of IFN-γ and cytotoxic activity (37). IL-18 is a member of IL-1 family, and involved not only in Th1 and NK cell activation, but also in Th2, IL-17-producing γδ T cells and macrophage activation (38). CXCL13 is critical for the recruitment of follicular Tfhs, and plasma CXCL13 can be a biomarker for germinal center activity (39,40). Thus, all these three cytokines may exert important regulatory effects on the immune responses in IgG4-RD. However, these three cytokines were not detected in our serum proteomic analysis, and these might result from their low concentration in serum, and the detection threshold for iTRAQ methods. Thus, more validations should be applied to detect the existence and levels of these cytokines in serum, and examine the availability as biomarkers in IgG4-RD.
Infectious agents can modulate the status of immune responses, and even be the triggers of systemic lupus erythematosus (SLE) and Sjögren's syndrome (41). Our analyses enriched a large amount of processes related to infectious agents in both proteomic (Figure 2A) and transcriptomic data (Tables S7, S8). We may infer that the most frequently involved organs in IgG4-RD including pancreas, salivary glands are mainly exocrine organs, which have more opportunity to contact infectious agents. Activation of toll-like receptors (TLRs) can involve in autoimmune responses indirectly via modulating innate immunity, and directly via modulating B cell signaling, and TLRs can even be the therapeutic targets for autoimmune connective tissue diseases (42,43). In our proteomic data, the elevation of TLR2 was detected in IgG4-RD samples (Table S1), which indicated the possibility, that TLR2 may have the potential to exert modulatory effects in IgG4-RD in vivo. Enrichment of platelet activation related processes in IgG4-RD patients was also observed in our study (Figures 2, 3,  5). As illustrated above, majority of the recent investigations focus on the pathogenetic roles of traditional immune cells like lymphocytes and macrophages (44). However, immune response is a complex process, and can be modulated by multiple factors. Recently, immunomodulatory effects of platelets have been observed in several immunoinflammatory conditions. Zhu et al. showed platelets first promoted the activation of Th1, Th17, and Treg cells, while only suppressed the immune response of Th1 and Th17 cells secondarily (45). Platelets can also express receptors including FcγRIIA, TLR4, and TLR9, which equips platelets the capacity to receive the stimuli, like immune complexes, PAMPs, and DAMPs, from microenvironments and response to them (46). Granules in platelets contain various growth factors, chemokines, and proinflammatory factors including TGF-β, EGF, CXCL12, HMGB1, and sCD40L etc., which can function in the modulation of immunoinflammatory processes (46,47).
Genes (or proteins) related to several signaling pathways were also detected in IgG4-RD samples (tissue or PBMC), including MAPK signaling pathway, PI3K-Akt signaling pathway, Ras signaling pathway, TGF-β signaling pathway, NF-κB signaling pathway and Rap1 signaling pathway (Figures 2A, 3, Table 4, Table S7). All these signaling pathway can regulate variety of biological process. Beside biological process described above, genes related to autophagy, necroptosis and adhesion molecules (e.g., "Focal adhesion, " "Cell adhesion molecules (CAMs), " etc.) were also detected in IgG4-RD samples. Interestingly, several terms related to malignancy (e.g., "Choline metabolism in cancer, " "Small cell lung cancer, " "Bladder cancer, " "Viral carcinogenesis, " and "Proteoglycans in cancer") were also enriched in IgG4-RD samples. All these clues can reflect some parts of the characteristics of IgG4-RD, and indicate the potential similarities of the pathogenesis mechanisms underlying IgG4-RD to the known mechanisms of the biological processes illustrated above.
How to treat IgG4-RD is another task we are facing. Despite glucocorticoids and rituximab are effective in the induction therapy, and application of conventional steroid sparing medication (e.g., azathioprine, mycophenolate mofetil, and methotrexate) can help control the disease, the problem of relapse and side effects after long-term usage of these medications are still remarkable (48)(49)(50). Based on the characteristics of IgG4-RD tissue proteomic data, we identified several hub proteins, which might be involved and play important roles in the pathogenesis of IgG4-RD. Most of these proteins are important signaling components in immune reactions, and some of these agents have been even applied into the clinical treatment of other diseases. Therefore, after the validation of the pathogenic roles of these proteins in IgG4-RD via laboratory experiments, these drugs have the potential to be therapeutic agents targeting IgG4-RD.
In summary, we provided the first integrative analysis of IgG4-RD via both proteomic and transcriptomic data, and described a landscape of biological processes of this mysterious disease, which indicated some potential pathogenic molecules and immunoinflammatory responses, and provided several potential therapeutic targets for the treatment of IgG4-RD. There are also limitations in our study. Firstly, despite the aim of our study is to provide the landscape at the level of mRNA and protein of IgG4-RD, the sample size of our analysis was relatively small, and we didn't enroll more IgG4-RD-like samples (e.g., tumors, and other rheumatic diseases), in our analyses; secondly, the transcriptomic data were originated from the published data of other centers, although analyses to them showed overlaps with our proteomic data, which may reflect the reliability of our analyses indirectly; thirdly, our analyses were based on the data of bulk samples, which cannot provide the information of specific cell types. Therefore, further studies with large sample size, and at single cell level are needed. Besides, other types of omics (e.g., lipidomics, metabolomics, and glycomics) and validation with laboratory experiments are also important to help us understand this mysterious disease more deeply.

ETHICS STATEMENT
The studies involving human participants were reviewed and approved by Tongji Hospital, Tongji Medical College, Huazhong University of Science, and Technology Institutional Review Board Approval. The ethics IRB ID is: TJ-C20151109. The patients/participants provided their written informed consent to participate in this study.

AUTHOR CONTRIBUTIONS
SC and YC made the sample collection, data analyses, and wrote this manuscript. SL and CY revised the manuscript and provided   Table S1 | Differentially expressed proteins between IgG4-RD patients and controls in tissue. Table S2 | Differentially expressed proteins between IgG4-RD patients and controls in serum.