Urinary peptidomics in early type 1 diabetes Title: Peptidomic analysis of urine from youths with early type 1 diabetes reveals novel bioactivity of uromodulin peptides in vitro

Chronic hyperglycemia is known to disrupt the proteolytic milieu, initiating compensatory and maladaptive pathways in the diabetic kidney. Such changes in intrarenal proteolysis are captured by the urinary peptidome. To elucidate the early kidney response to chronic hyperglycemia, we conducted a peptidomic investigation into urines from otherwise healthy youths with type 1 diabetes and their nondiabetic peers using unbiased and targeted mass spectrometry-based techniques. This cross-sectional study included two separate cohorts for the discovery (N = 30) and internal validation (N = 30) of differential peptide excretion. Peptide bioactivity was predicted using PeptideRanker and subsequently verified in vitro. Proteasix and the Nephroseq database were used to identify putative proteases responsible for peptide generation and examine their expression in diabetic nephropathy. A total of 6550 urinary peptides were identified in the discovery analysis. We chose to further examine a subset of 162 peptides that were quantified across all thirty samples. Of the 15 differentially excreted peptides (P < 0.05), seven derived from a C-terminal region (SGSVIDQSRVLNLGPITRK) of uromodulin, a kidney-specific protein. Increased excretion rates of five uromodulin peptides were replicated in the validation cohort using parallel reaction monitoring (P < 0.05). One of the validated peptides (SGSVIDQSRVLNLGPI) activated NFκB and AP-1 signalling, stimulated cytokine release, and enhanced neutrophil migration in vitro. In silico analyses highlighted several potential proteases such as hepsin, meprin A, and cathepsin B to be responsible for generating these peptides. In summary, we identified a urinary signature of uromodulin peptides associated with early type 1 diabetes before clinical manifestations of kidney disease and discovered novel bioactivity of uromodulin peptides in vitro. Our present findings lay the groundwork for future studies to validate peptide excretion in larger and broader populations, to investigate the role of bioactive uromodulin peptides in high glucose conditions, and to examine proteases that cleave uromodulin.


INTRODUCTION
Diabetes mellitus is a metabolic disorder characterized by chronic hyperglycemia. It arises from either an insufficient production of insulin as a result of an autoimmune attack on beta cells of the pancreas (referred to as type 1 diabetes) or insulin resistance as a result of decreased sensitivity (referred to as type 2 diabetes) (1). Type 1 diabetes typically emerges in childhood or adolescence; whereas type 2 diabetes typically appears later in life and is often complicated by hypertension, obesity, or cardiovascular disease at time of diagnosis (2). Despite differences in etiology, both types predispose individuals to serious complications that threaten survival and quality of life (3).
A common microvascular complication is diabetic kidney disease. It is in fact the leading cause of end-stage renal disease, accounting for 30-50% of incident cases in Canada and worldwide (4,5). The current standard of care is geared toward managing diabetes and delaying the onset of renal disease.

Landmark clinical trials in diabetes, namely the Diabetes Control and Complications Trial (DCCT) and
the Epidemiology of Diabetes Intervention and Complications (EDIC), have demonstrated that intensive glycemic control significantly reduces the risk of kidney disease (6); however, maintaining "normal" blood glucose levels may not be attainable for most individuals with diabetes. The search for better personalized and targeted therapies has been under way for decades; however, the greatest hurdle is our limited understanding of the early mechanisms responsible for initiating diabetic kidney disease.
Omics approaches have previously been employed to fill these gaps in knowledge. In a review of For the discovery phase, we performed power calculations to determine the appropriate total cohort size (N = 30) required to achieve a study power of 80% using the independent two-sample t-test and Benjamini-Hochberg (BH) multiple testing correction with a false detection rate (FDR) of 0.0001 and an effect size d of 2. We thus obtained urine samples from 15 youths with type 1 diabetes and 15 healthy age-and sex-matched controls for the discovery phase.
For the validation phase, we performed power calculations based on the observed effect sizes from the discovery phase. To achieve a study power of 80% and an FDR of 0.05, we determined that the minimum number of samples to validate the top peptide candidates as determined by the BH adjustment was 9 per group (Supplemental Table S1). To match the discovery analysis, we collected thirty urine samples for the validation cohort.
A total of 60 second-morning midstream urine samples were collected from 60 youths. All power calculations were performed with G*Power software.

Study Population
All sixty youths were under 19 years of age; had no history of hypertension, persistent microalbuminuria, renal disease, macrovascular disease, and chronic inflammatory disease; and had never used corticosteroid, anti-inflammatory, or anti-hypertensive drugs. In this study, youths with type 1 diabetes mellitus were considered to be in the earliest and uncomplicated stage of the natural history of diabetic kidney disease. Youths with diabetes were recruited and initially screen from multiple diabetes clinics in Urinary peptidomics in early type 1 diabetes for age and sex differences in youths (13). eGFR was determined using the Larsson's formula (eGFR = 77.24 × CysC -1.2623 ). For the discovery cohort, clinical characteristics were measured at the time of study enrolment for both groups ( Table 1) and at time of urine collection for youths with diabetes only (Supplemental Table S2). The mean time between the measurement of clinical characteristics and the urine collection was 2.7 ± 0.5 years and 2.5 ± 0.7 years for youths with diabetes and without diabetes, respectively. Groups were matched by age at time of urine collection and sex. For the validation cohort, we obtained clinical characteristics at or near the time of urine collection ( Table 2). Groups were matched based on age. The research ethics boards at the Hospital for Sick Children and Mount Sinai Hospital approved this study.

Collection, Handling, and Storage of Urines
We followed the Standard Protocol for Urine Collection and Storage, created by the Human Kidney and Urine Proteome Project (HKUPP) and the Human Proteome Organization (HUPO) (14). Following collection, all fresh urine samples were immediately kept at 4°C until further processing. Urine samples were centrifuged at 1000 g for 10 minutes to remove intact cells and debris. The supernatants were then aliquoted in either 15-or 50-mL Falcon tubes as well as 1.5-mL Eppendorf tubes and stored at -80°C. This initial processing step was completed within 3 hours of urine collection to obviate the need for urine preservatives and protease inhibitors. Samples were de-identified and randomized so that investigators were blinded to experimental groups during sample processing.
Plus hybrid quadrupole-orbitrap mass spectrometer (Thermo Fisher Scientific) using a nano-electrospray ionization source. Each fraction ran on a 60-minute gradient in data-dependent acquisition mode with full MS1 scan from 400 -1500 m/z with a resolution of 70,000 and subsequent MS2 scan of the top 12 parent ions with a resolution of 17,500. Xcalibur software (v. 3.0.63; Thermo Fisher Scientific) was utilized to generate RAW files of each run. Mass spectrometry data have been deposited onto the ProteomeXchange Consortium via the PRIDE (16) partner repository with the dataset identifier PXD012210 (http://www.ebi.ac.uk/pride/archive/login; username: reviewer85516@ebi.ac.uk; password: ht3PGbCg).
All mass spectrometry raw data were analyzed by MaxQuant software (version 1.5.3.28) and were searched against the human Uniprot FASTA database (version July 2016 with 42,158 protein entries) using the built-in Andromeda search engine (17). To study endogenous peptides, we selected the "unspecific" digestion mode. Cysteine carbamidomethylation was the only fixed modification. Variable modifications included methionine oxidation, proline oxidation, and N-terminal acetylation with a maximum number of 5 modifications per peptide. The false discovery rate was set to 1% for both proteins and peptides with a minimum length of 6 amino acids and was determined by searching a randomized database. Potential contaminants were allowed in the search. The initial peptide tolerance was set to 20 ppm against a small 'human-first-search' database. The main search peptide mass tolerance was 4.5 ppm, and the fragment mass MS/MS tolerance was set to 0.5 Da. Matching between runs was selected. Peptide intensities derived from extracted ion currents and were used as a measure of peptide quantification.
Peptidomic data analysis was done on Perseus software (version 1.5.5.3) (18). Non-human contaminants and reverse hits were manually checked and removed. We then filtered out peptides that were not quantified in all samples so that we could examine the most robust changes in peptide intensities between groups. Peptide intensities were log-transformed to approximate a normal distribution. Differential peptide intensities between groups were then determined using the independent t-test (P < 0.05), followed by the BH adjustment (Q < 0.05).

Downloaded from
We initially developed targeted methods for selected reaction monitoring to validate the observed differences in urinary peptide excretion (using Tier 2 analyses). The top seven transitions were chosen for each precursor on the basis of magnitude, order of intensity, co-elution, and peak shape. Linearity was determined for each of the detected crude heavy-labeled peptides (Supplemental Figure S1). We ultimately developed and performed a parallel reaction monitoring assay. Complete details about method development and optimization can be found in the Supplemental Methods. The master mix and list of transitions are summarized in Supplemental Table S3.
Overall, sample preparation of urine was similar to what was described for discovery peptidomics with some notable exceptions. After thawing, urine volumes were normalized to 20 μmol of creatinine.
Crude heavy-labeled peptide standards for each of the differentially excreted peptides were purchased from JPT Peptide Technologies and were spiked into the samples just prior to Vivaspin centrifugation.
Steps involving SCX-HPLC fractionation and OMIX tip desalting were not performed. Instead, peptides were analyzed at the SickKids Proteomics, Analytics, Robotics, and Chemical Biology Centre (SPARC Biocentre in Toronto, Canada) facility, where they were loaded onto Evotips and injected in duplicate onto an Evosep One nLC system coupled to a Thermo Scientific Q Exactive HF-X hybrid quadrupoleorbitrap mass spectrometer on a 21-minute unscheduled gradient. Raw files were imported into Skyline software (version 4.1). We then manually reviewed each precursor ion to ensure that all transitions were captured in the integration boundaries for relative quantification. Differential peptide excretion rates (deriving from peak area ratios of endogenous-to-heavy-labeled peptides) were determined using the Mann-Whitney test (P < 0.05) between groups in the validation phase.

Urinary peptidomics in early type 1 diabetes
uromodulin ELISA (DY5144) and ancillary kits (DY008) were purchased from R&D Systems. All plates were read at 450 and 540 nm using a Perkin Elmer 2103 EnVision Multilabel Plate Reader. Protein concentrations were extrapolated from optical densities using four-parameter logistic (4-PL) regression.
Urinary measurements were normalized to urinary creatinine concentrations.

Luciferase Activity Assay
HK-2 and human primary proximal tubular epithelial cells were transfected with Renilla luciferase control reporter vector pRL-TK and a luciferase reporter for either NF-κB or AP-1. Purified peptides (70%) were purchased from JPT Peptide Technologies and were used in all in vitro experiments.
Uromodulin protein (>95% purity) was purchased from Millipore Sigma (catalogue number AG733). 1 hour. Reporter activities were measured using the Promega dual-luciferase assay kit. The luciferase activity was normalized to the Renilla luciferase activity.

Neutrophil Chemotaxis Assay
Human primary neutrophils were isolated from peripheral blood from a healthy volunteer. Neutrophils were incubated with purified uromodulin peptide (1 µM) for 15 minutes at 37°C and then placed on a 5% bovine serum albumin-coated microscope cover glass (22 × 40 mm) for 10 minutes at 37°C. The cover glass was inverted onto Zigmond chambers: 100 μL of Hank's Balanced Salt Solution media in the right chamber, and 100 μL of media containing fMLP (1 µM) in the left chamber. We recorded neutrophil movement for 15 minutes at a rate of 0.05 frame per second using time-lapse video microscopy. Captured images were analyzed using cell-tracking Retrac software (version 2.1.01).

Cytokine Release from Kidney Tubular Cells
Human primary proximal tubular epithelial cells were serum-starved for 24 hours and subsequently treated with purified uromodulin peptides (0.1 to 10 µM) or LPS (1 µg/mL) for 24 hours. Conditioned media were collected, centrifuged at 2000 g for 10 minutes, and kept at either 4°C for immediate use or stored at -20°C. Human Interleukin (IL)-6 (ab178013) and IL-8 (ab214030) SimpleStep ELISA kits were purchased from Abcam. Conditioned media were analysed in duplicate according to kit instructions.

Bioinformatics
For the heatmap analysis, we converted log-transformed peptide intensities into z-scores and performed Euclidean hierarchical clustering using the R package pheatmap. PeptideRanker was used to predict bioactive peptides (23). Peptide sequences were aligned with precursor proteins using Peptide Extractor to visualize cleavage sites (9). We inferred protease activity using Proteasix (24). Structural information was extracted from the UniProtKB database. We used iceLogo to determine percent differences in amino acid proportion at each position and to visualize the consensus cleavage sites at the N-and C-termini of urinary peptides (25). The MEROPS database (version 12.1, https://www.ebi.ac.uk/merops/index.shtml) was manually searched to identify proteases with specific cleavage sites and to collate cleavage site sequence logo motifs (26, 27).

Statistics
Normal distribution of each variable was examined using the Shapiro-Wilks normality test. We compared clinical characteristics between groups using the t-test for continuous variables and the chi-squared test for categorical variables. Pearson correlations were performed on log-transformed urinary peptide excretion rates and uromodulin protein excretion. Differences in luciferase activities and protein concentrations were determined using the independent t-test. Differential excretion of cytokines and chemokines was determined using the Mann-Whitney test on positive, non-zero data only. Imputation of zero-value intensities was performed with Perseus software from a normal distribution with a downshift of 1.8 and width of 0.3 standard deviations to determine fold changes for the cleavage site analysis. Plots were created with R software.

Downloaded from
Overall, 6550 peptides deriving from 751 protein precursor groups were identified using our urinary peptidomics workflow (Figure 2A, Supplemental Tables S4 and S5). After the removal of false hits, unquantified peptides, and contaminants, our MaxQuant search provided label-free quantitation (LFQ) for 6323 peptides from 731 unique protein precursors (Supplemental Tables S6 and S7): 5708 peptides were found in youths with diabetes, 5011 peptides originated from controls, and 4396 peptides were common to both groups ( Figure 2B). Of the 615 peptides found exclusively in the non-diabetic group, only ten peptides were found in at least ten samples-all deriving from the C-peptide region of insulin (Supplemental Table S6). Notably, one C-peptide fragment (EAEDLQVGQVELGGGPGAGSLQP) was identified in all fifteen urine samples from youths without diabetes and in none from those with diabetes.
On the other hand, no peptide was exclusively identified in urine from all youths with diabetes.
The most frequently observed protein precursors, associated with at least 100 urinary peptides, were serum albumin, collagen alpha-1(I) chain, hemoglobin subunit alpha and subunit beta, keratin type II cytoskeletal 6A, and apolipoprotein A1. Together, these six proteins accounted for 17% of the quantified peptidome (Supplemental Table S7).
A subset of 162 peptides were found in all thirty samples (Supplemental Table S8), representing 2.6% of the total peptidome. Although this subset mainly consisted of medium-to-high-abundance peptides, the range of intensities spanned four orders of magnitude ( Figure 2C). As expected, lowabundance peptides were not consistently detected across samples. We examined the tissue origins of peptides using the Human Protein Atlas database (28). Roughly 80% of unique peptides derived from proteins with known expression in the kidneys. Kidney-specific or -enriched proteins accounted for 15% of the 162 peptides, compared to 4% of the total peptidome ( Figure 2D). We thus decided to focus our further analyses on these 162 peptides as they reflected the most robust and consistent changes in the peptidome.
Fifteen peptides were differentially excreted between study groups (P < 0.05), deriving from alpha-2-HS-glycoprotein, albumin, apolipoprotein A1, clusterin, latent transforming growth factor beta binding protein 4, Na + /K + -transporting ATPase subunit gamma, and uromodulin ( Figure 2E, Table 3). All by guest on May 6, 2020 but three peptides were more abundant in urines from youths with diabetes ( Figure 2F, Figure 3). After BH adjustment, four uromodulin peptides and one clusterin peptide remained statistically significant (Q < 0.05). Unsupervised hierarchical clustering based on normalized peptide intensities segregated the two groups based on diabetes status, leaving one youth with diabetes among the controls (Figure 2F).

Validation of differentially excreted peptides
Parallel reaction monitoring takes advantage of mass spectrometry-based techniques to quantify multiple peptides simultaneously on quadrupole-orbitrap hybrid instrument using heavy-labeled internal standards (29). We successfully developed targeted methods for the relative quantification of six peptides: UMOD-1, UMOD-2, UMOD-3, UMOD-4, UMOD-5, UMOD-7. Each of the thirty urine samples were analyzed in duplicate injections onto the mass spectrometer with median technical coefficients of variation less than 12% for each peptide (Supplemental Table S9). All peptides except UMOD-2 were differentially excreted between groups in the validation cohort ( Figure 4A, Table 4). Fold changes were similar to those from the discovery phase. Notably, UMOD-2 was the lone peptide whose urinary excretion decreased in diabetes (fold change of 0.75, P = 0.0672). We also demonstrated that total uromodulin protein excretion was markedly reduced in diabetes (fold change of 0.28, P = 0.0061), similar to UMOD-2 ( Figure 4B), although there was no significant correlation between peptide and protein excretion ( Figure 4C, Supplemental Table S10). We observed that the uromodulin peptides strongly correlated with one another, except for UMOD-2 ( Figure 4C, Supplemental Table S10). These findings suggest UMOD-2 may derive separately from the other uromodulin peptides. To better understand this, we sought to examine whether the downward trend of UMOD-2 would be reflected when applying other urinary normalization methods (Supplemental Table S11). In contrast to urinary excretion, urinary concentration of UMOD-2 was mildly increased in diabetes (fold change of 1.1, P = 0.7718). A similar trend was observed when UMOD-2 was normalized to total uromodulin protein (fold change of 2.8, P = 0.0561).
The direction of change of the other uromodulin peptides remained consistent despite different normalization methods.

Downloaded from
We next investigated whether these uromodulin peptides might exert a functional or biological role in the kidney. PeptideRanker was used to predict bioactivity of the differentially excreted peptides ( Table 5).
To verify the in silico prediction, we first performed NFκB and AP-1 luciferase activity assays in HK-2 cells. Both transcription factors are activated downstream of toll-like receptor 4 (TLR4) activation as part of the pro-inflammatory response (33). Chronic activation of inflammation is thought to contribute to diabetic kidney disease progression (34, 35). As an initial screen, cells were treated for 24 hours with graded concentrations (0.1, 1, and 10 μM) of 70%-purified uromodulin peptides, UMOD-1, UMOD-2, and UMOD-3 ( Figure 5A). Ang II was used as a positive control based on previous findings that it promotes inflammation and oxidative stress in the kidney (36-38). At 0.1 μM, UMOD-1 and UMOD-2 significantly stimulated NFκB activity (P < 0.05). As predicted, UMOD-3 failed to elicit a response, even at the highest tested dose of 10 μM.
Next, we conducted a neutrophil chemotaxis assay. Human blood neutrophils were first treated with one of the three previously tested uromodulin peptides, then placed in a Zigmond chamber, and finally subjected to a N-formyl-methionyl-leucyl-phenylalanine (fMLP) attractant gradient. Compared to untreated neutrophils or neutrophils treated with UMOD-3, UMOD-1 and UMOD-2 significantly enhanced neutrophil migration toward fMLP (P < 0.05) (Figure 5E).
To bridge NFκB signaling and neutrophil chemotaxis, we examined whether UMOD-1 might induce cytokine release from primary proximal tubular epithelial cells. Interleukin (IL)-8 was specifically chosen as a neutrophil chemokine; while IL-6 was selected for its pro-inflammatory action. Both cytokines are downstream products of TLR-4 mediated NFκB signaling. We demonstrated that UMOD-1 promotes IL-6 and IL-8 secretion into serum-starved conditioned media of primary proximal tubular epithelial cells at 0.1 to 10 μM (P < 0.05) (Figures 5F-G.

Urinary excretion of cytokines and chemokines
Our in vitro findings suggested a role for uromodulin peptides in inflammation. Accordingly, we speculated whether a pro-inflammatory signal on the basis of increased urinary excretion of cytokines and chemokines could be detected in early diabetes. To do so, we used a Luminex bead-based immunoassay and adjusted urinary analyte concentrations to urinary creatinine concentrations. Eight analytes were selected a priori based on the literature on urinary markers of inflammation in diabetic kidney disease: IL-1b, IL-6, IL-8, IL-18, IP-10, MCP-1, MIP-1b, and TNFa (Figure 6). IL-1b and TNFa, however, were omitted because more than two-thirds of the samples in each group could not be reliably quantified. Of the six remaining analytes, IL-6 (median fold change of 2.1), IL-18 (median fold change of 2.5), and IP-10 (median fold change of 2.8) were excreted to a greater extent by youths with diabetes, compared to their non-diabetic peers (P < 0.05).

Cleavage site analysis of urinary peptides
by guest on May 6, 2020 https://www.mcponline.org

Downloaded from
To investigate the proteolytic mechanisms altered in early diabetes, we first examined the cleavage sites associated with the N-termini and C-termini of urinary peptides. For this particular analysis, we initially included all 4396 peptides that were common to both groups and subsequently excluded peptides that were quantified in fewer than five samples per group. Missing LFQ data was imputed using the Perseus software. Of the 1941 eligible peptides, 549 peptides were increased (with a fold change of at least 1.5) and 331 peptides were decreased (with a fold change of 0.67 and below) in diabetes. We also added 94 peptides that were exclusively identified in one group if they were quantified in at least five samples: 58 from youths with diabetes and 36 from their non-diabetic peers. Thus, a total of 974 peptides was examined in our cleavage site analysis using iceLogo (25), in which the experimental set (N = 607 peptides that were elevated in diabetes) was compared to the reference set (N = 367 downregulated peptides) (Figure 7A).
We noted a preponderance of proline residues at or near the N-and C-termini of urinary peptides, especially peptides that were downregulated in diabetes (Figure 7B-C). In fact, proline residues accounted for 25% of all C-termini of downregulated peptides, compared to 14% of C-termini of peptides elevated in diabetes. Additionally, there was a higher percentage of trypsin-like cleavages (i.e., at the carboxyl side of an arginine or lysine residues) resulting in the N-terminus of urinary peptides downregulated in diabetes; the opposite trend seems to be true at the C-terminus of peptides elevated in diabetes (P < 0.05). Unlike the non-diabetic profile, the sequence motifs for diabetes were more diverse and highlighted a broader range of residues such as leucine, phenylalanine, and isoleucine. Notably, the top motif for the N-terminus cleavage site (L*LV) corresponded to renin (26), which is secreted by juxtaglomerular cells of the kidney.
We next focused on the proteolytic mechanisms responsible for the uromodulin peptides. Using the Peptide Extractor(9), we mapped all of the differentially excreted peptides to a specific region near the C-terminus of uromodulin: 589 SGSVIDQSRVLNLGPITRK 607 (Figure 7D). Adjacent regions include the hepsin consensus cleavage site at 586 RFRS 589 and a glycosylphosphatidylinositol (GPI) anchor site at S 614 , which anchors the uromodulin protein onto the cell membrane of the loop of Henle and distal tubular by guest on  cells. This peptide region contains the external hydrophobic patch (EHP) motif, 598 VLNLGPITRK 607 , which inhibits uromodulin polymerization and function (43). Given that urinary protein excretion of uromodulin was decreased in youths with diabetes ( Figure 4B), we hypothesize that changes in peptide excretion could reflect changes in proteolytic activity.
Proteasix (24) was used to predict the proteases responsible for generating the uromodulin peptides. A total of 156 combinations of predicted proteases and cleavage sites met the MEROPS threshold score. For each cleavage site, we selected the proteases according to the predicted probability as follows: 1) predicted probability fell above the 99 th percentile of the population distribution of all sequences, 2) predicted probability was below, but was the closest to the 99 th percentile, and 3) predicted probability was the furthest away from the MEROPS threshold score. After applying the criteria, nineteen protease/cleavage sites combinations remained ( Table 6), and the predicted cleavage sites are illustrated in Figure 7E. Of these predicted proteases, eight predominantly cleave at the carboxyl side of arginine or lysine (Supplemental Figure S5). According to our earlier sequence logo motif analyses, we speculate that the N-terminal cleavage of the peptide region ( 589 SGSVIDQSRVLNLGPITRK 607 ) occurs in normal, physiologic conditions (Figure 7B), whereas the C-terminal cleavage events are increased in diabetes ( Figure 7C). This is supported by an increased percentage of cleavage events that are at the carboxyl side of isoleucine (UMOD-1, UMOD-3, and UMOD-7), lysine (UMOD-6), threonine (UMOD-4) of uromodulin peptides elevated in diabetes.
We also searched the Nephroseq v5 database (www.nephroseq.org, November 2018, University of Michigan, Ann Arbor, MI) for expression data in diabetic nephropathy and found nine proteases whose expression was significantly altered in cases of diabetic nephropathy (P < 0.05): a disintegrin and metallo-proteinase domain-containing protein 10, cathepsin B granzyme A, hepsin, kallikrein-6, matrix metalloproteinase-2, matrix metalloproteinase-9, neprilysin, and plasmin. Hepsin has previously been identified as a putative protease responsible for the N-terminal cleavage of UMOD-1 and UMOD-2.
Meprin A and cathepsin B were predicted to cleave at carboxyl side of G 590 , resulting in the inert peptide UMOD-3. Neutrophil elastase was mapped to the C-terminus of UMOD-1, UMOD-3, and UMOD-7.

DISCUSSION
It is generally accepted that chronic hyperglycemia disrupts the proteolytic milieu (10-12). The resulting peptide products could thus be used to examine the upstream interactions between proteolytic enzymes and protein substrates (9). In this study, we report that otherwise healthy youths with type 1 diabetes excrete a urinary peptide profile that is distinct from their non-diabetic peers and that these differences precede the onset of microalbuminuria or renal function decline. Seven of the fifteen differentially excreted peptides originated from a small region ( 589 SGSVIDQSRVLNLGPITRK 607 ) near the C terminus of uromodulin; and five were validated in a second cohort by parallel reaction monitoring. Furthermore, we discovered that UMOD-1, one of the validated peptides, displays novel bioactivity in vitro, by stimulating TLR4-dependent NFκB luciferase activity in HK-2 and primary proximal tubular epithelial cells, by inducing cytokine release from primary proximal tubular epithelial cells, and by enhancing neutrophil chemotaxis. A second peptide, UMOD-2, also exhibited similar bioactivity, but was not differentially excreted in the validation cohort. Interestingly, uromodulin protein failed to elicit a proinflammatory response in tubular epithelial cells. In silico analyses identified several proteases that could be responsible for generating these uromodulin peptides, such as hepsin, meprin A, cathepsin B, and neutrophil elastase. Our findings present new insights into the early kidney response to chronic hyperglycemia long before the clinical manifestation of renal dysfunction.
Uromodulin, also known as Tamm-Horsfall glycoprotein, is exclusively expressed in the loop of Henle and distal tubules of the kidney (44,45). Naturally, its peptides have gained notoriety as potential kidney-specific biomarkers for diabetic nephropathy (46), hypertension (47, 48), chronic kidney disease (49), systematic lupus erythematosus(50), IgA nephropathy(51), and acute rejection following renal transplantation (52). Even though the exact peptide sequences differ between peptidomic studies, they typically originate from the same region as our peptides: 589 SGSVIDQSRVLNLGPITRK 607 . These observations suggest that perturbations in specific proteolytic events associated with kidney injury may not be unique to diabetes after all. Still, our findings demonstrate that the urinary peptidome may reveal by guest on May 6, 2020 early signs of kidney injury in youths with diabetes, even when they present with no clinical symptoms of renal stress or disease. Additional validation studies are needed to examine how urinary excretion of uromodulin peptides changes over the course of diabetic kidney disease.
Beyond their potential utility as indicators of injury, very little is known about these uromodulin peptides. We thus sought to investigate whether they may play a direct mechanistic role in the kidney; and, for the first time to our knowledge, we demonstrated that two peptides, UMOD-1 and UMOD-2, were capable of stimulating pro-inflammatory responses in kidney tubular epithelial cells and enhanced neutrophil migration. These findings significantly enrich our current understanding of uromodulin biology because the glycoprotein is thought to be immunologically inert in the tubules (40,53,54). Namely, uromodulin was shown to induce the maturation of myeloid dendritic cells by activating TLR4-dependent NFκB signalling pathways (39). Uromodulin nanoparticles stimulated the NLRP3-inflammasome in antigen-presenting cells, triggering IL-1β secretion (40). In addition, uromodulin protein expression on membranes of cultured tubular epithelial cells enhanced neutrophil adherence and trans-epithelial migration (55). Neutrophils also preferentially attached to uromodulin protein immobilized on microtiters plates; however, no such interaction was observed with soluble uromodulin (56). A recent study reported that uromodulin abated renal and systemic oxidative stress by inhibiting the transient receptor potential cation channel, subfamily M, member 2 (TRPM2), as a protective mechanism against acute kidney injury (57). Altogether, uromodulin protein seems to rely on interactions with immune cells to trigger inflammation and may be protective in the kidney. In contrast, our in vitro data suggest that uromodulin peptides could mediate pro-inflammatory processes in the kidney directly.
We next examined how these peptides might be generated. While basolateral expression has been observed, uromodulin is predominantly attached to the apical surface by a GPI anchor near the C terminus (58, 59). The bulk of the protein is subsequently shed into the urinary space following hepsin cleavage (60). More specifically, the type II transmembrane serine protease separates the active portion of the protein from our peptide region, which remains attached to the membrane and contains the inhibitory external hydrophobic patch (60-62). This particular cleavage event accounts for the N-terminus of the Two other notable cleavage sites emerged from our in silico analyses: the N-terminus of UMOD-3, UMOD-4, UMOD-5, and UMOD-6 by meprin A and cathepsin B; and the C-terminus of UMOD-1, UMOD-3, and UMOD-7 by neutrophil elastase. Meprin A is a metalloproteinase found in the brush border of the proximal tubules, where it is constitutively shed into the urinary space by furin (69). In its active form, meprin A stimulates pro-inflammatory cytokine release and fibrosis in the diabetic kidney (70)(71)(72)(73). Studies have also hinted at a protective role for meprin A, in which diminished expression and activity in the kidney were associated with later stages of diabetic kidney disease (74, 75). Urinary levels of meprin A were markedly elevated in microalbuminuric and late-stage disease, but were negligible in normoalbuminuric individuals with diabetes and non-diabetic healthy control subjects (76). Cathepsin B, a lysosome-associated enzyme, is involved in autophagy, extracellular matrix remodeling, inflammasome activation, and sodium reabsorption in the kidney (77)(78)(79)(80). Studies have shown that cathepsin B activity is decreased following in vitro stimulation of HK-2 cells with advanced glycation end-product (82) and in the kidneys of rats with streptozotocin-induced diabetes (81). Interestingly, lysosome-derived cathepsin B expression was elevated in glomeruli of a murine model of podocyte injury; and cathepsin B knockout mice were more resistant to injury and recovered more rapidly following injection of nephrotoxic serum (83). Urinary proteomic profiling in early diabetes has demonstrated increased urinary excretion of cathepsin B and other lysosome-associated enzymes (84,85). Lastly, neutrophil elastase is secreted from neutrophil granules during acute inflammation and promotes reactive oxygen species generation (86,87).

Urinary peptidomics in early type 1 diabetes
Neutrophil infiltration into the kidney has been well-documented in later stages of diabetic kidney disease and contributes to disease progression and fibrosis (88)(89)(90).
While it is possible that each peptide is generated de novo from uromodulin protein, we posit that they derive from peptide intermediates, similar to angiotensinogen and its peptide products. The 452residue precursor is first shortened into the inactive Ang I decapeptide ( 1 DRVYIHPFHL 10 ) by renin.
Subsequent angiotensin-converting enzyme activity produces the bioactive Ang II octapeptide ( 1 DRVYIHPF 8 ), which promotes sodium reabsorption via aldosterone secretion (91). Alternate processing of the Ang I peptide results in a myriad of bioactive and inert peptide sequences (92)(93)(94)(95). Notably, Ang (1-7) and Ang II (1-8) have opposing biological effects (96,97), underscoring the significance of a single residue. Likewise, the bioactive UMOD-1 (SGSVIDQSRVLNLGPI) and the inert UMOD-3 (SVIDQSRVLNLGPI) peptides differ in two residues at the N-terminus, suggesting that the loss of the serine-glycine residues at the N-terminus may render a peptide inactive and could act as a means to "turn off" peptide bioactivity. Interestingly, UMOD-2 excretion did not correlate with any of the other uromodulin peptides and may thus originate from a separate proteolytic mechanism. Although UMOD-2 was not differentially excreted in the validation cohort, it displayed a downward trend in diabetes. Its presence may be related to normal renal physiology or reflect decreased total uromodulin protein expression, rather than dysregulated processing in chronic hyperglycemia. Future studies will be necessary to define the proteolytic pathways and putative enzymes responsible for uromodulin processing.
Our study has several strengths. First, we carefully selected our study population to avoid possible confounding effects of concurrent medications and comorbidities. This allowed us to examine the early kidney response to chronic hyperglycemia in a homogenous cohort of youths with and without diabetes. Second, we conducted an internal validation of uromodulin peptides in a second, independent cohort. Finally, we provided three pieces of evidence of novel bioactivity of UMOD-1, whose urinary excretion was elevated in early type 1 diabetes, using three different cell systems. Our findings may open by guest on May 6, 2020 https://www.mcponline.org

Downloaded from
Urinary peptidomics in early type 1 diabetes a new field of investigation into the mechanisms responsible for early kidney injury in type 1 diabetes and into therapies that target prevention, rather than management, of diabetic kidney disease.
Our study also has some important limitations. First, the cause-effect relationships between uromodulin peptides and chronic hyperglycemia cannot be fully discerned in a cross-sectional study.
Second, the study population did not include individuals at later stages of diabetic kidney disease, and it is thus difficult to infer whether these peptides are associated with progressive disease. Finally, we were unable to experimentally verify the protease prediction due to the lack of a model system, which mimics the appropriate kidney cell type and naturally expresses uromodulin protein. Ultimately, an in-depth investigation into these proteolytic events will better define the pathophysiologic role of uromodulin peptides in the diabetic kidney.
In summary, we identified and validated a signature of uromodulin peptides associated with early, uncomplicated type 1 diabetes in two separate cohorts. Our discovery of novel bioactivity in uromodulin peptides suggests a potential role in mediating early changes in the diabetic kidney.        Fifteen peptides were differentially excreted (independent two-sample Student's t-test, P < 0.05, purple); however, only five survived the Benjamini-Hochberg adjustment (Q < 0.05, red). (F) Heatmap representation of the 15 differentially excreted peptides with unsupervised clustering of samples. Logtransformed peptide intensities were converted into z-scores using means, wherein higher scores are indicated in red and lower scores in blue. Figure 3. Boxplots for the fifteen differentially excreted peptides found in the discovery cohort (P < 0.05). Differential excretion was determined using the independent two-sample t-test after logtransforming peptide intensities (P < 0.05). HC, youths without T1D; T1D, youths with type 1 diabetes.   Figure 2B, peptides were eligible for inclusion if they were identified in at least five of the fifteen samples per group, where appropriate. Missing LFQ data was imputed using the Perseus software. Of the 1941 eligible peptides common to both groups, 549 peptides were increased (with a fold change of at least 1.5) in diabetes, and 331 peptides were decreased (with a fold change of 0.67 and below). A total of 974 peptides was examined in subsequent cleavage site analyses. (B-C) Sequence logo of N-and C-termini of urinary peptides using iceLogo. Significant differences in percentage of amino acids at each position in the cleavage are shown (P < 0.05). Positive values indicate increased percentages in type 1 diabetes. (D) Proteolytic map of uromodulin and the predicted proteases. Differentially excreted peptides originate from a region near the C-terminus of the precursor uromodulin protein, which is shown in bold at the bottom with predicted cleavage sites and proteases. Full-length uromodulin consists of 640 residues.