The Evolution of Virulence in RNA Viruses under a Competition–Colonization Trade-Off

RNA viruses exist in large intra-host populations which display great genotypic and phenotypic diversity. We analyze a model of viral competition between two viruses infecting a constantly replenished cell pool. We assume a trade-off between the ability of the virus to colonize new cells (cell killing rate or virulence) and its local competitiveness (replicative success within coinfected cells). We characterize the conditions that allow for viral spread by means of the basic reproductive number and show that a local coexistence equilibrium exists, which is asymptotically stable. At this equilibrium, the less virulent competitor has a reproductive advantage over the more virulent colonizer reflected by a larger equilibrium population size of the competitor. The equilibria at which one virus outcompetes the other one are unstable, i.e., a second virus is always able to permanently invade. We generalize the two-virus model to multiple viral strains, each displaying a different virulence. To account for the large phenotypic diversity in viral populations, we consider a continuous spectrum of virulences and present a continuum limit of this multiple viral strains model that describes the time evolution of an initial continuous distribution of virulence without mutations. We provide a proof of the existence of solutions of the model equations, analytically assess the properties of stationary solutions, and present numerical approximations of solutions for different initial distributions. Our simulations suggest that initial continuous distributions of virulence evolve toward a distribution that is extremely skewed in favor of competitors. At equilibrium, only the least virulent part of the population survives. The discrepancy of this finding in the continuum limit with the two-virus model is attributed to the skewed equilibrium subpopulation sizes and to the transition to a continuum. Consequently, in viral quasispecies with high virulence diversity, the model predicts collective virulence attenuation. This result may contribute to understanding virulence attenuation, which has been reported in several experimental studies.     

Force–Velocity Curves of Motor Proteins Cooperating In Vivo

Motor proteins convert chemical energy into work, thereby generating persistent motion of cellular and subcellular objects. The velocities of motor proteins as a function of opposing loads have been previously determined in vitro for single motors. These single molecule “force–velocity curves” have been useful for elucidating motor kinetics and for estimating motor performance under physiological loads due to, for example, the cytoplasmic drag force on transported organelles. Here we report force–velocity curves for single and multiple motors measured in vivo. Using motion enhanced differential interference contrast (MEDIC) movies of living NT2 (neuron-committed teratocarcinoma) cells at 37°C, three parameters were measured—velocity (v), radius (a), and effective cytoplasmic viscosity (η′)—as they applied to moving vesicles. These parameters were combined in Stokes’ equation, F = 6πaη′v, to determine the force, F, required to transport a single intracellular particle at velocity, v. In addition, the number of active motors was inferred from the multimodal pattern seen in a normalized velocity histogram. Using this inference, the resulting in vivo force–velocity curve for a single motor agrees with previously reported in vitro single motor force–velocity curves. Interestingly, however, the curves for two and three motors lie significantly higher in both measured velocity and computed force, which suggests that motors can work cooperatively to attain higher transport forces and velocities. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Tree fine root Ca/Al molar ratio – Indicator of Al and acidity stress

Many Aloe species are exploited as natural products. Generally, the leaves are unsustainably picked from wild plants to meet the market demand. Basic scientific information on seed biology and the ways of increasing levels of secondary metabolites in seedlings is still lacking for Aloe species. This study investigated seed germination requirements and evaluated levels of secondary metabolites in seedlings of Aloe arborescens, an important species in traditional medicine. The highest percentage germination (78%) and the fastest germination rate (GR) (10% d^? 1) with a mean germination time (MGT) of 9 days were achieved at 20°C under a 16-h photoperiod. At 25°C, maximum percentage germination (67%) (P < 0.05), higher GR (13% d^? 1) and shorter MGT (6 days) were obtained under constant light. These results indicate that temperature and light play a significant role in germination of A. arborescens seeds. Increasing osmotic pressure on seeds decreased percentage germination, whereas buffering the solution to a range of pH values (4–10) did not significantly affect germination. Smoke–water (1:500 v/v), smoke-isolated karrikinolide (10^? 8 and 10^? 9 M) and potassium nitrate (10^? 3 and 10^? 4 M) significantly promoted germination compared with the control at 25°C (supra-optimal temperature) under a 16-h photoperiod. These treatments were also effective in increasing secondary metabolite levels (flavonoids and phenolics) in A. arborescens seedlings.     

Conditions for Global Dynamic Stability of a Class of Resource-Bounded Model Ecosystems

This paper studies a class of dynamical systems that model multi-species ecosystems. These systems are ‘resource bounded’ in the sense that species compete to utilize an underlying limiting resource or substrate. This boundedness means that the relevant state space can be reduced to a simplex, with coordinates representing the proportions of substrate utilized by the various species. If the vector field is inward pointing on the boundary of the simplex, the state space is forward invariant under the system flow, a requirement that can be interpreted as the presence of non-zero exogenous recruitment. We consider conditions under which these model systems have a unique interior equilibrium that is globally asymptotically stable. The systems we consider generalize classical multi-species Lotka–Volterra systems, the behaviour of which is characterized by properties of the community (or interaction) matrix. However, the more general systems considered here are not characterized by a single matrix, but rather a family of matrices. We develop a set of ‘explicit conditions’ on the basis of a notion of ‘uniform diagonal dominance’ for such a family of matrices, that allows us to extract a set of sufficient conditions for global asymptotic stability based on properties of a single, derived matrix. Examples of these explicit conditions are discussed.     

Searching for Spatial Patterns in a Pollinator–Plant–Herbivore Mathematical Model

This paper deals with the spatio-temporal dynamics of a pollinator–plant–herbivore mathematical model. The full model consists of three nonlinear reaction–diffusion–advection equations defined on a rectangular region. In view of analyzing the full model, we firstly consider the temporal dynamics of three homogeneous cases. The first one is a model for a mutualistic interaction (pollinator–plant), later on a sort of predator–prey (plant–herbivore) interaction model is studied. In both cases, the interaction term is described by a Holling response of type II. Finally, by considering that the plant population is the unique feeding source for the herbivores, a mathematical model for the three interacting populations is considered. By incorporating a constant diffusion term into the equations for the pollinators and herbivores, we numerically study the spatiotemporal dynamics of the first two mentioned models. For the full model, a constant diffusion and advection terms are included in the equation for the pollinators. For the resulting model, we sketch the proof of the existence, positiveness, and boundedness of solution for an initial and boundary values problem. In order to see the separated effect of the diffusion and advection terms on the final population distributions, a set of numerical simulations are included. We used homogeneous Dirichlet and Neumann boundary conditions.     

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.     

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.     

Serological identification and molecular characterization of B (A) 02 subtype in patients and blood donors from Eastern China

This study was designed to identify the rare?type?ABO?blood?groups, B(A) 02, from Eastern China. Three samples with discordant serological results during routine blood type identification and four samples from one sample’ family were selected. All of them were detected by serological method. The exon 6 and 7 of the ABO genes were amplified by PCR and sequenced. They were typed as AsubB by serology and as BO by genotype. In AsubB samples, nt 700C>G mutation was detected in B gene, which was previously defined as B(A)02 alleles. In these seven samples, six showed B(A)02/O01 and one showed B(A)02/O02.B(A)02 allele was found to be more common in this study than B(A)04 which is considered to be more frequent than B(A)02. The careful identification of rare blood types is important for the safety of clinical blood transfusion.     

Genetic diversity and structure of Lilium pumilum DC. in southeast of Qinghai–Tibet plateau

Lilium pumilum DC. is a valuable species not only for its showy flowers but also for its edible and medicinal values. As one of the distribution areas of L. pumilum, Qinghai–Tibet plateau has unique environmental features which have high impact on the evolution of the species. No population genetic studies have been done for L. pumilum so far. To provide the first reference data for evolutionary study and understanding the influence of eco-geographic factors on the distribution of genetic variation in L. pumilum, interspecific simple sequence repeat markers were used to investigate genetic diversity and population structure of 28 populations sampled from southeast of Qinghai–Tibet plateau. Fifteen selected primers generated a total of 147 polymorphic bands. The genetic diversity was low within populations (average He = 0.173), but higher at the species level (He = 0.392). A clear population structure and high level of genetic differentiation (F _ST = 0.518) were detected by unweighted pair group method for arithmetic averages, principle coordinate analysis and Bayesian clustering. All clustering approaches supported a division of the 28 populations into 4 major groups for which analysis of molecular variance confirmed a significant variation among groups (34.3 %). These population genetic parameters suggest limited gene flow among populations and evidence for isolation by distance (r = 0.272, P < 0.0001) was found in this study. Altitude, AMT and AMP explained 9.5, 11.5 and 14.0 % of the total variance among populations indicating that eco-geographic factors have a significant effect. Considering the low within-population genetic diversity, high differentiation among populations and the increasing anthropogenic pressure on the species, in situ conservation measures were recommended to preserve L. pumilum in Qinghai–Tibet plateau.     

Antimicrobial Activities of Bacteriocins E 50–52 and B 602 Against Antibiotic-Resistant Strains Involved in Nosocomial Infections

The antimicrobial spectra of previously published bacteriocin E 50–52 (39 a.a.; 3,932 Da; pI = 8.5) and bacteriocin B 602 (29 a.a.; 3,864 Da; pI = 7.2) were determined. Named peptides were related to class IIa (pediocin-like) bacteriocins. Minimal inhibitory concentrations (MICs) of bacteriocins have been determined for bacterial isolates that were causative agents of nosocomial infections collected from Russian hospitals in 2003–2007, namely methicillin-resistant Staphylococcus aureus (MRSA) (n = 10); Acinetobacter baumannii (n = 11); Citrobacter freundii (n = 8); Escherichia coli (n = 9); Klebsiella pneumoniae (n = 10); Proteus spp. (n = 6); and Pseudomonas aeruginosa (n = 10). The majority of these tested isolates have been shown to be multidrug resistant and carry genetic determinants of antimicrobial resistance that were detected using polymerase chain reaction (PCR). The MICs of bacteriocin B 602 ranged from ≤0.025–1.56 μg/ml, and for bacteriocin E 50–52 from 0.05 to 6.25 μg/ml for all of 64 bacterial clinical isolates tested. Interestingly, the bacteriocins studied demonstrate activity on both Gram-positive and Gram-negative bacteria. Bacteriocins E 50–52 and B 602 show good activity against nosocomial bacterial agents resistant to many classes of modern antibacterials used in clinical practice. These bacteriocins should be examined as an alternative in treating infections caused by such agents.     

Structural investigation of human wild-type Cu,Zn superoxide dismutase – a biomolecular case study using DL_POLY_3

This article demonstrates the capability of DL_POLY_3 to study protein molecules at long time scales. The example shown here is the study of early stages of structural destabilisation of human superoxide dismutase (SOD1) as a result of metal depletions. Simulations of up to 4 ns have been carried out on wild-type SOD1. It is shown that the removal of metal ions from the active site essentially destabilises loop structures and provides a potential precursor for unintended protein–protein interactions to take place, as observed experimentally. The results reveal early stages of the destabilisation of structures that may enhance subsequent protein aggregations. It is found that the removal of Zn ion at the active sites of the protein enzyme causes a disruption to the Zn-loop and subsequently the electrostatic-loop (e-loops) with an increase in flexibility and mobility of the structures. The Zn ion is important to maintain structural integrity of the channel loops, whereas the Cu ions apparently play a lesser role in this aspect.     

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.     

Unravelling the past 1,000 years of history of human–climate–landscape interactions at the Lindu plain,Sulawesi, Indonesia

The Lindu plain, located in the northern mountainous region of the Lore Lindu National Park in Sulawesi, Indonesia, provides many ecosystem services for the population inhabiting the area and harbours a unique biodiversity. Palynological, charcoal and diatom analyses of a lake sediment core from Lake Lindu (Danau Lindu) reveal that during the last 1,000 years the Lindu plain has been modified by human activities. Evidence of frequent burning and possible shifting cultivation from an earlier phase from ca. ad 1000 to 1200 might be related to the metal age population which erected the megaliths in the province of Central Sulawesi. From ca. ad 1200–1700 there followed 500 years of wetter climate conditions, corresponding to the southward movement of the Intertropical Convergence Zone. At the same time, decreases of macro-charcoal concentrations and pioneer vegetation indicators show that the use of the landscape of Lindu plain had become more permanent. Following a phase of forest recovery from ca. ad 1730 to 1910, the most recent part of the Lake Lindu record shows a trend towards deforestation that started in the late 20th century, lasting until now. The lake level started to fall at the beginning of the 20th century, as shown by the increase of sedimentation rate and supported by low pollen concentration and palaeomagnetic data. Such a change was unprecedented for the last 1,000 years covered by the record, and it has no link to the climate variability as reconstructed for the last hundred years. If deforestation increases and a larger amount of water is channelled away from the lake for irrigation purposes, the lake level will continue to fall. This suggests that there is a need for better management of the forests surrounding the plain and of the irrigation systems in the area open for cultivation.     

Bistability Analysis of an Apoptosis Model in the Presence of Nitric Oxide

Bistability in apoptosis, or programmed cell death, is crucial for the healthy functioning of multicellular organisms. The aim in this study is to show the presence of bistability in a mitochondria-dependent apoptosis model under nitric oxide effects using chemical reaction network theory. The model equations are a set of coupled ordinary differential equations arising from the assumed mass action kinetics. Whether these equations have a capacity for bistability (cell survival and apoptosis) is determined using a modular approach in which the model is decomposed into modules. Each module contains only a subset of the whole model and is analyzed separately. It is seen that bistability in a module is preserved throughout the whole model after adding the remaining reactions in the pathway on these modules. It is also found that inhibitor effect of some proteins and the appearance of a reacting protein in a later stage as a product is a desired feature but not sufficient for bistability (in the absence of cooperativity effects). On the whole model, two apoptotic and two cell survival states are obtained depending on the initial cell conditions. The results suggest that the antiapoptotic effects of nitric oxide species are responsible for the bistable character of the apoptotic pathway when cooperativity is not assumed in the apoptosome formation.     

13C–13C NOESY spectra of a 480 kDa protein: solution NMR of ferritin

Molecular size has limited solution NMR analyses of proteins. We report ¹³C–¹³C NOESY experiments on a 480 kDa protein, the multi-subunit ferritin nanocage with gated pores. By exploiting ¹³C-resonance-specific chemical shifts and spin diffusion effects, we identified 75% of the amino acids, with intraresidue C–C connectivities between nuclei separated by 1–4 bonds. These results show the potential of ¹³C–¹³C NOESY for solution studies of molecular assemblies >100 kDa.     

The life – history of the peach silver mite,Aculus fockeui (Acari: Eriophyidae) in Egypt

The peach silver mite Aculus fockeui (Nalepa and Trouessart) successfully developed from egg to adult stage when reared on a nectarine or peach leaflet attached with soft lateral divided branch dipped in test tube at different constant temperatures and 70% r.h. The effect of temperature on the development, reproduction and population growth was investigated. Below 20°C all activity ceased and by 33°C the adult began to slow down and cease all activity. At least of 34% of the generation time was spent in the egg stage at 29°C. Adult longevity decreased with increasing temperature. Fecundity was highest at 29°C with 43.16 eggs per female. Life table parameters showed that the population of A. fockeui on nectarine leaves multiplied 30.04 times in a generation time of 18.49 days at 29°C, while the population on peach leaves increased 13.87 times in a generation time of 19.18 days under the same conditions.