key: cord-0985327-7ovsux3p
authors: Baratchian, Mehdi; McManus, Jeffrey M.; Berk, Mike P.; Nakamura, Fumihiko; Mukhopadhyay, Sanjay; Xu, Weiling; Erzurum, Serpil; Drazba, Judy; Peterson, John; Klein, Eric A.; Gaston, Benjamin; Sharifi, Nima
title: Androgen regulation of pulmonary AR, TMPRSS2 and ACE2 with implications for sex-discordant COVID-19 outcomes
date: 2021-05-27
journal: Sci Rep
DOI: 10.1038/s41598-021-90491-1
sha: 9fec61c3c6a62e55d12d644b4bdb8c2a9f51a68e
doc_id: 985327
cord_uid: 7ovsux3p

The sex discordance in COVID-19 outcomes has been widely recognized, with males generally faring worse than females and a potential link to sex steroids. A plausible mechanism is androgen-induced expression of TMPRSS2 and/or ACE2 in pulmonary tissues that may increase susceptibility or severity in males. This hypothesis is the subject of several clinical trials of anti-androgen therapies around the world. Here, we investigated the sex-associated TMPRSS2 and ACE2 expression in human and mouse lungs and interrogated the possibility of pharmacologic modification of their expression with anti-androgens. We found no evidence for increased TMPRSS2 expression in the lungs of males compared to females in humans or mice. Furthermore, in male mice, treatment with the androgen receptor antagonist enzalutamide did not decrease pulmonary TMPRSS2. On the other hand, ACE2 and AR expression was sexually dimorphic and higher in males than females. ACE2 was moderately suppressible with enzalutamide administration. Our work suggests that sex differences in COVID-19 outcomes attributable to viral entry are independent of TMPRSS2. Modest changes in ACE2 could account for some of the sex discordance.

www.nature.com/scientificreports/ hormonal therapy for prostate cancer have yielded differing results on a possible protective effect of androgen suppression on risk of COVID-19 7, 8 . Androgen regulation of TMPRSS2 raises the possibility that the physiological roles of androgens may, at least partially, account for the sex-specific clinical outcomes 9, 10 . Utilizing a high-throughput drug screening strategy, a recent study found that ACE2 levels in human alveolar epithelial cells can be downregulated by 5α-reductase inhibitors, suggesting an androgen-driven mode of expression 11 . Furthermore, due to its androgenregulated nature in the prostate and its essential role in SARS-CoV-2 etiology, TMPRSS2 expression has been postulated to follow a similar pattern of regulation in pulmonary cells by the potent androgens testosterone and dihydrotestosterone 12 . If this link proves correct, it could pave the path to novel strategies, including repurposing of FDA-approved potent androgen synthesis inhibitors or AR antagonists, such as enzalutamide (Enz) and apalutamide, for the treatment of COVID-19. These strategies are the subject of several clinical trials (e.g., NCT04374279, NCT04475601, NCT04509999, NCT04397718) 5, 13 .

Here, we show that the expression of pulmonary AR and ACE2 follows a sex-discordant pattern with males expressing considerably higher levels of protein than females. In humans, there is no difference in ACE2 expression between non-smoking men and women, while in contrast, ACE2 expression is significantly higher in the lungs of male smokers. We provide in vivo evidence in mice that neither mRNA nor protein levels of TMPRSS2 vary by sex or treatment with the potent AR-antagonist Enz. ACE2 expression however is modestly modifiable by anti-AR treatment and may to some extent explain the sex disparities in susceptibility to SARS-CoV-2.

Sexually dimorphic AR expression and ACE2 dimorphism in smokers. Certain pulmonary disease outcomes, including asthma, are sex steroid-associated 14 . Considering the poorer clinical outcome of COVID-19 in men, underlying androgen-related causes are suspected but not presently known. The SARS-CoV-2 co-receptor TMPRSS2 harbors an AR-responsive enhancer that is induced by androgens in prostate tissue 15 , raising the possibility of a similar mode of regulation in the respiratory system. We first asked whether, similar to TMPRSS2, ACE2 was also regulated by AR signaling in the human prostate adenocarcinoma cell line originally derived from a lymph node metastasis (LNCaP), which is AR expressing and androgen responsive. Indeed, both mRNA and protein expression of ACE2 were strongly induced by the synthetic androgen R1881 and suppressed by Enz-mediated AR blockade (Fig. 1A ,B, S1 and S2). The specificity of the changes in ACE2 protein levels was confirmed with knockdown experiments, using two different antibodies (Fig. S1 ). Moreover, ChIP-seq analysis of AR cistrome revealed multiple AR-binding sites upstream of the ACE2 region that were lost upon Enz treatment (Fig. 1C) . These findings collectively indicate that ACE2 is indeed an androgen-driven gene in prostate cells.

We next sought to investigate whether male sex was associated with higher expression of ACE2, TMPRSS2 or AR in human lung. To this end, we acquired the publicly available expression datasets in non-cancerous lung and associated respiratory tissues from the Genomic Expression Omnibus (GEO). Across all tissue-type comparisons, we found no evidence for elevated ACE2, TMPRSS2 or AR mRNA expression in males compared with females ( Fig. 2 ). Next, we performed immunohistochemical (IHC) analyses to explore the possibility of a sex-specific pattern of protein expression within distinct pulmonary cell types. The IHC-stained mouse and human lung specimens were assessed by an expert pulmonary pathologist (SM) and semi-quantified for protein expression using the H-score, a common and well-known semiquantitative scoring system that takes into consideration staining intensity as well as the percentage of cells staining positively [16] [17] [18] . H-scores can range from 0 (no expression) to 300 (3 + staining intensity in 100% of cells).

Pulmonary TMPRSS2, predominantly expressed in airways (bronchial or bronchiolar epithelial cells), was neither associated with sex in humans ( Fig. 3A ) nor in mice (Fig. 3D) . Furthermore, smoking did not appear to affect TMPRSS2 expression (Fig. 3A) . In human lungs, alveolar epithelial cells also stained positive for TMPRSS2; expression was unchanged by smoking status or sex.

In contrast, we discovered that AR and ACE2 demonstrate sex-discordant expression; AR showed limited expression in airway cells of human males which was absent in human females (Fig. 3B ). Smoking further elevated AR levels to a modest degree in men but the difference did not reach statistical significance (Fig. 3B,H) . Consistent with the overall pattern observed in humans, lungs of male mice expressed markedly higher levels of AR protein than females (Figs. 3E and 4D). Alveolar pneumocytes (epithelial cells) were the dominant ARexpressing population in mouse lungs (Fig. 3G ). Since the AR transcript levels did not differ by sex (Figs. 2B and 4C), we speculated that the high concentrations of circulating androgens in males may stabilize pulmonary AR expression 19 .

In the case of ACE2, we detected no significant difference between non-smoker men and women. Nevertheless, its expression was elevated in the lungs of male smokers compared with female smokers (Fig. 3C,I) . Thus, we observed sexually dimorphic ACE2 expression in smokers only. In mice, the sex disparity was readily detectable, with males expressing significantly higher ACE2 protein in airways than females (Fig. 3F ). These data, together, suggest that pulmonary ACE2 and AR expression is sexually dimorphic. www.nature.com/scientificreports/ TMPRSS2 expression is unaffected but ACE2 is modestly suppressed by potent AR blockade. To further explore the possibility of inhibiting TMPRSS2 expression by means of AR blockade, we harvested lungs from female and male mice treated with control diet or Enz for > 10 days and analyzed them for protein and mRNA expression. The animal studies performed on bulk lungs yielded results consistent with the human expression analysis: i.e., there were no sex-specific TMPRSS2 changes. Additionally, Enz treatment failed to detectably downregulate TMPRSS2 (Fig. 4A,D) . It has been previously noted that, due to glycosylation, fulllength TMPRSS2 may be detected by SDS-PAGE at a higher molecular weight of approximately 70 kDa 20 . In immunoblots performed with maximum-sensitivity substrate and extended exposures, we did not detect distinguishable bands for TMPRSS2 above 50 kDa (Fig. 4D ). Given the unchanged AR protein levels (despite a modest but statistically significant mRNA increase) in Enz-treated males, we infer that TMPRSS2 is not regulated by AR in the lung. Nevertheless, high resolution techniques such as single cell-sequencing or -proteomics may be employed to explore otherwise undetectable alterations of expression in minority cell populations. www.nature.com/scientificreports/ Similar to the human data, we found no evidence of sex-specific changes in AR transcript levels in mice (Fig. 4C ). The lungs of male mice, however, expressed significantly higher amounts of AR protein compared with females ( Fig. 4D) . As mentioned previously, this may be due to AR stabilizing effects of abundant circulating androgens in males. Similar to AR, ACE2 expression also showed clear sexual dimorphism, with males expressing substantially higher levels of protein (Fig. 4D ) and modestly but significantly higher amounts of mRNA (Fig. 4B ) in bulk lungs compared with females. Treatment with Enz lowered the transcript levels in males down to female levels ( Fig. 4B ) and also partially reduced protein quantities (Fig. 4D) as evidenced by immunoblotting. Some subtle ACE2 protein levels were apparent in the control groups that might be due to biologic variations. Nevertheless, the modest Enz-mediated suppression of ACE2 was not clearly captured in our IHC analyses. This could be explained by possible marginal changes of ACE2 within pulmonary cells that may fall below the IHC detection limit, however, accumulatively can be observed through total protein detection methods performed on bulk tissue. Overall, these data indicate sex-discordant AR and ACE2 regulation, and a potential androgen-regulated mode of expression for pulmonary ACE2. Therefore, androgen-AR-mediated mechanisms could explain sexspecific differences in COVID-19 outcomes by TMPRSS2-independent mechanisms. www.nature.com/scientificreports/ ACE2 and TMPRSS2 mRNA expression increases in current smokers whereas in former smokers, expression returns to levels found in never-smokers. In addition to male sex, smoking is a risk factor for COVID-19 susceptibility and poor clinical outcomes 21 . One recent study of 1,099 COVID-19-positive patients reported a more than two-fold increased risk for intensive care unit admission and death in smokers as compared with non-smokers 22 . We identified human expression GEO datasets of bronchial/airway epithelial cells containing subject smoking status and asked whether smoking is associated with TMPRSS2 expression. Our analysis indicated a consistent pattern whereby expression of both TMPRSS2 (Fig. 5A ) and the primary SARS-CoV-2 receptor ACE2 (Fig. 5B ) was modestly but significantly increased in smokers compared with nonsmokers. Interestingly, the levels were downregulated to never-smoker levels in former smokers. The results of our analysis are in keeping with several recent reports on ACE2 and smoking [23] [24] [25] . Although the p-values range widely across different data sets, almost all show changes in a consistently increased direction for current smokers including multiple data sets with small p-values (Fig. 5A,B) . Finally, for both TMPRSS2 (Fig. 5C) and ACE2 (Fig. 5D) , there was no correlation between smoking pack-years and mRNA expression in either current or former smokers, suggesting that the change does not build up over time but is instead a rapid process akin to a switch.

Sex-associated clinical outcomes have been long observed in a variety of infectious and inflammatory conditions. Sex steroids (i.e., androgens and estrogens) are possible mediators of these biologic differences. For COVID-19, potential androgen-mediated biologic differences include (1) SARS-CoV2 cellular receptor regulation and (2) immune modulation 5, 13 . Our study addresses the first of these possibilities. In both humans and mice, AR protein is clearly more highly expressed in the lungs of males compared with females. Detection of AR protein in male compared with female lung might be surprising in the absence of a detectable difference in AR transcript. This might be explained by stabilization of AR protein in the presence of androgens, which is recognized to occur in prostatic tissues 19 . To our knowledge, there are no prior reports of this sexually dimorphic pulmonary AR expression. The specific presence of AR protein in male lungs raises the question of its transcriptional program-namely expression of TMPRSS2 and ACE2.

We find no evidence for androgen regulation of TMPRSS2 in pulmonary tissues. This evidence includes (1) the absence of any TMPRSS2 increase in male compared with female human lung (2) no TMPRSS2 increase in male compared with female mouse lung (3) no evidence for TMPRSS2 suppression with next-generation AR antagonist treatment. Our observations are in agreement with a new study that provides evidence for the absence of AR-binding and open chromatin state within the TMPRSS2 locus in lung as compared with prostate cells 26 . We do find evidence for an increase in TMPRSS2 transcript but not protein expression with smoking. This difference in findings for transcript vs. protein expression might be attributable to differences in the biospecimens sampled across studies or differential expression in cellular subtypes.

In contrast to TMPRSS2, we do find evidence for sexually dimorphic ACE2 expression. Specifically, protein expression is higher in the lungs of male smokers compared with female smokers and in male mice compared with female mice. Others studies have suggested a slight male bias in lung ACE2 expression 27 . Importantly, source tissues and/or methods can account for some differences among studies. In terms of pharmacologic intervention, suppression of ACE2 expression with a potent AR antagonist in the lungs of male mice is significant albeit appears to be modest. Whether this apparent magnitude of androgen regulation of ACE2 expression is meaningful and accounts for the sex-associated differences in SARS-CoV-2 infection or COVID-19 severity is unclear.

Nevertheless, this will be tested in several ongoing clinical trials (including ClinicalTrials.gov NCT04475601, NCT04397718, NCT045009999 and NCT04374279).

Sexually dimorphic AR expression in an organ not associated with sexual differentiation-the male lungraises the question of function. This is reminiscent of sexually dimorphic AR protein expression in the male human kidney, in which AR function includes regulating glucocorticoid metabolism and downstream steroid receptor activity 28 . Whether this physiology also occurs in the lung has yet to be determined. This may have implications for regulation of inflammatory processes including asthma.

In conclusion, we find no evidence for androgen regulation of TMPRSS2 in the male lung. Therefore, TMPRSS2 regulation in the lung appears to fundamentally differ from clear androgen-dependent regulation in prostatic tissues. In contrast, there is a sex discordance in AR and ACE2 expression in lungs of mice and humans. In humans, elevated ACE2 expression is apparent in the lungs of male smokers compared with female smokers. Pulmonary TMPRSS2 regulation appears not to account for the sex-discordance in COVID-19 clinical outcomes. In contrast, ACE2 expression is elevated in smokers and particularly in males. The magnitude of ACE2 suppression with enzalutamide in mouse is modest. The ultimate effects of anti-androgens on human pulmonary ACE2 expression and COVID-19 outcomes are not yet known. Smokers could partly mitigate their increased risk by quitting smoking. 

Mice, treatments and lung harvest. A cohort of adult NSG mice (> 6 weeks old) were obtained from Cleveland Clinic Biological Resources Unit. The male mice were arbitrarily divided between two groups receiving control chow or Enz diet 62.5 mg/kg. Following 11 days on diet, the mice were sacrificed using a lethal dose of Nembutal followed by cardiac puncture. Once sacrificed, the abdominal and thoracic cavities of the mice were opened, the inferior vena cava was cut, and the lungs were gently perfused with warm saline via the right ventricle. Next, the lungs were removed, and the individual lobes were either fixed in 10% formalin, or embedded in paraffin, or snap frozen for subsequent RNA or protein analysis. Total RNA was harvested by homogenizing 25 mg lung tissue in 350 µl RLT buffer (RNeasy kit, Qiagen) following the manufacturer's instructions. cDNA synthesis were then carried out with the iScript cDNA Synthesis Kit (Bio-Rad). Quantitative PCR (qPCR) analysis was conducted in triplicate in an ABI 7500 Real-Time PCR machine (Applied Biosystems) using iTaq Fast SYBR Green Supermix with ROX (Bio-Rad) with the following primer sets:

Tmprss2-forwardgtcatccacacacatcccaagtc; reversetcccagaacctccaaagcaaga; Ace2-forwardactatgaagcagagggagcagatg; reverseggctgatgtaggaagggtaggtat; Ar-forward ggcagcagtgaagcaggtag; reversecggacagagccagcggaa; Rplp0 -forward gacctccttcttccaggctttg; reversectcccaccttgtctccagtcttta; TMPRSS2-forward atcggtgtgttcgcctctacg; reverse atccgctgtcatccactattcctt; ACE2-forward ggaggatgtgcgagtggcta; reversetaggctgttgtcattcagacgg; RPLP0 -forward attacaccttcccacttgctg; reverse actcttccttggcttcaacctta.

ChIP-seq. ChIP-Seq analysis was performed in LNCaP cells treated with vehicle (DMSO) or Enz 10 μM.

Briefly, 106 cells were cross-linked using 1% formaldehyde (reconstituted in 1X PBS) at room temperature for 10 min, followed by quenching with Glycine (final concentration 125 mM) and further incubation at room temperature for 5 min. Fixed cells were lysed in SDS buffer (50 mM Tris-HCl pH 8, 1% SDS, 10 mM EDTA), and sonicated at 4 degrees using Bioruptor sonicator (Diagenode) for 20 cycles on high setting: 10″ on/10″ off. The lysates were then immunoprecipitated in ChIP Dilution Buffer (20 mM Tris-HCl pH 8, 1% Triton-X 100, 2 mM EDTA, 150 mM NaCl + PIC), using 3 μg chromatin, 10 μl of anti-AR (Millipore, EMD Millipore: PG-21), and 50 μl of blocked A/G beads. Recovered DNA was used to prepare libraries using the Illumina Nextera library prep method, subsequently sequenced on a NextSeq 500 and analyzed using the ChiLin analytical pipeline. Finally, genomic AR occupancy was visualized and compared using the Integrated Genome Viewer.

Immunohistochemistry. Immunohistochemistry staining was performed using the Discovery ULTRA automated stainer from Roche Diagnostics (Indianapolis, IN). In brief, antigen retrieval was performed using a tris/borate/EDTA buffer (Discovery CC1, 06414575001; Roche), pH 8.0 to 8.5, at 95 °C for 32 min. For AR staining only, 64 min of antigen retrieval time was applied. The slides were then incubated with primary antibodies for 1 h at room temperature with the following dilutions: TMPRSS2 (ab92323), 1:3000; TMPRSS2 (ab214462), 1:200; Androgen Receptor (ab133273), 1:100; ACE2 (R&D Systems, #AF933), 1:400. The antibodies were visualized using the OmniMap anti-Rabbit HRP (05269679001; Roche), and OmniMap anti-Goat HRP (06607233001; Roche) in conjunction with the ChromoMap DAB detection kit (05266645001; Roche). Lastly, the slides were counterstained with hematoxylin and bluing. The specificity of each antibody was first tested on appropriate Figure 5 . Expression of both TMPRSS2 and ACE2 transcript in human bronchial epithelia increases in current smokers compared to both former and never smokers and the increases do not depend on smoking pack years. (A) Box plot of TMPRSS2 expression (normalized; mean within each data set equals 1) from GEO data sets. N for each data set: GSE16008 = 13 current/13 never smoker; GSE18385 large airway = 32 current/20 never smoker; GSE18385 small airway = 58 current/51 never smoker; GSE37147 = 99 current/139 former smoker; GSE61327 = 112 current/71 never smoker; GSE7895 = 52 current/31 former/21 never smoker; GSE994 = 34 current/18 former/23 never smoker. Data sets 18,385, 63,127, 7895, and 994 contained multiple TMPRSS2 reference sequences whose expression values were summed. In all data sets containing four reference sequences (18,385 large and small airway and 63,127), the difference between groups was largest for sequence AI660243, so expression using that sequence alone is also shown. For data sets with two groups, p-values were obtained from t-tests. For data sets with three groups, Tukey HSD p-values were obtained after one-way ANOVA (ANOVA p-values: GSE7895 0.009, GSE994 0.033). (B) Box plot of ACE2 expression from the same data sets as in (A). Data sets 18,385 and 63,127 contained multiple ACE2 reference sequences whose expression values were summed. For data sets with two vs. three groups, p-values were obtained as in (A) (ANOVA p-values for three group sets: GSE7895 0.82, GSE994 0.010). (C) Scatter plots of normalized TMPRSS2 expression vs. smoking pack years for current and former smokers in data sets containing pack year data. GSE37147 (former smokers) had adjusted R 2 = 0.03 and p = 0.03; no other linear regression had p < 0. www.nature.com/scientificreports/ control tissues (Fig. S3 ) before proceeding to staining of the lung sections. Controls without primary antibody were included in all experiments (Fig. S4 ).

Gene expression in human lung. The public genomics data repository Gene Expression Omnibus (GEO, ncbi.nlm.nih.gov/geo) was searched for data sets containing expression profiling of samples from non-cancerous human lung and bronchial/airway epithelial cells with samples identified by gender and/or smoking status of subjects. The following data sets were identified: GSE994 (airway epithelial cells from current/former/never smokers), GSE4115 (histologically normal bronchial epithelial cells from smokers with and without lung cancer), GSE7895 (airway epithelial cells from current/former/never smokers), GSE16008 (bronchial epithelial cells from healthy current and never smokers), GSE18385 (large and small airway epithelial cells from healthy current and never smokers), GSE37147 (bronchial epithelial cells from current and former smokers with and without COPD), GSE43696 (bronchial epithelial cells from asthma patients and healthy controls), GSE63127 (small airway epithelial cells from healthy current and never smokers), GSE103174 (lung tissue from smokers and nonsmokers with and without COPD), and GSE123352 (non-involved lung parenchyma from ever and never smokers with lung adenocarcinoma). TMPRSS2, AR, and ACE2 gene expression values were obtained from each data set and analyses for comparisons between groups (for each data set for which gender and/or current smoking status information was available) were performed using R.

Statistics. For comparisons involving more than two groups, ANOVA with post hoc tests as indicated in figure legends was performed, or Kruskal-Wallis with post hoc tests for IHC results. For comparisons between two groups, t-tests were performed. For correlations between gene expression and smoking pack years, linear models were fit to data. All statistical tests were two-sided. Analyses were performed in R or GraphPad Prism.

Study approval. All methods on animals were performed in accordance with a protocol approved by the Institutional Animal Care and Use Committee (IACUC) of the Cleveland Clinic Lerner Research Institute (Protocol number: 2018-2093), the relevant guidelines and regulations, and are compliant with the ARRIVE guidelines. Studies using human tissues are deidentified and were deemed to be exempt by the Cleveland Clinic Institutional Review Board (IRB) and are consistent with the Declaration of Helsinki. Therefore, informed consent was not required.

All data generated or analyzed during this study are included in this published article.

Received: 10 November 2020; Accepted: 29 April 2021

A pneumonia outbreak associated with a new coronavirus of probable bat origin

A new coronavirus associated with human respiratory disease in China

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

Angiotensin-converting enzyme 2 (ACE2) as a SARS-CoV-2 receptor: Molecular mechanisms and potential therapeutic target

Serendipity or opportunity for intervention

Recurrent fusion of TMPRSS2 and ETS transcription factor genes in prostate cancer

Androgen-deprivation therapies for prostate cancer and risk of infection by SARS-CoV-2: A population-based study (N = 4532)

Androgen deprivation therapy in men with prostate cancer does not affect risk of infection with SARS-CoV-2

Baseline characteristics and outcomes of 1591 patients infected with SARS-CoV-2 admitted to ICUs of the Lombardy region

Why are so many more men dying from coronavirus

Androgen signaling regulates SARS-CoV-2 receptor levels and is associated with severe COVID-19 symptoms in men

Androgen receptor and androgen-dependent gene expression in lung

Androgen hazards with COVID-19

HSD3B1 genotype identifies glucocorticoid responsiveness in severe asthma

A hierarchical network of transcription factors governs androgen receptor-dependent prostate cancer growth

A 'quickscore' method for immunohistochemical semiquantitation: Validation for oestrogen receptor in breast carcinomas

A new immunohistochemical antibody for the assessment of estrogen receptor status on routine formalin-fixed tissue samples

Sex steroid hormone receptors in human thymoma

Specificity of ligand-dependent androgen receptor stabilization: Receptor domain interactions influence ligand dissociation and receptor stability

Phenotypic analysis of mice lacking the Tmprss2-encoded protease

The impact of COPD and smoking history on the severity of Covid-19: A systemic review and meta-analysis

Clinical characteristics of coronavirus disease 2019 in China

Smoking upregulates angiotensin-converting enzyme-2 receptor: A potential adhesion site for novel coronavirus SARS-CoV-2 (Covid-19)

ACE-2 expression in the small airway epithelia of smokers and COPD patients: Implications for COVID-19

Cigarette smoke exposure and inflammatory signaling increase the expression of the SARS-CoV-2 receptor ACE2 in the respiratory tract

Distinct mechanisms for TMPRSS2 expression explain organ-specific inhibition of SARS-CoV-2 infection by enzalutamide

Landscape of X chromosome inactivation across human tissues

Deep androgen receptor suppression in prostate cancer exploits sexually dimorphic renal expression for systemic glucocorticoid exposure

This work is supported in part by grants from the Prostate Cancer Foundation and the National Cancer Institute (R01CA172382 and R01CA236780) and the National Heart Lung and Blood Institute (HL60917).

The authors declare no competing interests.

The online version contains supplementary material available at https:// doi. org/ 10. 1038/ s41598-021-90491-1.Correspondence and requests for materials should be addressed to N.S.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.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:// creat iveco mmons. org/ licen ses/ by/4. 0/.