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.

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.

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.

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.

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.

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.

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.

A study of multidrug-resistant tuberculosis in risk groups in the city of Santos, Sao Paulo, Brazil

Monitoring the extent of and trends in multidrug-resistant tuberculosis (MDR-TB) is a priority of the Brazilian National Tuberculosis Control Programme. The current study aimed to estimate the incidence of MDR-TB, describe the profile of TB drug resistance in risk groups and examine whether screening for MDR-TB adhered to the recommended guidelines. A descriptive study that examined diagnosed cases of pulmonary TB was conducted in the city of Santos, Brazil, between 2000-2004. Of the 2,176 pulmonary TB cases studied, 671 (30.8%) met the criteria for drug sensitivity testing and, of these cases, 31.7% (213/671) were tested. Among the tested cases, 9.4% were resistant to one anti-TB drug and 15% were MDR. MDR was observed in 11.6% of 86 new TB cases and 17.3% of 127 previously treated cases. The average annual incidence of MDR-TB was 1.9 per 100,000 inhabitants-years. The extent of known MDR-TB in the city of Santos is high, though likely to be underestimated. Our study therefore indicates an inadequate adherence to the guidelines for MDR-TB screening and suggests the necessity of alternative strategies of MDR-TB surveillance.

Whole genome sequencing analysis of Plasmodium vivax using whole genome capture

Background: Malaria caused by Plasmodium vivax is an experimentally neglected severe disease with a substantial burden on human health. Because of technical limitations, little is known about the biology of this important human pathogen. Whole genome analysis methods on patient-derived material are thus likely to have a substantial impact on our understanding of P. vivax pathogenesis and epidemiology. For example, it will allow study of the evolution and population biology of the parasite, allow parasite transmission patterns to be characterized, and may facilitate the identification of new drug resistance genes. Because parasitemias are typically low and the parasite cannot be readily cultured, on-site leukocyte depletion of blood samples is typically needed to remove human DNA that may be 1000X more abundant than parasite DNA. These features have precluded the analysis of archived blood samples and require the presence of laboratories in close proximity to the collection of field samples for optimal pre-cryopreservation sample preparation. Results: Here we show that in-solution hybridization capture can be used to extract P. vivax DNA from human contaminating DNA in the laboratory without the need for on-site leukocyte filtration. Using a whole genome capture method, we were able to enrich P. vivax DNA from bulk genomic DNA from less than 0.5% to a median of 55% (range 20%-80%). This level of enrichment allows for efficient analysis of the samples by whole genome sequencing and does not introduce any gross biases into the data. With this method, we obtained greater than 5X coverage across 93% of the P. vivax genome for four P. vivax strains from Iquitos, Peru, which is similar to our results using leukocyte filtration (greater than 5X coverage across 96% of the genome). Conclusion: The whole genome capture technique will enable more efficient whole genome analysis of P. vivax from a larger geographic region and from valuable archived sample collections.

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

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.

Phylogeographic Analysis of HIV-1 Subtype C Dissemination in Southern Brazil

The HIV-1 subtype C has spread efficiently in the southern states of Brazil (Rio Grande do Sul, Santa Catarina and Parana). Phylogeographic studies indicate that the subtype C epidemic in southern Brazil was initiated by the introduction of a single founder virus population at some time point between 1960 and 1980, but little is known about the spatial dynamics of viral spread. A total of 135 Brazilian HIV-1 subtype C pol sequences collected from 1992 to 2009 at the three southern state capitals (Porto Alegre, Florianopolis and Curitiba) were analyzed. Maximum-likelihood and Bayesian methods were used to explore the degree of phylogenetic mixing of subtype C sequences from different cities and to reconstruct the geographical pattern of viral spread in this country region. Phylogeographic analyses supported the monophyletic origin of the HIV-1 subtype C clade circulating in southern Brazil and placed the root of that clade in Curitiba (Parana state). This analysis further suggested that Florianopolis (Santa Catarina state) is an important staging post in the subtype C dissemination displaying high viral migration rates from and to the other cities, while viral flux between Curitiba and Porto Alegre (Rio Grande do Sul state) is very low. We found a positive correlation (r(2) = 0.64) between routine travel and viral migration rates among localities. Despite the intense viral movement, phylogenetic intermixing of subtype C sequences from different Brazilian cities is lower than expected by chance. Notably, a high proportion (67%) of subtype C sequences from Porto Alegre branched within a single local monophyletic sub-cluster. These results suggest that the HIV-1 subtype C epidemic in southern Brazil has been shaped by both frequent viral migration among states and in situ dissemination of local clades.

Molecular identification of Pichia guilliermondii, Debaryomyces hansenii and Candida palmioleophila

Traditional phenotypic methods and commercial kits based on carbohydrate assimilation patterns are unable to consistently distinguish among isolates of Pichia guilliermondii, Debaryomyces hansenii and Candida palmioleophila. As result, these species are often misidentified. In this work, we established a reliable method for the identification/differentiation of these species. Our assay was validated by DNA sequencing of the polymorphic region used in a real-time PCR assay driven by species-specific probes targeted to the fungal ITS 1 region. This assay provides a new tool for pathogen identification and for epidemiological, drug resistance and virulence studies of these organisms.