CorMut: an R/Bioconductor package for computing correlated mutations based on selection pressure

Natural polymorphisms in HIV-1 CRF01_AE strain and profile of acquired drug resistance mutations in a long-term combination treatment cohort in northeastern China

BACKGROUND

The impacts of genetic polymorphisms on drug resistance mutations (DRMs) among various HIV-1 subtypes have long been debated. In this study, we aimed to analyze the natural polymorphisms and acquired DRM profile in HIV-1 CRF01_AE-infected patients in a large first-line antiretroviral therapy (ART) cohort in northeastern China.

METHODS

The natural polymorphisms of CRF01_AE were analyzed in 2034 patients from a long-term ART cohort in northeastern China. The polymorphisms in 105 treatment failure (TF) patients were compared with those in 1148 treatment success (TS) patients. The acquired DRM profile of 42 patients who experienced TF with tenofovir/lamivudine/efavirenz (TDF/3TC/EFV) treatment was analyzed by comparing the mutations at TF time point to those at baseline. The Stanford HIVdb algorithm was used to interpret the DRMs. Binomial distribution, McNemar test, Wilcoxon test and CorMut package were used to analyze the mutation rates and co-variation. Deep sequencing was used to analyze the evolutionary dynamics of co-variation.

RESULTS

Before ART, there were significantly more natural polymorphisms of 31 sites on reverse transcriptase (RT) in CRF01_AE than subtype B HIV-1 (|Z value| ≥ 3), including five known drug resistance-associated sites (238, 118, 179, 103, and 40). However, only the polymorphism at site 75 was associated with TF (|Z value| ≥ 3). The mutation rate at 14 sites increased significantly at TF time point compared to baseline, with the most common DRMs comprising G190S/C, K65R, K101E/N/Q, M184 V/I, and V179D/I/A/T/E, ranging from 66.7 to 45.2%. Moreover, two unknown mutations (V75 L and L228R) increased by 19.0 and 11.9% respectively, and they were under positive selection (Ka/Ks > 1, log odds ratio [LOD] > 2) and were associated with several other DRMs (cKa/Ks > 1, LOD > 2). Deep sequencing of longitudinal plasma samples showed that L228R occurred simultaneously or followed the appearance of Y181C.

CONCLUSION

The high levels of natural polymorphisms in CRF01_AE had little impact on treatment outcomes. The findings regarding potential new CRF01_AE-specific minor DRMs indicate the need for more studies on the drug resistance phenotype of CRF01_AE.

Chloroplast genome analysis of box-ironbark Eucalyptus

Eucalyptus L'Hérit. (Myrtaceae) is a taxonomically complex and highly speciose genus that dominates much of Australia's woody vegetation. However, very little information is available about the molecular biology and chloroplast diversity of certain groups, such as Eucalyptus section Adnataria, which is found in many woodland habitats of eastern Australia. We report four new complete chloroplast genomes of Eucalyptus, including three genomes from species previously lacking any chloroplast reference sequences. Plastomes of E. albens, E. conica, E. crebra and E. melliodora assembled using a de novo approach were shown to be largely identical to each other, and similar in size and structure to previously published chloroplast genomes from Eucalyptus. A total of 132 genes (114 single-copy genes and 18 duplicated genes in the IR regions) were identified, and shown to be highly conserved in terms of gene order, content and organization. Slightly higher divergence in the intergenic spacers was identified through comparative genomic analyses. Chloroplast sequences of 35 additional individuals representing 12 species were assembled using a reference guided approach. Rates of nucleotide substitution varied among the protein coding genes, with 17 genes under possible positive selection, and 29 invariant genes. Phylogenetic analysis of either the whole reconstructed plastome sequences or the individual genes revealed extreme discordance with expected species boundaries or higher-level relationships. Plastome relationships were better predicted by geography than by nuclear DNA or taxonomic relationships, suggesting a substantial influence of gene flow over and above the effects of incomplete lineage sorting. These results provide resources for future research and valuable insights into the prevalence of interspecific gene flow among Eucalyptus species.

Distinct Campylobacter fetus lineages adapted as livestock pathogens and human pathobionts in the intestinal microbiota

Campylobacter fetus is a venereal pathogen of cattle and sheep, and an opportunistic human pathogen. It is often assumed that C. fetus infection occurs in humans as a zoonosis through food chain transmission. Here we show that mammalian C. fetus consists of distinct evolutionary lineages, primarily associated with either human or bovine hosts. We use whole-genome phylogenetics on 182 strains from 17 countries to provide evidence that C. fetus may have originated in humans around 10,500 years ago and may have "jumped" into cattle during the livestock domestication period. We detect C. fetus genomes in 8% of healthy human fecal metagenomes, where the human-associated lineages are the dominant type (78%). Thus, our work suggests that C. fetus is an unappreciated human intestinal pathobiont likely spread by human to human transmission. This genome-based evolutionary framework will facilitate C. fetus epidemiology research and the development of improved molecular diagnostics and prevention schemes for this neglected pathogen.

Escape from humoral immunity is associated with treatment failure in HIV-1-infected patients receiving long-term antiretroviral therapy

Interindividual heterogeneity in the disease progression of HIV-1-infected patients receiving long-term antiretroviral therapy suggests that some host-related factors may have limited treatment efficacy. To understand the nature of factors contributing to treatment failure, we performed a retrospective cohort study of 45 chronically HIV-1-infected individuals sharing a similar demographics and route of infection, compared the differences between virologically suppressed (VS) and treatment failure (TF) patients with respect to clinical, immunological and virological characteristics. We found that the baseline diversity of HIV-1 env quasispecies was the major difference between VS and TF group, and higher baseline diversity in TF patients. We further predicted TF-related env mutations using a selection pressure-based approach, followed by an analysis of these mutations based on the available three-dimensional structures of gp120/gp41 or their complexes with neutralizing antibodies. Notably, almost all of the identified residues could be mapped to the epitopes of known HIV-1 neutralizing antibodies, especially the epitopes of broadly neutralizing antibodies, and these mutations tended to compromise antibody-antigen interactions. These results indicate that the escape of HIV-1 from host humoral immunity may play a direct role in TF in long-term antiretroviral-experienced patients and that based on env gene sequence of the viruses in the patients.

Emerging Computational Methods for the Rational Discovery of Allosteric Drugs

Allosteric drug development holds promise for delivering medicines that are more selective and less toxic than those that target orthosteric sites. To date, the discovery of allosteric binding sites and lead compounds has been mostly serendipitous, achieved through high-throughput screening. Over the past decade, structural data has become more readily available for larger protein systems and more membrane protein classes (e.g., GPCRs and ion channels), which are common allosteric drug targets. In parallel, improved simulation methods now provide better atomistic understanding of the protein dynamics and cooperative motions that are critical to allosteric mechanisms. As a result of these advances, the field of predictive allosteric drug development is now on the cusp of a new era of rational structure-based computational methods. Here, we review algorithms that predict allosteric sites based on sequence data and molecular dynamics simulations, describe tools that assess the druggability of these pockets, and discuss how Markov state models and topology analyses provide insight into the relationship between protein dynamics and allosteric drug binding. In each section, we first provide an overview of the various method classes before describing relevant algorithms and software packages.

HIV-1 transmitted drug resistance-associated mutations and mutation co-variation in HIV-1 treatment-naïve MSM from 2011 to 2013 in Beijing, China

Background

Transmitted drug resistance (TDR) is an important public health issue, because TDR-associated mutation may affect the outcome of antiretroviral treatment potentially or directly. Men who have sex with men (MSM) constitute a major risk group for HIV transmission. However, current reports are scarce on HIV TDR-associated mutations and their co-variation among MSM.

Methods

Blood samples from 262 newly diagnosed HIV-positive, antiretroviral therapy (ART)-naïve MSM, were collected from January 2011 and December 2013 in Beijing. The polymerase viral genes were sequenced to explore TDR-associated mutations and mutation co-variation.

Results

A total of 223 samples were sequenced and analyzed. Among them, HIV-1 CRF01_AE are accounted for 60.5%, followed by CRF07_BC (27.8%), subtype B (9.9%), and others. Fifty-seven samples had at least one TDR-associated mutation, mainly including L10I/V (6.3%), A71L/T/V (6.3%), V179D/E (5.4%), and V106I (2.7%), with different distributions of TDR-associated mutations by different HIV-1 subtypes and by each year. Moreover, eight significant co-variation pairs were found between TDR-associated mutations (V179D/E) and seven overlapping polymorphisms in subtype CRF01_AE.

Conclusions

To date, this work consists the most comprehensive genetic characterization of HIV-1 TDR-associated mutations prevalent among MSM. It provides important information for understanding TDR and viral evolution among Chinese MSM, a population currently at particularly high risk of HIV transmission.

Electronic supplementary material

The online version of this article (doi:10.1186/s12879-014-0689-7) contains supplementary material, which is available to authorized users.

Single genome analysis reveals genetic characteristics of Neuroadaptation across HIV-1 envelope

Background

The widespread use of highly effective, combination antiretroviral therapy (cART) has led to a significant reduction in the incidence of HIV-associated dementia (HAD). Despite these advances, the prevalence of HIV-1 associated neurocognitive disorders (HANDs) has been estimated at approximately 40%-50%. In the cART era, the majority of this disease burden is represented by asymptomatic neurocognitive impairment and mild neurocognitive disorder (ANI and MND respectively). Although less severe than HAD, these diagnoses carry with them substantial morbidity.

Results

In this cross-sectional study, single genome amplification (SGA) was used to sequence 717 full-length HIV-1 envelope (env) clade B variants from the paired cerebrospinal fluid (CSF) and blood plasma samples of fifteen chronically infected HIV-positive individuals with normal neurocognitive performance (NCN), ANI and MND. Various degrees of compartmentalization were found across disease states and history of cART utilization. In individuals with compartmentalized virus, mean HIV-1 env population diversity was lower in the CSF than plasma-derived variants. Overall, mean V1V2 loop length was shorter in CSF-derived quasispecies when compared to contemporaneous plasma populations, and this was found to correlate with a lower mean number of N-linked glycosylation sites in this region. A number of discrete amino acid positions that correlate strongly with compartmentalization in the CSF were identified in both variable and constant regions of gp120 as well as in gp41. Correlated mutation analyses further identified that a subset of amino acid residues in these compartmentalization “hot spot” positions were strongly correlated with one another, suggesting they may play an important, definable role in the adaptation of viral variants to the CSF. Analysis of these hot spots in the context of a well-supported crystal structure of HIV-1 gp120 suggests mechanisms through which amino acid differences at the identified residues might contribute to viral compartmentalization in the CSF.

Conclusions

The detailed analyses of SGA-derived full length HIV-1 env from subjects with both normal neurocognitive performance and the most common HAND diagnoses in the cART era allow us to identify novel and confirm previously described HIV-1 env genetic determinants of neuroadaptation and relate potential motifs to HIV-1 env structure and function.

Electronic supplementary material

The online version of this article (doi:10.1186/s12977-014-0065-0) contains supplementary material, which is available to authorized users.