Skip to main content
medRxiv
  • Home
  • About
  • Submit
  • ALERTS / RSS
Advanced Search

Integration of urine proteomics and renal single-cell genomics identifies an interferon-gamma response gradient in lupus nephritis

View ORCID ProfileAndrea Fava, Jill Buyon, Chandra Mohan, Ting Zhang, H. Michael Belmont, Peter Izmirly, Robert Clancy, Jose Monroy Trujillo, Derek Fine, Yuji Zhang, Laurence Magder, Deepak A. Rao, Arnon Arazi, Celine C. Berthier, Anne Davidson, Betty Diamond, Nir Hacohen, David Wofsy, the Accelerating Medicines Partnership in SLE network, Soumya Raychaudhuri, Michelle Petri
doi: https://doi.org/10.1101/2020.03.13.20034348
Andrea Fava
1Division of Rheumatology, Johns Hopkins University School of Medicine, Baltimore, MD, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Andrea Fava
  • For correspondence: afava1@jh.edu
Jill Buyon
2New York University School of Medicine, New York, New York, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Chandra Mohan
3University of Houston, Houston, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ting Zhang
3University of Houston, Houston, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
H. Michael Belmont
2New York University School of Medicine, New York, New York, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Peter Izmirly
2New York University School of Medicine, New York, New York, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Robert Clancy
2New York University School of Medicine, New York, New York, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jose Monroy Trujillo
4Division of Nephrology, Johns Hopkins University, Baltimore, MD, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Derek Fine
4Division of Nephrology, Johns Hopkins University, Baltimore, MD, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Yuji Zhang
5Department of Epidemiology and Public Health, University of Maryland, Baltimore, MD, USA
6University of Maryland Marlene and Stewart Greenebaum Comprehensive Cancer Center, Baltimore, Maryland
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Laurence Magder
5Department of Epidemiology and Public Health, University of Maryland, Baltimore, MD, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Deepak A. Rao
7Division of Rheumatology, Immunology, Allergy, Department of Medicine, Brigham and Women’s Hospital, Harvard Medical School, Boston, MA, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Arnon Arazi
8Broad Institute of MIT and Harvard, Cambridge, MA, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Celine C. Berthier
9Internal Medicine, Department of Nephrology, University of Michigan, Ann Arbor, MI, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Anne Davidson
10Center for Autoimmune and Musculoskeletal Diseases, The Feinstein Institute for Medical Research, Northwell Health, Manhasset, NY, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Betty Diamond
10Center for Autoimmune and Musculoskeletal Diseases, The Feinstein Institute for Medical Research, Northwell Health, Manhasset, NY, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Nir Hacohen
8Broad Institute of MIT and Harvard, Cambridge, MA, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
David Wofsy
11Rheumatology Division, University of California San Francisco, San Francisco, CA, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
12Multiple Organizations
Soumya Raychaudhuri
13Center for Data Sciences, Brigham and Women’s Hospital, Boston, MA 02115, USA
14Division of Rheumatology and Genetics, Department of Medicine, Brigham and Women’s Hospital, Boston, MA 02115, USA
15Department of Biomedical Informatics, Harvard Medical School, Boston, MA 02115 USA
16Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, Cambridge, MA, USA
17Centre for Genetics and Genomics Versus Arthritis, Centre for Musculoskeletal Research, Manchester Academic Health Science Centre, The University of Manchester, Oxford Road, Manchester, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Michelle Petri
1Division of Rheumatology, Johns Hopkins University School of Medicine, Baltimore, MD, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Abstract
  • Full Text
  • Info/History
  • Metrics
  • Supplementary material
  • Data/Code
  • Preview PDF
Loading

Abstract

Lupus nephritis, one of the most serious manifestations of systemic lupus erythematosus (SLE), has both a heterogeneous clinical and pathological presentation. For example, proliferative nephritis identifies a more aggressive disease class that requires immunosuppression. However, the current classification system relies on the static appearance of histopathological morphology which does not capture differences in the inflammatory response. Therefore, a biomarker grounded in the disease biology is needed to understand the molecular heterogeneity of lupus nephritis and identify immunologic mechanism and pathways. Here, we analyzed the patterns of 1000 urine protein biomarkers in 30 patients with active lupus nephritis. We found that patients stratify over a chemokine gradient inducible by interferon-gamma. Higher values identified patients with proliferative lupus nephritis. After integrating the urine proteomics with the single-cell transcriptomics of kidney biopsies, it was observed that the urinary chemokines defining the gradient were predominantly produced by infiltrating CD8 T cells, along with natural killer and myeloid cells. The urine chemokine gradient significantly correlated with the number of kidney-infiltrating CD8 cells. These findings suggest that urine proteomics can capture the complex biology of the kidney in lupus nephritis. Patient-specific pathways may be noninvasively tracked in the urine in real time, enabling diagnosis and personalized treatment.

INTRODUCTION

Lupus nephritis is one of the most severe manifestations of systemic lupus erythematosus (SLE) (1). Lupus nephritis affects up to 50% of patients leading to end-stage renal disease in 10% and carries an 8-fold increase in mortality (2–8). Since nephritis is usually asymptomatic, patients with SLE are serially screened for the presence of abnormal urine protein (9). When proteinuria is elevated, a renal biopsy is obtained to confirm the diagnosis and guide treatment. Histopathology classifications subgroup lupus nephritis into six classes based on the presence, amount, and location of inflammation and fibrosis (10). For example, the presence of intraglomerular immune cell infiltration (class III and IV) identifies the most aggressive proliferative form of SLE nephritis, which is treated with stronger immunosuppression. Membranous glomerulonephritis (class V) is characterized by thickened glomerular basal membranes, subepithelial immune complex deposition and absence of intraglomerular infiltration. Endocapillary proliferation and membranous disease frequently overlap in the “mixed” phenotype.

Although the current treatment approach to lupus nephritis is loosely grounded in the morphological classification, findings from recent studies have challenged this paradigm. Kidney biopsies repeated after one year showed a transition in class in up to 70% of cases, from non-proliferative to proliferative and vice versa, with the second biopsy having better prognostic value (11, 12). These findings indicate that lupus nephritis is a dynamic process and therefore disease state inferred from kidney biopsies at one time point may have a limited value. In addition, recent studies have underscored the prognostic importance of interstitial disease which is inadequately addressed by current histologic classification schemes (13). Finally, routine histopathology of lupus nephritis cannot evaluate the underlying molecular pathways that may inform personalized treatment choices. Recent single-cell RNA-seq studies have revealed the complex network of immune cells, and their diversity, in the context of lupus nephritis (14–16). These studies offer promise for more personalized strategies for therapeutics and for prevention of renal fibrosis.

As patients with the same histological class may have dramatically different outcomes (17), there is an unmet need for a new biology-based classification of lupus nephritis that can dissect the heterogeneity. Ideally, such immunologic classification needs to be easily queried during the course of the disease to dynamically assess changes in response to intervention over months (rather than year). Several previous studies explored potential biomarkers to assess and predict lupus nephritis by using targeted and systematic proteomic approaches (18–21). However, these have not yet been standardized for clinical practice and do not consider lupus nephritis biological endotypes (18). In this study, we sought to identify patterns in the urine proteome that can identify distinct groups of SLE patients and infer the ongoing renal molecular pathology. We quantified 1000 urine protein biomarkers in active lupus nephritis and identified that patients stratify on an immune activation gradient. Integration with renal single-cell RNA sequencing revealed that the urine signature reflected intrarenal secretion by immune infiltrating cells and that interferon-gamma is a main driver of such increased immune activation.

RESULTS

Urine proteomics stratifies lupus nephritis patients on an immune activation gradient

This study was initiated to address whether we could classify lupus nephritis patients based on distinct biological patterns in the urine. We quantified 1000 analytes including cytokines, growth factors, and other soluble markers in the urine from 30 SLE patients with proteinuria and a kidney biopsy on the same day confirming lupus nephritis (class III, IV, or V) as part of Phase 1 of the Accelerating Medicine Partnership (AMP; Table 1 and figure S1) (14). To agnostically determine patterns of inflammation, we used principal component analysis (PCA), a dimension reduction technique that identifies major axes of variation in the urine proteome. For example, the first two principal components (PCs) clearly separated healthy donors from SLE patients based on distinct patterns of urinary protein excretion (Figure S2). Within lupus nephritis, the first principal component (PC1) stratified patients on a gradient rather than distinct clusters (Figure 1A). We noted that patients with higher PC1 values were more likely to have proliferative lupus nephritis, while the patients with pure membranous lupus nephritis had PC1 values close to 0 or negative (Figure 1A-B) suggesting that PC1 detected a biological response shared by many patients with proliferative lupus nephritis. No association was detected between PC1 and demographics or technical confounders such as batch, proteinuria, age or gender. PC2 was associated with site of urine collection (Figure S3) and therefore it was not evaluated further.

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 1.

Demographics and clinical characteristics.

Figure 1.
  • Download figure
  • Open in new tab
Figure 1. Lupus nephritis patients stratify over an immune activation gradient.

(A-B) PCA (n=30) of the first 2 principal components of the urine proteome (% variance explained is indicated). Patients with higher PC1 value have almost exclusively proliferative lupus nephritis (class III, IV, or mixed). P value calculated by t test. (C) Top and bottom 10 PC1 loading values of the measured urine protein. (D) Top 10 enriched pathways PC1 using Gene Ontology Biological Process indicating the biological significance of PC1.

Next, we sought to characterize the biological significance of the PC1 gradient. As principal components are determined by the weighted sum of all the analytes, we tested whether a molecular pathway was enriched using the weight (“loading”) of each analyte on PC1 (Figure 1C). Gene set enrichment analysis (GSEA) revealed that PC1 detected chemotaxis pathways (FDR < 0.01) (Figure 1D and S4): in particular, a pattern of chemokines secreted in response to IFN-gamma, IL1-beta, and TNF and directed to attract monocytes, NK, and CD8 T cells.

Together, these findings indicate that urine proteomics analyses identified a differential expression of a specific immune activation signature in patients with lupus nephritis and that stronger signals were associated with proliferative lupus nephritis.

Urine cytokines reflect intrarenal production by myeloid and cytotoxic lymphocytes

We asked whether the chemokines detected in the urine were indicative of intrarenal production by kidney infiltrating immune cells. To answer this question, we analyzed single-cell transcriptomics from 24 renal biopsies of patients with lupus nephritis (Figure 2A-B) (14). We carried out single-cell transcriptomic analysis on genes coding for the proteins found in two of the most enriched pathways in the urine PC1 (Chemokine-mediated signaling pathway GO:0070098 and Cellular response to interferon-gamma (GO:0071346)) (Figure S5) in the renal cells. A chemokine score was defined as the sum of the normalized expression of the aforementioned genes. While these chemokines were expressed by most infiltrating immune cells as well as epithelial cells, the dominant expression was observed in myeloid, NK, and CD8 T cells (Figure 2C-F).

Figure 2.
  • Download figure
  • Open in new tab
Figure 2. The chemokine pathway detected by urine proteomics reflect intrarenal production by myeloid, CD8, and NK cells.

(A) UMAP plot of single-cell RNA-sequencing of 24 renal biopsies of patients with active lupus nephritis (medium resolution clustering). (B) The identity of the cell clusters was determined by the expression of lineage markers. (C-F) A chemokine and an IFNγ response score based on the expression of the markers in the top urine PC1 pathways (“Chemokine-mediated signaling pathway GO:0070098” and “cellular response to interferon-gamma (GO:0071346)”) identified kidney infiltrating myeloid, CD8, and NK cells as the main producers of the chemokines in response to IFNγ detected in the urine. (G) Distribution of IFNγ, TNFα, and IL1β expressing cells. (H) Distribution of the total UMIs mapped to IFNγ, TNFα, and IL1β as a proxy for the amount of cytokine produced by each cell type.

To evaluate whether the proteins found in the urine derived from serum as a consequence of glomerular or tubular damage, we repeated the analysis after adjustments for albumin or beta 2 microglobulin, respectively. Results were virtually identical to the unadjusted analysis (Figure S6).

These findings suggest that the chemokines identified by urine proteomics derive from intrarenal chemokine production and, in particular, by myeloid, NK, and CD8 T cells.

IFN-gamma is mostly produced by infiltrating CD8 and NK cells in all patients with lupus nephritis

Pathway analysis indicated a chemokine pattern inducible by IFN-gamma, IL1-beta, or TNF. We quantified the kidney-infiltrating cells expressing such cytokines. Out of the three cytokines, IFN-gamma-positive cells were the most abundant, followed by TNF and IL1-beta (Figure 2G). IFN-gamma was mostly produced by CD8 T and NK cells, TNF by myeloid, CD8 and CD4 T cells, and IL1-beta almost exclusively by myeloid cells (Figure 2B, 2F, and 2I). Local transcription of these cytokines suggests that their transcriptional signatures in the kidney and protein levels in the urine are a direct consequence of intrarenal production. Since patients contributed a different number of cells to the dataset, we evaluated whether these findings were biased by a handful of patients contributing more cells. IFN-gamma+ cells were identified in all patients ranging between 4% and 45% of the total immune cells (mean 18%), TNF+ cells in 4-32% (mean 16%), and IL1-beta+ in 0-13% (mean 7%) (Table 2 and Figure S7). In contrast, we did not detect any transcription of type 1 interferons (Figure 3) or type 2 immunity (Th2) cytokines such as IL4, IL5, or IL13 (Figure S8), suggesting that a type 1 (Th1) immune response (22) is ubiquitous in all patients with lupus nephritis including non-proliferative histological classes and that IFN-gamma (type 2 interferon) is the main interferon being produced by immune cells in lupus kidney disease.

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 2.

Prevalence cytokine positive cells based on single cell RNA sequencing.

Figure 3.
  • Download figure
  • Open in new tab
Figure 3. Type 2 interferon (IFN-gamma), but not type 1 (IFN-alpha or IFN-beta), is expressed in lupus nephritis kidney.

Violin plots of the expression of type (IFNA1 and IFNB1) and type 2 (IFNG interferons by single cell RNA sequencing. No other type 1 interferon (i.e. IFNA2-21, IFNB3, IFNK, etc were detected (not shown).

Urine proteomics may predict the intrarenal cellular infiltrate

We then hypothesized that urine proteomics might be used to infer the composition of kidney infiltrating cell types in lupus nephritis. We analyzed six patients with matching urine proteomics and single cell transcriptomics on a renal biopsy collected the same day. Of these, four patients were diagnosed with proliferative lupus nephritis (including mixed) and two with pure membranous lupus nephritis based on histopathology (Figure 4A). Single-cell transcriptomics revealed a heterogenous cellular infiltrate regardless of histopathological class (Figure 4B). We found that urine proteomics PC1 values strongly correlated with the relative abundance of kidney infiltrating CD8 T cells (r2 = 0.84, p<0.01; Figure 4C). There was no statistically significant correlation with the other cell populations. While needing validation in a larger cohort, these findings suggest that urinary protein patterns would predict the cell composition of the renal immune infiltrate in lupus nephritis and that CD8 T cells, the major producer of IFN-gamma, may be the key contributor of the immune activation signature in the urine.

Figure 4.
  • Download figure
  • Open in new tab
Figure 4. Urine proteomics may predict the intrarenal cellular infiltrate.

Six patients had matching urine proteomic and single cell RNA sequencing data. Their coordinates on the urine proteomics PCA shown in Figure 1 are plotted in (A). (B) Distribution of the kidney infiltrating immune cells determined by single cell RNA sequencing. Patients are ordered based on their urine proteomics PC1 value as they appear in (A). (C) Correlation of the urine PC1 values and the frequency of kidney infiltrating CD8.

DISCUSSION

Urine is an ideal non-invasive specimen to repetitively assess lupus nephritis. It contains the byproducts of the ongoing renal biological processes in real time, is highly non-invasive, and can be easily monitored over time. Here, we show that urine proteomics in patients with lupus nephritis may detect biologically relevant pathways that are mirrored in the kidney tissue. We found that 1) patients with higher expression of urinary chemokine and immune activation pathways have almost exclusively proliferative lupus nephritis; 2) these chemokines are renally produced by infiltrating myeloid, NK, and CD8 T cells; and 3) urine proteomics may predict the composition of the renal cellular infiltrate.

The current classification of lupus nephritis relies on morphological features that may not fully capture the complex, dynamic, and heterogeneous biology. We showed that systematic, unbiased analysis of the urine proteome may reveal the active immunological pathways in the kidney. Not surprisingly, patients with a stronger chemokine signal in the urine have glomerular immune cell infiltration as seen in proliferative lupus nephritis. Strikingly, the chemokines detected are part of a specific immune response inducible by IFN-gamma, IL1-beta, and TNF, that matches the inflammatory environment in the kidney. The significant presence of IFN-gamma positive cells suggests a strong type 1 (Th1) immune response (22). In fact, the analysis of the single-cell transcriptome of kidney infiltrating immune cells revealed IFN-gamma production in all patients, but not type 2 (Th2) response cytokines (IL4, IL5, and IL13). Furthermore, the IFN-gamma inducible urinary PC1 strongly correlated with the predominance of infiltrating CD8 T cells, the major source of IFN-gamma in lupus nephritis. Altogether, these findings suggest an ongoing IFN-gamma driven response in the kidneys of patients with lupus nephritis, which can be detected in real time in the urine by measurements of chemokines. These chemokines are directed to attract mostly macrophages, neutrophils and lymphocytes, and especially other IFN-gamma producing cells via ligation of CXCR3, CXCR4, and CCR5 that may amplify the immune response in a feedforward loop.

We deliberately performed an unsupervised analysis to detect biology-driven differences rather than looking for markers that associated with pre-specified outcomes such as histological class. While a stronger immune activation signal in the urine identified patients with proliferative lupus nephritis, a substantial group of patients classified as proliferative were grouped with patients with non-proliferative disease. Compared to those with an elevated immune activation signal, these patients may represent those whose nephritis differs in the molecular pathways involved, or those in a different stage of the inflammatory process. It is also conceivable that investigation of the relevant molecular pathways would provide a better assessment of the disease process in lupus nephritis compared to the current classification system based on glomerular endocapillary infiltration in histology. Since histological class may vary in >50% of patients on re-biopsy, the dynamic identification of precise disease driving pathways may better inform treatment selection, development of new treatments, and patient selection for clinical trials. Urine proteomics might also obviate biopsy sampling error. The next phase of AMP studies is aimed to molecularly dissect lupus nephritis and likely reclassify it based on disease biology.

Our results highlight the importance of IFN-gamma in lupus nephritis as a potential driver and coverable therapeutic target. For example, an increase in IFN-gamma activity is the first abnormality detected years before the diagnosis of SLE, preceding the appearance of autoantibodies and the dysregulation of type 1 interferons (23). IFN-gamma receptor polymorphisms are linked to the risk of developing SLE and IFN-gamma polymorphisms resulting in higher IFN-gamma production that predisposes SLE patients to develop proliferative lupus nephritis rather than membranous (24, 25). The immune response in proliferative lupus nephritis has in fact been shown to be skewed toward the Th1 axis, which is IFN-gamma predominant (24, 26, 27). In mouse models, IFN-gamma is a crucial cytokine for the development of lupus nephritis and its blockage prevents and ameliorates kidney disease as well as reduces mortality (28–30). Strikingly, type 1 interferons are not required for the generation of germinal centers and the development of lupus nephritis in some mouse models, in contrast to IFN-gamma signaling (30). IFN-gamma has a very pleiotropic effect involving cellular events implicated in lupus nephritis, including the recruitment and activation of neutrophils, CD8 and CD4 T cells, NK cells, and macrophages (31). Importantly, this molecule induces the formation of germinal centers, Tfh development and class switch in B cells (32). In the gut, IFN-gamma mediates the death of intestinal stem cells suggesting that a similar mechanism may have implications for irreversible kidney damage and fibrosis (33). In humans, encouragingly positive phase II trials showed the efficacy of blocking the IFN-gamma pathway in extrarenal lupus (34– 36). In the ustekinumab trial, for example, clinical response was associated with reduction of the IFN-gamma signature, rather than type I interferon (37). In aggregate, these findings suggest that IFN-gamma is a central cytokine in lupus nephritis and, given the acceptable safety profile of its direct blockage (38), further studies in lupus nephritis should pursue IFN-gamma inhibition.

An interferon gene response, especially type 1, is the archetypal signature of SLE (39– 41). We and others have previously showed a strong type 1 interferon signature in lupus nephritis (14–16). There is a large overlap between the gene signatures induced by type 1 and type 2 interferons (42–44), as indicated by shared signaling pathways. Each interferon type can induce the production of the other one (43) eventually leading to stimulation from both sides and thus a mixed signature. Here, we showed the unequivocal presence of IFN-gamma-producing cells in lupus nephritis along a typical chemokine signature in the urine, suggesting that IFN-gamma is central to the pathogenesis of lupus nephritis. Type 1 interferon intrarenal transcripts were not detected in our studies, suggesting that, assuming the absence of technical limitations, the source of type 1 interferons may be outside the kidney or precedes the clinical events leading to a kidney biopsy.

We acknowledge the limitations of our study. We did not detect an association with clinical variables such as renal activity/chronicity or medication use. There are likely other major detectable molecular pathways that could identify distinct patient groups. However, the relatively limited sample size did not allow for further exploration without the risk of a type 1 error. Nevertheless, our findings are confirmed by other smaller and independent studies that assessed the concentration of a fraction of our analytes in the urine of patients with lupus nephritis (45). Further studies, as those ongoing with the Phase 2 of the Accelerating Medicines Partnership, are needed to extend and validate these results. Second, while our platform did not assay the whole urine proteome, it allowed for interpretability, a broad assessment of the immune response, and limited major confounders from batch effect and unalignable peptides as commonly seen in classical mass spectrometry experiments. To limit this bias, we utilized a self-contained enrichment strategy (GSEA). In order to identify protein patterns independently of the amount of proteinuria, we have scaled the analytes concentration within and between patients. As this might introduce a bias related to this dataset, a larger study will be needed to validate the findings and to develop a convenient assay to apply in clinical practice. However, the significant infiltrate of IFN-gamma-producing cells sustains the biological plausibility of our results. Finally, the scRNA-seq did not include certain cell types such as neutrophils, likely because they did not survive the freeze-thaw process.

In summary, we demonstrated that the complex molecular biology in the kidney biopsy could be captured by the urine in patients with lupus nephritis. These processes provide insight into patient-specific disease pathogenesis that can address the issue of heterogeneity in SLE. These pathways can be tracked in the urine in real time and may not only improve diagnosis but also guide selection, dosing and monitoring of treatment, thus paving the way to develop a biologically grounded liquid biopsy.

METHODS

Patients and sample collection

This study enrolled SLE patient with urine protein to creatinine ratio > 0.5 undergoing clinically indicated renal biopsy, under protocols approved institutional review boards at each site and with informed consents. Only patients with a pathology report confirming lupus nephritis were included in the study. Renal biopsies were scored by a renal pathologist according to the International Society of Nephrology/Renal Pathology Society (ISN/RPS) guidelines and NIH activity and chronicity scales. Urine specimens were acquired at two clinical sites in the United States (Johns Hopkins University and New York University). The total urine volume was split into 2 50-ml Falcon tubes. Urine cells were pelleted by centrifugation at 200g for 10 min and the supernatant was aliquoted and stored at -80 °C.

Urine Quantibody assay

The Kiloplex Quantibody protein array platform (QAH-CAA-X00, Raybiotech Life, Norcross, GA) was used for screening urine samples. The array was spotted with 1000 capture antibodies specific for 1000 different proteins in quadruplicate. All urine samples were clarified by centrifugation, and then diluted to yield a total protein concentration within the working range (500 ug/m - 1 mg/mL) before application to the arrays. Briefly, protein standards and urine samples were incubated on the array for 2 hours to allow the proteins to bind to the capture antibodies. A biotinylated antibody cocktail comprised of 1000 detection antibodies was subsequently added for incubation for 2 hours. Finally, streptavidin-Cy3 was added and left to incubate for 1 hour. Washing was performed between each step to remove the unbound reagents. After a final wash and dry, the slides were read with a fluorescent scanner, and data were extracted from the image using vendor-provided GAL file with a suitable microarray analysis software. Creatinine was measured for each urine sample (KGE005, R&D Systems, Inc., Minneapolis, MN), and all data were creatinine normalized before analysis. The complete list of the analytes measured in the urine is in the supplementary file “Kiloplex targets.csv”.

Renal tissue single cell RNA sequencing

Renal tissue was collected, stored and processed as previously described (14). Briefly, research biopsy cores were collected from consented subjects as an additional biopsy pass or tissue from routine clinical passes. Only biopsies with confirmed lupus nephritis were included. Kidney tissue was frozen on site and shipped to a central processing location where it was thawed and disaggregated. Individual cells were retrieved and sorted by flow cytometry enriching for CD45+ leukocytes. For each single cell, the whole gene expression profile was sequenced using the CEL-Seq2 method (46). Data are publicly available (https://www.immport.org/resources/amp-ra-sle).

Statistics

Principal component analysis

Urine protein concentrations were scaled (Z-normalized) within and across patients after log normalization. Scaling within samples favors the identification of the sample-specific pattern of expression regardless of the total urine protein concentration. Principal components were computed using the R Stats package version 3.5.2.

Pathway enrichment analysis

GSEA (47) was performed using the genes coding for the protein assessed by the Quantibody assay using the Gene Ontology Biological Process (48, 49). Pathways with at least 10 (but less than 500) genes represented in the Quantibody assay were included in the analysis. Enrichment p-values were computed based on 10,000 permutations and adjusted for multiple comparisons using the Benjamini-Hochberg procedure (i.e., false discovery rate) (50). The pathways captured by PC1 were determined based on the loading of each analyte on PC1.

Data Availability

Data are publicly available.

https://www.immport.org/resources/amp-ra-sle

AUTHOR CONTRIBUTION

A.F., M.P., J.B., and C.M. ideated the project. A.F, M.P, D.F, J.M., H.M.B., P.I., R.C, and B.D. acquired samples and supervised the processing. T.Z. and C.M. generated the proteomic data. A.F., M.P., J.B., D.F, J.M., H.M.B., P.I., R.C, A.A., D.A.R., C.C.B, N.H, D.W., and B.D contributed to the sample collection, processing, analysis and interpretation of the renal biopsy scRNAseq. A.F. devised the analytical plan, interpreted the results, and wrote the manuscript with the guidance of M.P, S.R., Y.Z., and L.M, and the contribution of J.B, D.A.R., R.C., H.M.B, and P.I. See Supplemental Acknowledgments for consortium details.

SUPPLEMENTARY FILES

Kiloplex targets.csv

Supplementary Figure 1.
  • Download figure
  • Open in new tab
Supplementary Figure 1. Overview of the urine biomarkers.

Heatmap showing the relative concentration of the 1000 measure analytes (scaled by subject and analyte). Hierarchical clustering based on Ward’s rule and Euclidean distance.

Supplementary Figure 2.
  • Download figure
  • Open in new tab
Supplementary Figure 2. Principal component analysis of all urine analytes separates healthy subjects from lupus nephritis.

PCA plot (n=37) of the first 2 principal components with % variance explained. Principal components are calculated independently in each analysis: PC1 in this figure does not correspond to PC in Figure 1.

Supplementary Figure 3.
  • Download figure
  • Open in new tab
Supplementary Figure 3. Principal component 2 detects collection site.

(A-B) PCA (n=30) of the first 2 principal components of the urine proteome (% variance explained is indicated). PC2, but not PC1, was associated with the urine collection site. There was no biological pathway significantly enriched in PC2.

Supplementary Figure 4.
  • Download figure
  • Open in new tab
Supplementary Figure 4. Complete list of pathways enriched in the urine proteome PC1.

2018 Gene Ontology Biological Process pathways with a FDR <5% are ordered by their normalized enrichment score.

Supplementary Figure 5.
  • Download figure
  • Open in new tab
Supplementary Figure 5. Definition of a chemokine score.

First, the genes coding for the proteins defining the leading edge of “Chemokine-mediated signaling pathway GO:0070098”, the most enriched pathway in urine PC1 were identified. Then, a score for each cell was defined based on the sum of the normalized expression at the single cell level in the kidney. The same operation was repeated for “Cellular response to interferon-gamma (GO:0071346)”.

Supplementary Figure 6.
  • Download figure
  • Open in new tab
Supplementary Figure 6. Adjustment for glomerular (albumin) or tubular (B2M) dysfunction.

Urine analyte concentrations were normalized for the concentration of albumin (A-D), beta-2 microglobulin (E-H), or both (I-M) in order to account for non-specific protein excretion through the glomerulus, tubule, or both, respectively. Panels A,B, E, F, I, and J show the PCA (n=30) of the first 2 principal components of the urine proteome (% variance explained is indicated). Panels C, G, and L showed the top and bottom 10 PC1 loading values of the measured urine protein. Panels D, H, and M the top 10 enriched pathways PC1 using Gene Ontology Biological Process indicating the biological significance of PC1.

Supplementary Figure 6.
  • Download figure
  • Open in new tab
Supplementary Figure 6. Distribution of cytokine positive cells among patients.

Cells with positive cytokine UMIs were quantified in the renal biopsy single cell transcriptomics. The bar plot summarizes the individual contribution from samples with at least 10 high quality cells. As detailed in table S1, cytokine positive cells were found in all patients.

Supplementary Figure 8.
  • Download figure
  • Open in new tab
Supplementary Figure 8. No expression of type 2 (Th2) immunity cytokines was detected in lupus nephritis kidney.

Violin plots showing that there was abundant expression of the type 1 (Th1) immunity cytokine IFNG but none of the type 2 (Th2) immunity cytokines (IL4, IL5, and IL13) were detected in lupus nephritis kidney by single cell RNA sequencing.

ACKNOWLEDGMENTS

We thank the participating Accelerating Medicines Partnership network clinical sites and participants. We thank Ilya Korsunski, Maria Gutierrez, Fan Zhang, and Felipe Andrade for their critical review of the analytical plan and help with the interpretation of the results.

Footnotes

  • Funding: This work was supported by the Accelerating Medicines Partnership (AMP) Rheumatoid Arthritis and Lupus Network. AMP is a public-private partnership (AbbVie Inc., Arthritis Foundation, Bristol-Myers Squibb Company, Lupus Foundation of America, Lupus Research Alliance, Merck Sharp & Dohme Corp., National Institute of Allergy and Infectious Diseases, National Institute of Arthritis and Musculoskeletal and Skin Diseases, Pfizer Inc., Rheumatology Research Foundation, Sanofi and Takeda Pharmaceuticals International, Inc.) created to develop new ways of identifying and validating promising biological targets for diagnostics and drug development Funding was provided through grants from the National Institutes of Health (UH2-AR067676, UH2-AR067677, UH2-AR067679, UH2-AR067681, UH2-AR067685, UH2- AR067688, UH2-AR067689, UH2-AR067690, UH2-AR067691, UH2-AR067694, UM2- AR067678, and AR074096). The Hopkins Lupus Cohort is funded by NIH AR 69572.

  • Disclosures: A. Fava, None; Y. Zhang, None; J. Buyon, Bristol Myers Squibb, Exagen Diagnostics; H. Belmont Exagen and Ra Pharmaceuticals; P. Izmirly, None; R. Clancy, None; C. Mohan, Equillium, Inc; T. Zhang, None; J. Monroy Trujillo, None; D. Fine, None; A. Arazi, None; D.A. Rao, None; C.C Berthier, None; N Hacohen, None; B. Diamond, None; the Accelerating Medicines Partnership, None; S Raychaudhuri, None; M. Petri, AbbVie, Aleon, Amgen, Annenberg Center for Health Sciences, AstraZeneca, Blackrock, BMS, Decision Resources, Exagen, Glenmark, GSK, INOVA, IQVIA, Janssen, Lilly, Merck EMD Serono, Novartis, Quintiles, Sanofi Japan, Thermofisher, UCB, AstraZeneca, Exagen, GSK, Eli Lilly.

REFERENCES

  1. 1.↵
    Fava A, Petri M. Systemic lupus erythematosus: Diagnosis and clinical management [Internet]. J. Autoimmun. 2019;96(November 2018):1–13.
    OpenUrl
  2. 2.↵
    Alarcón GS et al. Systemic lupus erythematosus in three ethnic groups: II. Features predictive of disease activity early in its course. LUMINA Study Group. Lupus in minority populations, nature versus nurture. [Internet]. Arthritis Rheum. 1998;41(7):1173–80.
    OpenUrlCrossRefPubMedWeb of Science
  3. 3.
    Mok CC, Kwok RCL, Yip PSF. Effect of renal disease on the standardized mortality ratio and life expectancy of patients with systemic lupus erythematosus. Arthritis Rheum. 2013;65(8):2154–2160.
    OpenUrlCrossRefPubMedWeb of Science
  4. 4.
    Lim SS et al. The incidence and prevalence of systemic lupus erythematosus, 2002-2004: The Georgia lupus registry. Arthritis Rheumatol. 2014;66(2):357–368.
    OpenUrl
  5. 5.
    Ferucci ED et al. Prevalence and incidence of systemic lupus erythematosus in a population- based registry of American Indian and Alaska Native people, 2007-2009. Arthritis Rheumatol. 2014;66(9):2494–2502.
    OpenUrl
  6. 6.
    Dall’Era M et al. The Incidence and Prevalence of Systemic Lupus Erythematosus in San Francisco County, California: The California Lupus Surveillance Project. Arthritis Rheumatol. 2017;69(10):1996–2005.
    OpenUrl
  7. 7.
    Izmirly PM et al. The Incidence and Prevalence of Systemic Lupus Erythematosus in New York County (Manhattan), New York: The Manhattan Lupus Surveillance Program. Arthritis Rheumatol. 2017;69(10):2006–2017.
    OpenUrl
  8. 8.↵
    Somers EC et al. Population-based incidence and prevalence of systemic lupus erythematosus: The Michigan lupus epidemiology and surveillance program. Arthritis Rheumatol. 2014;66(2):369–378.
    OpenUrl
  9. 9.↵
    Hahn BH et al. American College of Rheumatology guidelines for screening, treatment, and management of lupus nephritis. Arthritis Care Res. 2012;64(6):797–808.
    OpenUrlCrossRefPubMedWeb of Science
  10. 10.↵
    Weening JJ et al. The Classification of Glomerulonephritis in Systemic Lupus Erythematosus Revisited. J. Am. Soc. Nephrol. 2004;15(2):241–250.
    OpenUrlAbstract/FREE Full Text
  11. 11.↵
    Greloni G et al. Value of repeat biopsy in lupus nephritis flares. Lupus Sci. Med. 2014;1(1):1–6.
    OpenUrl
  12. 12.↵
    De Rosa M et al. A prospective observational cohort study highlights kidney biopsy findings of lupus nephritis patients in remission who flare following withdrawal of maintenance therapy [Internet]. Kidney Int. 2018;94(4):788–794.
    OpenUrlPubMed
  13. 13.↵
    Clark MR, Trotter K, Chang A. The Pathogenesis and Therapeutic Implications of Tubulointerstitial Inflammation in Human Lupus Nephritis [Internet]. Semin. Nephrol. 2015;35(5):455–464.
    OpenUrlCrossRefPubMed
  14. 14.↵
    Arazi A et al. The immune cell landscape in kidneys of patients with lupus nephritis. Nat. Immunol. 2019;20(7):902–914.
    OpenUrlCrossRefPubMed
  15. 15.
    Der E et al. Tubular cell and keratinocyte single-cell transcriptomics applied to lupus nephritis reveal type I IFN and fibrosis relevant pathways [Internet]. Nat. Immunol. [published online ahead of print: 2019]; doi:10.1038/s41590-019-0386-1
    OpenUrlCrossRef
  16. 16.↵
    Der E et al. Single cell RNA sequencing to dissect the molecular heterogeneity in lupus nephritis [Internet]. JCI Insight 2017;2(9):1–12.
    OpenUrl
  17. 17.↵
    Ginzler EM et al. Mycophenolate mofetil or intravenous cyclophosphamide for lupus nephritis. [Internet]. N. Engl. J. Med. 2005;353(21):2219–28.
    OpenUrlCrossRefPubMedWeb of Science
  18. 18.↵
    Soliman S, Mohan C. Lupus nephritis biomarkers [Internet]. Clin. Immunol. 2017;185:10–20.
    OpenUrlPubMed
  19. 19.
    Kiani AN et al. Urinary vascular cell adhesion molecule, but not neutrophil gelatinase- associated lipocalin, is associated with lupus nephritis. J. Rheumatol. 2012;39(6):1231–1237.
    OpenUrlAbstract/FREE Full Text
  20. 20.
    Brunner HI et al. Association of noninvasively measured renal protein biomarkers with histologic features of lupus nephritis. Arthritis Rheum. 2012;64(8):2687–2697.
    OpenUrlCrossRefPubMedWeb of Science
  21. 21.↵
    Kiani AN et al. Serum osteoprotegrin (OPG) in subclinical atherosclerosis in systemic lupus erythematosus. Lupus 2017;26(8):865–870.
    OpenUrlPubMed
  22. 22.↵
    Annunziato F, Romagnani C, Romagnani S. The 3 major types of innate and adaptive cell- mediated effector immunity [Internet]. J. Allergy Clin. Immunol. 2015;135(3):626–635.
    OpenUrlCrossRefPubMed
  23. 23.↵
    Munroe ME et al. Altered type II interferon precedes autoantibody accrual and elevated type I interferon activity prior to systemic lupus erythematosus classification. [Internet]. Ann. Rheum. Dis. 2016;75(11):2014–2021.
    OpenUrlAbstract/FREE Full Text
  24. 24.↵
    Miyake K et al. Genetically determined interferon-gamma production influences the histological phenotype of lupus nephritis. [Internet]. Rheumatology (Oxford). 2002;41(5):518–24.
    OpenUrlCrossRefPubMedWeb of Science
  25. 25.↵
    Nakashima H et al. The combination of polymorphisms within interferon-γ receptor I and receptor 2 associated with the risk of systemic lupus erythematosus. FEBS Lett. 1999;453(1– 2):187–190.
    OpenUrlCrossRefPubMedWeb of Science
  26. 26.↵
    Uhm WS et al. Cytokine balance in kidney tissue from lupus nephritis patients. Rheumatology 2003;42(8):935–938.
    OpenUrlCrossRefPubMedWeb of Science
  27. 27.↵
    Masutani K et al. Predominance of Th1 Immune Response in Diffuse Proliferative Lupus Nephritis. Arthritis Rheum. 2001;44(9):2097–2106.
    OpenUrlCrossRefPubMedWeb of Science
  28. 28.↵
    Schwarting A, Wada T, Kinoshita K, Tesch G, Kelley VR. IFN-gamma receptor signaling is essential for the initiation, acceleration, and destruction of autoimmune kidney disease in MRL- Fas(lpr) mice. [Internet]. J. Immunol. 1998;161(1):494–503.
    OpenUrlAbstract/FREE Full Text
  29. 29.
    Haas C, Le Hir M, Ryffel B. IFN-γ receptor deletion prevents autoantibody production and glomerulonephritis in lupus-prone (NZB x NZW)F1 mice. J. Immunol. 1998;160(8):3713–3718.
    OpenUrlAbstract/FREE Full Text
  30. 30.↵
    Jackson SW et al. B cell IFN-γ receptor signaling promotes autoimmune germinal centers via cell-intrinsic induction of BCL-6. J. Exp. Med. 2016;213(5):733–750.
    OpenUrlAbstract/FREE Full Text
  31. 31.↵
    Pollard KM, Cauvi DM, Toomey CB, Morris K V., Kono DH. Interferon-gamma and systemic autoimmunity. Discov. Med. 2013;16(87):123–131.
    OpenUrlPubMed
  32. 32.↵
    Domeier PP et al. IFN-γ receptor and STAT1 signaling in B cells are central to spontaneous germinal center formation and autoimmunity. J. Exp. Med. 2016;213(5):715–732.
    OpenUrlAbstract/FREE Full Text
  33. 33.↵
    Takashima S et al. T cell-derived interferon-γ programs stem cell death in immune-mediated intestinal damage. [Internet]. Sci. Immunol. 2019;4(42):1–15.
    OpenUrl
  34. 34.
    van Vollenhoven RF et al. Efficacy and safety of ustekinumab, an IL-12 and IL-23 inhibitor, in patients with active systemic lupus erythematosus: results of a multicentre, double-blind, phase 2, randomised, controlled study. [Internet]. Lancet (London, England) 2018;392(10155):1330–1339.
    OpenUrlPubMed
  35. 35.
    Wallace DJ et al. Baricitinib for systemic lupus erythematosus: a double-blind, randomised, placebo-controlled, phase 2 trial [Internet]. Lancet 2018;392(10143):222–231.
    OpenUrlCrossRefPubMed
  36. 36.
    Gadina M et al. Translational and clinical advances in JAK-STAT biology: The present and future of jakinibs. [Internet]. J. Leukoc. Biol. 2018;(May):1–16.
  37. 37.↵
    Jordan J et al. 251 Type II but not type I interferon signifies clinical response to ustekinumab in patients with systemic lupus erythematosus [Internet]. In: Abstracts. Lupus Foundation of America; 2019:pA185.1-A185
  38. 38.↵
    DeBenedetti F et al. Interferon-gamma (IFN-γ) Neutralization with Emapalumab and Time to Response in Patients with Macrophage Activation Syndrome (MAS) Complicating Systemic Juvenile Idiopathic Arthritis (s-JIA) who failed High-Dose Glucocorticoids. Arthritis Rheumatol. 2019;71 (suppl.
  39. 39.
    Baechler EC et al. Interferon-inducible gene expression signature in peripheral blood cells of patients with severe lupus. [Internet]. Proc. Natl. Acad. Sci. U. S. A. 2003;100(5):2610–5.
    OpenUrlAbstract/FREE Full Text
  40. 40.
    Kirou KA et al. Coordinate overexpression of interferon-?-induced genes in systemic lupus erythematosus [Internet]. Arthritis Rheum. 2004;50(12):3958–3967.
    OpenUrlCrossRefPubMedWeb of Science
  41. 41.
    Bennett L et al. Interferon and Granulopoiesis Signatures in Systemic Lupus Erythematosus Blood [Internet]. J. Exp. Med. 2003;197(6):711–723.
    OpenUrlAbstract/FREE Full Text
  42. 42.↵
    Hall JC et al. Precise probes of type II interferon activity define the origin of interferon signatures in target tissues in rheumatic diseases. Proc. Natl. Acad. Sci. 2012;109(43):17609–17614.
    OpenUrlAbstract/FREE Full Text
  43. 43.↵
    Rusinova I et al. INTERFEROME v2.0: An updated database of annotated interferon- regulated genes. Nucleic Acids Res. 2013;41(D1):1040–1046.
    OpenUrl
  44. 44.↵
    Catalina MD, Bachali P, Geraci NS, Grammer AC, Lipsky PE. Gene expression analysis delineates the potential roles of multiple interferons in systemic lupus erythematosus [Internet]. Commun. Biol. 2019;2(1). doi:10.1038/s42003-019-0382-x
    OpenUrlCrossRef
  45. 45.↵
    Klocke J et al. Mapping urinary chemokines in human lupus nephritis: Potentially redundant pathways recruit CD4+ and CD8+ T cells and macrophages. Eur. J. Immunol. 2017;47(1):180–192.
    OpenUrl
  46. 46.↵
    Hashimshony T et al. CEL-Seq2: Sensitive highly-multiplexed single-cell RNA-Seq [Internet]. Genome Biol. 2016;17(1):1–7.
    OpenUrlCrossRefPubMed
  47. 47.↵
    Subramanian A et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. U. S. A. 2005;102(43):15545–15550.
    OpenUrlAbstract/FREE Full Text
  48. 48.↵
    Ontology CTG et al. The Gene Ontology Consortium, Michael Ashburner1, Catherine A. Ball3, Judith A. Blake4, David Botstein3, Heather Butler1, J. Michael Cherry3, Allan P. Davis4, Kara Dolinski3, Selina S. Dwight3, Janan T. Eppig4, Midori A. Harris3, David P. Hill4, Laurie Is. Nat. Genet. 2000;25(1):25–29.
    OpenUrlCrossRefPubMedWeb of Science
  49. 49.↵
    Carbon S et al. The Gene Ontology Resource: 20 years and still GOing strong. Nucleic Acids Res. 2019;47(D1):D330–D338.
    OpenUrlCrossRefPubMed
  50. 50.↵
    Sergushichev AA. An algorithm for fast preranked gene set enrichment analysis using cumulative statistic calculation [Internet]. bioRxiv 2016;060012.
Back to top
PreviousNext
Posted March 17, 2020.
Download PDF

Supplementary Material

Data/Code
Email

Thank you for your interest in spreading the word about medRxiv.

NOTE: Your email address is requested solely to identify you as the sender of this article.

Enter multiple addresses on separate lines or separate them with commas.
Integration of urine proteomics and renal single-cell genomics identifies an interferon-gamma response gradient in lupus nephritis
(Your Name) has forwarded a page to you from medRxiv
(Your Name) thought you would like to see this page from the medRxiv website.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Integration of urine proteomics and renal single-cell genomics identifies an interferon-gamma response gradient in lupus nephritis
Andrea Fava, Jill Buyon, Chandra Mohan, Ting Zhang, H. Michael Belmont, Peter Izmirly, Robert Clancy, Jose Monroy Trujillo, Derek Fine, Yuji Zhang, Laurence Magder, Deepak A. Rao, Arnon Arazi, Celine C. Berthier, Anne Davidson, Betty Diamond, Nir Hacohen, David Wofsy, the Accelerating Medicines Partnership in SLE network, Soumya Raychaudhuri, Michelle Petri
medRxiv 2020.03.13.20034348; doi: https://doi.org/10.1101/2020.03.13.20034348
Digg logo Reddit logo Twitter logo Facebook logo Google logo LinkedIn logo Mendeley logo
Citation Tools
Integration of urine proteomics and renal single-cell genomics identifies an interferon-gamma response gradient in lupus nephritis
Andrea Fava, Jill Buyon, Chandra Mohan, Ting Zhang, H. Michael Belmont, Peter Izmirly, Robert Clancy, Jose Monroy Trujillo, Derek Fine, Yuji Zhang, Laurence Magder, Deepak A. Rao, Arnon Arazi, Celine C. Berthier, Anne Davidson, Betty Diamond, Nir Hacohen, David Wofsy, the Accelerating Medicines Partnership in SLE network, Soumya Raychaudhuri, Michelle Petri
medRxiv 2020.03.13.20034348; doi: https://doi.org/10.1101/2020.03.13.20034348

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Subject Area

  • Rheumatology
Subject Areas
All Articles
  • Addiction Medicine (160)
  • Allergy and Immunology (413)
  • Anesthesia (90)
  • Cardiovascular Medicine (855)
  • Dentistry and Oral Medicine (159)
  • Dermatology (97)
  • Emergency Medicine (248)
  • Endocrinology (including Diabetes Mellitus and Metabolic Disease) (393)
  • Epidemiology (8550)
  • Forensic Medicine (4)
  • Gastroenterology (383)
  • Genetic and Genomic Medicine (1744)
  • Geriatric Medicine (167)
  • Health Economics (371)
  • Health Informatics (1237)
  • Health Policy (619)
  • Health Systems and Quality Improvement (467)
  • Hematology (196)
  • HIV/AIDS (372)
  • Infectious Diseases (except HIV/AIDS) (10281)
  • Intensive Care and Critical Care Medicine (552)
  • Medical Education (192)
  • Medical Ethics (51)
  • Nephrology (210)
  • Neurology (1670)
  • Nursing (97)
  • Nutrition (248)
  • Obstetrics and Gynecology (325)
  • Occupational and Environmental Health (450)
  • Oncology (926)
  • Ophthalmology (263)
  • Orthopedics (100)
  • Otolaryngology (172)
  • Pain Medicine (112)
  • Palliative Medicine (40)
  • Pathology (252)
  • Pediatrics (534)
  • Pharmacology and Therapeutics (247)
  • Primary Care Research (207)
  • Psychiatry and Clinical Psychology (1764)
  • Public and Global Health (3832)
  • Radiology and Imaging (622)
  • Rehabilitation Medicine and Physical Therapy (319)
  • Respiratory Medicine (520)
  • Rheumatology (207)
  • Sexual and Reproductive Health (165)
  • Sports Medicine (157)
  • Surgery (190)
  • Toxicology (36)
  • Transplantation (101)
  • Urology (74)