VIPERA: Viral Intra-Patient Evolution Reporting and Analysis.

Viral mutations within patients nurture the adaptive potential of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) during chronic infections, which are a potential source of variants of concern. However, there is no integrated framework for the evolutionary analysis of intra-patient SARS-CoV-2 serial samples. Herein, we describe Viral Intra-Patient Evolution Reporting and Analysis (VIPERA), a new software that integrates the evaluation of the intra-patient ancestry of SARS-CoV-2 sequences with the analysis of evolutionary trajectories of serial sequences from the same viral infection. We have validated it using positive and negative control datasets and have successfully applied it to a new case, which revealed population dynamics and evidence of adaptive evolution. VIPERA is available under a free software license at https://github.com/PathoGenOmics-Lab/VIPERA.

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.

Role of the first WHO mutation catalogue in the diagnosis of antibiotic resistance in Mycobacterium tuberculosis in the Valencia Region, Spain: a retrospective genomic analysis.

BACKGROUND: In June, 2021, WHO published the most complete catalogue to date of resistance-conferring mutations in Mycobacterium tuberculosis. Here, we aimed to assess the performance of genome-based antimicrobial resistance prediction using the catalogue and its potential for improving diagnostics in a real low-burden setting. METHODS: In this retrospective population-based genomic study M tuberculosis isolates were collected from 25 clinical laboratories in the low-burden setting of the Valencia Region, Spain. Culture-positive tuberculosis cases reported by regional public health authorities between Jan 1, 2014, and Dec 31, 2016, were included. The drug resistance profiles of these isolates were predicted by the genomic identification, via whole-genome sequencing (WGS), of the high-confidence resistance-causing variants included in the catalogue and compared with the phenotype. We determined the minimum inhibitory concentration (MIC) of the isolates with discordant resistance profiles using the resazurin microtitre assay. FINDINGS: WGS was performed on 785 M tuberculosis complex culture-positive isolates, and the WGS resistance prediction sensitivities were: 85·4% (95% CI 70·8-94·4) for isoniazid, 73·3% (44·9-92·2) for rifampicin, 50·0% (21·1-78·9) for ethambutol, and 57·1% (34·0-78·2) for pyrazinamide; all specificities were more than 99·6%. Sensitivity values were lower than previously reported, but the overall pan-susceptibility accuracy was 96·4%. Genotypic analysis revealed that four phenotypically susceptible isolates carried mutations (rpoB Leu430Pro and rpoB Ile491Phe for rifampicin and fabG1 Leu203Leu for isoniazid) known to give borderline resistance in standard phenotypic tests. Additionally, we identified three putative resistance-associated mutations (inhA Ser94Ala, katG Leu48Pro, and katG Gly273Arg for isoniazid) in samples with substantially higher MICs than those of susceptible isolates. Combining both genomic and phenotypic data, in accordance with the WHO diagnostic guidelines, we could detect two new multidrug-resistant cases. Additionally, we detected 11 (1·6%) of 706 isolates to be monoresistant to fluoroquinolone, which had been previously undetected. INTERPRETATION: We showed that the WHO catalogue enables the detection of resistant cases missed in phenotypic testing in a low-burden region, thus allowing for better patient-tailored treatment. We also identified mutations not included in the catalogue, relevant at the local level. Evidence from this study, together with future updates of the catalogue, will probably lead in the future to the partial replacement of culture testing with WGS-based drug susceptibility testing in our setting. FUNDING: European Research Council and the Spanish Ministerio de Ciencia.

Monitoring of First-line Drug Resistance Mutations Outside the Scope of Xpert MTB/RIF Ultra is Needed for Successful Control of DR-TB in Southern Mozambique.

Multidrug-resistant(MDR) tuberculosis in Southern Africa is of great concern, exacerbated by the spread of a clone harboring a mutation missed by Xpert Ultra. In Southern Mozambique, the presence of such mutation and rising cases of non-MDR isoniazid resistance highlights the need to ensure accurate detection of antimicrobial-resistance in the country.

Tracking SARS-CoV-2 introductions in Mozambique using pandemic-scale phylogenies: a retrospective observational study.

BACKGROUND: From the start of the SARS-CoV-2 outbreak, global sequencing efforts have generated an unprecedented amount of genomic data. Nonetheless, unequal sampling between high-income and low-income countries hinders the implementation of genomic surveillance systems at the global and local level. Filling the knowledge gaps of genomic information and understanding pandemic dynamics in low-income countries is essential for public health decision making and to prepare for future pandemics. In this context, we aimed to discover the timing and origin of SARS-CoV-2 variant introductions in Mozambique, taking advantage of pandemic-scale phylogenies. METHODS: We did a retrospective, observational study in southern Mozambique. Patients from Manhiça presenting with respiratory symptoms were recruited, and those enrolled in clinical trials were excluded. Data were included from three sources: (1) a prospective hospital-based surveillance study (MozCOVID), recruiting patients living in Manhiça, attending the Manhiça district hospital, and fulfilling the criteria of suspected COVID-19 case according to WHO; (2) symptomatic and asymptomatic individuals with SARS-CoV-2 infection recruited by the National Surveillance system; and (3) sequences from SARS-CoV-2-infected Mozambican cases deposited on the Global Initiative on Sharing Avian Influenza Data database. Positive samples amenable for sequencing were analysed. We used Ultrafast Sample placement on Existing tRees to understand the dynamics of beta and delta waves, using available genomic data. This tool can reconstruct a phylogeny with millions of sequences by efficient sample placement in a tree. We reconstructed a phylogeny (~7·6 million sequences) adding new and publicly available beta and delta sequences. FINDINGS: A total of 5793 patients were recruited between Nov 1, 2020, and Aug 31, 2021. During this time, 133 328 COVID-19 cases were reported in Mozambique. 280 good quality new SARS-CoV-2 sequences were obtained after the inclusion criteria were applied and an additional 652 beta (B.1.351) and delta (B.1.617.2) public sequences were included from Mozambique. We evaluated 373 beta and 559 delta sequences. We identified 187 beta introductions (including 295 sequences), divided in 42 transmission groups and 145 unique introductions, mostly from South Africa, between August, 2020 and July, 2021. For delta, we identified 220 introductions (including 494 sequences), with 49 transmission groups and 171 unique introductions, mostly from the UK, India, and South Africa, between April and November, 2021. INTERPRETATION: The timing and origin of introductions suggests that movement restrictions effectively avoided introductions from non-African countries, but not from surrounding countries. Our results raise questions about the imbalance between the consequences of restrictions and health benefits. This new understanding of pandemic dynamics in Mozambique can be used to inform public health interventions to control the spread of new variants. FUNDING: European and Developing Countries Clinical Trials, European Research Council, Bill & Melinda Gates Foundation, and Agència de Gestió d'Ajuts Universitaris i de Recerca.

Lack of methoxy-mycolates characterizes the geographically restricted lineage 7 of Mycobacterium tuberculosis complex.

Lineage 7 (L7) emerged in the phylogeny of the Mycobacterium tuberculosis complex (MTBC) subsequent to the branching of 'ancient' lineage 1 and prior to the Eurasian dispersal of 'modern' lineages 2, 3 and 4. In contrast to the major MTBC lineages, the current epidemiology suggests that prevalence of L7 is highly confined to the Ethiopian population, or when identified outside of Ethiopia, it has mainly been in patients of Ethiopian origin. To search for microbiological factors that may contribute to its restricted distribution, we compared the genome of L7 to the genomes of globally dispersed MTBC lineages. The frequency of predicted functional mutations in L7 was similar to that documented in other lineages. These include mutations characteristic of modern lineages - such as constitutive expression of nitrate reductase - as well as mutations in the VirS locus that are commonly found in ancient lineages. We also identified and characterized multiple lineage-specific mutations in L7 in biosynthesis pathways of cell wall lipids, including confirmed deficiency of methoxy-mycolic acids due to a stop-gain mutation in the mmaA3 gene that encodes a methoxy-mycolic acid synthase. We show that the abolished biosynthesis of methoxy-mycolates of L7 alters the cell structure and colony morphology on selected growth media and impacts biofilm formation. The loss of these mycolic acid moieties may change the host-pathogen dynamic for L7 isolates, explaining the limited geographical distribution of L7 and contributing to further understanding the spread of MTBC lineages across the globe.

Epidemiological cluster identification using multiple data sources: an approach using logistic regression.

In the management of infectious disease outbreaks, grouping cases into clusters and understanding their underlying epidemiology are fundamental tasks. In genomic epidemiology, clusters are typically identified either using pathogen sequences alone or with sequences in combination with epidemiological data such as location and time of collection. However, it may not be feasible to culture and sequence all pathogen isolates, so sequence data may not be available for all cases. This presents challenges for identifying clusters and understanding epidemiology, because these cases may be important for transmission. Demographic, clinical and location data are likely to be available for unsequenced cases, and comprise partial information about their clustering. Here, we use statistical modelling to assign unsequenced cases to clusters already identified by genomic methods, assuming that a more direct method of linking individuals, such as contact tracing, is not available. We build our model on pairwise similarity between cases to predict whether cases cluster together, in contrast to using individual case data to predict the cases' clusters. We then develop methods that allow us to determine whether a pair of unsequenced cases are likely to cluster together, to group them into their most probable clusters, to identify which are most likely to be members of a specific (known) cluster, and to estimate the true size of a known cluster given a set of unsequenced cases. We apply our method to tuberculosis data from Valencia, Spain. Among other applications, we find that clustering can be predicted successfully using spatial distance between cases and whether nationality is the same. We can identify the correct cluster for an unsequenced case, among 38 possible clusters, with an accuracy of approximately 35 %, higher than both direct multinomial regression (17 %) and random selection (< 5 %).

Deciphering the Tangible Spatio-Temporal Spread of a 25-Year Tuberculosis Outbreak Boosted by Social Determinants.

Outbreak strains of Mycobacterium tuberculosis are promising candidates as targets in the search for intrinsic determinants of transmissibility, as they are responsible for many cases with sustained transmission; however, the use of low-resolution typing methods and restricted geographical investigations represent flaws in assessing the success of long-lived outbreak strains. We can now address the nature of outbreak strains by combining large genomic data sets and phylodynamic approaches. We retrospectively sequenced the whole genome of representative samples assigned to an outbreak circulating in the Canary Islands (the GC strain) since 1993, which accounts for ~20% of local tuberculosis cases. We selected a panel of specific single nucleotide polymorphism (SNP) markers for an in-silico search for additional outbreak-related sequences within publicly available tuberculosis genomic data. Using this information, we inferred the origin, spread, and epidemiological parameters of the GC strain. Our approach allowed us to accurately trace the historical and more recent dispersion of the GC strain. We provide evidence of a highly successful nature within the Canarian archipelago but limited expansion abroad. Estimation of epidemiological parameters from genomic data disagree with a distinctive biology of the GC strain. With the increasing availability of genomic data allowing for the accurate inference of strain spread and critical epidemiological parameters, we can now revisit the link between Mycobacterium tuberculosis genotypes and transmission, as is routinely carried out for SARS-CoV-2 variants of concern. We demonstrate that social determinants rather than intrinsically higher bacterial transmissibility better explain the success of the GC strain. Importantly, our approach can be used to trace and characterize strains of interest worldwide. IMPORTANCE Infectious disease outbreaks represent a significant problem for public health. Tracing outbreak expansion and understanding the main factors behind emergence and persistence remain critical to effective disease control. Our study allows researchers and public health authorities to use Whole-Genome Sequencing-based methods to trace outbreaks, and shows how available epidemiological information helps to evaluate the factors underpinning outbreak persistence. Taking advantage of all the freely available information placed in public repositories, researchers can accurately establish the expansion of an outbreak beyond original boundaries, and determine the potential risk of a strain to inform health authorities which, in turn, can define target strategies to mitigate expansion and persistence. Finally, we show the need to evaluate strain transmissibility in different geographic contexts to unequivocally associate spread to local or pathogenic factors, an important lesson taken from genomic surveillance of SARS-CoV-2.

Evaluation of two different concentration methods for surveillance of human viruses in sewage and their effects on SARS-CoV-2 sequencing.

During the current COVID-19 pandemic, wastewater-based epidemiology (WBE) emerged as a reliable strategy both as a surveillance method and a way to provide an overview of the SARS-CoV-2 variants circulating among the population. Our objective was to compare two different concentration methods, a well-established aluminum-based procedure (AP) and the commercially available Maxwell® RSC Enviro Wastewater TNA Kit (TNA) for human enteric virus, viral indicators and SARS-CoV-2 surveillance. Additionally, both concentration methods were analyzed for their impact on viral infectivity, and nucleic acids obtained from each method were also evaluated by massive sequencing for SARS-CoV-2. The percentage of SARS-CoV-2 positive samples using the AP method accounted to 100 %, 83.3 %, and 33.3 % depending on the target region while 100 % positivity for these same three target regions was reported using the TNA procedure. The concentrations of norovirus GI, norovirus GII and HEV using the TNA method were significantly greater than for the AP method while no differences were reported for rotavirus, astrovirus, crAssphage and PMMoV. Furthermore, TNA kit in combination with the Artic v4 primer scheme yields the best SARS-CoV-2 sequencing results. Regarding impact on infectivity, the concentration method used by the TNA kit showed near-complete lysis of viruses. Our results suggest that although the performance of the TNA kit was higher than that of the aluminum procedure, both methods are suitable for the analysis of enveloped and non-enveloped viruses in wastewater by molecular methods.

The mineralocorticoid receptor modulates timing and location of genomic binding by glucocorticoid receptor in response to synthetic glucocorticoids in keratinocytes.

Glucocorticoids (GCs) exert potent antiproliferative and anti-inflammatory properties, explaining their therapeutic efficacy for skin diseases. GCs act by binding to the GC receptor (GR) and the mineralocorticoid receptor (MR), co-expressed in classical and non-classical targets including keratinocytes. Using knockout mice, we previously demonstrated that GR and MR exert essential nonoverlapping functions in skin homeostasis. These closely related receptors may homo- or heterodimerize to regulate transcription, and theoretically bind identical GC-response elements (GRE). We assessed the contribution of MR to GR genomic binding and the transcriptional response to the synthetic GC dexamethasone (Dex) using control (CO) and MR knockout (MREKO ) keratinocytes. GR chromatin immunoprecipitation (ChIP)-seq identified peaks common and unique to both genotypes upon Dex treatment (1 h). GREs, AP-1, TEAD, and p53 motifs were enriched in CO and MREKO peaks. However, GR genomic binding was 35% reduced in MREKO , with significantly decreased GRE enrichment, and reduced nuclear GR. Surface plasmon resonance determined steady state affinity constants, suggesting preferred dimer formation as MR-MR > GR-MR ~ GR-GR; however, kinetic studies demonstrated that GR-containing dimers had the longest lifetimes. Despite GR-binding differences, RNA-seq identified largely similar subsets of differentially expressed genes in both genotypes upon Dex treatment (3 h). However, time-course experiments showed gene-dependent differences in the magnitude of expression, which correlated with earlier and more pronounced GR binding to GRE sites unique to CO including near Nr3c1. Our data show that endogenous MR has an impact on the kinetics and differential genomic binding of GR, affecting the time-course, specificity, and magnitude of GC transcriptional responses in keratinocytes.

Cumulative incidence of SARS-CoV-2 infection in the general population of the Valencian Community (Spain) after the surge of the Omicron BA.1 variant.

Studies investigating the cumulative incidence of and immune status against SARS-CoV-2 infection provide valuable information for shaping public health decision-making. A cross-sectional study on 935 participants, conducted in the Valencian Community (VC), measuring anti-SARS-CoV-2-receptor binding domain-RBD-total antibodies and anti-Nucleocapsid (N)-IgGs via electrochemiluminescence assays. Quantitation of neutralizing antibodies (NtAb) against ancestral and Omicron BA.1 and BA.2 variants and enumeration of SARS-CoV-2-S specific-IFNγ-producing CD4+ and CD8+ T cells was performed in 100 and 137 participants, respectively. The weighted cumulative incidence was 51.9% (95% confidence interval [CI]: 48.7-55.1) and was inversely related to age. Anti-RBD total antibodies were detected in 97% of participants; vaccinated and SARS-CoV-2-experienced (VAC-ex; n = 442) presented higher levels (p < 0.001) than vaccinated/naïve (VAC-n; n = 472) and nonvaccinated/experienced (UNVAC-ex; n = 63) subjects. Antibody levels correlated inversely with time elapsed since last vaccine dose in VAC-n (Rho, -0.52; p < 0.001) but not in VAC-ex (rho -0.02; p = 0.57). Heterologous booster shots resulted in increased anti-RBD antibody levels compared with homologous schedules in VAC-n, but not in VAC-ex. NtAbs against Omicron BA.1 were detected in 94%, 75%, and 50% of VAC-ex, VAC-n and UNVAC-ex groups, respectively. For Omicron BA.2, the figures were 97%, 84%, and 40%, respectively. SARS-CoV-2-S-reactive IFN-γ T cells were detected in 73%, 75%, and 64% of VAC-ex, VAC-n and UNVAC-ex, respectively. Median frequencies for both T-cell subsets were comparable across groups. In summary, by April 2022, around half of the VC population had been infected with SARS-CoV-2 and, due to extensive vaccination, displayed hybrid immunity.

Signatures of transmission in within-host M. tuberculosis variation.

Background: Because M. tuberculosis evolves slowly, transmission clusters often contain multiple individuals with identical consensus genomes, making it difficult to reconstruct transmission chains. Finding additional sources of shared M. tuberculosis variation could help overcome this problem. Previous studies have reported M. tuberculosis diversity within infected individuals; however, whether within-host variation improves transmission inferences remains unclear. Methods: To evaluate the transmission information present in within-host M. tuberculosis variation, we re-analyzed publicly available sequence data from three household transmission studies, using household membership as a proxy for transmission linkage between donor-recipient pairs. Findings: We found moderate levels of minority variation present in M. tuberculosis sequence data from cultured isolates that varied significantly across studies (mean: 6, 7, and 170 minority variants above a 1% minor allele frequency threshold, outside of PE/PPE genes). Isolates from household members shared more minority variants than did isolates from unlinked individuals in the three studies (mean 98 shared minority variants vs. 10; 0.8 vs. 0.2, and 0.7 vs. 0.2, respectively). Shared within-host variation was significantly associated with household membership (OR: 1.51 [1.30,1.71], for one standard deviation increase in shared minority variants). Models that included shared within-host variation improved the accuracy of predicting household membership in all three studies as compared to models without within-host variation (AUC: 0.95 versus 0.92, 0.99 versus 0.95, and 0.93 versus 0.91). Interpretation: Within-host M. tuberculosis variation persists through culture and could enhance the resolution of transmission inferences. The substantial differences in minority variation recovered across studies highlights the need to optimize approaches to recover and incorporate within-host variation into automated phylogenetic and transmission inference. Funding: NIAID: 5K01AI173385.

Prevalence of SARS-CoV-2 and co-infection with malaria during the first wave of the pandemic (the Burkina Faso case).

Africa accounts for 1.5% of the global coronavirus disease 2019 (COVID-19) cases and 2.7% of deaths, but this low incidence has been partly attributed to the limited testing capacity in most countries. In addition, the population in many African countries is at high risk of infection with endemic infectious diseases such as malaria. Our aim is to determine the prevalence and circulation of SARS-CoV-2 variants, and the frequency of co-infection with the malaria parasite. We conducted serological tests and microscopy examinations on 998 volunteers of different ages and sexes in a random and stratified population sample in Burkina-Faso. In addition, nasopharyngeal samples were taken for RT-qPCR of SARS-CoV-2 and for whole viral genome sequencing. Our results show a 3.2 and a 2.5% of SARS-CoV-2 seroprevalence and PCR positivity; and 22% of malaria incidence, over the sampling period, with marked differences linked to age. Importantly, we found 8 cases of confirmed co-infection and 11 cases of suspected co-infection mostly in children and teenagers. Finally, we report the genome sequences of 13 SARS-CoV-2 isolates circulating in Burkina Faso at the time of analysis, assigned to lineages A.19, A.21, B.1.1.404, B.1.1.118, B.1 and grouped into clades; 19B, 20A, and 20B. This is the first population-based study about SARS-CoV-2 and malaria in Burkina Faso during the first wave of the pandemic, providing a relevant estimation of the real prevalence of SARS-CoV-2 and variants circulating in this Western African country. Besides, it highlights the non-negligible frequency of co-infection with malaria in African communities.

Epidemiological and Genomic Analysis of a Large SARS-CoV-2 Outbreak in a Long-Term Care Facility in Catalonia, Spain.

Limiting outbreaks in long-term care facilities (LTCFs) is a cornerstone strategy to avoid an excess of COVID-19-related morbidity and mortality and to reduce its burden on the health system. We studied a large outbreak that occurred at an LTCF, combining methods of classical and genomic epidemiology analysis. The outbreak lasted for 31 days among residents, with an attack rate of 98% and 57% among residents and staff, respectively. The case fatality rate among residents was 16% (n = 15). Phylogenetic analysis of 59 SARS-CoV-2 isolates revealed the presence of two closely related viral variants in all cases (B.1.177 lineage), revealing a far more complex outbreak than initially thought and suggesting an initial spread driven by staff members. In turn, our results suggest that resident relocations to mitigate viral spread might have increased the risk of infection for staff members, creating secondary chains of transmission that were responsible for prolonging the outbreak. Our results highlight the importance of considering unnoticed chains of transmission early during an outbreak and making an adequate use and interpretation of diagnostic tests. Outbreak containment measures should be carefully tailored to each LTCF. IMPORTANCE The impact of COVID-19 on long-term care facilities (LTCFs) has been disproportionately large due to the high frailty of the residents. Here, we report epidemiological and genomic findings of a large outbreak that occurred at an LTCF, which ultimately affected almost all residents and nearly half of staff members. We found that the outbreak was initially driven by staff members; however, later resident relocation to limit the outbreak resulted in transmission from residents to staff members, evidencing the complexity and different phases of the outbreak. The phylogenetic analysis of SARS-CoV-2 isolates indicated that two closely related variants were responsible for the large outbreak. Our study highlights the importance of combining methods of classical and genomic epidemiology to take appropriate outbreak containment measures in LTCFs.

Correction: The structural role of SARS-CoV-2 genetic background in the emergence and success of spike mutations: The case of the spike A222V mutation.

[This corrects the article DOI: 10.1371/journal.ppat.1010631.].

Fine-grain population structure and transmission patterns of Mycobacterium tuberculosis in southern Mozambique, a high TB/HIV burden area.

Genomic studies of the Mycobacterium tuberculosis complex (MTBC) might shed light on the dynamics of its transmission, especially in high-burden settings, where recent outbreaks are embedded in the complex natural history of the disease. To this end, we conducted a 1 year prospective surveillance-based study in Mozambique. We applied whole-genome sequencing (WGS) to 295 positive cultures. We fully characterized MTBC isolates by phylogenetics and dating evaluation, and carried out a molecular epidemiology analysis to investigate further associations with pre-defined transmission risk factors. The majority of strains (49.5%, 136/275) belonged to lineage (L) 4; 57.8 % of them (159/275) were in genomic transmission clusters (cut-off 5 SNPs), and a strikingly high proportion (45.5%) shared an identical genotype (0 SNP pairwise distance). We found two 'likely endemic' clades, comprising 67 strains, belonging to L1.2, which dated back to the late 19th century and were associated with recent spread among people living with human immunodeficiency virus (PLHIV). We describe for the first time the population structure of MTBC in our region, a high tuberculosis (TB)/HIV burden area. Clustering analysis revealed an unforeseen pattern of spread and high rates of progression to active TB, suggesting weaknesses in TB control activities. The long-term presence of local strains in Mozambique, which were responsible for large transmission among HIV/TB-coinfected patients, calls into question the role of HIV in TB transmission.

The structural role of SARS-CoV-2 genetic background in the emergence and success of spike mutations: The case of the spike A222V mutation.

The S:A222V point mutation, within the G clade, was characteristic of the 20E (EU1) SARS-CoV-2 variant identified in Spain in early summer 2020. This mutation has since reappeared in the Delta subvariant AY.4.2, raising questions about its specific effect on viral infection. We report combined serological, functional, structural and computational studies characterizing the impact of this mutation. Our results reveal that S:A222V promotes an increased RBD opening and slightly increases ACE2 binding as compared to the parent S:D614G clade. Finally, S:A222V does not reduce sera neutralization capacity, suggesting it does not affect vaccine effectiveness.

Population-based sequencing of Mycobacterium tuberculosis reveals how current population dynamics are shaped by past epidemics.

Transmission is a driver of tuberculosis (TB) epidemics in high-burden regions, with assumed negligible impact in low-burden areas. However, we still lack a full characterization of transmission dynamics in settings with similar and different burdens. Genomic epidemiology can greatly help to quantify transmission, but the lack of whole genome sequencing population-based studies has hampered its application. Here, we generate a population-based dataset from Valencia region and compare it with available datasets from different TB-burden settings to reveal transmission dynamics heterogeneity and its public health implications. We sequenced the whole genome of 785 Mycobacterium tuberculosis strains and linked genomes to patient epidemiological data. We use a pairwise distance clustering approach and phylodynamic methods to characterize transmission events over the last 150 years, in different TB-burden regions. Our results underscore significant differences in transmission between low-burden TB settings, i.e., clustering in Valencia region is higher (47.4%) than in Oxfordshire (27%), and similar to a high-burden area as Malawi (49.8%). By modeling times of the transmission links, we observed that settings with high transmission rate are associated with decades of uninterrupted transmission, irrespective of burden. Together, our results reveal that burden and transmission are not necessarily linked due to the role of past epidemics in the ongoing TB incidence, and highlight the need for in-depth characterization of transmission dynamics and specifically tailored TB control strategies.