ABSTRACT
The emergence and spread of new HIV-1 variants pose a challenge for the effectiveness of antiretrovirals (ARV) targeting Pol proteins. During viral evolution, non-synonymous mutations have fixed along the viral genome, leading to amino acid (aa) changes that can be variant-specific (V-markers). Those V-markers fixed in positions associated with drug resistance mutations (DRM), or R-markers, can impact drug susceptibility and resistance pathways. All available HIV-1 Pol sequences from ARV-naïve subjects were downloaded from the United States Los Alamos HIV Sequence Database, selecting 59,733 protease (PR), 6,437 retrotranscriptase (RT), and 6,059 integrase (IN) complete sequences ascribed to the four HIV-1 groups and group M subtypes and circulating recombinant forms (CRFs). Using a bioinformatics tool developed in our laboratory (EpiMolBio), we inferred the consensus sequences for each Pol protein and HIV-1 variant to analyze the aa conservation in Pol. We analyzed the Wu-Kabat protein variability coefficient (WK) in PR, RT, and IN group M to study the susceptibility of each site to evolutionary replacements. We identified as V-markers the variant-specific aa changes present in >75% of the sequences in variants with >5 available sequences, considering R-markers those V-markers that corresponded to DRM according to the IAS-USA2019 and Stanford-Database 9.0. The mean aa conservation of HIV-1 and group M consensus was 82.60%/93.11% in PR, 88.81%/94.07% in RT, and 90.98%/96.02% in IN. The median group M WK was 10 in PR, 4 in RT, and 5 in IN. The residues involved in binding or catalytic sites showed a variability <0.5%. We identified 106 V-markers: 31 in PR, 28 in RT, and 47 in IN, present in 11, 12, and 13 variants, respectively. Among them, eight (7.5%) were R-markers, present in five variants, being minor DRM with little potential effect on ARV susceptibility. We present a thorough analysis of Pol variability among all HIV-1 variants circulating to date. The relatively high aa conservation observed in Pol proteins across HIV-1 variants highlights their critical role in the viral cycle. However, further studies are needed to understand the V-markers' impact on the Pol proteins structure, viral cycle, or treatment strategies, and periodic variability surveillance studies are also required to understand PR, RT, and IN evolution.
ABSTRACT
Monitoring SARS-CoV-2's genetic diversity and emerging mutations in this ongoing pandemic is crucial to understanding its evolution and ensuring the performance of COVID-19 diagnostic tests, vaccines, and therapies. Spain has been one of the main epicenters of COVID-19, reaching the highest number of cases and deaths per 100,000 population in Europe at the beginning of the pandemic. This study aims to investigate the epidemiology of SARS-CoV-2 in Spain and its 18 Autonomous Communities across the six epidemic waves established from February 2020 to January 2022. We report on the circulating SARS-CoV-2 variants in each epidemic wave and Spanish region and analyze the mutation frequency, amino acid (aa) conservation, and most frequent aa changes across each structural/non-structural/accessory viral protein among the Spanish sequences deposited in the GISAID database during the study period. The overall SARS-CoV-2 mutation frequency was 1.24 × 10−5. The aa conservation was >99% in the three types of protein, being non-structural the most conserved. Accessory proteins had more variable positions, while structural proteins presented more aa changes per sequence. Six main lineages spread successfully in Spain from 2020 to 2022. The presented data provide an insight into the SARS-CoV-2 circulation and genetic variability in Spain during the first two years of the pandemic.
Subject(s)
COVID-19 , Pandemics , Amino Acids/genetics , COVID-19/epidemiology , COVID-19/genetics , Genome, Viral , Humans , Mutation , Phylogeny , SARS-CoV-2/genetics , Spain/epidemiologyABSTRACT
HIV envelope transmembrane glycoproteins gp41 (HIV-1) and gp36 (HIV-2) present high variability and play a key role in the HIV-host cell membrane's fusion, as a target for human broadly neutralizing antibodies (bnAbs) and drugs. Thus, a better knowledge of amino acid (aa) conservation across structural domains and HIV variants can help to identify conserved targets to direct new therapeutic and diagnostic strategies. All available gp41/gp36 nucleotide sequences were downloaded from Los Alamos National Laboratory (LANL) HIV Sequence Database, selecting 17,078 sequences ascribed to HIV-1 and HIV-2 variants with ≥3 sequences. After aligning and translating into aa with MEGAv6.0, an in-house bioinformatics program (EpiMolBio) was used to identify the most conserved aa and the aa changes that were specific for each variant (V-markers) vs. HXB2/BEN (HIV-1/HIV-2) reference sequence. We analyzed the presence of specific aa changes among V-markers affecting infectivity, gp41 structure, function, or resistance to the enfuvirtide viral fusion inhibitor (T-20). We also inferred the consensus sequences per HIV variant, describing in each HIV-1 group (M, N, O, P) the conservation level along the complete gp41 per structural domain and locating in each binding site the anti-gp41 human Abs (bnAbs and non bnAbs) described in LANL. We found 38.3/59.7% highly conserved aa present in ≥90% of the 16,803/275 gp41/gp36 sequences ascribed to 105/3 HIV-1/HIV-2 variants, with 9/12.6% of them showing complete conservation across LANL sequences. The fusion peptide, its proximal region, the N-heptad repeat, and the membrane-proximal external region were the gp41 domains with ≥84% of conserved aa in the HIV-1 consensus sequence, the target of most Abs. No natural major resistance mutations to T-20 were observed. Our results show, for the first time, a complete conservation study of gp41/gp36 per variant in the largest panel of HIV variants analyzed to date, providing useful information for a more rational design of drugs, vaccines, and molecular detection tests targeting the HIV transmembrane glycoprotein.
ABSTRACT
The HIV p24 capsid protein has an essential, structural, and functional role in the viral replication cycle, being an interesting target for vaccine design, diagnostic tests, and new antiretroviral drugs (ARVs). The HIV-1 variability poses a challenge for the accuracy and efficiency of diagnostic and treatment tools. This study analyzes p24 diversity among HIV-1 variants and within its secondary structure in HIV-1 M, O, P, and N groups. All available HIV-1 p24 nucleotide sequences were downloaded from the Los Alamos HIV Sequence Database, selecting 23,671 sequences belonging to groups O, N, P, and M (9 subtypes, 7 sub-sub types, and 109 circulating recombinant forms or CRFs). Using a bioinformatics tool developed in our laboratory (EpiMolBio program), we analyzed the amino acid conservation compared to the HXB2 subtype B reference sequence and the V-markers, or amino acid changes that were specific for each variant with at least 10 available sequences. We inferred the p24 consensus sequence for HIV-1 and for each group to analyze the overall conservation in p24 main structural regions, reporting the percentage of substitutions per variant affecting the capsid assembly and molecule-binding, including those associated with resistance to the new capsid-inhibitor lenacapavir, and the key residues involved in lenacapavir-p24 interaction, according to the bibliography. Although the overall structure of p24 was highly conserved, the conservation in the secondary structure varied between HIV-1 variants and the type of secondary structure. All HIV-1 variants presented >80% amino acid conservation vs. HXB2 reference sequence, except for group M sub-subtype F1 (69.27%). Mutants affecting the capsid assembly or lenacapavir capsid-binding were found in <1% of the p24 consensus sequence. Our study reports the HIV-1 variants carrying 14 unique single V-markers in 9/38 group M variants and the level of p24 conservation in each secondary structure region among the 4 HIV-1 groups and group M variants, revealing no natural resistance to lenacapavir in any HIV-1 variant. We present a thorough analysis of p24 variability among all HIV-1 variants circulating to date. Since p24 genetic variability can impact the viral replication cycle and the efficacy of new p24-based diagnostic, therapeutic, and vaccine strategies, conservation studies must consider all HIV-1 variants circulating worldwide.
ABSTRACT
HIV variants carry natural polymorphisms related to drug resistance (R-markers) fixed during viral evolution in the absence of antiretroviral therapy (ART) that may impact on drug susceptibility and resistance pathways. We aimed to identify the HIV type 2 (HIV-2) variant-specific R-markers at Pol in all available sequences from ART-naive subjects deposited in Los Alamos database according to reported HIV-2 drug resistance-associated mutations (DRMs) and report the performance of two online HIV-2 resistance interpretation tools (HIV2EU Tool and Stanford HIVdb Program for HIV-2) to detect them. From a total of 587 sequences, we found 23 R-markers in low frequency, in groups A, B, and G. Four were present in >10% of the sequences with no direct impact on antiretroviral susceptibility. HIV2EU Tool detected one, whereas Stanford program all four. Stanford new tool, although still under development, seems effective in detecting HIV-2 DRMs and may prove a useful tool for HIV-2 resistance interpretation when fully developed.
Subject(s)
Anti-HIV Agents , Drug Resistance, Viral , HIV-2 , Anti-HIV Agents/pharmacology , Drug Resistance, Viral/genetics , Genotype , HIV Infections/drug therapy , HIV-2/drug effects , HIV-2/genetics , Humans , MutationABSTRACT
Monitoring acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genetic diversity and emerging mutations in this ongoing pandemic is crucial for understanding its evolution and assuring the performance of diagnostic tests, vaccines, and therapies against coronavirus disease (COVID-19). This study reports on the amino acid (aa) conservation degree and the global and regional temporal evolution by epidemiological week for each residue of the following four structural SARS-CoV-2 proteins: spike, envelope, membrane, and nucleocapsid. All, 105,276 worldwide SARS-CoV-2 complete and partial sequences from 117 countries available in the Global Initiative on Sharing All Influenza Data (GISAID) from 29 December 2019 to 12 September 2020 were downloaded and processed using an in-house bioinformatics tool. Despite the extremely high conservation of SARS-CoV-2 structural proteins (>99%), all presented aa changes, i.e., 142 aa changes in 65 of the 75 envelope aa, 291 aa changes in 165 of the 222 membrane aa, 890 aa changes in 359 of the 419 nucleocapsid aa, and 2671 changes in 1132 of the 1273 spike aa. Mutations evolution differed across geographic regions and epidemiological weeks (epiweeks). The most prevalent aa changes were D614G (81.5%) in the spike protein, followed by the R203K and G204R combination (37%) in the nucleocapsid protein. The presented data provide insight into the genetic variability of SARS-CoV-2 structural proteins during the pandemic and highlights local and worldwide emerging aa changes of interest for further SARS-CoV-2 structural and functional analysis.