Environmental Fluctuations and Level of Density-Compensation Strongly Affects the Probability of Fixation and Fixation Times

The probability of, and time to, fixation of a mutation in a population has traditionally been studied by the classic Wright–Fisher model where population size is constant. Recent theoretical expansions have covered fluctuating populations in various ways but have not incorporated models of how the environment fluctuates in combination with different levels of density-compensation affecting fecundity. We tested the hypothesis that the probability of, and time to, fixation of neutral, advantageous and deleterious mutations is dependent on how the environment fluctuates over time, and on the level of density-compensation. We found that fixation probabilities and times were dependent on the pattern of autocorrelation of carrying capacity over time and interacted with density-compensation. The pattern found was most pronounced at small population sizes. The patterns differed greatly depending on whether the mutation was neutral, advantageous, or disadvantageous. The results indicate that the degree of mismatch between carrying capacity and population size is a key factor, rather than population size per se, and that effective population sizes can be very low also when the census population size is far above the carrying capacity. This study highlights the need for explicit population dynamic models and models for environmental fluctuations for the understanding of the dynamics of genes in populations.     

The Final Size of an Epidemic and Its Relation to the Basic Reproduction Number

We study the final size equation for an epidemic in a subdivided population with general mixing patterns among subgroups. The equation is determined by a matrix with the same spectrum as the next generation matrix and it exhibits a threshold controlled by the common dominant eigenvalue, the basic reproduction number R₀{\mathcal{R}_{0}}: There is a unique positive solution giving the size of the epidemic if and only if R₀{\mathcal{R}_{0}} exceeds unity. When mixing heterogeneities arise only from variation in contact rates and proportionate mixing, the final size of the epidemic in a heterogeneously mixing population is always smaller than that in a homogeneously mixing population with the same basic reproduction number R₀{\mathcal{R}_{0}}. For other mixing patterns, the relation may be reversed.     

Prey Selection,Vertical Migrations and the Impacts of Harvesting Upon the Population Dynamics of a Predator-Prey System

A model is developed to describe the interaction between a predator and two prey types located in different regions. Conditions for stability and persistence are analysed. The effects of harvesting the predators are investigated by making the predator mortality rate habitat dependent. Results demonstrate that for any given set of parameter values there is a value of the intrinsic preference of the predator for each prey type at which the system undergoes a Hopf bifurcation. Above this critical value the system evolves towards a stable equilibrium, whereas below it, stable limit cycles arise by Hopf bifurcations. Simulations demonstrate that the presence of demographic stochasticity may destabilise oscillatory populations, thereby causing population extinctions. An application of the model to the foraging behaviour of North Sea cod is described. It is shown that if the preferred prey is more productive, it is likely that the equilibrium will be stable, whereas if the less preferred prey is more productive, populations are likely to display cycles and in the stochastic case become extinct. As cod fishing mortality is increased, the point of bifurcation and region of parameter space for which the system is unstable decreases. An increased understanding of how cod behave may enable fish stocks to be managed more successfully, for example by indicating where marine reserves should be placed.     

Age and Sex Structured Model for Assessing the Demographic Impact of Mother-to-Child Transmission of HIV/AIDS

Age and sex structured HIV/AIDS model with explicit incubation period is proposed as a system of delay differential equations. The model consists of two age groups that are children (0–14 years) and adults (15–49 years). Thus, the model considers both mother-to-child transmission (MTCT) and heterosexual transmission of HIV in a community. MTCT can occur prenatally, at labour and delivery or postnatally through breastfeeding. In the model, we consider the children age group as a one-sex formulation and divide the adult age group into a two-sex structure consisting of females and males. The important mathematical features of the model are analysed. The disease-free and endemic equilibria are found and their stabilities investigated. We use the Lyapunov functional approach to show the local stability of the endemic equilibrium. Qualitative analysis of the model including positivity and boundedness of solutions, and persistence are also presented. The basic reproductive number (ℛ₀) for the model shows that the adult population is responsible for the spread HIV/AIDS epidemic, thus up-to-date developed HIV/AIDS models to assess intervention strategies have focused much on heterosexual transmission by the adult population and the children population has received little attention. We numerically analyse the HIV/AIDS model to assess the community benefits of using antiretroviral drugs in reducing MTCT and the effects of breastfeeding in settings with high HIV/AIDS prevalence ratio using demographic and epidemiological parameters for Zimbabwe.     

Mathematical Study of the Role of Gametocytes and an Imperfect Vaccine on Malaria Transmission Dynamics

A mathematical model is developed to assess the role of gametocytes (the infectious sexual stage of the malaria parasite) in malaria transmission dynamics in a community. The model is rigorously analysed to gain insights into its dynamical features. It is shown that, in the absence of disease-induced mortality, the model has a globally-asymptotically stable disease-free equilibrium whenever a certain epidemiological threshold, known as the basic reproduction number (denoted by ℛ₀), is less than unity. Further, it has a unique endemic equilibrium if ℛ₀>1. The model is extended to incorporate an imperfect vaccine with some assumed therapeutic characteristics. Theoretical analyses of the model with vaccination show that an imperfect malaria vaccine could have negative or positive impact (in reducing disease burden) depending on whether or not a certain threshold (denoted by ∇) is less than unity. Numerical simulations of the vaccination model show that such an imperfect anti-malaria vaccine (with a modest efficacy and coverage rate) can lead to effective disease control if the reproduction threshold (denoted by ℛ_vac) of the disease is reasonably small. On the other hand, the disease cannot be effectively controlled using such a vaccine if ℛ_vac is high. Finally, it is shown that the average number of days spent in the class of infectious individuals with higher level of gametocyte is critically important to the malaria burden in the community.     

In vivo response to the different LET regions of a 55 MeV helium ion beam

We evaluated the in vivo response to heavy particle irradiation in rat tail epidermis using silver-stained nucleolar organizer regions (AgNOR) as the end-point. The energy degradation of the beam across the circular section of the tail allowed us to study the damage elicited by two different LET regions of a helium beam, i.e. non-Bragg peak (NBP) and Bragg peak (BP), at different sites on the same sample. The tails were locally irradiated with a helium ion beam at different fluences. AgNOR exhibited marked variations between tissue areas only a few micrometers apart within each tail exposed to a given beam fluence. An analysis of the AgNOR variations in NBP and BP areas of tails exposed to different beam fluences revealed a dose-dependent effect. The AgNOR provide quantitative evidence of differential damage in neighboring tissue areas exposed to different LET regions of a helium-ion beam.     

Erratum to: Attractor Dynamics and Thermodynamic Analogies in the Cerebral Cortex: Synchronous Oscillation,the Background EEG,and the Regulation of Attention

Ongoing changes in attention and cognition depend upon cortical/subcortical interactions, which select sequences of different spatial patterns of activation in the cortex.     

The Role of Cohesiveness in the Permeability of the Spatial Assemblies of FG Nucleoporins

Nuclear pore complexes (NPCs) conduct selective, bidirectional transport across the nuclear envelope. The NPC passageway is lined by intrinsically disordered proteins that contain hydrophobic phenylalanine-glycine (FG) motifs, known as FG nucleoporins (FG nups), that play the key role in the NPC transport mechanism. Cohesive interactions among the FG nups, which arise from the combination of hydrophobic, electrostatic, and other forces, have been hypothesized to control the morphology of the assemblies of FG nups in the NPC, as well as their permeability with respect to the transport proteins. However, the role of FG nup cohesiveness is still vigorously debated. Using coarse-grained polymer theory and numerical simulations, we study the effects of cohesiveness on the selective permeability of in vitro FG nup assemblies in different geometries that have served as proxies for the morphological and transport properties of the NPC. We show that in high-density FG nup assemblies, increase in cohesiveness leads to the decrease in their permeability, in accordance with the accepted view. On the other hand, the permeability of low-density assemblies is a nonmonotonic function of the cohesiveness, and a moderate increase in cohesiveness can enhance permeability. The density- and cohesiveness-dependent effects on permeability are explained by considering the free-energy cost associated with penetrating the FG nup assemblies. We discuss the implications of these findings for the organization and function of the NPC.     

Physical mapping and evolution of the centromeric class I gene-containing region of the rat MHC

We physically mapped the centromeric part of the BN rat MHC (RT1n haplotype) in a contig of overlapping P1-derived artificial chromosome (PAC) clones encompassing about 300 kb. The following genes were identified and ordered as: (Syngap, Hset, Daxx, Bing1)-Tapbp-Rgl2-Ke2-Bing4-B3galt4- Rps18-Sacm2l-RT1-A1-RT1-A2-RT1-A3-Ring1-Ring2-++ +Ke4-Rxrb-Col11a2-RT1-Hb-Ring3-RT1-DMb. Thus, in contrast to other RT1 haplotypes, RT1n contains three class I genes, RT1-A1, RT1-A2, and RT1-A3, mapping between the Sacm2l and Ring1 genes. Comparisons of the sequences flanking the Sacm2L and Ring1 genes in rat, human, and mouse suggest that the class I gene-containing region was inserted between these genes in rat and mouse at a similar position. Thus, this insertion is likely to have occurred in a common ancestor of these rodents, although the presence of a site particularly permissive for insertions cannot be excluded.     

Modeling and Estimation of Kinetic Parameters and Replicative Fitness of HIV-1 from Flow-Cytometry-Based Growth Competition Experiments

Growth competition assays have been developed to quantify the relative fitness of HIV-1 mutants. In this article, we develop mathematical models to describe viral/cellular dynamic interactions in the assay system from which the competitive fitness indices or parameters are defined. In our previous HIV-viral fitness experiments, the concentration of uninfected target cells was assumed to be constant (Wu et al. 2006). But this may not be true in some experiments. In addition, dual infection may frequently occur in viral fitness experiments and may not be ignorable. Here, we relax these two assumptions and extend our earlier viral fitness model (Wu et al. 2006). The resulting models then become nonlinear ODE systems for which closed-form solutions are not achievable. In the new model, the viral relative fitness is a function of time since it depends on the target cell concentration. First, we studied the structure identifiability of the nonlinear ODE models. The identifiability analysis showed that all parameters in the proposed models are identifiable from the flow-cytometry-based experimental data that we collected. We then employed a global optimization approach (the differential evolution algorithm) to directly estimate the kinetic parameters as well as the relative fitness index in the nonlinear ODE models using nonlinear least square regression based on the experimental data. Practical identifiability was investigated via Monte Carlo simulations.     

Elite swimmers and the D allele of the ACE I/D polymorphism

A polymorphism of the human angiotensin-1-converting enzyme (ACE) gene has been identified in which the presence (insertion, I allele) of a 287-bp fragment rather than the absence (deletion, D allele) is associated with lower ACE activity. Several recent studies have shown an association of the I allele with endurance performance, it being found with excess frequency in elite distance runners, rowers and mountaineers. Other workers using heterogeneous cohorts of athletes from mixed sporting disciplines have found no such association. An increasing linear trend of I allele frequency with the distance run amongst Olympic runners and an excess of the D allele amongst sprinters led us to examine whether the ratio of I and D alleles in swimmers competing over different distances would also vary. Swimmers (n=120) from the European and Commonwealth championships and an American college team had their ACE genotype determined and their gene and allele frequencies compared with several control groups, the most closely age-matched of which were 1,248 military recruits. Of the 103 Caucasians, there was a significant excess of the D allele compared with this control group only in the truly elite swimmers of the European and Commonwealth championships (P=0.004). This association remained in those competing over shorter distances (P=0.005 for 400 m and below) but not in the longer events. These findings were confirmed in three further large control groups. A population association study testing whether a genetic marker (the ACE I/D polymorphism) occurs more frequently in cases (elite athletes) than in controls therefore requires a homogeneous cohort of subjects from the same sporting discipline.     

Emergence and Spread of Basal Lineages of Yersinia pestis during the Neolithic Decline

1. Download high-res image (316KB)
2. Download full-size image

     

Involvement of sulfoquinovosyl diacylglycerol in the structural integrity and heat-tolerance of photosystem II

To examine the role of sulfoquinovosyl diacylglycerol (SQDG) in thylakoid membranes, we compared the structural and functional properties of photosystem II (PSII) between a mutant of Chlamydomonas reinhardtii defective in SQDG ( hf-2) and the wild type. The PSII core complex of hf-2, as compared with that of the wild type, showed structural fragility when solubilized with a detergent, dodecyl beta- d-maltoside, suggesting that the physical properties of the PSII complex were altered by the loss of SQDG. On the other hand, exposure of the cells to 41 degrees C for 120 min in the dark decreased the PSII activity to 70% and 50% of the initial levels in the wild type and hf-2, respectively, which implies that the PSII activity, in the absence of SQDG, becomes less stable under heat-stress conditions. PSII inactivated to 60% of the initial level by dark incubation at 41 degrees C was reactivated by following illumination even at 41 degrees C to more than 90% in the wild type, but only to 70% in hf-2. These results suggest that PSII inactivated by heat recovers through some mechanism dependent on light, and that SQDG participates in functioning of the mechanism. The conformational disorder of PSII caused by the defect in SQDG might be correlated with the increased susceptibility of its activity to heat-stress.     

Mucolipidosis type IV and the mucolipins

MLIV (mucolipidosis type?IV) is a neurodegenerative lysosomal storage disorder caused by mutations in MCOLN1, a gene that encodes TRPML1 (mucolipin-1), a member of the TRPML (transient receptor potential mucolipin) cation channels. Two additional homologues are TRPML2 and TRPML3 comprising the TRPML subgroup in the TRP superfamily. The three proteins play apparently key roles along the endocytosis process, and thus their cellular localization varies among the different group members. Thus TRPML1 is localized exclusively to late endosomes and lysosomes, TRPML2 is primarily located in the recycling clathrin-independent GPI (glycosylphosphatidylinositol)-anchored proteins and early endosomes, and TRPML3 is primarily located in early endosomes. Apparently, all three proteins' main physiological function underlies Ca(2+) channelling, regulating the endocytosis process. Recent findings also indicate that the three TRPML proteins form heteromeric complexes at least in some of their cellular content. The physiological role of these complexes in lysosomal function remains to be elucidated, as well as their effect on the pathophysiology of MLIV. Another open question is whether any one of the TRPMLs bears additional function in channel activity.     

Cloning and expression of Penicillium minioluteum dextranase in Saccharomyces cerevisiae and its exploitation as a reporter in the detection of mycotoxins

A dextranase gene from Penicillium minioluteum (strain IMI068219) has been cloned, sequenced and expressed in Saccharomyces cerevisiae via fusion of the DNA segment encoding the mature dextranase protein with α-factor signal sequence, and insertion into the GAL1–controlled expression vector pYES2/CT. Galactose-induced expression yielded extracellular dextranase activity of 0.63 units/ml and cell-associated dextranase activity of 0.48 units/ml, after 24 h incubation. The dextranase construct was introduced into a strain of S. cerevisiae expressing the human cytochrome P450 3A4 (CYP3A4) and the cognate reductase, which was then used to develop a microplate toxicity bioassay. Toxicity was signalled as inhibition of dextranase activity, assayed fluorimetrically. This novel bioassay was assessed using six economically significant mycotoxins.     

Resonance assignments of the 56 kDa chimeric avidin in the biotin-bound and free forms

Avidin is a homotetrameric ~56 kDa protein found in chicken egg white. Avidin’s ability to bind biotin with a very high affinity has widely been exploited in biotechnological applications. Protein engineering has further diversified avidin’s feasibility. ChiAVD(I117Y) is a product of rational protein engineering. It is a hyperthermostable synthetic hybrid of avidin and avidin-related protein 4 (AVR4). In this chimeric protein a 23-residue segment in avidin has been replaced with the corresponding sequence found in AVR4, and a point mutation at subunit interface 1–3 (and 2–4) has been introduced. Here we report the backbone and sidechain resonance assignments of the biotin-bound form of ChiAVD(I117Y) as well as the backbone resonance assignments of the free form.     

Density functional theory studies on the generation of trimethylenemethanes from the ring opening of dialkoxymethylenecyclopropanes and methylenecyclopropanethioacetals and follow-up reactions

This work reports a detailed DFT study on the generation of trimethylenemethanes (TMMs) from the ring opening of dialkoxymethylenecyclopropane (DMCP), methylenecyclopropanethioacetal (MCPT), and substituted derivatives of DMCP and MCPT, as well as follow-up reactions of the TMMs. The singlet DMCP and MCPT were found to be 51.32 and 53.77 kcal mol^?1 more stable than the triplet DMCP and MCPT respectively, corresponding to triplet:singlet population ratios of 1:10³⁸ and 1:10⁴⁰, respectively, at 25 °C using Boltzmann distribution, implying that the proportion of the triplet species is negligible at 25 °C. The ring-opening reactions occur through singlet transition states with barriers of 40.68 and 42.27 kcal mol^?1 for DMCP and MCPT, respectively, and yield TMMs that are very unstable compared to the precursors, with the triplet TMM being far more stable than the singlet. Whereas the singlet TMMs readily undergo cycloaddition reactions with olefins to form five-membered carbocyclic rings, the triplet species do not. The selectivity of the reactions of the DMCP TMMs is very sensitive to temperature; at 25°C, cycloaddition with olefins and ring-closure to form ketenes have very comparable barriers while temperatures above 150 °C favor the exclusive formation of a ketene followed by dimerization. In MCPT, ring closure to form ketenes is the favored reaction at all temperatures studied.

Open image in new window

Graphical abstract

Pathways for the generation of trimethylenemethanes from their precursors and follow-up reactions