Dispersal history of SARS-CoV-2 in Galicia, Spain

The dynamics of SARS-CoV-2 transmission are influenced by a variety of factors, including social restrictions and the emergence of distinct variants. In this study, we delve into the origins and dissemination of the Alpha, Delta, and Omicron variants of concern in Galicia, northwest Spain. For this, we leveraged genomic data collected by the EPICOVIGAL Consortium and from the GISAID database, along with mobility information from other Spanish regions and foreign countries. Our analysis indicates that initial introductions during the Alpha phase were predominantly from other Spanish regions and France. However, as the pandemic progressed, introductions from Portugal and the USA became increasingly significant. Notably, Galicia's major coastal cities emerged as critical hubs for viral transmission, highlighting their role in sustaining and spreading the virus. This research emphasizes the critical role of regional connectivity in the spread of SARS-CoV-2 and offers essential insights for enhancing public health strategies and surveillance measures.

Influence of perinatal and childhood exposure to tobacco and mercury in children's gut microbiota

Background

Early life determinants of the development of gut microbiome composition in infants have been widely investigated; however, if early life pollutant exposures, such as tobacco or mercury, have a persistent influence on the gut microbial community, its stabilization at later childhood remains largely unknown.

Objective

In this exposome-wide study, we aimed at identifying the contribution of exposure to tobacco and mercury from the prenatal period to childhood, to individual differences in the fecal microbiome composition of 7-year-old children, considering co-exposure to a width of established lifestyle and clinical determinants.

Methods

Gut microbiome was studied by 16S rRNA amplicon sequencing in 151 children at the genus level. Exposure to tobacco was quantified during pregnancy through questionnaire (active tobacco consumption, second-hand smoking -SHS) and biomonitoring (urinary cotinine) at 4 years (urinary cotinine, SHS) and 7 years (SHS). Exposure to mercury was quantified during pregnancy (cord blood) and at 4 years (hair). Forty nine other potential environmental determinants (12 at pregnancy/birth/infancy, 15 at 4 years and 22 at 7 years, such as diet, demographics, quality of living/social environment, and clinical records) were registered. We used multiple models to determine microbiome associations with pollutants including multi-determinant multivariate analysis of variance and linear correlations (wUnifrac, Bray-Curtis and Aitchison ß-diversity distances), single-pollutant permutational multivariate analysis of variance adjusting for co-variates (Aitchison), and multivariable association model with single taxa (MaAsLin2; genus). Sensitivity analysis was performed including genetic data in a subset of 107 children.

Results

Active smoking in pregnancy was systematically associated with microbiome composition and ß-diversity (R2 2-4%, p < 0.05, Aitchison), independently of other co-determinants. However, in the adjusted single pollutant models (PERMANOVA), we did not find any significant association. An increased relative abundance of Dorea and decreased relative abundance of Akkermansia were associated with smoking during pregnancy (q < 0.05).

Discussion

Our findings suggest a long-term sustainable effect of prenatal tobacco exposure on the children's gut microbiota. This effect was not found for mercury exposure or tobacco exposure during childhood. Assessing the role of these exposures on the children's microbiota, considering multiple environmental factors, should be further investigated.

Anal Dysplasia Screening in People Living with HIV: Long-Term Follow-Up in a Large Cohort from Northwest Spain

Anal squamous cell carcinoma (SCC) is not a common disease in the general population, although its incidence is higher in people living with human immunodeficiency virus (PLWH). Anal SCC is caused by human papillomavirus (HPV) infection and arises from premalignant lesions termed squamous intraepithelial lesions (SILs). SIL surveillance programs are based on the early detection and treatment of SILs, especially those with a higher risk of transforming into cancer. An anal surveillance program has been under development in our institution since 2011. In this context, we performed a retrospective cohort study at the anal dysplasia unit of Álvaro-Cunqueiro Hospital (Spain). Epidemiological and clinical data were gathered from our Infectious Diseases Sample Collection (an open sample cohort including PLWH) from January 2011 to January 2022. A total of 493 PLWH were considered, 122 (24.7%) of whom were diagnosed with anal dysplasia at baseline, including 2 cases of anal SCC. Briefly, most of individuals were young men (median age, 38 years old) born in Spain (76%), whose vaccination rate before their inclusion in the program was scarce (<3%). Throughout the study period, 81 (16.4%) cases were diagnosed with high-grade squamous-intraepithelial lesions (HSILs) and 3 with anal SCC. At the baseline, severe immunosuppression (i.e., nadir CD4+ lymphocyte count below 200 cell/μL), and prior diagnosis of condyloma acuminata were more frequent within the group with SILs. Conversely, the baseline CD4+ lymphocyte count was similar among both groups. HPV-16 was related to a higher risk of HSILs (odds ratio: 2.76). At the end of the follow-up, 385 PLWH had been retained in care; one patient had died of anal cancer. Anal dysplasia was common (25% of cases), especially among patients infected by HPV-16, diagnosed with condyloma acuminata, and who were severely immunosuppressed. HPV-16 was the main risk factor for the presentation of HSILs.

Case report: BA.1 subvariant showing a BA.2-like pattern using a variant-specific PCR assay due to a single point mutation downstream the spike 69/70 deletion

Background:

SARS-CoV-2 variant tracking is key to the genomic surveillance of the COVID-19 pandemic. While next-generation sequencing (NGS) is commonly used for variant determination, it is expensive and time-consuming. Variant-specific PCR (vsPCR) is a faster, cheaper method that detects specific mutations that are considered variant-defining. These tests usually rely on specific amplification when a mutation is present or a specific melting temperature peak after amplification.

Case presentation:

A discrepant result between vsPCR and NGS was found in seventeen SARS-CoV-2 samples from Galicia, Spain. A cluster of BA.1 Omicron SARS-CoV-2 variant showed a BA.2-like melting temperature pattern due to a point mutation (C21772T) downstream the deletion of the spike amino acids 69/70. As the 69/70 deletion is widely used for differentiation between BA.1 and BA.2 by vsPCR, C21772T can cause BA.1 samples to be misinterpreted as BA.2. Over a thousand BA.1 sequences in the EpiCoV database contain this mutation.

Conclusions:

To our knowledge, this is the first case of a point mutation causing a vsPCR algorithm to misclassify BA.1 samples as BA.2. This is an example of how mutations in the probe target area of vsPCR tests based on melting curve analysis can lead to variant misclassification. NGS confirmation of vsPCR results is relevant for the accuracy of the epidemiological surveillance. In order to overcome the possible impact of novel mutations, diagnostic tools must be constantly updated.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12985-022-01883-2.

Gastrointestinal colonization by OXA-48-producing Enterobacterales: risk factors for persistent carriage

Carbapenem-resistant Enterobacterales (CRE) infections are a major health problem. Intestinal colonization is a key factor in developing infection. However, factors associated with persistent colonization by CRE are unknown. The aim of the study was to identify factors associated with persistent CRE gut colonization. This is a retrospective, single-centre, observational study of adult patients with CRE gut colonization between January 2015 and January 2020. Epidemiologic characteristics, comorbidities, infectious events, duration of hospitalization and antimicrobial treatment received in the follow-up period were collected. Colonization was defined as isolation in at least 2 rectal swab culture samples of CRE. Decolonization was defined as 3 negative rectal swab cultures or 2 negative cultures and a negative molecular test. A cohort of 86 patients with CRE gut colonization was selected: 44 patients with spontaneous decolonization (DC) and 42 patients with persistent colonization (PC). The mean follow-up period was 24 months (IQR 14-33) in the DC group vs. 25 months (IQR 16-36) in the PC group (p = 0.478). Patient characteristics were similar between both groups. Colonization by other MDR microorganisms was high (44 patients, 51%) and slightly more common in the PC group (PC 60% vs. DC 43%, p = 0.139). The use of ceftazidime-avibactam was more common among the PC group (PC 33% vs. DC 14%, p = 0.041). We observed a higher percentage of antimicrobial therapy in the previous 30 days (PC 68% vs. DC 57%, p = 0.371) and 90 days (PC 81% vs. DC 82%, p = 0.353) in the PC group. Multivariable analysis showed that patients that have received ceftazidime-avibactam therapy (OR 4.9 95% CI [1.45-16.39], p = 0.010), and those colonized by other MDR microorganisms (OR 2.5, 95% CI [0.96-6.25], p = 0.060) presented a higher risk of PC. Ceftazidime-avibactam use and colonization by other MDR microorganisms might be associated with CRE persistent gut colonization.

In vitro neutralizing activity of BNT162b2 mRNA-induced antibodies against full B.1.351 SARS-CoV-2 variant

SARS-CoV-2 variation represents a serious challenge to current COVID-19 vaccines. Recent reports suggest that B.1.351 and other variants may escape the neutralization activity of the antibodies generated by current vaccines. Ninety-nine healthcare workers undertaking BNT162b2 mRNA vaccination were sampled at baseline, on the day of the second dose, and 14 days after the latter. Neutralization activity against SARS-CoV-2 B.1, B.1.1.7 and B.1.351 was investigated using a Vero-E6 model. Eleven of the study participants had prior infection with SARS-CoV-2. Neutralization titers against the B.1 and the B.1.1.7 variants were not statistically different and were significantly higher than titers against the B.1.351 variant across pre-exposed and non-pre-exposed vaccinated individuals (p < .01). While all vaccinated individuals presented neutralizing antibodies against B.1 and B 1.1.7 after the second dose, 14% were negative against B.1.351 and 76% had low titers (1/201/80). Pre-exposed vaccinated individuals showed higher titers than non-pre-exposed after the first (median titers of 1/387 versus 1/28, respectively) and the second doses (1/995 versus 1/703, respectively). As high as 72% of the pre-exposed vaccines presented titers >1/80 after a single dose, while only 11% of non-exposed vaccinated individuals had titers >1/80. BNT162b2 mRNA-induced antibodies show a lower in vitro neutralizing activity against B.1.351 variant compared to neutralization against B.1.1.7 or B.1 variants. Interestingly, for individuals pre-exposed to SARS-CoV-2, one dose of BNT162b2 mRNA may be adequate to produce neutralizing antibodies against B.1.1.7 and B.1, while two doses of BNT162b2 mRNA provide optimal neutralizing antibody response against B.1.351 too.

Enterovirus D68-associated respiratory and neurological illness in Spain, 2014-2018

During 2014, enterovirus D68 (EV-D68) outbreaks were described globally, causing severe respiratory diseases in children and, in some cases, subsequent paralysis. In this study, the type characterization of enterovirus (EV) detected in respiratory illnesses and the epidemiology and clinical association of EV-D68 infections in Spain over a five-year period were described. A total of 546 EV-positive samples from hospitalized patients with respiratory infections were included. EV-D68 was the most frequently detected type (46.6%, 191/410 typed EV). Other EV from species A (25.1%), B (27.8%) and C (0.5%) were also identified. EV-D68 infections were more associated with bronchitis while EV-A/B types were more frequent in upper respiratory illness (p < 0.01). EV-D68 was also detected in patients with neurological symptoms (nine meningitis/meningoencephalitis and eight acute flaccid paralysis cases). Phylogenetic analysis of 3′-VP1 region showed most Spanish EV-D68 sequences from 2014 to 2016 belonged to subclades B2/B3, as other American and European strains circulating during the same period. However, those detected in 2017 and 2018 clustered to the emerged subclade D1. In summary, different EV can cause respiratory infections but EV-D68 was the most prevalent, with several strains circulating in Spain at least since 2014. Association between EV-D68 infection and neurological disease was also described.

Prevalence and distribution of hepatitis C virus genotypes in Spain during the 2000-2015 period (the GEHEP 005 study)

We report the largest study on the prevalence and distribution of HCV genotypes in Spain (2000-2015), and we relate them with clinical, epidemiological and virological factors. Patients from 29 hospitals in 10 autonomous communities (Andalusia, Aragon, Castilla-Leon, Catalonia, Galicia, Canary Islands, Madrid Community, Valencian Community, Murcia Region and Basque Country) have been studied. Annual distribution of HCV genotypes and subtypes, as well as gender, age, transmission route, HIV and/or HBV coinfection, and treatment details were recorded. We included 48595 chronically HCV-infected patients with the following characteristics: median age 51 years (IQR, 44-58), 67.9% male, 19.1% HIV-coinfected, 23.5% HBV-coinfected. Parenteral transmission route was the most frequent (58.7%). Genotype distribution was 66.9% GT1 (24.9% subtype 1a and 37.9% subtype 1b), 2.8% GT2, 17.3% GT3, 11.4% GT4 and 0.1% GT5 and 0.02% GT6. LiPA was the most widely HCV genotyping test used (52.4%). HCV subtype 1a and genotypes 3 and 4 were closely associated with male gender, parenteral route of infection and HIV and HBV coinfection; in contrast, subtype 1b and genotype 2 were associated with female gender, nonparenteral route and mono-infection. Age was related to genotype distribution, and different patterns of distribution and biodiversity index were observed between different geographical areas. Finally, we describe how treatment and changes in transmission routes may have affected HCV genotype prevalence and distribution patterns. We present the most recent data on molecular epidemiology of hepatitis C virus in Spain. This study confirms that genotype distributions vary with age, sex, HIV and HBV coinfection and within geographical areas and epidemiological groups.

Sequence Analysis of In Vivo-Expressed HIV-1 Spliced RNAs Reveals the Usage of New and Unusual Splice Sites by Viruses of Different Subtypes

HIV-1 RNAs are generated through a complex splicing mechanism, resulting in a great diversity of transcripts, which are classified in three major categories: unspliced, singly spliced (SS), and doubly spliced (DS). Knowledge on HIV-1 RNA splicing in vivo and by non-subtype B viruses is scarce. Here we analyze HIV-1 RNA splice site usage in CD4+CD25+ lymphocytes from HIV-1-infected individuals through pyrosequencing. HIV-1 DS and SS RNAs were amplified by RT-PCR in 19 and 12 samples, respectively. 13,108 sequences from HIV-1 spliced RNAs, derived from viruses of five subtypes (A, B, C, F, G), were identified. In four samples, three of non-B subtypes, five 3' splice sites (3'ss) mapping to unreported positions in the HIV-1 genome were identified. Two, designated A4i and A4j, were used in 22% and 25% of rev RNAs in two viruses of subtypes B and A, respectively. Given their close proximity (one or two nucleotides) to A4c and A4d, respectively, they could be viewed as variants of these sites. Three 3'ss, designated A7g, A7h, and A7i, located 20, 32, and 18 nucleotides downstream of A7, respectively, were identified in a subtype C (A7g, A7h) and a subtype G (A7i) viruses, each in around 2% of nef RNAs. The new splice sites or variants of splice sites were associated with the usual sequence features of 3'ss. Usage of unusual 3'ss A4d, A4e, A5a, A7a, and A7b was also detected. A4f, previously identified in two subtype C viruses, was preferentially used by rev RNAs of a subtype C virus. These results highlight the great diversity of in vivo splice site usage by HIV-1 RNAs. The fact that four of five newly identified splice sites or variants of splice sites were detected in non-subtype B viruses allows anticipating an even greater diversity of HIV-1 splice site usage than currently known.

Identification of an HIV-1 BG Intersubtype Recombinant Form (CRF73_BG), Partially Related to CRF14_BG, Which Is Circulating in Portugal and Spain

HIV-1 exhibits a characteristically high genetic diversity, with the M group, responsible for the pandemic, being classified into nine subtypes, 72 circulating recombinant forms (CRFs) and numerous unique recombinant forms (URFs). Here we characterize the near full-length genome sequence of an HIV-1 BG intersubtype recombinant virus (X3208) collected in Galicia (Northwest Spain) which exhibits a mosaic structure coincident with that of a previously characterized BG recombinant virus (9601_01), collected in Germany and epidemiologically linked to Portugal, and different from currently defined CRFs. Similar recombination patterns were found in partial genome sequences from three other BG recombinant viruses, one newly derived, from a virus collected in Spain, and two retrieved from databases, collected in France and Portugal, respectively. Breakpoint coincidence and clustering in phylogenetic trees of these epidemiologically-unlinked viruses allow to define a new HIV-1 CRF (CRF73_BG). CRF73_BG shares one breakpoint in the envelope with CRF14_BG, which circulates in Portugal and Spain, and groups with it in a subtype B envelope fragment, but the greatest part of its genome does not appear to derive from CRF14_BG, although both CRFs share as parental strain the subtype G variant circulating in the Iberian Peninsula. Phylogenetic clustering of partial pol and env segments from viruses collected in Portugal and Spain with X3208 and 9691_01 indicates that CRF73_BG is circulating in both countries, with proportions of around 2–3% Portuguese database HIV-1 isolates clustering with CRF73_BG. The fact that an HIV-1 recombinant virus characterized ten years ago as a URF has been shown to represent a CRF suggests that the number of HIV-1 CRFs may be much greater than currently known.

Phylogeny and Phylogeography of a Recent HIV-1 Subtype F Outbreak among Men Who Have Sex with Men in Spain Deriving from a Cluster with a Wide Geographic Circulation in Western Europe

We recently reported the rapid expansion of an HIV-1 subtype F cluster among men who have sex with men (MSM) in the region of Galicia, Northwest Spain. Here we update this outbreak, analyze near full-length genomes, determine phylogenetic relationships, and estimate its origin. For this study, we used sequences of HIV-1 protease-reverse transcriptase and env V3 region, and for 17 samples, near full-length genome sequences were obtained. Phylogenetic analyses were performed via maximum likelihood. Locations and times of most recent common ancestors were estimated using Bayesian inference. Among samples analyzed by us, 100 HIV-1 F1 subsubtype infections of monophyletic origin were diagnosed in Spain, including 88 in Galicia and 12 in four other regions. Most viruses (n = 90) grouped in a subcluster (Galician subcluster), while 7 from Valladolid (Central Spain) grouped in another subcluster. At least 94 individuals were sexually-infected males and at least 71 were MSM. Seventeen near full-length genomes were uniformly of F1 subsubtype. Through similarity searches and phylogenetic analyses, we identified 18 viruses from four other Western European countries [Switzerland (n = 8), Belgium (n = 5), France (n = 3), and United Kingdom (n = 2)] and one from Brazil, from samples collected in 2005–2011, which branched within the subtype F cluster, outside of both Spanish subclusters, most of them corresponding to recently infected individuals. The most probable geographic origin and age of the Galician subcluster was Ferrol, Northwest Galicia, around 2007, while the Western European cluster probably emerged in Switzerland around 2002. In conclusion, a recently expanded HIV-1 subtype F cluster, the largest non-subtype B cluster reported in Western Europe, continues to spread among MSM in Spain; this cluster is part of a larger cluster with a wide geographic circulation in diverse Western European countries.

Cervical intraepithelial neoplasia grade 2 or worse in Galicia, Spain: HPV 16 prevalence and vaccination impact

INTRODUCTION

The etiology of cervical intraepithelial neoplasia grade 2 or worse (CIN2+) can influence the efficacy of Public Health preventive strategies. This study aimed to determine the high-risk papillomavirus (HR-HPV) prevalence in CIN2+ cases in unvaccinated women in Galicia (Spain), the expected impact of bivalent vaccination, and the distribution of HPV 16 in squamous lesions.

MATERIAL AND METHODS

Ninety-four histologically confirmed cases of CIN2+ (2009-2010) were retrospectively studied: 23 CIN2, 58 CIN3- squamous carcinoma in situ (CIN3-CIS), 5 adenocarcinoma in situ (AIS), and 8 invasive squamous cervical cancer (SCC). Linear Array HPV Genotyping Test (Roche Diagnostics, Mannheim, Germany) was performed on the cervical specimens. Bivalent vaccination impact was calculated, based on regional vaccination coverage data, local HR-HPV prevalence, and reported efficacy (direct and cross-protection) of the vaccine.

RESULTS

HR-HPV prevalence was 96.8%. The most frequent genotypes were HPV 16 (48.8-58.2%) and HPV 31 (9.3%-12.1%), considering single infections or single-multiple infections, respectively (hierarchical attribution). In squamous lesions, HPV 16 prevalence in women younger than 45 years of age increased in severe lesions (CIN3-CIS/SCC, OR 4.2), and was higher than in older women (OR 5.5). The vaccine could reduce the cumulative incidence of CIN2+ by 50.6% (direct protection), or by 62.7% (direct and cross-protection).

CONCLUSION

HPV vaccination could have a great impact in women younger than 45 years of age due to the high prevalence of HPV 16 in their lesions.

Predominance of CXCR4 tropism in HIV-1 CRF14_BG strains from newly diagnosed infections

OBJECTIVES

R5-tropic viruses are associated with HIV-1 transmission and predominate during the early stages of infection. X4-tropic populations have been detected in ~50% of patients with late-stage disease infected with subtype B viruses. In this study, we compared the frequency of X4 tropism in individuals infected with HIV-1 CRF14_BG viruses, which have a V3 loop of subtype B, with a control group of individuals infected with subtype B viruses.

METHODS

Sixty-three individuals infected with HIV-1 CRF14_BG (n = 31) or subtype B (n = 32) were studied. Similar proportions of newly diagnosed and chronically infected individuals were included in the subtype B and CRF14_BG groups. V3 sequences were obtained and coreceptor tropism was predicted using the Geno2pheno[coreceptor] algorithm. V3 net charge and 11/25 rules were also used for coreceptor prediction.

RESULTS

Overall, X4 tropism was more frequent among individuals infected with CRF14_BG viruses (87.1%) than subtype B viruses (34.3%), a difference that was statistically highly significant (P = 0.00001). Importantly, the frequencies among newly diagnosed individuals were 90% and 13.3%, respectively (P = 0.0007). Characteristic amino acids in the V3 loop (T13, M14, V19 and W20) were identified at higher frequencies in CRF14_BG viruses (54%) than subtype B viruses (0%; P < 0.000001).

CONCLUSIONS

CRF14_BG is the genetic form with the highest proportion of X4-tropic viruses reported to date in newly diagnosed and chronic infections. This suggests high pathogenicity for CRF14_BG viruses, potentially leading to rapid disease progression. CCR5 antagonists will be ineffective in most CRF14_BG-infected patients, even at early stages of infection.

Human papillomavirus infection in a male population attending a sexually transmitted infection service

OBJECTIVE

Human Papillomavirus (HPV) infection in men may produce cancer and other major disorders. Men play an important role in the transmission of the virus and act as a reservoir. The aim of this study was to determine the HPV-genotypes and their prevalence in a group of men attending a Sexually Transmitted Infection service. PATIENTS AND SAMPLES: Between July 2002 and June 2011, 1392 balanopreputial, 435 urethral, 123 anal, and 67 condyloma lesions from 1551 men with a mean age of 35.8±11.3 years old (range: 17-87) were collected for HPV-DNA testing.

METHODS

A fragment of the L1-gene and a fragment of the E6/E7-genes were amplified by PCR. Positive samples were typed by hybridization.

RESULTS

The HPV genome was detected in 36.9% (486/1318) balanopreputial and in 24.9% (101/405) urethral (p<0.0001) swabs from 38.1% (538) of 1469 men. Co-infections were present in 5.4% (80/1469) of cases. HPV was found in 43.9% (373/850) of men younger than 35 vs. 31.7% (187/589) of men aged >35. HPV was found in 59.4% (104) of 165 men with lesions (macroscopic or positive peniscopy), and in 22.8% (61/267) without clinical alterations. HPV was also detected in 71.4% (40/56) men with condylomata and in 58.7% (64/109) of men with positive peniscopy.

CONCLUSIONS

HPV prevalence in men was high and decreased with age. HPV was found more frequently in balanopreputial than in urethral swabs. There was a low rate of co-infections. Low-risk HPV vaccine genotypes were the most recurrent especially in younger. Although HPV has been associated with clinical alterations, it was also found in men without any clinical presentation. Inclusion of men in the national HPV vaccination program may reduce their burden of HPV-related disease and reduce transmission of the virus to non-vaccinated women.

Identification of a new HIV type 1 circulating BF intersubtype recombinant form (CRF47_BF) in Spain

We report the identification of a new HIV-1 circulating recombinant form (CRF47_BF) derived from subtypes B and F. It was initially identified in protease-reverse transcriptase sequences from nine individuals from three separate regions of Spain who acquired HIV-1 infection via sexual contact. All nine sequences formed a strongly supported phylogenetic cluster, branching apart from all known CRFs, and in bootscan analyses were BF mosaics with two coincident breakpoints. Two epidemiologically unlinked viruses were sequenced in near full-length genomes, which exhibited identical mosaic structures, with 16 intersubtype breakpoints in a genome predominantly of subtype B. Subtype F segments of the new CRF failed to cluster with any of the near full-length genome subtype F sequences available in public databases. Recent dates of HIV-1 diagnoses and short genetic distances suggest a recent origin of this CRF. This is the tenth reported CRF_BF, the first apparently having originated outside of South America.