The PDZ domain of the SpoIVB serine peptidase facilitates multiple functions

During spore formation in Bacillus subtilis, the SpoIVB protein is a critical component of the sigma (K) regulatory checkpoint. SpoIVB has been shown to be a serine peptidase that is synthesized in the spore chamber and which self-cleaves, releasing active forms. These forms can signal proteolytic processing of the transcription factor sigma (K) in the outer mother cell chamber of the sporulating cell. This forms the basis of the sigma (K) checkpoint and ensures accurate sigma (K)-controlled gene expression. SpoIVB has also been shown to activate a second distinct process, termed the second function, which is essential for the formation of heat-resistant spores. In addition to the serine peptidase domain, SpoIVB contains a PDZ domain. We have altered a number of conserved residues in the PDZ domain by site-directed mutagenesis and assayed the sporulation phenotype and signaling properties of mutant SpoIVB proteins. Our work has revealed that the SpoIVB PDZ domain could be used for up to four distinct processes, (i) targeting of itself for trans proteolysis, (11) binding to the protease inhibitor BofC, (iii) signaling of pro-sigma (K) processing, and (iv) signaling of the second function of SpoIVB

Non-L\'evy mobility patterns of Mexican Me'Phaa peasants searching for fuelwood

We measured mobility patterns that describe walking trajectories of individual Me'Phaa peasants searching and collecting fuelwood in the forests of "La Monta\~na de Guerrero" in Mexico. These one-day excursions typically follow a mixed pattern of nearly-constant steps when individuals displace from their homes towards potential collecting sites and a mixed pattern of steps of different lengths when actually searching for fallen wood in the forest. Displacements in the searching phase seem not to be compatible with L\'evy flights described by power-laws with optimal scaling exponents. These findings however can be interpreted in the light of deterministic searching on heavily degraded landscapes where the interaction of the individuals with their scarce environment produces alternative searching strategies than the expected L\'evy flights. These results have important implications for future management and restoration of degraded forests and the improvement of the ecological services they may provide to their inhabitants.Comment: 15 pages, 4 figures. First version submitted to Human Ecology. The final publication will be available at http://www.springerlink.co

Evaluation of a candidate breast cancer associated SNP in ERCC4 as a risk modifier in BRCA1 and BRCA2 mutation carriers. Results from the Consortium of Investigators of Modifiers of BRCA1/BRCA2 (CIMBA)

Background: In this study we aimed to evaluate the role of a SNP in intron 1 of the ERCC4 gene (rs744154), previously reported to be associated with a reduced risk of breast cancer in the general population, as a breast cancer risk modifier in BRCA1 and BRCA2 mutation carriers. Methods: We have genotyped rs744154 in 9408 BRCA1 and 5632 BRCA2 mutation carriers from the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA) and assessed its association with breast cancer risk using a retrospective weighted cohort approach. Results: We found no evidence of association with breast cancer risk for BRCA1 (per-allele HR: 0.98, 95% CI: 0.93–1.04, P=0.5) or BRCA2 (per-allele HR: 0.97, 95% CI: 0.89–1.06, P=0.5) mutation carriers. Conclusion: This SNP is not a significant modifier of breast cancer risk for mutation carriers, though weak associations cannot be ruled out. A Osorio1, R L Milne2, G Pita3, P Peterlongo4,5, T Heikkinen6, J Simard7, G Chenevix-Trench8, A B Spurdle8, J Beesley8, X Chen8, S Healey8, KConFab9, S L Neuhausen10, Y C Ding10, F J Couch11,12, X Wang11, N Lindor13, S Manoukian4, M Barile14, A Viel15, L Tizzoni5,16, C I Szabo17, L Foretova18, M Zikan19, K Claes20, M H Greene21, P Mai21, G Rennert22, F Lejbkowicz22, O Barnett-Griness22, I L Andrulis23,24, H Ozcelik24, N Weerasooriya23, OCGN23, A-M Gerdes25, M Thomassen25, D G Cruger26, M A Caligo27, E Friedman28,29, B Kaufman28,29, Y Laitman28, S Cohen28, T Kontorovich28, R Gershoni-Baruch30, E Dagan31,32, H Jernström33, M S Askmalm34, B Arver35, B Malmer36, SWE-BRCA37, S M Domchek38, K L Nathanson38, J Brunet39, T Ramón y Cajal40, D Yannoukakos41, U Hamann42, HEBON37, F B L Hogervorst43, S Verhoef43, EB Gómez García44,45, J T Wijnen46,47, A van den Ouweland48, EMBRACE37, D F Easton49, S Peock49, M Cook49, C T Oliver49, D Frost49, C Luccarini50, D G Evans51, F Lalloo51, R Eeles52, G Pichert53, J Cook54, S Hodgson55, P J Morrison56, F Douglas57, A K Godwin58, GEMO59,60,61, O M Sinilnikova59,60, L Barjhoux59,60, D Stoppa-Lyonnet61, V Moncoutier61, S Giraud59, C Cassini62,63, L Olivier-Faivre62,63, F Révillion64, J-P Peyrat64, D Muller65, J-P Fricker65, H T Lynch66, E M John67, S Buys68, M Daly69, J L Hopper70, M B Terry71, A Miron72, Y Yassin72, D Goldgar73, Breast Cancer Family Registry37, C F Singer74, D Gschwantler-Kaulich74, G Pfeiler74, A-C Spiess74, Thomas v O Hansen75, O T Johannsson76, T Kirchhoff77, K Offit77, K Kosarin77, M Piedmonte78, G C Rodriguez79, K Wakeley80, J F Boggess81, J Basil82, P E Schwartz83, S V Blank84, A E Toland85, M Montagna86, C Casella87, E N Imyanitov88, A Allavena89, R K Schmutzler90, B Versmold90, C Engel91, A Meindl92, N Ditsch93, N Arnold94, D Niederacher95, H Deißler96, B Fiebig97, R Varon-Mateeva98, D Schaefer99, U G Froster100, T Caldes101, M de la Hoya101, L McGuffog49, A C Antoniou49, H Nevanlinna6, P Radice4,5 and J Benítez1,3 on behalf of CIMB

Search for the Proton Decay Mode proton to neutrino K+ in Soudan 2

We have searched for the proton decay mode proton to neutrino K+ using the one-kiloton Soudan 2 high resolution calorimeter. Contained events obtained from a 3.56 kiloton-year fiducial exposure through June 1997 are examined for occurrence of a visible K+ track which decays at rest into mu+ nu or pi+ pi0. We found one candidate event consistent with background, yielding a limit, tau/B > 4.3 10^{31} years at 90% CL with no background subtraction.Comment: 13 pages, Latex, 3 tables and 3 figures, Accepted by Physics Letters

GENOMEPOP: A program to simulate genomes in populations

<p>Abstract</p> <p>Background</p> <p>There are several situations in population biology research where simulating DNA sequences is useful. Simulation of biological populations under different evolutionary genetic models can be undertaken using backward or forward strategies. Backward simulations, also called coalescent-based simulations, are computationally efficient. The reason is that they are based on the history of lineages with surviving offspring in the current population. On the contrary, forward simulations are less efficient because the entire population is simulated from past to present. However, the coalescent framework imposes some limitations that forward simulation does not. Hence, there is an increasing interest in forward population genetic simulation and efficient new tools have been developed recently. Software tools that allow efficient simulation of large DNA fragments under complex evolutionary models will be very helpful when trying to better understand the trace left on the DNA by the different interacting evolutionary forces. Here I will introduce GenomePop, a forward simulation program that fulfills the above requirements. The use of the program is demonstrated by studying the impact of intracodon recombination on global and site-specific <it>dN/dS </it>estimation.</p> <p>Results</p> <p>I have developed algorithms and written software to efficiently simulate, forward in time, different Markovian nucleotide or codon models of DNA mutation. Such models can be combined with recombination, at inter and intra codon levels, fitness-based selection and complex demographic scenarios.</p> <p>Conclusion</p> <p>GenomePop has many interesting characteristics for simulating SNPs or DNA sequences under complex evolutionary and demographic models. These features make it unique with respect to other simulation tools. Namely, the possibility of forward simulation under General Time Reversible (GTR) mutation or GTR×MG94 codon models with intra-codon recombination, arbitrary, user-defined, migration patterns, diploid or haploid models, constant or variable population sizes, etc. It also allows simulation of fitness-based selection under different distributions of mutational effects. Under the 2-allele model it allows the simulation of recombination hot-spots, the definition of different frequencies in different populations, etc. GenomePop can also manage large DNA fragments. In addition, it has a scaling option to save computation time when simulating large sequences and population sizes under complex demographic and evolutionary situations. These and many other features are detailed in its web page <abbrgrp><abbr bid="B1">1</abbr></abbrgrp>.</p

The Probability of a Gene Tree Topology within a Phylogenetic Network with Applications to Hybridization Detection

Gene tree topologies have proven a powerful data source for various tasks, including species tree inference and species delimitation. Consequently, methods for computing probabilities of gene trees within species trees have been developed and widely used in probabilistic inference frameworks. All these methods assume an underlying multispecies coalescent model. However, when reticulate evolutionary events such as hybridization occur, these methods are inadequate, as they do not account for such events. Methods that account for both hybridization and deep coalescence in computing the probability of a gene tree topology currently exist for very limited cases. However, no such methods exist for general cases, owing primarily to the fact that it is currently unknown how to compute the probability of a gene tree topology within the branches of a phylogenetic network. Here we present a novel method for computing the probability of gene tree topologies on phylogenetic networks and demonstrate its application to the inference of hybridization in the presence of incomplete lineage sorting. We reanalyze a Saccharomyces species data set for which multiple analyses had converged on a species tree candidate. Using our method, though, we show that an evolutionary hypothesis involving hybridization in this group has better support than one of strict divergence. A similar reanalysis on a group of three Drosophila species shows that the data is consistent with hybridization. Further, using extensive simulation studies, we demonstrate the power of gene tree topologies at obtaining accurate estimates of branch lengths and hybridization probabilities of a given phylogenetic network. Finally, we discuss identifiability issues with detecting hybridization, particularly in cases that involve extinction or incomplete sampling of taxa

Evolution of Disease Response Genes in Loblolly Pine: Insights from Candidate Genes

BACKGROUND: Host-pathogen interactions that may lead to a competitive co-evolution of virulence and resistance mechanisms present an attractive system to study molecular evolution because strong, recent (or even current) selective pressure is expected at many genomic loci. However, it is unclear whether these selective forces would act to preserve existing diversity, promote novel diversity, or reduce linked neutral diversity during rapid fixation of advantageous alleles. In plants, the lack of adaptive immunity places a larger burden on genetic diversity to ensure survival of plant populations. This burden is even greater if the generation time of the plant is much longer than the generation time of the pathogen. METHODOLOGY/PRINCIPAL FINDINGS: Here, we present nucleotide polymorphism and substitution data for 41 candidate genes from the long-lived forest tree loblolly pine, selected primarily for their prospective influences on host-pathogen interactions. This dataset is analyzed together with 15 drought-tolerance and 13 wood-quality genes from previous studies. A wide range of neutrality tests were performed and tested against expectations from realistic demographic models. CONCLUSIONS/SIGNIFICANCE: Collectively, our analyses found that axr (auxin response factor), caf1 (chromatin assembly factor) and gatabp1 (gata binding protein 1) candidate genes carry patterns consistent with directional selection and erd3 (early response to drought 3) displays patterns suggestive of a selective sweep, both of which are consistent with the arm-race model of disease response evolution. Furthermore, we have identified patterns consistent with diversifying selection at erf1-like (ethylene responsive factor 1), ccoaoemt (caffeoyl-CoA-O-methyltransferase), cyp450-like (cytochrome p450-like) and pr4.3 (pathogen response 4.3), expected under the trench-warfare evolution model. Finally, a drought-tolerance candidate related to the plant cell wall, lp5, displayed patterns consistent with balancing selection. In conclusion, both arms-race and trench-warfare models seem compatible with patterns of polymorphism found in different disease-response candidate genes, indicating a mixed strategy of disease tolerance evolution for loblolly pine, a major tree crop in southeastern United States

The Genomics of Speciation in Drosophila: Diversity, Divergence, and Introgression Estimated Using Low-Coverage Genome Sequencing

In nature, closely related species may hybridize while still retaining their distinctive identities. Chromosomal regions that experience reduced recombination in hybrids, such as within inversions, have been hypothesized to contribute to the maintenance of species integrity. Here, we examine genomic sequences from closely related fruit fly taxa of the Drosophila pseudoobscura subgroup to reconstruct their evolutionary histories and past patterns of genic exchange. Partial genomic assemblies were generated from two subspecies of Drosophila pseudoobscura (D. ps.) and an outgroup species, D. miranda. These new assemblies were compared to available assemblies of D. ps. pseudoobscura and D. persimilis, two species with overlapping ranges in western North America. Within inverted regions, nucleotide divergence among each pair of the three species is comparable, whereas divergence between D. ps. pseudoobscura and D. persimilis in non-inverted regions is much lower and closer to levels of intraspecific variation. Using molecular markers flanking each of the major chromosomal inversions, we identify strong crossover suppression in F1 hybrids extending over 2 megabase pairs (Mbp) beyond the inversion breakpoints. These regions of crossover suppression also exhibit the high nucleotide divergence associated with inverted regions. Finally, by comparison to a geographically isolated subspecies, D. ps. bogotana, our results suggest that autosomal gene exchange between the North American species, D. ps. pseudoobscura and D. persimilis, occurred since the split of the subspecies, likely within the last 200,000 years. We conclude that chromosomal rearrangements have been vital to the ongoing persistence of these species despite recent hybridization. Our study serves as a proof-of-principle on how whole genome sequencing can be applied to formulate and test hypotheses about species formation in lesser-known non-model systems

Levels and Patterns of Nucleotide Variation in Domestication QTL Regions on Rice Chromosome 3 Suggest Lineage-Specific Selection

Oryza sativa or Asian cultivated rice is one of the major cereal grass species domesticated for human food use during the Neolithic. Domestication of this species from the wild grass Oryza rufipogon was accompanied by changes in several traits, including seed shattering, percent seed set, tillering, grain weight, and flowering time. Quantitative trait locus (QTL) mapping has identified three genomic regions in chromosome 3 that appear to be associated with these traits. We would like to study whether these regions show signatures of selection and whether the same genetic basis underlies the domestication of different rice varieties. Fragments of 88 genes spanning these three genomic regions were sequenced from multiple accessions of two major varietal groups in O. sativa—indica and tropical japonica—as well as the ancestral wild rice species O. rufipogon. In tropical japonica, the levels of nucleotide variation in these three QTL regions are significantly lower compared to genome-wide levels, and coalescent simulations based on a complex demographic model of rice domestication indicate that these patterns are consistent with selection. In contrast, there is no significant reduction in nucleotide diversity in the homologous regions in indica rice. These results suggest that there are differences in the genetic and selective basis for domestication between these two Asian rice varietal groups