Fast and Robust Characterization of Time-Heterogeneous Sequence Evolutionary Processes Using Substitution Mapping

Genes and genomes do not evolve similarly in all branches of the tree of life. Detecting and characterizing the heterogeneity in time, and between lineages, of the nucleotide (or amino acid) substitution process is an important goal of current molecular evolutionary research. This task is typically achieved through the use of non-homogeneous models of sequence evolution, which being highly parametrized and computationally-demanding are not appropriate for large-scale analyses. Here we investigate an alternative methodological option based on probabilistic substitution mapping. The idea is to first reconstruct the substitutional history of each site of an alignment under a homogeneous model of sequence evolution, then to characterize variations in the substitution process across lineages based on substitution counts. Using simulated and published datasets, we demonstrate that probabilistic substitution mapping is robust in that it typically provides accurate reconstruction of sequence ancestry even when the true process is heterogeneous, but a homogeneous model is adopted. Consequently, we show that the new approach is essentially as efficient as and extremely faster than (up to 25 000 times) existing methods, thus paving the way for a systematic survey of substitution process heterogeneity across genes and lineages

BASE-STEP: A transportable antiproton reservoir for fundamental interaction studies

Currently, the only worldwide source of low-energy antiprotons is the AD/ELENA facility located at CERN. To date, all precision measurements on single antiprotons have been conducted at this facility and provide stringent tests of the fundamental interactions and their symmetries. However, the magnetic field fluctuations from the facility operation limit the precision of upcoming measurements. To overcome this limitation, we have designed the transportable antiproton trap system BASE-STEP to relocate antiprotons to laboratories with a calm magnetic environment. We anticipate that the transportable antiproton trap will facilitate enhanced tests of CPT invariance with antiprotons, and provide new experimental possibilities of using transported antiprotons and other accelerator-produced exotic ions. We present here the technical design of the transportable trap system. This includes the transportable superconducting magnet, the cryogenic inlay consisting of the trap stack and the detection systems, and the differential pumping section to suppress the residual gas flow into the cryogenic trap chamber.Comment: To be submitted to Rev. Sci. Instrument

Lakeside Cemeteries in the Sahara: 5000 Years of Holocene Population and Environmental Change

Background: Approximately two hundred human burials were discovered on the edge of a paleolake in Niger that providea uniquely preserved record of human occupation in the Sahara during the Holocene (,8000 B.C.E. to the present). CalledGobero, this suite of closely spaced sites chronicles the rapid pace of biosocial change in the southern Sahara in response tosevere climatic fluctuation.Methodology/Principal Findings: Two main occupational phases are identified that correspond with humid intervals in theearly and mid-Holocene, based on 78 direct AMS radiocarbon dates on human remains, fauna and artifacts, as well as 9 OSLdates on paleodune sand. The older occupants have craniofacial dimensions that demonstrate similarities with mid-Holocene occupants of the southern Sahara and Late Pleistocene to early Holocene inhabitants of the Maghreb. Theirhyperflexed burials compose the earliest cemetery in the Sahara dating to ,7500 B.C.E. These early occupants abandon thearea under arid conditions and, when humid conditions return ,4600 B.C.E., are replaced by a more gracile people withelaborated grave goods including animal bone and ivory ornaments.Conclusions/Significance: The principal significance of Gobero lies in its extraordinary human, faunal, and archaeologicalrecord, from which we conclude the following:(1) The early Holocene occupants at Gobero (7700–6200 B.C.E.) were largely sedentary hunter-fisher-gatherers withlakeside funerary sites that include the earliest recorded cemetery in the Sahara.(2) Principal components analysis of craniometric variables closely allies the early Holocene occupants at Gobero with askeletally robust, trans-Saharan assemblage of Late Pleistocene to mid-Holocene human populations from the Maghreband southern Sahara.(3) Gobero was abandoned during a period of severe aridification possibly as long as one millennium (6200–5200 B.C.E).(4) More gracile humans arrived in the mid-Holocene (5200–2500 B.C.E.) employing a diversified subsistence economybased on clams, fish, and savanna vertebrates as well as some cattle husbandry.(5) Population replacement after a harsh arid hiatus is the most likely explanation for the occupational sequence at Gobero.(6) We are just beginnin

Bio++: Efficient Extensible Libraries and Tools for Computational Molecular Evolution

Efficient algorithms and programs for the analysis of the ever-growing amount of biological sequence data are strongly needed in the genomics era. The pace at which new data and methodologies are generated calls for the use of pre-existing, optimized—yet extensible—code, typically distributed as libraries or packages. This motivated the Bio++ project, aiming at developing a set of C++ libraries for sequence analysis, phylogenetics, population genetics, and molecular evolution. The main attractiveness of Bio++ is the extensibility and reusability of its components through its object-oriented design, without compromising the computer-efficiency of the underlying methods. We present here the second major release of the libraries, which provides an extended set of classes and methods. These extensions notably provide built-in access to sequence databases and new data structures for handling and manipulating sequences from the omics era, such as multiple genome alignments and sequencing reads libraries. More complex models of sequence evolution, such as mixture models and generic n-tuples alphabets, are also included

Evolutionary Modeling of Rate Shifts Reveals Specificity Determinants in HIV-1 Subtypes

A hallmark of the human immunodeficiency virus 1 (HIV-1) is its rapid rate of evolution within and among its various subtypes. Two complementary hypotheses are suggested to explain the sequence variability among HIV-1 subtypes. The first suggests that the functional constraints at each site remain the same across all subtypes, and the differences among subtypes are a direct reflection of random substitutions, which have occurred during the time elapsed since their divergence. The alternative hypothesis suggests that the functional constraints themselves have evolved, and thus sequence differences among subtypes in some sites reflect shifts in function. To determine the contribution of each of these two alternatives to HIV-1 subtype evolution, we have developed a novel Bayesian method for testing and detecting site-specific rate shifts. The RAte Shift EstimatoR (RASER) method determines whether or not site-specific functional shifts characterize the evolution of a protein and, if so, points to the specific sites and lineages in which these shifts have most likely occurred. Applying RASER to a dataset composed of large samples of HIV-1 sequences from different group M subtypes, we reveal rampant evolutionary shifts throughout the HIV-1 proteome. Most of these rate shifts have occurred during the divergence of the major subtypes, establishing that subtype divergence occurred together with functional diversification. We report further evidence for the emergence of a new sub-subtype, characterized by abundant rate-shifting sites. When focusing on the rate-shifting sites detected, we find that many are associated with known function relating to viral life cycle and drug resistance. Finally, we discuss mechanisms of covariation of rate-shifting sites

Prevalence of Epistasis in the Evolution of Influenza A Surface Proteins

The surface proteins of human influenza A viruses experience positive selection to escape both human immunity and, more recently, antiviral drug treatments. In bacteria and viruses, immune-escape and drug-resistant phenotypes often appear through a combination of several mutations that have epistatic effects on pathogen fitness. However, the extent and structure of epistasis in influenza viral proteins have not been systematically investigated. Here, we develop a novel statistical method to detect positive epistasis between pairs of sites in a protein, based on the observed temporal patterns of sequence evolution. The method rests on the simple idea that a substitution at one site should rapidly follow a substitution at another site if the sites are positively epistatic. We apply this method to the surface proteins hemagglutinin and neuraminidase of influenza A virus subtypes H3N2 and H1N1. Compared to a non-epistatic null distribution, we detect substantial amounts of epistasis and determine the identities of putatively epistatic pairs of sites. In particular, using sequence data alone, our method identifies epistatic interactions between specific sites in neuraminidase that have recently been demonstrated, in vitro, to confer resistance to the drug oseltamivir; these epistatic interactions are responsible for widespread drug resistance among H1N1 viruses circulating today. This experimental validation demonstrates the predictive power of our method to identify epistatic sites of importance for viral adaptation and public health. We conclude that epistasis plays a large role in shaping the molecular evolution of influenza viruses. In particular, sites with , which would normally not be identified as positively selected, can facilitate viral adaptation through epistatic interactions with their partner sites. The knowledge of specific interactions among sites in influenza proteins may help us to predict the course of antigenic evolution and, consequently, to select more appropriate vaccines and drugs

A Phylogenetic Perspective on the Evolution of Mediterranean Teleost Fishes

The Mediterranean Sea is a highly diverse, highly studied, and highly impacted biogeographic region, yet no phylogenetic reconstruction of fish diversity in this area has been published to date. Here, we infer the timing and geographic origins of Mediterranean teleost species diversity using nucleotide sequences collected from GenBank. We assembled a DNA supermatrix composed of four mitochondrial genes (12S ribosomal DNA, 16S ribosomal DNA, cytochrome c oxidase subunit I and cytochrome b) and two nuclear genes (rhodopsin and recombination activating gene I), including 62% of Mediterranean teleost species plus 9 outgroups. Maximum likelihood and Bayesian phylogenetic and dating analyses were calibrated using 20 fossil constraints. An additional 124 species were grafted onto the chronogram according to their taxonomic affinity, checking for the effects of taxonomic coverage in subsequent diversification analyses. We then interpreted the time-line of teleost diversification in light of Mediterranean historical biogeography, distinguishing non-endemic natives, endemics and exotic species. Results show that the major Mediterranean orders are of Cretaceous origin, specifically ∼100–80 Mya, and most Perciformes families originated 80–50 Mya. Two important clade origin events were detected. The first at 100–80 Mya, affected native and exotic species, and reflects a global diversification period at a time when the Mediterranean Sea did not yet exist. The second occurred during the last 50 Mya, and is noticeable among endemic and native species, but not among exotic species. This period corresponds to isolation of the Mediterranean from Indo-Pacific waters before the Messinian salinity crisis. The Mediterranean fish fauna illustrates well the assembly of regional faunas through origination and immigration, where dispersal and isolation have shaped the emergence of a biodiversity hotspot

A primary electron beam facility at CERN -- eSPS Conceptual design report

The design of a primary electron beam facility at CERN is described. The study has been carried out within the framework of the wider Physics Beyond Colliders study. It re-enables the Super Proton Synchrotron (SPS) as an electron accelerator, and leverages the development invested in Compact Linear Collider (CLIC) technology for its injector and as an accelerator research and development infrastructure. The facility would be relevant for several of the key priorities in the 2020 update of the European Strategy for Particle Physics, such as an electron-positron Higgs factory, accelerator R\&D, dark sector physics, and neutrino physics. In addition, it could serve experiments in nuclear physics. The electron beam delivered by this facility would provide access to light dark matter production significantly beyond the targets predicted by a thermal dark matter origin, and for natures of dark matter particles that are not accessible by direct detection experiments. It would also enable electro-nuclear measurements crucial for precise modelling the energy dependence of neutrino-nucleus interactions, which is needed to precisely measure neutrino oscillations as a function of energy. The implementation of the facility is the natural next step in the development of X-band high-gradient acceleration technology, a key technology for compact and cost-effective electron/positron linacs. It would also become the only facility with multi-GeV drive bunches and truly independent electron witness bunches for plasma wakefield acceleration. A second phase capable to deliver positron witness bunches would make it a complete facility for plasma wakefield collider studies. [...