key: cord-0910332-a3e9cy4v
authors: Lin, Hongchun; Ma, Xinxin; Xiao, Fang; Su, Hua; Shi, Yaling; Liu, Yuntao; Song, Lan; Zhang, Zhongde; Zhang, Chun; Peng, Hui
title: Identification of a special cell type as a determinant of the kidney tropism of SARS‐CoV‐2
date: 2021-07-29
journal: FEBS J
DOI: 10.1111/febs.16114
sha: 05a64972d0b960c2edea87d50f26d82fdeae86ee
doc_id: 910332
cord_uid: a3e9cy4v

The kidney tropism of severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) has been well‐validated clinically and often leads to various forms of renal damage in coronavirus disease‐2019 (COVID‐19) patients. However, the underlying mechanisms and diagnostic approaches remain to be determined. We interrogated the expression of virus‐related host factors in single‐cell RNA sequencing (scRNA‐seq) datasets of normal human kidneys and kidneys with pre‐existing diseases and validated the results with urinary proteomics of COVID‐19 patients and healthy individuals. We also assessed the effects of genetic variants on kidney susceptibility using expression quantitative trait loci (eQTLs) databases. We identified a subtype of tubular cells, which we named PT‐3 cells, as being vulnerable to SARS‐CoV‐2 infections in the kidneys. PT‐3 cells were enriched in viral entry factors and replication and assembly machinery but lacked antiviral restriction factors. Immunohistochemistry confirmed positive staining of PT‐3 cell marker SCL36A2 on kidney sections from COVID‐19 patients. Urinary proteomic analyses of COVID‐19 patients revealed that markers of PT‐3 cells were significantly increased, along with elevated viral receptor angiotensin‐converting enzyme 2. We further found that the proportion of PT‐3 cells increased in diabetic nephropathy but decreased in kidney allografts and lupus nephropathy, suggesting that kidney susceptibility varied among these diseases. We finally identified several eQTLs that regulate the expression of host factors in kidney cells. PT‐3 cells may represent a key determinant for the kidney tropism of SARS‐CoV‐2, and detection of PT‐3 cells may be used to assess the risk of renal infection during COVID‐19.

Hongchun Lin, Xinxin Ma and Fang Xiao contributed equally to this work.

The kidney tropism of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been well-validated clinically and often leads to various forms of renal damage in coronavirus disease-2019 (COVID- 19) patients. However, the underlying mechanisms and diagnostic approaches remain to be determined. We interrogated the expression of virus-related host factors in single-cell RNA sequencing (scRNA-seq) datasets of normal human kidneys and kidneys with pre-existing diseases and validated the results with urinary proteomics of COVID-19 patients and healthy individuals. We also assessed the effects of genetic variants on kidney susceptibility using expression quantitative trait loci (eQTLs) databases. We identified a subtype of tubular cells, which we named PT-3 cells, as being vulnerable to SARS-CoV-2 infections in the kidneys. PT-3 cells were enriched in viral entry factors and replication and assembly machinery but lacked antiviral restriction factors. Immunohistochemistry confirmed positive staining of PT-3 cell marker SCL36A2 on kidney sections from COVID-19 patients. Urinary proteomic analyses of COVID-19 patients revealed that markers of PT-3 cells were significantly increased, along with elevated viral receptor angiotensin-converting enzyme 2. We further found that the proportion of PT-3 cells increased in diabetic nephropathy but decreased in kidney allografts and lupus nephropathy, suggesting that kidney susceptibility varied among these diseases. We finally identified several eQTLs that regulate the expression of host factors in kidney cells. PT-3 cells may represent a key determinant for the kidney tropism of SARS-CoV-2, and detection of PT-3 cells may be used to assess the risk of renal infection during COVID-19.

Coronavirus disease-2019 (COVID- 19) is an infectious disease caused by a novel discovered coronavirus, the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). It has been reported that the kidney is one of the most commonly affected extrapulmonary organs of SARS-CoV-2 infection and varying degrees of renal damage is often observed in COVID-19 patients [1, 2] . In addition to indirect factors (e.g., cytokine storm, hypoxia), the direct kidney infection of SARS-CoV-2 could also lead to kidney injury [3] [4] [5] [6] . Autopsy reports have provided evidence that the virus could directly infect kidney cells of COVID-19 patients and patient-derived SARS-CoV-2 were able to replicate in nonhuman primate kidney tubular epithelial cells [7] . However, why kidney cells are permissive to SARS-CoV-2 infection and which renal cell types are more susceptible remain poorly understood.

The successful infection of SARS-CoV-2 depends critically on its interaction with target cells, which involves various host factors such as cell surface receptors, cellular proteases, antiviral restriction factors, factors participate viral replication, and assembly [8] [9] [10] . Angiotensin-converting enzyme 2 (ACE2) has been identified as the SARS-CoV-2 receptor, and the transmembrane serine protease 2 (TMPRSS2) has been confirmed to cleave the viral spike (S) protein, which facilitates viral entry [11] . However, several candidate receptors and proteases may serve similar functions. For example, basign (BSG, also known as CD147), dipeptidyl peptidase 4 (DPP4), and alanyl aminopeptidase (ANPEP) reportedly interact with the S protein to promote SARS-CoV-2 or other human coronavirus entry into cells [12] [13] [14] . Other members of the TMPRSS family, such as TMPRSS4, may function similarly to TMPRSS2 [15] . In addition, cellular proteases, including cathepsins B and L (CTSB/L) and furin, can mediate S protein cleavage [11, 16] . In contrast, lymphocyte antigen 6 family member E (LY6E) and interferon-induced transmembrane proteins had been identified as potent coronavirus restriction factors [9, 17] . Thus, comprehensive characterizations of these host factors in kidney cells will greatly improve our understanding on the renal tropism and pathogenicity of SARS-CoV-2.

Single-cell RNA sequencing (scRNA-seq) analysis is a powerful approach to characterize complex tissues at single-cell resolution and reveal diverse cellular responses in different cell types. Many studies have applied scRNA-seq to predict the tissue tropism of SARS-CoV-2, including the kidney. The transcripts of ACE2 were primarily detected in proximal tubular epithelial cells (PTECs), while the mRNA of TMPRSS2 was robustly expressed in the distal convoluted tubule (DCT) and intercalated cells rather than the PTECs [18, 19] . Few studies have comprehensively evaluated the expression patterns of other host factors (such as restriction factors and replication/assembly factors) in different kidney cells.

The expression of host factors determines SARS-CoV-2 susceptibility, which is partially regulated by genetic variations, especially expression quantitative trait loci (eQTLs) [20] . eQTLs refer to genetic variants that regulate transcript expression levels. Although several studies have described eQTLs associated with the expression of ACE2 in multiple tissues [20, 21] , little is known regarding the association of eQTLs with other host factors that underlie the SARS-CoV-2 kidney tropism.

In this study, we surveyed the expression patterns of host factors in published scRNA-seq datasets of healthy human kidneys and clarified the distribution of ACE2 and TMPRSS2 in PTECs. We also identified and characterized a PTEC subtype that was highly susceptible to SARS-CoV-2 infection, which we named PT-3 cells. Due to the limitation of scRNA-seq data, such as dropout events, and possible differences between the transcriptome and the proteome, urinary proteomics was applied to validate the results obtained via single-cell analysis. We confirmed that the marker proteins of PT-3 cells were significantly elevated in the urine of COVID-19 patients. Moreover, we assessed the expression pattern of host factors in scRNA-seq datasets of kidneys with kidney diseases, to predict the tropism of SARS-CoV-2. Finally, we explored the eQTLs regulating the expression of host factors and identified the eQTLs associated with host proteases and restriction factors, providing a potential interpretation for the differences in kidney susceptibility among individuals.

Identification of a specific SARS-CoV-2susceptible cell type in the kidney We first evaluated the autopsy kidney specimens from 26 COVID-19 patients and found that spherical particles surrounded by crown-like projections, which were highly suspected as SARS-CoV-2 virus, scattered within the apical side of renal PTECs by electron microscopy (Fig. 1A) . This finding led us to the question that there could be differences in susceptibility between distinct kidney cell types to SARS-CoV-2 infection. To answer this question, we systematically surveyed the expression patterns of host factors required for the SARS-CoV-2 replication cycle in publicly available scRNA-seq datasets of normal human kidneys (GSE131882) [22] , which included 15 728 cells from three normal human kidney samples. Our analysis identified 11 distinct kidney cell types (Fig. 1B) , in concordance with the original report. Although ACE2 mRNA was predominantly expressed in PTECs, it was undetectable in other cell types, such as podocytes ACE2  BSG  CLEC4M  DPP4  ANPEP  TMPRSS2  TMPRSS4   CTSB  CTSL  FURIN  IFITM2  IFITM3  LY6E  ZCRB1  TOP3B  RAB10  RAB1A  RAB14  TAPT1  AP2A2  AP2M1 (PODO) and endothelial cells. In contrast, the expression of BSG, a potential SARS-CoV-2 receptor, was broadly detected in kidney cells, with the highest expression levels observed in collecting duct (CD) cells (Fig. 1B) . We also surveyed the receptors of other known coronaviruses, such as C-type lectin domain family 4 member M (CLEC4M; SARS-CoV) [23] , DPP4 (Middle East respiratory syndrome coronavirus), and ANPEP (human coronavirus 229E) and found that the expression of these proteins was predominantly observed on PTECs, except for CLEC4M (Fig. 1B) . When examining cellular proteases, we found that CTSB was broadly expressed in different cell types, including PTECs; in contrast, TMPRSS2 was strongly expressed in CD and DCT cells, rather than PTECs. Other proteases, such as CTSL and furin, could also be detected in various kidney cells but were almost undetectable in PTECs (Fig. 1B) . These results indicated that PTECs were the only kidney cell type that showed appreciable ACE2 mRNA levels coexpressed with varying levels of proteases. We next investigated the gene expression patterns of restriction factors. As shown in Fig. 1B , LY6E and IFITM2/3 were highly expressed in glomerular cells, while only minimal transcripts for these genes were detected in renal tubular cells, especially in PTECs. Finally, we surveyed the expression patterns of factors involved in viral replication and assembly. The factors required for viral replication (ZCRB1, zinc finger CCHC-type and RNA-binding motif-containing 1; TOP3B, DNA topoisomerase III beta) could be detected in most kidney cell types. In addition, the factors engaged in viral assembly and trafficking were detected across the whole dataset, with moderate to high levels of expression (Fig. 1b) .

Since only a proportion of PTECs are ACE2positive, we speculated that there might be some subtypes of PTECs that were more vulnerable to SARS-CoV-2 infection. Hence, we performed an unsupervised clustering analysis of 6071 PTECs from the same dataset and obtained seven subclusters of PTECs (PT-1 to PT-7; Fig. 1C ). We then delineated the distribution of double-positive cells for various combinations of ACE2 with TMPRSS2, CTSB, and CTSL in PTECs (Fig. 1D) . Intriguingly, we indeed identified that only the PT-3 subcluster demonstrated a higher proportion of ACE2 + /CTSB + double-positive cells compared with other PTECs (20.54% vs 2.24%; Fig. 1E ). In addition, the ACE2 + /TMPRSS2 + and ACE2 + /CTSL + doublepositive cells also primarily appear among PT-3 cells, with lower proportions than the ACE2 + /CTSB + cells (Fig. 1D,E) . Then, we characterized the expression of host factors in these PTEC subclusters ( Fig. 2A) .

Overall, PT-3 cells showed the highest expression of AEC2 and co-expressed with two other human coronavirus receptors, DPP4 and ANPEP. Compared with other PTEC subtypes, PT-3 cells exhibited higher levels of proteases (e.g., CTSB) and cellular machinery required for viral replication (ZCRB1, TOP3B) and assembly, but lacked viral restriction factors (LY6E, IFITM2, and IFITM3; Fig. 2A ). Meanwhile, gene set enrichment analysis (GSEA) revealed that PT-3 cells were not only significantly enriched in various pathways related to virus infection, but also enriched in endocytosis and autophagy pathways (Fig. 2B ). The cell identities and distribution of proximal tubule segment-specific markers [24] in these PTECs were described in Fig. 2C ,D. PT-3 cells primarily expressed specific markers for Segment 3, and SLC36A2, PEPD were identified as the signature genes of PT-3 cells.

Considering the optimal environment for coronavirus replication provided by PT-3 cells, it strongly indicated that these cells were the target cells of SARS-CoV-2 infection in the kidney.

To verify our findings derived from scRNA-seq data analysis, we first performed immunohistochemistry (IHC) staining for PT-3 cell marker SCL36A2 on kidney sections from COVID-19 patients and controls. In control kidney, SLC36A2 was detected in some of the proximal tubules. This focal or scattered staining pattern indicated that PT-3 cells were dispersed among the proximal tubules. In kidney from COVID-19 patients, light microscopy showed decrease in SLC36A2-positive cells in the proximal tubules which had typical features of acute kidney injury (AKI), including loss of brush border, flattening, and focal loss of tubular epithelial cells (Fig. 3A) .

We next performed urinary proteomic analysis and focused on protein abundance changes in various host factors, as well as that of specific markers for kidney cells. Urine samples were collected from 97 COVID-19 patients and 350 age-and gender-matched healthy individuals ( Table 1 ). The median of proteins identified were 1370 and 1024 in healthy individuals and COVID-19 patients, respectively [false discovery rate (FDR) < 1%]. We observed a significantly increased level of ACE2 in the urine of COVID-19 patients, compared with healthy individuals (Fig. 3B ). Since ACE2 is predominantly expressed in PTECs in kidney, we speculated that the increased ACE2 protein abundance may be due to the shedding of PTECs. The level of CTSB was also significantly increased. In contrast, TMPRSS2 and CTSL appeared unchanged. Moreover, the protein abundance of the restriction factors IFITM1 and IFITM3 was increased in COVID-19 group compared with the healthy group. And the protein levels of cellular factors involved in virus assembly, such as AP2M1, RHOA, and RAB14, were also increased in COVID-19 patients (Fig. 3B ).

As for cell type-specific markers (Fig. 3C ), the levels of PTEC marker LRP2 and ALDOB in the urine of COVID-19 patients were significantly higher than that of healthy individuals. Whereas markers associated with complement factor H cells, PODOs, and mesangial cells (MESs) decreased, markers for loop of Henle cells, DCT cells, and CD cells remained unchanged. No markers for endothelial cells were detected in either group. Importantly, a significantly elevated level of PT-3 cell-specific markers (PEPD, SLC36A2) was also detected (Fig. 3D) . Thus, our urinary proteomic results provided supporting evidence that PT-3 cells represent the preferred targeting cell type for SARS-CoV-2 in the kidney.

Although COVID-19 patients with comorbidities, especially those with chronic kidney disease (CKD) are at higher risk of kidney injury, it is still unclear whether underlying kidney disease, per se, increases the vulnerability of kidney cells to SARS-CoV-2 infection. To CLEC4M  DPP4  TMPRSS2  TMPRSS4   CTSB  CTSL  FURIN  IFITM2  IFITM3  LY6E  ZCRB1  TOP3B  RAB10  RAB1A  RAB14  TAPT1  AP2A2  AP2M1  CHMP2A   -1  0  1   20  40  60  10  20  30  40  5  10  15  20  2  3  4  5  6  20  40  60 Percent Expressed

Receptor Protease Restriction Replication Assembly

Average expression Assembly   ACE2  BSG  DPP4  ANPEP  TMPRSS2  CTSL  CTSB  FURIN  IFITM1  IFITM3  AP2A2  AP2M1  RHOA  RAB1A  RAB14 RAB10 CHMP2A predict kidney susceptibility to SARS-CoV-2 infection, we assessed changes in the expression patterns of host factors and proportions of PT-3 cells in kidneys with pre-existing kidney diseases, including diabetic nephropathy (DN), kidney allografts, and lupus nephropathy (LN; Table 2 ). We first analyzed the scRNA-seq dataset (GSE131882), which included 9883 cells from 3 DN biopsy samples and 15 728 cells from three normal human kidney samples [22] . The expression of ACE2 was predominantly detected in PTECs, and its level was significantly increased in diabetic kidneys (Fig. 4A,B) . This result was in line with previous studies at the mRNA and protein levels [25, 26] . Compared with controls, the expression levels of protease CTSB and factors related to virus assembly were also significantly increased in diabetic PTECs (Fig. 4A,B) . In diabetic kidneys, although the overall number of PTECs decreased, the proportion of PT-3 cells was significantly higher compared with controls (28.5% vs 11.1%; Fig. 4C ). Concomitantly, the proportions of ACE2 + /CTSB + , ACE2 + /TMPRSS2 + , and ACE2 + / CTSL + cells among PTECs were higher than those in controls (Fig. 4D) . We also surveyed changes in the expression levels of host factors in PT-3 cells (Fig. 4E) . There was no upregulation of ACE2 mRNA in diabetic PT-3 cells, although a significant increase in DPP4 expression was observed. Except furin, the expression of major cellular proteases was also upregulated in diabetic PT-3 cells. Compared with control kidneys, the expression levels of IFITM2 and IFITM3 in PT-3 cells from diabetic kidneys were increased, while LY6E remained unchanged (Fig. 4E) .

In summary, the proportion of PT-3 cells in diabetic kidney was increased, which was accompanied by upregulation of protease expression. These findings indicated a higher risk of SARS-CoV-2 entry into the kidneys of DN patients, which consequently increases the incidence of kidney infection.

We next interrogated a public scRNA-seq dataset GSE140989, which included 4690 cells from five kidney allograft biopsy samples and 17 578 cells from 19 normal kidney tissue samples [27] . In PTECs from kidney allografts, the proportions of cells expressing receptors (ACE2, DPP4, and ANPEP) and cells expressing CTSB were lower than that in healthy kidneys. Compared with healthy kidneys, restriction factors, especially LY6E, in kidney allografts were increased, while there was no difference in the expression of factors involved in virus assembly (Fig. 5A) . Then, we found the fraction of PT-3 cells in kidney allografts was significantly decreased (Fig. 5B) , and the expression of ACE2 was downregulated in these PT-3 cells, while proteases CTSB/CTSL and restriction factor LY6E were upregulated (Fig. 5C) .

We further re-analyzed a scRNA-seq dataset for LN (Immport, SDY997), which included 2629 cells from 21 LN kidney biopsy samples and 160 cells from three healthy donor kidney samples [28] . We observed that both ACE2 expression and the proportion of ACE2 + cells decreased in LN PTECs. Compared with controls, we found a higher proportion of cells expressing CTSB, restriction factors, and factors for viral replication (ZCRB1) and assembly in LN PTECs (Fig. 5D) . We also found that PT-3 cells decreased slightly in LN kidneys compared with healthy kidneys (28.6% vs ACE2  BSG  CLEC4M  DPP4  ANPEP  TMPRSS2  TMPRSS4   CTSB  CTSL  FURIN  IFITM2  IFITM3  LY6E  ZCRB1  TOP3B  RAB10 RAB1A RAB14 TAPT1 AP2M1 CHMP2A   10  20  30  40  50  20  40  60  80  10  20  30  40 Fig. 5E ). And in these lupus PT-3 cells, ACE2 mRNA level decreased while the levels of proteases and restriction factors increased (Fig. 5F) .

Overall, the decreased proportion of PT-3 cells in kidney allografts and LN suggests that kidney susceptibility varied among different kidney diseases.

Since a considerable proportion of COVID-19 patients present no evidence of kidney infection, even among those with multiple organ infections, suggesting the effects of genetic variants on kidney susceptibility. To delineate this genetic basis, we examined eQTLs of genes encoding host factors, using the Human Kidney eQTL Atlas [29] and GTEx database. We did not identify any significant eQTL for ACE2 gene in the kidney (Table S1 ), consistent with previous study [21] . However, we did find some eQTLs associated with the kidney expression of proteases and restriction factors, such as IFITM2 and LY6E. Among these, five eQTLsingle nucleotide polymorphisms (SNPs; rs3809111, rs111797074, rs75117940, rs909097, and rs909098) were significantly associated with the low expression level of IFITM2. And these associations were only detected in the tubules, suggesting individual differences in IFITM2 expression in renal tubules. In kidney cortex samples from GTEx database, another significant eQTL rs111817998 was found to contribute to high IFITM2 expression.

Among the six eQTLs identified, three of them (rs3809111, rs75117940, and rs111817998) located in the promoter region, and two located in intronic regions (rs909097 and rs909098) (Fig. 6A) . A single nucleotide change in transcription factor-binding sites (TFBSs) can affect interactions with transcription factors, altering downstream gene expression. To explore the regulatory mechanisms associated with these eQTL-SNPs, we examined the potential TFBSs in the promoter region of IFITM2, using PROMO [30] . The (Fig. 6B ). STAT4 and TFII-I have been reported to participate in the host innate immune response against viral infection [31, 32] . Although the association and regulatory mechanism need to be further confirmed, special attention should be paid to these genetic variants in future studies of COVID-19 patients.

Although the high prevalence of kidney involvement is well-documented in patients with COVID-19, limited information exists regarding host factors that determine this kidney tropism of SARS-CoV-2. In the current study, we comprehensively characterized the expression patterns and expression levels of host factors in kidney cells, by mining human kidney scRNAseq datasets, and identified PT-3 cells, a special subtype of PTECs that may function as a "gateway" cell type for the SARS-CoV-2 invasion of the kidney. Our analysis revealed that PT-3 cells not only enriched for ACE2 expression, but also expressed two other human coronavirus receptors, DPP4 and ANPEP, at high levels than other kidney cells. Additionally, the expression levels of proteases and host factors involved in virus replication and assembly in PT-3 cells were higher than other PTECs. PT-3 cells were also significantly enriched in various virus infection-related pathways, as well as endocytosis and autophagy pathways. These pathways play important roles in viral infection [8, 33] . Although the expression of restriction factors was also observed in PT-3 cells, the basal expression level was lower than that in other PTECs. Based on the transcription signature of PT-3 cells, we predict that these cells are highly permissive to infection by SARS-CoV-2 and other human coronaviruses. Urinary proteomics confirmed that the protein levels of PT-3 markers, in addition to ACE2, CTSB, and restriction factors, were significantly increased in the urine of COVID-19 patients. This result further supported the findings derived from our scRNA-seq analyses, which suggested that PT-3 cells were the target cells of SARS-CoV-2 and provided new insights to understand the pathogenesis of kidney infections. Moreover, the detection of PT-3 cell-specific markers in the urine could be utilized as a noninvasive and dynamic method for kidney infection risk assessment. TMPRSS2 transcripts were most abundantly expressed in CD cells and DCT cells and rarely expressed in the PTECs that highly expressed ACE2. This expression pattern is distinct from the pattern observed in nasal epithelial cells, where the protease TMPRSS2 is highly expressed on the cell membranes of goblet cells and ciliated cells and is co-expressed with ACE2 [19] . In contrast, the proportion of ACE2 and TMPRSS2 co-expressing cells was much lower than that of ACE2 and CTSB co-expressing cells among PT-3 cells. For both SARS-CoV and SARS-CoV-2, S protein priming by TMPRSS2 has been shown to be crucial for entry into target cells [11] . Compared with SARS-CoV, the S protein of SARS- CoV-2 contains a four-amino acid insertion at the S1/ S2 site, which allows it to be proteolytically cleaved by a variety of cellular proteases, including CTSB and CTSL [34] . Previous studies also verified that CTSB had high proteolytic activity against SARS-CoV-2, whereas CTSL showed lower proteolytic activity against SARS-CoV-2 than SARS-CoV [34] . These findings indicated that SARS-CoV-2 might use both TMPRSS2 and CTSB/L for S protein priming. Thus, the prominent co-expression of ACE2 and CTSB in PT-3 cells suggests that the invasion of SARS-CoV-2 into these cells is likely mediated by CTSB and that the routes of viral entry may be cell type-specific. However, the specific mechanism requires further study.

Our survey identified an increased proportion of PT-3 cells in diabetic kidneys, which could lead to enhanced kidney susceptibility to SARS-CoV-2 infection. This conclusion is consistent with a previous study and clinical observations [2, 3, 25] . We also found lower proportion of PT-3 cells in kidney allografts than healthy controls, accompanied by the upregulation of restriction factors. However, it has been reported that kidney transplant recipients with COVID-19 are at high risk of acute kidney injury [35] . The kidney abnormalities in kidney transplant patients could be due to other mechanisms such as immunemediated effects. We also predicted the decreased viral susceptibility of kidneys with LN, based on the decreased proportion of PT-3 cells observed in LN kidneys compared with healthy kidneys. Whether patients with systemic lupus erythematosus have a higher risk of SARS-CoV-2 infection than the general population remains unclear [36] . These findings provided useful information regarding the SARS-CoV-2 infection of kidneys in LN patients. In addition, evidence on the kidney infectivity of SARS-CoV-2 virus has been largely derived from autopsy studies [4, 7, 37, 38] . Some studies also reported the detection of SARS-CoV-2 RNA or protein in the urine of COVID-19 patients [39, 40] . Despite these positive findings, many other studies did not find the presence of SARS-CoV-2 virus in the kidneys of COVID-19 patients [41] . Therefore, caution is needed to draw conclusions about the direct kidney infection by SARS-CoV-2 [38] .

Indeed, there are discrepancies between transcriptome and proteome, and this issue could be more significant in single-cell transcriptomics due to shallow sequencing depth and low mRNA capture efficiency of this approach. For instance, at single-cell transcriptomics level, the mRNA levels of ACE2 in major glomeruli cell types such as PODOs and MESs were hardly detected. While at protein level, low level of ACE2 was confirmed in the glomeruli [26] . Regarding TMPRSS2, it was hardly detected in the ACE2positive PTECs based on the results of single-cell analysis, whereas immunostaining confirmed that it also colocalized with ACE2 in the proximal tubules of kidney organoids [42] . Meanwhile, posttranscriptional modification is an important factor contributing to this discrepancy [43] . In this context, accurate prediction based on results from scRNA-seq analysis alone could be challenging, and a multi-omics data integrative strategy can provide a more definitive conclusion.

Overall, our analysis identified PT-3 cells, a PTEC subtype, as key determinant for SARS-CoV-2 infection in the kidney. These findings provide new insights into the mechanism underlying the SARS-CoV-2 infection of kidney cells and present a potential strategy for risk assessment of renal infection during COVID-19.

We first collected autopsy kidney specimens from 26 cases of patients died from COVID-19. We also collected the urine samples of 97 COVID-19 patients. Laboratory confirmed cases of COVID-19 were defined as a positive result on viral gene sequencing highly similar with known SARS-CoV-2 or real-time polymerase chain reaction assay of respiratory tract specimens or blood specimens. The control group consisted of 350 healthy individuals who were matched on age, gender with the 97 COVID-19 patients.

Autopsies were performed in the Union Hospital of Tongji Medical College of Huazhong University of Science and Technology. The study was approved by the local ethics committee. Urinary proteomics study of healthy individuals and COVID-19 patients was approved by institutional review boards of Beijing Proteome Research Center and ethics committee of Guangdong Provincial Hospital of Chinese Medicine (BF2020-052-01). All clinical investigations were conducted according to the Declaration of Helsinki. Written informed consent was obtained from enrolled subjects or next of kin (for deceased subjects).

Kidney tissue specimens were fixed in 2.5% glutaraldehyde or 10% formalin for 48-72 h before transmission electron microscope (TEM) examination or histological staining. For evidence of SARS-CoV-2 infection in kidney samples, we performed electron microscope examination. After the osmium tetroxide postfixation and gradient dehydration, the samples were embedded in Epon for TEM (HT-7800; Hitachi, Tokyo, Japan) observation, as previously described [5] .

To determine the distribution and localization of SLC36A2 in kidney, immunohistochemical staining was performed in autopsied kidney specimens obtained from four COVID-19 patients and four healthy controls, using SLC36A2 antibody from NOVUS (Catalog Number: NBP1-59418; Centennial, CO, USA).

The midstream of the first morning urine was collected and put in 56°C water bath for 30 min. The urine samples were then prepared, sterilized, and stored at À80°C. Briefly, one milliliter urine sample was centrifuged at 176 000 g for 70 min to save the pellets. The samples were then processed according to a published streamlined workflow with minor modification, in which the centrifuged pellets were digested in solution with 1 lg of trypsin at 37°C for 4 h [44] .

Tryptic peptides were separated using a self-packed capillary column packed with C18 reverse-phase particles, and the treated samples were analyzed with Thermo Fisher Orbitrap mass spectrometer (MS) coupled with online Easy-nLC 1000 nano (Nano HPLC) system (Thermo Fisher Scientific, Waltham, MA, USA). Raw LC-MS/MS files were uploaded onto FIRMIANA platform for further processing and analysis, and the MASCOT v.2.5.1 (Matrix Science, Boston, MA, USA) was used for peptide identification with peptide FDR < 1%.

We used the intensity-based absolute quantification (iBAQ) algorithm for protein quantification. For batch-tobatch comparisons, iBAQ was further transformed into iFOT (intensity-based Fraction of Total), which represented the normalized intensity of the protein identified in the LC-MS/MS analysis. For visualization purposes, the value of iFOTs is multiplied by 10 5 . Trypsin digests of human HEK293T cells were used as control quality samples and routinely assessed by the LC-MS/MS to ensure instrument reproducibility.

Mann-Whitney U-test was performed to compare protein abundance between two groups.

We obtained publicly available scRNA-seq/single-nucleus RNA sequencing (snRNA-seq) datasets of human kidney with CKDs from Gene Expression Omnibus database or ImmPort database (Table 2) . Additional data related to this paper may be requested from the authors.

The original gene-barcode matrix was downloaded and processed by Seurat3 R package (https://satijalab.org/seurat/). Cells/nuclei with < 500 or more than 5000 expressed genes were filtered out, and genes expressed in < 3 cells or nuclei were removed [22] . Then, a log-normalize step was performed using the method "LogNormalize."

We used the top 2000 most variable genes to generate clusters at multiple resolutions and selected the best matching clusters according to the instruction in Seurat [45] . Then, we performed principal component analysis for dimensionality reduction and the Louvain method for detection of cell subgroups. The first 10 principal components were selected and used for the Uniform Manifold Approximation and Projection (UMAP). Specific marker genes of different cell types were identified by likelihood ratio test in Seurat package. PTECs which enriched in the expression of ACE2 were extracted for further clustering analysis, and PTEC subclusters were identified using the FindCluster function in Seurat with the resolution parameter of 0.3. DEGs within clusters were identified using FindMarkers function with logfc.threshold = 0.25 and min.pct = 0.25.

Gene set enrichment analysis was conducted using clus-terProfiler package (v3.5) with default parameters [46] . Differential expressed genes in PT-3 subtype (compared to other PTEC subtypes) were identified with the limma R package using a log-fold-change threshold of 1.25 and P value threshold of 0.05. The preranked gene list according to fold change values was input for GSEA on the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. P value < 0.05 and adjusted P value < 0.05 were set up as the significant cutoff for GSEA.

We inquired the eQTLs which might regulate virus-related host factor gene expression in the Human Kidney eQTL Atlas (http://susztaklab.com/eqtl) and Genotype-Tissue Expression (GTEx) database.

All statistical analyses were conducted using R or PYTHON (https://www.python.org/). For comparison between two independent groups, we performed two-tailed, unpaired Student's t-test for data that were normally distributed. Otherwise, we performed Mann-Whitney U-test. P < 0.05 was considered statistically significant. The level of significance was assigned as *P < 0.05, **P < 0.01, ***P < 0.001, or ****P < 0.0001.

Additional Supporting Information may be found online in the Supporting Information section at the end of the article. Table S1 . eQTLs associated with gene expression of host factors in kidney.

Histopathology and ultrastructural findings of fatal COVID-19 infections in Washington State: a case series

Presenting characteristics, comorbidities, and outcomes among 5700 patients hospitalized with COVID-19 in the

Renal involvement and early prognosis in patients with COVID-19 pneumonia

Multiorgan and renal tropism of SARS-CoV-2

Renal histopathological analysis of 26 postmortem findings of patients with COVID-19 in China

Ultrastructural evidence for direct renal infection with SARS-CoV-2

SARS-CoV-2 renal tropism associates with acute kidney injury

Human coronavirus: hostpathogen interaction

Distinct patterns of IFITMmediated restriction of filoviruses, SARS coronavirus, and influenza A virus

Host cell proteases: critical determinants of coronavirus tropism and pathogenesis

SARS-CoV-2 cell entry depends on ACE2 and TMPRSS2 and is blocked by a clinically proven protease inhibitor

Dipeptidyl peptidase 4 is a functional receptor for the emerging human coronavirus-EMC

CD147-spike protein is a novel route for SARS-CoV-2 infection to host cells

Look AT & Holmes KV

Human aminopeptidase N is a receptor for human coronavirus 229E

TMPRSS2 and TMPRSS4 promote SARS-CoV-2 infection of human small intestinal enterocytes

Structure, function, and antigenicity of the SARS-CoV-2 spike glycoprotein

LY6E impairs coronavirus fusion and confers immune control of viral disease

Swaminathan S; COVID-19 and ACE2 in Cardiovascular, Lung, and Kidney Working Group (2020) Acute kidney injury in COVID-19: emerging evidence of a distinct pathophysiology

SARS-CoV-2 entry factors are highly expressed in nasal epithelial cells together with innate immune genes

Comparative genetic analysis of the novel coronavirus (2019-nCoV/ SARS-CoV-2) receptor ACE2 in different populations

Genetic roadmap for kidney involvement of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Infection

The single-cell transcriptomic landscape of early human diabetic nephropathy

Homozygous L-SIGN (CLEC4M) plays a protective role in SARS coronavirus infection

A single-nucleus RNAsequencing pipeline to decipher the molecular anatomy and pathophysiology of human kidneys

SARS-CoV-2 receptor networks in diabetic and COVID-19-associated kidney disease

Glomerular localization and expression of Angiotensin-converting enzyme 2 and Angiotensinconverting enzyme: implications for albuminuria in diabetes

Single cell transcriptomics identifies focal segmental glomerulosclerosis remission endothelial biomarker

Tubular cell and keratinocyte single-cell transcriptomics applied to lupus nephritis reveal type I IFN and fibrosis relevant pathways

Renal compartment-specific genetic variation analyses identify new pathways in chronic kidney disease

PROMO: detection of known transcription regulatory elements using speciestailored searches

Critical role for STAT4 activation by type 1 interferons in the interferon-gamma response to viral infection

Biochemistry and biology of the inducible multifunctional transcription factor TFII-I: 10 years later

Targeting the endocytic pathway and autophagy process as a novel therapeutic strategy in COVID-19

Proteolytic cleavage of the SARS-CoV-2 spike protein and the role of the novel S1/S2 site

Clinical characteristics and risk factors for severe COVID-19 in hospitalized kidney transplant recipients: a multicentric cohort study

COVID-19 and systemic lupus erythematosus: a case series

Histopathological findings and viral tropism in UK patients with severe fatal COVID-19: a post-mortem study

Evidence for and against direct kidney infection by SARS-CoV-2 in patients with COVID-19

SARS-CoV-2 can be detected in urine, blood, anal swabs, and oropharyngeal swabs specimens

Evidence for SARS-CoV-2 spike protein in the urine of COVID-19 patients

Kidney biopsy findings in patients with COVID-19

A novel soluble ACE2 variant with prolonged duration of action neutralizes SARS-CoV-2 infection in human kidney organoids

ACE2, the kidney and the emergence of COVID-19 two decades after ACE2 discovery

Proof-of-concept workflow for establishing reference intervals of human urine proteome for monitoring physiological and pathological changes

Integrating single-cell transcriptomic data across different conditions, technologies, and species

clusterProfiler: an R package for comparing biological themes among gene clusters

The authors declare no conflict of interest.

HP conceived the project and designed the study. HL, FX, and XM conducted all single-cell analysis and prepared the figures. HS carried out the histological examination and immunohistochemical staining. LS carried out experiments and analyzed the urinary proteomics data. CZ and ZZ supervised these works. XM, HL, FX, YS, and YL participated in results interpretation. XM, HL, and HP wrote the manuscript, and all the authors approved for the manuscript submission.

The peer review history for this article is available at https://publons.com/publon/10.1111/febs.16114.