Stochastic processes strongly influence HIV-1 evolution during suboptimal protease-inhibitor therapy

It has long been assumed that HIV-1 evolution is best described by deterministic evolutionary models because of the large population size. Recently, however, it was suggested that the effective population size (Ne) may be rather small, thereby allowing chance to influence evolution, a situation best described by a stochastic evolutionary model. To gain experimental evidence supporting one of the evolutionary models, we investigated whether the development of resistance to the protease inhibitor ritonavir affected the evolution of the env gene. Sequential serum samples from five patients treated with ritonavir were used for analysis of the protease gene and the V3 domain of the env gene. Multiple reverse transcription–PCR products were cloned, sequenced, and used to construct phylogenetic trees and to calculate the genetic variation and Ne. Genotypic resistance to ritonavir developed in all five patients, but each patient displayed a unique combination of mutations, indicating a stochastic element in the development of ritonavir resistance. Furthermore, development of resistance induced clear bottleneck effects in the env gene. The mean intrasample genetic variation, which ranged from 1.2% to 5.7% before treatment, decreased significantly (P < 0.025) during treatment. In agreement with these findings, Ne was estimated to be very small (500–15,000) compared with the total HIV-1 RNA copy number. This study combines three independent observations, strong population bottlenecking, small Ne, and selection of different combinations of protease-resistance mutations, all of which indicate that HIV-1 evolution is best described by a stochastic evolutionary model.

Genetic drift and within-host metapopulation dynamics of HIV-1 infection

Most HIV replication occurs in solid lymphoid tissue, which has prominent architecture at the histological level, which separates groups of productively infected CD4+ cells. Nevertheless, current population models of HIV assume panmixis within lymphoid tissue. We present a simple “metapopulation” model of HIV replication, where the population of infected cells is comprised of a large number of small populations, each of which is established by a few founder viruses and undergoes turnover. To test this model, we analyzed viral genetic variation of infected cell subpopulations within the spleen and demonstrated the action of founder effects as well as significant variation in the extent of genetic differentiation between subpopulations among patients. The combination of founder effects and subpopulation turnover can result in an effective population size much lower than the actual population size and may contribute to the importance of genetic drift in HIV evolution despite a large number of infected cells.

A pre-Columbian Y chromosome-specific transition and its implications for human evolutionary history.

A polymorphic C-->T transition located on the human Y chromosome was found by the systematic comparative sequencing of Y-specific sequence-tagged sites by denaturing high-performance liquid chromatography. The results of genotyping representative global indigenous populations indicate that the locus is polymorphic exclusively within the Western Hemisphere. The pre-Columbian T allele occurs at > 90% frequency within the native South and Central American populations examined, while its occurrence in North America is approximately 50%. Concomitant genotyping at the polymorphic tetranucleotide microsatellite DYS19 locus revealed that the C-->T mutation displayed significant linkage disequilibrium with the 186-bp allele. The data suggest a single origin of linguistically diverse native Americans with subsequent haplotype differentiation within radiating indigenous populations as well as post-Columbian European and African gene flow. The mutation may have originated either in North America at a very early time during the expansion or before it, in the ancestral population(s) from which all Americans may have originated. The analysis of linkage of the DYS199 and the DYS19 tetranucleotide loci suggests that the C-->T mutation may have occurred around 30,000 years ago. We estimate the nucleotide diversity over 4.2 kb of the nonrecombining portion of the Y chromosome to be 0.00014. compared to autosomes, the majority of variation is due to the smaller effective population size of the Y chromosome rather than selective sweeps. There begins to emerge a pattern of pronounced geographical localization of Y-specific nucleotide substitution polymorphisms.

Deriva genética de caracteres quantitativos em milho

A obtenção de genótipos superiores no melhoramento de plantas depende da existência de variabilidade genética. A existência de coleções de germoplasma representativas e a utilização de um tamanho adequado de amostra são fundamentais para a preservação das frequências alélicas e genotípicas, diminuindo a perda de variabilidade genética e postergando o aparecimento dos efeitos da deriva genética. Assim, teve-se como objetivo avaliar os efeitos da deriva genética em caracteres quantitativos em subpopulações de milho. Este estudo foi realizado a partir das populações originais BR-105 e BR-106, das quais 10 subpopulações foram obtidas em cada um dos cinco ciclos sucessivos de amostragem com tamanho efetivo reduzido, totalizando 50 subpopulações para cada população original, as quais foram posteriormente autofecundadas, gerando um nível a mais de endogamia. Os tratamentos foram constituídos de 10 amostras da população original sem autofecundação, 10 amostras com autofecundação, 50 subpopulações obtidas da população original e 50 subpopulações autofecundadas, totalizando 120 tratamentos para cada população, avaliados separadamente. Utilizou-se o delineamento em blocos casualizados no esquema de parcelas subdivididas em faixas hierárquico, em quatro ambientes com duas repetições por ambiente. Os caracteres avaliados foram produção de grãos (PG), prolificidade (PROL), comprimento e diâmetro de espigas (CE e DE), número de fileiras por espiga (NFE), número de grãos por fileira (NGF), altura de planta e espiga (AP e AE), florescimento masculino e feminino (FM e FF) e número de ramificações do pendão (NRP). Foram estimados os efeitos da deriva genética entre as médias das subpopulações nos dois níveis de endogamia e os efeitos da depressão por endogamia nas subpopulações dentro dos ciclos. Posteriormente, realizaram-se análises de regressão linear para as subpopulações nos dois níveis de endogamia, separadamente, e em conjunto. Foi verificada uma grande variação nas médias das subpopulações ao longo dos ciclos, indicando que a deriva genética causou diferenciação entre as mesmas e que estas se diferenciaram das populações originais. Detectaram-se efeitos significativos da deriva genética nas populações não autofecundadas para todos os caracteres avaliados, em maior número para PG, já que este caráter é mais sensível à deriva genética por possuir maior grau de dominância que os demais. Houve diminuição no número de estimativas de deriva significativas para as populações autofecundadas, incluindo mudanças na magnitude e no sinal das mesmas em relação às populações não autofecundadas. Para as estimativas de depressão por endogamia, os caracteres PG, NGF, FM e FF apresentaram maior quantidade de estimativas significativas que os demais. Para a maioria dos caracteres, a regressão linear explicou a maior parte da variação encontrada com o aumento dos coeficientes de endogamia. As populações BR-105 e BR-106, por terem estruturas genéticas distintas, apresentaram performances diferentes quanto aos efeitos da deriva genética. Enfim, como a deriva genética interfere na integridade genética das populações, torna-se importante considerar seus efeitos na coleta e manutenção dos bancos de germoplasma e nas populações utilizadas no melhoramento genético de plantas.

Genetic variation and inference of demographic histories in non-model species

Both long-term environmental changes such as those driven by the glacial cycles and more recent anthropogenic impacts have had major effects on the past demography in wild organisms. Within species, these changes are reflected in the amount and distribution of neutral genetic variation. In this thesis, mitochondrial and microsatellite DNA was analysed to investigate how environmental and anthropogenic factors have affected genetic diversity and structure in four ecologically different animal species. Paper I describes the post-glacial recolonisation history of the speckled-wood butterfly (Pararge aegeria) in Northern Europe. A decrease in genetic diversity with latitude and a marked population structure were uncovered, consistent with a hypothesis of repeated founder events during the postglacial recolonisation. Moreover, Approximate Bayesian Computation analyses indicate that the univoltine populations in Scandinavia and Finland originate from recolonisations along two routes, one on each side of the Baltic. Paper II aimed to investigate how past sea-level rises affected the population history of the convict surgeonfish (Acanthurus triostegus) in the Indo-Pacific. Assessment of the species’ demographic history suggested a population expansion that occurred approximately at the end of the last glaciation. Moreover, the results demonstrated an overall lack of phylogeographic structure, probably due to the high dispersal rates associated with the species’ pelagic larval stage. Populations at the species’ eastern range margin were significantly differentiated from other populations, which likely is a consequence of their geographic isolation. In Paper III, we assessed the effect of human impact on the genetic variation of European moose (Alces alces) in Sweden. Genetic analyses revealed a spatial structure with two genetic clusters, one in northern and one in southern Sweden, which were separated by a narrow transition zone. Moreover, demographic inference suggested a recent population bottleneck. The inferred timing of this bottleneck coincided with a known reduction in population size in the 19th and early 20th century due to high hunting pressure. In Paper IV, we examined the effect of an indirect but well-described human impact, via environmental toxic chemicals (PCBs), on the genetic variation of Eurasian otters (Lutra lutra) in Sweden. Genetic clustering assignment revealed differentiation between otters in northern and southern Sweden, but also in the Stockholm region. ABC analyses indicated a decrease in effective population size in both northern and southern Sweden. Moreover, comparative analyses of historical and contemporary samples demonstrated a more severe decline in genetic diversity in southern Sweden compared to northern Sweden, in agreement with the levels of PCBs found.

Bayesian inferences on the recent island colonization history by the bird Zosterops lateralis lateralis

The founding of new populations by small numbers of colonists has been considered a potentially important mechanism promoting evolutionary change in island populations. Colonizing species, such as members of the avian species complex Zosterops lateralis, have been used to support this idea. A large amount of background information on recent colonization history is available for one Zosterops subspecies, Z. lateralis lateralis, providing the opportunity to reconstruct the population dynamics of its colonization sequence. We used a Bayesian approach to combine historical and demographic information available on Z. l. lateralis with genotypic data from six microsatellite loci, and a rejection algorithm to make simultaneous inferences on the demographic parameters describing the recent colonization history of this subspecies in four southwest Pacific islands. Demographic models assuming mutation–drift equilibrium or a large number of founders were better supported than models assuming founder events for three of four recently colonized island populations. Posterior distributions of demographic parameters supported (i) a large stable effective population size of several thousands individuals with point estimates around 4000–5000; (ii) a founder event of very low intensity with a large effective number of founders around 150–200 individuals for each island in three of four islands, suggesting the colonization of those islands by one flock of large size or several flocks of average size; and (iii) a founder event of higher intensity on Norfolk Island with an effective number of founders around 20 individuals, suggesting colonization by a single flock of moderate size. Our inferences on demographic parameters, especially those on the number of founders, were relatively insensitive to the precise choice of prior distributions for microsatellite mutation processes and demographic parameters, suggesting that our analysis provides a robust description of the recent colonization history of the subspecies.

Genetic analysis of complex demographic scenarios: Spatially expanding populations of the cane toad, Bufo marinus

Inferring the spatial expansion dynamics of invading species from molecular data is notoriously difficult due to the complexity of the processes involved. For these demographic scenarios, genetic data obtained from highly variable markers may be profitably combined with specific sampling schemes and information from other sources using a Bayesian approach. The geographic range of the introduced toad Bufo marinus is still expanding in eastern and northern Australia, in each case from isolates established around 1960. A large amount of demographic and historical information is available on both expansion areas. In each area, samples were collected along a transect representing populations of different ages and genotyped at 10 microsatellite loci. Five demographic models of expansion, differing in the dispersal pattern for migrants and founders and in the number of founders, were considered. Because the demographic history is complex, we used an approximate Bayesian method, based on a rejection-regression algorithm. to formally test the relative likelihoods of the five models of expansion and to infer demographic parameters. A stepwise migration-foundation model with founder events was statistically better supported than other four models in both expansion areas. Posterior distributions supported different dynamics of expansion in the studied areas. Populations in the eastern expansion area have a lower stable effective population size and have been founded by a smaller number of individuals than those in the northern expansion area. Once demographically stabilized, populations exchange a substantial number of effective migrants per generation in both expansion areas, and such exchanges are larger in northern than in eastern Australia. The effective number of migrants appears to be considerably lower than that of founders in both expansion areas. We found our inferences to be relatively robust to various assumptions on marker. demographic, and historical features. The method presented here is the only robust, model-based method available so far, which allows inferring complex population dynamics over a short time scale. It also provides the basis for investigating the interplay between population dynamics, drift, and selection in invasive species.

High nuclear genetic diversity, high levels of outcrossing and low differentiation among remnant populations of Quercus petraea at the margin of its range in Ireland

Background and Aims Quercus petraea colonized Ireland after the last glaciation from refugia on mainland Europe. Deforestation. however. beginning in Neolithic times, has resulted in small, scattered forest fragments, now covering less than 12 000 ha. Methods Plastid (three fragments) and microsatellite variation (13 loci) were characterized in seven Irish populations sampled along a north-south gradient. Using Bayesian approaches and Wright's F-statistics, the effects of colonization and fragmentation on the genetic structure and mating patterns of extant oak populations were investigated. Key-Results All Populations possessed cytotypes common to the Iberian Peninsula. Despite the distance from the refugial core and the extensive deforestation in Ireland, nuclear genetic variation was high and comparable to mainland Europe. Low population differentiation was observed within Ireland and populations showed no evidence for isolation by distance. As expected of a marker with an effective Population size of one-quarter relative to the nuclear genome, plastid variation indicated higher differentiation. Individual inbreeding coefficients indicated high levels of outcrossing. Conclusions Consistent with a large effective Population size in the historical migrant gene pool and/or with high gene flow among populations, high within-population diversity and low population differentiation was observred within Ireland. It is proposed that native Q. petraea populations in Ireland share a common phylogeographic history and that the present genetic structure does not reflect founder effects. (C) 2004 Annals of Botany Company.

The effects of habitat fragmentation due to forestry plantation establishment on the demography and genetic variation of a marsupial carnivore, Antechinus agilis

We conducted a demographic and genetic study to investigate the effects of fragmentation due to the establishment of an exotic softwood plantation on populations of a small marsupial carnivore, the agile antechinus (Antechinus agilis), and the factors influencing the persistence of those populations in the fragmented habitat. The first aspect of the study was a descriptive analysis of patch occupancy and population size, in which we found a patch occupancy rate of 70% among 23 sites in the fragmented habitat compared to 100% among 48 sites with the same habitat characteristics in unfragmented habitat. Mark-recapture analyses yielded most-likely population size estimates of between 3 and 85 among the 16 occupied patches in the fragmented habitat. Hierarchical partitioning and model selection were used to identify geographic and habitat-related characteristics that influence patch occupancy and population size. Patch occupancy was primarily influenced by geographic isolation and habitat quality (vegetation basal area). The variance in population size among occupied sites was influenced primarily by forest type (dominant Eucalyptus species) and, to a lesser extent, by patch area and topographic context (gully sites had larger populations). A comparison of the sex ratios between the samples from the two habitat contexts revealed a significant deficiency of males in the fragmented habitat. We hypothesise that this is due to male-biased dispersal in an environment with increased dispersal-associated mortality. The population size and sex ratio data were incorporated into a simulation study to estimate the proportion of genetic diversity that would have been lost over the known timescale since fragmentation if the patch populations had been totally isolated. The observed difference in genetic diversity (gene diversity and allelic richness at microsatellite and mitochondrial markers) between 16 fragmented and 12 unfragmented sites was extremely low and inconsistent with the isolation of the patch populations. Our results show that although the remnant habitat patches comprise approximately 2% of the study area, they can support non-isolated populations. However, the distribution of agile antechinus populations in the fragmented system is dependent on habitat quality and patch connectivity. (C) 2004 Elsevier Ltd. All rights reserved.

Sociality and the rate of molecular evolution

The molecular clock does not tick at a uniform rate in all taxa but maybe influenced by species characteristics. Eusocial species (those with reproductive division of labor) have been predicted to have faster rates of molecular evolution than their nonsocial relatives because of greatly reduced effective population size; if most individuals in a population are nonreproductive and only one or few queens produce all the offspring, then eusocial animals could have much lower effective population sizes than their solitary relatives, which should increase the rate of substitution of nearly neutral mutations. An earlier study reported faster rates in eusocial honeybees and vespid wasps but failed to correct for phylogenetic nonindependence or to distinguish between potential causes of rate variation. Because sociality has evolved independently in many different lineages, it is possible to conduct a more wide-ranging study to test the generality of the relationship. We have conducted a comparative analysis of 25 phylogenetically independent pairs of social lineages and their nonsocial relatives, including bees, wasps, ants, termites, shrimps, and mole rats, using a range of available DNA sequences (mitochondrial and nuclear DNA coding for proteins and RNAs, and nontranslated sequences). By including a wide range of social taxa, we were able to test whether there is a general influence of sociality on rates of molecular evolution and to test specific predictions of the hypothesis: (1) that social species have faster rates because they have reduced effective population sizes; (2) that mitochondrial genes would show a greater effect of sociality than nuclear genes; and (3) that rates of molecular evolution should be correlated with the degree of sociality. We find no consistent pattern in rates of molecular evolution between social and nonsocial lineages and no evidence that mitochondrial genes show faster rates in social taxa. However, we show that the most highly eusocial Hymenoptera do have faster rates than their nonsocial relatives. We also find that social parasites (that utilize the workers from related species to produce their own offspring) have faster rates than their social relatives, which is consistent with an effect of lower effective population size on rate of molecular evolution. Our results illustrate the importance of allowing for phylogenetic nonindependence when conducting investigations of determinants of variation in rate of molecular evolution.

A multilocus perspective on refugial isolation and divergence in rainforest skinks (Carlia)

To explore the evolutionary consequences of climate-induced fluctuations in distribution of rainforest habitat we contrasted demographic histories of divergence among three lineages of Australian rainforest endemic skinks. The red-throated rainbow skink, Carlia rubrigularis, consists of morphologically indistinguishable northern and southern mitochondrial DNA (mtDNA) lineages that are partially reproductively isolated at their parapatric boundary. The third lineage (C. rhomboidalis) inhabits rainforests just to the south of C. rubrigularis, has blue, rather than red-throated males, and for mtDNA is more closely related to southern C. rubrigularis than is northern C. rubrigularis. Multigene coalescent analyses supported more recent divergence between morphologically distinct lineages than between morphologically conservative lineages. There was effectively no migration and therefore stronger isolation between southern C. rubrigularis and C. rhomboidalis, and low unidirectional migration between morphologically conservative lineages of C. rubrigularis. We found little or no evidence for strong differences in effective population size, and hence different contributions of genetic drift in the demographic history of the three lineages. Overall the results suggest contrasting responses to long-term fluctuations in rainforest habitats, leading to varying opportunities for speciation.

The genetic structure of Australasian green turtles (Chelonia mydas): exploring the geographical scale of genetic exchange

Ecological and genetic studies of marine turtles generally support the hypothesis of natal homing, but leave open the question of the geographical scale of genetic exchange and the capacity of turtles to shift breeding sites. Here we combine analyses of mitochondrial DNA (mtDNA) variation and recapture data to assess the geographical scale of individual breeding populations and the distribution of such populations through Australasia. We conducted multiscale assessments of mtDNA variation among 714 samples from 27 green turtle rookeries and of adult female dispersal among nesting sites in eastern Australia. Many of these rookeries are on shelves that were flooded by rising sea levels less than 10 000 years (c. 450 generations) ago. Analyses of sequence variation among the mtDNA control region revealed 25 haplotypes, and their frequency distributions indicated 17 genetically distinct breeding stocks (Management Units) consisting either of individual rookeries or groups of rookeries in general that are separated by more than 500 km. The population structure inferred from mtDNA was consistent with the scale of movements observed in long-term mark-recapture studies of east Australian rookeries. Phylogenetic analysis of the haplotypes revealed five clades with significant partitioning of sequence diversity (Phi = 68.4) between Pacific Ocean and Southeast Asian/Indian Ocean rookeries. Isolation by distance was indicated for rookeries separated by up to 2000 km but explained only 12% of the genetic structure. The emerging general picture is one of dynamic population structure influenced by the capacity of females to relocate among proximal breeding sites, although this may be conditional on large population sizes as existed historically across this region.

Genetic structure and male-mediated gene flow in the ghost bat (Macroderma gigas)

The Australian ghost bat is a large, opportunistic carnivorous species that has undergone a marked range contraction toward more mesic, tropical sites over the past century. Comparison of mitochondrial DNA (mtDNA) control region sequences and six nuclear microsatellite loci in 217 ghost bats from nine populations across subtropical and tropical Australia revealed strong population subdivision (mtDNA phi(ST) = 0.80; microsatellites URST = 0.337). Low-latitude (tropical) populations had higher heterozygosity and less marked phylogeographic structure and lower subdivision among sites within regions (within Northern Territory [NT] and within North Queensland [NQ]) than did populations at higher latitudes (subtropical sites; central Queensland [CQ]), although sampling of geographically proximal breeding sites is unavoidably restricted for the latter. Gene flow among populations within each of the northern regions appears to be male biased in that the difference in population subdivision for mtDNA and microsatellites (NT phi(ST) = 0.39, URST = 0.02; NQ phi(ST) = 0.60, URST = -0.03) is greater than expected from differences in the effective population size of haploid versus diploid loci. The high level of population subdivision across the range of the ghost bat contrasts with evidence for high gene flow in other chiropteran species and may be due to narrow physiological tolerances and consequent limited availability of roosts for ghost bats, particularly across the subtropical and relatively arid regions. This observation is consistent with the hypothesis that the contraction of the species' range is associated with late Holocene climate change. The extreme isolation among higher-latitude populations may predispose them to additional local extinctions if the processes responsible for the range contraction continue to operate.

Biogeographical distribution of Rimicaris exoculata resident gut epibiont communities along the Mid-Atlantic Ridge hydrothermal vent sites

Rimicaris exoculata is a deep-sea hydrothermal vent shrimp which enlarged gill chamber houses a complex trophic epibiotic community. Its gut harbours an autochthonous and distinct microbial community. This species dominates hydrothermal ecosystems megafauna along the Mid-Atlantic Ridge, regardless of contrasted geochemical conditions prevailing in them. Here, the resident gut epibiont community at four contrasted hydrothermal vent sites (Rainbow/TAG/Logatchev/Ashadze) was analysed and compiled with previous data to evaluate the possible influence of site location, using 16S rRNA surveys and microscopic observations (TEM, SEM and FISH analyses). Filamentous epibionts inserted between the epithelial cells microvilli were observed on all examined samples. Results confirmed resident gut community affiliation to Deferribacteres, Mollicutes, Epsilonproteobacteria and to a lesser extent Gammaproteobacteria lineages. Still a single Deferribacteres phylotype was retrieved at all sites. Four Mollicutes-related OTUs were distinguished, one being only identified on Rainbow specimens. The topology of ribotypes median-joining networks illustrated a community diversification possibly following demographic expansions, suggesting a more ancient evolutionary history and/or a larger effective population size at Rainbow. Finally, the gill chamber community distribution was also analysed through ribotypes networks based on sequences from R. exoculata collected at Rainbow/Snake Pit/TAG/Logatchev/Ashadze sites. Results allow refining hypotheses on the epibiont role and transmission pathways.

Phylodynamic modelling of foot-and-mouth disease virus sequence data

The under-reporting of cases of infectious diseases is a substantial impediment to the control and management of infectious diseases in both epidemic and endemic contexts. Information about infectious disease dynamics can be recovered from sequence data using time-varying coalescent approaches, and phylodynamic models have been developed in order to reconstruct demographic changes of the numbers of infected hosts through time. In this study I have demonstrated the general concordance between empirically observed epidemiological incidence data and viral demography inferred through analysis of foot-and-mouth disease virus VP1 coding sequences belonging to the CATHAY topotype over large temporal and spatial scales. However a more precise and robust relationship between the effective population size (