Uncovering genetic causes of hypophosphatemia

Background Chronic hypophosphatemia can result from a variety of acquired disorders, such as malnutrition, intestinal malabsorption, hyperparathyroidism, vitamin D deficiency, excess alcohol intake, some drugs, or organ transplantation. Genetic disorders can be a cause of persistent hypophosphatemia, although they are less recognized. We aimed to better understand the prevalence of genetic hypophosphatemia in the population. Methods By combining retrospective and prospective strategies, we searched the laboratory database of 815,828 phosphorus analyses and included patients 17?55 years old with low serum phosphorus. We reviewed the charts of 1287 outpatients with at least 1 phosphorus result ?2.2 mg/dL. After ruling out clear secondary causes, 109 patients underwent further clinical and analytical studies. Among them, we confirmed hypophosphatemia in 39 patients. After excluding other evident secondary causes, such as primary hyperparathyroidism and vitamin D deficiency, we performed a molecular analysis in 42 patients by sequencing the exonic and flanking intronic regions of a panel of genes related to rickets or hypophosphatemia (CLCN5, CYP27B1, dentin matrix acidic phosphoprotein 1, ENPP1, FAM20C, FGFR1, FGF23, GNAS, PHEX, SLC34A3, and VDR). Results We identified 14 index patients with hypophosphatemia and variants in genes related to phosphate metabolism. The phenotype of most patients was mild, but two patients with X-linked hypophosphatemia (XLH) due to novel PHEX mutations had marked skeletal abnormalities. Conclusion Genetic causes should be considered in children, but also in adult patients with hypophosphatemia of unknown origin. Our data are consistent with the conception that XLH is the most common cause of genetic hypophosphatemia with an overt musculoskeletal phenotype

Guía para la aplicación de la Ley 11/2018, Información no financiera: versión 1.0

Ministerio de Economía y Competitividad, FEDER y la Consejería de Educación, Junta de Castilla y León (ECO2015- 65782-P y BU058P17). Asimismo, agradecen la financiación recibida de la Junta de Castilla y León para la Contratación de personal técnico de apoyo a la investigación por las universidades públicas de Castilla y León (Orden de 14 de diciembre de 2017, de la Consejería de Educación, cofinanciadas por el Fondo Social Europeo y la Iniciativa de Empleo Juvenil

Predictors of Loss of Functional Independence in Parkinson’s Disease: Results from the COPPADIS Cohort at 2-Year Follow-Up and Comparison with a Control Group

COPPADIS Study Group.[Background and objective] The aim of this study was to compare the progression of independence in activities of daily living (ADL) in Parkinson’s disease (PD) patients versus a control group, as well as to identify predictors of disability progression and functional dependency (FD).[Patients and Methods] PD patients and control subjects, who were recruited from 35 centers of Spain from the COPPADIS cohort between January 2016 and November 2017 (V0), were included. Patients and subjects were then evaluated again at the 2-year follow-up (V2). Disability was assessed with the Schwab & England Activities of Daily Living Scale (S&E-ADLS) at V0 and V2. FD was defined as an S&E-ADLS score less than 80%.[Results] In the PD group, a significant decrease in the S&E-ADLS score from V0 to V2 (N = 507; from 88.58 ± 10.19 to 84.26 ± 13.38; p < 0.0001; Cohen’s effect size = −0.519) was observed but not in controls (N = 124; from 98.87 ± 6.52 to 99.52 ± 2.15; p = 0.238). When only patients considered functional independent at baseline were included, 55 out of 463 (11.9%) converted to functional dependent at V2. To be a female (OR = 2.908; p = 0.009), have longer disease duration (OR = 1.152; p = 0.002), have a non-tremoric motor phenotype at baseline (OR = 3.574; p = 0.004), have a higher score at baseline in FOGQ (OR = 1.244; p < 0.0001) and BDI-II (OR = 1.080; p = 0.008), have a lower score at baseline in PD-CRS (OR = 0.963; p = 0.008), and have a greater increase in the score from V0 to V2 in UPDRS-IV (OR = 1.168; p = 0.0.29), FOGQ (OR = 1.348; p < 0.0001) and VAFS-Mental (OR = 1.177; p = 0.013) (adjusted R-squared 0.52; Hosmer and Lemeshow test = 0.94) were all found to be independent predictors of FD at V2.[Conclusions] In conclusion, autonomy for ADL worsens in PD patients compared to controls. Cognitive impairment, gait problems, fatigue, depressive symptoms, more advanced disease, and a non-tremor phenotype are independent predictors of FD in the short-term.Fundación Curemos el Parkinson (www.curemoselparkinson.org).Peer reviewe

Staging Parkinson’s Disease Combining Motor and Nonmotor Symptoms Correlates with Disability and Quality of Life

COPPADIS Study Group.[Introduction] In a degenerative disorder such as Parkinson’s disease (PD), it is important to establish clinical stages that allow to know the course of the disease. Our aim was to analyze whether a scale combining Hoehn and Yahr’s motor stage (H&Y) and the nonmotor symptoms burden (NMSB) (assessed by the nonmotor symptoms scale (NMSS)) provides information about the disability and the patient’s quality of life (QoL) with regard to a defined clinical stage.[Materials and Methods] Cross-sectional study in which 603 PD patients from the COPPADIS cohort were classified according to H&Y (1, stage I; 2, stage II; 3, stage III; 4, stage IV/V) and NMSB (A: NMSS = 0–20; B: NMSS = 21–40; C: NMSS = 41–70; D: NMSS ≥ 71) in 16 stages (HY.NMSB, from 1A to 4D). QoL was assessed with the PDQ-39SI, PQ-10, and EUROHIS-QOL8 and disability with the Schwab&England ADL (Activities of Daily Living) scale.[Results] A worse QoL and greater disability were observed at a higher stage of H&Y and NMSB (). Combining both (HY.NMSB), patients in stages 1C and 1D and 2C and 2D had significantly worse QoL and/or less autonomy for ADL than those in stages 2A and 2B and 3A and 3B, respectively (; e.g., PDQ-39SI in 1D [n = 15] vs 2A [n = 101]: 28.6 ± 17.1 vs 7.9 ± 5.8; ).[Conclusion] The HY.NMSB scale is simple and reflects the degree of patient involvement more accurately than the H&Y. Patients with a lower H&Y stage may be more affected if they have a greater NMS burden.Peer reviewe

Predictors of Global Non-Motor Symptoms Burden Progression in Parkinson’s Disease. Results from the COPPADIS Cohort at 2-Year Follow-Up

COPPADIS Study Group.[Background and Objective] Non-motor symptoms (NMS) progress in different ways between Parkinson’s disease (PD) patients. The aim of the present study was to (1) analyze the change in global NMS burden in a PD cohort after a 2-year follow-up, (2) to compare the changes with a control group, and (3) to identify predictors of global NMS burden progression in the PD group.[Material and Methods] PD patients and controls, recruited from 35 centers of Spain from the COPPADIS cohort from January 2016 to November 2017, were followed-up with after 2 years. The Non-Motor Symptoms Scale (NMSS) was administered at baseline (V0) and at 24 months ± 1 month (V2). Linear regression models were used for determining predictive factors of global NMS burden progression (NMSS total score change from V0 to V2 as dependent variable).[Results] After the 2-year follow-up, the mean NMS burden (NMSS total score) significantly increased in PD patients by 18.8% (from 45.08 ± 37.62 to 53.55 ± 42.28; p < 0.0001; N = 501; 60.2% males, mean age 62.59 ± 8.91) compared to no change observed in controls (from 14.74 ± 18.72 to 14.65 ± 21.82; p = 0.428; N = 122; 49.5% males, mean age 60.99 ± 8.32) (p < 0.0001). NMSS total score at baseline (β = −0.52), change from V0 to V2 in PDSS (Parkinson’s Disease Sleep Scale) (β = −0.34), and change from V0 to V2 in NPI (Neuropsychiatric Inventory) (β = 0.25) provided the highest contributions to the model (adjusted R-squared 0.41; Durbin-Watson test = 1.865).[Conclusions] Global NMS burden demonstrates short-term progression in PD patients but not in controls and identifies worsening sleep problems and neuropsychiatric symptoms as significant independent predictors of this NMS progression.This research was funded by Fundación Española de Ayuda a la Investigación en Parkinson y otras Enfermedades Neuro-degenerativas (Curemos el Parkinson; www.curemoselparkinson.org).Peer reviewe

Predictors of clinically significant quality of life impairment in Parkinson’s disease

COPPADIS Study Group.Quality of life (QOL) plays an important role in independent living in Parkinson’s disease (PD) patients, being crucial to know what factors impact QoL throughout the course of the disease. Here we identified predictors of QoL impairment in PD patients from a Spanish cohort. PD patients recruited from 35 centers of Spain from the COPPADIS cohort from January 2016, to November 2017, were followed up during 2 years. Health-related QoL (HRQoL) and global QoL (GQoL) were assessed with the 39-item Parkinson’s disease Questionnaire (PDQ-39) and the EUROHIS-QOL 8-item index (EUROHIS-QOL8), respectively, at baseline (V0) and at 24 months ± 1 month (V2). Clinically significant QoL impairment was defined as presenting an increase (PDQ-39SI) or decrement (EUROHIS-QOL8) at V2 ≥ 10% of the score at baseline (V0). A comparison with a control group was conducted for GQoL. GQoL did not change significantly in PD patients (N = 507; p = 0.686) or in the control group (N = 119; p = 0.192). The mean PDQ-39SI was significantly increased in PD patients (62.7 ± 8.5 years old; 58.8% males; N = 500) by 21.6% (from 16.7 ± 13 to 20.3 ± 16.4; p < 0.0001) at V2. Ninety-three patients (18.6%) presented a clinically significant HRQoL impairment at V2. To be younger (OR = 0.896; 95% CI 0.829–0.968; p = 0.006), to be a female (OR = 4.181; 95% CI 1.422–12.290; p = 0.009), and to have a greater increase in BDI-II (Beck Depression Inventory-II) (OR = 1.139; 95% CI 1.053–1.231; p = 0.001) and NMSS (Non-Motor Symptoms Scale) (OR = 1.052; 95% CI 1.027–1.113; p < 0.0001) total scores from V0 to V2 were associated with clinically significant HRQoL impairment at the 2-year follow-up (Hosmer–Lemeshow test, p = 0.665; R 2 = 0.655). An increase in ≥5 and ≥10 points of BDI-II and NMSS total score at V2 multiplied the probability of presenting clinically significant HRQoL impairment by 5 (OR = 5.453; 95% CI 1.663–17.876; p = 0.005) and 8 (OR = 8.217; 95% CI, 2.975–22.696; p = 0.002), respectively. In conclusion, age, gender, mood, and non-motor impairment were associated with clinically significant HRQoL impairment after the 2-year follow-up in PD patients.Mir P. has received honoraria from AbbVie, Abbott, Allergan, Bial, Merz, UCB and Zambon and have received grants from the Spanish Ministry of Economy and Competitiveness [PI16/01575] co-founded by ISCIII (Subdirección General de Evaluación y Fomento de la Investigación) and by Fondo Europeo de Desarrollo Regional (FEDER), the Consejería de Economía, Innovación, Ciencia y Empleo de la Junta de Andalucía [CVI-02526, CTS-7685], the Consejería de Salud y Bienestar Social de la Junta de Andalucía [PI-0437-2012, PI-0471-2013], the Sociedad Andaluza de Neurología, the Jacques and Gloria Gossweiler Foundation, the Fundación Alicia Koplowitz, the Fundación Mutua Madrileña.Peer reviewe

Singular Location and Signaling Profile of Adenosine A2A-Cannabinoid CB1 Receptor Heteromers in the Dorsal Striatum

The dorsal striatum is a key node for many neurobiological processes such as motor activity, cognitive functions, and affective processes. The proper functioning of striatal neurons relies critically on metabotropic receptors. Specifically, the main adenosine and endocannabinoid receptors present in the striatum, ie, adenosine A2A receptor (A2AR) and cannabinoid CB1 receptor (CB1R), are of pivotal importance in the control of neuronal excitability. Facilitatory and inhibitory functional interactions between striatal A2AR and CB1R have been reported, and evidence supports that this cross-talk may rely, at least in part, on the formation of A2AR-CB1R heteromeric complexes. However, the specific location and properties of these heteromers have remained largely unknown. Here, by using techniques that allowed a precise visualization of the heteromers in situ in combination with sophisticated genetically-modified animal models, together with biochemical and pharmacological approaches, we provide a high resolution expression map and a detailed functional characterization of A2AR-CB1R heteromers in the dorsal striatum. Specifically, our data unveil that the A2AR-CB1R heteromer (i) is essentially absent from corticostriatal projections and striatonigral neurons, and, instead, is largely present in striatopallidal neurons, (ii) displays a striking G protein-coupled signaling profile, where co-stimulation of both receptors leads to strongly reduced downstream signaling, and (iii) undergoes an unprecedented dysfunction in Huntington’s disease, an archetypal disease that affects striatal neurons. Altogether, our findings may open a new conceptual framework to understand the role of coordinated adenosine-endocannabinoid signaling in the indirect striatal pathway, which may be relevant in motor function and neurodegenerative diseases

Diagnóstico Territorial del Litoral de Cantabria. Volumen I: Informe.

Estudios de base para la redacción del Plan de Ordenación del Litoral (POL) de Cantabria.Este proyecto de investigación aplicada se ha realizado gracias al Convenio de Colaboración entre el Gobierno de Cantabria y la Universidad de Cantabria titulado “Diagnóstico Territorial del Litoral de Cantabria”

Falls Predict Acute Hospitalization in Parkinson's Disease

[Background] There is a need for identifying risk factors for hospitalization in Parkinson’s disease (PD) and also interventions to reduce acute hospital admission.[Objective] To analyze the frequency, causes, and predictors of acute hospitalization (AH) in PD patients from a Spanish cohort.[Methods] PD patients recruited from 35 centers of Spain from the COPPADIS-2015 (COhort of Patients with PArkinson’s DIsease in Spain, 2015) cohort from January 2016 to November 2017, were included in the study. In order to identify predictors of AH, Kaplan-Meier estimates of factors considered as potential predictors were obtained and Cox regression performed on time to hospital encounter 1-year after the baseline visit.[Results] Thirty-five out of 605 (5.8%) PD patients (62.5±8.9 years old; 59.8% males) presented an AH during the 1-year follow-up after the baseline visit. Traumatic falls represented the most frequent cause of admission, being 23.7% of all acute hospitalizations. To suffer from motor fluctuations (HR [hazard ratio] 2.461; 95% CI, 1.065–5.678; p = 0.035), a very severe non-motor symptoms burden (HR [hazard ratio] 2.828; 95% CI, 1.319–6.063; p = 0.008), falls (HR 3.966; 95% CI 1.757–8.470; p = 0.001), and dysphagia (HR 2.356; 95% CI 1.124–4.941; p = 0.023) was associated with AH after adjustment to age, gender, disease duration, levodopa equivalent daily dose, total number of non-antiparkinsonian drugs, and UPDRS-IIIOFF. Of the previous variables, only falls (HR 2.998; 95% CI 1.080–8.322; p = 0.035) was an independent predictor of AH.[Conclusion] Falls is an independent predictor of AH in PD patients.Peer reviewe

The first wave of the COVID-19 epidemic in Spain was associated with early introductions and fast spread of a dominating genetic variant

SeqCOVID-Spain consortium: Álvaro Chiner-Oms, Irving Cancino-Muñoz, Mariana G. López, Manuela Torres-Puente, Inmaculada Gómez-Navarro, Santiago Jiménez-Serrano, Jordi Pérez-Tur, Darío García de Viedma, Laura Pérez-Lago, Marta Herranz, Jon Sicilia, Pilar Catalán-Alonso, Julia Suárez González, Patricia Muñoz, Mireia Coscolla, Paula Ruiz-Rodríguez, Fernando González-Candelas, Iñaki Comas, Lidia Ruiz-Roldán, María Alma Bracho, Neris García-González, Llúcia Martínez Priego, Inmaculada Galán-Vendrell, Paula Ruiz-Hueso, Griselda De Marco, María Loreto Ferrús-Abad, Sandra Carbó-Ramírez, Giuseppe D’Auria, Galo Adrian Goig, Juan Alberola, Jose Miguel Nogueira, Juan José Camarena, David Navarro, Eliseo Albert, Ignacio Torres, Maitane Aranzamendi Zaldumbide, Óscar Martínez Expósito, Nerea Antona Urieta, María de Toro, María Pilar Bea-Escudero, Jose Antonio Boga, Cristian Castelló-Abietar, Susana Rojo-Alba, Marta Elena Álvarez-Argüelles, Santiago Melón, Elisa Martró, Antoni E. Bordoy, Anna Not, Adrián Antuori, Anabel Fernández-Navarro, Andrés Canut-Blasco, Silvia Hernáez Crespo, Maria Luz Cordón Rodríguez, Maria Concepción Lecaroz Agara, Carmen Gómez-González, Amaia Aguirre-Quiñonero, José Israel López-Mirones, Marina Fernández-Torres, Maria Rosario Almela-Ferrer, Ana Carvajal, Juan Miguel Fregeneda-Grandes, Héctor Argüello, Gustavo Cilla Eguiluz, Milagrosa Montes Ros, Luis Piñeiro Vázquez, Ane Sorarrain, José María Marimón, José J. Costa-Alcalde, Rocío Trastoy, Gema Barbeito Castiñeiras, Amparo Coira, María Luisa Pérez del Molino, Antonio Aguilera, Begoña Palop-Borrás, Inmaculada de Toro Peinado, Maria Concepción Mediavilla Gradolph, Mercedes Pérez-Ruiz, Mirian Fernández-Alonso, Jose Luis del Pozo, Oscar González-Recio, Mónica Gutiérrez-Rivas, Jovita Fernández-Pinero, Miguel Ángel Jiménez Clavero, Begoña Fuster Escrivá, Concepción Gimeno Cardona, María Dolores Ocete Mochón, Rafael Medina-Gonzalez, José Antonio Lepe, Verónica González Galán, Ángel Rodríguez-Villodres, Nieves Gonzalo Jiménez, Jordi Reina, Carla López-Causapé, Maria Dolores Gómez-Ruiz, Eva M. Gonzalez-Barbera, José Luis López-Hontangas, Vicente Martín, Antonio J. Molina, Tania Fernandez-Villa, Ana Milagro Beamonte, Nieves Felisa Martínez-Cameo, Yolanda Gracia-Grataloup, Rosario Moreno-Muñoz, Maria Dolores Tirado Balaguer, José María Navarro-Marí, Irene Pedrosa-Corral, Sara Sanbonmatsu-Gámez, Antonio Oliver, Mónica Parra Grande, Bárbara Gómez Alonso, Francisco José Arjona Zaragozí, Maria Carmen Pérez González, Francisco Javier Chamizo López, Ana Bordes-Benítez, Núria Rabella, Ferran Navarro, Elisenda Miró, Antonio Rezusta, Alexander Tristancho, Encarnación Simarro Córdoba, Julia Lozano-Serra, Lorena Robles Fonseca, Álex Soriano, Francisco Javier Roig Sena, Hermelinda Vanaclocha Luna, Isabel Sanmartín, Daniel García-Souto, Ana Pequeño-Valtierra, Jose M. C. Tubio, Javier Temes, Jorge Rodríguez-Castro, Martín Santamarina García, Manuel Rodríguez-Iglesias, Fátima Galán-Sanchez, Salud Rodríguez-Pallares, José Manuel Azcona-Gutiérrez, Miriam Blasco-Alberdi, Alfredo Mayor, Alberto L. García-Basteiro, Gemma Moncunill, Carlota Dobaño, Pau Cisteró, Oriol Mitjà, Camila González-Beiras, Martí Vall-Mayans, Marc Corbacho-Monné, Andrea Alemany, Cristina Muñoz-Cuevas, Guadalupe Rodríguez-Rodríguez, Rafael Benito, Sonia Algarate, Jessica Bueno, Andrea Vergara-Gómez, Miguel J. Martínez, Jordi Vila, Elisa Rubio, Aida Peiró-Mestres, Jessica Navero-Castillejos, David Posada, Diana Valverde, Nuria Estévez, Iria Fernández-Silva, Loretta de Chiara, Pilar Gallego-García, Nair Varela, Ulises Gómez-Pinedo, Mónica Gozalo-Margüello, Maria Eliecer Cano García, José Manuel Méndez-Legaza, Jesus Rodríguez-Lozano, María Siller, Daniel Pablo-Marcos, Maria Montserrat Ruiz-García, Antonio Galiana, Judith Sánchez-Almendro, Maria Isabel Gascón Ros, Cristina Juana Torregrosa-Hetland, Eva María Pastor Boix, Paloma Cascales Ramos, Pedro Luis Garcinuño Enríquez, Salvador Raga Borja, Julia González Cantó, Olalla Martínez Macias, Adolfo de Salazar, Laura Viñuela González, Natalia Chueca, Federico García, Cristina Gómez-Camarasa, Amparo Farga Martí, Rocío Falcón, Victoria Domínguez-Márquez, Anna M. Planas, Israel Fernández-Cádenas, Maria Ángeles Marcos, Carmen Ezpeleta, Ana Navascués, Ana Miqueleiz Zapatero, Manuel Segovia, Antonio Moreno-Docón, Esther Viedma, Raúl Recio Martínez, Irene Muñoz-Gallego, Sara Gonzalez-Bodi, Maria Dolores Folgueira, Jesús Mingorance, Elias Dahdouh, Fernando Lázaro-Perona, María Rodríguez-Tejedor, María Pilar Romero-Gómez, Julio García-Rodríguez, Juan Carlos Galán, Mario Rodríguez-Dominguez, Laura Martínez-García, Melanie Abreu Di Berardino, Manuel Ponce-Alonso, Jose Maria González-Alba, Ivan Sanz-Muñoz, Diana Pérez San José, Maria Gil Fortuño, Juan B. Bellido-Blasco, Alberto Yagüe Muñoz, Noelia Hernández Pérez, Helena Buj Jordá, Óscar Pérez Olaso, Alejandro González Praetorius, Nora Mariela Martínez Ramírez, Aida Ramírez Marinero, Eduardo Padilla León, Alba Vilas Basil, Mireia Canal Aranda, Albert Bernet Sánchez, Alba Bellés Bellés, Eric López González, Iván Prats Sánchez, Mercè García-González, Miguel José Martínez-Lirola, Manuel Ángel Rodríguez Maresca, Maria Teresa Cabezas Fernández, María Eugenia Carrillo Gil, Maria Paz Ventero Martín, Carmen Molina Pardines, Nieves Orta Mira, María Navarro Cots, Inmaculada Vidal Catalá, Isabel García Nava, Soledad Illescas Fernández-Bermejo, José Martínez-Alarcón, Marta Torres-Narbona, Cristina Colmenarejo, Lidia García-Agudo, Jorge A. Pérez García, Martín Yago López, María Ángeles Goberna Bravo, Victoria Simón García, Gonzalo Llop Furquet, Agustín Iranzo Tatay, Sandra Moreno-Marro, Noelia Lozano Rodríguez, Amparo Broseta Tamarit, Juan José Badiola Díez, Amparo Martínez-Ramírez, Ana Dopazo, Sergio Callejas, Alberto Benguría, Begoña Aguado, Antonio Alcamí, Marta Bermejo Bermejo, Ricardo Ramos-Ruíz, Víctor Manuel Fernández Soria, Fernando Simón Soria & Mercedes Roig CardellsThe coronavirus disease 2019 (COVID-19) pandemic has affected the world radically since 2020. Spain was one of the European countries with the highest incidence during the first wave. As a part of a consortium to monitor and study the evolution of the epidemic, we sequenced 2,170 samples, diagnosed mostly before lockdown measures. Here, we identified at least 500 introductions from multiple international sources and documented the early rise of two dominant Spanish epidemic clades (SECs), probably amplified by superspreading events. Both SECs were related closely to the initial Asian variants of SARS-CoV-2 and spread widely across Spain. We inferred a substantial reduction in the effective reproductive number of both SECs due to public-health interventions (Re < 1), also reflected in the replacement of SECs by a new variant over the summer of 2020. In summary, we reveal a notable difference in the initial genetic makeup of SARS-CoV-2 in Spain compared with other European countries and show evidence to support the effectiveness of lockdown measures in controlling virus spread, even for the most successful genetic variants.This work was mainly funded by the Instituto de Salud Carlos III project COV20/00140, with additional funding by Spanish National Research Council project CSIC-COV19-021, Ministerio de Ciencia project PID2019-104477RB-100, ERC StG 638553 and ERC CoG 101001038 to I.C., and BFU2017-89594R to F.G.C. M.C. is supported by Ramón y Cajal program from Ministerio de Ciencia and grants RTI2018-094399-A-I00 and Generalitat Valenciana (Regional Government) project SEJI/2019/011. We gratefully acknowledge Hospital Universitari Vall d’Hebron, Instituto de Salud Carlos III, IrsiCaixa AIDS Research Lab and all the international researchers and institutions that submitted sequenced SARS-CoV-2 genomes to the GISAID’s EpiCov Database (Supplementary Table 1), as an important part of our analyses has been made possible by the sharing of their work. We also thank Unidad de Bioinformática y Estadística, Centro de Investigación Príncipe Felipe, for allowing us to use the Computer Cluster to perform some of the bioinformatic analysis.Peer reviewe