Genomic epidemiology of the main SARS-CoV-2 variants in Italy between summer 2020 and winter 2021.

Since the beginning of the pandemic, SARS-CoV-2 has shown a great genomic variability, resulting in the continuous emergence of new variants that has made their global monitoring and study a priority. This work aimed to study the genomic heterogeneity, the temporal origin, the rate of viral evolution and the population dynamics of the main circulating variants (20E.EU1, Alpha and Delta) in Italy, in August 2020-January 2022 period. For phylogenetic analyses, three datasets were set up, each for a different main lineage/variant circulating in Italy in that time including other Italian and International sequences of the same lineage/variant, available in GISAID sampled in the same times. The international dataset showed 26 (23% Italians, 23% singleton, 54% mixed), 40 (60% mixed, 37.5% Italians, 1 singleton) and 42 (85.7% mixed, 9.5% singleton, 4.8% Italians) clusters with at least one Italian sequence, in 20E.EU1  clade, Alpha and Delta variants, respectively. The estimation of tMRCAs in the Italian clusters (including >70% of genomes from Italy) showed that in all the lineage/variant, the earliest clusters were the largest in size and the most persistent in time and frequently mixed. Isolates from the major Italian Islands tended to segregate in clusters more frequently than those from other part of Italy. The study of infection dynamics showed a positive correlation between the trend in the effective number of infections estimated by BSP model and the Re curves estimated by birth-death skyline plot. The present work highlighted different evolutionary dynamics of studied lineages with high concordance between epidemiological parameters estimation and phylodynamic trends suggesting that the mechanism of replacement of the SARS-CoV-2 variants must be related to a complex of factors involving the transmissibility, as well as the implementation of control measures, and the level of cross-immunization within the population.

No need to modify treatment within the first month after rapid start of a tailored antiretroviral therapy: the TWODAY Study.

The aim of the TWODAY Study was to investigate the frequency of early treatment change after rapid start of a tailored ART regimen (a 2-drug regimen - 2DR, when clinically feasible or a 3-drug regimen - 3DR, otherwise). TWODAY was an open-label, prospective, proof-of-concept, single center study. ART-naïve patients started their first-line regimen within a few days from the first laboratory testing with a 2DR of dolutegravir (DTG) and lamivudine (3TC) if CD4+ count >200 cells/mL, HIVRNA <500,000 copies/mL, no transmitted drug resistance to DTG or 3TC and HBsAg undetectable; otherwise, ART was started with a 3DR. The primary endpoint was the proportion of patients who needed to change ART within four week from start, for any reason. Thirty-two patients were enrolled; 19 (59.3%) were deemed eligible for a 2DR. Median time from laboratory testing to ART start was 5 days (5; 5). No regimen modification occurred within one month. In conclusion, no regimen modification was needed within the first month of treatment. Starting a 2DR within a few days after HIV diagnosis was feasible, relying upon complete results of the needed laboratory tests (including resistance testing). A 2DR can be safely proposed provided full laboratory tests are readily available.

Detection of low-level HCV variants in DAA treated patients: comparison amongst three different NGS data analysis protocols.

BACKGROUND: Notwithstanding the efforts of direct-acting antivirals (DAAs) for the treatment of chronically infected hepatitis C virus (HCV) patients, concerns exist regarding the emergence of resistance-associated substitutions (RAS) related to therapy failure. Sanger sequencing is still the reference technique used for the detection of RAS and it detects viral variants present up to 15%, meaning that minority variants are undetectable, using this technique. To date, many studies are focused on the analysis of the impact of HCV low variants using next-generation sequencing (NGS) techniques, but the importance of these minority variants is still debated, and importantly, a common data analysis method is still not defined. METHODS: Serum samples from four patients failing DAAs therapy were collected at baseline and failure, and amplification of NS3, NS5A and NS5B genes was performed on each sample. The genes amplified were sequenced using Sanger and NGS Illumina sequencing and the data generated were analyzed with different approaches. Three different NGS data analysis methods, two homemade in silico pipeline and one commercially available certified user-friendly software, were used to detect low-level variants. RESULTS: The NGS approach allowed to infer also very-low level virus variants. Moreover, data processing allowed to generate high accuracy data which results in reduction in the error rates for each single sequence polymorphism. The results improved the detection of low-level viral variants in the HCV quasispecies of the analyzed patients, and in one patient a low-level RAS related to treatment failure was identified. Importantly, the results obtained from only two out of the three data analysis strategies were in complete agreement in terms of both detection and frequency of RAS. CONCLUSIONS: These results highlight the need to find a robust NGS data analysis method to standardize NGS results for a better comprehension of the clinical role of low-level HCV variants. Based on the extreme importance of data analysis approaches for wet-data interpretation, a detailed description of the used pipelines and further standardization of the in silico analysis could allow increasing diagnostic laboratory networking to unleash true potentials of NGS.

Gonadal Function Recovery and Fertility in Women Treated with Chemo- and/or Radiotherapy for Hodgkin's and Non-Hodgkin Lymphoma.

BACKGROUND: Fertility and gonadal function represent one of the most important aspects for long-term lymphoma survivors. AIMS: The aim of our study was to determine possible risk factors, such as age at treatment, chemotherapeutic regimen, protection with oral contraceptives (OCs), and gonadotropin-releasing hormone (GnRH) analogues in female patients treated for Hodgkin's lymphoma (HL) or non-Hodgkin lymphoma (NHL) at a reproductive age. METHODS: Patients between the age of 16 and 50 years at the time of HL or NHL diagnosis were selected. Eligible patients were requested to respond to a questionnaire by phone interview about fertility, menstrual status, sexual aspects, and treatment with OCs or GnRH analogues during chemotherapy. RESULTS: The resumption of menstrual activity was associated with the use of the OCs and GnRH analogues during chemotherapy (p = 0.008 and 0.034, respectively). At univariate analysis, the use of OCs during chemotherapy was associated with a lower risk of amenorrhea (prevalence ratio [PR] = 0.37; 95% CI 0.17-0.82). A higher age at the time of treatment correlated positively with therapy-induced amenorrhea, with a difference of 12.8 years between the mean age at diagnosis of the women with therapy-induced amenorrhea and those who resumed their menses. Amenorrhea was significantly higher in women receiving R-CHOP than in women treated with ABVD (PR = 6.00; 95% CI 2.32-15.54). Moreover, NHL had an infertility PR of 1.51 (95% CI 0.86-2.45) at multivariate analysis compared to HL. CONCLUSIONS: This study suggests a possible role of pharmacological prophylaxis with OCs and GnRH analogues.

48 week outcomes of maraviroc-containing regimens following the genotypic or Trofile assay in HIV-1 failing subjects: the OSCAR Study.

This study assessed the 48-week efficacy of an antiretroviral therapy including maraviroc following the assessment of co-receptor tropism by use of Geno2Pheno algorithm or the Trofile phenotypic assay in failing treatment-experienced HIV-1 patients. This was a multicenter, randomized, open-label, non-inferiority trial. Treatment-experienced subjects with HIV-RNA ≥500 copies/mL were randomized (1:1) to undergo co-receptor tropism testing by the Geno- 2Pheno algorithm (with a false positive rate >10%) or the Trofile assay before starting a new antiretroviral treatment which included maraviroc. The primary endpoint was the 48 week proportion of patients with treatment success (TS). Intention-to-treat analyses are also reported. One hundred and fifty-five experienced patients were analysed: 77 patients in the Trofile arm and 78 in the Genotype arm. The 48-week proportion of TS was 87% in the Trofile arm and 89% in the Genotype arm (difference: 1.5%, 95%CI: -8.9% to 11.8%) suggesting non-inferiority. In the Trofile arm, 10 patients had treatment failure: 5 viral rebound, 5 discontinuations. In the Genotype arm, 9 patients had treatment failure: 7 viral rebound, 2 lost to follow-up. CD4+ significantly increased from baseline to week 48 in both arms. 48-week treatment success was similar for maraviroc-including therapy prescribed following the Trofile phenotypic assay or Geno2Pheno algorithm.

Performance of commonly used genotypic assays and comparison with phenotypic assays of HIV-1 coreceptor tropism in acutely HIV-1-infected patients.

OBJECTIVES: Although founder viruses in primary HIV-1 infections (PHIs) typically use the CCR5 coreceptor (R5-tropic), 3%-19% of subjects also harbour CXCR4-using viruses (X4-tropic), making tropism determination before CCR5 antagonist usage mandatory. Genotypic methods can be used to accurately determine HIV-1 tropism in chronically infected patients. METHODS: We compared the results of genotypic methods [geno2pheno, PSSMx4r5 including a novel nucleotide-input version (ntPSSM) and distant segments (ds)Kernel] to predict coreceptor usage in a cohort of 67 PHIs. Specimens with discrepant results were phenotypically tested after cloning the V3 gene region into proviral backbones. Recombinant viruses were used to infect U87 indicator cell lines bearing CD4 and either CCR5 or CXCR4. RESULTS: Geno2pheno10%, PSSMx4r5 and (ds)Kernel gave identical predictions in 85% of cases. Geno2pheno10% predicted the presence of CXCR4 viruses in 18% of patients. Two patients were predicted to carry X4-tropic viruses by all algorithms and X4-tropic viruses were detected in at least one of the recombinant AD8 or NL4-3 backbone-based assays. Ten samples resulted in discordant predictions with at least one algorithm. Full concordance between tropism prediction by using population sequencing and phenotypic assays was observed only with ntPSSM. Geno2pheno prediction and the phenotypic assay gave the same results in a minority of 'discordant' patients. CONCLUSIONS: Compared with both PSSMx4r5 versions, (ds)Kernel and our phenotypic assay, geno2pheno10% overestimated the frequency of X4-tropic viruses (18% versus 3%). ntPSSM was able to detect one additional X4 virus compared with (ds)Kernel that was confirmed with the phenotypic assay.

2-Aminothiazolones as anti-HIV agents that act as gp120-CD4 inhibitors.

We report here the synthesis of 2-aminothiazolones along with their biological properties as novel anti-HIV agents. Such compounds have proven to act through the inhibition of the gp120-CD4 protein-protein interaction that occurs at the very early stage of the HIV-1 entry process. No cytotoxicity was found for these compounds, and broad antiviral activities against laboratory strains and pseudotyped viruses were documented. Docking simulations have also been applied to predict the mechanism, at the molecular level, by which the inhibitors were able to interact within the Phe43 cavity of HIV-1 gp120. Furthermore, a preliminary absorption, distribution, metabolism, and excretion (ADME) evaluation was performed. Overall, this study led the basis for the development of more potent HIV entry inhibitors.

Clinico-Virological Outcomes and Mutational Profile of SARS-CoV-2 in Adults Treated with Ribavirin Aerosol for COVID-19 Pneumonia.

The emergence of new SARS-CoV-2 variants can affect vaccine efficacy, laboratory diagnosis and the therapies already available, triggering interest in the search for antiviral agents for SARS-CoV-2 infections. Ribavirin (RBV) is a broad-spectrum antiviral with demonstrated in vitro activity against multiple viruses, including SARS-CoV-2. This retrospective study evaluated the dynamics and viral clearance of SARS-CoV-2 in hospitalised adult participants (PTs) with COVID-19 pneumonia who received an RBV aerosol within a compassionate use study. The impact of RBV on the clinical outcome and the mutational profile of SARS-CoV-2 was also assessed. The median RNA values measured in nine PTs included in this study decreased from baseline to discharge (at BL, threshold cycle (Ct) = 22.4, IQR 19.84-5.07; at discharge, Ct = 27.92, IQR 26.43-36.11), with a significant decline in the Ct value evaluated by Friedman rank ANOVA analysis, p = 0.032. Seven out of nine PTs experienced a clinical improvement, while two PTs deceased during hospitalisation. In PTs with a favourable outcome, the virus clearance rate at discharge was 28.6%. The cumulative clearance rate was 71.4% within 14 days from discharge. A mutational pattern after RBV was detected in three out of five PTs in whom whole-genome sequencing was available. Our findings suggest that RBV limits SARS-CoV-2 replication, possibly resulting in a favourable clinical outcome. Ribavirin may also contribute to the mutational spectrum of SARS-CoV-2.

Infection and coinfection of human rhinovirus C in stem cell transplant recipients.

In 54 adult stem cell transplant recipients, the presence and persistence of human rhinoviruses (including the novel lineage C) were evaluated by molecular detection and phylogenetic analysis, independently from respiratory symptoms. In the same group of patients, the presence of other coinfecting respiratory pathogens, including the novel enterovirus 109, was also evaluated.

Case Report: Favorable outcome of allogeneic hematopoietic stem cell transplantation in SARSCoV2 positive recipient, risk-benefit balance between infection and leukemia.

Allogeneic hematopoietic stem cell transplantation (allo-HSCT) in SARS-CoV-2 positive candidates is usually delayed until the clinical resolution of the infection's symptoms and a negative nasopharyngeal molecular test. However, prolonged SARS-CoV-2 positivity has been frequently observed in haematological malignancies, thus representing a challenge for the timing of transplant procedures. Here, we report on the case of a 34-year-old patient with recent pauci-symptomatic COVID-19 undergoing transplant for high-risk acute B-lymphoblastic leukemia before achieving viral clearance. Shortly before their scheduled allogeneic HSCT from a matched unrelated donor, the patient developed mild Omicron BA.5 infection receiving nirmatrelvir/ritonavir with fever resolution within 72 hours. Twenty-three days after COVID-19 diagnosis, because of increasing minimal residual disease values in the context of high-risk refractory leukemia and clinical resolution of SARS-2-CoV infection with reduction of viral load at surveillance nasopharyngeal swabs, it was decided not to delay further allo-HSCT. During myelo-ablative conditioning, the nasopharyngeal SARS-CoV-2 viral load increased while the patient remained asymptomatic. Consequently, two days before the transplant, intra-muscular tixagevimab/cilgavimab 300/300 mg and a 3-day course of intravenous remdesivir were administered. During the pre-engraftment phase, veno-occlusive disease (VOD) occurred at day +13, requiring defibrotide treatment to obtain a slow but complete recovery. The post-engraftment phase was characterized by mild COVID-19 at day +23 (cough, rhino-conjunctivitis, fever) that spontaneously resolved, achieving viral clearance at day +28. At day +32, she experienced grade I acute graft-versus host disease (a-GVHD, skin grade II) treated with steroids and photo-apheresis, without further complications during follow-up until day +180. Addressing the issue of allo-HSCT timing in patients recovering from SARS-CoV-2 infection with high-risk malignant diseases is challenging because of 1] the high risk of COVID-19 clinical progression, 2] the impact of transplant delay on leukemia prognosis and 3] the occurrence of endothelial complications such as VOD, a-GVHD, and transplant associated thrombotic micro-angiopathy. Our report describes the favourable outcome of allo-HSCT in a recipient with active SARS-CoV2 infection and high-risk leukemia thanks to timely anti-SARS-CoV-2 preventive therapies and prompt management of transplant-related complications.

Epidemiological, molecular and clinical features of enterovirus 109 infection in children and in adult stem cell transplant recipients.

BACKGROUND: A novel human enterovirus (HEV) type within the species HEV-C, named EV109, was discovered from cases of respiratory illness in Nicaragua in September 2010. The aim of this study, was to retrospectively examine the presence and the role of EV109 in respiratory samples from two patients populations; infants below the age of 2 years, hospitalized for acute respiratory diseases (ARDs) and adult hematopoietic stem cell transplantation recipients. RESULTS: A total of 1149 nasopharingeal aspirates were collected and tested for the presence of EV109 by reverse transcription-PCR (RT-PCR). In positive samples, the presence of the most common respiratory viruses was also assayed and clinical symptoms were evaluated. Samples from 2 of the 974 infants tested positive for EV109 RNA (0.2%) and belonged to patients with lower ARDs; co-infection with other viral pathogens under study was observed in both cases. In transplant recipients, one out of the 175 samples analyzed, from a patients with upper respiratory simptoms tested positive for HEV 109 in the absence of co-infecting viruses. Sequence analysis of amplified EV109 genomic regions, showed only a few nucleotide differences when compared with the Nicaraguan strains. CONCLUSIONS: Overall these results indicate that HEV109 variants have circulated and differentiated in different lineages worldwide. Although more cases and larger studies are needed, HEV109 infection may be associated to ARDs both in infants and in hematopoietic stem cell transplantation recipients. If these preliminary observations will be confirmed, improved molecular methods with a wider panel of potential pathogens will be useful for monitoring these categories of patients.

Neutrophils predominate the immune signature of cerebral thrombi in COVID-19 stroke patients.

Coronavirus disease 2019 (COVID-19) is associated with an increased risk of thrombotic events. Ischemic stroke in COVID-19 patients entails high severity and mortality rates. Here we aimed to analyze cerebral thrombi of COVID-19 patients with large vessel occlusion (LVO) acute ischemic stroke to expose molecular evidence for SARS-CoV-2 in the thrombus and to unravel any peculiar immune-thrombotic features. We conducted a systematic pathological analysis of cerebral thrombi retrieved by endovascular thrombectomy in patients with LVO stroke infected with COVID-19 (n = 7 patients) and non-covid LVO controls (n = 23). In thrombi of COVID-19 patients, the SARS-CoV-2 docking receptor ACE2 was mainly expressed in monocytes/macrophages and showed higher expression levels compared to controls. Using polymerase chain reaction and sequencing, we detected SARS-CoV-2 Clade20A, in the thrombus of one COVID-19 patient. Comparing thrombus composition of COVID-19 and control patients, we noted no overt differences in terms of red blood cells, fibrin, neutrophil extracellular traps (NETs), von Willebrand Factor (vWF), platelets and complement complex C5b-9. However, thrombi of COVID-19 patients showed increased neutrophil density (MPO+ cells) and a three-fold higher Neutrophil-to-Lymphocyte Ratio (tNLR). In the ROC analysis both neutrophils and tNLR had a good discriminative ability to differentiate thrombi of COVID-19 patients from controls. In summary, cerebral thrombi of COVID-19 patients can harbor SARS-CoV2 and are characterized by an increased neutrophil number and tNLR and higher ACE2 expression. These findings suggest neutrophils as the possible culprit in COVID-19-related thrombosis.

Characterization of a Lineage C.36 SARS-CoV-2 Isolate with Reduced Susceptibility to Neutralization Circulating in Lombardy, Italy.

SARS-CoV-2 spike is evolving to maximize transmissibility and evade the humoral response. The massive genomic sequencing of SARS-CoV-2 isolates has led to the identification of single-point mutations and deletions, often having the recurrence of hotspots, associated with advantageous phenotypes. We report the isolation and molecular characterization of a SARS-CoV-2 strain, belonging to a lineage (C.36) not previously associated with concerning traits, which shows decreased susceptibility to vaccine sera neutralization.

Genotypic/phenotypic patterns of HIV-1 integrase resistance to raltegravir.

OBJECTIVES: To understand the dynamic viral evolution observed during failure on raltegravir-containing regimens, we studied the genotypic and phenotypic patterns of resistance to raltegravir and the residual replication capacity (rRC) of HIV-1 variants selected in vivo. METHODS: Clonal genotypic analyses were performed on sequential HIV-1 integrase sequences amplified from 11 failing patients and sampled every 4-24 weeks for up to 64 weeks. Fully replicating recombinant viruses were generated using modified vectors in which selected viral integrase genes amplified from patients' plasma were cloned. rRC was measured by a novel multiple cycle competition assay. Resistance to raltegravir and the rRC of resistant HIV-1 variants selected in vivo were evaluated in purified CD4+ T cells. RESULTS: In all of the patients but one, failure was associated with selection of mutations in positions 143, 148 or 155. Unlike mutations at position 143 (Y143S/K/R), identified alone or in combination with others, mutations at position 148 and 155 were always found in combination. A wide range of resistance levels to raltegravir [from 10- to 770-fold change in 50% inhibitory concentration (IC(50)) compared with baseline] was observed using recombinant viral clones. Finally, rRC was not significantly altered in highly resistant variants. DISCUSSION: Two patterns of viral evolution were observed in the resistant viral populations, driving the variants towards a fast (most of them with G140S + Q148H mutations) or progressive increase in resistance to raltegravir. These results may have implications either for the evaluation of genotypic results, or for the correct clinical use of the compound.

Molecular epidemiology of KI and WU polyomaviruses in infants with acute respiratory disease and in adult hematopoietic stem cell transplant recipients.

Polyomaviruses KI (KIPyV) and WU (WUPyV) were described recently in children with acute respiratory disease. The pathogenic potential of these human viruses has not been determined completely, but a correlation between immunosuppression and virus reactivation has been suggested. In the present study, the association between KI/WUPyV infection and immunosuppression was investigated using sequential nasopharyngeal aspirates from asymptomatic adult hematopoietic stem cell transplant recipients. In parallel, an investigation on the WU/KIPyV prevalence in children with acute respiratory disease was also carried out. Two of the 126 samples obtained from the 31 hematopoietic transplant recipients were positive for KIPyV (1 sample, 0.79%) and WUPyV (1 sample, 0.79%). Both samples were obtained 15 days after allogeneic transplantation and virus persistence was not observed in subsequent samples. In symptomatic children, 7 of the 486 nasopharyngeal aspirates were positive for WUPyV (1.4%) and 1 for KIPyV (0.2%). Single polyomavirus infection was detected in four patients, whereas the remaining patients were co-infected with respiratory syncityal virus (three patients) or adenovirus (one patient). The results suggest that WU/KIPyVs have a limited circulation in Italy and a low pathogenic potential in young children. Brief and asymptomatic infection can occur in hematopoietic transplant recipients.

Dynamic features of the selective pressure on the human immunodeficiency virus type 1 (HIV-1) gp120 CD4-binding site in a group of long term non progressor (LTNP) subjects.

The characteristics of intra-host human immunodeficiency virus type 1 (HIV-1) env evolution were evaluated in untreated HIV-1-infected subjects with different patterns of disease progression, including 2 normal progressor [NP], and 5 Long term non-progressor [LTNP] patients. High-resolution phylogenetic analysis of the C2-C5 env gene sequences of the replicating HIV-1 was performed in sequential samples collected over a 3-5 year period; overall, 301 HIV-1 genomic RNA sequences were amplified from plasma samples, cloned, sequenced and analyzed. Firstly, the evolutionary rate was calculated separately in the 3 codon positions. In all LTNPs, the 3rd codon mutation rate was equal or even lower than that observed at the 1st and 2nd positions (p = 0.016), thus suggesting strong ongoing positive selection. A Bayesian approach and a maximum-likelihood (ML) method were used to estimate the rate of virus evolution within each subject and to detect positively selected sites respectively. A great number of N-linked glycosylation sites under positive selection were identified in both NP and LTNP subjects. Viral sequences from 4 of the 5 LTNPs showed extensive positive selective pressure on the CD4-binding site (CD4bs). In addition, localized pressure in the area of the IgG-b12 epitope, a broad neutralizing human monoclonal antibody targeting the CD4bs, was documented in one LTNP subject, using a graphic colour grade 3-dimensional visualization. Overall, the data shown here documenting high selective pressure on the HIV-1 CD4bs of a group of LTNP subjects offers important insights for planning novel strategies for the immune control of HIV-1 infection.

Persistent replication of severe acute respiratory syndrome coronavirus in human tubular kidney cells selects for adaptive mutations in the membrane protein.

Severe acute respiratory syndrome (SARS) is a systemic disease characterized by both lung pathology and widespread extrapulmonary virus dissemination causing multiple organ injuries. In this regard, renal dysfunction is an ominous sign in patients with SARS. Indeed, clusters of SARS coronavirus (SARS-CoV) particles have been detected in the cytoplasm of renal tubular epithelial cells in postmortem studies, explaining the presence of infectious virus in the urine of SARS patients. In order to investigate the potential SARS-CoV kidney tropism, we have evaluated the susceptibility of human renal cells of tubular and glomerular origin to in vitro SARS-CoV infection. Immortalized cultures of differentiated proximal tubular epithelial cells (PTEC), glomerular mesangial cells (MC), and glomerular epithelial cells (podocytes) were found to express the SARS-CoV receptor angiotensin-converting enzyme 2 on their surface. Productive infection, however, occurred only in PTEC but not in glomerular cells. A transient infection with poor virus production was observed in MC, whereas podocytes were not permissive to SARS-CoV infection. In contrast to the cytopathic infection of the Vero E6 cell line, SARS-CoV did not cause overt cytopathic effects in PTEC or MC. Of interest, PTEC, but not MC, maintained stable levels of SARS-CoV production in serial subcultures, suggesting a persistent state of infection. In this regard, a SARS-CoV variant with increased replication capacity in PTEC was selected after four serial subculture passages. This SARS-CoV variant acquired a single nonconservative amino acid change from glutamic acid (E) to alanine (A) at position 11 in the viral membrane (M) protein. The E11A point mutation was sufficient for enhanced SARS-CoV replication and persistence in PTEC when introduced in a SARS-CoV recombinant infectious clone. These findings indicate that human PTEC may represent a site of SARS-CoV productive and persistent replication favoring the emergence of viral variants with increased replication capacity, at least in these kidney cells.