Systematic review of type 1 diabetes biomarkers reveals regulation in circulating proteins related to complement, lipid metabolism, and immune response

Pre-onset proteomic profiles

Our literature search identified 6 papers that investigated the temporal protein abundance changes in individuals with T1D. Studies by Moulder et al. [8], Fronhert et al. [9], Nakayasu et al. [10], and Webb-Robertson et al. (unpublished) looked at protein abundance changes at both pre- and post-seroconversion stages. In contrast, von Toerne et al. and Lui et al. examined the protein profile only after seroconversion [11, 12]. Moulder et al., von Toerne et al., Lui et al., and Nakayasu et al. used untargeted proteomics and identified 68, 26, 30, and 59 proteins, respectively, that were significantly different in post-seroconversion vs controls. Nakayasu et al., used targeted proteomics and validated 83 proteins. Webb-Robertson et al. used targeted proteomics, investigating the expression of 24 complement proteins in plasma from pre-seroconversion, and post-seroconversion subjects. In addition, von Toerne et al., and Fronhert et al. used targeted proteomics as a validation method to look at 3 and 5 unique proteins, respectively, whereas Liu et al. used ELISA as its validation step. In conclusion, these studies have identified many biomarker candidates, but with limited validation.

Post-diagnosis proteomic profiles

One of the first plasma/serum proteomics studies of individuals with T1D was performed by Metz et al. in 2008 [13]. They identified 5 differentially abundant proteins in recently diagnosed T1D patients compared to controls. Similarly, studies by Zhi et al. and Chen et al. utilized untargeted proteomics and identified 21 and 36 differentially abundant proteins, respectively, in sera from individuals with T1D compared to controls [14, 15]. Zhang et al. and Oliveira et al. performed untargeted proteomics of serum/plasma samples from individuals with T1D and identified 24, and 8 differentially abundant proteins, respectively [16, 17]. Zhang et al. tested the 24 proteins using targeted proteomics in 50 T1D vs. 100 controls, validating 16 proteins with high discriminating power. A subsequent blinded experiment in an independent cohort of 10 individuals with T1D and 10 controls identified the chemokine proplatelet basic factor (PPBP/CXCL7) and C1 inhibitor with 100% sensitivity and specificity to discriminate between the groups. Manjunatha et al. and Gourgari et al. performed untargeted proteomic analysis on high-density lipoproteins (HDL) and found a compositional but not level change of the HDL proteome in T1D individuals with high risk of cardiovascular complications [18, 19]. Overall, these studies showed proteomic changes in plasma profiles after T1D onset, which have the potential to be developed as diagnostic biomarkers.

Complement system

The complement system is a cascade of proteases making up a humoral extension of the innate immune system. Dysregulation of the complement pathway is linked to chronic and autoimmune diseases. Complement deficiencies are either inherited or acquired. Inherited deficiencies of C1-C4 are strongly associated with bacterial infection and systemic lupus erythematosus, while inherited deficiencies of C5-C9 are associated with bacterial infection and sepsis [20]. Acquired deficiency or factor level changes arise when activation or inflammation-related acute phase responses exist, resulting in either up or downstream exhaustion of some factor [21]. They also interact with adaptive immunity by forming complexes with antibodies, including islet autoantibodies [22]. A cytotoxic effect [23] of such complexes has been controversial, but this concept has gained a renewed interest with recent discovery of β-cell surface autoantibodies [24]. In individuals with T1D, complement components C3 and C4 are highly expressed in the pancreas, including the islets [25, 26]. Studies of pancreata obtained from cadaveric T1D donors have reported C4D immunostaining in blood vessel endothelium and exocrine ducts [25], a change typically associated with activation at that site, and significant upregulation of the complement cascade (C1QA, C1QB, C1QC, C1R, C1S, C3, C4B, C5, C5AR1, C6, C7, C8A, C8B, C8G, and especially C9) [27]. Their tissue compartment localization was extrapolated from transcriptomic data but remains uncertain. This information may be of limited value for our project as the tissue or plasma status of the complement system at the time of death in a person with long-standing T1D may not resemble that at the time of the appearance of islet autoantibodies in an otherwise healthy young child.

Importantly for studies of pre-diagnosis T1D, a functional relationship has been demonstrated between activation of components C3 and C5 and improved β-cell function in mice and humans [28], suggesting direct effects on β cells. The pro-inflammatory cytokines interleukin-1β and interferon-γ increase C3 expression in rodent and human β cells [29]. C3 silencing exacerbates β-cell apoptosis. On the other hand, upregulation of the complement system may improve β-cell autophagic response - a protective homeostatic response to the β-cell stress [29] that is impaired in T1D [30]. Exogenously added C3 protects against cytokine-induced β-cell death via AKT activation and c-Jun N-terminal kinase inhibition. While locally produced C3 is an important survival mechanism in β-cells under a pro-inflammatory assault. However, it is not known if a C3-focused therapy could slow or abort the progression of diabetes in humans. In addition, a variant of C3 gene is associated with T1D [31], which may suggest that the response to such intervention could be genetically modified. In sum, although changes in the complement system are clearly linked to the risk of T1D development and its rapid progression, the direction of the effect and the therapeutic implications are uncertain. Thus, we need to determine whether the system should be activated or modulated, what components of the pathway are most relevant to T1D development, and at what point in the evolution of the disease should a specific change in the pathway be introduced. Regarding complement therapeutics designed to block or modulate activation, there is a range of drugs that are either available or in clinical development. For example, these therapeutics will modulate the C3 and C5 convertases, thus dampening overall activation, or be more specific to target C5, C5a, C3a, or complement receptors for activation fragments (reviewed in [32]). Conversely, activation of the pathways is being explored for the treatment of infectious diseases, cancer, and disorders of metabolism [33].

The 13 papers we examined reported some aspect of the complement system as dysregulated, with 44 out of the 251 biomarker candidates identified from our KEGG analysis (Figure 4). C4 was the most identified protein, followed by C3 in this systematic review, and was reported to be primarily downregulated in post-seroconversion and post-diagnosis compared to controls [11, 14, 16, 19]. This was further corroborated by Webb-Robertson et al., where C3 levels were consistently low in pre- and post-seroconverted subjects. This unaltered level of C3 expression throughout the course of the disease was also reported by longitudinal studies conducted by Moulder et al. and Lui et al. These results are consistent with observations of deficiencies in downstream complement components coupled with increased abundance of the MAC inhibitor clusterin pre-seroconversion. However, there are conflicting reports of abundance following seroconversion. Previous ex-vivo characterization of the T1D pancreas corroborated our identification of complement dysregulation but instead found an increased abundance of complement markers following diagnosis [25, 27, 34]. Pre-clinical models of C3 knockout or receptor blockage leading to reduction in T1D development and other findings also suggest a role for the complement system in the T1D development [35-37]. While at first, the serum and tissue complement abundances appear to be at odds with one another, the low complement levels pre-seroconversion may be due to consumption through C3, C4, and MAC depositing in the pancreas throughout T1D development. Low C3 and C4 have been seen in both COVID-19 and lipodystrophy cases, and low C3 with high C5-C9 are common in glomerulopathy [38]. Evidence of lifelong complement deposits in the pancreas matching the parallel findings of lifelong low C3 levels provides further evidence of complement deposition in the pancreas as the driving force behind low blood levels of complement and, therefore, progressive loss of β cells due to increased immune response.

• Download figure
• Open in new tab

Figure 4. Complement cascade.

The diagram represents all the complement pathway proteins denoted by orange color identified in our systematic review. The Image was modified from “Complement cascade pathway” on the Reactome website (https://reactome.org/PathwayBrowser/#/R-HSA-166658 with StableID: R-HSA-166658).

In the context of T1D, the complement system may play a part in modulating adaptive immune response against islets, with HLA class II genes being associated with T1D risk [39]. Activation of the complement cascade and deposition of complement factors into the pancreas have been reported during insulitis [22, 27]. In late-stage diabetes, elevated complement levels in serum are linked to diabetic nephropathy [40]. Overall, the disruption of complement pathway seems to be a characteristic trait of T1D, however, further studies are warranted to help navigate the path to consistent biomarker or drug development.

Immune pathways

A recurrent theme among the T1D biomarker candidates is the enrichment of proteins related to antigen processing and presentation. For instance, the complement pathway can opsonize pathogens and dead cells to be presented to the antigen-presenting cells. Other proteins related to antigen opsonization are also regulated in T1D, such as antibodies and opsonization/scavenger receptors. Once the antigen is opsonized, it is phagocyted into phagosomes. In fact, β cells infected with coxsackievirus are efficiently phagocyted by dendritic cells, making coxsackievirus infection a potential trigger of the islet autoimmune response [41, 42]. The phagocytic process requires an extensive cell cytoskeleton remodeling [43], which was another pathway enriched in our analysis. The phagosome can be next fused to lysosomes to initially process the antigens, which are further processed in the proteasome and loaded into HLA for presenting to T cells [44]. HLA alleles represent the main risk factor of T1D development, further supporting that this pathway is involved in the autoimmune response [44].

Another process that occurs in parallel is the cytokine and chemokine signaling [45]. Among the biomarker candidates, the chemokine PPBP/CXCL7 has been identified along with 3 other immunoregulatory molecules (Figure 5). Despite all the cytokines/chemokines regulated in T1D, little is known about their mechanistic roles in disease development. Cytokine/chemokine and even phagocytosis can trigger signaling cascades in the cells that further regulate these processes but also leads to the expression of other effector molecules. Among signaling proteins, 4 kinases, phosphatases, and phospho-binding proteins have been described (Figure 5). In addition, the transcription factor TNIP1 has also been shown to be regulated in T1D (Figure 5). Regarding the effector molecules, oxidative stress proteins myeloperoxidase, glutathione peroxidase 3, peroxiredoxin-1, and sulfhydryl oxidase 1, were also shown to be regulated in T1D. Oxidative stress has been shown to induce β-cell dysfunction and death and has also been proposed as a potential therapeutic target [46].

• Download figure
• Open in new tab

Figure 5. Immune pathways.

The diagram represents proteins identified as the extracellular matrix, cytoskeleton/actin filament, oxidative stress, gene expression, inflammatory signaling, antigen presentation, cytokine/chemokine, opsonization, antibodies, and other immune receptors/regulators.

A common feature in the plasma of individuals with T1D is the regulation of extracellular matrix proteins, of which 27 have been described to be regulated (Figure 5). The extracellular matrix is an integral part of the immune response regulated by cytokines and chemokines [47, 48]. For instance, the extracellular matrix peri-islet basement membrane serves as a barrier, protecting islets from immune cell infiltration in insulitis in mouse models of T1D [49]. In addition, preservation of the extracellular matrix by administration of dextran sulfate, a mimic molecule of the extracellular matrix proteoglycans, has been shown to protect β cells and to be a potential treatment for T1D in mice [50].

Plasma lipoproteins

The DAVID analysis showed ∼5% of the reported proteins to be key players in lipid metabolism. Most circulating plasma lipids are packaged into lipoproteins which are traditionally classified into four common subfractions based on particle density: chylomicrons (CM), very low-density (VLDL), low-density (LDL), and high-density lipoproteins (HDL). Structurally dynamic apolipoprotein scaffolds reside at the water-lipid interface of all subclasses where they modulate particle interactions with plasma enzymes, cofactors, and cell surface receptors that continuously remodel the lipoproteins throughout their lifespan. Though traditionally defined based on their “cholesterol” content, proteomics studies over the last decade have revealed significant compositional heterogeneity exists within the lipoprotein subfractions which contain upwards of 273 different proteins [51] with the HDL subfraction accounting for >250 of these proteins. Thus, lipoproteins are thought to consist of a variety of compositionally distinct subspecies which have now been shown to modulate a diverse array of metabolic pathways [52, 53].

Though DAVID analysis implicated “cholesterol metabolism,” and by proxy lipoproteins, the analysis fails to capture the full lipoproteome due to the recency of lipoprotein molecular profiling studies in the literature. When compared to a lipoprotein-specific database [51], we found nearly half (N = 112) of the protein biomarkers identified in individuals with T1D are associated with lipoproteins (Figure 6A, ESM3), and 66 were more specifically associated to HDL (ESM4). Approximately 85% of these changes were unique to post-seroconversion and post-diagnosis indicating the most profound changes in lipoprotein metabolism occurs later in the disease process (Figure 6B). A total of 17 members were altered pre-seroconversion with most overlapping with previously discussed immune response and complement cascade (Figure 6C). Outside of the immune and complement proteins, we noted a few well-studied APOCs and clusterin were altered pre-seroconversion (Figure 6D), hinting some changes occur in traditional lipoprotein metabolism pathways prior to onset of dysglycemia or hyperglycemia.

• Download figure
• Open in new tab

Figure 6. Overlap of T1D-relevant proteins with HDL/LDL lipoproteome. A.

Venn diagram showing common protein hits (grey) between the HDL/LDL proteome [51] (black) and those reported in our systematic review (white). B. Venn diagram detailing the T1D stage as pre-seroconversion (pre-sero), post-seroconversion (post-sero) and post-diagnosis of the 112 HDL/LDL/T1D shared proteins of figure A. C. Gene Ontology of HDL lipoprotein-known functions of pre-seroconversion reported proteins. D. Heatmap of apolipoproteins with altered levels reported in at least one stage of T1D development. Panel C was modified from a published figure by Davidson et al. (copyright permission was obtained from the publisher, license number: 5433910818199) [51].

Most of the lipoproteome members altered post-seroconversion and post-diagnosis have documented roles in triacylglycerol metabolism. Perhaps the most robust of these observations was the associated with changes in plasma APOA4 altered in 8 studies post-seroconversion and post-T1D diagnosis [8, 11, 14, 18, 19]. APOA4 is well-documented to modulate the triacylglycerol packaging in the triacylglycerol-rich lipoproteins (CMs and VLDL). [54]. Additionally, APOA4 is reported to play key roles in satiety, gastric function, and glucose homeostasis [55, 56], all of which have been reported altered in individuals with T1D [56-60]. Two post-seroconversion and post-diagnosis studies reported increased plasma APOA2; a well-known HDL scaffold protein. While HDL which has little triacylglycerol, APOA2 has been shown to be implicated in triacylglycerol metabolism [61] though the mechanism is still poorly understood. Several studies report changes in the APOCs and APOE [8, 11, 18] which are also thought to modulate triacylglycerol lipolysis on VLDLs [62, 63]. These observations are in-line with elevated triacylglycerol levels and inhibition of lipoprotein lipase associated with the innate immune response [64].

Most of the biomarker studies were performed on whole plasma. As most apolipoproteins are exchangeable, more detailed lipoprotein speciation studies are required to determine the subclass on which these particles are located and how they are modulating particle function in the context of T1D. Though two studies attempted to speciate lipoproteins from individuals with T1D, both limited their analysis to HDLs (21, 22) thus missing most of the changes involved with the triacylglycerol-rich particles. Future studies that examine the temporal changes in the lipoproteome across all subclasses will better inform on the role of these pleiotropic particles and how they cross-communicate with the aforementioned complement system and immune pathways in T1D.