Genetic diversity of microsatellite loci in hierarchically structured populations

Microsatellite loci are widely used for investigating patterns of genetic variation within and among populations. Those patterns are in turn determined by population sizes, migration rates, and mutation rates. We provide exact expressions for the first two moments of the allele frequency distribution in a stochastic model appropriate for studying microsatellite evolution with migration, mutation, and drift under the assumption that the range of allele sizes is bounded. Using these results, we study the behavior of several measures related to Wright’s F_ST, including Slatkin’s R_ST. Our analytical approximations for F_ST and R_ST show that familiar relationships between N_em and F_ST or R_ST hold when the migration and mutation rates are small. Using the exact expressions for F_ST and R_ST, our numerical results show that, when the migration and mutation rates are large, these relationships no longer hold. Our numerical results also show that the diversity measures most closely related to F_ST depend on mutation rates, mutational models (stepwise versus two-phase), migration rates, and population sizes. Surprisingly, R_ST is relatively insensitive to the mutation rates and mutational models. The differing behaviors of R_ST and F_ST suggest that properties of the among-population distribution of allele frequencies may allow the roles of mutation and migration in producing patterns of diversity to be distinguished, a topic of continuing investigation.     

Evolutionary stability of first-order-information indirect reciprocity in sizable groups

Indirect reciprocity is considered as a key mechanism for explaining the evolution of cooperation in situations where the same individuals interact only a few times. Under indirect reciprocity, an individual who helps others gets returns indirectly from others who know her good reputation. Recently, many studies have discussed the effect of reputation criteria based only on the former actions of the others (first-order information) and of those based also on the former reputation of opponents of the others (second-order information) on the evolution of indirect reciprocity. In this study, we investigate the evolutionary stability of the indirectly reciprocal strategy (discriminating strategy: DIS), which cooperates only with opponents who have good reputations, in -person games where more than two individuals take part in a single group (interaction). We show that in n-person games, DIS is an evolutionarily stable strategy (ESS) even under the image-scoring reputation criterion, which is based only on first-order information and where cooperations (defections) are judged to be good (bad). This result is in contrast to that of 2-person games where DIS is not an ESS under reputation criteria based only on first-order information.     

The shared reward dilemma

One of the most direct human mechanisms of promoting cooperation is rewarding it. We study the effect of sharing a reward among cooperators in the most stringent form of social dilemma, namely the prisoner's dilemma (PD). Specifically, for a group of players that collect payoffs by playing a pairwise PD game with their partners, we consider an external entity that distributes a fixed reward equally among all cooperators. Thus, individuals confront a new dilemma: on the one hand, they may be inclined to choose the shared reward despite the possibility of being exploited by defectors; on the other hand, if too many players do that, cooperators will obtain a poor reward and defectors will outperform them. By appropriately tuning the amount to be shared a vast variety of scenarios arises, including the traditional ones in the study of cooperation as well as more complex situations where unexpected behavior can occur. We provide a complete classification of the equilibria of the n-player game as well as of its evolutionary dynamics.     

The role of low avidity T cells in the protection against type 1 diabetes: a modeling investigation

Cytotoxic T lymphocytes (CTLs) play a dominant role in the pathogenesis of autoimmune diabetes, commonly denoted Type 1 Diabetes (T1D). These CTLs (notably CD8⁺ T cells) recognize and kill insulin-secreting pancreatic β cells, reducing their number by ∼90%. The resulting reduction of insulin secretion causes the defective regulation of glucose metabolism, leading to the characteristic symptoms of diabetes. Recognition of β cells as targets by CTLs depends on the interactions between MHC-peptide complexes on the surface of β cells and receptors (TCRs) on T cells. Those CTLs with high affinity TCRs (also called high avidity T cells) cause most of the harm, while those with low affinity TCRs (also called low avidity T cells) play a more mysterious role. Recent experimental evidence suggests that low avidity T cells accumulate as memory T cells during the disease and may be protective in NOD mice (a strain prone to developing T1D), delaying disease progression. It has been hypothesized that such low avidity T cells afford disease protection either by crowding the islets of Langerhans, where β cells reside, or by killing antigen presenting cells (APCs).In this paper, we explore the hypothesized mechanisms for this protective effect in the context of a series of models for (1) the interactions of low and high avidity T cells, (2) the effect of APCs and (3) the feedback from β cell killing to autoantigen-induced T cell proliferation. We analyze properties of these models, noting consistency of predictions with observed behaviour. We then use the models to examine the influence of various treatment strategies on the progression of the disease. The model reveals that progressive accumulation of memory low avidity autoreactive T cells during disease progression makes treatments aimed at expanding these protective T cell types more effective close to, or at the onset of clinical disease. It also provides evidence for the hypothesis that low avidity T cells kill APCs (rather than the alternate hypothesis that they crowd the islets).     

Transmission dynamics of Toxoplasma gondii along an urban-rural gradient

Recently, several authors have proposed that the availability of intermediate hosts (IHs) for definitive hosts (DHs) may contribute to determining the dynamics and evolutionary ecology of parasites with facultative complex life cycles. The protozoa Toxoplasma gondii may be transmitted to DHs either via predation of infected IHs through a complex life cycle (CLC) or directly from a contaminated environment through a simple life cycle (SLC). This parasite is also present in contrasting host density environments. We tested the hypothesis that the relative contributions of the CLC and SLC along an urban-rural gradient depend on the IH supply. We built and analysed a deterministic model of the T. gondii transmission cycle. The SLC relative contribution is important only in urban-type environments, i.e., with low predation rate on IHs. In contrast, the parasite is predominantly transmitted through a CLC in suburban and rural environments. The association of the two cycles enables the parasite to spread in situations of low IH availability and low DH population size for which each cycle alone is insufficient.     

Melanoblast proliferation dynamics during mouse embryonic development. Modeling and validation

In this paper, we are looking for mathematical modeling of mouse embryonic melanoblast proliferation dynamics, taking into account, the expression level of β‐catenin. This protein plays an important role into the whole signal pathway process. Different assumptions on some unobservable features lead to different candidate models. From real data measured, from biological experiments and from a priori biological knowledge, it was able to validate or invalidate some of the candidate models. Data assimilation and parameter identification allowed us to derive a mathematical model that is in very good agreement with biological data. As a result, the produced model can give tracks for biologists into their biological investigations and experimental evidence. Another interest is the use of this model for robust hidden parameter identification like double times or number of founder melanoblasts.     

Modeling species-specific diacylglycerol dynamics in the RAW 264.7 macrophage

A mathematical model of the G protein signaling pathway in RAW 264.7 macrophages downstream of P2Y₆ receptors activated by the ubiquitous signaling nucleotide uridine 5’-diphosphate is developed. The model, which is based on time-course measurements of inositol trisphosphate, cytosolic calcium, and diacylglycerol, focuses particularly on differential dynamics of multiple chemical species of diacylglycerol. When using the canonical pathway representation, the model predicted that key interactions were missing from the current network structure. Indeed, the model suggested that accurate depiction of experimental observations required an additional branch to the signaling pathway. An intracellular pool of diacylglycerol is immediately phosphorylated upon stimulation of an extracellular receptor for uridine 5’-diphosphate and subsequently used to aid replenishment of phosphatidylinositol. As a result of sensitivity analysis of the model parameters, key predictions can be made regarding which of these parameters are the most sensitive to perturbations and are therefore most responsible for output uncertainty.     

Uncovering operational interactions in genetic networks using asynchronous Boolean dynamics

Biological networks of large dimensions, with their diagram of interactions, are often well represented by a Boolean model with a family of logical rules. The state space of a Boolean model is finite, and its asynchronous dynamics are fully described by a transition graph in the state space. In this context, a model reduction method will be developed for identifying the active or operational interactions responsible for a given dynamic behaviour. The first step in this procedure is the decomposition of the asynchronous transition graph into its strongly connected components, to obtain a “reduced” and hierarchically organized graph of transitions. The second step consists of the identification of a partial graph of interactions and a sub-family of logical rules that remain operational in a given region of the state space. This model reduction method and its usefulness are illustrated by an application to a model of programmed cell death. The method identifies two mechanisms used by the cell to respond to death-receptor stimulation and decide between the survival and apoptotic pathways.     

Population differentiation and migration: coalescence times in a two-sex island model for autosomal and X-linked loci

Evolutionists have debated whether population-genetic parameters, such as effective population size and migration rate, differ between males and females. In humans, most analyses of this problem have focused on the Y chromosome and the mitochondrial genome, while the X chromosome has largely been omitted from the discussion. Past studies have compared F_ST values for the Y chromosome and mitochondrion under a model with migration rates that differ between the sexes but with equal male and female population sizes. In this study we investigate rates of coalescence for X-linked and autosomal lineages in an island model with different population sizes and migration rates for males and females, obtaining the mean time to coalescence for pairs of lineages from the same deme and for pairs of lineages from different demes. We apply our results to microsatellite data from the Human Genome Diversity Panel, and we examine the male and female migration rates implied by observed F_ST values.     

Effects of agronomic application of olive mill wastewater in a field of olive trees on carbohydrate profiles, chlorophyll a fluorescence and mineral nutrient content

Olive mill wastewater (OMW) management is a serious environmental issue for the Mediterranean area where there is the most production of olive oil. OMW contains a high organic load, substantial amounts of plant nutrients but also several compounds with recognized toxicity towards living organisms. Moreover, OMW may represent a low cost source of water. We studied the influence of irrigation with OMW (amounts applied: 30, 60, 100 and 150 m³ h⁻¹) in a field of olive trees on root colonization, photosynthesis, chlorophyll fluorescence, leaf nutrient concentration and soluble carbohydrate. The soil fatty acid methyl ester (FAME) 16:1ω5 was used to quantify biomass of arbuscular mycorrhizal (AM) fungi and the root FAME 16:1ω5 analysis was used as index for the development of colonization in the roots. Agronomic application of OMW decreased significantly the abundance of the soil FAME 16:1ω5 and the root FAME 16:1ω5 in the soil amended with 60, 100 and 150 m³ ha⁻¹ OMW. Decreased root FAME 16:1ω5 due to OMW amendment was associated with a significant reduction of tissue nutrient concentrations in the olive trees. The highest application of OMW to the soil reduced significantly the olive trees uptake of N, P, K, Ca, Mg, Fe, Cu, Mn and Zn. Land spreading of OMW increased concentration of soluble carbohydrate in the olive leaves, mostly due to decreased sink demand for carbon by the root. In the olive trees amended with 150 m³ ha⁻¹ OMW, net CO₂ uptake rate (A), quantum yield of photosystem II electron transport (Φ_PSII), maximal photochemical efficiency of photosystem II (Fv/Fm), photochemical quenching (q_p) and the electron transport rate (ETR) were significantly depressed, whereas non-photochemical quenching (NPQ) was found to increase. Taken with data from experiments in field conditions, our results suggest that agronomic application of OMW alters the functioning of arbuscular mycorrhizas and can even disrupt the relationship between AM fungi and olive trees.     

Compressible gas gills of diving insects: Measurements and models

Many diving insects collect a bubble of air from the surface to supply their oxygen requirements while submerged. It has been theorised that these air bubbles may also act as compressible gas gills, as the low oxygen partial pressure (PO₂) within the bubble caused by the insect's respiration creates a gradient capable of driving the diffusion of oxygen from the water into the bubble. Under these conditions nitrogen diffuses in the opposite direction, resulting in a situation where the volume of the bubble is continually shrinking while oxygen is obtained. This study measures changes in volume and PO₂ within the gas gills held by a tethered water bug, Agraptocorixa eurynome. Both gill volume and PO₂ drop rapidly at the beginning of a dive, but eventually the PO₂ reaches an apparently stable level while volume continually declines at a slower rate. Active ventilation of the gill is crucial to maintaining oxygen uptake. These measurements are used to calculate oxygen flux into the gas gill and the oxygen consumption rate of the bug. The effectiveness of a gas gill as a respiratory organ is also demonstrated by determining the critical PO₂ of the water bug and comparing this with measured gas gill PO₂ and calculated .     

The naturally split Npu DnaE intein exhibits an extraordinarily high rate in the protein trans-splicing reaction

We have studied the naturally split α subunit of the DNA polymerase III (DnaE) intein from Nostoc punctiforme PCC73102 (Npu) using purified proteins and determined an apparent first-order rate constant of (1.1±0.2)×10^-2 s⁻¹ at 37 °C. This represents the highest rate reported for the protein trans-splicing reaction so far (t_1/2 of ∼ 60 s). Furthermore, the reaction was very robust and high-yielding with respect to different extein sequences, temperatures from 6 to 37 °C, and the presence of up to 6 M urea. Given these outstanding properties, the Npu DnaE intein appears to be the intein of choice for many applications in protein and cellular chemistry.     

KatG from Salmonella Typhimurium is a peroxynitritase

Pathogenic bacteria elicit protective responses to oxidative and nitrosative stresses. Although such responses are generally distinct, it was recently reported in Mycobacterium tuberculosis that catalase-peroxidase (KatG), a classical defence against peroxides, also exhibits peroxynitritase activity. Here, the katG gene from Salmonella Typhimurium was cloned and protein purified and characterised. An increase in the rate of decomposition of peroxynitrite was observed for KatG from the enterobacterium with a second-order rate constant of 4.2 × 10⁴ M⁻¹ s⁻¹ at pH 7.4, 25 °C. This enzyme was able to reduce dihydrorhodamine oxidation by peroxynitrite to ∼83%. Given the peroxynitritase activity demonstrated here it is likely that KatG may play a wider role in the detoxification of oxidative stresses than previously thought.     

E versus Z geometry in β-d-arabino-hexopyranosidulose oximes

Koenigs-Knorr-type glycosidations of peracylated 2Z-benzoyloxyimino-glycopyranosyl bromides invariably proceed with retention of the Z-geometry. Accordingly, the many β-d-hexosidulose oximes in literature which were prepared in this way and for which the oxime geometry has not been addressed explicitly, are the Z-oximes throughout. By contrast, oximation of β-d-hexopyranosid-2-uloses leads to mixtures of E and Z oximes readily separable and structurally verifiable by ¹H and ¹³C NMR. Configurational assignments rested on comparative evaluation of NMR data of E and Z isomers, and, most notably on an X-ray structural analysis of the pivaloylated isopropyl 2E-benzoyloxyimino-2-deoxy-β-d-arabino-hexopyranoside revealing the unusual ¹S₅?^1,4B conformation for the pyranoid ring.     

Metabolic and energetic aspects of biohydrogen production of Clostridium tyrobutyricum: The effects of hydraulic retention time and peptone addition

This study evaluates the microbial metabolism and energy demand in fermentative biohydrogen production using Clostridium tyrobutyricum FYa102 at different hydraulic retention times (HRT) over a period of 1-18 h. The hydrogen yield shows a positive correlation with the butyrate yield, the B/A ratio, and the Y_H2/2(Y_HAc+Y_HBu) ratio, but a negative correlation with the lactate yield. A decrease in HRT, which is accompanied by an increased biomass growth, tends to decrease the B/A ratio, due presumably to a higher energy demand for microbial growth. The production of lactate at a low HRT, however, may involve an unfavorable change in e⁻ equiv distribution to result in a reduced hydrogen production. Finally, the relatively high hydrogen yields observed in the bioreactor with the peptone addition may be ascribed to the utilization of peptone as an additional energy and/or amino-acid source, thus reducing the glucose demand for biomass growth during the hydrogen production process.     

The xanthophyll cycle and antioxidative defense system are enhanced in the wheat hybrid subjected to high light stress

Although the wheat hybrids have often shown higher grain yields, the physiological basis of the higher yields remains unknown. Previous studies suggest that tolerance to photoinhibition in the hybrid may be one of the physiological bases (Yang et al., 2006, Plant Sci 171:389-97). The objective of this study was to further investigate the possible mechanism responsible for tolerance to photoinhibition in the hybrid. Photosystem II (PSII) photochemistry, the xanthophyll cycle, and antioxidative defense system were compared between the hybrid and its parents subjected to high light stress (1500 μmol m⁻² s⁻¹). The analyses of oxygen-evolving activity, chlorophyll fluorescence, and protein blotting demonstrated that the higher tolerance in the hybrid than in its parents was associated with its higher tolerance of PSII to photoinhibition. High light induced an increase in non-photochemical quenching, and this increase was greater in the hybrid than in its parents. There were no differences in the pool size of the xanthophyll cycle between the hybrid and its parents. The content of violaxanthin decreased significantly, whereas the content of zeaxanthin + antherxanthin increased considerably during high light treatments. However, the decrease in violaxanthin content and the increase in zeaxanthin + antherxanthin content were greater in the hybrid than in its parents. High light resulted in a significant accumulation of H₂O₂, O₂⁻ and catalytic Fe, and this accumulation was less in the hybrid than in its parents. High light induced a significant increase in the activities of superoxide dismutase, catalase, ascorbate peroxidase, glutathione reductase, dehydroascorbate reductase, and monodehydroascorbate reductase, and these increases were greater in the hybrid than its parents. These results suggest that the higher tolerance to photoinhibition in the hybrid may be associated with its higher capacity for antioxidative defense metabolism and the xanthophyll cycle.     

Cytochrome c oxidase loses catalytic activity and structural integrity during the aging process in Drosophila melanogaster

The hypothesis, that structural deterioration of cytochrome c oxidase (CcO) is a causal factor in the age-related decline in mitochondrial respiratory activity and an increase in H₂O₂ generation, was tested in Drosophila melanogaster. CcO activity and the levels of seven different nuclear DNA-encoded CcO subunits were determined at three different stages of adult life, namely, young-, middle-, and old-age. CcO activity declined progressively with age by 33%. Western blot analysis, using antibodies specific to Drosophila CcO subunits IV, Va, Vb, VIb, VIc, VIIc, and VIII, indicated that the abundance these polypeptides decreased, ranging from 11% to 40%, during aging. These and previous results suggest that CcO is a specific intra-mitochondrial site of age-related deterioration, which may have a broad impact on mitochondrial physiology.     

Molecular characterization of β1,4-galactosyltransferase 7 genetic mutations linked to the progeroid form of Ehlers-Danlos syndrome (EDS)

β1,4-Galactosyltransferase 7 (β4GalT7) is a key enzyme initiating glycosaminoglycan (GAG) synthesis. Based on in vitro and ex vivo kinetics studies and structure-based modelling, we molecularly characterized β4GalT7 mutants linked to the progeroid form of Ehlers-Danlos syndrome (EDS), a severe connective tissue disorder. Our results revealed that loss of activity upon L206P substitution due to altered protein folding is the primary cause for the GAG synthesis defect in patients carrying the compound A186D and L206P mutations. We showed that R270C substitution strongly reduced β4GalT7 affinity towards xyloside acceptor, thus affecting GAG chains formation. This study establishes the molecular basis for β4GalT7 defects associated with altered GAG synthesis in EDS.