The Carpathian range represents a weak genetic barrier in South-East Europe

Background: In the present study we have assessed whether the Carpathian Mountains represent a genetic barrier in East Europe. Therefore, we have analyzed the mtDNA of 128 native individuals of Romania: 62 of them from the North of Romania, and 66 from South Romania. Results: We have analyzed their mtDNA variability in the context of other European and Near Eastern populations through multivariate analyses. The results show that regarding the mtDNA haplogroup and haplotype distributions the Romanian groups living outside the Carpathian range (South Romania) displayed some degree of genetic differentiation compared to those living within the Carpahian range (North Romania). Conclusion: The main differentiation between the mtDNA variability of the groups from North and South Romania can be attributed to the demographic movements from East to West (prehistoric or historic) that differently affected in these regions, suggesting that the Carpathian mountain range represents a weak genetic barrier in South-East Europe.

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.

Timing the Emergence of Resistance to Anti-HIV Drugs with Large Genetic Barriers

New antiretroviral drugs that offer large genetic barriers to resistance, such as the recently approved inhibitors of HIV-1 protease, tipranavir and darunavir, present promising weapons to avert the failure of current therapies for HIV infection. Optimal treatment strategies with the new drugs, however, are yet to be established. A key limitation is the poor understanding of the process by which HIV surmounts large genetic barriers to resistance. Extant models of HIV dynamics are predicated on the predominance of deterministic forces underlying the emergence of resistant genomes. In contrast, stochastic forces may dominate, especially when the genetic barrier is large, and delay the emergence of resistant genomes. We develop a mathematical model of HIV dynamics under the influence of an antiretroviral drug to predict the waiting time for the emergence of genomes that carry the requisite mutations to overcome the genetic barrier of the drug. We apply our model to describe the development of resistance to tipranavir in in vitro serial passage experiments. Model predictions of the times of emergence of different mutant genomes with increasing resistance to tipranavir are in quantitative agreement with experiments, indicating that our model captures the dynamics of the development of resistance to antiretroviral drugs accurately. Further, model predictions provide insights into the influence of underlying evolutionary processes such as recombination on the development of resistance, and suggest guidelines for drug design: drugs that offer large genetic barriers to resistance with resistance sites tightly localized on the viral genome and exhibiting positive epistatic interactions maximally inhibit the emergence of resistant genomes.

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.

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.

Is life history a barrier to dispersal? Contrasting patterns of genetic differentiation along an oceanographically complex coast

Extreme variation in early life-history strategies is considered a moderately good predictor of genetic subdivision and hence dispersal for a range of marine species. In reality, however, a good deal of population differentiation must reflect historical effects, more subtle variation in life histories, and, particularly, the interaction of larvae with oceanographic processes. Using a combination of allozyme and microsatellite markers, we show that the large-scale genetic structure of populations of three species (direct and planktonically developing cushion stars and a planktonic developing sea anemone that is also asexually viviparous) varies consistently, in line with the predicted capacity for dispersal within three geographic regions. We detected high levels of genetic subdivision for the direct developing cushion star (FST = 0.6), low levels for the planktonically developing cushion star (FST = 0.009), and intermediate levels for the sexual/asexual sea anmone (FST = 0.19). These patterns are exhibited despite the highly variable patterns of current movement and the presence of biogeographic barriers. Our results suggest that, although there is large scale genetic differentiation for two species, patterns of population connectivity are remarkably consistent within major regions and do not reflect variation in major oceanographic processes or genetic discontinuity coincident with biogeographic boundaries.

Ocean currents influence the genetic structure of an intertidal mollusc in southeastern Australia – implications for predicting the movement of passive dispersers across a marine biogeographic barrier

Major disjunctions among marine communities in southeastern Australia have been well documented, although explanations for biogeographic structuring remain uncertain. Converging ocean currents, environmental gradients, and habitat discontinuities have been hypothesized as likely drivers of structuring in many species, although the extent to which species are affected appears largely dependent on specific life histories and ecologies. Understanding these relationships is critical to the management of native and invasive species, and the preservation of evolutionary processes that shape biodiversity in this region. In this study we test the direct influence of ocean currents on the genetic structure of a passive disperser across a major biogeographic barrier. Donax deltoides (Veneroida: Donacidae) is an intertidal, soft-sediment mollusc and an ideal surrogate for testing this relationship, given its lack of habitat constraints in this region, and its immense dispersal potential driven by year-long spawning and long-lived planktonic larvae. We assessed allele frequencies at 10 polymorphic microsatellite loci across 11 sample locations spanning the barrier region and identified genetic structure consistent with the major ocean currents of southeastern Australia. Analysis of mitochondrial DNA sequence data indicated no evidence of genetic structuring, but signatures of a species range expansion corresponding with historical inundations of the Bassian Isthmus. Our results indicate that ocean currents are likely to be the most influential factor affecting the genetic structure of D. deltoides and a likely physical barrier for passive dispersing marine fauna generally in southeastern Australia.

Differing impact of a major biogeographic barrier on genetic structure in two large kangaroos from the monsoon tropics of Northern Australia

Tropical savannas cover 20-30% of the world's land surface and exhibit high levels of regional endemism, but the evolutionary histories of their biota remain poorly studied. The most extensive and unmodified tropical savannas occur in Northern Australia, and recent studies suggest this region supports high levels of previously undetected genetic diversity. To examine the importance of barriers to gene flow and the environmental history of Northern Australia in influencing patterns of diversity, we investigated the phylogeography of two closely related, large, vagile macropodid marsupials, the antilopine wallaroo (Macropus antilopinus; n = 78), and the common wallaroo (Macropus robustus; n = 21). Both species are widespread across the tropical savannas of Australia except across the Carpentarian Barrier (CB) where there is a break in the distribution of M. antilopinus. We determined sequence variation in the hypervariable Domain I of the mitochondrial DNA control region and genotyped individuals at 12 polymorphic microsatellite loci to assess the historical and contemporary influence of the CB on these species. Surprisingly, we detected only limited differentiation between the disjunct Northern Territory and QueenslandM. antilopinus populations. In contrast, the continuously distributedM. robustus was highly divergent across the CB. Although unexpected, these contrasting responses appear related to minor differences in species biology. Our results suggest that vicariance may not explain well the phylogeographic patterns in Australia's dynamic monsoonal environments. This is because Quaternary environmental changes in this region have been complex, and diverse individual species' biologies have resulted in less predictable and idiosyncratic responses.

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.