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.

Gene genealogies and population variation in plants

Early in the development of plant evolutionary biology, genetic drift, fluctuations in population size, and isolation were identified as critical processes that affect the course of evolution in plant species. Attempts to assess these processes in natural populations became possible only with the development of neutral genetic markers in the 1960s. More recently, the application of historically ordered neutral molecular variation (within the conceptual framework of coalescent theory) has allowed a reevaluation of these microevolutionary processes. Gene genealogies trace the evolutionary relationships among haplotypes (alleles) with populations. Processes such as selection, fluctuation in population size, and population substructuring affect the geographical and genealogical relationships among these alleles. Therefore, examination of these genealogical data can provide insights into the evolutionary history of a species. For example, studies of Arabidopsis thaliana have suggested that this species underwent rapid expansion, with populations showing little genetic differentiation. The new discipline of phylogeography examines the distribution of allele genealogies in an explicit geographical context. Phylogeographic studies of plants have documented the recolonization of European tree species from refugia subsequent to Pleistocene glaciation, and such studies have been instructive in understanding the origin and domestication of the crop cassava. Currently, several technical limitations hinder the widespread application of a genealogical approach to plant evolutionary studies. However, as these technical issues are solved, a genealogical approach holds great promise for understanding these previously elusive processes in plant evolution.

The vesicle size of DDA:TDB liposomal adjuvants plays a role in the cell-mediated immune response but has no significant effect on antibody production

The use of cationic liposomes as experimental adjuvants for subunit peptide of protein vaccines is well documented. Recently the cationic liposome CAF01, composed of dimethyldioctadecylammonium (DDA) and trehalose dibehenate (TDB), has entered Phase I clinical trials for use in a tuberculosis (TB) vaccine. CAF01 liposomes are a heterogeneous population with a mean vesicle size of 500 nm; a strong retention of antigen at the injection site and a Th1-biassed immune response are noted. The purpose of this study was to investigate whether CAF01 liposomes of significantly different vesicle sizes exhibited altered pharmacokinetics in vivo and cellular uptake with activation in vitro. Furthermore, the immune response against the TB antigen Ag85B-ESAT-6 was followed when various sized CAF01 liposomes were used as vaccine adjuvants. The results showed no differences in vaccine (liposome or antigen) draining from the injection site, however, significant differences in the movement of liposomes to the popliteal lymph node were noted. Liposome uptake by THP-1 vitamin D3 stimulated macrophage-like cells did not show a liposome size-dependent pattern of uptake. Finally, whilst there were no significant differences in the IgG1/2 regardless of the liposome size used as a delivery vehicle for Ag85B-ESAT-6, vesicle size has a size dependent effect on cell proliferation and IL-10 production with larger liposomes (in excess of 2 µm) promoting the highest proliferation and lowest IL-10 responses, yet vesicles of ~500 nm promoting higher IFN-? cytokine production from splenocytes and higher IL-1ß at the site of injection.

A New Class of Processes for Formalizing and Generalizing Individual-Based Models: The Semi-Semi-Markov Processes

2000 Mathematics Subject Classification: 60K15, 60K20, 60G20,60J75, 60J80, 60J85, 60-08, 90B15.

Data compilation on the biological response to ocean acidification: environmental and experimental context of data sets and related literature

The exponential growth of studies on the biological response to ocean acidification over the last few decades has generated a large amount of data. To facilitate data comparison, a data compilation hosted at the data publisher PANGAEA was initiated in 2008 and is updated on a regular basis (doi:10.1594/PANGAEA.149999). By January 2015, a total of 581 data sets (over 4 000 000 data points) from 539 papers had been archived. Here we present the developments of this data compilation five years since its first description by Nisumaa et al. (2010). Most of study sites from which data archived are still in the Northern Hemisphere and the number of archived data from studies from the Southern Hemisphere and polar oceans are still relatively low. Data from 60 studies that investigated the response of a mix of organisms or natural communities were all added after 2010, indicating a welcomed shift from the study of individual organisms to communities and ecosystems. The initial imbalance of considerably more data archived on calcification and primary production than on other processes has improved. There is also a clear tendency towards more data archived from multifactorial studies after 2010. For easier and more effective access to ocean acidification data, the ocean acidification community is strongly encouraged to contribute to the data archiving effort, and help develop standard vocabularies describing the variables and define best practices for archiving ocean acidification data.

Seawater carbonate chemistry, calcification and shell size of hard clam Mercenaria mercenaria during experiments, 2011

Increasing atmospheric carbon dioxide threatens to decrease pH in the world's oceans. Coastal and estuarine calcifying organisms of significant ecological and economical importance are at risk; however, several biogeochemical processes drive pH in these habitats. In particular, coastal and estuarine sediments are frequently undersaturated with respect to calcium carbonate due to high rates of organic matter remineralization, even when overlying waters are saturated. As a result, the post-larval stages of infaunal marine bivalves must be able to deposit new shell material in conditions that are corrosive to shell. We measured calcification rates on the hard clam, Mercenaria spp.,in 5 post-larval size classes (0.39, 0.56, 0.78, 0.98, and 2.90 mm shell height) using the alkalinity anomaly method. Acidity of experimental water was controlled by bubbling with air-CO2 blends to obtain pH values of 8.02, 7.64, and 7.41, corresponding to pCO2 values of 424, 1120, and 1950 µatm. These pH values are typical of those found in many near-shore terrigenous marine sediments. Our results show that calcification rate decreased with lower pH in all 5 size classes measured. We also found a significant effect of size on calcification rate, with the smaller post-larval sizes unable to overcome dissolution pressure. Increased calcification rate with size allowed the larger sizes to overcome dissolution pressure and deposit new shell material under corrosive conditions. Size dependency of pH effects on calcification is likely due to organogenesis and developmental shifts in shell mineralogy occurring through the post-larval stage. Furthermore, we found significantly different calcification rates between the 2 sources of hard clams we used for these experiments, most likely due to genotypic differences. Our findings confirm the susceptibility of the early life stages of this important bivalve to decreasing pH and reveal mechanisms behind the increased mortality in post-larval juvenile hard clams related to dissolution pressure, that has been found in previous studies.

Applied evolution: An integrated approach to studying life history traits in response to drug selection

The use of chemical control measures to reduce the impact of parasite and pest species has frequently resulted in the development of resistance. Thus, resistance management has become a key concern in human and veterinary medicine, and in agricultural production. Although it is known that factors such as gene flow between susceptible and resistant populations, drug type, application methods, and costs of resistance can affect the rate of resistance evolution, less is known about the impacts of density-dependent eco-evolutionary processes that could be altered by drug-induced mortality. The overall aim of this thesis was to take an experimental evolution approach to assess how life history traits respond to drug selection, using a free-living dioecious worm (Caenorhabditis remanei) as a model. In Chapter 2, I defined the relationship between C. remanei survival and Ivermectin dose over a range of concentrations, in order to control the intensity of selection used in the selection experiment described in Chapter 4. The dose-response data were also used to appraise curve-fitting methods, using Akaike Information Criterion (AIC) model selection to compare a series of nonlinear models. The type of model fitted to the dose response data had a significant effect on the estimates of LD50 and LD99, suggesting that failure to fit an appropriate model could give misleading estimates of resistance status. In addition, simulated data were used to establish that a potential cost of resistance could be predicted by comparing survival at the upper asymptote of dose-response curves for resistant and susceptible populations, even when differences were as low as 4%. This approach to dose-response modeling ensures that the maximum amount of useful information relating to resistance is gathered in one study. In Chapter 3, I asked how simulations could be used to inform important design choices used in selection experiments. Specifically, I focused on the effects of both within- and between-line variation on estimated power, when detecting small, medium and large effect sizes. Using mixed-effect models on simulated data, I demonstrated that commonly used designs with realistic levels of variation could be underpowered for substantial effect sizes. Thus, use of simulation-based power analysis provides an effective way to avoid under or overpowering a study designs incorporating variation due to random effects. In Chapter 4, I 3 investigated how Ivermectin dosage and changes in population density affect the rate of resistance evolution. I exposed replicate lines of C. remanei to two doses of Ivermectin (high and low) to assess relative survival of lines selected in drug-treated environments compared to untreated controls over 10 generations. Additionally, I maintained lines where mortality was imposed randomly to control for differences in density between drug treatments and to distinguish between the evolutionary consequences of drug treatment versus ecological processes affected by changes in density-dependent feedback. Intriguingly, both drug-selected and random-mortality lines showed an increase in survivorship when challenged with Ivermectin; the magnitude of this increase varied with the intensity of selection and life-history stage. The results suggest that interactions between density-dependent processes and life history may mediate evolved changes in susceptibility to control measures, which could result in misleading conclusions about the evolution of heritable resistance following drug treatment. In Chapter 5, I investigated whether the apparent changes in drug susceptibility found in Chapter 4 were related to evolved changes in life-history of C. remanei populations after selection in drug-treated and random-mortality environments. Rapid passage of lines in the drug-free environment had no effect on the measured life-history traits. In the drug-free environment, adult size and fecundity of drug-selected lines increased compared to the controls but drug selection did not affect lifespan. In the treated environment, drug-selected lines showed increased lifespan and fecundity relative to controls. Adult size of randomly culled lines responded in a similar way to drug-selected lines in the drug-free environment, but no change in fecundity or lifespan was observed in either environment. The results suggest that life histories of nematodes can respond to selection as a result of the application of control measures. Failure to take these responses into account when applying control measures could result in adverse outcomes, such as larger and more fecund parasites, as well as over-estimation of the development of genetically controlled resistance. In conclusion, my thesis shows that there may be a complex relationship between drug selection, density-dependent regulatory processes and life history of populations challenged with control measures. This relationship could have implications for how resistance is monitored and managed if life histories of parasitic species show such eco-evolutionary responses to drug application.

Small-scale genetic structure of Cerastoderma glaucum in a lagoonal environment: potential significance of habitat discontinuity and unstable population dynamics

Environmental heterogeneity in coastal lagoons is expected to facilitate local adaptation in response to different ecological conditions, causing significant genetic structuring within lagoon populations at a small scale and also differentiation between lagoons. However, these patterns and processes of genetic structuring are still poorly understood. The aims of our study were (1) to seek genetic structure at a small scale in Cerastoderma glaucum inside the Mar Menor coastal lagoon using a mitochondrial DNA marker (COI) that has previously detected genetic differentiation inside the lagoon in other species and (2) to evaluate the influence of extreme environmental conditions and habitat discontinuity on its genetic composition. The results indicate high levels of haplotype diversity and low values of nucleotide diversity. COI data provide evidence of significant population differentiation among some localities within the lagoon. Limited gene flow and unstable population dynamics (i.e. fluctuations in population size caused by local extinction and recolonization), probably due to the high environmental heterogeneity, could generate the small-scale genetic divergence detected between populations within the lagoon.

Survival, population pize, and gonotrophic cycle duration of Nyssomyia neivai (Diptera: Psychodidae) at an endemic area of American Cutaneous Leishmaniasis in Southeastern Brazil

The survival, absolute population size, gonotrophic cycle duration, and temporal and spatial abundance of Nyssomyia neivai (Pinto) were studied in a rural area endemic for American cutaneous leishmaniasis (ACL) in Conchal, Sõo Paulo State, southeastern Brazil, using mark-release-recapture techniques and by monitoring population fluctuation. The monthly abundance exhibited a unimodal pattern, with forest and domicile habitats having the highest relative abundances. A total of 1,873 males and 3,557 females were marked and released during the six experiments, of which 4.1-13.0 per cent of males and 4.1-11.8 per cent of females were recaptured. Daily survivorship estimated from the decline in recaptures per day was 0.681 for males and 0.667 for females. Gonotrophic cycle duration was estimated to be 4.0 d. Absolute population size was calculated using the Lincoln Index and ranged from 861 to 4,612 males and from 2,187 to 19,739 females. The low proportion of females that reach the age when they are potentially infective suggests that N. neivai has a low biological capacity to serve as a vector and that factors such as high biting rates and opportunistic feeding behavior would be needed to enable Leishmania (Viannia) braziliensis Vianna transmission. This agreed with the epidemiological pattern of ACL in southeastern Brazil that is characterized by low incidence, with isolated cases acquired principally within domiciliary habitats

POPREP: a generic report for population management

Genetic variation provides a basis upon which populations can be genetically improved. Management of animal genetic resources in order to minimize loss of genetic diversity both within and across breeds has recently received attention at different levels, e. g., breed, national and international levels. A major need for sustainable improvement and conservation programs is accurate estimates of population parameters, such as rate of inbreeding and effective population size. A software system (POPREP) is presented that automatically generates a typeset report. Key parameters for population management, such as age structure, generation interval, variance in family size, rate of inbreeding, and effective population size form the core part of this report. The report includes a default text that describes definition, computation and meaning of the various parameters. The report is summarized in two pdf files, named Population Structure and Pedigree Analysis Reports. In addition, results (e. g., individual inbreeding coefficients, rate of inbreeding and effective population size) are stored in comma-separate-values files that are available for further processing. Pedigree data from eight livestock breeds from different species and countries were used to describe the potential of POPREP and to highlight areas for further research.

Structure of a Population of Hydromedusa maximiliani (Testudines, Chelidae) from Parque Estadual da Serra do Mar, an Atlantic Rainforest Preserve in Southeastern Brazil

We studied the structure of a population of Hydromedusa maximiliani associated with a stream in Parque Estadual da Serra do Mar, Nucleo Itutinga-Piloes, southeastern Brazil, between October 2004 and October 2005. Twenty-five individuals were captured, and a population size of 43.72 +/- 23.7 individuals was estimated. This value is similar to that of the population of Parque Estadual Carlos Botelho, another Atlantic forest reserve of southeastern Brazil. Males were recaptured more frequently than females, suggesting higher activity and/or greater movement of males.

Understanding the genetic diversity, spatial genetic structure and mating system at the hierarchical levels of fruits and individuals of a continuous Theobroma cacao population from the Brazilian Amazon

Understanding the mating patterns of populations of tree species is a key component of ex situ genetic conservation. In this study, we analysed the genetic diversity, spatial genetic structure (SGS) and mating system at the hierarchical levels of fruits and individuals as well as pollen dispersal patterns in a continuous population of Theobroma cacao in Para State, Brazil. A total of 156 individuals in a 0.56 ha plot were mapped and genotyped for nine microsatellite loci. For the mating system analyses, 50 seeds were collected from nine seed trees by sampling five fruits per tree (10 seeds per fruit). Among the 156 individuals, 127 had unique multilocus genotypes, and the remaining were clones. The population was spatially aggregated; it demonstrated a significant SGS up to 15m that could be attributed primarily to the presence of clones. However, the short seed dispersal distance also contributed to this pattern. Population matings occurred mainly via outcrossing, but selfing was observed in some seed trees, which indicated the presence of individual variation for self-incompatibility. The matings were also correlated, especially within ((r) over cap (p(m)) = 0.607) rather than among the fruits ((r) over cap (p(m)) = 0.099), which suggested that a small number of pollen donors fertilised each fruit. The paternity analysis suggested a high proportion of pollen migration (61.3%), although within the plot, most of the pollen dispersal encompassed short distances (28m). The determination of these novel parameters provides the fundamental information required to establish long-term ex situ conservation strategies for this important tropical species. Heredity (2011) 106, 973-985; doi:10.1038/hdy.2010.145; published online 8 December 2010

The population of Brazil

In 1970 the population of Brazil with 94,508,554 inhabitants was extreme youth, since 42.67% was composed of children under 14 years old. In that year the proportion of female was 50.2%. The population density increased from 1.17 inhabitants /km² in 1872 to 11.18 in 1970, and in this last year the range was 1.03 in the North region and 43.90 in the South-East region. The urban population increased from 31.24% in 1940 to 55.98% in 1970 and for the first time the rural population was smaller than the urban population. In 1950 concerning with marital status 39% of the population 15 years old and over was single and 54% married. In 1970 this rate was respectively 35.4% and 56.6%. The population economically inactive increased from 49.17% in 1940 to 52.24% in 1970. The literacy ratio increased from 43% in 1940, to 48% in 1950 and 68.04% in 1970. The crude birth rate was 43/1000 live births in 1950 and fell to 37.7/1000 in 1970. The fertility rate decreased from 179.3/1000 women (15-49 years old} to 156.7/1000 in 1960/70. The crude death rate decreased from 20.60/1000 inhabitants in 1940/50 to 9.4/1000 in 1960/70. The infant mortality rate still remains high: 171/1000 live births in 1940/50 and 170/1000 in 1971. Concerning with the size of the cities, 8 in 1940 had 100,000 or more inhabitants and in 1970 this number increased to 94 cities. The population growth increased from 2.38% in 1940/50, to 2.99% in 1950/60 and 2.83% in 1960./70. Brazil is the first country in population size in Latin America and the eighth in the world. Concerning his area, Brazil is the fifth country in size.

Evaluation of dog population in an urban area of Southeastern Brazil

Planning control programs, for diseases such as rabies requires information on the size and structure of the dog and cat population. In order to evaluate the dog population of the urban area of Araçatuba city, S. Paulo State, Brazil, a survey was conducted using a questionnaire to interview members of households. Eighty-eight districts were visited (37,778 houses) and the interview was possible at 77.93% of these. Human population size evaluated was 113,157 inhabitants. Houses that owned animals represented 55.2%, 26,926 of the animals concerned were dogs and 5,755 were cats. Of the dogs, 56.64% were 1-4 year olds and males represented 56.2% of the total population. Dog: person ratio was estimated at 2.8 dogs to every 10 persons, almost 3 times the ratio hitherto estimated and used in the planning of rabies vaccination campaigns.

Effects of density versus size on grazing by a key exploited herbivore

Fishing decreases the biomass of target species via reduction in the numbers and/or size of individuals. In natural systems, the strength of biological interactions, including predator-prey dynamics, are often density or size-dependent. Hence, changes in the numbers or size of key taxa may be expected to influence biological interactions but their effects do not need to be identical. Here we compare the effects of biomass reduction in populations of the exploited limpet Patella candei. Biomass removal was experimentally achieved by either removing individuals (density reduction) or by replacing large by small individuals (size reduction), while controlling for total limpet biomass in a laboratory-based experiment. At the experiment’s termination, biomass reduction led to proportional changes in area grazed. However, there was no difference whether this was achieved via changes in density or in size. Furthermore, no discernible effects of treatments were evident on different components of the algal assemblage. A field survey also revealed that P. candei biomass explained a greater proportion in variation in the area free of algae than density alone. Our results suggest that loss of biomass in populations of P. candei has quantitatively and qualitatively similar effects on algal cover regardless of whether it is caused by an equivalent (biomass) reduction in the numbers or size of individuals.