Mojgan Haddad

A systems analysis of mutational effects in HIV-1 protease and reverse transcriptase

The development of a quantitative understanding of viral evolution and the fitness landscape in HIV-1 drug resistance is a formidable challenge given the large number of available drugs and drug resistance mutations. We analyzed a dataset measuring the in vitro fitness of 70,081 virus samples isolated from HIV-1 subtype B infected individuals undergoing routine drug resistance testing. We assayed virus samples for in vitro replicative capacity in the absence of drugs as well as in the presence of 15 individual drugs. We employed a generalized kernel ridge regression to estimate main fitness effects and epistatic interactions of 1,859 single amino acid variants found within the HIV-1 protease and reverse transcriptase sequences. Models including epistatic interactions predict an average of 54.8% of the variance in replicative capacity across the 16 different environments and substantially outperform models based on main fitness effects only. We find that the fitness landscape of HIV-1 protease and reverse transcriptase is characterized by strong epistasis.

Exploring the complexity of the HIV-1 fitness landscape.

Although fitness landscapes are central to evolutionary theory, so far no biologically realistic examples for large-scale fitness landscapes have been described. Most currently available biological examples are restricted to very few loci or alleles and therefore do not capture the high dimensionality characteristic of real fitness landscapes. Here we analyze large-scale fitness landscapes that are based on predictive models for in vitro replicative fitness of HIV-1. We find that these landscapes are characterized by large correlation lengths, considerable neutrality, and high ruggedness and that these properties depend only weakly on whether fitness is measured in the absence or presence of different antiretrovirals. Accordingly, adaptive processes on these landscapes depend sensitively on the initial conditions. While the relative extent to which mutations affect fitness on their own (main effects) or in combination with other mutations (epistasis) is a strong determinant of these properties, the fitness landscape of HIV-1 is considerably less rugged, less neutral, and more correlated than expected from the distribution of main effects and epistatic interactions alone. Overall this study confirms theoretical conjectures about the complexity of biological fitness landscapes and the importance of the high dimensionality of the genetic space in which adaptation takes place.

HER3, p95HER2, and HER2 protein expression levels define multiple subtypes of HER2-positive metastatic breast cancer.

Trastuzumab is effective in the treatment of HER2/neu over-expressing breast cancer, but not all patients benefit from it. In vitro data suggest a role for HER3 in the initiation of signaling activity involving the AKT–mTOR pathway leading to trastuzumab insensitivity. We sought to investigate the potential of HER3 alone and in the context of p95HER2 (p95), a trastuzumab resistance marker, as biomarkers of trastuzumab escape. Using the VeraTag® assay platform, we developed a dual antibody proximity-based assay for the precise quantitation of HER3 total protein (H3T) from formalin-fixed paraffin-embedded (FFPE) breast tumors. We then measured H3T in 89 patients with metastatic breast cancer treated with trastuzumab-based therapy, and correlated the results with progression-free survival and overall survival using Kaplan–Meier and decision tree analyses that also included HER2 total (H2T) and p95 expression levels. Within the sub-population of patients that over-expressed HER2, high levels of HER3 and/or p95 protein expression were significantly associated with poor clinical outcomes on trastuzumab-based therapy. Based on quantitative H3T, p95, and H2T measurements, multiple subtypes of HER2-positive breast cancer were identified that differ in their outcome following trastuzumab therapy. These data suggest that HER3 and p95 are informative biomarkers of clinical outcomes on trastuzumab therapy, and that multiple subtypes of HER2-positive breast cancer may be defined by quantitative measurements of H3T, p95, and H2T.

Assessing predicted HIV-1 replicative capacity in a clinical setting

HIV-1 replicative capacity (RC) provides a measure of within-host fitness and is determined in the context of phenotypic drug resistance testing. However it is unclear how these in-vitro measurements relate to in-vivo processes. Here we assess RCs in a clinical setting by combining a previously published machine-learning tool, which predicts RC values from partial pol sequences with genotypic and clinical data from the Swiss HIV Cohort Study. The machine-learning tool is based on a training set consisting of 65000 RC measurements paired with their corresponding partial pol sequences. We find that predicted RC values (pRCs) correlate significantly with the virus load measured in 2073 infected but drug naïve individuals. Furthermore, we find that, for 53 pairs of sequences, each pair sampled in the same infected individual, the pRC was significantly higher for the sequence sampled later in the infection and that the increase in pRC was also significantly correlated with the increase in plasma viral load and with the length of the time-interval between the sampling points. These findings indicate that selection within a patient favors the evolution of higher replicative capacities and that these in-vitro fitness measures are indicative of in-vivo HIV virus load.

Correlation Between Quantitative HER-2 Protein Expression and Risk for Brain Metastases in HER-2+ Advanced Breast Cancer Patients Receiving Trastuzumab-Containing Therapy

BACKGROUND Patients with human epidermal growth factor receptor (HER)-2+ breast cancer are at particularly high risk for brain metastases; however, the biological basis is not fully understood. Using a novel HER-2 assay, we investigated the correlation between quantitative HER-2 expression in primary breast cancers and the time to brain metastasis (TTBM) in HER-2+ advanced breast cancer patients treated with trastuzumab. METHODS The study group included 142 consecutive patients who were administered trastuzumab-based therapy for HER-2+ metastatic breast cancer. HER-2/neu gene copy number was quantified as the HER-2/centromeric probe for chromosome 17 (CEP17) ratio by central laboratory fluorescence in situ hybridization (FISH). HER-2 protein was quantified as total HER-2 protein expression (H2T) by the HERmark® assay (Monogram Biosciences, Inc., South San Francisco, CA) in formalin-fixed, paraffin-embedded tumor samples. HER-2 variables were correlated with clinical features and TTBM was measured from the initiation of trastuzumab-containing therapy. RESULTS A higher H2T level (continuous variable) was correlated with shorter TTBM, whereas HER-2 amplification by FISH and a continuous HER-2/CEP17 ratio were not predictive (p = .013, .28, and .25, respectively). In the subset of patients that was centrally determined by FISH to be HER-2+, an above-the-median H2T level was significantly associated with a shorter TTBM (hazard ratio, [HR], 2.4; p = .005), whereas this was not true for the median HER-2/CEP17 ratio by FISH (p = .4). Correlation between a continuous H2T level and TTBM was confirmed on multivariate analysis (HR, 3.3; p = .024). CONCLUSIONS These data reveal a strong relationship between the quantitative HER-2 protein expression level and the risk for brain relapse in HER-2+ advanced breast cancer patients. Consequently, quantitative assessment of HER-2 protein expression may inform and facilitate refinements in therapeutic treatment strategies for selected subpopulations of patients in this group.

Correlation of HER2, p95HER2 and HER3 expression and treatment outcome of lapatinib plus capecitabine in her2-positive metastatic breast cancer.

Background Lapatinib plus capecitabine is an effective treatment option for trastuzumab-refractory HER2-positive metastatic breast cancer. We have investigated the correlation between quantitative measures of HER2, p95HER2, and HER3 and treatment outcomes using lapatinib and capecitabine. Methods Total HER2 (H2T), p95HER2 (p95), and total HER3 (H3T) expression were quantified in formalin-fixed paraffin-embedded samples using the VeraTag assays. Patients received lapatinib and capecitabine treatment following trastuzumab failure according to the Lapatinib Expanded Access Program. The association between the protein expression levels and clinical outcomes was analyzed. Results A total of 52 patients were evaluable. H2T level was significantly higher in responders (median 93.49 in partial response, 47.66 in stable disease, and 17.27 in progressive disease; p = 0.020). Longer time-to-progression (TTP) was observed in patients with high H2T [p = 0.018, median 5.2 months in high (>14.95) vs. 1.8 in low (<14.95)] and high H3T [p = 0.017, median 5.0 months in high (>0.605) vs. 2.2 in low (<0.605)]. Patients having both high H2T and high H3T had significantly longer TTP [adjusted hazard ratio (HR) 0.38 (95% CI 0.20–0.73), p = 0.004] and overall survival [adjusted HR 0.46 (95% CI 0.24–0.89), p = 0.020]. No significant association between p95 and response or survival was observed. Conclusions These data suggest a correlation between high HER2 and high HER3 expression and treatment outcome, while no significant difference was observed between clinical outcome and p95 expression level in this cohort of HER2-positive, trastuzumab-refractory metastatic breast cancer patients treated with lapatinib and capecitabine.

Quantitative Measurements of Tumoral p95HER2 Protein Expression in Metastatic Breast Cancer Patients Treated with Trastuzumab: Independent Validation of the p95HER2 Clinical Cutoff

Purpose: P95HER2 (p95) is a truncated form of the HER2, which lacks the trastuzumab-binding site and contains a hyperactive kinase domain. Previously, an optimal clinical cutoff of p95 expression for progression-free survival (PFS) and overall survival (OS) was defined using a quantitative VeraTag assay (Monogram Biosciences) in a training set of trastuzumab-treated metastatic breast cancer (MBC) patients. Experimental Design: In the current study, the predictive value of the p95 VeraTag assay cutoff established in the training set was retrospectively validated for PFS and OS in an independent series of 240 trastuzumab-treated MBC patients from multiple institutions. Results: In the subset of 190 tumors assessed as HER2-total (H2T)-positive using the quantitative HERmark assay (Monogram Biosciences), p95 VeraTag values above the predefined cutoff correlated with shorter PFS (HR = 1.43; P = 0.039) and shorter OS (HR = 1.94; P = 0.0055) where both outcomes were stratified by hormone receptor status and tumor grade. High p95 expression correlated with shorter PFS (HR = 2.41; P = 0.0003) and OS (HR = 2.57; P = 0.0025) in the hormone receptor-positive subgroup of patients (N = 78), but not in the hormone receptor-negative group. In contrast with the quantitative p95 VeraTag measurements, p95 immunohistochemical expression using the same antibody was not significantly correlated with outcomes. Conclusions: The consistency in the p95 VeraTag cutoff across different cohorts of patients with MBC treated with trastuzumab justifies additional studies using blinded analyses in larger series of patients. Clin Cancer Res; 20(10); 2805–13. ©2014 AACR.

Recombination Accelerates Adaptation on a Large-Scale Empirical Fitness Landscape in HIV-1

Recombination has the potential to facilitate adaptation. In spite of the substantial body of theory on the impact of recombination on the evolutionary dynamics of adapting populations, empirical evidence to test these theories is still scarce. We examined the effect of recombination on adaptation on a large-scale empirical fitness landscape in HIV-1 based on in vitro fitness measurements. Our results indicate that recombination substantially increases the rate of adaptation under a wide range of parameter values for population size, mutation rate and recombination rate. The accelerating effect of recombination is stronger for intermediate mutation rates but increases in a monotonic way with the recombination rates and population sizes that we examined. We also found that both fitness effects of individual mutations and epistatic fitness interactions cause recombination to accelerate adaptation. The estimated epistasis in the adapting populations is significantly negative. Our results highlight the importance of recombination in the evolution of HIV-I.

Prevalence in the USA of rilpivirine resistance-associated mutations in clinical samples and effects on phenotypic susceptibility to rilpivirine and etravirine.

BACKGROUND The prevalence of rilpivirine resistance-associated mutations (RAMs) in the USA, and their effect on phenotypic susceptibility to rilpivirine and etravirine, was evaluated in clinical samples from HIV-1-infected patients. METHODS In total, 15,991 samples submitted to Monogram Biosciences (South San Francisco, CA, USA) for routine resistance testing between January 2010 and June 2011 were assessed for the presence of known rilpivirine RAMs K101E/P, E138A/G/K/Q/R, V179L, Y181C/I/V, Y188L, H221Y, F227C and M230I/L; non-nucleoside reverse transcriptase inhibitor (NNRTI) RAMs K103N, L100I and L100I+K103N; and the nucleoside reverse transcriptase inhibitor (NRTI) RAMs M184I/V and their combinations with rilpivirine RAMs. Phenotypic susceptibility (PhenoSenseGT(®) assay; Monogram Biosciences) was evaluated, with reduced susceptibility defined as fold change (FC) in 50% inhibitory concentration (IC50)>2.0 for rilpivirine and FC>2.9 for etravirine. RESULTS Of the 15,991 samples, 17% harboured ≥1 rilpivirine RAMs. The prevalence of most rilpivirine RAMs and combinations of NNRTI RAMs of interest was low (≤3%), except for Y181C (7%). Rilpivirine RAMs were often associated with reduced rilpivirine phenotypic susceptibility. Median FC values >2.0 were observed for clinical isolates with rilpivirine RAMs K101P, E138Q/R, Y181C/I/V, Y188L or M230L, and for the combination of E138K with M184I/V, and K101E with M184I. Most rilpivirine FC values >2.0 were associated with etravirine FC values >2.9 for individual rilpivirine RAMs and those combined with M184I/V. There was no relationship between the presence of K103N and rilpivirine FC. However, the L100I+K103N combination (without rilpivirine RAMs), at <2% prevalence, was associated with a rilpivirine FC>2.0. CONCLUSIONS Based on 15,991 US clinical samples from HIV-1-infected patients, the frequency of most known rilpivirine RAMs apart from Y181C was low.

A decade of HIV-1 drug resistance in the United States: trends and characteristics in a large protease/reverse transcriptase and co-receptor tropism database from 2003 to 2012.

BACKGROUND Drug resistance testing and co-receptor tropism determination are key components of the management of antiretroviral therapy for HIV-1-infected individuals. The purpose of this study was to examine trends of HIV-1 resistance and viral evolution in the past decade by surveying a large commercial patient testing database. METHODS Temporal trends of drug resistance, viral fitness and co-receptor usage among samples submitted for routine phenotypic and genotypic resistance testing to protease inhibitors (PIs), nucleoside reverse transcriptase inhibitors (NRTIs) and non-nucleoside reverse transcriptase inhibitors (NNRTIs), as well as for tropism determination were investigated. RESULTS Within 62,397 resistant viruses reported from 2003 to 2012, we observed a decreasing trend in the prevalence of three-class resistance (from 25% to 9%) driven by decreased resistance to PIs (43% to 21%) and NRTIs (79% to 57%), while observing a slight increase in NNRTI resistance (68% to 75%). The prevalence of CXCR4-mediated entry among tropism testing samples (n=52,945) declined over time from 47% in 2007 to 40% in 2012. A higher proportion of CXCR4-tropic viruses was observed within samples with three-class resistance (50%) compared with the group with no resistance (36%). CONCLUSIONS Decreased prevalence of three-class resistance and increased prevalence of one-class resistance was observed within samples reported between 2003 and 2012. The fraction of CXCR4-tropic viruses has decreased over time; however, CXCR4 usage was more prevalent among multi-class-resistant samples, which may be due to the more advanced disease stage of treatment-experienced patients. These trends have important implications for clinical practice and future drug discovery and development.