Redefining the treponemal history through pre-Columbian genomes from Brazil.

The origins of treponemal diseases have long remained unknown, especially considering the sudden onset of the first syphilis epidemic in the late 15th century in Europe and its hypothesized arrival from the Americas with Columbus' expeditions1,2. Recently, ancient DNA evidence has revealed various treponemal infections circulating in early modern Europe and colonial-era Mexico3-6. However, there has been to our knowledge no genomic evidence of treponematosis recovered from either the Americas or the Old World that can be reliably dated to the time before the first trans-Atlantic contacts. Here, we present treponemal genomes from nearly 2,000-year-old human remains from Brazil. We reconstruct four ancient genomes of a prehistoric treponemal pathogen, most closely related to the bejel-causing agent Treponema pallidum endemicum. Contradicting the modern day geographical niche of bejel in the arid regions of the world, the results call into question the previous palaeopathological characterization of treponeme subspecies and showcase their adaptive potential. A high-coverage genome is used to improve molecular clock date estimations, placing the divergence of modern T. pallidum subspecies firmly in pre-Columbian times. Overall, our study demonstrates the opportunities within archaeogenetics to uncover key events in pathogen evolution and emergence, paving the way to new hypotheses on the origin and spread of treponematoses.

Spread of a SARS-CoV-2 variant through Europe in the summer of 2020.

Following its emergence in late 2019, the spread of SARS-CoV-21,2 has been tracked by phylogenetic analysis of viral genome sequences in unprecedented detail3-5. Although the virus spread globally in early 2020 before borders closed, intercontinental travel has since been greatly reduced. However, travel within Europe resumed in the summer of 2020. Here we report on a SARS-CoV-2 variant, 20E (EU1), that was identified in Spain in early summer 2020 and subsequently spread across Europe. We find no evidence that this variant has increased transmissibility, but instead demonstrate how rising incidence in Spain, resumption of travel, and lack of effective screening and containment may explain the variant's success. Despite travel restrictions, we estimate that 20E (EU1) was introduced hundreds of times to European countries by summertime travellers, which is likely to have undermined local efforts to minimize infection with SARS-CoV-2. Our results illustrate how a variant can rapidly become dominant even in the absence of a substantial transmission advantage in favourable epidemiological settings. Genomic surveillance is critical for understanding how travel can affect transmission of SARS-CoV-2, and thus for informing future containment strategies as travel resumes.

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

The dynamics of severe acute respiratory syndrome coronavirus 2 (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-BA.1 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 United States became increasingly significant. The number of detected introductions varied from 96 and 101 for Alpha and Delta to 39 for Omicron-BA.1. Most of these introductions left a low number of descendants (<10), suggesting a limited impact on the evolution of the pandemic in Galicia. 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.

Evolutionary Processes in the Emergence and Recent Spread of the Syphilis Agent, Treponema pallidum.

The incidence of syphilis has risen worldwide in the last decade in spite of being an easily treated infection. The causative agent of this sexually transmitted disease is the bacterium Treponema pallidum subspecies pallidum (TPA), very closely related to subsp. pertenue (TPE) and endemicum (TEN), responsible for the human treponematoses yaws and bejel, respectively. Although much focus has been placed on the question of the spatial and temporary origins of TPA, the processes driving the evolution and epidemiological spread of TPA since its divergence from TPE and TEN are not well understood. Here, we investigate the effects of recombination and selection as forces of genetic diversity and differentiation acting during the evolution of T. pallidum subspecies. Using a custom-tailored procedure, named phylogenetic incongruence method, with 75 complete genome sequences, we found strong evidence for recombination among the T. pallidum subspecies, involving 12 genes and 21 events. In most cases, only one recombination event per gene was detected and all but one event corresponded to intersubspecies transfers, from TPE/TEN to TPA. We found a clear signal of natural selection acting on the recombinant genes, which is more intense in their recombinant regions. The phylogenetic location of the recombination events detected and the functional role of the genes with signals of positive selection suggest that these evolutionary processes had a key role in the evolution and recent expansion of the syphilis bacteria and significant implications for the selection of vaccine candidates and the design of a broadly protective syphilis vaccine.

Enhancing SARS-CoV-2 Surveillance through Regular Genomic Sequencing in Spain: The RELECOV Network.

Millions of SARS-CoV-2 whole genome sequences have been generated to date. However, good quality data and adequate surveillance systems are required to contribute to meaningful surveillance in public health. In this context, the network of Spanish laboratories for coronavirus (RELECOV) was created with the main goal of promoting actions to speed up the detection, analyses, and evaluation of SARS-CoV-2 at a national level, partially structured and financed by an ECDC-HERA-Incubator action (ECDC/GRANT/2021/024). A SARS-CoV-2 sequencing quality control assessment (QCA) was developed to evaluate the network's technical capacity. QCA full panel results showed a lower hit rate for lineage assignment compared to that obtained for variants. Genomic data comprising 48,578 viral genomes were studied and evaluated to monitor SARS-CoV-2. The developed network actions showed a 36% increase in sharing viral sequences. In addition, analysis of lineage/sublineage-defining mutations to track the virus showed characteristic mutation profiles for the Delta and Omicron variants. Further, phylogenetic analyses strongly correlated with different variant clusters, obtaining a robust reference tree. The RELECOV network has made it possible to improve and enhance the genomic surveillance of SARS-CoV-2 in Spain. It has provided and evaluated genomic tools for viral genome monitoring and characterization that make it possible to increase knowledge efficiently and quickly, promoting the genomic surveillance of SARS-CoV-2 in Spain.

Systematic Genomic and Clinical Analysis of Severe Acute Respiratory Syndrome Coronavirus 2 Reinfections and Recurrences Involving the Same Strain.

Estimates of the burden of severe acute respiratory syndrome coronavirus 2 reinfections are limited by the scarcity of population-level studies incorporating genomic support. We conducted a systematic study of reinfections in Madrid, Spain, supported by genomic viral analysis and host genetic analysis, to cleanse laboratory errors and to discriminate between reinfections and recurrences involving the same strain. Among the 41,195 cases diagnosed (March 2020-March 2021), 93 (0.23%) had 2 positive reverse transcription PCR tests (55-346 days apart). After eliminating cases with specimens not stored, of suboptimal sequence quality, or belonging to different persons, we obtained valid data from 22 cases. Of those, 4 (0.01%) cases were recurrences involving the same strain; case-patients were 39-93 years of age, and 3 were immunosuppressed. Eighteen (0.04%) cases were reinfections; patients were 19-84 years of age, and most had no relevant clinical history. The second episode was more severe in 8 cases.

SARS-CoV-2 RNA load in nasopharyngeal specimens from outpatients with breakthrough COVID-19 due to Omicron BA.1 and BA.2.

This retrospective observational study compared severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA load in nasopharyngeal specimens (NPs) from patients with breakthrough coronavirus disease 2019 (COVID-19) caused by the Omicron BA.1 or BA.2 sublineages. The convenience sample was composed of 277 outpatients (176 female/112 male; median age, 48 years; range, 12-97) with breakthrough COVID-19 (n = 130 due to BA.1 and n = 147 due to BA.2). All participants had completed a full vaccination schedule and 56% had received a booster vaccine dose at the time of COVID-19 breakthrough microbiological diagnosis. NPs were collected within 7 days (median 2 days) after symptom onset. The TaqPath COVID-19 Combo Kit (Thermo Fisher Scientific) was used to estimate viral loads in NPs. Overall, viral RNA loads in NPs were comparable (p = 0.31) for BA.1 (median, 7.1 log10 copies/ml; range, 2.7-10.6) and BA.2 (median, 7.5 log10 copies/ml; range, 2.7-10.6), yet peak viral load appeared to be reached sooner for BA.2 than for BA.1 (Day 1 vs. Days 3-5; p = 0.002). Time elapsed since last vaccine dose had no significant impact on SARS-CoV-2 RNA loads in the upper respiratory tract (URT) for either BA.1 or BA.2. The data presented do not support that the transmissibility advantage of BA.2 over BA.1 is related to generation of higher viral loads in the URT early after infection.

SARS-CoV-2 Delta-variant breakthrough infections in nursing home residents at midterm after Comirnaty® COVID-19 vaccination.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Delta variant breakthrough infections in nursing home residents following vaccination with Comirnaty® COVID-19 vaccine were characterized. In total, 201 participants (median age, 87 years; range, 64-100; 133 female) from two nursing homes in the Valencian community (Spain) were included. SARS-CoV-2-Spike (S) antibody responses were determined by a lateral flow immunocromatography (LFIC) assay and by quantitative electrochemiluminescent assay in LFIC-negative participants. SARS-CoV-2-S-IFNγ T cells were enumerated by flow cytometry in 10 participants. Nasopharyngeal SARS-CoV-2 RNA loads were quantified by real-time polymerase chain reaction assays. Vaccine breakthrough COVID-19 due to the Delta variant occurred in 39 residents (median age, 87 years; range, 69-96; 31 female) at a median of 6.5 months after vaccination (nine requiring hospitalization). Breakthrough infections occurred at a higher rate (p < 0.0001) in residents who had not been previously infected with SARS-CoV-2 (naïve) (33/108; 18%) than in those with prior diagnosis of SARS-CoV-2 infection (experienced) (6/93; 6.4%), and were more likely (p < 0.0001) to develop in residents who tested negative by LFIC (20/49) at 3 months after vaccination as compared to their LFIC-positive counterparts (19/142). Among LFIC-negative residents, a trend towards lower plasma anti-RBD antibody levels was noticed in those developing breakthrough infection (p = 0.16). SARS-CoV-2 RNA loads in nasopharyngeal specimens were lower in SARS-CoV-2-experienced residents (p < 0.001) and in those testing positive by LFIC (p = 0.13). The frequency of SARS-CoV-2-S-reactive T cells at 3 months was similar in LFIC-negative residents with (n = 7) or without (n = 3) breakthrough infection. Prior history of SARS-CoV-2 infection and detection of S-reactive antibodies by LFIC at 3 months is associated with a lower risk of Delta-variant breakthrough infection in nursing home residents at midterm after Comirnaty® COVID-19 vaccination.

SARS-CoV-2 Omicron BA.1 variant breakthrough infections in nursing home residents after an homologous third dose of the Comirnaty® COVID-19 vaccine: Looking for correlates of protection.

We investigated whether peripheral blood levels of SARS-CoV-2 Spike (S) receptor binding domain antibodies (anti-RBD), neutralizing antibodies (NtAb) targeting Omicron S, and S-reactive-interferon (IFN)-γ-producing CD4+ and CD8+ T cells measured after a homologous booster dose (3D) with the Comirnaty® vaccine was associated with the likelihood of subsequent breakthrough infections due to the Omicron variant. An observational study including 146 nursing home residents (median age, 80 years; range, 66-99; 109 female) evaluated for an immunological response after 3D (at a median of 16 days). Anti-RBD total antibodies were measured by chemiluminescent immunoassay. NtAb were quantified by an Omicron S pseudotyped virus neutralization assay. SARS-CoV-2-S specific-IFNγ-producing CD4+ and CD8+ T cells were enumerated by whole-blood flow cytometry for intracellular cytokine staining. In total, 33/146 participants contracted breakthrough Omicron infection (symptomatic in 30/33) within 4 months after 3D. Anti-RBD antibody levels were comparable in infected and uninfected participants (21 123 vs. 24 723 BAU/ml; p = 0.34). Likewise, NtAb titers (reciprocal IC50 titer, 157 vs. 95; p = 0.32) and frequency of virus-reactive CD4+ (p = 0.82) and CD8+ (p = 0.91) T cells were similar across participants in both groups. anti-RBD antibody levels and NtAb titers estimated at around the time of infection were also comparable (3445 vs. 4345 BAU/ml; p = 0.59 and 188.5 vs. 88.9; p = 0.70, respectively). Having detectable NtAb against Omicron or SARS-CoV-2-S-reactive-IFNγ-producing CD4+ or CD8+ T cells after 3D was not correlated with increased protection from breakthrough infection (OR, 1.50; p = 0.54; OR, 0.0; p = 0.99 and OR 3.70; p = 0.23, respectively). None of the immune parameters evaluated herein, including NtAb titers against the Omicron variant, may reliably predict at the individual level the risk of contracting COVID-19 due to the Omicron variant in nursing home residents.

One is not enough: On the effects of reference genome for the mapping and subsequent analyses of short-reads.

Mapping of high-throughput sequencing (HTS) reads to a single arbitrary reference genome is a frequently used approach in microbial genomics. However, the choice of a reference may represent a source of errors that may affect subsequent analyses such as the detection of single nucleotide polymorphisms (SNPs) and phylogenetic inference. In this work, we evaluated the effect of reference choice on short-read sequence data from five clinically and epidemiologically relevant bacteria (Klebsiella pneumoniae, Legionella pneumophila, Neisseria gonorrhoeae, Pseudomonas aeruginosa and Serratia marcescens). Publicly available whole-genome assemblies encompassing the genomic diversity of these species were selected as reference sequences, and read alignment statistics, SNP calling, recombination rates, dN/dS ratios, and phylogenetic trees were evaluated depending on the mapping reference. The choice of different reference genomes proved to have an impact on almost all the parameters considered in the five species. In addition, these biases had potential epidemiological implications such as including/excluding isolates of particular clades and the estimation of genetic distances. These findings suggest that the single reference approach might introduce systematic errors during mapping that affect subsequent analyses, particularly for data sets with isolates from genetically diverse backgrounds. In any case, exploring the effects of different references on the final conclusions is highly recommended.

Proper Assignation of Reactivation in a COVID-19 Recurrence Initially Interpreted as a Reinfection.

A 77-year-old man (case R) with previous diagnosis of a mild COVID-19 episode was hospitalized 35 days later. On day 23 postadmission, he developed a second COVID-19 episode, now severe, and finally died. Initially, case R's COVID-19 recurrence was interpreted as a reinfection due to the exposure to a SARS-CoV-2 RT-PCR-positive roommate. However, whole-genome sequencing indicated that case R's recurrence corresponded to a reactivation of the strain involved in his first episode. Case R's reactivation had major consequences, leading to a more severe episode, and causing subsequent transmission to another 2 hospitalized patients, 1 of them with fatal outcome.

An unusually high substitution rate in transplant-associated BK polyomavirus in vivo is further concentrated in HLA-C-bound viral peptides.

Infection with human BK polyomavirus, a small double-stranded DNA virus, potentially results in severe complications in immunocompromised patients. Here, we describe the in vivo variability and evolution of the BK polyomavirus by deep sequencing. Our data reveal the highest genomic evolutionary rate described in double-stranded DNA viruses, i.e., 10(-3)-10(-5) substitutions per nucleotide site per year. High mutation rates in viruses allow their escape from immune surveillance and adaptation to new hosts. By combining mutational landscapes across viral genomes with in silico prediction of viral peptides, we demonstrate the presence of significantly more coding substitutions within predicted cognate HLA-C-bound viral peptides than outside. This finding suggests a role for HLA-C in antiviral immunity, perhaps through the action of killer cell immunoglobulin-like receptors. The present study provides a comprehensive view of viral evolution and immune escape in a DNA virus.

Isolation of Four Lytic Phages Infecting Klebsiella pneumoniae K22 Clinical Isolates from Spain.

The emergence of multi-drug-resistant bacteria represents a major public-health threat. Phages constitute a promising alternative to chemical antibiotics due to their high host specificity, abundance in nature, and evolvability. However, phage host specificity means that highly diverse bacterial species are particularly difficult to target for phage therapy. This is the case of Klebsiella pneumoniae, which presents a hypervariable extracellular matrix capsule exhibiting dozens of variants. Here, we report four novel phages infecting K. pneumoniae capsular type K22 which were isolated from environmental samples in Valencia, Spain. Full genome sequencing showed that these phages belong to the Podoviridae family and encode putative depolymerases that allow digestion of specific K22 K. pneumoniae capsules. Our results confirm the capsular type-specificity of K. pneumoniae phages, as indicated by their narrow infectivity in a panel of K. pneumoniae clinical isolates. Nonetheless, this work represents a step forward in the characterization of phage diversity, which may culminate in the future use of large panels of phages for typing and/or for combating multi-drug-resistant K. pneumoniae.

Isolation and Characterization of Two Klebsiella pneumoniae Phages Encoding Divergent Depolymerases.

The emergence of multidrug-resistant bacteria is a major global health concern. The search for new therapies has brought bacteriophages into the spotlight, and new phages are being described as possible therapeutic agents. Among the bacteria that are most extensively resistant to current antibiotics is Klebsiella pneumoniae, whose hypervariable extracellular capsule makes treatment particularly difficult. Here, we describe two new K. pneumoniae phages, πVLC5 and πVLC6, isolated from environmental samples. These phages belong to the genus Drulisvirus within the family Podoviridae. Both phages encode a similar tail spike protein with putative depolymerase activity, which is shared among other related phages and probably determines their ability to specifically infect K. pneumoniae capsular types K22 and K37. In addition, we found that phage πVLC6 also infects capsular type K13 and is capable of striping the capsules of K. pneumoniae KL2 and KL3, although the phage was not infectious in these two strains. Genome sequence analysis suggested that the extended tropism of phage πVLC6 is conferred by a second, divergent depolymerase. Phage πVLC5 encodes yet another putative depolymerase, but we found no activity of this phage against capsular types other than K22 and K37, after testing a panel of 77 reference strains. Overall, our results confirm that most phages productively infected one or few Klebsiella capsular types. This constitutes an important challenge for clinical applications.

Hepatitis C virus intrinsic molecular determinants may contribute to the development of cholestatic hepatitis after liver transplantation.

Cholestatic hepatitis C (CHC) is a severe form of hepatitis C virus (HCV) infection recurrence that leads to high graft loss rates early after liver transplantation (LT). To investigate the pathogenic mechanisms of CHC, we analysed HCV quasispecies in CHC patients compared to a control group (mild hepatitis C recurrence) by deep pyrosequencing. At the time of LT, NS5B quasispecies complexity was similar between the two groups but, after LT, it decreased more sharply in CHC patients than in the control group. Interestingly, the major variant before LT propagated efficiently and remained as the dominant sequence after LT in 62 % of CHC patients versus 11 % of controls (P=0.031). Sequence analysis of the complete non-structural region in a limited number of patients revealed a potential 12 aa signature specific to the CHC group. These data suggest that intrinsic molecular determinants in the circulating HCV quasispecies may provide a fitness advantage, contributing to the development of CHC.

Phylogenomics of Enterococcus faecalis from wild birds: new insights into host-associated differences in core and accessory genomes of the species.

Wild birds have been suggested to be reservoirs of antimicrobial resistant and/or pathogenic Enterococcus faecalis (Efs) strains, but the scarcity of studies and available sequences limit our understanding of the population structure of the species in these hosts. Here, we analysed the clonal and plasmid diversity of 97 Efs isolates from wild migratory birds. We found a high diversity, with most sequence types (STs) being firstly described here, while others were found in other hosts including some predominant in poultry. We found that pheromone-responsive plasmids predominate in wild bird Efs while 35% of the isolates entirely lack plasmids. Then, to better understand the ecology of the species, the whole genome of fivestrains with known STs (ST82, ST170, ST16 and ST55) were sequenced and compared with all the Efs genomes available in public databases. Using several methods to analyse core and accessory genomes (AccNET, PLACNET, hierBAPS and PANINI), we detected differences in the accessory genome of some lineages (e.g. ST82) demonstrating specific associations with birds. Conversely, the genomes of other Efs lineages exhibited divergence in core and accessory genomes, reflecting different adaptive trajectories in various hosts. This pangenome divergence, horizontal gene transfer events and occasional epidemic peaks could explain the population structure of the species.

Hepatitis C virus early kinetics and resistance-associated substitution dynamics during antiviral therapy with direct-acting antivirals.

The emergence of resistance-associated substitutions (RASs) can compromise the high efficacy of direct-acting antivirals (DAAs). Little is known about RASs selection at very early time points during DAA treatment. Therefore, we analyzed the potential emergence of RASs immediately after therapy initiation. Samples of 71 patients treated with different DAAs were collected at baseline, during therapy (hours 4 and 8; days 1-7; weeks 2-4) or until target not detected. HCV-RNA levels were determined by qPCR, and RASs were detected by deep sequencing. Sixty-three (89%) patients achieved a sustained virological response (SVR), 7 (10%) relapsed, and 1 (1%) experienced a breakthrough. Almost all non-SVR (7/8, 88%) showed RASs either at baseline or relapse. High-frequency RASs detected at baseline (Y93H and L159F+C316N) remained detectable at early time points during therapy and reappeared as most prevalent substitutions at relapse. Conversely, emergent RASs at relapse (Q80R, D168E/V, R155K and L31V) were not observed during the first hours-days, before HCV-RNA became undetectable. HCV-RNA decay and genetic evolution of the quasispecies followed a similar pattern during the first hours of therapy in SVR and non-SVR patients. In conclusion, the absence of early RASs selection and the similar dynamics of HCV kinetics and quasispecies in SVR and non-SVR patients after therapy initiation suggest that RASs selection may occur at later stages in the remaining reservoir, where viral populations persist hidden at very low replication levels. Nevertheless, we cannot completely exclude very early selection, when RASs are present below the sensitivity limit of deep sequencing.

Genotypes and subgenotypes of hepatitis B virus circulating in an endemic area in Peru.

Although hepatitis B virus (HBV) infection is still endemic in Abancay, Peru, two decades after vaccination against hepatitis B started in the area, little is known about the diversity and circulation of genotypes and subgenotypes of the virus. To identify the genotypes and subtypes of HBV circulating in Abancay, complete genome sequences of 11 treatment-naive HBV-infected patients were obtained, and phylogenetic analysis was conducted with these and additional sequences from GenBank. Genotyping revealed the presence of genotype F in all the samples from Abancay. Subgenotype F1b was dominant and only one isolate belonged to subgenotype F4, which represents the first description of this subgenotype in Peru. Phylogenetic analysis revealed that most subgenotype F1b isolates from Peru clustered in a subgroup along with two sequences from Argentina, whereas two clusters with two HBV/F1b sequences each were indicative of recent epidemiological linkage, but only one could be verified by independent data. These results suggest that the HBV subgenotype F1b seems to be the predominant subgenotype in Abancay, Peru.

The impact of virus population diversity on the dynamics of cytomegalovirus DNAemia in allogeneic stem cell transplant recipients.

Mixed cytomegalovirus (CMV) infections are associated with delayed viral clearance in solid organ transplant recipients. We investigated whether this could be extrapolated to allogeneic stem cell transplant (allo-SCT) recipients. A total of 48 plasma specimens, obtained during 29 episodes of active CMV infection in 25 non-consecutive allo-SCT patients, were analysed. Baseline blood specimens, drawn shortly prior to the inception of pre-emptive antiviral therapy (pre-treatment specimen; n=29), as well as follow-up samples obtained either after the initiation of antiviral therapy (post-treatment specimen; n=15) or during recurrent episodes (n=4) were analysed. Plasma CMV DNA loads were quantified by real-time PCR and the CMV genotyping was performed by ultra-deep sequencing of hypervariable regions in the genes coding for glycoproteins N (gN) and O (gO). A trend towards higher CMV DNA peak loads, longer CMV DNAemia episode durations and slower CMV DNAemia decay rates was observed for episodes with mixed CMV genotype populations compared to those caused by single CMV variants, although the differences did not reach statistical significance. The length of the treatment course required to clear DNAemia was significantly longer in these mixed episodes (P=0.002). Significant changes in the number or frequency of CMV gN or gO genetic variants were documented following the initiation of antiviral therapy or in recurrent episodes. CMV diversity may have a major impact on the kinetics of CMV DNAemia clearance during the treatment of active CMV infection episodes in allo-SCT recipients.

Comprehensive Screening for Naturally Occurring Hepatitis C Virus Resistance to Direct-Acting Antivirals in the NS3, NS5A, and NS5B Genes in Worldwide Isolates of Viral Genotypes 1 to 6.

There is no comprehensive study available on the natural hepatitis C virus (HCV) polymorphism in sites associated with resistance including all viral genotypes which may present variable susceptibilities to particular direct-acting antivirals (DAAs). This study aimed to analyze the frequencies, genetic barriers, and evolutionary histories of naturally occurring resistance-associated variants (RAVs) in the six main HCV genotypes. A comprehensive analysis of up to 103 RAVs was performed in 2,901, 2,216, and 1,344 HCV isolates for the NS3, NS5A, and NS5B genes, respectively. We report significant intergenotypic differences in the frequencies of natural RAVs for these three HCV genes. In addition, we found a low genetic barrier for the generation of new RAVs, irrespective of the viral genotype. Furthermore, in 1,126 HCV genomes, including sequences spanning the three genes, haplotype analysis revealed a remarkably high frequency of viruses carrying more than one natural RAV to DAAs (53% of HCV-1a, 28.5% of HCV-1b, 67.1% of HCV-6, and 100% of genotype 2, 3, 4, and 5 haplotypes). With the exception of HCV-1a, the most prevalent haplotypes showed RAVs in at least two different viral genes. Finally, evolutionary analyses revealed that, while most natural RAVs appeared recently, others have been efficiently transmitted over time and cluster in well-supported clades. In summary, and despite the observed high efficacy of DAA-based regimens, we show that naturally occurring RAVs are common in all HCV genotypes and that there is an overall low genetic barrier for the selection of resistance mutations. There is a need for natural DAA resistance profiling specific for each HCV genotype.