The Individualized Genetic Barrier Predicts Treatment Response in a Large Cohort of HIV-1 Infected Patients.

The success of combination antiretroviral therapy is limited by the evolutionary escape dynamics of HIV-1. We used Isotonic Conjunctive Bayesian Networks (I-CBNs), a class of probabilistic graphical models, to describe this process. We employed partial order constraints among viral resistance mutations, which give rise to a limited set of mutational pathways, and we modeled phenotypic drug resistance as monotonically increasing along any escape pathway. Using this model, the individualized genetic barrier (IGB) to each drug is derived as the probability of the virus not acquiring additional mutations that confer resistance. Drug-specific IGBs were combined to obtain the IGB to an entire regimen, which quantifies the virus' genetic potential for developing drug resistance under combination therapy. The IGB was tested as a predictor of therapeutic outcome using between 2,185 and 2,631 treatment change episodes of subtype B infected patients from the Swiss HIV Cohort Study Database, a large observational cohort. Using logistic regression, significant univariate predictors included most of the 18 drugs and single-drug IGBs, the IGB to the entire regimen, the expert rules-based genotypic susceptibility score (GSS), several individual mutations, and the peak viral load before treatment change. In the multivariate analysis, the only genotype-derived variables that remained significantly associated with virological success were GSS and, with 10-fold stronger association, IGB to regimen. When predicting suppression of viral load below 400 cps/ml, IGB outperformed GSS and also improved GSS-containing predictors significantly, but the difference was not significant for suppression below 50 cps/ml. Thus, the IGB to regimen is a novel data-derived predictor of treatment outcome that has potential to improve the interpretation of genotypic drug resistance tests.

The Individualized Genetic Barrier Predicts Treatment Response in a Large Cohort of HIV-1 Infected Patients

The success of combination antiretroviral therapy is limited by the evolutionary escape dynamics of HIV-1. We used Isotonic Conjunctive Bayesian Networks (I-CBNs), a class of probabilistic graphical models, to describe this process. We employed partial order constraints among viral resistance mutations, which give rise to a limited set of mutational pathways, and we modeled phenotypic drug resistance as monotonically increasing along any escape pathway. Using this model, the individualized genetic barrier (IGB) to each drug is derived as the probability of the virus not acquiring additional mutations that confer resistance. Drug-specific IGBs were combined to obtain the IGB to an entire regimen, which quantifies the virus' genetic potential for developing drug resistance under combination therapy. The IGB was tested as a predictor of therapeutic outcome using between 2,185 and 2,631 treatment change episodes of subtype B infected patients from the Swiss HIV Cohort Study Database, a large observational cohort. Using logistic regression, significant univariate predictors included most of the 18 drugs and single-drug IGBs, the IGB to the entire regimen, the expert rules-based genotypic susceptibility score (GSS), several individual mutations, and the peak viral load before treatment change. In the multivariate analysis, the only genotype-derived variables that remained significantly associated with virological success were GSS and, with 10-fold stronger association, IGB to regimen. When predicting suppression of viral load below 400 cps/ml, IGB outperformed GSS and also improved GSS-containing predictors significantly, but the difference was not significant for suppression below 50 cps/ml. Thus, the IGB to regimen is a novel data-derived predictor of treatment outcome that has potential to improve the interpretation of genotypic drug resistance tests.

Comparison of different methods to detect genetic barriers in a small mammal population

Habitat fragmentation and the consequently the loss of connectivity between populations can reduce the individuals interchange and gene flow, increasing the chances of inbreeding, and the increase the risk of local extinction. Landscape genetics is providing more and better tools to identify genetic barriers.. To our knowledge, no comparison of methods in terms of consistency has been made with observed data and species with low dispersal ability. The aim of this study is to examine the consistency of the results of five methods to detect barriers to gene flow in a Mediterranean pine vole population Microtus duodecimcostatus: F-statistics estimations, Non-Bayesian clustering, Bayesian clustering, Boundary detection and Simple/Partial Mantel tests. All methods were consistent in detecting the stream as a non-genetic barrier. However, no consistency in results among the methods were found regarding the role of the highway as a genetic barrier. Fst, Bayesian clustering assignment test and Partial Mantel test identifyed the highway as a filter to individual interchange. The Mantel tests were the most sensitive method. Boundary detection method (Monmonier’s Algorithm) and Non-Bayesian approaches did not detect any genetic differentiation of the pine vole due to the highway. Based on our findings we recommend that the genetic barrier detection in low dispersal ability populations should be analyzed with multiple methods such as Mantel tests, Bayesian clustering approaches because they show more sensibility in those scenarios and with boundary detection methods by having the aim of detect drastic changes in a variable of interest between the closest individuals. Although simulation studies highlight the weaknesses and the strengths of each method and the factors that promote some results, tests with real data are needed to increase the effectiveness of genetic barrier detection.

Differential patterns of Male and Female mtDNA exchange across the Atlantic Ocean in the blue mussel, Mytilus edulis

Comparisons among loci with differing modes of inheritance can reveal unexpected aspects of population history. We employ a multilocus approach to ask whether two types of independently assorting mitochondrial DNAs (maternally and paternally inherited: F- and M-mtDNA) and a nuclear locus (ITS) yield concordant estimates of gene flow and population divergence. The blue mussel, Mytilus edulis, is distributed on both North American and European coastlines and these populations are separated by the waters of the Atlantic Ocean. Gene flow across the Atlantic Ocean differs among loci, with F-mtDNA and ITS showing an imprint of some genetic interchange and M-mtDNA showing no evidence for gene flow. Gene flow of F-mtDNA and ITS causes trans-Atlantic population divergence times to be greatly underestimated for these loci, although a single trans-Atlantic population divergence time (1.2 MYA) can be accommodated by considering all three loci in combination in a coalescent framework. The apparent lack of gene flow for M-mtDNA is not readily explained by different dispersal capacities of male and female mussels. A genetic barrier to M-mtDNA exchange between North American and European mussel populations is likely to explain the observed pattern, perhaps associated with the double uniparental system of mitochondrial DNA inheritance.

The Role of Second Generation Antiretroviral Drugs in HIV-1 Subtype B and non-B Variants Harboring Natural Polymorphisms and Drug Resistance Mutations.

Cette thèse traite de la résistance du VIH-1 aux antirétroviraux, en particulier de l'activité antivirale de plusieurs inhibiteurs non nucléosidiques de la transcriptase inverse (INNTI) ainsi que des inhibiteurs de protéase (IP). Nous avons exploré l’émergence et la spécificité des voies de mutations qui confèrent la résistance contre plusieurs nouveaux INNTI (étravirine (ETR) et rilpivirine (RPV)) (chapitres 2 et 3). En outre, le profil de résistance et le potentiel antirétroviral d'un nouvel IP, PL-100, est présenté dans les chapitres 4 et 5. Pour le premier projet, nous avons utilisé des sous-types B et non-B du VIH-1 pour sélectionner des virus résistants à ETR, et ainsi montré que ETR favorise l’émergence des mutations V90I, K101Q, E138K, V179D/E/F, Y181C, V189I, G190E, H221H/Y et M230L, et ce, en 18 semaines. Fait intéressant, E138K a été la première mutation à émerger dans la plupart des cas. Les clones viraux contenant E138K ont montré un faible niveau de résistance phénotypique à ETR (3,8 fois) et une diminution modeste de la capacité de réplication (2 fois) par rapport au virus de type sauvage. Nous avons également examiné les profils de résistance à ETR et RPV dans les virus contenant des mutations de résistance aux INNTI au début de la sélection. Dans le cas du virus de type sauvage et du virus contenant la mutation unique K103N, les premières mutations à apparaître en présence d’ETR ou de RPV ont été E138K ou E138G suivies d’autres mutations de résistance aux INNTI. À l’inverse, dans les mêmes conditions, le virus avec la mutation Y181C a évolué pour produire les mutations V179I/F ou A62V/A, mais pas E138K/G. L'ajout de mutations à la position 138 en présence de Y181C n'augmente pas les niveaux de résistance à ETR ou RPV. Nous avons également observé que la combinaison de Y181C et E138K peut conduire à un virus moins adapté par rapport au virus contenant uniquement Y181C. Sur la base de ces résultats, nous suggérons que les mutations Y181C et E138K peuvent être antagonistes. L’analyse de la résistance au PL-100 des virus de sous-type C et CRF01_AE dans les cellules en culture est décrite dans le chapitre 4. Le PL-100 sélectionne pour des mutations de résistance utilisant deux voies distinctes, l'une avec les mutations V82A et L90M et l'autre avec T80I, suivi de l’addition des mutations M46I/L, I54M, K55R, L76F, P81S et I85V. Une accumulation d'au moins trois mutations dans le rabat protéique et dans le site actif est requise dans chaque cas pour qu’un haut niveau de résistance soit atteint, ce qui démontre que le PL-100 dispose d'une barrière génétique élevée contre le développement de la résistance. Dans le chapitre 5, nous avons évalué le potentiel du PL-100 en tant qu’inhibiteur de protéase de deuxième génération. Les virus résistants au PL-100 émergent en 8-48 semaines alors qu’aucune mutation n’apparaît avec le darunavir (DRV) sur une période de 40 semaines. La modélisation moléculaire montre que la haute barrière génétique du DRV est due à de multiples interactions avec la protéase dont des liaison hydrogènes entre les groupes di-tétrahydrofuranne (THF) et les atomes d'oxygène des acides aminés A28, D29 et D30, tandis que la liaison de PL-100 est principalement basée sur des interactions polaires et hydrophobes délocalisées à travers ses groupes diphényle. Nos données suggèrent que les contacts de liaison hydrogène et le groupe di-THF dans le DRV, ainsi que le caractère hydrophobe du PL-100, contribuent à la liaison à la protéase ainsi qu’à la haute barrière génétique contre la résistance et que la refonte de la structure de PL-100 pour inclure un groupe di-THF pourrait améliorer l’activité antivirale et le profil de résistance.

Mutations in H5N1 influenza virus hemagglutinin that confer binding to human tracheal airway epithelium

The emergence in 2009 of a swine-origin H1N1 influenza virus as the first pandemic of the 21st Century is a timely reminder of the international public health impact of influenza viruses, even those associated with mild disease. The widespread distribution of highly pathogenic H5N1 influenza virus in the avian population has spawned concern that it may give rise to a human influenza pandemic. The mortality rate associated with occasional human infection by H5N1 virus approximates 60%, suggesting that an H5N1 pandemic would be devastating to global health and economy. To date, the H5N1 virus has not acquired the propensity to transmit efficiently between humans. The reasons behind this are unclear, especially given the high mutation rate associated with influenza virus replication. Here we used a panel of recombinant H5 hemagglutinin (HA) variants to demonstrate the potential for H5 HA to bind human airway epithelium, the predominant target tissue for influenza virus infection and spread. While parental H5 HA exhibited limited binding to human tracheal epithelium, introduction of selected mutations converted the binding profile to that of a current human influenza strain HA. Strikingly, these amino-acid changes required multiple simultaneous mutations in the genomes of naturally occurring H5 isolates. Moreover, H5 HAs bearing intermediate sequences failed to bind airway tissues and likely represent mutations that are an evolutionary "dead end." We conclude that, although genetic changes that adapt H5 to human airways can be demonstrated, they may not readily arise during natural virus replication. This genetic barrier limits the likelihood that current H5 viruses will originate a human pandemic.

Restrictions to the adaptation of influenza A virus H5 hemagglutinin to the human host

The binding specificities of a panel of avian influenza virus subtype H5 hemagglutinin (RA) proteins bearing mutations at key residues in the receptor binding site were investigated. The results demonstrate that two simultaneous mutations in the receptor binding site resulted in H5 RA binding in a pattern similar to that shown by human viruses. Coexpression of the ion channel protein, M2, from most avian and human strains tested protected H5 RA conformation during trafficking, indicating that no genetic barrier to the reassortment of the H5 surface antigen gene with internal genes of human viruses existed at this level.

Darunavir: A Critical Review of Its Properties, Use and Drug Interactions

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Avaliação do perfil de mutações e resistência aos inibidores de transcriptase reversa e protease em pacientes pediátricos infectados pleo HIV-1

Pós-graduação em Pesquisa e Desenvolvimento (Biotecnologia Médica) - FMB

Effects of HIV type-1 immune selection on susceptability to integrase inhibitor resistance

BACKGROUND: All site-specific interactions between HIV type-1 (HIV-1) subtype, human leukocyte antigen (HLA)-associated immune selection and integrase inhibitor resistance are not completely understood. We examined naturally occurring polymorphisms in HIV-1 integrase sequences from 342 antiretroviral-naive individuals from the Western Australian HIV Cohort Study and the Swiss HIV Cohort Study. METHODS: Standard bulk sequencing and sequence-based typing were used to generate integrase sequences and high-resolution HLA genotypes, respectively. Viral residues were examined with respect to drug resistance mutations and CD8(+) T-cell escape mutations. RESULTS: In both predominantly subtype B cohorts, 12 of 38 sites that mediate integrase inhibitor resistance mutations were absolutely conserved, and these included the primary resistance mutations. There were 18 codons with non-primary drug resistance-associated substitutions at rates of up to 58.8% and eight sites with alternative polymorphisms. Five viral residues were potentially subject to dual-drug and HLA-associated immune selection in which both selective pressures either drove the same amino acid substitution (codons 72, 157 and 163) or HLA alleles were associated with an alternative polymorphism that would alter the genetic barrier to resistance (codons 125 and 193). The common polymorphism T125A, which was characteristic of non-subtype B and was also associated with carriage of HLA-B*57/*5801, increased the mutational barrier to the resistance mutation T125K. CONCLUSIONS: Primary integrase inhibitor resistance mutations were not detected in the absence of drug exposure in keeping with sites of high constraint. Viral polymorphisms caused by immune selection and/or associated with non-subtype B might alter the genetic barrier to some non-primary resistance-associated mutations.

Tight Genetic Linkage of Prezygotic Barrier Loci Creates a Multifunctional Speciation Island in Petunia

Most flowering plants depend on animal vectors for pollination and seed dispersal. Differential pollinator preferences lead to premating isolation and thus reduced gene flow between interbreeding plant populations [1, 2, 3 and 4]. Sets of floral traits, adapted to attract specific pollinator guilds, are called pollination syndromes [5]. Shifts in pollination syndromes have occurred surprisingly frequently [6], considering that they must involve coordinated changes in multiple genes affecting multiple floral traits. Although the identification of individual genes specifying single pollination syndrome traits is in progress in many species, little is known about the genetic architecture of coadapted pollination syndrome traits and how they are embedded within the genome [7]. Here we describe the tight genetic linkage of loci specifying five major pollination syndrome traits in the genus Petunia: visible color, UV absorption, floral scent production, pistil length, and stamen length. Comparison with other Solanaceae indicates that, in P. exserta and P. axillaris, loci specifying these floral traits have specifically become clustered into a multifunctional “speciation island” [ 8 and 9]. Such an arrangement promotes linkage disequilibrium and avoids the dissolution of pollination syndromes by recombination. We suggest that tight genetic linkage provides a mechanism for rapid switches between distinct pollination syndromes in response to changes in pollinator availabilities.

Genetic structure of a reef-building coral from thermally distinct environments on the Great Barrier Reef

Adaptation to localised thermal regimes is facilitated by restricted gene flow, ultimately leading to genetic divergence among populations and differences in their physiological tolerances. Allozyme analysis of six polymorphic loci was used to assess genetic differentiation between nine populations of the reef-building coral Acropora millepora over a latitudinal temperature gradient on the inshore regions of the Great Barrier Reef (GBR). Small but significant genetic differentiation indicative of moderate levels of gene flow (pairwise F-ST 0.023 to 0.077) was found between southern populations of A. millepora in cooler regions of the GBR and the warmer, central or northern GBR populations. Patterns of genetic differentiation at these putatively neutral allozyme loci broadly matched experimental variation in thermal tolerance and were consistent with local thermal regimes (warmest monthly-averages) for the A. millepora populations examined. It is therefore hypothesized that natural selection has influenced the thermal tolerance of the A. millepora populations examined and greater genetic divergence is likely to be revealed by examination of genetic markers under the direct effects of natural selection.

Macrogeographic Genetic Variation in Broad-Snouted Caiman (Caiman latirostris)

Broad-snouted caiman`s (Caiman latirostris) geographic distribution comprises one of the widest latitudinal ranges among all crocodilians. In this study we analyzed the relationship between geographic distance (along the species latitudinal range) and genetic differentiation using DNA microsatellite loci developed for C. latirostris and Alligator mississippiensis. The results suggest that there is a consistent relationship between geographic distance and genetic differentiation; however, other biogeographical factors seem to be relevant. The Atlantic Chain (Serra do Mar) seems to be an effective geographic barrier, as well as the relatively narrow (<= 1.5 km) sea channel between Cardoso Island and the continent. In addition, coastal populations seem to have been well connected in recent geological time (Pleistocene 16,000 years ago) all along the eastern Brazilian coast. Further studies should focus on the Sao Francisco River drainage, which is still poorly known for this species. J. Exp. Zool. 30.9A:628-636, 2008. (C) 2008 Wiley-Liss, Inc.

Low symbiont diversity in southern Great Barrier Reef corals, relative to those of the Caribbean

The specific identity of endosymbiotic dinoflagellates (Symbiodinium spp.) from most zooxanthellate corals is unknown. In a survey of symbiotic cnidarians from the southern Great Barrier Reef (GBR), 23 symbiont types were identified from 86 host species representing 40 genera. A majority (>85%) of these symbionts belong to a single phylogenetic clade or subgenus (C) composed of closely related (as assessed by sequence data from the internal transcribed spacer region and the ribosomal large subunit gene), yet ecologically and physiologically distinct, types. A few prevalent symbiont types, or generalists, dominate the coral community of the southern GBR, whereas many rare and/or specific symbionts, or specialists, are found uniquely within certain host taxa. The comparison of symbiont diversity between southern GBR and Caribbean reefs shows an inverse relationship between coral diversity and symbiont diversity, perhaps as a consequence of more-rapid diversification of Caribbean symbionts. Among clade C types, generalists C1 and C3 are common to both Caribbean and southern GBR symbiont assemblages, whereas the rest are regionally endemic. Possibly because of environmental changes in the Caribbean after geographic isolation through the Quaternary period, a high proportion of Caribbean fauna associate with symbiont taxa from two other distantly related Symbiodinium clades (A and B) that rarely occur in Pacific hosts. The resilience of Porites spp. and the resistance of Montipora digitata to thermal stress and bleaching are partially explained by their association with a thermally tolerant symbiont type, whereas the indiscriminant widespread bleaching and death among certain Pacific corals, during El Nino Southern Oscillation events, are influenced by associations with symbionts possessing higher sensitivity to thermal stress.

Genetic population structure of the catadromous perciform: Macquaria novemaculeata (Percichthyidae)

Genetic population structure in the catadromous Australian bass Macquaria novemaculeata was investigated using samples from four locations spanning 600 km along the eastern Australian coastline. Both allozymes and mtDNA control region sequences were examined. Population subdivision estimates based on allozymes revealed low levels of population structuring (G(st)=0.043, P<0.05). However, mtDNA indicated moderate levels of geographic population structure (G(st)=0.146, P<0.01). Phylogenetic analysis of mtDNA control region sequences (mean sequence divergence 1.9%) indicated little phylogeographic structuring. Results suggested that genotypic variation within each river population, while bring affected primarily by genetic drift, was also prevented from more significant divergence by homogenizing levels of gene flow-synonymous with a one-dimensional stepping-stone model of population structure. The catadromous life history of Macquaria novemaculeata was considered to br influential on the pattern of population structure displayed. Results were compared to the few population genetic studies involving catadromous fishes, indicating that catadromy alone is unlikely to be a good predictor of population structure. A more comprehensive suite of biological characteristics than simple life-history traits must be considered fully to allow reliable predictive models of population structure to be formulated. (C) 1997 The Fisheries Society of the British Isles.