Skip to main content
  • Research article
  • Open access
  • Published:

Possible kidney-lung cross-talk in COVID-19: in silico modeling of SARS-CoV-2 infection

Abstract

Background

Publicly available genomics datasets have grown drastically during the past decades. Although most of these datasets were initially generated to answer a pre-defined scientific question, their repurposing can be useful when new challenges such as COVID-19 arise. While the establishment and use of experimental models of COVID-19 are in progress, the potential hypotheses for mechanisms of onset and progression of COVID-19 can be generated by using in silico analysis of known molecular changes during COVID-19 and targets for SARS-CoV-2 invasion.

Methods

Selecting condition: COVID-19 infection leads to pneumonia and mechanical ventilation (PMV) and associated with acute kidney injury (AKI). There is increasing data demonstrating mechanistic links between AKI and lung injury caused by mechanical ventilation.

Selecting targets: SARS-CoV-2 uses angiotensin-converting enzyme 2 (ACE2) and transmembrane protease serine 2 (TMPRSS2) for cell entry. We hypothesized that expression of ACE2 and TMPRSS2 would be affected in models of AKI and PMV. We therefore evaluated expression of ACE2 and TMPRSS2 as well as other novel molecular players of AKI and AKI-lung cross-talk in the publicly available microarray datasets GSE6730 and GSE60088, which represent gene expression of lungs and kidneys in mouse models of AKI and PMV, respectively.

Results

Expression of COVID-19 related genes ACE2 and TMPRSS2 was downregulated in lungs after 6 h of distant AKI effects. The expression of ACE2 decreased further after 36 h, while expression of TMPRSS2 recovered. In kidneys, both genes were downregulated by AKI, but not by distant lung injury. We also identified 53 kidney genes upregulated by PMV; and 254 lung genes upregulated by AKI, 9 genes of which were common to both organs. 3 of 9 genes were previously linked to kidney-lung cross-talk: Lcn2 (Fold Change (FC)Lung (L) = 18.6, FCKidney (K) = 6.32), Socs3 (FCL = 10.5, FCK = 10.4), Inhbb (FCL = 6.20, FCK = 6.17). This finding validates the current approach and reveals 6 new candidates, including Maff (FCL = 7.21, FCK = 5.98).

Conclusions

Using our in silico approach, we identified changes in COVID-19 related genes ACE2 and TMPRSS2 in traditional mouse models of AKI and kidney-lung cross-talk. We also found changes in new candidate genes, which could be involved in the combined kidney-lung injury during COVID-19.

Peer Review reports

Background

Publicly available genomics datasets have grown drastically during the past decades. Although most datasets were initially generated to answer a pre-defined scientific question, the enormous amount of accumulated data makes genomics data repurposing possible when new challenges such as COVID-19 arise and need to be rapidly dealt with. These data sets are also very useful to study when traditional wet-lab research is slowed down during the pandemic, and experimental models of COVID-19 are still under development.

Clinical data for COVID-19 suggested that its mortality and morbidity can be linked not only to respiratory complications but also to other organ dysfunction [1]. Acute kidney injury (AKI) is associated with COVID-19. This association is thought to be secondary to hemodynamic effects of the disease related septic condition, and the poor outcome is secondary to sepsis, where AKI is often a surrogate marker of severity of illness in these patients [2, 3]. There is also a strong association between AKI and dysfunction of extra-renal organs, and more recent animal research has shown a significant causal effect of AKI on a distant organ dysfunction [48]. Despite frustrating outcomes, little is known about the potential pathophysiological interactions between the kidney and the extrarenal organs in critically ill patients. The critical illness during COVID-19 warrants study COVID-19 related genes and kidney-lung cross-talk to discover potential therapeutic targets.

We previously used combined genomics data analysis of kidney and lung to discover new molecular pathways of kidney-lung cross-talk during lung injury [9], and demonstrated that inflammation is a major component of the initiation and exacerbation of AKI. To relate our findings to COVID-19, we re-visited the global gene expression profiles of kidney in mice with pneumonia and mechanical ventilation [10, 11] with the main focus on changes in two important COVID-19 related genes angiotensin-converting enzyme 2 (ACE2) and transmembrane protease serine 2 (TMPRSS2). These genes were not previously considered during prior studies of kidney-lung cross-talk. However, SARS-CoV-2 uses ACE2 to enter both kidney and lung, and is highly expressed in proximal tubule cells of the kidney [10]. While virus enters the lung directly by inhalation the kidney becomes infected via blood (viremia), that occurs in ~ 15% of patients and is reported in association with severe infection [12]. The rate of coronoviral infection is multiplied by thousand times in the presence of co-expressed proteases [13]. TMPRSS2 gene codes for such a protease, which cleaves (primes) S protein of SARS-CoV-2 bound to ACE2 receptor [12]. TMPRSS2 is abundantly expressed in the distal nephron [14]. The data from Human Protein Atlas (https://www.proteinatlas.org) also demonstrates high co-expression of ACE2 and TMPRSS2 in microvilli of proximal tubules.

In addition to detected significant changes in expression of ACE2 and TMPRSS2 in AKI and PMV models, we also identified novel candidates that mediate cross-talk between the lung and kidney, which could be pertinent to COVID-19.

Methods

Study aims

This expression data repurposing study aims to extrapolate the molecular changes during COVID-19 to cross-talk between the injured kidney and lung using publicly available genomic data. Mouse models were selected based upon the phenotypes similar to those in patients with COVID-19. Two major COVID-19 related clinical conditions were selected: PMV and AKI. We hypothesized that genes involved in SARS-CoV-2 invasion will be affected in these models and might be potentially involved in COVID-19 outcome.

The microarray Genomic Series GSE60088 and GSE6730 were downloaded from the Gene Expression Omnibus (GEO).

Brief description of gene expression series 60,088: The S. aureus (SA) ~ 107 cfu was deposited in mice oropharynx (n = 5) and mice were intubated via tracheostomy with a 20-gauge blunt metal catheter. Although the injurious pathogen in the model was bacterium the recent clinical observations showed that 37% of critically ill COVID-19 patients acquired bacterial pneumonia [15], which supports the model selection. Intubated mice were connected to a MiniVent rodent ventilator and mechanically ventilated (MV) with a tidal volume of 10 mL/kg, a respiratory rate of 150 breaths per minute. Control mice (n = 5) were maintained in their cages. After 6 h the lungs and kidneys were collected for RNA isolation. The Office of Animal Welfare at the University of Washington approved all experiments.

Brief description of gene expression series 6730: Male 6–8-week-old mice (C57BL6/J), weighing approximately 25-30 g were placed on a heating blanket and underwent midline laparotomy with isolation of bilateral renal pedicles. For mice assigned to experimental ischemia-reperfusion injury (IRI), a non-traumatic microvascular clamp was applied across both renal pedicles for 60 min. After the allotted ischemia time, the clamps were gently removed, the animals were administered 1 ml of sterile saline intraperitoneally, and the incision was closed in two layers with 4-0 silk suture. The animals were then allowed to recover with free access to food and water. Sham animals underwent the identical procedure without placement of the vascular clamps. At 6 or 36 h following the experimental procedure, the mice were euthanized by exsanguination under pentobarbital anesthesia and lung tissues were collected for the analysis. All procedures were approved by the Johns Hopkins Animal Care and Use Committee, and were consistent with the National Institutes of Health (NIH) Guide for the Care and Use of Laboratory Animals.

Mouse Genome 430 2.0 microarray (Affymetrix) was used for both models.

Statistical analysis

Differentially expressed genes were identified using GEO-build-in GEO2R tool. For PMV model the groups were defined as follow: Control Lung Genomic Sample GSM1464822-GSM1464825, MV + SA Lung GSM1464826-GSM14648230. Control Kidney GSM1464839-GSM1464843, MV + SA Kidney GSM1464844-GSM1464848.

For AKI model the groups were defined as follow: Control 36 h GSM155092- GSM155096, AKI 36 h GSM155100- GSM155102. Control 6 h GSM155086- GSM155088, AKI 6 h GSM155089- GSM155091. Then GEO2R outputs comparing all groups were downloaded.

The main SARS-CoV-2 targets: angiotensin-converting enzyme 2 (ACE2) and transmembrane protease serine 2 (TMPRSS2) were selected for the detailed expression analysis.

Results

We identified changes in lungs and kidneys of COVID-19 related genes ACE2 and TMPRSS2 in pneumonia and mechanical ventilation (PMV) and ischemic AKI-induced lung dysfunction (Table 1). Expression of ACE2 in lungs was significantly downregulated after 6 h in PMV model, and was similar to the gene downregulation in initially normal lungs at 6 h after the induction of ischemic injury to kidneys in AKI model. Lung ACE2 was further downregulated 36 h after ischemic AKI. (Table 1). Expression of TMPRSS2 in the lung was also significantly decreased after PMV, as it was 6 h after ischemic AKI. However, at 36 h of AKI lung TMPRSS2 was similar to the baseline.

Table 1 Response of ACE2 and TMPRSS2 to direct injurious stimuli in kidney and lung

In kidneys ACE2 sharply decreased at 6 h after AKI and was further downregulated 36 h after the induction of ischemic AKI. However, kidney ACE2 was unaffected by PMV. Kidney TMPRSS2 was unaffected in PMV model, but decreased 6 h after ischemic AKI. (Table 1). We then explored other molecular targets in both organs by filtering for the known genes with P-value < 0.01 and fold change (FC) > 4. This approach identified 53 genes upregulated in the kidneys from mice with PMV; and 254 genes upregulated in the lungs from mice with AKI (Fig. 1). Cross-referencing of these candidates identified 9 common genes (Table 2). Three of them Lcn2, Socs3, and Inhbb had already been associated with kidney-lung cross-talk [9], which validate by our current approach. We also identify 6 novel molecular targets for kidney lung-cross-talk.

Fig. 1
figure 1

Theoretical model of kidney-lung interaction during COVID-19. We postulate that SARS-CoV2 enters a lung and perhaps a kidney using ACE2 as a main entry point and TMPRSS2 for proteolitic cleavage (priming) of SARS-CoV2-ACE2 complex to increase the infection rate. Ensuing injury (“second hit”), such as ventilation injury and AKI (local tissue inflammation) can alter the expression of these two genes and change response to the infection. The left panel, therefore, represents ventilation injury and right panel represents AKI. Expression values for ACE2 and TMPRSS2 genes are depicted in closed boxes and are taken from Table 1. Solid boxes/lines represent direct effects of an injury, and dashed boxes/lines represent distant genomic effects at 6 h after the injury. Aside from ACE2/TMPRSS2 axis additional pathways can be altered in a lung and kidney that might modify survival and tissue injury/repair. Genes from these key pathways could be involved in lung-kidney cross-talk in COVID-19 patient. We have identified 9 genes, expression of which was significantly changed during experimental PMV and AKI; theses changes were common to both organs during injurious response. Out of common 9 genes 3 were previously described as AKI-related genes (Lcn2, Socs3, and Inhbb); other six genes were novel (Mt2, Maff, Junb, Hmgcs2, Tnfrsf12a and Ch25h)

Table 2 Gene candidates distantly upregulated in kidney or lung by injury to distant organ

Discussion

Severe COVID-19 is usually characterized by lung injury and often associated with acute kidney injury. There has been mechanistic studies of combined lung and kidney injury in animal models, but not related to COVID-19. Prior studies focusing on mechanisms of kidney-lung cross-talk did not examine COVID-19 specific genes because they were not relevant at the time. Publicly available databases can be re-analyzed while labs are slowed down and animal models of COVID-19 are being developed. We used this approach and found unexpected decreases in key SARS-CoV2 molecules ACE2 and TMPRSS2 in two models of kidney lung cross-talk: 1) pneumonia and mechanical ventilation as a trigger of kidney injury and 2) AKI as a trigger of lung injury. Consequently, we found significant changes in expression of 3 genes that were previously associated with kidney-lung cross-talk, and also identified 6 novel targets.

ACE2 and TMPRSS2 are important for SARS-CoV2 entry into lungs via inhalation and kidneys via blood stream. It is believed that expression level of these molecules predisposes to viral invasion, especially the co-expression of both genes in targeted cells, which significantly accelerates infection and disease development. This could explain why some patients get more serious injury while others do not [16]. There are also “second hits” that can alter patient outcome, such as bacterial pneumonia, mechanical ventilation induced injury, and AKI. We examined expression of these two molecules in two traditional published models and expected an increase in mRNA expression of the receptors. Unexpectedly, both ACE2 and TMPRSS2 decreased in lung and kidney during PMV and AKI, though with different intensities and kinetics. This decrease is a real suppression of expression, rather than a reflection of a cell death after the injury, given that 53 kidney genes and 254 lung genes demonstrated significant increase in response to the injury. The meaning of this is unknown and could be tested in the future studies when experimental models of COVID-19 are better developed. The matter is also complicated by the known difference in expression of TPMRSS2 in the mouse and human [17]. However, one still can speculate that “second hits” could make it more difficult for SARS-CoV-2 to enter lung and kidney cells and might be protective. Based on the level of gene expression in Table 2, we directly extrapolated the gene expression values to receptor expression on the surface of a cell and built our protective model against SARS-CoV-2 invasion (Fig. 2). This model might suggest a straightforward reduction in the severity of COVID-19 infection, which was not noted during this pandemic. However, the future comparison of virus burden between patients with and without AKI can address feasibility of this model.

Fig. 2
figure 2

Protective model of a “second hit” during distant organ injury. The schematic cells represent proximal tubule cells of kidney. It has been shown in the Human Protein Atlas that ACE2 and TMPRSS2 are co-expressed in proximal tubules, especially in its microvilli (https://www.proteinatlas.org/ENSG00000130234-ACE2/tissue/kidney and https://www.proteinatlas.org/ENSG00000184012-TMPRSS2/tissue/kidney). The virus invasion flow was adopted from Parmar et al. [12]. The proposed model demonstrates that with increased time after the initial injury it becomes more difficult for SARS-CoV-2 to enter the injured lung or kidney cell. 6 h after the injury the both targets of SARS-CoV-2 (ACE2 and TMPRSS2) are downregulated. 36 h after the injury, while the expression of TMPRSS2 is recovered, the expression of ACE2 is downregulated even more than after 6 h

While focusing on the important COVID-19 targets ACE2 and TMPRSS2 molecules, we re-analyzed these data sets and confirmed that our approach was valid when we identified 3 previously known candidates of lung-kidney cross-talk during AKI: lipocalin 2 (NGAL), suppressor of cytokine signaling 3 (member of STAT family), and inhibin beta B. NGAL is a bio-marker of AKI but also has important role in iron metabolism, epithelial cell and immune cell functions. Socs3 in stressed proximal tubules plays an important role during AKI by inhibition of reparative proliferation [18] and its anti-inflammatory properties might also play protective role in kidney injury. Inhibin beta B is a member of the transforming growth factor β (TGF-β) superfamily. It has been reported that this gene is significantly up-regulated in renal endothelial cells from male mice with ischemia reperfusion injury [19]. We also identified 6 new candidate kidney-lung cross-talk genes (Table 2). We found Maff gene particularly interesting. This gene is abundantly expressed in the kidney and its up-regulation has been linked to deleterious effects on podocytes [20]. Maff also binds the oxytocin receptor promoter [21], which could contribute to sex differences in disease outcome – seen in COVID-19 [22]. Moreover, it had been demonstrated in the highly related mouse model of SARS-CoV that ovariectomised or antiestrogen-treated female mice had more severe CoV infections than their respective controls [23].

Conclusions

Study of expression patterns of ACE2 and TMPRSS2 genes during PMV and AKI could be important in understanding the pathogenesis of COVID-19 in critically ill patients. These two genes may also interplay with other key genes during kidney-lung cross-talk. Analytical methods of repurposing available genomics data allowed us to generate new COVID-19 related hypotheses, which can be studied and validated once reliable COVID-19 models are established. Moreover, this in silico repurposing method is generalizable and can be applied to different studies.

Availability of data and materials

The data supporting the results reported in this manuscript can be found at GEO:

https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE60088

https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE6730

Abbreviations

ACE2:

Angiotensin-converting enzyme 2

AKI:

Acute kidney injury

CoV:

Corona virus

COVID:

Corona virus disease

FC:

Fold change

GEO:

Gene Expression Omnibus

IRI:

Ischemia reperfusion injury

Maff:

Musculoaponeurotic fibrosarcoma oncogene family F

NGAL:

Neutrophil gelatinase-associated lipocalin

PMV:

Pneumonia and mechanical ventilation

SARS:

Severe acute respiratory syndrome

Socs3:

Suppressor of cytokine signaling 3

TGF-β:

Transforming growth factor β

TMPRSS2:

Transmembrane protease serine 2

References

  1. Kow CS, Hasan SS. Possible protective effect of renin-angiotensin system inhibitors in COVID-19 induced acute kidney injury. JASN. 2020;31:1917–8.

    Article  CAS  Google Scholar 

  2. Kellum JA, Nadim MK, Forni LG. Sepsis-associated acute kidney injury: is COVID-19 different? Kidney Int. 2020;98:1370–2.

    Article  CAS  Google Scholar 

  3. Parmar MS. Acute kidney injury associated with COVID-19-cumulative evidence and rationale supporting against direct kidney injury (infection). Nephrology (Carlton). 2021;26:239–47.

    Article  CAS  Google Scholar 

  4. Hoke TS, Douglas IS, Klein CL, He Z, Fang W, Thurman JM, et al. Acute renal failure after bilateral nephrectomy is associated with cytokine-mediated pulmonary injury. JASN. 2007;18:155–64.

    Article  CAS  Google Scholar 

  5. Hassoun HT, Grigoryev DN, Lie ML, Liu M, Cheadle C, Tuder RM, et al. Ischemic acute kidney injury induces a distant organ functional and genomic response distinguishable from bilateral nephrectomy. Am J Physiol Renal Physiol. 2007;293:F30–40.

    Article  CAS  Google Scholar 

  6. Kelly KJ. Distant effects of experimental renal ischemia/reperfusion injury. JASN. 2003;14:1549–58.

    Article  CAS  Google Scholar 

  7. Rabb H, Wang Z, Postler G, Soleimani M. Possible molecular basis for changes in potassium handling in acute renal failure. Am J Kidney Dis. 2000;35:871–7.

    Article  CAS  Google Scholar 

  8. Zarbock A, Schmolke M, Spieker T, Jurk K, Van Aken H, Singbartl K. Acute uremia but not renal inflammation attenuates aseptic acute lung injury: a critical role for uremic neutrophils. JASN. 2006;17:3124–31.

    Article  CAS  Google Scholar 

  9. Grigoryev DN, Liu M, Hassoun HT, Cheadle C, Barnes KC, Rabb H. The local and systemic inflammatory Transcriptome after acute kidney injury. JASN. 2008;19:547–58.

    Article  CAS  Google Scholar 

  10. Gu C, Qiao W, Wang L, Li M, Song K. Identification of genes and pathways associated with multiple organ dysfunction syndrome by microarray analysis. Mol Med Report. 2018. https://doi.org/10.3892/mmr.2018.8973.

  11. Gharib SA, Mar D, Bomsztyk K, Denisenko O, Dhanireddy S, Liles WC, et al. System-wide mapping of activated circuitry in experimental systemic inflammatory response syndrome. Shock. 2016;45:148–56.

    Article  CAS  Google Scholar 

  12. Parmar MS. TMPRSS2: an equally important protease as ACE2 in the pathogenicity of SARS-CoV-2 infection. Mayo Clin Proc. 2021;96:2748–52.

    Article  CAS  Google Scholar 

  13. Matsuyama S, Ujike M, Morikawa S, Tashiro M, Taguchi F. Protease-mediated enhancement of severe acute respiratory syndrome coronavirus infection. Proc Natl Acad Sci U S A. 2005;102:12543–7.

    Article  CAS  Google Scholar 

  14. Ni W, Yang X, Yang D, Bao J, Li R, Xiao Y, et al. Role of angiotensin-converting enzyme 2 (ACE2) in COVID-19. Crit Care. 2020;24:422.

    Article  Google Scholar 

  15. Dudoignon E, Caméléna F, Deniau B, Habay A, Coutrot M, Ressaire Q, et al. Bacterial pneumonia in COVID-19 critically ill patients: a case series. Clin Infect Dis. 2021;72:905–6.

    Article  CAS  Google Scholar 

  16. Batlle D, Soler MJ, Sparks MA, Hiremath S, South AM, Welling PA, et al. Acute kidney injury in COVID-19: emerging evidence of a distinct pathophysiology. JASN. 2020;31:1380–3.

    Article  CAS  Google Scholar 

  17. Vaarala MH, Porvari KS, Kellokumpu S, Kyllönen AP, Vihko PT. Expression of transmembrane serine protease TMPRSS2 in mouse and human tissues. J Pathol. 2001;193:134–40.

    Article  CAS  Google Scholar 

  18. Susnik N, Sörensen-Zender I, Rong S, von Vietinghoff S, Lu X, Rubera I, et al. Ablation of proximal tubular suppressor of cytokine signaling 3 enhances tubular cell cycling and modifies macrophage phenotype during acute kidney injury. Kidney Int. 2014;85:1357–68.

    Article  CAS  Google Scholar 

  19. Viñas JL, Porter CJ, Douvris A, Spence M, Gutsol A, Zimpelmann JA, et al. Sex diversity in proximal tubule and endothelial gene expression in mice with ischemic acute kidney injury. Clin Sci. 2020;134:1887–909.

    Article  Google Scholar 

  20. Okabe M, Motojima M, Miyazaki Y, Pastan I, Yokoo T, Matsusaka T. Global polysome analysis of normal and injured podocytes. Am J Physiol Ren Physio. 2019;316:F241–52.

    Article  Google Scholar 

  21. Kimura T, Ivell R, Rust W, Mizumoto Y, Ogita K, Kusui C, et al. Molecular cloning of a human MafF homologue, which specifically binds to the oxytocin receptor gene in term myometrium. Biochem Biophys Res Commun. 1999;264:86–92.

    Article  CAS  Google Scholar 

  22. Jin J-M, Bai P, He W, Wu F, Liu X-F, Han D-M, et al. Gender differences in patients with COVID-19: focus on severity and mortality. Front Public Health. 2020;8:152.

    Article  Google Scholar 

  23. Channappanavar R, Fett C, Mack M, Ten Eyck PP, Meyerholz DK, Perlman S. Sex-based differences in susceptibility to severe acute respiratory syndrome coronavirus infection. JI. 2017;198:4046–53.

    CAS  Google Scholar 

Download references

Acknowledgments

The study investigators thank the owners of the raw expression data, who made it available for other investigators.

Funding

This in silico study was conducted outside of any funding agencies.

Author information

Authors and Affiliations

Authors

Contributions

HR proposed and guided the study, revised the manuscript. DNG participated in protocol editing and conduct the study. All authors have approved the final manuscript and certify that this manuscript represents valid work and has not been previously published.

Corresponding author

Correspondence to Dmitry N. Grigoryev.

Ethics declarations

Ethics approval and consent to participate

Given that the data for this study was obtained from publicly available repository, the authors do not have a detailed ethics approval information. Therefore, ethical statements are reproduced from the data submission page of GEO and pertinent manuscripts.

GSE6730: All procedures were approved by the Johns Hopkins Animal Care and Use Committee, and were consistent with the National Institutes of Health (NIH) Guide for the Care and Use of Laboratory Animals.

GSE60088: The Office of Animal Welfare at the University of Washington approved all experiments.

Consent for publication

Both authors are consenting for the publication.

Competing interests

The authors have no competing interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Grigoryev, D.N., Rabb, H. Possible kidney-lung cross-talk in COVID-19: in silico modeling of SARS-CoV-2 infection. BMC Nephrol 23, 57 (2022). https://doi.org/10.1186/s12882-022-02682-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12882-022-02682-1

Keywords