Control of time-dependent biological processes by temporally patterned input

Temporal patterning of biological variables, in the form of oscillations and rhythms on many time scales, is ubiquitous. Altering the temporal pattern of an input variable greatly affects the output of many biological processes. We develop here a conceptual framework for a quantitative understanding of such pattern dependence, focusing particularly on nonlinear, saturable, time-dependent processes that abound in biophysics, biochemistry, and physiology. We show theoretically that pattern dependence is governed by the nonlinearity of the input–output transformation as well as its time constant. As a result, only patterns on certain time scales permit the expression of pattern dependence, and processes with different time constants can respond preferentially to different patterns. This has implications for temporal coding and decoding, and allows differential control of processes through pattern. We show how pattern dependence can be quantitatively predicted using only information from steady, unpatterned input. To apply our ideas, we analyze, in an experimental example, how muscle contraction depends on the pattern of motorneuron firing.

Stochastic processes strongly influence HIV-1 evolution during suboptimal protease-inhibitor therapy

It has long been assumed that HIV-1 evolution is best described by deterministic evolutionary models because of the large population size. Recently, however, it was suggested that the effective population size (Ne) may be rather small, thereby allowing chance to influence evolution, a situation best described by a stochastic evolutionary model. To gain experimental evidence supporting one of the evolutionary models, we investigated whether the development of resistance to the protease inhibitor ritonavir affected the evolution of the env gene. Sequential serum samples from five patients treated with ritonavir were used for analysis of the protease gene and the V3 domain of the env gene. Multiple reverse transcription–PCR products were cloned, sequenced, and used to construct phylogenetic trees and to calculate the genetic variation and Ne. Genotypic resistance to ritonavir developed in all five patients, but each patient displayed a unique combination of mutations, indicating a stochastic element in the development of ritonavir resistance. Furthermore, development of resistance induced clear bottleneck effects in the env gene. The mean intrasample genetic variation, which ranged from 1.2% to 5.7% before treatment, decreased significantly (P < 0.025) during treatment. In agreement with these findings, Ne was estimated to be very small (500–15,000) compared with the total HIV-1 RNA copy number. This study combines three independent observations, strong population bottlenecking, small Ne, and selection of different combinations of protease-resistance mutations, all of which indicate that HIV-1 evolution is best described by a stochastic evolutionary model.

Gene genealogies and population variation in plants

Early in the development of plant evolutionary biology, genetic drift, fluctuations in population size, and isolation were identified as critical processes that affect the course of evolution in plant species. Attempts to assess these processes in natural populations became possible only with the development of neutral genetic markers in the 1960s. More recently, the application of historically ordered neutral molecular variation (within the conceptual framework of coalescent theory) has allowed a reevaluation of these microevolutionary processes. Gene genealogies trace the evolutionary relationships among haplotypes (alleles) with populations. Processes such as selection, fluctuation in population size, and population substructuring affect the geographical and genealogical relationships among these alleles. Therefore, examination of these genealogical data can provide insights into the evolutionary history of a species. For example, studies of Arabidopsis thaliana have suggested that this species underwent rapid expansion, with populations showing little genetic differentiation. The new discipline of phylogeography examines the distribution of allele genealogies in an explicit geographical context. Phylogeographic studies of plants have documented the recolonization of European tree species from refugia subsequent to Pleistocene glaciation, and such studies have been instructive in understanding the origin and domestication of the crop cassava. Currently, several technical limitations hinder the widespread application of a genealogical approach to plant evolutionary studies. However, as these technical issues are solved, a genealogical approach holds great promise for understanding these previously elusive processes in plant evolution.

The interspecific mass–density relationship and plant geometry

We present an a priori theoretical framework for the interspecific allometric relationship between stand mass and plant population density. Our model predicts a slope of −\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}\frac{1}{3}\end{equation*}\end{document} between the logarithm of stand mass and the logarithm of stand density, thus conflicting with a previously assumed slope of −½. Our model rests on a heuristic separation of resource-limited living mass and structural mass in the plant body. We point out that because of similar resource requirements among plants of different sizes, a nonzero plant mass–density slope is primarily defined by structural mass. Specifically, the slope is a result of (i) the physical size-dependent relationship between stem width and height, (ii) foliage-dependent demands of conductance, and (iii) the cumulative nature of structural mass. The data support our model, both when the potential sampling bias of taxonomic relatedness is accounted for and when it is not. Independent contrasts analyses show that observed relationships among variables are not significantly different from the assumptions made to build the model or from its a priori predictions. We note that the dependence of the plant mass–density slope on the functions of structural mass provides a cause for the difference from the zero slope found in the animal population mass–density relationship; for the most part, animals do not have a comparable cumulative tissue type.

Different population dynamics of human T cell lymphotropic virus type II in intravenous drug users compared with endemically infected tribes

The phylogeny of human T cell lymphotropic virus type II (HTLV-II) was investigated by using strains isolated from Amerindian and Pygmy tribes, in which the virus is maintained primarily through mother-to-child transmission via breast-feeding, and strains from intravenous drug users (IDUs), in which spread is mainly blood-borne via needle sharing. Molecular clock analysis showed that HTLV-II has two different evolutionary rates with the molecular clock for the virus in IDUs ticking 150–350 times faster than the one in endemically infected tribes: 2.7 × 10−4 compared with 1.71/7.31 × 10−7 nucleotide substitutions per site per year in the long terminal repeat region. This dramatic acceleration of the evolutionary rate seems to be related with the mode of transmission. Mathematical models showed the correlation of these two molecular clocks with an endemic spread of HTLV-II in infected tribes compared with the epidemic spread in IDUs. We also noted a sharp increase in the population size of the virus among IDUs during the last decades probably caused by the worldwide increase in intravenous drug use.

Social transmission of reproductive behavior increases frequency of inherited disorders in a young-expanding population

The observation of high frequencies of certain inherited disorders in the population of Saguenay–Lac Saint Jean can be explained in terms of the variance and the correlation of effective family size (EFS) from one generation to the next. We have shown this effect by using the branching process approach with real demographic data. When variance of EFS is included in the model, despite its profound effect on mutant allele frequency, any mutant introduced in the population never reaches the known carrier frequencies (between 0.035 and 0.05). It is only when the EFS correlation between generations is introduced into the model that we can explain the rise of the mutant alleles. This correlation is described by a c parameter that reflects the dependency of children’s EFS on their parents’ EFS. The c parameter can be considered to reflect social transmission of demographic behavior. We show that such social transmission dramatically reduces the effective population size. This could explain particular distributions in allele frequencies and unusually high frequency of certain inherited disorders in some human populations.

Population dynamics, demographic stochasticity, and the evolution of cooperation

A basic evolutionary problem posed by the Iterated Prisoner’s Dilemma game is to understand when the paradigmatic cooperative strategy Tit-for-Tat can invade a population of pure defectors. Deterministically, this is impossible. We consider the role of demographic stochasticity by embedding the Iterated Prisoner’s Dilemma into a population dynamic framework. Tit-for-Tat can invade a population of defectors when their dynamics exhibit short episodes of high population densities with subsequent crashes and long low density periods with strong genetic drift. Such dynamics tend to have reddened power spectra and temporal distributions of population size that are asymmetric and skewed toward low densities. The results indicate that ecological dynamics are important for evolutionary shifts between adaptive peaks.

The evolution of RNA viruses: A population genetics view

RNA viruses are excellent experimental models for studying evolution under the theoretical framework of population genetics. For a proper justification of this thesis we have introduced some properties of RNA viruses that are relevant for studying evolution. On the other hand, population genetics is a reductionistic theory of evolution. It does not consider or make simplistic assumptions on the transformation laws within and between genotypic and phenotypic spaces. However, such laws are minimized in the case of RNA viruses because the phenotypic space maps onto the genotypic space in a much more linear way than on higher DNA-based organisms. Under experimental conditions, we have tested the role of deleterious and beneficial mutations in the degree of adaptation of vesicular stomatitis virus (VSV), a nonsegmented virus of negative strand. We also have studied how effective population size, initial genetic variability in populations, and environmental heterogeneity shapes the impact of mutations in the evolution of vesicular stomatitis virus. Finally, in an integrative attempt, we discuss pros and cons of the quasispecies theory compared with classic population genetics models for haploid organisms to explain the evolution of RNA viruses.

Protein misfolding and temperature up-shift cause G1 arrest via a common mechanism dependent on heat shock factor in Saccharomyces cerevisiae

Accumulation of misfolded proteins in the cell at high temperature may cause entry into a nonproliferating, heat-shocked state. The imino acid analog azetidine 2-carboxylic acid (AZC) is incorporated into cellular protein competitively with proline and can misfold proteins into which it is incorporated. AZC addition to budding yeast cells at concentrations sufficient to inhibit proliferation selectively activates heat shock factor (HSF). We find that AZC treatment fails to cause accumulation of glycogen and trehalose (Msn2/4-dependent processes) or to induce thermotolerance (a protein kinase C-dependent process). However, AZC-arrested cells can accumulate glycogen and trehalose and can acquire thermotolerance in response to a subsequent heat shock. We find that AZC treatment arrests cells in a viable state and that this arrest is reversible. We find that cells at high temperature or cells deficient in the ubiquitin-conjugating enzymes Ubc4 and Ubc5 are hypersensitive to AZC-induced proliferation arrest. We find that AZC treatment mimics temperature up-shift in arresting cells in G1 and represses expression of CLN1 and CLN2. Mutants with reduced G1 cyclin-Cdc28 activity are hypersensitive to AZC-induced proliferation arrest. Expression of the hyperstable Cln3–2 protein prevents G1 arrest upon AZC treatment and temperature up-shift. Finally, we find that the EXA3–1 mutation, encoding a defective HSF, prevents efficient G1 arrest in response to both temperature up-shift and AZC treatment. We conclude that nontoxic levels of misfolded proteins (induced by AZC treatment or by high temperature) selectively activate HSF, which is required for subsequent G1 arrest.

Investigation of the bottleneck leading to the domestication of maize

Maize (Zea mays ssp. mays) is genetically diverse, yet it is also morphologically distinct from its wild relatives. These two observations are somewhat contradictory: the first observation is consistent with a large historical population size for maize, but the latter observation is consistent with strong, diversity-limiting selection during maize domestication. In this study, we sampled sequence diversity, coupled with simulations of the coalescent process, to study the dynamics of a population bottleneck during the domestication of maize. To do this, we determined the DNA sequence of a 1,400-bp region of the Adh1 locus from 19 individuals representing maize, its presumed progenitor (Z. mays ssp. parviglumis), and a more distant relative (Zea luxurians). The sequence data were used to guide coalescent simulations of population bottlenecks associated with domestication. Our study confirms high genetic diversity in maize—maize contains 75% of the variation found in its progenitor and is more diverse than its wild relative, Z. luxurians—but it also suggests that sequence diversity in maize can be explained by a bottleneck of short duration and very small size. For example, the breadth of genetic diversity in maize is consistent with a founding population of only 20 individuals when the domestication event is 10 generations in length.

Deleterious mutations destabilize ribosomal RNA in endosymbiotic bacteria

In populations that are small and asexual, mutations with slight negative effects on fitness will drift to fixation more often than in large or sexual populations in which they will be eliminated by selection. If such mutations occur in substantial numbers, the combined effects of long-term asexuality and small population size may result in substantial accumulation of mildly deleterious substitutions. Prokaryotic endosymbionts of animals that are transmitted maternally for very long periods are effectively asexual and experience smaller effective population size than their free-living relatives. The contrast between such endosymbionts and related free-living bacteria allows us to test whether a population structure imposing frequent bottlenecks and asexuality does lead to an accumulation of slightly deleterious substitutions. Here we show that several independently derived insect endosymbionts, each with a long history of maternal transmission, have accumulated destabilizing base substitutions in the highly conserved 16S rRNA. Stabilities of Domain I of this subunit are 15–25% lower in endosymbionts than in closely related free-living bacteria. By mapping destabilizing substitutions onto a reconstructed phylogeny, we show that decreased ribosomal stability has evolved separately in each endosymbiont lineage. Our phylogenetic approach allows us to demonstrate statistical significance for this pattern: becoming endosymbiotic predictably results in decreased stability of rRNA secondary structure.

A negative regulator mediates quorum-sensing control of exopolysaccharide production in Pantoea stewartii subsp. stewartii

Classical quorum-sensing (autoinduction) regulation, as exemplified by the lux system of Vibrio fischeri, requires N-acyl homoserine lactone (AHL) signals to stimulate cognate transcriptional activators for the cell density-dependent expression of specific target gene systems. For Pantoea stewartii subsp. stewartii, a bacterial pathogen of sweet corn and maize, the extracellular polysaccharide (EPS) stewartan is a major virulence factor, and its production is controlled by quorum sensing in a population density-dependent manner. Two genes, esaI and esaR, encode essential regulatory proteins for quorum sensing. EsaI is the AHL signal synthase, and EsaR is the cognate gene regulator. esaI, ΔesaR, and ΔesaI-esaR mutations were constructed to establish the regulatory role of EsaR. We report here that strains containing an esaR mutation produce high levels of EPS independently of cell density and in the absence of the AHL signal. Our data indicate that quorum-sensing regulation in P. s. subsp. stewartii, in contrast to most other described systems, uses EsaR to repress EPS synthesis at low cell density, and that derepression requires micromolar amounts of AHL. In addition, derepressed esaR strains, which synthesize EPS constitutively at low cell densities, were significantly less virulent than the wild-type parent. This finding suggests that quorum sensing in P. s. subsp. stewartii may be a mechanism to delay the expression of EPS during the early stages of infection so that it does not interfere with other mechanisms of pathogenesis.

In vitro selection of a 7-methyl-guanosine binding RNA that inhibits translation of capped mRNA molecules

Using systematic evolution of ligands by exponential enrichment (SELEX), an RNA molecule was isolated that displays a 1,000-fold higher affinity for guanosine residues that carry an N-7 methyl group than for nonmethylated guanosine residues. The methylated guanosine residue closely resembles the 5′ terminal cap structure present on all eukaryotic mRNA molecules. The cap-binding RNA specifically inhibited the translation of capped but not uncapped mRNA molecules in cell-free lysates prepared from either human HeLa cells or from Saccharomyces cerevisiae. These findings indicate that the cap-binding RNA will also be useful in studies of other cap-dependent processes such as pre-mRNA splicing and nucleocytoplasmic mRNA transport.

Persistence of neutral polymorphisms in Lake Victoria cichlid fish

Phylogenetic trees for groups of closely related species often have different topologies, depending on the genes used. One explanation for the discordant topologies is the persistence of polymorphisms through the speciation phase, followed by differential fixation of alleles in the resulting species. The existence of transspecies polymorphisms has been documented for alleles maintained by balancing selection but not for neutral alleles. In the present study, transspecific persistence of neutral polymorphisms was tested in the endemic haplochromine species flock of Lake Victoria cichlid fish. Putative noncoding region polymorphisms were identified at four randomly selected nuclear loci and tested on a collection of 12 Lake Victoria species and their putative riverine ancestors. At all loci, the same polymorphism was found to be present in nearly all the tested species, both lacustrine and riverine. Different polymorphisms at these loci were found in cichlids of other East African lakes (Malawi and Tanganyika). The Lake Victoria polymorphisms must have therefore arisen after the flocks now inhabiting the three great lakes diverged from one another, but before the riverine ancestors of the Lake Victoria flock colonized the Lake. Calculations based on the mtDNA clock suggest that the polymorphisms have persisted for about 1.4 million years. To maintain neutral polymorphisms for such a long time, the population size must have remained large throughout the entire period.

RacF1, a Novel Member of the Rho Protein Family in Dictyostelium discoideum, Associates Transiently with Cell Contact Areas, Macropinosomes, and Phagosomes

Using a PCR approach we have isolated racF1, a novel member of the Rho family in Dictyostelium. The racF1 gene encodes a protein of 193 amino acids and is constitutively expressed throughout the Dictyostelium life cycle. Highest identity (94%) was found to a RacF2 isoform, to Dictyostelium Rac1A, Rac1B, and Rac1C (70%), and to Rac proteins of animal species (64–69%). To investigate the role of RacF1 in cytoskeleton-dependent processes, we have fused it at its amino-terminus with green fluorescent protein (GFP) and studied the dynamics of subcellular redistribution using a confocal laser scanning microscope and a double-view microscope system. GFP–RacF1 was homogeneously distributed in the cytosol and accumulated at the plasma membrane, especially at regions of transient intercellular contacts. GFP–RacF1 also localized transiently to macropinosomes and phagocytic cups and was gradually released within <1 min after formation of the endocytic vesicle or the phagosome, respectively. On stimulation with cAMP, no enrichment of GFP–RacF1 was observed in leading fronts, from which it was found to be initially excluded. Cell lines were obtained using homologous recombination that expressed a truncated racF1 gene lacking sequences encoding the carboxyl-terminal region responsible for membrane targeting. These cells displayed normal phagocytosis, endocytosis, and exocytosis rates. Our results suggest that RacF1 associates with dynamic structures that are formed during pinocytosis and phagocytosis. Although RacF1 appears not to be essential, it might act in concert and/or share functions with other members of the Rho family in the regulation of a subset of cytoskeletal rearrangements that are required for these processes.