Stability analysis of mixtures of mutagenetic trees

<p>Abstract</p> <p>Background</p> <p>Mixture models of mutagenetic trees are evolutionary models that capture several pathways of ordered accumulation of genetic events observed in different subsets of patients. They were used to model HIV progression by accumulation of resistance mutations in the viral genome under drug pressure and cancer progression by accumulation of chromosomal aberrations in tumor cells. From the mixture models a genetic progression score (GPS) can be derived that estimates the genetic status of single patients according to the corresponding progression along the tree models. GPS values were shown to have predictive power for estimating drug resistance in HIV or the survival time in cancer. Still, the reliability of the exact values of such complex markers derived from graphical models can be questioned.</p> <p>Results</p> <p>In a simulation study, we analyzed various aspects of the stability of estimated mutagenetic trees mixture models. It turned out that the induced probabilistic distributions and the tree topologies are recovered with high precision by an EM-like learning algorithm. However, only for models with just one major model component, also GPS values of single patients can be reliably estimated.</p> <p>Conclusion</p> <p>It is encouraging that the estimation process of mutagenetic trees mixture models can be performed with high confidence regarding induced probability distributions and the general shape of the tree topologies. For a model with only one major disease progression process, even genetic progression scores for single patients can be reliably estimated. However, for models with more than one relevant component, alternative measures should be introduced for estimating the stage of disease progression.</p

Mapping HIV-1 Vaccine Induced T-Cell Responses: Bias towards Less-Conserved Regions and Potential Impact on Vaccine Efficacy in the Step Study

T cell directed HIV vaccines are based upon the induction of CD8+ T cell memory responses that would be effective in inhibiting infection and subsequent replication of an infecting HIV-1 strain, a process that requires a match or near-match between the epitope induced by vaccination and the infecting viral strain. We compared the frequency and specificity of the CTL epitope responses elicited by the replication-defective Ad5 gag/pol/nef vaccine used in the Step trial with the likelihood of encountering those epitopes among recently sequenced Clade B isolates of HIV-1. Among vaccinees with detectable 15-mer peptide pool ELISpot responses, there was a median of four (one Gag, one Nef and two Pol) CD8 epitopes per vaccinee detected by 9-mer peptide ELISpot assay. Importantly, frequency analysis of the mapped epitopes indicated that there was a significant skewing of the T cell response; variable epitopes were detected more frequently than would be expected from an unbiased sampling of the vaccine sequences. Correspondingly, the most highly conserved epitopes in Gag, Pol, and Nef (defined by presence in >80% of sequences currently in the Los Alamos database www.hiv.lanl.gov) were detected at a lower frequency than unbiased sampling, similar to the frequency reported for responses to natural infection, suggesting potential epitope masking of these responses. This may be a generic mechanism used by the virus in both contexts to escape effective T cell immune surveillance. The disappointing results of the Step trial raise the bar for future HIV vaccine candidates. This report highlights the bias towards less-conserved epitopes present in the same vaccine used in the Step trial. Development of vaccine strategies that can elicit a greater breadth of responses, and towards conserved regions of the genome in particular, are critical requirements for effective T-cell based vaccines against HIV-1

Kinetic Pathway of Pyrophosphorolysis by a Retrotransposon Reverse Transcriptase

DNA and RNA polymerases use a common phosphoryl transfer mechanism for base addition that requires two or three acidic amino acid residues at their active sites. We previously showed, for the reverse transcriptase (RT) encoded by the yeast retrotransposon Ty1, that one of the three conserved active site aspartates (D211) can be substituted by asparagine and still retain in vitro polymerase activity, although in vivo transposition is lost. Transposition is partially restored by second site suppressor mutations in the RNAse H domain. The novel properties of this amino acid substitution led us to express the WT and D211N mutant enzymes, and study their pre-steady state kinetic parameters. We found that the kpol was reduced by a factor of 223 in the mutant, although the Kd for nucleotide binding was unaltered. Further, the mutant enzyme had a marked preference for Mn2+ over Mg2+. To better understand the functions of this residue within the Ty1 RT active site, we have now examined the in vitro properties of WT and D211N mutant Ty1 RTs in carrying out pyrophosphorolysis, the reverse reaction to polymerization, where pyrophosphate is the substrate and dNTPs are the product. We find that pyrophosphorolysis is efficient only when the base-paired primer template region is >14 bases, and that activity increases when the primer end is blunt-ended or recessed by only a few bases. Using pre-steady state kinetic analysis, we find that the rate of pyrophosphorolysis (kpyro) in the D211N mutant is nearly 320 fold lower than the WT enzyme, and that the mutant enzyme has an ∼170 fold lower apparent Kd for pyrophosphate. These findings indicate that subtle substrate differences can strongly affect the enzyme's ability to properly position the primer-end to carry out pyrophosphorolysis. Further the kinetic data suggests that the D211 residue has a role in pyrophosphate binding and release, which could affect polymerase translocation, and help explain the D211N mutant's transposition defect

Cross-validated stepwise regression for identification of novel non-nucleoside reverse transcriptase inhibitor resistance associated mutations

<p>Abstract</p> <p>Background</p> <p>Linear regression models are used to quantitatively predict drug resistance, the phenotype, from the HIV-1 viral genotype. As new antiretroviral drugs become available, new resistance pathways emerge and the number of resistance associated mutations continues to increase. To accurately identify which drug options are left, the main goal of the modeling has been to maximize predictivity and not interpretability. However, we originally selected linear regression as the preferred method for its transparency as opposed to other techniques such as neural networks. Here, we apply a method to lower the complexity of these phenotype prediction models using a 3-fold cross-validated selection of mutations.</p> <p>Results</p> <p>Compared to standard stepwise regression we were able to reduce the number of mutations in the reverse transcriptase (RT) inhibitor models as well as the number of interaction terms accounting for synergistic and antagonistic effects. This reduction in complexity was most significant for the non-nucleoside reverse transcriptase inhibitor (NNRTI) models, while maintaining prediction accuracy and retaining virtually all known resistance associated mutations as first order terms in the models. Furthermore, for etravirine (ETR) a better performance was seen on two years of unseen data. By analyzing the phenotype prediction models we identified a list of forty novel NNRTI mutations, putatively associated with resistance. The resistance association of novel variants at known NNRTI resistance positions: 100, 101, 181, 190, 221 and of mutations at positions not previously linked with NNRTI resistance: 102, 139, 219, 241, 376 and 382 was confirmed by phenotyping site-directed mutants.</p> <p>Conclusions</p> <p>We successfully identified and validated novel NNRTI resistance associated mutations by developing parsimonious resistance prediction models in which repeated cross-validation within the stepwise regression was applied. Our model selection technique is computationally feasible for large data sets and provides an approach to the continued identification of resistance-causing mutations.</p

Adaptation of HIV-1 Depends on the Host-Cell Environment

Many viruses have the ability to rapidly develop resistance against antiviral drugs and escape from the host immune system. To which extent the host environment affects this adaptive potential of viruses is largely unknown. Here we show that for HIV-1, the host-cell environment is key to the adaptive potential of the virus. We performed a large-scale selection experiment with two HIV-1 strains in two different T-cell lines (MT4 and C8166). Over 110 days of culture, both virus strains adapted rapidly to the MT4 T-cell line. In contrast, when cultured on the C8166 T-cell line, the same strains did not show any increase in fitness. By sequence analyses and infections with viruses expressing either yellow or cyan fluorescent protein, we were able to show that the absence of adaptation was linked to a lower recombination rate in the C8166 T-cell line. Our findings suggest that if we can manipulate the host-cellular factors that mediate viral evolution, we may be able to significantly retard viral adaptability

Deciphering the Complex Distribution of Human Immunodeficiency Virus Type 1 Subtypes among Different Cohorts in Northern Tanzania.

Increased understanding of the genetic diversity of HIV-1 is challenging but important in the development of an effective vaccine. We aimed to describe the distribution of HIV-1 subtypes in northern Tanzania among women enrolled in studies preparing for HIV-1 prevention trials (hospitality facility-worker cohorts), and among men and women in an open cohort demographic surveillance system (Kisesa cohort). The polymerase encompassing partial reverse transcriptase was sequenced and phylogenetic analysis performed and subtype determined. Questionnaires documented demographic data. We examined factors associated with subtype using multinomial logistic regression, adjusted for study, age, and sex. Among 140 individuals (125 women and 15 men), subtype A1 predominated (54, 39%), followed by C (46, 33%), D (25, 18%) and unique recombinant forms (URFs) (15, 11%). There was weak evidence to suggest different subtype frequencies by study (for example, 18% URFs in the Kisesa cohort versus 5-9% in the hospitality facility-worker cohorts; adjusted relative-risk ratio (aRR) = 2.35 [95% CI 0.59,9.32]; global p = 0.09). Compared to men, women were less likely to have subtype D versus A (aRR = 0.12 [95% CI 0.02,0.76]; global p = 0.05). There was a trend to suggest lower relative risk of subtype D compared to A with older age (aRR = 0.44 [95% CI 0.23,0.85] per 10 years; global p = 0.05). We observed multiple subtypes, confirming the complex genetic diversity of HIV-1 strains circulating in northern Tanzania, and found some differences between cohorts and by age and sex. This has important implications for vaccine design and development, providing opportunity to determine vaccine efficacy in diverse HIV-1 strains

Mechanisms of relapse in acute leukaemia: involvement of p53 mutated subclones in disease progression in acute lymphoblastic leukaemia

Mutations of the p53 tumour suppressor gene are infrequent at presentation of both acute myeloblastic leukaemia (AML) and acute lymphoblastic leukaemia (ALL), being found in between 5–10% of AML and 2–3% of ALL. Here we have studied the frequency of detection of p53 mutations at relapse of both AML and B-precursor ALL. In those patients with detectable mutations at relapse we investigated whether the mutation was detectable at presentation and was thus an early initiating event or whether it had arisen as a late event associated with relapse. Bone marrow samples from 55 adults and children with relapsed AML (n = 41) or ALL (n = 14) were analysed for p53 gene alterations by direct sequencing of exons 5–9. For samples where a p53 mutation was found at relapse, analysis of presentation samples was carried out by direct sequencing of the exon involved, or by allele-specific polymerase chain reaction (PCR) if the mutation could not be detected using direct sequencing. A p53 mutated gene was found at relapse in seven out of 55 cases. The frequency was higher in relapsed ALL (four out of 14 cases; 28.6%) compared to AML (three out of 41 cases; 7.3%). In five out of the seven cases presentation samples were available to study for the presence of the mutation. In two out of two AML patients the p53 mutation was detectable in the presentation sample by direct sequencing. In three ALL patients analysis of presentation material by direct sequencing showed a small mutant peak in one case, the other two being negative despite the sample analysed containing > 90% blast cells. However in both of these patients, the presence of p53 mutation was confirmed in the presentation sample using allele-specific PCR. In one of these patients the emergence of a subclone at relapse was confirmed by clonality analysis using IgH fingerprinting. Our results confirm that in ALL p53 mutations are present in a proportion of patients at relapse. Furthermore cells carrying the mutation are detectable at presentation in a minor clone suggesting that p53 mutations in ALL may be a mechanism contributing to disease relapse. © 1999 Cancer Research Campaig