Prelife catalysts and replicators

Life is based on replication and evolution. But replication cannot be taken for granted. We must ask what there was prior to replication and evolution. How does evolution begin? We have proposed prelife as a generative system that produces information and diversity in the absence of replication. We model prelife as a binary soup of active monomers that form random polymers. ‘Prevolutionary’ dynamics can have mutation and selection prior to replication. Some sequences might have catalytic activity, thereby enhancing the rates of certain prelife reactions. We study the selection criteria for these prelife catalysts. Their catalytic efficiency must be above certain critical values. We find a maintenance threshold and an initiation threshold. The former is a linear function of sequence length, and the latter is an exponential function of sequence length. Therefore, it is extremely hard to select for prelife catalysts that have long sequences. We compare prelife catalysis with a simple model for replication. Assuming fast template-based elongation reactions, we can show that replicators have selection thresholds that are independent of their sequence length. Our calculation demonstrates the efficiency of replication and provides an explanation of why replication was selected over other forms of prelife catalysis.     

Originator dynamics

We study the origin of evolution. Evolution is based on replication, mutation, and selection. But how does evolution begin? When do chemical kinetics turn into evolutionary dynamics? We propose "prelife" and "prevolution" as the logical precursors of life and evolution. Prelife generates sequences of variable length. Prelife is a generative chemistry that proliferates information and produces diversity without replication. The resulting "prevolutionary dynamics" have mutation and selection. We propose an equation that allows us to investigate the origin of evolution. In one limit, this "originator equation" gives the classical selection equation. In the other limit, we obtain "prelife." There is competition between life and prelife and there can be selection for or against replication. Simple prelife equations with uniform rate constants have the property that longer sequences are exponentially less frequent than shorter ones. But replication can reverse such an ordering. As the replication rate increases, some longer sequences can become more frequent than shorter ones. Thus, replication can lead to "reversals" in the equilibrium portraits. We study these reversals, which mark the transition from prelife to life in our model. If the replication potential exceeds a critical value, then life replicates into existence.     

The distribution of the coalescence time and the number of pairwise nucleotide differences in a model of population divergence or speciation with an initial period of gene flow

This paper is concerned with a model of “isolation with an initial period of migration”, where a panmictic ancestral population split into n descendant populations which exchanged migrants symmetrically at a constant rate for a period of time and subsequently became completely isolated. In the limit as the population split occurred an infinitely long time ago, the model becomes an “isolation after migration” model, describing completely isolated descendant populations which arose from a subdivided ancestral population. The probability density function of the coalescence time of a pair of genes and the probability distribution of the number of pairwise nucleotide differences are derived for both models. Whilst these are theoretical results of interest in their own right, they also give an exact analytical expression for the likelihood, for data consisting of the numbers of nucleotide differences between pairs of DNA sequences where each pair is at a different, independent locus. The behaviour of the distribution of the number of pairwise nucleotide differences under these models is illustrated and compared to the corresponding distributions under the “isolation with migration” and “complete isolation” models. It is shown that the distribution of the number of nucleotide differences between a pair of DNA sequences from different descendant populations in the model of “isolation with an initial period of migration” can be quite different from that under the “isolation with migration model”, even if the average migration rate over time (and hence the total number of migrants) is the same in both scenarios. It is also illustrated how the results can be extended to other demographic scenarios that can be described by a combination of isolated panmictic populations and “symmetric island” models.     

“Protoisochores” in certain archaeal species are formed by replication-associated mutational pressure

This report shows that isochore-like structures can be found not only in warm-blooded animals, some reptiles, fishes and yeast, but also in certain archaeal species. In perfectly shaped isochore-like structures (in “protoisochores”) from Sulfolobus acidocaldarius and Thermofilum pendens genomes the difference in 3GC levels between genes from different “protoisochores” is about 30%. In these archaeal species GC-poor “protoisochores” are situated near the origin of replication, while GC-rich “protoisochores” are situated near the terminus of replication. There is a strong linear dependence between position of a gene and its 3GC level in S. acidocaldarius (an average difference in 3GC per 100?000 base pairs is equal to 3.6%). Detailed analyses of nucleotide usage biases in genes from leading and lagging strands led us to the suggestion that 3GC in genes situated near terminus of replication grows due to higher rates of thymine oxidation producing T to C transitions in lagging strands.     

Sequence variations of Env signal peptide alleles in different clinical stages of HIV infection

The human immunodeficiency virus has been shown to increase its infectivity throughout the course of infection. This virus selection property has been associated with genome mutations and recombinations among virus variants, causing amino acid residue alterations in important viral proteins. In order to explore the contribution of Env signal peptide (Env-sp) to Env glycoprotein expression and its possible relationship to increased virus infectivity observed at late stages of infection, we characterized Env-sp sequences derived from twelve patients at “early” and “late” stages of HIV infection without antiretroviral therapy use. In spite of the remarkable overall similarity between both stages, we observed the deletion of a sequence of neutral and basic residues at the Env-sp amino terminus in virus from early stage specimens and the insertion of basic residues in the hydrophobic region on late-stage viral isolates. The Env-sp sequence alterations may have viral adaptive functions during HIV infection.     

The state transition mechanism—simply depending on light-on and -off in Spirulina platensis

The state transition in cyanobacteria is a long-discussed topic of how the photosynthetic machine regulates the excitation energy distribution in balance between the two photosystems. In the current work, whether the state transition is realized by “mobile phycobilisome (PBS)” or “energy spillover” has been clearly answered by monitoring the spectral responses of the intact cells of the cyanobacterium Spirulina platensis. Firstly, light-induced state transition depends completely on a movement of PBSs toward PSI or PSII while the redox-induced one on not only the “mobile PBS” but also an “energy spillover”. Secondly, the “energy spillover” is triggered by dissociation of PSI trimers into the monomers which specially occurs under a case from light to dark, while the PSI monomers will re-aggregate into the trimers under a case from dark to light, i.e., the PSI oligomerization is reversibly regulated by light switch on and off. Thirdly, PSI oligomerization is regulated by the local H⁺ concentration on the cytosol side of the thylakoid membranes, which in turn is regulated by light switch on and off. Fourthly, PSI oligomerization change is the only mechanism for the “energy spillover”. Thus, it can be concluded that the “mobile PBS” is a common rule for light-induced state transition while the “energy spillover” is only a special case when dark condition is involved.     

Characterisation of physiological and immunological differences between Pacific oysters (Crassostrea gigas) genetically selected for high or low survival to summer mortalities and fed different rations under controlled conditions

Within the framework of a national scientific program named “MORtalités ESTivales de l'huître creuse Crassostrea gigas” (MOREST), a family-based experiment was developed to study the genetic basis of resistance to summer mortality in the Pacific oyster, Crassostrea gigas. As part of the MOREST project, the second generation of three resistant families and two susceptible families were chosen and pooled into two respective groups: “R” and “S”. These two groups of oysters were conditioned for 6 months on two food levels (4% and 12% of oyster soft-tissue dry weight in algal dry weight per day) with a temperature gradient that mimicked the Marennes-Oléron natural cycle during the oyster reproductive period. Oyster mortality remained low for the first two months, but then rapidly increased in July when seawater temperature reached 19 °C and above. Mortality was higher in “S” oysters than in “R” oysters, and also higher in oysters fed the 12% diet than those fed 4%, resulting in a decreasing, relative order in cumulative mortality as follows; 12% “S” > 12% “R” > 4% “S” > 4% “R”. Although the observed mortality rates were lower than those previously observed in the field, the mortality differential between “R” and “S” oysters was similar. Gonadal development, estimated by tissue lipid content, followed a relative order yielding a direct, positive relationship between reproductive effort and mortality as we reported precedently by quantitative histology. Regarding hemocyte parameters, one of the most striking observations was that reactive oxygen species (ROS) production was significantly higher in “S” oysters than in “R” oysters in May and June, regardless of food level. The absence of known environmental stress under these experimental conditions suggests that the ROS increase in “S” oyster could be related to their higher reproductive activity. Finally, a higher increase in hyalinocyte counts was observed for”S” oysters, compared to “R” oysters, in July, just before mortality. Taken together, our results suggest an association of genetically based resistance to summer mortality, reproductive strategy and hemocyte parameters.     

Question 5: Does the RNA-World Still Retain its Appeal After 40 Years of Research?

Forty years after its formulation, the hypothesis of the RNA-World remains rather controversial even though studies of RNA catalysis in cellular processes (for example, in the ubiquitous ribosomal peptide-bond formation) have clearly lent increased plausibility to the idea that an RNA-World existed at some point in the evolution leading to the emergence of cellular life. Indeed, several issues remain that weaken the concept: the synthesis of the RNA monomers under prebiotic conditions, their subsequent, efficient polymerization to yield ribozymes that specifically catalyze their own replication. This communication summarizes existing studies of the RNA polymerization from monomers. In our opinion, the recent developments show that given time plausible answers to some of the issues facing the RNA-World hypothesis will be found. Presented at: International School of Complexity – 4th Course: Basic Questions on the Origins of Life; “Ettore Majorana” Foundation and Centre for Scientific Culture, Erice, Italy, 1–6 October 2006.     

Spatially directed assembly of a heterotetrameric Cre-Lox synapse restricts recombination specificity

The pseudo-fourfold homotetrameric synapse formed by Cre protein and target DNA restricts site-specific recombination to sequences containing dyad-symmetric Cre-binding repeats. Mixtures of engineered altered-specificity Cre monomers can form heterotetramers that recombine nonidentical asymmetric sequences, allowing greater flexibility for target site selection in the genome of interest. However, the variety of tetramers allowed by random subunit association increases the chances of unintended reactivity at nontarget sites. This problem can be circumvented by specifying a unique spatial arrangement of heterotetramer subunits. By reconfiguring intersubunit protein-protein contacts, we directed the assembly of two different Cre monomers, each having a distinct DNA sequence specificity, in an alternating (ABAB) configuration. This designed heterotetramer preferentially recombined a particular pair of asymmetric Lox sites over other pairs, whereas a mixture of freely associating subunits showed little bias. Alone, the engineered monomers had reduced reactivity towards both dyad-symmetric and asymmetric sites. Specificity arose because the organization of Cre-binding repeats of the preferred substrate matched the programmed arrangement of the subunits in the heterotetrameric synapse. When this “spatial matching” principle is applied, Cre-mediated recombination can be directed to asymmetric DNA sequences with greater fidelity.     

Regulatory Roles of the Nuclear Envelope

Roles of the nuclear envelope are considered in the regulation of nuclear protein import, ribonucleoprotein export, and coupling of DNA replication to the cell cycle. First, evidence is discussed that indicates that neutral and acidic amino acids can be important in nuclear localization signals as well as the widely acknowledged basic amino acids. Second, the recognition of nuclear localization signals by their receptor “importin” is discussed, focusing on the different roles of the two subunits of importin. Third, a role for the α subunit of importin in RNP export is considered together with the question of how the direction of traffic through nuclear pores is determined. The final part of this article considers evidence that the nuclear membrane prevents reinitiation of DNA replication in Xenopus eggs, by excluding a “licensing factor” that is essential for DNA replication. Replication licensing in Xenopus appears to involve several proteins including the MCM (minichromosome maintenance) complex and ORC, the origin recognition complex, which must bind before the MCM complex can bind to chromatin.     

Functional Characterization of Nuclear Localization and Export Signals in Hepatitis C Virus Proteins and Their Role in the Membranous Web

The hepatitis C virus (HCV) is a positive strand RNA virus of the Flavivirus family that replicates in the cytoplasm of infected hepatocytes. Previously, several nuclear localization signals (NLS) and nuclear export signals (NES) have been identified in HCV proteins, however, there is little evidence that these proteins travel into the nucleus during infection. We have recently shown that nuclear pore complex (NPC) proteins (termed nucleoporins or Nups) are present in the membranous web and are required during HCV infection. In this study, we identify a total of 11 NLS and NES sequences in various HCV proteins. We show direct interactions between HCV proteins and importin α5 (IPOA5/kapα1), importin β3 (IPO5/kap β3), and exportin 1 (XPO1/CRM1) both in-vitro and in cell culture. These interactions can be disrupted using peptides containing the specific NLS or NES sequences of HCV proteins. Moreover, using a synchronized infection system, we show that these peptides inhibit HCV infection during distinct phases of the HCV life cycle. The inhibitory effects of these peptides place them in two groups. The first group binds IPOA5 and inhibits infection during the replication stage of HCV life cycle. The second group binds IPO5 and is active during both early replication and early assembly. This work delineates the entire life cycle of HCV and the active involvement of NLS sequences during HCV replication and assembly. Given the abundance of NLS sequences within HCV proteins, our previous finding that Nups play a role in HCV infection, and the relocation of the NLS double-GFP reporter in HCV infected cells, this work supports our previous hypothesis that NPC-like structures and nuclear transport factors function in the membranous web to create an environment conducive to viral replication.     

Paramyosin and the filaments of molluscan “catch” muscles : I. Paramyosin: Structure and assembly

Aggregates of paramyosin precipitated with divalent cations show a variety of forms with simple staining patterns and a period of 725 Å in the electron microscope. All the forms may be generated from a basic polar array where the molecules do not bond end-to-end: dark staining regions represent “gap” areas where stain can penetrate the paracrystal; light regions are “overlap” areas where stain is largely excluded. The various polar and dihedral arrays can be accounted for by a superposition of two of the basic polar arrays. This interpretation is supported by the appearance of “corrugation” at the lateral edge of the overlap regions of the fibers, by positive staining experiments, and by the “fringe” length observed at the ends of tactoids.     

Survivor's dilemma: Defend the group or flee?

We consider a survival game of gregarious individuals, in which the aim of the players is survival to reproductive age under predator attacks. The survivor’s dilemma (shortly: SVD) game consists in the following: a group member either surely survives alone by fleeing, while its defensive mate may be killed; or tries to save its mate’s life, risking to get killed. The dilemma is that, in every single attack, fleeing ensures maximal survival probability, but if its mate survives by fighting both, and they remain together, its risk to be killed at the next attack will be lower. We show that, if defense is successful enough, then the one-attack game is a prisoner’s dilemma (PD), where fleeing is the strict ESS. We have additively decomposed the SVD game, according to the survival of the group mate of the focal prey, into two games: the aim of the “collective game” is survival of the group of prey. Counter-wise, the aim of the “hostile game” is survival alone (focal prey survives and its mate is killed by the predator). We obtain the following results: if the attack number is large enough, the multi-attack SVD game is dominated by the “collective game” in the sense that each individual can ensure its own maximal survival probability by maximizing the group survival probability in each attack. In the hostile game, the only strict ESS is the fleeing strategy. In the collective game there are two different cases: either defense is a unique strict ESS, or the collective game is bistable, i.e. fleeing and defense are local strict ESS’s. If defense is the only strict ESS in the collective game, and the attack number is large enough, defense replaces fleeing strategy in the multi-attack SVD game. However, in the bistable case, defense cannot invade into the fleeing population. It is shown that, if the interaction between relatives is frequent enough, than defense can replace fleeing strategy, in spite of the fact that in the well-mixed population the collective game is bistable.     

Taxonomical status and phylogenetic relations between the “thomasi” and “atticus” chromosomal races of the underground vole Microtus thomasi (Rodentia, Arvicolinae)

Laboratory crosses among wild caught individuals of the chromosomal races “atticus” and “thomasi”, were performed to analyze the degree of interracial reproductive isolation. The fertility of the studied specimens was evaluated by taking into consideration the reproductive success, the litter size and performing comparative histological examination of the testicular material. All studied populations were submitted to classical cytogenetic and mitochondrial analysis (cytochrome b gene), providing new evidences to the potential phylogenetic relations and taxonomical status of the two chromosomal races. The previously described “atticus” populations are divided in two genetically distinct, geographically and reproductively isolated lineages (2.9% total and 2.4% net divergence), which probably derived from different glacial refugia of Southern Greece. Here, we suggest that the lineage, consisting of the populations from Attiki and Evia Island, should be distinguished as a valid species, named Microtus atticus, including the two chromosomal races “atticus” and “evia”. On the contrary, the ex-“atticus” populations from North Peloponnesus belong to the same mitochondrial lineage with the other Microtus thomasi populations and should be considered as a chromosomal polymorphism inside the chromosomal race “thomasi”.     

ODE models for oncolytic virus dynamics

Replicating oncolytic viruses are able to infect and lyse cancer cells and spread through the tumor, while leaving normal cells largely unharmed. This makes them potentially useful in cancer therapy, and a variety of viruses have shown promising results in clinical trials. Nevertheless, consistent success remains elusive and the correlates of success have been the subject of investigation, both from an experimental and a mathematical point of view. Mathematical modeling of oncolytic virus therapy is often limited by the fact that the predicted dynamics depend strongly on particular mathematical terms in the model, the nature of which remains uncertain. We aim to address this issue in the context of ODE modeling, by formulating a general computational framework that is independent of particular mathematical expressions. By analyzing this framework, we find some new insights into the conditions for successful virus therapy. We find that depending on our assumptions about the virus spread, there can be two distinct types of dynamics. In models of the first type (the “fast spread” models), we predict that the viruses can eliminate the tumor if the viral replication rate is sufficiently high. The second type of models is characterized by a suboptimal spread (the “slow spread” models). For such models, the simulated treatment may fail, even for very high viral replication rates. Our methodology can be used to study the dynamics of many biological systems, and thus has implications beyond the study of virus therapy of cancers.     

Clonal polymerase chain reaction-single-strand conformational polymorphism analysis: an effective approach for identifying cloned sequences

The amplification of target sequences from genomic DNA can result in more than one amplicon sequence being produced even when highly specific primers are used. Here we present a clonal polymerase chain reaction-single-strand conformational polymorphism (PCR-SSCP) approach for screening cloned amplicons and identifying particular clones prior to sequence determination. Comparison of the PCR-SSCP patterns of the cloned amplicons with the PCR-SSCP patterns observed for the DNA templates from which the clones were derived allows PCR artifacts, different alleles, and even different loci to be differentiated prior to sequencing. Using this approach, the number of clones required for reliable sequence determination is minimized, and complex “mixed” amplicons can be resolved easily, cost-effectively, and reliably.     

Pitfalls in comparisons of genetic distances: a case study of the avian family Acrocephalidae

Genetic distances are increasingly being used for identification and species delimitation, especially since the introduction of “barcoding”. While for phylogenetic inferences great care is generally taken to choose the best-fit evolutionary model, this is usually neglected in calculating genetic distances. Moreover, distances obtained from others than best-fit models, different lengths of sequences, and even different loci are often freely compared. We examined the influence of different methods on calculating genetic distances using mitochondrial cytochrome b sequences for the passerine family Acrocephalidae.We found substantial differences between: (1) corrected distances based on the best-fit model (TrN + Γ) vs. uncorrected p-distances; (2) distances calculated based on different parts of the same gene; and (3) distances calculated using the methods of “complete deletion” vs. “pairwise deletion” for sequences that included uncertain nucleotides. All these methodological differences affected comparisons between species and potential taxonomical conclusions.We suggest that (1) different loci are incomparable. (2) Only perfectly homologous regions (same length, same part of locus) should be compared. (3) In the case of sequences with some uncertain nucleotides, only distances calculated by the method of “complete deletion” are fully comparable. (4) Only distances based on the optimal substitution model should be used. (5) Even within the same locus, corrected genetic distances are unique to the study in which they are calculated, as they are conditional on the particular dataset and model selected for that dataset.     

Universal sequence replication, reversible polymerization and early functional biopolymers: a model for the initiation of prebiotic sequence evolution

Many models for the origin of life have focused on understanding how evolution can drive the refinement of a preexisting enzyme, such as the evolution of efficient replicase activity. Here we present a model for what was, arguably, an even earlier stage of chemical evolution, when polymer sequence diversity was generated and sustained before, and during, the onset of functional selection. The model includes regular environmental cycles (e.g. hydration-dehydration cycles) that drive polymers between times of replication and functional activity, which coincide with times of different monomer and polymer diffusivity. Template-directed replication of informational polymers, which takes place during the dehydration stage of each cycle, is considered to be sequence-independent. New sequences are generated by spontaneous polymer formation, and all sequences compete for a finite monomer resource that is recycled via reversible polymerization. Kinetic Monte Carlo simulations demonstrate that this proposed prebiotic scenario provides a robust mechanism for the exploration of sequence space. Introduction of a polymer sequence with monomer synthetase activity illustrates that functional sequences can become established in a preexisting pool of otherwise non-functional sequences. Functional selection does not dominate system dynamics and sequence diversity remains high, permitting the emergence and spread of more than one functional sequence. It is also observed that polymers spontaneously form clusters in simulations where polymers diffuse more slowly than monomers, a feature that is reminiscent of a previous proposal that the earliest stages of life could have been defined by the collective evolution of a system-wide cooperation of polymer aggregates. Overall, the results presented demonstrate the merits of considering plausible prebiotic polymer chemistries and environments that would have allowed for the rapid turnover of monomer resources and for regularly varying monomer/polymer diffusivities.     

Ovarian parameters and fertility of dairy cows selected for one QTL located on BTA3

Recently, one Quantitative Trait Locus (QTL) of female fertility located on Bos Taurus chromosome 3 (BTA3), QTL-F-Fert-BTA3, has been identified in Holstein breed. It is implied in the success rate after the first AI (AI1) in cow. The failure of pregnancy can be due to several factors involved in the different steps of the reproductive process. The aim of our study was to finely phenotype heifers and primiparous cows selected for their haplotype at the QTL-F-Fert-BTA3. We specifically studied the ovarian follicular dynamic and several fertility parameters. Females carrying the favourable haplotype “fertil+” or unfavourable haplotype “fertil−” were monitored by transrectal ultrasonography during their cycle before the first AI (AI1). Follicular dynamic was similar between the two groups. However, the length of the estrus cycle was shorter in heifers than in primiparous cows and two-wave cycles were shorter than three-wave cycles, regardless of the age and the haplotype. The concentration of plasma anti-Müllerian hormone was correlated with the number of small antral follicles. It was higher in heifers than in primiparous cows, independently of their haplotype. The success rate at the AI1 was significantly higher in “fertil+” than in “fertil−” primiparous cows, 35 d after the AI1 (70% vs 39%). In both haplotypes, pregnancy failure occurred mainly before 21 d after AI1. The commencement of luteal activity after calving was significantly earlier in “fertil+” than in “fertil−” primiparous cows. Calving-AI1 and calving-calving intervals were similar between “fertil+” and “fertil−” primiparous cows. Taken together, “fertil+” and “fertil−” primiparous cows present a difference in the success rate after AI1 that is not explained by variations of ovarian dynamics.     

Requirements for Rapid Plasmid ColE1 Copy Number Adjustments: A Mathematical Model of Inhibition Modes and RNA Turnover Rates

The random distribution of ColE1 plasmids between the daughter cells at cell division introduces large copy number variations. Statistic variation associated with limited copy number in single cells also causes fluctuations to emerge spontaneously during the cell cycle. Efficient replication control out of steady state is therefore important to tame such stochastic effects of small numbers. In the present model, the dynamic features of copy number control are divided into two parts: first, how sharply the replication frequency per plasmid responds to changes in the concentration of the plasmid-coded inhibitor, RNA I, and second, how tightly RNA I and plasmid concentrations are coupled. Single (hyperbolic)- and multiple (exponential)-step inhibition mechanisms are compared out of steady state and it is shown how the response in replication frequency depends on the mode of inhibition. For both mechanisms, sensitivity of inhibition is “bought” at the expense of a rapid turnover of a replication preprimer, RNA II. Conventional, single-step, inhibition kinetics gives a sloppy replication control even at high RNA II turnover rates, whereas multiple-step inhibition has the potential of working with unlimited precision. When plasmid concentration changes rapidly, RNA I must be degraded rapidly to be “up to date” with the change. Adjustment to steady state is drastically impaired when the turnover rate constants of RNA I decrease below certain thresholds, but is basically unaffected for a corresponding increase. Several features of copy number control that are shown to be crucial for the understanding of ColE1-type plasmids still remain to be experimentally characterized. It is shown how steady-state properties reflect dynamics at the heart of regulation and therefore can be used to discriminate between fundamentally different copy number control mechanisms. The experimental tests of the predictions made require carefully planned assays, and some suggestions for suitable experiments arise naturally from the present work. It is also discussed how the presence of the Rom protein may affect dynamic qualities of copy number control.