Low-Cost Ultra-Wide Genotyping Using Roche/454 Pyrosequencing for Surveillance of HIV Drug Resistance

Background: Great efforts have been made to increase accessibility of HIV antiretroviral therapy (ART) in low and middle-income countries. The threat of wide-scale emergence of drug resistance could severely hamper ART scale-up efforts. Population-based surveillance of transmitted HIV drug resistance ensures the use of appropriate first-line regimens to maximize efficacy of ART programs where drug options are limited. However, traditional HIV genotyping is extremely expensive, providing a cost barrier to wide-scale and frequent HIV drug resistance surveillance. Methods/Results: We have developed a low-cost laboratory-scale next-generation sequencing-based genotyping method to monitor drug resistance. We designed primers specifically to amplify protease and reverse transcriptase from Brazilian HIV subtypes and developed a multiplexing scheme using multiplex identifier tags to minimize cost while providing more robust data than traditional genotyping techniques. Using this approach, we characterized drug resistance from plasma in 81 HIV infected individuals collected in Sao Paulo, Brazil. We describe the complexities of analyzing next-generation sequencing data and present a simplified open-source workflow to analyze drug resistance data. From this data, we identified drug resistance mutations in 20% of treatment naive individuals in our cohort, which is similar to frequencies identified using traditional genotyping in Brazilian patient samples. Conclusion: The developed ultra-wide sequencing approach described here allows multiplexing of at least 48 patient samples per sequencing run, 4 times more than the current genotyping method. This method is also 4-fold more sensitive (5% minimal detection frequency vs. 20%) at a cost 3-5 x less than the traditional Sanger-based genotyping method. Lastly, by using a benchtop next-generation sequencer (Roche/454 GS Junior), this approach can be more easily implemented in low-resource settings. This data provides proof-of-concept that next-generation HIV drug resistance genotyping is a feasible and low-cost alternative to current genotyping methods and may be particularly beneficial for in-country surveillance of transmitted drug resistance.

Modeling the emergence of HIV-1 drug resistance resulting from antiretroviral therapy: Insights from theoretical and numerical studies

The use of antiretroviral therapy has proven to be remarkably effective in controlling the progression of human immunodeficiency virus (HIV) infection and prolonging patient's survival. Therapy however may fail and therefore these benefits can be compromised by the emergence of HIV strains that are resistant to the therapy. In view of these facts, the question of finding the reason for which drug-resistant strains emerge during therapy has become a worldwide problem of great interest. This paper presents a deterministic HIV-1 model to examine the mechanisms underlying the emergence of drug-resistance during therapy. The aim of this study is to determine whether, and how fast, antiretroviral therapy may determine the emergence of drug resistance by calculating the basic reproductive numbers. The existence, feasibility and local stability of the equilibriums are also analyzed. By performing numerical simulations we show that Hopf bifurcation may occur. The model suggests that the individuals with drug-resistant infection may play an important role in the epidemic of HIV. (C) 2011 Elsevier Ireland Ltd. All rights reserved.

Molecular evolution of a malaria resistance gene (DARC) in primates

Genes involved in host-pathogen interactions are often strongly affected by positive natural selection. The Duffy antigen, coded by the Duffy antigen receptor for chemokines (DARC) gene, serves as a receptor for Plasmodium vivax in humans and for Plasmodium knowlesi in some nonhuman primates. In the majority of sub-Saharan Africans, a nucleic acid variant in GATA-1 of the gene promoter is responsible for the nonexpression of the Duffy antigen on red blood cells and consequently resistance to invasion by P. vivax. The Duffy antigen also acts as a receptor for chemokines and is expressed in red blood cells and many other tissues of the body. Because of this dual role, we sequenced a 3,000-bp region encompassing the entire DARC gene as well as part of its 5' and 3' flanking regions in a phylogenetic sample of primates and used statistical methods to evaluate the nature of selection pressures acting on the gene during its evolution. We analyzed both coding and regulatory regions of the DARC gene. The regulatory analysis showed accelerated rates of substitution at several sites near known motifs. Our tests of positive selection in the coding region using maximum likelihood by branch sites and maximum likelihood by codon sites did not yield statistically significant evidence for the action of positive selection. However, the maximum likelihood test in which the gene was subdivided into different structural regions showed that the known binding region for P. vivax/P. knowlesi is under very different selective pressures than the remainder of the gene. In fact, most of the gene appears to be under strong purifying selection, but this is not evident in the binding region. We suggest that the binding region is under the influence of two opposing selective pressures, positive selection possibly exerted by the parasite and purifying selection exerted by chemokines.

Exploration of the parameter space in an agent-based model of tuberculosis spread: emergence of drug resistance in developing vs developed countries

In this work we present an agent-based model for the spread of tuberculosis where the individuals can be infected with either drug-susceptible or drug-resistant strains and can also receive a treatment. The dynamics of the model and the role of each one of the parameters are explained. The whole set of parameters is explored to check their importance in the numerical simulation results. The model captures the beneficial impact of the adequate treatment on the prevalence of tuberculosis. Nevertheless, depending on the treatment parameters range, it also captures the emergence of drug resistance. Drug resistance emergence is particularly likely to occur for parameter values corresponding to less efficacious treatment, as usually found in developing countries.

Production of human factor VIII-FL in 293T cells using the bicistronic MGMT(P140K)-retroviral vector

Hemophilia A is the most common X-linked bleeding disorder; it is caused by deficiency of coagulation factor VIII (FVIII). Replacement therapy with rFVIII produced from human cell line is a major goal for treating hemophilia patients. We prepared a full-length recombinant FVIII (FVIII-FL), using the pMFG-P140K retroviral vector. The IRES DNA fragment was cloned upstream to the P140K gene, providing a 9.34-kb bicistronic vector. FVIII-FL cDNA was then cloned upstream to IRES, resulting in a 16.6-kb construct. In parallel, an eGFP control vector was generated, resulting in a 10.1-kb construct. The 293T cells were transfected with these constructs, generating the 293T-FVIII-FL/P140K and 293T-eGFP/P140K cell lines. In 293T-FVIII-FL/P140K cells, FVIII and P140K mRNAs levels were 4,410 (+/- 931.7)- and 295,400 (+/- 75,769)-fold higher than in virgin cells. In 293T-eGFP/P140K cells, the eGFP and P140K mRNAs levels were 1,501,000 (+/- 493,700)- and 308,000 (+/- 139,300)-fold higher than in virgin cells. The amount of FVIII-FL was 0.2 IU/mL and 45 ng/mL FVIII cells or 4.4 IU/mu g protein. These data demonstrate the efficacy of the bicistronic retroviral vector expressing FVIII-FL and MGMT(P140K), showing that it could be used for producing the FVIII-FL protein in a human cell line.

Pretherapeutic Expression of the hOCT1 Gene Predicts a Complete Molecular Response to Imatinib Mesylate in Chronic-Phase Chronic Myeloid Leukemia

In this retrospective study we evaluated the pretherapeutic mRNA expression of the hOCT1 (human organic cation transporter 1) gene in patients with chronic-phase (CP) chronic myeloid leukemia (CML) who varied in terms of their response to imatinib (IM). hOCT1 mRNA was quantified by real-time PCR. Patients were classified as expressing either high (n = 44) or low hOCT1 mRNA (n = 44). The complete cytogenetic response rates observed at 6, 12 and 18 months were 47.7, 84.1 and 91%, respectively, in patients with high hOCT1 mRNA and 47.5, 81.8 and 86.3%, respectively, in patients with low hOCT1 transcripts. The major molecular response rates were not significantly different between patients with high and low hOCT1 mRNA after 6 months of therapy (22.7 vs. 9.1%; p = 0.07), but they were significantly different after 12 months (54.5 vs. 31.8%; p = 0.026) and 18 months (77.2 vs. 56.8%; p = 0.034). Complete molecular responses were observed in 5 patients with low and 17 patients with high hOCT1 mRNA (p = 0.003). The 5-year event-free and overall survival analyses revealed no significant differences between the groups. These data imply that knowledge of the pretherapeutic level of hOCT1 could be a useful marker to predict IM therapy outcome in treatment-naive CP CML patients. Copyright (C) 2012 S. Karger AG, Basel

Evaluation of HIV-1 resistance to antiretroviral drugs among 150 patients after six months of therapeutic interruption

Most of the antiretroviral (ARV) studies in Brazil have been reported in treatment-experienced and naive patients rather than in the setting of treatment interruption (TI). In this study, we analysed reasons given for TI and resistance mutations occurring in 150 HIV-1-infected patients who underwent TI. Of the patients analysed, 110 (73.3%) experienced TI following medical advice, while the remaining patients stopped antiretroviral therapy (ART) of their own accord. The main justifications for TI were: ARV-related toxicities (38.7%), good laboratory parameters (30%) and poor adherence (20%). DNA sequencing of the partial pol gene was successful in 137 (91.3%) patients, of whom 38 (27.7%) presented mutations conferring ARV resistance. A higher viral load prior to TI correlated with drug resistance (P < 0.05). Our results demonstrate that there are diverse rationales for TI and that detection of resistant strains during TI most likely indicates a fitter virus than the wild type. High viral loads coupled with unprotected sex in this group could increase the likelihood of transmission of drug-resistant virus. Thus, treating physicians should be alerted to this problem when the use of ARVs is interrupted.

Single-Nucleotide Polymorphism and Copy Number Variation of the Multidrug Resistance-1 Locus of Plasmodium vivax: Local and Global Patterns

Emerging resistance to chloroquine (CQ) poses a major challenge for Plasmodium vivax malaria control, and nucleotide substitutions and copy number variation in the P. vivax multidrug resistance 1 (pvmdr-1) locus, which encodes a digestive vacuole membrane transporter, may modulate this phenotype. We describe patterns of genetic variation in pvmdr-1 alleles from Acre and Amazonas in northwestern Brazil, and compare then with those reported in other malaria-endemic regions. The pvmdr-1 mutation Y976F, which is associated with CQ resistance in Southeast Asia and Oceania, remains rare in northwestern Brazil (1.8%) and its prevalence mirrors that of CO resistance worldwide. Gene amplification of pvmdr-1, which is associated with mefloquine resistance but increased susceptibility to CO, remains relatively rare in northwestern Brazil (0.9%) and globally (< 4%), but became common (> 10%) in Tak Province, Thailand, possibly because of drug-mediated selection. The global database we have assembled provides a baseline for further studies of genetic variation in pvmdr-1 and drug resistance in P. vivax malaria.

IGFBP7 participates in the reciprocal interaction between acute lymphoblastic leukemia and BM stromal cells and in leukemia resistance to asparaginase

The interaction of acute lymphoblastic leukemia (ALL) blasts with bone marrow (BM) stromal cells (BMSCs) has a positive impact on ALL resistance to chemotherapy. We investigated the modulation of a series of putative asparaginase-resistance/sensitivity genes in B-precursor ALL cells upon coculture with BMSCs. Coculture with stromal cells resulted in increased insulin-like growth factor (IGF)-binding protein 7 (IGFBP7) expression by ALL cells. Assays with IGFBP7 knockdown ALL and stromal cell lines, or with addition of recombinant rIGFBP7 (rIGFBP7) to the culture medium, showed that IGFBP7 acts as a positive regulator of ALL and stromal cells growth, and significantly enhances in-vitro resistance of ALL to asparaginase. In these assays, IGFBP7 function occurred mainly in an insulin-and stromal-dependent manner. ALL cells were found to contribute substantially to extracellular IGFBP7 levels in the conditioned coculture medium. Diagnostic BM plasma from children with ALL had higher levels of IGFBP7 than controls. IGFBP7, in an insulin/IGF-dependent manner, enhanced asparagine synthetase expression and asparagine secretion by BMSCs, thus providing a stromal-dependent mechanism by which IGFBP7 protects ALL cells against asparaginase in this coculture system. Importantly, higher IGFBP7 mRNA levels were associated with lower leukemia-free survival (Cox regression model, P = 0.003) in precursor B-cell Ph(-) ALL patients (n = 147) treated with a contemporary polychemotherapy protocol.

In silico identification of conserved intercoding sequences in Leishmania genomes: Unraveling putative cis-regulatory elements

In silico analyses of Leishmania spp. genome data are a powerful resource to improve the understanding of these pathogens' biology. Trypanosomatids such as Leishmania spp. have their protein-coding genes grouped in long polycistronic units of functionally unrelated genes. The control of gene expression happens by a variety of posttranscriptional mechanisms. The high degree of synteny among Leishmania species is accompanied by highly conserved coding sequences (CDS) and poorly conserved intercoding untranslated sequences. To identify the elements involved in the control of gene expression, we conducted an in silico investigation to find conserved intercoding sequences (CICS) in the genomes of L major, L infantum, and L braziliensis. We used a combination of computational tools, such as Linux-Shell, PERL and R languages, BLAST, MSPcrunch, SSAKE, and Pred-A-Term algorithms to construct a pipeline which was able to: (i) search for conservation in target-regions, (ii) eliminate CICS redundancy and mask repeat elements, (iii) predict the mRNA's extremities, (iv) analyze the distribution of orthologous genes within the generated LeishCICS-clusters, (v) assign GO terms to the LeishCICS-clusters. and (vi) provide statistical support for the gene-enrichment annotation. We associated the LeishCICS-cluster data, generated at the end of the pipeline, with the expression profile oft. donovani genes during promastigote-amastigote differentiation, as previously evaluated by others (GEO accession: GSE21936). A Pearson's correlation coefficient greater than 0.5 was observed for 730 LeishCICS-clusters containing from 2 to 17 genes. The designed computational pipeline is a useful tool and its application identified potential regulatory cis elements and putative regulons in Leishmania. (C) 2012 Elsevier B.V. All rights reserved.

BF Integrase Genes of HIV-1 Circulating in Sao Paulo, Brazil, with a Recurrent Recombination Region

Although some studies have shown diversity in HIV integrase (IN) genes, none has focused particularly on the gene evolving in epidemics in the context of recombination. The IN gene in 157 HIV-1 integrase inhibitor-naive patients from the Sao Paulo State, Brazil, were sequenced tallying 128 of subtype B (23 of which were found in non-B genomes), 17 of subtype F (8 of which were found in recombinant genomes), 11 integrases were BF recombinants, and 1 from subtype C. Crucially, we found that 4 BF recombinant viruses shared a recurrent recombination breakpoint region between positions 4900 and 4924 (relative to the HXB2) that includes 2 gRNA loops, where the RT may stutter. Since these recombinants had independent phylogenetic origin, we argue that these results suggest a possible recombination hotspot not observed so far in BF CRF in particular, or in any other HIV-1 CRF in general. Additionally, 40% of the drug-naive and 45% of the drug-treated patients had at least 1 raltegravir (RAL) or elvitegravir (EVG) resistance-associated amino acid change, but no major resistance mutations were found, in line with other studies. Importantly, V151I was the most common minor resistance mutation among B, F and BF IN genes. Most codon sites of the IN genes had higher rates of synonymous substitutions (dS) indicative of a strong negative selection. Nevertheless, several codon sites mainly in the subtype B were found under positive selection. Consequently, we observed a higher genetic diversity in the B portions of the mosaics, possibly due to the more recent introduction of subtype F on top of an ongoing subtype B epidemics and a fast spread of subtype F alleles among the B population.

Clinical, epidemiological and molecular features of the HIV-1 subtype C and recombinant forms that are circulating in the city of Sao Paulo, Brazil

Background: The city of Sao Paulo has the highest AIDS case rate, with nearly 60% in Brazil. Despite, several studies involving molecular epidemiology, lack of data regarding a large cohort study has not been published from this city. Objectives: This study aimed to describe the HIV-1 subtypes, recombinant forms and drug resistance mutations, according to subtype, with emphasis on subtype C and BC recombinants in the city of Sao Paulo, Brazil. Study design: RNA was extracted from the plasma samples of 302 HIV-1-seropositive subjects, of which 211 were drug-naive and 82 were exposed to ART. HIV-1 partial pol region sequences were used in phylogenetic analyses for subtyping and identification of drug resistance mutations. The envelope gene of subtype C and BC samples was also sequenced. Results: From partial pol gene analyses, 239 samples (79.1%) were assigned as subtype B, 23 (7.6%) were F1, 16 (5.3%) were subtype C and 24 (8%) were mosaics (3 CRF28/CRF29-like). The subtype C and BC recombinants were mainly identified in drug-naive patients (72.7%) and the heterosexual risk exposure category (86.3%), whereas for subtype B, these values were 69.9% and 57.3%, respectively (p = 0.97 and p = 0.015, respectively). An increasing trend of subtype C and BC recombinants was observed (p < 0.01). Conclusion: The HIV-1 subtype C and CRFs seem to have emerged over the last few years in the city of Sao Paulo, principally among the heterosexual population. These findings may have an impact on preventive measures and vaccine development in Brazil.

Analysis of Transmitted Resistance to Raltegravir and Selective Pressure among HIV-1-Infected Patients on a Failing HAART in Sao Paulo, Brazil

We studied the presence of primary resistance to raltegravir (RAL), natural polymorphisms, and selection pressure on HIV-1 integrase. We found a high frequency of integrase polymorphisms related to the resistance to RAL and sequence stability. Further studies are needed to determine the importance of these polymorphisms to RAL resistance.

Clinical, epidemiological and molecular features of the HIV-1 subtype C and recombinant forms that are circulating in the city of São Paulo, Brazil

Abstract Background The city of Sao Paulo has the highest AIDS case rate, with nearly 60% in Brazil. Despite, several studies involving molecular epidemiology, lack of data regarding a large cohort study has not been published from this city. Objectives This study aimed to describe the HIV-1 subtypes, recombinant forms and drug resistance mutations, according to subtype, with emphasis on subtype C and BC recombinants in the city of São Paulo, Brazil. Study design RNA was extracted from the plasma samples of 302 HIV-1-seropositive subjects, of which 211 were drug-naive and 82 were exposed to ART. HIV-1 partial pol region sequences were used in phylogenetic analyses for subtyping and identification of drug resistance mutations. The envelope gene of subtype C and BC samples was also sequenced. Results From partial pol gene analyses, 239 samples (79.1%) were assigned as subtype B, 23 (7.6%) were F1, 16 (5.3%) were subtype C and 24 (8%) were mosaics (3 CRF28/CRF29-like). The subtype C and BC recombinants were mainly identified in drug-naïve patients (72.7%) and the heterosexual risk exposure category (86.3%), whereas for subtype B, these values were 69.9% and 57.3%, respectively (p = 0.97 and p = 0.015, respectively). An increasing trend of subtype C and BC recombinants was observed (p < 0.01). Conclusion The HIV-1 subtype C and CRFs seem to have emerged over the last few years in the city of São Paulo, principally among the heterosexual population. These findings may have an impact on preventive measures and vaccine development in Brazil.

Evaluation of the in vitro Activity of Dermaseptin 01, a Cationic Antimicrobial Peptide, against Schistosoma mansoni

Schistosomiasis is a neglected tropical disease that remains a considerable public health problem worldwide. Since the mainstay of schistosomiasis control is chemotherapy with a single drug, praziquantel, drug resistance is a concern. Here, we examined the in vitro effects of dermaseptin 01 (DS 01), an antimicrobial peptide found in the skin secretion of frogs of the genus Phyllomedusa, on Schistosoma mansoni adult worms. DS 01 at a concentration of 100 mg/ml reduced the worm motor activity and caused the death of all worms within 48 h in RPMI 1640 medium. At the highest sublethal concentration of antimicrobial peptide (75 mg/ml), a 100% reduction in egg output of paired female worms was observed. Additionally, DS 01 induced morphological alterations on the tegument of S. mansoni, and a quantitative analysis carried out by confocal microscopy revealed extensive destruction of the tubercles in a dose-dependent manner over the concentration range of 50-200 mu g/ml. It was the first time that an anthelmintic activity towards schistosomes has been reported for a dermaseptin.