Stochastic slowdown in evolutionary processes

We examine birth--death processes with state dependent transition probabilities and at least one absorbing boundary. In evolution, this describes selection acting on two different types in a finite population where reproductive events occur successively. If the two types have equal fitness the system performs a random walk. If one type has a fitness advantage it is favored by selection, which introduces a bias (asymmetry) in the transition probabilities. How long does it take until advantageous mutants have invaded and taken over? Surprisingly, we find that the average time of such a process can increase, even if the mutant type always has a fitness advantage. We discuss this finding for the Moran process and develop a simplified model which allows a more intuitive understanding. We show that this effect can occur for weak but non--vanishing bias (selection) in the state dependent transition rates and infer the scaling with system size. We also address the Wright-Fisher model commonly used in population genetics, which shows that this stochastic slowdown is not restricted to birth-death processes.Comment: 8 pages, 3 figures, accepted for publicatio

Cryptic MHC Polymorphism Revealed but Not Explained by Selection on the Class IIB Peptide-Binding Region

The immune genes of the major histocompatibility complex (MHC) are characterized by extraordinarily high levels of nucleotide and haplotype diversity. This variation is maintained by pathogen-mediated balancing selection that is operating on the peptide-binding region (PBR). Several recent studies have found, however, that some populations possess large clusters of alleles that are translated into virtually identical proteins. Here, we address the question of how this nucleotide polymorphism is maintained with little or no functional variation for selection to operate on. We investigate circa 750–850 bp of MHC class II DAB genes in four wild populations of the guppy Poecilia reticulata. By sequencing an extended region, we uncovered 40.9% more sequences (alleles), which would have been missed if we had amplified the exon 2 alone. We found evidence of several gene conversion events that may have homogenized sequence variation. This reduces the visible copy number variation (CNV) and can result in a systematic underestimation of the CNV in studies of the MHC and perhaps other multigene families. We then focus on a single cluster, which comprises 27 (of a total of 66) sequences. These sequences are virtually identical and show no signal of selection. We use microsatellites to reconstruct the populations' demography and employ simulations to examine whether so many similar nucleotide sequences can be maintained in the populations. Simulations show that this variation does not behave neutrally. We propose that selection operates outside the PBR, for example, on linked immune genes or on the “sheltered load” that is thought to be associated to the MHC. Future studies on the MHC would benefit from extending the amplicon size to include polymorphisms outside the exon with the PBR. This may capture otherwise cryptic haplotype variation and CNV, and it may help detect other regions in the MHC that are under selection

Phylogenetic relationships of Indian caecilians (Amphibia: Gymnophiona) inferred from mitochondrial rRNA gene sequences

India has a diverse caecilian fauna, including representatives of three of the six currently recognized families, the Caeciliidae, Ichthyophiidae, the endemic Uraeotyphlidae, but previous molecular phylogenetic studies of caecilians have not included sequences for any Indian caecilians. Partial 12S and 16S mitochondrial gene sequences were obtained for a single representative of each of the caecilian families found in India and aligned against previously reported sequences for 13 caecilian species. The resulting alignment (16 taxa, 1200 sites, of which 288 cannot be aligned unambiguously) was analyzed using parsimony, maximum-likelihood, and distance methods. As judged by bootstrap proportions, decay indices, and leaf stabilities, well-supported relationships of the Indian caecilians are recovered from the alignment. The data (1) corroborate the hypothesis, based on morphology, that the Uraeotyphlidae and Ichthyophiidae are sister taxa, (2) recover a monophyletic Ichthyophiidae, including Indian and South East Asian representatives, and (3) place the Indian caeciliid Gegeneophis ramaswamii as the sister group of the caeciliid caecilians of the Seychelles. Rough estimates of divergence times suggest an origin of the Uraeotyphlidae and Ichthyophiidae while India was isolated from Laurasia and Africa and are most consistent with an Indian origin of these families and subsequent dispersal of ichthyophiids into South East Asia

Groundwater quality comparison between rural farms and riparian wells in the western Amazon, Brazil

Groundwater quality of a riparian forest is compared to wells in surrounding rural areas at Urupá River basin. Groundwater types were calcium bicarbonated at left margin and sodium chloride at right, whereas riparian wells exhibited a combination of both (sodium bicarbonate). Groundwater was mostly solute-depleted with concentrations within permissible limits for human consumption, except for nitrate. Isotopic composition suggests that inorganic carbon in Urupá River is mostly supplied by runoff instead of riparian groundwater. Hence, large pasture areas in addition to narrow riparian forest width in this watershed may have an important contribution in the chemical composition of this river.FAPESPCNPq - CT-HIDROCNPq - Milêni

X-ray Structure And Molecular Dynamics Simulations Of Endoglucanase 3 From Trichoderma Harzianum: Structural Organization And Substrate Recognition By Endoglucanases That Lack Cellulose Binding Module.

Plant biomass holds a promise for the production of second-generation ethanol via enzymatic hydrolysis, but its utilization as a biofuel resource is currently limited to a large extent by the cost and low efficiency of the cellulolytic enzymes. Considerable efforts have been dedicated to elucidate the mechanisms of the enzymatic process. It is well known that most cellulases possess a catalytic core domain and a carbohydrate binding module (CBM), without which the enzymatic activity can be drastically reduced. However, Cel12A members of the glycosyl hydrolases family 12 (GHF12) do not bear a CBM and yet are able to hydrolyze amorphous cellulose quite efficiently. Here, we use X-ray crystallography and molecular dynamics simulations to unravel the molecular basis underlying the catalytic capability of endoglucanase 3 from Trichoderma harzianum (ThEG3), a member of the GHF12 enzymes that lacks a CBM. A comparative analysis with the Cellulomonas fimi CBM identifies important residues mediating interactions of EG3s with amorphous regions of the cellulose. For instance, three aromatic residues constitute a harboring wall of hydrophobic contacts with the substrate in both ThEG3 and CfCBM structures. Moreover, residues at the entrance of the active site cleft of ThEG3 are identified, which might hydrogen bond to the substrate. We advocate that the ThEG3 residues Asn152 and Glu201 interact with the substrate similarly to the corresponding CfCBM residues Asn81 and Arg75. Altogether, these results show that CBM motifs are incorporated within the ThEG3 catalytic domain and suggest that the enzymatic efficiency is associated with the length and position of the substrate chain, being higher when the substrate interact with the aromatic residues at the entrance of the cleft and the catalytic triad. Our results provide guidelines for rational protein engineering aiming to improve interactions of GHF12 enzymes with cellulosic substrates.8e5906

Poisson structures for reduced non-holonomic systems

Borisov, Mamaev and Kilin have recently found certain Poisson structures with respect to which the reduced and rescaled systems of certain non-holonomic problems, involving rolling bodies without slipping, become Hamiltonian, the Hamiltonian function being the reduced energy. We study further the algebraic origin of these Poisson structures, showing that they are of rank two and therefore the mentioned rescaling is not necessary. We show that they are determined, up to a non-vanishing factor function, by the existence of a system of first-order differential equations providing two integrals of motion. We generalize the form of that Poisson structures and extend their domain of definition. We apply the theory to the rolling disk, the Routh's sphere, the ball rolling on a surface of revolution, and its special case of a ball rolling inside a cylinder.Comment: 22 page

Bayesian paternity analysis and mating patterns in a parasitic nematode, Trichostrongylus tenuis

Mating behaviour is a fundamental aspect of the evolutionary ecology of sexually reproducing species, but one that has been under-researched in parasitic nematodes. We analysed mating behaviour in the parasitic nematode Trichostrongylus tenuis by performing a paternity analysis in a population from a single red grouse host. Paternity of the 150 larval offspring of 25 mothers (sampled from one of the two host caeca) was assigned among 294 candidate fathers (sampled from both caeca). Each candidate father's probability of paternity of each offspring was estimated from 10-locus microsatellite genotypes. Seventy-six (51%) offspring were assigned a father with a probability of >0.8, and the estimated number of unsampled males was 136 (95% credible interval (CI) 77-219). The probability of a male from one caecum fathering an offspring in the other caecum was estimated as 0.024 (95% CI 0.003-0.077), indicating that the junction of the caeca is a strong barrier to dispersal. Levels of promiscuity (defined as the probability of two of an adult's offspring sharing only one parent) were high for both sexes. Variance in male reproductive success was moderately high, possibly because of a combination of random mating and high variance in post-copulatory reproductive success. These results provide the first data on individual mating behaviour among parasitic nematodes