Phylogeography of SARS-CoV-2 pandemic in Spain: a story of multiple introductions, micro-geographic stratification, founder effects, and super-spreaders

Spain has been one of the main global pandemic epicenters for coronavirus disease 2019 (COVID-19). Here, we analyzed >41 000 genomes (including >26 000 high-quality (HQ) genomes) downloaded from the GISAID repository, including 1 245 (922 HQ) sampled in Spain. The aim of this study was to investigate genome variation of novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and reconstruct phylogeographic and transmission patterns in Spain. Phylogeographic analysis suggested at least 34 independent introductions of SARS-CoV-2 to Spain at the beginning of the outbreak. Six lineages spread very successfully in the country, probably favored by super-spreaders, namely, A2a4 (7.8%), A2a5 (38.4%), A2a10 (2.8%), B3a (30.1%), and B9 (8.7%), which accounted for 87.9% of all genomes in the Spanish database. One distinct feature of the Spanish SARS-CoV-2 genomes was the higher frequency of B lineages (39.3%, mainly B3a+B9) than found in any other European country. While B3a, B9, (and an important sub-lineage of A2a5, namely, A2a5c) most likely originated in Spain, the other three haplogroups were imported from other European locations. The B3a strain may have originated in the Basque Country from a B3 ancestor of uncertain geographic origin, whereas B9 likely emerged in Madrid. The time of the most recent common ancestor (TMRCA) of SARS-CoV-2 suggested that the first coronavirus entered the country around 11 February 2020, as estimated from the TMRCA of B3a, the first lineage detected in the country. Moreover, earlier claims that the D614G mutation is associated to higher transmissibility is not consistent with the very high prevalence of COVID-19 in Spain when compared to other countries with lower disease incidence but much higher frequency of this mutation (56.4% in Spain vs. 82.4% in rest of Europe). Instead, the data support a major role of genetic drift in modeling the micro-geographic stratification of virus strains across the country as well as the role of SARS-CoV-2 super-spreaders

Increased Serum Levels of sCD14 and sCD163 Indicate a Preponderant Role for Monocytes in COVID-19 Immunopathology

Background: Emerging evidence indicates a potential role for monocytes in COVID-19 immunopathology. We investigated two soluble markers of monocyte activation, sCD14 and sCD163, in COVID-19 patients, with the aim of characterizing their potential role in monocyte-macrophage disease immunopathology. To the best of our knowledge, this is the first study of its kind. Methods: Fifty-nine SARS-Cov-2 positive hospitalized patients, classified according to ICU or non-ICU admission requirement, were prospectively recruited and analyzed by ELISA for levels of sCD14 and sCD163, along with other laboratory parameters, and compared to a healthy control group. Results: sCD14 and sCD163 levels were significantly higher among COVID-19 patients, independently of ICU admission requirement, compared to the control group. We found a significant correlation between sCD14 levels and other inflammatory markers, particularly Interleukin-6, in the non-ICU patients group. sCD163 showed a moderate positive correlation with the time lapsed from admission to sampling, independently of severity group. Treatment with corticoids showed an interference with sCD14 levels, whereas hydroxychloroquine and tocilizumab did not. Conclusions: Monocyte-macrophage activation markers are increased and correlate with other inflammatory markers in SARS-Cov-2 infection, in association to hospital admission. These data suggest a preponderant role for monocyte-macrophage activation in the development of immunopathology of COVID-19 patients

Biogeographical origin and timing of the founder ichthyosis TGM1 c.1187G > A mutation in an isolated Ecuadorian population

An unusually high frequency of the lamellar ichthyosis TGM1 mutation, c.1187G > A, has been observed in the Ecuadorian province of Manabi. Recently, the same mutation has been detected in a Galician patient (Northwest of Spain). By analyzing patterns of genetic variation around this mutation in Ecuadorian patients and population matched controls, we were able to estimate the age of c.1187G > A and the time to their most recent common ancestor (TMRCA) of c.1187G > A Ecuadorian carriers. While the estimated mutation age is 41 generations ago (~1,025 years ago [ya]), the TMRCA of Ecuadorian c.1187G > A carrier haplotypes dates to just 17 generations (~425 ya). Probabilistic-based inferences of local ancestry allowed us to infer a most likely European origin of a few (16% to 30%) Ecuadorian haplotypes carrying this mutation. In addition, inferences on demographic historical changes based on c.1187G > A Ecuadorian carrier haplotypes estimated an exponential population growth starting ~20 generations, compatible with a recent founder effect occurring in Manabi. Two main hypotheses can be considered for the origin of c.1187G > A: (i) the mutation could have arisen in Spain >1,000 ya (being Galicia the possible homeland) and then carried to Ecuador by Spaniards in colonial times ~400 ya, and (ii) two independent mutational events originated this mutation in Ecuador and Galicia. The geographic and cultural characteristics of Manabi could have favored a founder effect that explains the high prevalence of TGM1 c.1187G > A in this region