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

Leptin levels were negatively associated with lumbar spine bone mineral content in children with overweight or obesity

Aim: Adipokines seem to play a role in bone morphogenesis, although this also depends on the mechanical forces applied to the skeleton. The aim was to assess the relationships of resting leptin and adiponectin with bone parameters and whether high muscular fitness levels affect these relationships in children with overweight or obesity. Methods: This cross-sectional study took part from 2014 to 2016 in Granada, Spain. Participants were recruited from University Hospitals, and we also used advertisements in local media and school contacts in the city. Adipokines were analysed in plasma. Muscular fitness was assessed by one repetition maximum in bench and leg press tests. Dual-energy X-ray absorptiometry was used to measure bone parameters. Results: We included 84 children (10.0 ± 1.2y; 63% boys) in this analysis. Leptin was negatively associated with lumbar spine bone mineral content (β = −0.162, p = 0.053)

Author Correction: A ferroptosis–based panel of prognostic biomarkers for Amyotrophic Lateral Sclerosis

Correction to: Scientific Reports https://doi.org/10.1038/s41598-019-39739-5, published online 27 February 201

Inference of Co-Evolving Site Pairs: an Excellent Predictor of Contact Residue Pairs in Protein 3D structures

Residue-residue interactions that fold a protein into a unique three-dimensional structure and make it play a specific function impose structural and functional constraints on each residue site. Selective constraints on residue sites are recorded in amino acid orders in homologous sequences and also in the evolutionary trace of amino acid substitutions. A challenge is to extract direct dependences between residue sites by removing indirect dependences through other residues within a protein or even through other molecules. Recent attempts of disentangling direct from indirect dependences of amino acid types between residue positions in multiple sequence alignments have revealed that the strength of inferred residue pair couplings is an excellent predictor of residue-residue proximity in folded structures. Here, we report an alternative attempt of inferring co-evolving site pairs from concurrent and compensatory substitutions between sites in each branch of a phylogenetic tree. First, branch lengths of a phylogenetic tree inferred by the neighbor-joining method are optimized as well as other parameters by maximizing a likelihood of the tree in a mechanistic codon substitution model. Mean changes of quantities, which are characteristic of concurrent and compensatory substitutions, accompanied by substitutions at each site in each branch of the tree are estimated with the likelihood of each substitution. Partial correlation coefficients of the characteristic changes along branches between sites are calculated and used to rank co-evolving site pairs. Accuracy of contact prediction based on the present co-evolution score is comparable to that achieved by a maximum entropy model of protein sequences for 15 protein families taken from the Pfam release 26.0. Besides, this excellent accuracy indicates that compensatory substitutions are significant in protein evolution.Comment: 17 pages, 4 figures, and 4 tables with supplementary information of 5 figure

Substitution of adeno-associated virus Rep protein binding and nicking sites with human Chromosome 19 sequences

<p>Abstract</p> <p>Background</p> <p>Adeno-associated virus type 2 (AAV2) preferentially integrates its DNA at a ~2 kb region of human chromosome 19, designated <it>AAVS1 </it>(also known as <it>MBS85</it>). Integration at <it>AAVS1 </it>requires the AAV2 replication (Rep) proteins and a DNA sequence within <it>AAVS1 </it>containing a 16 bp Rep recognition sequence (RRS) and closely spaced Rep nicking site (also referred to as a terminal resolution site, or <it>trs</it>). The AAV2 genome is flanked by inverted terminal repeats (ITRs). Each ITR contains an RRS and closely spaced <it>trs</it>, but the sequences differ from those in <it>AAVS1</it>. These ITR sequences are required for replication and packaging.</p> <p>Results</p> <p>In this study we demonstrate that the <it>AAVS1 </it>RRS and <it>trs </it>can function in AAV2 replication, packaging and integration by replacing a 61 bp region of the AAV2 ITR with a 49 bp segment of <it>AAVS1 </it>DNA. Modifying one or both ITRs did not have a large effect on the overall virus titers. These modifications did not detectably affect integration at <it>AAVS1</it>, as measured by semi-quantitative nested PCR assays. Sequencing of integration junctions shows the joining of the modified ITRs to <it>AAVS1 </it>sequences.</p> <p>Conclusions</p> <p>The ability of these <it>AAVS1 </it>sequences to substitute for the AAV2 RRS and <it>trs </it>provides indirect evidence that the stable secondary structure encompassing the <it>trs </it>is part of the AAV2 packaging signal.</p

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

BMJ Open

INTRODUCTION: The prevalence of postnatal depression (PND) is significant: reaching up to 20% in the general population. In mechanistic terms, the risk of PND lies in an interaction between a maternal psychophysiological vulnerability and a chronic environmental context of stress. On the one hand, repetition of stressor during pregnancy mimics a chronic stress model that is relevant to the study of the allostatic load and the adaptive mechanisms. On the other hand, vulnerability factors reflect a psychological profile mirroring mindfulness functioning (psychological quality that involves bringing one's complete and non-judgemental attention to the present experience on a moment-to-moment basis). This psychological resource is linked to protective and resilient psychic functioning. Thus, PND appears to be a relevant model for studying the mechanisms of chronic stress and vulnerability to psychopathologies.In this article, we present the protocol of an ongoing study (started in May 2017). METHODS AND ANALYSIS: The study is being carried out in five maternities and will involve 260 women. We aim to determine the predictive psychobiological factors for PND emergence and to provide a better insight into the mechanisms involved in chronic stress during pregnancy. We use a multidisciplinary approach that encompasses psychological resources and biophysiological and genetic profiles in order to detect relevant vulnerability biomarkers for chronic stress and the development of PND. To do so, each woman will be involved in the study from her first trimester of pregnancy until 12 months postdelivery. ETHICS AND DISSEMINATION: Ethics approval was obtained from the Ile de France III Ethics Committee, France (2016-A00887-44). We aim to disseminate the findings through international conferences and international peer-reviewed journals. TRIAL REGISTRATION NUMBER: NCT03088319; Pre-results

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

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

Pooled analysis of iron-related genes in Parkinson's disease: Association with transferrin

Pathologic features of Parkinson's disease (PD) include death of dopaminergic neurons in the substantia nigra, presence of α-synuclein containing Lewy bodies, and iron accumulation in PD-related brain regions. The observed iron accumulation may be contributing to PD etiology but it also may be a byproduct of cell death or cellular dysfunction. To elucidate the possible role of iron accumulation in PD, we investigated genetic variation in 16 genes related to iron homeostasis in three case-control studies from the United States, Australia, and France. After screening 90 haplotype tagging single nucleotide polymorphisms (SNPs) within the genes of interest in the US study population, we investigated the five most promising gene regions in two additional independent case-control studies. For the pooled data set (1289 cases, 1391 controls) we observed a protective association (OR. = 0.83, 95% CI: 0.71-0.96) between PD and a haplotype composed of the A allele at rs1880669 and the T allele at rs1049296 in transferrin (TF; GeneID: 7018). Additionally, we observed a suggestive protective association (OR. = 0.87, 95% CI: 0.74-1.02) between PD and a haplotype composed of the G allele at rs10247962 and the A allele at rs4434553 in transferrin receptor 2 (TFR2; GeneID: 7036). We observed no associations in our pooled sample for haplotypes in SLC40A1, CYB561, or HFE. Taken together with previous findings in model systems, our results suggest that TF or a TF- TFR2 complex may have a role in the etiology of PD, possibly through iron misregulation or mitochondrial dysfunction within dopaminergic neurons