The Impact of COVID-19 on Treatment Practices for Patients With Early Breast Cancer: A Cross-Sectional Study From a Large Cancer Center in Italy.

INTRODUCTION: The Coronavirus Disease 2019 (COVID-19) has disrupted health services worldwide. The evidence on the impact of the pandemic on cancer care provision, however, is conflicting. We aimed to audit the management of patients diagnosed with early breast cancer (EBC) during the pandemic in a large, tertiary-level cancer center in Italy. METHODS: We conducted a cross-sectional study to track the route to first treatment for patients diagnosed with EBC during 2019, 2020, and 2021. We abstracted data for all consecutive patients referred to the Veneto Institute of Oncology (Padua, Italy). We defined as point of contact (POC) the date of the first consultation with a breast cancer specialist of the breast unit. First treatment was defined as either upfront surgery or neoadjuvant chemotherapy (NACT). RESULTS: We reviewed medical records for 878 patients for whom an MDT report during 2019-2021 (April through June) was available. Of these, 431 (49%) were eligible. The proportion of screen-detected tumors was larger in 2019 and 2021 than in 2020 (59%). Conversely, the proportion of screen-detected tumors was offset by the proportion of palpable tumors in 2020 (P = .004). Distribution of tumor and nodal stage was unchanged over time, but in situ tumors were slightly fewer in 2020 than in 2019 or 2021. The adjusted odds ratio for treatment delay (45 days or more) was 0.87 for 2020 versus 2019 (95% CI, 0.5-1.53) and 0.9 for 2021 versus 2019 (95% CI, 0.52-1.55). CONCLUSIONS: There was no evidence for major changes in the management of patients with EBC during 2019-2021 and no treatment delays were observed. Our findings suggest that more women presented with palpable nodules at diagnosis, but the stage distribution did not change over time. Validation on a larger cohort of patients is warranted to robustly assess the impact of the COVID-19 pandemic on treatment practices for patients with EBC.

Characterization of Gut Microbiome Composition in Patients with Triple-Negative Breast Cancer Treated with Neoadjuvant Chemotherapy.

INTRODUCTION: Patients with triple-negative breast cancer (TNBC) achieving a pathological complete response (pCR) after neoadjuvant chemotherapy have a better event-free survival. The role of gut microbiome in early TNBC is underexplored. METHODS: Microbiome was analyzed by 16SrRNA sequencing. RESULTS: Twenty-five patients with TNBC treated with neoadjuvant anthracycline/taxane-based chemotherapy were included. Fifty-six percent achieved a pCR. Fecal samples were collected before (t0), at 1 (t1), and 8 weeks (t2) from chemotherapy. Overall, 68/75 samples (90.7%) were suitable for microbiome analysis. At t0, pCR group showed a significantly higher α-diversity as compared with no-pCR, (P = .049). The PERMANOVA test on ß-diversity highlighted a significant difference in terms of BMI (P = 0.039). Among patients with available matched samples at t0 and t1, no significant variation in microbiome composition was reported over time. CONCLUSIONS: Fecal microbiome analysis in early TNBC is feasible and deserves further investigation in order to unravel its complex correlation with immunity and cancer.

Validation of Residual Proliferative Cancer Burden as a Predictor of Long-Term Outcome Following Neoadjuvant Chemotherapy in Patients with Hormone Receptor-Positive/Human Epidermal Growth Receptor 2-Negative Breast Cancer.

BACKGROUND: The integration of residual cancer burden (RCB) and post-treatment Ki67 as residual proliferative cancer burden (RPCB) has been proposed as a stronger predictor of long-term outcome in unselected patients with breast cancer (BC) undergoing neoadjuvant chemotherapy (NACT), as compared with RCB. However, no specific analysis in hormone-receptor-positive (HR+) human epidermal growth receptor 2-negative (HER2-) BC is available so far. MATERIALS AND METHODS: A cohort of 130 patients with HR+/HER2- BC who underwent NACT between 2000 and 2014 was included. Archival surgical specimens were evaluated for RCB. RPCB was calculated by combining RCB and Ki67 as previously described. Patients were categorized in four RCB and RPCB categories (pathological complete response and tertiles). Disease-free survival (DFS) and overall survival (OS) estimates were determined by Kaplan-Meier analysis and compared using the log-rank test. Overall change of χ2 and c-indexes were used to compare the performance of the prognostic models. RESULTS: RPCB was calculated for 85 patients. After a median follow up of 8.5 years, RCB was associated with OS (p = .048) but not with DFS (p = .152); RPCB was instead significantly associated with both DFS and OS (p = .034 and p < .001, respectively). In terms of OS, RPCB provided a significant amount of prognostic information beyond RCB (∆χ2 5.73, p < .001). In addition, c-index for OS prediction was significantly higher for RPCB as compared with RCB (0.79 vs. 0.61, p = .03). CONCLUSION: This is the first study evaluating RPCB in patients with HR+/HER2- BC treated with NACT. In this independent cohort, RPCB was a strong predictor of DFS and OS. The better performance of RPCB versus RCB was in part due to the ability of RPCB to discriminate a subgroup of patients with a particularly worse prognosis after NACT, who may be candidates for clinical trials evaluating novel adjuvant strategies. IMPLICATIONS FOR PRACTICE: The present work validated residual proliferative cancer burden (RPCB) as a strong predictor of long-term outcome in patients with hormone receptor-positive human epidermal growth receptor 2-negative (HR+/HER2-) breast cancer (BC) treated with neoadjuvant chemotherapy. In addition, results from the present study suggest RPCB as a promising tool to identify patients with HR+/HER2- BC who might potentially benefit from the inclusion in clinical trials evaluating novel or escalated postneoadjuvant treatment strategies because it allowed to discriminate a subgroup of patients with particularly poor prognosis despite having received subsequent endocrine therapy in the adjuvant setting.

Use of scalp cooling device to prevent alopecia for early breast cancer patients receiving chemotherapy: A prospective study.

Chemotherapy-induced alopecia (CIA) affects the majority of patients receiving chemotherapy (CT) for early breast cancer. It is a highly distressing side effect of CT, with psychological and social impact. Primary aim of the present analysis was to assess the efficacy of scalp cooling with DigniCap® in preventing CIA. Success rate was defined as patients' self-reported hair loss <50% according to Dean scale. In this analysis, we reported success rate at 3 weeks after the first CT course and at 3 weeks after the last CT course. Secondary endpoints included self-reported tolerability and patients' judgment on scalp cooling performance. Consecutive early breast cancer patients admitted to Istituto Oncologico Veneto who were recommended to receive neoadjuvant or adjuvant CT, were eligible to undergo scalp cooling during the CT administration within this study. 135 patients were included: 74% received adjuvant CT and 26% neoadjuvant CT (P < .001). The type of CT was: docetaxel-cyclophosphamide (26%), paclitaxel (23%), epirubicin-cyclophosphamide followed by paclitaxel (32%), and paclitaxel followed by epirubicincyclophosphamide (19%). The rate of success in preventing alopecia was 77% (104/135) at 3 weeks from the start of CT and 60% (81/135) at 3 weeks from the end of treatment. Higher success rates were reported in non-anthracycline (71%) compared to anthracycline-containing CT regimens (54%; P < 0.001). Premature discontinuation of scalp cooling was reported in 29/135 patients (21.5%), including withdrawal for alopecia (16/29), for low scalp cooling tolerability (8/29) or both (5/29). Scalp cooling was generally well tolerated. These results overall suggest that the use of scalp cooling is effective in preventing alopecia in the majority of early breast cancer patients receiving neoadjuvant or adjuvant CT, especially for patients undergoing a taxane-based non-anthracycline regimen.

Whole genome and phylogenetic analysis of two SARS-CoV-2 strains isolated in Italy in January and February 2020: additional clues on multiple introductions and further circulation in Europe.

Whole genome sequences of SARS-CoV-2 obtained from two patients, a Chinese tourist visiting Rome and an Italian, were compared with sequences from Europe and elsewhere. In a phylogenetic tree, the Italian patient's sequence clustered with sequences from Germany while the tourist's sequence clustered with other European sequences. Some additional European sequences in the tree segregated outside the two clusters containing the patients' sequences. This suggests multiple SARS-CoV-2 introductions in Europe or virus evolution during circulation.

Coronavirus disease (COVID-19) in a paucisymptomatic patient: epidemiological and clinical challenge in settings with limited community transmission, Italy, February 2020.

Data concerning the transmission of the novel severe acute respiratory syndrome coronavirus (SARS-CoV-2) in paucisymptomatic patients are lacking. We report an Italian paucisymptomatic case of coronavirus disease 2019 with multiple biological samples positive for SARS-CoV-2. This case was detected using the World Health Organization protocol on cases and contact investigation. Current discharge criteria and the impact of extra-pulmonary SARS-CoV-2 samples are discussed.

Immune characterization of breast cancer metastases: prognostic implications.

BACKGROUND: Tumor-infiltrating lymphocytes (TILs) evaluated in primary breast cancer (BC) convey prognostic information. Limited data in the metastatic setting are available. METHODS: Secondary lesions from 94 BC patients, 43 triple-negative (TN) and 51 HER2-positive, were evaluated for TILs and expression of CD8, FOXP3, and PD-L1 by immunohistochemistry. RESULTS: TILs levels on metastasis were generally low (median 5%) and did not differ between TN and HER2+ tumors. Younger patients showed significantly lower TILs (p = 0.002). In HER2+ patients, TILs were higher in lung metastases as compared to other sites (p = 0.038). TILs composition was different across metastatic sites: skin metastases presented higher FOXP3 (p = 0.002) and lower CD8/FOXP3 ratio (p = 0.032). Patients treated for metastatic BC prior to biopsy had lower CD8 (overall: p = 0.005, HER2+: p = 0.011, TN: p = 0.075). In TN patients, median overall survival (OS) was 11.8 and 62.9 months for patients with low and high TILs, respectively (HR 0.29, 95%CI 0.11-0.76, log-rank p = 0.008). CD8/FOXP3 ratio was also prognostic in TN patients (median OS 8.0, 13.2, and 54.0 months in 1st, 2nd and 3th tertile, log-rank p = 0.019). Both TILs and CD8/FOXP3 ratio were independent factors at multivariate analysis. Counterintuitively, in HER2+ BC, low TILs tumors showed better prognosis (median OS 53.7 vs 39.9 months in TILs low and TILs high, not statistically significant). CONCLUSIONS: Our findings indicate the relevance of TILs as prognostic biomarker for TNBC even in the advanced setting and provide novel hypothesis-generating data on potential sources of immune heterogeneity of metastatic BC.

Prevalence of Usutu and West Nile virus antibodies in human sera, Modena, Italy, 2012.

A collection of 3069 human sera collected in the area of the municipality of Modena, Emilia Romagna, Italy, was retrospectively investigated for specific antibodies against Usutu (USUV) and West Nile viruses (WNV). All the samples resulting positive using a preliminary screening test were analyzed with the plaque reduction neutralization test. Overall, 24 sera were confirmed as positive for USUV (0.78%) and 13 for WNV (0.42%). The results suggest that in 2012, USUV was circulating more than WNV in North-eastern Italy.

Molecular evidence of Plasmodium vivax infection in Duffy negative symptomatic individuals from Dschang, West Cameroon.

BACKGROUND: Plasmodium vivax infection is known to be rare in West/Central Africa, the most accepted explanation being the lack of expression of erythroid Duffy antigen in the local human populations. Duffy negativity prevents the parasite to exploit the entry mechanism on the red blood cell surface. However, there are a growing number of reported vivax infections in Duffy-negative individuals. Data on P. vivax circulation in Cameroon are limited. The aim of the study was to evaluate the P. vivax presence, and its association with the Duffy genotype in West Cameroon. RESULTS: Overall, 484 blood samples were collected consecutively from febrile outpatients attending the Dschang's Hospital (West Cameroon) during a 3-months period. Plasmodium vivax infection was detected by PCR in 5.6% (n = 27/484) of the cases, representing 38.6% (n = 27/70) of all Plasmodium infections detected. All P. vivax infected individuals showed a Duffy-negative genotype, and the frequency of Duffy-positive individuals in the whole tested population was 1.7%. CONCLUSIONS: The results of this study confirm the circulation of P. vivax in Cameroon, as well as that the lack of expression of Duffy-antigen does not confer full protection against vivax malaria acquisition.

Natural humoral immune response to ribosomal P0 protein in colorectal cancer patients.

BACKGROUND: Tumor associated antigens are useful in colorectal cancer (CRC) management. The ribosomal P proteins (P0, P1, P2) play an important role in protein synthesis and tumor formation. The immunogenicity of the ribosomal P0 protein in head and neck, in breast and prostate cancer patients and the overexpression of the carboxyl-terminal P0 epitope (C-22 P0) in some tumors were reported. METHODS: Sera from 72 colorectal tumor patients (67 malignant and 5 benign tumors) were compared with 73 healthy donor sera for the presence of antibodies to CEA, EGFR, ErbB2 and ribosomal P proteins by western blotting or ELISA. Expression of the C-22 P0 epitope on tissues and colon cancer cells was determined by immunoperoxidase staining and indirect immunofluorescence/western blotting, respectively, employing MAb 2B2. Biological effects of MAb 2B2 on colon cancer cells were assessed by the Sulforhodamine B cell proliferation assay, trypan blue exclusion test and cleaved caspase-3 detection. Fisher's exact test was used to compare the number of auto-antibodies positive patients with healthy donors. Variation in the C-22 P0 expression, and in the number of apoptotic cells was evaluated by Student's t-test. Variation in cell survival and cell death was evaluated by Newman-Keuls test. RESULTS: No significant humoral response was observed to CEA, EGFR and ErbB2 in CRC patients. Conversely, 7 out of 67 CRC patient sera reacted to ribosomal P proteins. The prevalence of P proteins auto-antibodies in CRC patients was significant. Five patients showed restricted P0 immunoreactivity, while two patients reacted simultaneously to all P proteins. The C-22 P0 epitope was homogenously expressed both in malignant tumors and the adjacent mucosa, but the intensity of expression was higher in the tumor. Starved colon cancer cells showed a higher C-22 P0 epitope plasma membrane expression compared to control cells. MAb 2B2 inhibited colon cancer cell growth and induced cell death in a dose dependent manner. CONCLUSIONS: Our study shows a spontaneous humoral immune response to ribosomal P0 protein in CRC patients and the inhibition of in vitro cancer cell growth after C-22 P0 epitope targeting. The ribosomal P0 protein might be a useful immunological target in CRC patients.

Co-circulation of Dengue and Chikungunya Viruses, Al Hudaydah, Yemen, 2012.

We investigated 400 cases of dengue-like illness in persons hospitalized during an outbreak in Al Hudaydah, Yemen, in 2012. Overall, 116 dengue and 49 chikungunya cases were diagnosed. Dengue virus type 2 was the predominant serotype. The co-circulation of these viruses indicates that mosquitoborne infections represent a public health threat in Yemen.

Virus investigation in ticks from migratory birds in Italy.

The role of migratory birds in circulation tick-borne viruses needs to be better defined. In order to assess the potential role of migratory birds in exotic virus spread, we conducted a study to identify ticks collected from migratory birds in the Central Region of Italy, and performed molecular investigation for Crimea-Congo Hemorrhagic Fever (CCHFV), West Nile fever (WNFV) and Usutu (USUV) in the vectors. A total of 137 competent ticks were collected with predominance of Hyalomma species. Although, negative results were obtained for all viruses considered, the high proportion of Hyalomma ticks highlights the potential risk for the dissemination of tick-borne viruses through infested migratory birds.

Diketo acid inhibitors of nsp13 of SARS-CoV-2 block viral replication.

For RNA viruses, RNA helicases have long been recognized to play critical roles during virus replication cycles, facilitating proper folding and replication of viral RNAs, therefore representing an ideal target for drug discovery. SARS-CoV-2 helicase, the non-structural protein 13 (nsp13) is a highly conserved protein among all known coronaviruses, and, at the moment, is one of the most explored viral targets to identify new possible antiviral agents. In the present study, we present six diketo acids (DKAs) as nsp13 inhibitors able to block both SARS-CoV-2 nsp13 enzymatic functions. Among them four compounds were able to inhibit viral replication in the low micromolar range, being active also on other human coronaviruses such as HCoV229E and MERS CoV. The experimental investigation of the binding mode revealed ATP-non-competitive kinetics of inhibition, not affected by substrate-displacement effect, suggesting an allosteric binding mode that was further supported by molecular modelling calculations predicting the binding into an allosteric conserved site located in the RecA2 domain.

Durability of neutralizing antibodies against yellow fever virus after vaccination in healthy adults.

In accordance with the World Health Organization, one dose of yellow fever vaccine may guarantee protection lifelong in healthy adults. However, relatively little information is still available from ad hoc studies. We evaluated the persistence of neutralizing antibodies, which are considered to be an immune correlate of protection, in a large number of military personnel vaccinated up to 47 years before. Overall, 322 individuals were studied. The median time from vaccination to blood collection for neutralizing antibody evaluation was 9 years, ranging from <1 to 47 years. Of the 322 participants, 319 had neutralizing antibodies (99.1 %). The highest median PRNT50 value was observed in those vaccinated ≤1 year before (median PRNT50 = 320). In conclusion, our study confirms on a larger scale that, in healthy adults, neutralizing antibodies may persist as long as 47 years after a single yellow fever vaccines dose.

Pan-yellow fever virus detection and lineage assignment by real-time RT-PCR and amplicon sequencing.

Yellow fever disease is a viral zoonosis that may result in a severe hemorrhagic disease. A safe and effective vaccine used in mass immunization campaigns has allowed control and mitigation against explosive outbreaks in endemic areas. Since the 1960's, re-emergent of the yellow fever virus has been observed. The timely implementation of control measures, to avoid or contain an ongoing outbreak requires rapid specific viral detection methods. Here a novel molecular assay, expected to detect all known yellow fever virus strains, is described. The method has demonstrated high sensitivity and specificity in real-time RT-PCR as well as in an endpoint RT-PCR set-up. Sequence alignment and phylogenetic analysis reveal that the amplicon resulting from the novel method covers a genomic region whose mutational profile is completely associated to the yellow fever viral lineages. Therefore, sequencing analysis of this amplicon allows for assignment of the viral lineage.

West Nile alternative open reading frame (N-NS4B/WARF4) is produced in infected West Nile Virus (WNV) cells and induces humoral response in WNV infected individuals.

BACKGROUND: West Nile Virus (WNV) is a flavivirus that requires an efficient humoral and cellular host response for the control of neuroinvasive infection. We previously reported the existence of six alternative open reading frame proteins in WNV genome, one of which entitled WARF4 is exclusively restricted to the lineage I of the virus. WARF4 is able to elicit antibodies in WNV infected horses; however, there was no direct experimental proof of the existence of this novel protein. The purpose of this study was to demonstrate the in vitro production of WARF4 protein following WNV infection of cultured VERO cells and its immunity in WNV infected individuals. RESULTS: We produced a monoclonal antibody against WARF4 protein (MAb 3A12) which detected the novel protein in WNV lineage I-infected, cultured VERO cells while it did not react with WNV lineage II infected cells. MAb 3A12 specificity to WARF4 protein was confirmed by its reactivity to only one peptide among four analyzed that cover the full WARF4 amino acids sequence. In addition, WARF4 protein was expressed in the late phase of WNV lineage I infection. Western blotting and bioinformatics analyses strongly suggest that the protein could be translated by programmed -1 ribosomal frameshifting process. Since WARF4 is embedded in the NS4B gene, we rename this novel protein N-NS4B/WARF4. Furthermore, serological analysis shows that N-NS4B/WARF4 is able to elicit antibodies in WNV infected individuals. CONCLUSIONS: N-NS4B/WARF4 is the second Alternative Reading Frame (ARF) protein that has been demonstrated to be produced following WNV infection and might represent a novel tool for a better characterization of immune response in WNV infected individuals. Further serological as well as functional studies are required to characterize the function of the N-NS4B/WARF4 protein. Since the virus might actually make an extensive use of ARFs, it appears important to investigate the novel six ARF putative proteins of WNV.

Real-World Clinical Outcomes of Ribociclib in Combination with a Non-Steroidal Aromatase Inhibitor and a Luteinizing Hormone-Releasing Hormone Agonist in Premenopausal HR+/HER2- Advanced Breast Cancer Patients: An Italian Managed Access Program.

Ribociclib plus an aromatase inhibitor and ovarian function suppression is the preferred first-line option for pre-/perimenopausal women with hormone receptor-positive/human epidermal growth factor receptor-2-negative advanced or metastatic breast cancer. We opened an italian managed access program (MAP) that permitted access to ribociclib to selected patients and allowed to collect informative results on the clinical impact of the therapy. The MAP (April 2018-May 2020) included 64 premenopausal patients, with characteristics similar to those of the MONALEESA-7 trial. Of 57 patients with a known response, 48 (84.2%) achieved a clinical benefit (i.e., complete response, N = 7 (12.3%); partial response, N = 17 (29.8%); stable disease, N = 24 (42.1%)), while 9 (15.8%) experienced tumor progression. Some patients (N = 15-23.4%) needed ribociclib dose reduction because of adverse events. Thereafter, the treatment was well tolerated, and no new safety signals emerged. Our study is the first reported Italian real-world evidence of ribociclib effectiveness in premenopausal HR+/HER2- advanced breast cancer patients. Response and clinical benefit rates were particularly encouraging compared with those of the ribociclib group of MONALEESA-7. Our work confirms that ribociclib in combination with endocrine therapy is highly effective in the treatment of premenopausal HR+/HER2- advanced breast cancer patients with an expected safety profile.

Analysis of Genomic Characteristics of SARS-CoV-2 in Italy, 29 January to 27 March 2020.

We performed next-generation sequencing (NGS), phylogenetic analysis, gene flows, and N- and O-glycosylation prediction on SARS-CoV-2 genomes collected from lab-confirmed cases from different Italian regions. To this end, a total of 111 SARS-CoV-2 genomes collected in Italy between 29 January and 27 March 2020 were investigated. The majority of the genomes belonged to lineage B.1, with some descendant lineages. The gene flow analysis showed that the spread occurred mainly from the north to the center and to the south of Italy, as confirmed by epidemiological data. The mean evolutionary rate estimated here was 8.731 × 10-4 (95% highest posterior density, HPD intervals 5.809 × 10-4 to 1.19 × 10-3), in line with values reported by other authors. The dated phylogeny suggested that SARS-CoV-2 lineage B.1 probably entered Italy between the end of January and early February 2020. Continuous molecular surveillance is needed to trace virus circulation and evolution.

High throughput MLVA-16 typing for Brucella based on the microfluidics technology.

BACKGROUND: Brucellosis, a zoonosis caused by the genus Brucella, has been eradicated in Northern Europe, Australia, the USA and Canada, but remains endemic in most areas of the world. The strain and biovar typing of Brucella field samples isolated in outbreaks is useful for tracing back source of infection and may be crucial for discriminating naturally occurring outbreaks versus bioterrorist events, being Brucella a potential biological warfare agent. In the last years MLVA-16 has been described for Brucella spp. genotyping. The MLVA band profiles may be resolved by different techniques i.e. the manual agarose gels, the capillary electrophoresis sequencing systems or the microfluidic Lab-on-Chip electrophoresis. In this paper we described a high throughput system of MLVA-16 typing for Brucella spp. by using of the microfluidics technology. RESULTS: The Caliper LabChip 90 equipment was evaluated for MLVA-16 typing of sixty-three Brucella samples. Furthermore, in order to validate the system, DNA samples previously resolved by sequencing system and Agilent technology, were de novo genotyped. The comparison of the MLVA typing data obtained by the Caliper equipment and those previously obtained by the other analysis methods showed a good correlation. However the outputs were not accurate as the Caliper DNA fragment sizes showed discrepancies compared with real data and a conversion table from observed to expected data was created. CONCLUSION: In this paper we described the MLVA-16 using a rapid, sophisticated microfluidics technology for detection of amplification product sizes. The comparison of the MLVA typing data produced by Caliper LabChip 90 system with the data obtained by different techniques showed a general concordance of the results. Furthermore this platform represents a significant improvement in terms of handling, data acquiring, computational efficiency and rapidity, allowing to perform the strain genotyping in a time equal to one sixth respect to other microfluidics systems as e.g. the Agilent 2100 bioanalyzer.Finally, this platform can be considered a valid alternative to standard genotyping techniques, particularly useful dealing with a large number of samples in short time. These data confirmed that this technology represents a significative advancement in high-throughput accurate Brucella genotyping.

Identification and characterization of SARS-CoV-2 clusters in the EU/EEA in the first pandemic wave: additional elements to trace the route of the virus.

A high-quality dataset of 3289 complete SARS-CoV-2 genomes collected in Europe and European Economic Area (EAA) in the early phase of the first wave of the pandemic was analyzed. Among all single nucleotide mutations, 41 had a frequency ≥ 1%, and the phylogenetic analysis showed at least 6 clusters with a specific mutational profile. These clusters were differentially distributed in the EU/EEA, showing a statistically significant association with the geographic origin. The analysis highlighted that the mutations C14408T and C14805T played an important role in clusters selection and further virus spread. Moreover, the molecular analysis suggests that the SARS-CoV-2 strain responsible for the first Italian confirmed COVID-19 case was already circulating outside the country.