Quasispecies distribution of Eigen model

We study sharp peak landscapes (SPL) of Eigen model from a new perspective about how the quasispecies distribute in the sequence space. To analyze the distribution more carefully, we bring forth two tools. One tool is the variance of Hamming distance of the sequences at a given generation. It not only offers us a different avenue for accurately locating the error threshold and illustrates how the configuration of the distribution varies with copying fidelity $q$ in the sequence space, but also divides the copying fidelity into three distinct regimes. The other tool is the similarity network of a certain Hamming distance $d_{0}$, by which we can get a visual and in-depth result about how the sequences distribute. We find that there are several local optima around the center (global optimum) in the distribution of the sequences reproduced near the threshold. Furthermore, it is interesting that the distribution of clustering coefficient $C(k)$ follows lognormal distribution and the curve of clustering coefficient $C$ of the network versus $d_{0}$ appears as linear behavior near the threshold.Comment: 13 pages, 6 figure

The Error and Repair Catastrophes: A Two-Dimensional Phase Diagram in the Quasispecies Model

This paper develops a two gene, single fitness peak model for determining the equilibrium distribution of genotypes in a unicellular population which is capable of genetic damage repair. The first gene, denoted by $\sigma_{via}$, yields a viable organism with first order growth rate constant $k > 1$ if it is equal to some target ``master'' sequence $\sigma_{via, 0}$. The second gene, denoted by $\sigma_{rep}$, yields an organism capable of genetic repair if it is equal to some target ``master'' sequence $\sigma_{rep, 0}$. This model is analytically solvable in the limit of infinite sequence length, and gives an equilibrium distribution which depends on \mu \equiv L\eps , the product of sequence length and per base pair replication error probability, and \eps_r , the probability of repair failure per base pair. The equilibrium distribution is shown to exist in one of three possible ``phases.'' In the first phase, the population is localized about the viability and repairing master sequences. As \eps_r exceeds the fraction of deleterious mutations, the population undergoes a ``repair'' catastrophe, in which the equilibrium distribution is still localized about the viability master sequence, but is spread ergodically over the sequence subspace defined by the repair gene. Below the repair catastrophe, the distribution undergoes the error catastrophe when $\mu$ exceeds \ln k/\eps_r , while above the repair catastrophe, the distribution undergoes the error catastrophe when $\mu$ exceeds $\ln k/f_{del}$, where $f_{del}$ denotes the fraction of deleterious mutations.Comment: 14 pages, 3 figures. Submitted to Physical Review

Evolutionary trajectories in rugged fitness landscapes

We consider the evolutionary trajectories traced out by an infinite population undergoing mutation-selection dynamics in static, uncorrelated random fitness landscapes. Starting from the population that consists of a single genotype, the most populated genotype \textit{jumps} from a local fitness maximum to another and eventually reaches the global maximum. We use a strong selection limit, which reduces the dynamics beyond the first time step to the competition between independent mutant subpopulations, to study the dynamics of this model and of a simpler one-dimensional model which ignores the geometry of the sequence space. We find that the fit genotypes that appear along a trajectory are a subset of suitably defined fitness \textit{records}, and exploit several results from the record theory for non-identically distributed random variables. The genotypes that contribute to the trajectory are those records that are not \textit{bypassed} by superior records arising further away from the initial population. Several conjectures concerning the statistics of bypassing are extracted from numerical simulations. In particular, for the one-dimensional model, we propose a simple relation between the bypassing probability and the dynamic exponent which describes the scaling of the typical evolution time with genome size. The latter can be determined exactly in terms of the extremal properties of the fitness distribution.Comment: Figures in color; minor revisions in tex

Characterization of the role of BGS13 in the secretory mechanism of Pichia pastoris

The methylotrophic yeast Pichia pastoris has been utilized for heterologous protein expression for over 30 years. Because P. pastoris secretes few of its own proteins, the exported recombinant protein is the major polypeptide in the extracellular medium, making purification relatively easy. Unfortunately, some recombinant proteins intended for secretion are retained within the cell. A mutant strain isolated in our laboratory, containing a disruption of the BGS13 gene, displayed elevated levels of secretion for a variety of reporter proteins. The Bgs13 peptide (Bgs13p) is similar to the Saccharomyces cerevisiae protein kinase C 1 protein (Pkc1p), but its specific mode of action is currently unclear. To illuminate differences in the secretion mechanism between the wild-type (wt) strain and the bgs13 strain, we determined that the disrupted bgs13 gene expressed a truncated protein that had reduced protein kinase C activity and a different location in the cell, compared to the wt protein. Because the Pkc1p of baker’s yeast plays a significant role in cell wall integrity, we investigated the sensitivity of the mutant strain’s cell wall to growth antagonists and extraction by dithiothreitol, determining that the bgs13 strain cell wall suffered from inherent structural problems although its porosity was normal. A proteomic investigation of the bgs13 strain secretome and cell wall-extracted peptides demonstrated that, compared to its wt parent, the bgs13 strain also displayed increased release of an array of normally secreted, endogenous proteins, as well as endoplasmic reticulum-resident chaperone proteins, suggesting that Bgs13p helps regulate the unfolded protein response and protein sorting on a global scale. IMPORTANCE The yeast Pichia pastoris is used as a host system for the expression of recombinant proteins. Many of these products, including antibodies, vaccine antigens, and therapeutic proteins such as insulin, are currently on the market or in late stages of development. However, one major weakness is that sometimes these proteins are not secreted from the yeast cell efficiently, which impedes and raises the cost of purification of these vital proteins. Our laboratory has isolated a mutant strain of Pichia pastoris that shows enhanced secretion of many proteins. The mutant produces a modified version of Bgs13p. Our goal is to understand how the change in the Bgs13p function leads to improved secretion. Once the Bgs13p mechanism is illuminated, we should be able to apply this understanding to engineer new P. pastoris strains that efficiently produce and secrete life-saving recombinant proteins, providing medical and economic benefits

Phylogeographic analysis reveals multiple international transmission events have driven the global emergence of Escherichia coli O157:H7

This work was supported by: Scotland by Food Standards Scotland [Grant Number FS102029] and University of Aberdeen; New Zealand, Institute of Environmental Science and Research; Canada, the Public Health Agency of Canada; United States, United States Department of AgriculturePeer reviewedPostprin

GrassPlot v. 2.00 – first update on the database of multi-scale plant diversity in Palaearctic grasslands

Abstract: GrassPlot is a collaborative vegetation-plot database organised by the Eurasian Dry Grassland Group (EDGG) and listed in the Global Index of Vegetation-Plot Databases (GIVD ID EU-00-003). Following a previous Long Database Report (Dengler et al. 2018, Phyto- coenologia 48, 331–347), we provide here the first update on content and functionality of GrassPlot. The current version (GrassPlot v. 2.00) contains a total of 190,673 plots of different grain sizes across 28,171 independent plots, with 4,654 nested-plot series including at least four grain sizes. The database has improved its content as well as its functionality, including addition and harmonization of header data (land use, information on nestedness, structure and ecology) and preparation of species composition data. Currently, GrassPlot data are intensively used for broad-scale analyses of different aspects of alpha and beta diversity in grassland ecosystems

Expression analysis and in silico characterization of intronic long noncoding RNAs in renal cell carcinoma: emerging functional associations

Abstract Background Intronic and intergenic long noncoding RNAs (lncRNAs) are emerging gene expression regulators. The molecular pathogenesis of renal cell carcinoma (RCC) is still poorly understood, and in particular, limited studies are available for intronic lncRNAs expressed in RCC Methods Microarray experiments were performed with custom-designed arrays enriched with probes for lncRNAs mapping to intronic genomic regions. Samples from 18 primary RCC tumors and 11 nontumor adjacent matched tissues were analyzed. Meta-analyses were performed with microarray expression data from three additional human tissues (normal liver, prostate tumor and kidney nontumor samples), and with large-scale public data for epigenetic regulatory marks and for evolutionarily conserved sequences. Results A signature of 29 intronic lncRNAs differentially expressed between RCC and nontumor samples was obtained (false discovery rate (FDR) <5%). A signature of 26 intronic lncRNAs significantly correlated with the RCC five-year patient survival outcome was identified (FDR <5%, p-value ≤0.01). We identified 4303 intronic antisense lncRNAs expressed in RCC, of which 22% were significantly (p <0.05) cis correlated with the expression of the mRNA in the same locus across RCC and three other human tissues. Gene Ontology (GO) analysis of those loci pointed to 'regulation of biological processes’ as the main enriched category. A module map analysis of the protein-coding genes significantly (p <0.05) trans correlated with the 20% most abundant lncRNAs, identified 51 enriched GO terms (p <0.05). We determined that 60% of the expressed lncRNAs are evolutionarily conserved. At the genomic loci containing the intronic RCC-expressed lncRNAs, a strong association (p <0.001) was found between their transcription start sites and genomic marks such as CpG islands, RNA Pol II binding and histones methylation and acetylation. Conclusion Intronic antisense lncRNAs are widely expressed in RCC tumors. Some of them are significantly altered in RCC in comparison with nontumor samples. The majority of these lncRNAs is evolutionarily conserved and possibly modulated by epigenetic modifications. Our data suggest that these RCC lncRNAs may contribute to the complex network of regulatory RNAs playing a role in renal cell malignant transformation.We would like to thank Helder I Nakaya and Tarik A El-Jundi for help with the 44 k microarray experiments and, together with Rodrigo Louro, for helpful discussions. This work was supported by grants from Fundação de Amparo a Pesquisa do Estado de São Paulo (FAPESP) and from Instituto Nacional de Ciência e Tecnologia em Oncogenômica to SVA and EMR, and by fellowships from FAPESP to AAF, ACT, SAVA and VMC. SVA and EMR received research fellowship awards from Conselho Nacional de Desenvolvimento Cientifico e Tecnologico (CNPq), Brazil

Modelling human visual navigation using multi-view scene reconstruction

It is often assumed that humans generate a 3D reconstruction of the environment, either in egocentric or world-based coordinates, but the steps involved are unknown. Here, we propose two reconstruction-based models, evaluated using data from two tasks in immersive virtual reality. We model the observer’s prediction of landmark location based on standard photogrammetric methods and then combine location predictions to compute likelihood maps of navigation behaviour. In one model, each scene point is treated independently in the reconstruction; in the other, the pertinent variable is the spatial relationship between pairs of points. Participants viewed a simple environment from one location, were transported (virtually) to another part of the scene and were asked to navigate back. Error distributions varied substantially with changes in scene layout; we compared these directly with the likelihood maps to quantify the success of the models. We also measured error distributions when participants manipulated the location of a landmark to match the preceding interval, providing a direct test of the landmark-location stage of the navigation models. Models such as this, which start with scenes and end with a probabilistic prediction of behaviour, are likely to be increasingly useful for understanding 3D vision

Pan-cancer Alterations of the MYC Oncogene and Its Proximal Network across the Cancer Genome Atlas

Although theMYConcogene has been implicated incancer, a systematic assessment of alterations ofMYC, related transcription factors, and co-regulatoryproteins, forming the proximal MYC network (PMN),across human cancers is lacking. Using computa-tional approaches, we define genomic and proteo-mic features associated with MYC and the PMNacross the 33 cancers of The Cancer Genome Atlas.Pan-cancer, 28% of all samples had at least one ofthe MYC paralogs amplified. In contrast, the MYCantagonists MGA and MNT were the most frequentlymutated or deleted members, proposing a roleas tumor suppressors.MYCalterations were mutu-ally exclusive withPIK3CA,PTEN,APC,orBRAFalterations, suggesting that MYC is a distinct onco-genic driver. Expression analysis revealed MYC-associated pathways in tumor subtypes, such asimmune response and growth factor signaling; chro-matin, translation, and DNA replication/repair wereconserved pan-cancer. This analysis reveals insightsinto MYC biology and is a reference for biomarkersand therapeutics for cancers with alterations ofMYC or the PMN