Guangdi Li

A new ensemble coevolution system for detecting HIV-1 protein coevolution

BackgroundA key challenge in the field of HIV-1 protein evolution is the identification of coevolving amino acids at the molecular level. In the past decades, many sequence-based methods have been designed to detect position-specific coevolution within and between different proteins. However, an ensemble coevolution system that integrates different methods to improve the detection of HIV-1 protein coevolution has not been developed.ResultsWe integrated 27 sequence-based prediction methods published between 2004 and 2013 into an ensemble coevolution system. This system allowed combinations of different sequence-based methods for coevolution predictions. Using HIV-1 protein structures and experimental data, we evaluated the performance of individual and combined sequence-based methods in the prediction of HIV-1 intra- and inter-protein coevolution. We showed that sequence-based methods clustered according to their methodology, and a combination of four methods outperformed any of the 27 individual methods. This four-method combination estimated that HIV-1 intra-protein coevolving positions were mainly located in functional domains and physically contacted with each other in the protein tertiary structures. In the analysis of HIV-1 inter-protein coevolving positions between Gag and protease, protease drug resistance positions near the active site mostly coevolved with Gag cleavage positions (V128, S373-T375, A431, F448-P453) and Gag C-terminal positions (S489-Q500) under selective pressure of protease inhibitors.ConclusionsThis study presents a new ensemble coevolution system which detects position-specific coevolution using combinations of 27 different sequence-based methods. Our findings highlight key coevolving residues within HIV-1 structural proteins and between Gag and protease, shedding light on HIV-1 intra- and inter-protein coevolution.ReviewersThis article was reviewed by Dr. Zoltán Gáspári.

High-Throughput Analysis of Human Cytomegalovirus Genome Diversity Highlights the Widespread Occurrence of Gene-Disrupting Mutations and Pervasive Recombination

ABSTRACT Human cytomegalovirus is a widespread pathogen of major medical importance. It causes significant morbidity and mortality in immunocompromised individuals, and congenital infections can result in severe disabilities or stillbirth. Development of a vaccine is prioritized, but no candidate is close to release. Although correlations of viral genetic variability with pathogenicity are suspected, knowledge about the strain diversity of the 235-kb genome is still limited. In this study, 96 full-length human cytomegalovirus genomes from clinical isolates were characterized, quadrupling the amount of information available for full-genome analysis. These data provide the first high-resolution map of human cytomegalovirus interhost diversity and evolution. We show that cytomegalovirus is significantly more divergent than all other human herpesviruses and highlight hot spots of diversity in the genome. Importantly, 75% of strains are not genetically intact but contain disruptive mutations in a diverse set of 26 genes, including the immunomodulatory genes UL40 and UL111A. These mutants are independent of culture passage artifacts and circulate in natural populations. Pervasive recombination, which is linked to the widespread occurrence of multiple infections, was found throughout the genome. The recombination density was significantly higher than those of other human herpesviruses and correlated with strain diversity. While the overall effects of strong purifying selection on virus evolution are apparent, evidence of diversifying selection was found in several genes encoding proteins that interact with the host immune system, including UL18, UL40, UL142, and UL147. These residues may present phylogenetic signatures of past and ongoing virus-host interactions. IMPORTANCE Human cytomegalovirus has the largest genome of all viruses that infect humans. Currently, there is a great interest in establishing associations between genetic variants and strain pathogenicity of this herpesvirus. Since the number of publicly available full-genome sequences is limited, knowledge about strain diversity is highly fragmented and biased toward a small set of loci. Combined with our previous work, we have now contributed 101 complete genome sequences. We have used these data to conduct the first high-resolution analysis of interhost genome diversity, providing an unbiased and comprehensive overview of cytomegalovirus variability. These data are of major value to the development of novel antivirals and a vaccine and to identify potential targets for genotype-phenotype experiments. Furthermore, these data have enabled a thorough study of the evolutionary processes that have shaped cytomegalovirus diversity.

An integrated map of HIV genome-wide variation from a population perspective

BackgroundThe HIV pandemic is characterized by extensive genetic variability, which has challenged the development of HIV drugs and vaccines. Although HIV genomes have been classified into different types, groups, subtypes and recombinants, a comprehensive study that maps HIV genome-wide diversity at the population level is still lacking to date. This study aims to characterize HIV genomic diversity in large-scale sequence populations, and to identify driving factors that shape HIV genome diversity.ResultsA total of 2996 full-length genomic sequences from 1705 patients infected with 16 major HIV groups, subtypes and circulating recombinant forms (CRFs) were analyzed along with structural, immunological and peptide inhibitor information. Average nucleotide diversity of HIV genomes was almost 50% between HIV-1 and HIV-2 types, 37.5% between HIV-1 groups, 14.7% between HIV-1 subtypes, 8.2% within individual HIV-1 subtypes and less than 1% within single patients. Along the HIV genome, diversity patterns and compositions of nucleotides and amino acids were highly similar across different groups, subtypes and CRFs. Current HIV-derived peptide inhibitors were predominantly derived from conserved, solvent accessible and intrinsically ordered structures in the HIV-1 subtype B genome. We identified these conserved regions in Capsid, Nucleocapsid, Protease, Integrase, Reverse transcriptase, Vpr and the GP41 N terminus as potential drug targets. In the analysis of factors that impact HIV-1 genomic diversity, we focused on protein multimerization, immunological constraints and HIV-human protein interactions. We found that amino acid diversity in monomeric proteins was higher than in multimeric proteins, and diversified positions were preferably located within human CD4 T cell and antibody epitopes. Moreover, intrinsic disorder regions in HIV-1 proteins coincided with high levels of amino acid diversity, facilitating a large number of interactions between HIV-1 and human proteins.ConclusionsThis first large-scale analysis provided a detailed mapping of HIV genomic diversity and highlighted drug-target regions conserved across different groups, subtypes and CRFs. Our findings suggest that, in addition to the impact of protein multimerization and immune selective pressure on HIV-1 diversity, HIV-human protein interactions are facilitated by high variability within intrinsically disordered structures.

Current therapy for chronic hepatitis C: The role of direct-acting antivirals

Abstract One of the most exciting developments in antiviral research has been the discovery of the direct‐acting antivirals (DAAs) that effectively cure chronic hepatitis C virus (HCV) infections. Based on more than 100 clinical trials and real‐world studies, we provide a comprehensive overview of FDA‐approved therapies and newly discovered anti‐HCV agents with a special focus on drug efficacy, mechanisms of action, and safety. We show that HCV drug development has advanced in multiple aspects: (i) interferon‐based regimens were replaced by interferon‐free regimens; (ii) genotype‐specific drugs evolved to drugs for all HCV genotypes; (iii) therapies based upon multiple pills per day were simplified to a single pill per day; (iv) drug potency increased from moderate (˜60%) to high (>90%) levels of sustained virologic responses; (v) treatment durations were shortened from 48 to 12 or 8 weeks; and (vi) therapies could be administered orally regardless of prior treatment history and cirrhotic status. However, despite these remarkable achievements made in HCV drug discovery, challenges remain in the management of difficult‐to‐treat patients. HighlightsHCV genotype‐specific drugs evolve to pan‐genotypic drugs.Drug potency increases from moderate (˜60%) to high (>90%) levels of sustained virologic response.Treatment durations are shortened from a 48‐week to 12‐week or 8‐week period.HCV therapies based upon multiple pills per day are simplified to a single pill per day.HCV therapies are administered orally regardless of prior treatment history and cirrhotic status.

Genetic Diversity and Selective Pressure in Hepatitis C Virus Genotypes 1–6: Significance for Direct-Acting Antiviral Treatment and Drug Resistance

Treatment with pan-genotypic direct-acting antivirals, targeting different viral proteins, is the best option for clearing hepatitis C virus (HCV) infection in chronically infected patients. However, the diversity of the HCV genome is a major obstacle for the development of antiviral drugs, vaccines, and genotyping assays. In this large-scale analysis, genome-wide diversity and selective pressure was mapped, focusing on positions important for treatment, drug resistance, and resistance testing. A dataset of 1415 full-genome sequences, including genotypes 1–6 from the Los Alamos database, was analyzed. In 44% of all full-genome positions, the consensus amino acid was different for at least one genotype. Focusing on positions sharing the same consensus amino acid in all genotypes revealed that only 15% was defined as pan-genotypic highly conserved (≥99% amino acid identity) and an additional 24% as pan-genotypic conserved (≥95%). Despite its large genetic diversity, across all genotypes, codon positions were rarely identified to be positively selected (0.23%–0.46%) and predominantly found to be under negative selective pressure, suggesting mainly neutral evolution. For NS3, NS5A, and NS5B, respectively, 40% (6/15), 33% (3/9), and 14% (2/14) of the resistance-related positions harbored as consensus the amino acid variant related to resistance, potentially impeding treatment. For example, the NS3 variant 80K, conferring resistance to simeprevir used for treatment of HCV1 infected patients, was present in 39.3% of the HCV1a strains and 0.25% of HCV1b strains. Both NS5A variants 28M and 30S, known to be associated with resistance to the pan-genotypic drug daclatasvir, were found in a significant proportion of HCV4 strains (10.7%). NS5B variant 556G, known to confer resistance to non-nucleoside inhibitor dasabuvir, was observed in 8.4% of the HCV1b strains. Given the large HCV genetic diversity, sequencing efforts for resistance testing purposes may need to be genotype-specific or geographically tailored.

Functional conservation of HIV-1 Gag: implications for rational drug design.

BackgroundHIV-1 replication can be successfully blocked by targeting gag gene products, offering a promising strategy for new drug classes that complement current HIV-1 treatment options. However, naturally occurring polymorphisms at drug binding sites can severely compromise HIV-1 susceptibility to gag inhibitors in clinical and experimental studies. Therefore, a comprehensive understanding of gag natural diversity is needed.FindingsWe analyzed the degree of functional conservation in 10862 full-length gag sequences across 8 major HIV-1 subtypes and identified the impact of natural variation on known drug binding positions targeted by more than 20 gag inhibitors published to date. Complete conservation across all subtypes was detected in 147 (29%) out of 500 gag positions, with the highest level of conservation observed in capsid protein. Almost half (41%) of the 136 known drug binding positions were completely conserved, but all inhibitors were confronted with naturally occurring polymorphisms in their binding sites, some of which correlated with HIV-1 subtype. Integration of sequence and structural information revealed one drug binding pocket with minimal genetic variability, which is situated at the N-terminal domain of the capsid protein.ConclusionsThis first large-scale analysis of full-length HIV-1 gag provided a detailed mapping of natural diversity across major subtypes and highlighted the considerable variation in current drug binding sites. Our results contribute to the optimization of gag inhibitors in rational drug design, given that drug binding sites should ideally be conserved across all HIV-1 subtypes.

Trends and Predictors of Transmitted Drug Resistance (TDR) and Clusters with TDR in a Local Belgian HIV-1 Epidemic

We aimed to study epidemic trends and predictors for transmitted drug resistance (TDR) in our region, its clinical impact and its association with transmission clusters. We included 778 patients from the AIDS Reference Center in Leuven (Belgium) diagnosed from 1998 to 2012. Resistance testing was performed using population-based sequencing and TDR was estimated using the WHO-2009 surveillance list. Phylogenetic analysis was performed using maximum likelihood and Bayesian techniques. The cohort was predominantly Belgian (58.4%), men who have sex with men (MSM) (42.8%), and chronically infected (86.5%). The overall TDR prevalence was 9.6% (95% confidence interval (CI): 7.7–11.9), 6.5% (CI: 5.0–8.5) for nucleoside reverse transcriptase inhibitors (NRTI), 2.2% (CI: 1.4–3.5) for non-NRTI (NNRTI), and 2.2% (CI: 1.4–3.5) for protease inhibitors. A significant parabolic trend of NNRTI-TDR was found (p = 0.019). Factors significantly associated with TDR in univariate analysis were male gender, Belgian origin, MSM, recent infection, transmission clusters and subtype B, while multivariate and Bayesian network analysis singled out subtype B as the most predictive factor of TDR. Subtype B was related with transmission clusters with TDR that included 42.6% of the TDR patients. Thanks to resistance testing, 83% of the patients with TDR who started therapy had undetectable viral load whereas half of the patients would likely have received a suboptimal therapy without this test. In conclusion, TDR remained stable and a NNRTI up-and-down trend was observed. While the presence of clusters with TDR is worrying, we could not identify an independent, non-sequence based predictor for TDR or transmission clusters with TDR that could help with guidelines or public health measures.

Impact of HCV genotype on treatment regimens and drug resistance: a snapshot in time.

The introduction of highly potent direct‐acting antivirals (DAAs) has revolutionized hepatitis C virus treatment. Nevertheless, viral eradication worldwide remains a challenge also in the era of DAA treatment, because of the high associated costs, high numbers of undiagnosed patients, high re‐infection rates in some risk groups and suboptimal drug efficacies associated with host and viral factors as well as advanced stages of liver disease. A correct determination of the HCV genotype allows administration of the most appropriate antiviral regimen. Additionally, HCV genetic sequencing improves our understanding of resistance‐associated variants, either naturally occurring before treatment, acquired by transmission at HCV infection, or emerging after virological failure. Because treatment response rates, and the prevalence and development of drug resistance variants differ for each DAA regimen and HCV genotype, this review summarizes treatment opportunities per HCV genotype, and focuses on viral genetic sequencing to guide clinical decision making. Although approval of the first pan‐genotypic DAA‐only regimen is expected soon, HCV genetic sequencing will remain important because when DAA therapies fail, genotyping and resistance testing to select a new active DAA combination will be essential. Copyright

A stably expressed llama single-domain intrabody targeting Rev displays broad-spectrum anti-HIV activity

The HIV Rev protein mediates the transport of partially and unspliced HIV mRNA from the nucleus to the cytoplasm. Rev multimerizes on a secondary stem-loop structure present in the viral intron-containing mRNA species and recruits the cellular karyopherin CRM1 to export viral mRNAs from the nucleus to the cytoplasm. Previously we have identified a single-domain intrabody (Nb(190)), derived from a llama heavy-chain antibody, which efficiently inhibits Rev multimerization and suppresses the production of infectious virus. We recently mapped the epitope of this nanobody and demonstrated that Rev residues K20 and Y23 are crucial for interaction while residues V16, H53 and L60 are important to a lesser extent. Here, we generated cell lines stably expressing Nb(190) and assessed the capacity of these cell lines to suppress the replication of different HIV-1 subtypes. These cells stably expressing the single-domain antibody are protected from virus-induced cytopathogenic effect even in the context of high multiplicity of infection. In addition, the replication of different subtypes of group M and one strain of group O is significantly suppressed in these cell lines. Next, we analysed the natural variations of Rev amino acids in sequence samples from HIV-1 infected patients worldwide and assessed the effect of Nb(190) on the most prevalent polymorphisms occurring at the key epitope positions (K20 and Y23) in Rev. We found that Nb(190) was able to suppress the function of these Rev variants except for the K20N mutant, which was present in only 0.7% of HIV-1 sequence populations (n = 4632). Cells stably expressing the single-domain intrabody Nb(190) are protected against virus-induced cytopathogenic effect and display a selective survival advantage upon infection. In addition, Nb(190) suppresses the replication of a wide range of different HIV-1 subtypes. Large-scale sequence analysis reveals that the Nb(190) epitope positions in Rev are well conserved across major HIV-1 subtypes and groups. Altogether, our results indicate that Nb(190) may have broad potential as a gene therapeutic agent against HIV-1.

HIV-1 Gag C-terminal amino acid substitutions emerging under selective pressure of protease inhibitors in patient populations infected with different HIV-1 subtypes

HIV-1 Gag amino acid substitutions associated with protease inhibitor (PI) treatment have mainly been reported in subtype B, while information on other subtypes is scarce. Using sequences from 11613 patients infected with different HIV-1 subtypes, we evaluated the prevalence of 93 Gag amino acid substitutions and their association with genotypic PI resistance. A significant association was found for 13 Gag substitutions, including A431V in both subtype B and CRF01_AE. K415R in subtype C and S451G in subtype B were newly identified. Most PI-associated Gag substitutions are located in the flexible C-terminal domain, revealing the key role this region plays in PI resistance.