Publications

The mycobiota are a critical part of the gut microbiome, but hostfungal interactions and specific functional contributions of commensal fungi to host fitness remain incompletely understood. Here, we report the identification of a new fungal commensal, Kazachstania heterogenica var. weizmannii, isolated from murine intestines. K. weizmannii exposure prevented Candida albicans colonization and significantly reduced the commensal C. albicans burden in colonized animals. Following immunosuppression of C. albicans colonized mice, competitive fungal commensalism thereby mitigated fatal candidiasis. Metagenome analysis revealed K. heterogenica or K. weizmannii presence among human commensals. Our results reveal competitive fungal commensalism within the intestinal microbiota, independent of bacteria and immune responses, that could bear potential therapeutic value for the management of C. albicansmediated diseases.

Some drugs increase the mutation rate of their target pathogen, a potentially concerning mechanism as the pathogen might evolve faster toward an undesired phenotype. We suggest a four-step assessment of evolutionary safety for the approval of such treatments.

The mutation rate plays an important role in adaptive evolution. It can be modified by mutator and anti-mutator alleles. Recent empirical evidence hints that the mutation rate may vary among genetically identical individuals: evidence from bacteria suggests that the mutation rate can be affected by expression noise of a DNA repair protein and potentially also by translation errors in various proteins. Importantly, this non-genetic variation may be heritable via a transgenerational epigenetic mode of inheritance, giving rise to a mutator phenotype that is independent from mutator alleles. Here, we investigate mathematically how the rate of adaptive evolution is affected by the rate of mutation rate phenotype switching. We model an asexual population with two mutation rate phenotypes, non-mutator and mutator. An offspring may switch from its parental phenotype to the other phenotype. We find that switching rates that correspond to so-far empirically described non-genetic systems of inheritance of the mutation rate lead to higher rates of adaptation on both artificial and natural fitness landscapes. These switching rates can maintain within the same individuals both a mutator phenotype and intermediary mutations, a combination that facilitates adaptation. Moreover, non-genetic inheritance increases the proportion of mutators in the population, which in turn increases the probability of hitchhiking of the mutator phenotype with adaptive mutations. This in turns facilitates the acquisition of additional adaptive mutations. Our results rationalize recently observed noise in the expression of proteins that affect the mutation rate and suggest that non-genetic inheritance of this phenotype may facilitate evolutionary adaptive processes.

Fitness landscape mapping and the prediction of evolutionary trajectories on these landscapes are major tasks in evolutionary biology research. Evolutionary dynamics is tightly linked to the landscape topography, but this relation is not straightforward. Here, we analyze a fitness landscape of a yeast tRNA gene, previously measured under four different conditions. We find that the wild type allele is sub-optimal, and 810% of its variants are fitter. We rule out the possibilities that the wild type is fittest on average on these four conditions or located on a local fitness maximum. Notwithstanding, we cannot exclude the possibility that the wild type might be fittest in some of the many conditions in the complex ecology that yeast lives at. Instead, we find that the wild type is mutationally robust (\u201cflat\u201d), while more fit variants are typically mutationally fragile. Similar observations of mutational robustness or flatness have been so far made in very few cases, predominantly in viral genomes.

Trypanosoma brucei, the parasite that causes sleeping sickness, cycles between an insect and a mammalian host. However, the effect of RNA modifications such as pseudouridinylation on its ability to survive in these two different host environments is unclear. Here, two genome-wide approaches were applied for mapping pseudouridinylation sites (Ψs) on small nucleolar RNA (snoRNA), 7SL RNA, vtRNA, and tRNAs from T. brucei. We show using HydraPsiSeq and RiboMeth-seq, that the Ψ on C/D snoRNA guiding 2-O-methylation increased the efficiency of the guided modification on its target, rRNA. We found differential levels of Ψs on these ncRNAs in the two life stages (insect host and mammalian host) of the parasite. Furthermore, tRNA isoform abundance and Ψ modifications were characterized in these two life stages demonstrating stage-specific regulation. We conclude that the differential Ψ modifications identified here may contribute to modulating the function of non-coding (nc)RNAs involved in rRNA processing, rRNA modification, protein synthesis, and protein translocation during cycling of the parasite between its two hosts.

The tRNA pool determines the efficiency, throughput, and accuracy of translation. Previous studies have identified dynamic changes in the tRNA (transfer RNA) supply and mRNA (messenger RNA) demand during cancerous proliferation. Yet dynamic changes may also occur during physiologically normal proliferation, and these are less well characterized. We examined the tRNA and mRNA pools of T cells during their vigorous proliferation and differentiation upon triggering their antigen receptor. We observed a global signature of switch in demand for codons at the early proliferation phase of the response, accompanied by corresponding changes in tRNA expression levels. In the later phase, upon differentiation, the response of the tRNA pool relaxed back to the basal level, potentially restraining excessive proliferation. Sequencing of tRNAs allowed us to evaluate their diverse base-modifications. We found that two types of tRNA modifications, wybutosine and ms2t6A, are reduced dramatically during T cell activation. These modifications occur in the anticodon loops of two tRNAs that decode \u201cslippery codons,\u201d which are prone to ribosomal frameshifting. Attenuation of these frameshift-protective modifications is expected to increase the potential for proteome-wide frameshifting during T cell proliferation. Indeed, human cell lines deleted of a wybutosine writer showed increased ribosomal frameshifting, as detected with an HIV gag-pol frameshifting site reporter. These results may explain HIVs specific tropism toward proliferating T cells since it requires ribosomal frameshift exactly on the corresponding codon for infection. The changes in tRNA expression and modifications uncover a layer of translation regulation during T cell proliferation and expose a potential tradeoff between cellular growth and translation fidelity.

Tracing evolutionary processes that lead to fixation of genomic variation in wild bacterial populations is a prime challenge in molecular evolution. In particular, the relative contribution of horizontal gene transfer (HGT) vs. de novo mutations during adaptation to a new environment is poorly understood. To gain a better understanding of the dynamics of HGT and its effect on adaptation, we subjected several populations of competent Bacillus subtilis to a serial dilution evolution on a high-salt-containing medium, either with or without foreign DNA from diverse pre-adapted or naturally salt tolerant species. Following 504 generations of evolution, all populations improved growth yield on the medium. Sequencing of evolved populations revealed extensive acquisition of foreign DNA from close Bacillus donors but not from more remote donors. HGT occurred in bursts, whereby a single bacterial cell appears to have acquired dozens of fragments at once. In the largest burst, close to 2% of the genome has been replaced by HGT. Acquired segments tend to be clustered in integration hotspots. Other than HGT, genomes also acquired spontaneous mutations. Many of these mutations occurred within, and seem to alter, the sequence of flagellar proteins. Finally, we show that, while some HGT fragments could be neutral, others are adaptive and accelerate evolution.

Splicing expands, reshapes, and regulates the transcriptome of eukaryotic organisms. Despite its importance, key questions remain unanswered, including the following: Can splicing evolve when organisms adapt to new challenges? How does evolution optimize inefficiency of introns' splicing and of the splicing machinery? To explore these questions, we evolved yeast cells that were engineered to contain an inefficiently spliced intron inside a gene whose protein product was under selection for an increased expression level. We identified a combination of mutations in Cis (within the gene of interest) and in Trans (in mRNA-maturation machinery). Surprisingly, the mutations in Cis resided outside of known intronic functional sites and improved the intron's splicing efficiency potentially by easing tight mRNA structures. One of these mutations hampered a protein's domain that was not under selection, demonstrating the evolutionary flexibility of multi-domain proteins as one domain functionality was improved at the expense of the other domain. The Trans adaptations resided in two proteins, Npl3 and Gbp2, that bind pre-mRNAs and are central to their maturation. Interestingly, these mutations either increased or decreased the affinity of these proteins to mRNA, presumably allowing faster spliceosome recruitment or increased time before degradation of the pre-mRNAs, respectively. Altogether, our work reveals various mechanistic pathways toward optimizations of intron splicing to ultimately adapt gene expression patterns to novel demands.

The localization of mRNAs encoding secreted/membrane proteins (mSMPs) to the endoplasmic reticulum (ER) likely facilitates the co-translational translocation of secreted proteins. However, studies have shown that mSMP recruitment to the ER in eukaryotes can occur in a manner that is independent of the ribosome, translational control, and the signal recognition particle, although the mechanism remains largely unknown. Here, we identify a cis-acting RNA sequence motif that enhances mSMP localization to the ER and appears to increase mRNA stability, and both the synthesis and secretion of secretome proteins. Termed SECReTE, for secretion-enhancing cis regulatory targeting element, this motif is enriched in mRNAs encoding secretome proteins translated on the ER in eukaryotes and on the inner membrane of prokaryotes. SECReTE consists of >= 10 nucleotide triplet repeats enriched with pyrimidine (C/U) every third base (i.e. NNY, where N = any nucleotide, Y = pyrimidine) and can be present in the untranslated as well as the coding regions of the mRNA. Synonymous mutations that elevate the SECReTE count in a given mRNA (e.g. SUC2, HSP150, and CCW12) lead to an increase in protein secretion in yeast, while a reduction in count led to less secretion and physiological defects. Moreover, the addition of SECReTE to the 3'UTR of an mRNA for an exogenously expressed protein (e.g. GFP) led to its increased secretion from yeast cells. Thus, SECReTE constitutes a novel RNA motif that facilitates ER-localized mRNA translation and protein secretion.

In experimental evolution, scientists evolve organisms in the lab, typically by challenging them to new environmental conditions. How best to evolve a desired trait? Should the challenge be applied abruptly, gradually, periodically, sporadically? Should one apply chemical mutagenesis, and do strains with high innate mutation rate evolve faster? What are ideal population sizes of evolving populations? There are endless strategies, beyond those that can be exposed by individual labs. We therefore arranged a community challenge, Evolthon, in which students and scientists from different labs were asked to evolve Escherichia coli or Saccharomyces cerevisiae for an abiotic stresslow temperature. About 30 participants from around the world explored diverse environmental and genetic regimes of evolution. After a period of evolution in each lab, all strains of each species were competed with one another. In yeast, the most successful strategies were those that used mating, underscoring the importance of sex in evolution. In bacteria, the fittest strain used a strategy based on exploration of different mutation rates. Different strategies displayed variable levels of performance and stability across additional challenges and conditions. This study therefore uncovers principles of effective experimental evolutionary regimens and might prove useful also for biotechnological developments of new strains and for understanding natural strategies in evolutionary arms races between species. Evolthon constitutes a model for community-based scientific exploration that encourages creativity and cooperation.

DNA harbors the blueprint for life. However, the instructions stored in the DNA could be altered at the RNA level before they are executed. One of these processes is RNA editing, which was shown to modify RNA sequences in many organisms. The most abundant modification is the deamination of adenosine (A) into inosine (I). In turn, inosine can be identified as a guanosine (G) by the ribosome and other cellular machineries such as reverse transcriptase. In multicellular organisms, enzymes from the ADAR (adenosine deaminase acting on RNA) family mediate RNA editing in mRNA, whereas enzymes from the ADAT family mediate A-to-I editing on tRNAs. In bacteria however, until recently, only one editing site was described, in tRNA(Arg), but never in mRNA. The tRNA site was shown to be modified by tadA (tRNA specific adenosine deaminase) which is believed to be the ancestral enzyme for the RNA editing family of enzymes. In our recent work, we have shown for the first time, editing on multiple sites in bacterial mRNAs and identified tadA as the enzyme responsible for this editing activity. Focusing on one of the identified targets - the self-killing toxin hokB, we found that editing is physiologically regulated and that it increases protein activity. Here we discuss possible modes of regulation on hokB editing, potential roles of RNA editing in bacteria, possible implications, and future research directions.

Although the genetic code is redundant, synonymous codons for the same amino acid are not used with equal frequencies in genomes, a phenomenon termed "codon usage bias." Previous studies have demonstrated that synonymous changes in a coding sequence can exert significant cis effects on the gene's expression level. However, whether the codon composition of a gene can also affect the translation efficiency of other genes has not been thoroughly explored. To study how codon usage bias influences the cellular economy of translation, we massively converted abundant codons to their rare synonymous counterpart in several highly expressed genes in Escherichia coli. This perturbation reduces both the cellular fitness and the translation efficiency of genes that have high initiation rates and are naturally enriched with the manipulated codon, in agreement with theoretical predictions. Interestingly, we could alleviate the observed phenotypes by increasing the supply of the tRNA for the highly demanded codon, thus demonstrating that the codon usage of highly expressed genes was selected in evolution to maintain the efficiency of global protein translation.

Codon usage bias affects protein translation because tRNAs that recognize synonymous codons differ in their abundance. Although the current dogma states that tRNA expression is exclusively regulated by intrinsic control elements (A- and B-box sequences), we revealed, using a reporter that monitors the levels of individual tRNA genes in Caenorhabditis elegans, that eight tryptophan tRNA genes, 100% identical in sequence, are expressed in different tissues and change their expression dynamically. Furthermore, the expression levels of the sup-7 tRNA gene at day 6 were found to predict the animals lifespan. We discovered that the expression of tRNAs that reside within introns of protein-coding genes is affected by the host genes promoter. Pairing between specific Pol II genes and the tRNAs that are contained in their introns is most likely adaptive, since a genome-wide analysis revealed that the presence of specific intronic tRNAs within specific orthologous genes is conserved across Caenorhabditis species.

Correlation does not imply causation. If two variables, say A and B, are correlated, it could be because A causes B, or that B causes A, or because a third factor affects them both. We suggest that in many cases in biology, the causal link might be bi-directional: A causes B through a fast-acting physio- logical process, while B causes A through a slowly accumulating evolutionary process. Furthermore, many trained biologists tend to consistently focus at first on the fast-acting direction, and overlook the slower process in the opposite direction. We analyse several examples from modern biology that dem- onstrate this bias (codon usage optimality and gene expression, gene duplication and genetic dispens- ability, stem cell division and cancer risk, and the microbiome and host metabolism) and also discuss an example from linguistics. These examples demonstrate mutual effects between the fast physiologi- cal processes and the slow evolutionary ones. We believe that building awareness of inference biases among biologists who tend to prefer one causal direction over another could improve scientific reasoning.

Which came first, the Chicken or the Egg?" We suggest this question is not a paradox. The Modern Synthesis envisions speciation through genetic changes in germ cells via random mutations, an "Egg first" scenario, but perhaps epigenetic inheritance mechanisms can transmit adaptive changes initiated in the soma ("Chicken first"). Reviewers: The article was reviewed by Dr. Eugene Koonin, Dr. Itai Yanai, Dr. Laura Landweber.

p53 is a transcription factor that governs numerous stress response pathways within the cell. Maintaining the right levels of p53 is crucial for cell survival and proper cellular homeostasis. The tight regulation of p53 involves many cellular components, most notably its major negative regulators Mdm2 and Mdm4, which maintain p53 protein amount and activity in tight check. microRNAs (miRNAs) are small non-coding RNAs that target specific mRNAs to translational arrest and degradation. miRNAs are also key components of the normal p53 pathway, joining forces with Mdm2 and Mdm4 to maintain proper p53 activity. Here we review the current knowledge of miRNAs targeting Mdm2 and Mdm4, and their importance in different tissues and in pathological states such as cancer. In addition, we address the role of Alu sequences-highly abundant retroelements spread throughout the human genome, and their impact on gene regulation via the miRNA machinery. Alus occupy a significant portion of genes' 3UTR, and as such they have the potential to impact mRNA regulation. Since Alus are primate-specific, they introduce a new regulatory layer into primate genomes. Alus can influence and alter gene regulation, creating primate-specific cancer-preventive regulatory mechanisms to sustain the transition to longer life span in primates. We review the possible influence of Alu sequences on miRNA functionality in general and specifically within the p53 network.

The p53 pathway is pivotal in tumor suppression. Cellular p53 activity is subject to tight regulation, in which the two related proteins Mdm2 and Mdm4 have major roles. The delicate interplay between the levels of Mdm2, Mdm4 and p53 is crucial for maintaining proper cellular homeostasis. microRNAs (miRNAs) are short non-coding RNAs that downregulate the level and translatability of specific target mRNAs. We report that miR-661, a primate-specific miRNA, can target both Mdm2 and Mdm4 mRNA in a cell type-dependent manner. miR-661 interacts with Mdm2 and Mdm4 RNA within living cells. The inhibitory effect of miR-661 is more prevalent on Mdm2 than on Mdm4. Interestingly, the predicted miR-661 targets in both mRNAs reside mainly within Alu elements, suggesting a primate-specific mechanism for regulatory diversification during evolution. Downregulation of Mdm2 and Mdm4 by miR-661 augments p53 activity and inhibits cell cycle progression in p53-proficient cells. Correspondingly, low miR-661 expression correlates with bad outcome in breast cancers that typically express wild-type p53. In contrast, the miR-661 locus tends to be amplified in tumors harboring p53 mutations, and miR-661 promotes migration of cells derived from such tumors. Thus, miR-661 may either suppress or promote cancer aggressiveness, depending on p53 status.

Motivation: The massive spread of repetitive elements in the human genome presents a substantial challenge to the organism, as such elements may accidentally contain seemingly functional motifs. A striking example is offered by the roughly one million copies of Alu repeats in the genome, of which ∼0.5% reside within genes' untranslated regions (UTRs), presenting ∼30 000 novel potential targets for highly conserved microRNAs (miRNAs). Here, we examine the functionality of miRNA targets within Alu elements in 30UTRs in the human genome. Results: Using a comprehensive dataset of miRNA overexpression assays, we show that mRNAs with miRNA targets within Alus are significantly less responsive to the miRNA effects compared with mRNAs that have the same targets outside Alus. Using Ago2-binding mRNA profiling, we confirm that the miRNA machinery avoids miRNA targets within Alus, as opposed to the highly efficient binding of targets outside Alus. We propose three features that prevent potential miRNA sites within Alus from being recognized by the miRNA machinery: (i) Alu repeats that contain miRNA targets and genuine functional miRNA targets appear to reside in distinct mutually exclusive territories within 3'UTRs; (ii) Alus have tight secondary structure that may limit access to the miRNA machinery; and (iii) A-to-I editing of Alu-derived mRNA sequences may divert miRNA targets. The combination of these features is proposed to allow toleration of Alu insertions into mRNAs. Nonetheless, a subset of miRNA targets within Alus appears not to possess any of the aforementioned features, and thus may represent cases where Alu insertion in the genome has introduced novel functional miRNA targets.

Translation of a gene is assumed to be efficient if the supply of the tRNAs that translate it is high. Yet high-abundance tRNAs are often also at high demand since they correspond to preferred codons in genomes. Thus to fully model translational efficiency one must gauge the supply-to-demand ratio of the tRNAs that are required by the transcriptome at a given time. The tRNAs' supply is often approximated by their gene copy number in the genome. Yet neither the demand for each tRNA nor the extent to which its concentration changes across environmental conditions has been extensively examined. Here we compute changes in the codon usage of the transcriptome across different conditions in several organisms by inspecting conventional mRNA expression data. We find recurring dynamics of codon usage in the transcriptome in multiple stressful conditions. In particular, codons that are translated by rare tRNAs become over-represented in the transcriptome in response to stresses. These results raise the possibility that the tRNA pool might dynamically change upon stress to support efficient translation of stress-transcribed genes. Alternatively, stress genes may be typically translated with low efficiency, presumably due to lack of sufficient evolutionary optimization pressure on their codon usage.

Effective translation of the viral genome during the infection cycle most likely enhances its fitness. In this study, we reveal two different strategies employed by cyanophages, viruses infecting cyanobacteria, to enhance their translation efficiency. Cyanophages of the T7-like Podoviridae family adjust their GC content and codon usage to those of their hosts. In contrast, cyanophages of the T4-like Myoviridae family maintain genomes with low GC content, thus sometimes differing from that of their hosts. By introducing their own specific set of tRNAs, they appear to modulate the tRNA pools of hosts with tRNAs that fit the viral low GC preferred codons. We assessed the possible effects of those viral tRNAs on cyanophages and cyanobacterial genomes using the tRNA adaptation index, which measures the extent to which a given pool of tRNAs translates efficiently particular genes. We found a strong selective pressure to gain and maintain tRNAs that will boost translation of myoviral genes when infecting a high GC host, contrasted by a negligible effect on the host genes. Thus, myoviral tRNAs may represent an adaptive strategy to enhance fitness when infecting high GC hosts, thereby potentially broadening the spectrum of hosts while alleviating the need to adjust global parameters such as GC content for each specific host.

Genetic regulatory circuits are often regarded as precise machines that accurately determine the level of expression of each protein. Most experimental technologies used to measure gene expression levels are incapable of testing and challenging this notion, as they often measure levels averaged over entire populations of cells. Yet, when expression levels are measured at the single cell level of even genetically identical cells, substantial cell-to-cell variation (or \u201cnoise\u201d) may be observed. Sometimes different genes in a given genome may display different levels of noise; even the same gene, expressed under different environmental conditions, may display greater cell-to-cell variability in specific conditions and more tight control in other situations. While at first glance noise may seem to be an undesired property of biological networks, it might be beneficial in some cases. For instance, noise will increase functional heterogeneity in a population of microorganisms facing variable, often unpredictable, environmental changes, increasing the probability that some cells may survive the stress. In that respect, we can speculate that the population is implementing a risk distribution strategy, long before genetic heterogeneity could be acquired. Organisms may have evolved to regulate not only the averaged gene expression levels but also the extent of allowed deviations from such an average, setting it at the desired level for every gene under each specific condition. Here we review the evolving understanding of noise, its molecular underpinnings, and its effect on phenotype and fitness when it can be detrimental, beneficial, or neutral and which regulatory tools eukaryotic cells may use to optimally control it.

Survival in natural habitats selects for microorganisms that are well-adapted to a wide range of conditions. Recent studies revealed that cells evolved innovative response strategies that extend beyond merely sensing a given stimulus and responding to it on encounter. A diversity of microorganisms, including Escherichia coli, Vibrio cholerae, and several yeast species, were shown to use a predictive regulation strategy that uses the appearance of one stimulus as a cue for the likely arrival of a subsequent one. A better understanding of such a predictive strategy requires elucidating the interplay between key biological and environmental forces. Here, we describe a mathematical framework to address this challenge. We base this framework on experimental systems featuring early preparation to either a stress or an exposure to improvement in the growth medium. Our model calculates the fitness advantage originating under each regulation strategy in a given habitat. We conclude that, although a predictive response strategy might by advantageous under some ecologies, its costs might exceed the benefit in others. The combined theoretical-experimental treatment presented here helps assess the potential of natural ecologies to support a predictive behavior.

Background: Early embryos contain mRNA transcripts expressed from two distinct origins; those expressed from the mother's genome and deposited in the oocyte (maternal) and those expressed from the embryo's genome after fertilization (zygotic). The transition from maternal to zygotic control occurs at different times in different animals according to the extent and form of maternal contributions, which likely reflect evolutionary and ecological forces. Maternally deposited transcripts rely on post-transcriptional regulatory mechanisms for precise spatial and temporal expression in the embryo, whereas zygotic transcripts can use both transcriptional and post-transcriptional regulatory mechanisms. The differences in maternal contributions between animals may be associated with gene regulatory changes detectable by the size and complexity of the associated regulatory regions.Results: We have used genomic data to identify and compare maternal and/or zygotic expressed genes from six different animals and find evidence for selection acting to shape gene regulatory architecture in thousands of genes. We find that mammalian maternal genes are enriched for complex regulatory regions, suggesting an increase in expression specificity, while egg-laying animals are enriched for maternal genes that lack transcriptional specificity.Conclusions: We propose that this lack of specificity for maternal expression in egg-laying animals indicates that a large fraction of maternal genes are expressed non-functionally, providing only supplemental nutritional content to the developing embryo. These results provide clear predictive criteria for analysis of additional genomes.

MicroRNAs (miRNAs) appear to be key players in the maintenance of genomic integrity. Recent evidence implies that cancers often avoid miRNA-mediated regulation, and global repression of miRNAs is associated with increased tumorigenicity. Here we suggest that miRNAs are directly involved in the maintenance of genomic integrity through global repression of transposable elements (TEs), whose expression and transposition are well-documented causes of genomic instability in mammalian somatic tissues. Hence, one outcome of the tumor's ability to avoid miRNA-mediated regulation might be the enhancement of genomic instability and mutability due to derepression of TEs. We outline possible mechanisms underlying TE repression by miRNAs, including post-transcriptional silencing and transcriptional silencing through DNA and histone methylation. This hypothesis calls into consideration the need to study the role of miRNAs and the RNAi machinery in the nucleus, and specifically their impact on the maintenance of genomic integrity in the context of cancer.

Single-cell experiments of simple regulatory networks can markedly differ from cell population experiments. Such differences arise from stochastic events in individual cells that are averaged out in cell populations. For instance, while individual cells may show sustained oscillations in the concentrations of some proteins, such oscillations may appear damped in the population average. In this paper we investigate the role of RNA stochastic fluctuations as a leading force to produce a sustained excitatory behavior at the single-cell level. As opposed to some previous models, we build a fully stochastic model of a negative feedback loop that explicitly takes into account the RNA stochastic dynamics. We find that messenger RNA random fluctuations can be amplified during translation and produce sustained pulses of protein expression. Motivated by the recent appreciation of the importance of noncoding regulatory RNAs in post-transcription regulation, we also consider the possibility that a regulatory RNA transcript could bind to the messenger RNA and repress translation. Our findings show that the regulatory transcript helps reducing gene expression variability both at the single-cell level and at the cell population level.

Injury to nerve axons induces diverse responses in neuronal cell bodies, some of which are influenced by the distance from the site of injury. This suggests that neurons have the capacity to estimate the distance of the injury site from their cell body. Recent work has shown that the molecular motor dynein transports importin-mediated retrograde signaling complexes from axonal lesion sites to cell bodies, raising the question whether dynein-based mechanisms enable axonal distance estimations in injured neurons? We used computer simulations to examine mechanisms that may provide nerve cells with dynein-dependent distance assessment capabilities. A multiple-signals model was postulated based on the time delay between the arrival of two or more signals produced at the site of injury-a rapid signal carried by action potentials or similar mechanisms and slower signals carried by dynein. The time delay between the arrivals of these two types of signals should reflect the distance traversed, and simulations of this model show that it can indeed provide a basis for distance measurements in the context of nerve injuries. The analyses indicate that the suggested mechanism can allow nerve cells to discriminate between distances differing by 10% or more of their total axon length, and suggest that dynein-based retrograde signaling in neurons can be utilized for this purpose over different scales of nerves and organisms. Moreover, such a mechanism might also function in synapse to nucleus signaling in uninjured neurons. This could potentially allow a neuron to dynamically sense the relative lengths of its processes on an ongoing basis, enabling appropriate metabolic output from cell body to processes.

The state of the transcriptome reflects a balance between mRNA production and degradation. Yet how these two regulatory arms interact in shaping the kinetics of the transcriptome in response to environmental changes is not known. We subjected yeast to two stresses, one that induces a fast and transient response, and another that triggers a slow enduring response. We then used microarrays following transcriptional arrest to measure genome-wide decay profiles under each condition. We found condition-specific changes in mRNA decay rates and coordination between mRNA production and degradation. In the transient response, most induced genes were surprisingly destabilized, whereas repressed genes were somewhat stabilized, exhibiting counteraction between production and degradation. This strategy can reconcile high steady-state level with short response time among induced genes. In contrast, the stress that induces the slow response displays the more expected behavior, whereby most induced genes are stabilized, and repressed genes are destabilized. Our results show genome-wide interplay between mRNA production and degradation, and that alternative modes of such interplay determine the kinetics of the transcriptome in response to stress.

microRNAs (miRs) are small RNAs that regulate gene expression at the posttranscriptional level. It is anticipated that, in combination with transcription factors (TFs), they span a regulatory network that controls thousands of mammalian genes. Here we set out to uncover local and global architectural features of the mammalian miR regulatory network. Using evolutionarily conserved potential binding sites of miRs in human targets, and conserved binding sites of TFs in promoters, we uncovered two regulation networks. The first depicts combinatorial interactions between pairs of miRs with many shared targets. The network reveals several levels of hierarchy, whereby a few miRs interact with many other lowly connected miR partners. We revealed hundreds of "target hubs" genes, each potentially subject to massive regulation by dozens of miRs. Interestingly, many of these target hub genes are transcription regulators and they are often related to various developmental processes. The second network consists of miR-TF pairs that coregulate large sets of common targets. We discovered that the network consists of several recurring motifs. Most notably, in a significant fraction of the miR-TF coregulators the TF appears to regulate the miR, or to be regulated by the miR, forming a diversity of feed-forward loops. Together these findings provide new insights on the architecture of the combined transcriptional-post transcriptional regulatory network.

This work addresses a central question in medical genetics - the distinction between disease-causing SNPs and neutral variations. Unlike previous studies that focused mainly on coding SNPs, our efforts were centered around variations in regulatory regions and specifically within transcription factor (TF) binding sites. We have compiled a comprehensive collection of genome wide TF binding sites and developed computational measures to estimate the effects of binding site variations on the expression profiles of the regulated genes. Applying these measures to binding sites of known TFs, we were able to make predictions that were in line with published experimental evidence and with structural data on DNA-protein interactions. We attempted to generalize the properties of expression-altering substitutions by accumulating statistics from many substitutions across multiple binding sites. We found that in the yeast genome substitutions that abolish a G or a C are on average more severe than substitutions that abolish an A or a T. This may be attributed to the low GC content of the yeast genome, in which G and C may be important for conferring specificity. We found additional factors that are correlated with the severity of a substitution. Such factors can be integrated in order to create a set of rules for the prioritization of regulatory SNPs according to their disease-causing potential.

Functional redundancies, generated by gene duplications, are highly widespread throughout all known genomes. One consequence of these redundancies is a tremendous increase to the robustness of organisms to mutations and other stresses. Yet, this very robustness also renders redundancy evolutionary unstable, and it is, thus, predicted to have only a transient lifetime. In contrast, numerous reports describe instances of functional overlaps that have been conserved throughout extended evolutionary periods. More interestingly, many such backed-up genes were shown to be transcriptionally responsive to the intactness of their redundant partner and are up-regulated if the latter is mutationally inactivated. By manual inspection of the literature, we have compiled a list of such "responsive backup circuits" in a diverse list of species. Reviewing these responsive backup circuits, we extract recurring principles characterizing their regulation. We then apply modeling approaches to explore further their dynamic properties. Our results demonstrate that responsive backup circuits may function as ideal devices for filtering nongenetic noise from transcriptional pathways and obtaining regulatory precision. We thus challenge the view that such redundancies are simply leftovers of ancient duplications and suggest they are an additional component to the sophisticated machinery of cellular regulation. In this respect, we suggest that compensation for gene loss is merely a side effect of sophisticated design principles using functional redundancy.

Deciphering regulatory events that drive malignant transformation represents a major challenge for systems biology. Here, we analyzed genome-wide transcription profiling of an in vitro cancerous transformation process. We focused on a cluster of genes whose expression levels increased as a function of p53 and p16(INK4A) tumor suppressors inactivation. This cluster predominantly consists of cell cycle genes and constitutes a signature of a diversity of cancers. By linking expression profiles of the genes in the cluster with the dynamic behavior of p53 and p16(INK4A), we identified a promoter architecture that integrates signals from the two tumor suppressive channels and that maps their activity onto distinct levels of expression of the cell cycle genes, which, in turn, correspond to different cellular proliferation rates. Taking components of the mitotic spindle as an example, we experimentally verified our predictions that p53-mediated transcriptional repression of several of these novel targets is dependent on the activities of p21, NFY, and E2F. Our study demonstrates how a well-controlled transformation process allows linking between gene expression, promoter architecture, and activity of upstream signaling molecules.

BACKGROUND: In recent years, intensive computational efforts have been directed towards the discovery of promoter motifs that correlate with mRNA expression profiles. Nevertheless, it is still not always possible to predict steady-state mRNA expression levels based on promoter signals alone, suggesting that other factors may be involved. Other genic regions, in particular 3' UTRs, which are known to exert regulatory effects especially through controlling RNA stability and localization, were less comprehensively investigated, and deciphering regulatory motifs within them is thus crucial. RESULTS: By analyzing 3' UTR sequences and mRNA decay profiles of Saccharomyces cerevisiae genes, we derived a catalog of 53 sequence motifs that may be implicated in stabilization or destabilization of mRNAs. Some of the motifs correspond to known RNA-binding protein sites, and one of them may act in destabilization of ribosome biogenesis genes during stress response. In addition, we present for the first time a catalog of 23 motifs associated with subcellular localization. A significant proportion of the 3' UTR motifs is highly conserved in orthologous yeast genes, and some of the motifs are strikingly similar to recently published mammalian 3' UTR motifs. We classified all genes into those regulated only at transcription initiation level, only at degradation level, and those regulated by a combination of both. Interestingly, different biological functionalities and expression patterns correspond to such classification. CONCLUSION: The present motif catalogs are a first step towards the understanding of the regulation of mRNA degradation and subcellular localization, two important processes which--together with transcription regulation--determine the cell transcriptome.

Ohno [Ohno, S. (1970) in Evolution by Gene Duplication, Springer, New York] proposed that gene duplication with subsequent divergence of paralogs could be a major force in the evolution of new gene functions. In practice the functional differences between closely related homologues produced by duplications can be subtle and difficult to separate experimentally. Here we show that DNA microarrays can distinguish the functions of two closely related homologues from the yeast Saccharomyces cerevisiae, Yap1p and Yap2p. Although Yap1p and Yap2p are both bZIP transcription factors involved in multiple stress responses and are 88% identical in their DNA binding domains, our work shows that these proteins activate nonoverlapping sets of genes. Yap1p controls a set of genes involved in detoxifying the effects of reactive oxygen species, whereas Yap2p controls a set of genes over represented for the function of stabilizing proteins. In addition we show that the binding sites in the promoters of the Yap1p-dependent genes differ from the sites in the promoters of Yap2p-dependent genes and we validate experimentally that these differences are important for regulation by Yap1p. We conclude that while Yap1p and Yap2p may have some overlapping functions they are clearly not redundant and, more generally, that DNA microarray analysis will be an important tool for distinguishing the functions of the large numbers of highly conserved genes found in all eukaryotic genomes.

The olfactory receptor (OR) subgenome harbors the largest known gene family in mammals, disposed in clusters on numerous chromosomes. One of the best characterized OR clusters, located at human chromosome 17p13.3, has previously been studied by us in human and in other primates, revealing a conserved set of 17 OR genes. Here, we report the identification of a syntenic OR cluster in the mouse and the partial DNA sequence of many of its OR genes. A probe for the mouse M5 gene, orthologous to one of the OR genes in the human cluster (OR17-25), was used to isolate six PAC clones, all mapping by in situ hybridization to mouse chromosome 11B3-11B5, a region of shared synteny with human chromosome 17p13.3. Thirteen mouse OR sequences amplified and sequenced from these PACs allowed us to construct a putative physical map of the OR gene cluster at the mouse Olfr1 locus. Several points of evidence, including a strong similarity in subfamily composition and at least four cases of gene orthology, suggest that the mouse Olfr1 and the human 17p13.3 clusters are orthologous. A detailed comparison of the OR sequences within the two clusters helps trace their independent evolutionary history in the two species. Two types of evolutionary scenarios are discerned: cases of "true orthologous genes" in which high sequence similarity suggests a shared conserved function, as opposed to instances in which orthologous genes may have undergone independent diversification in the realm of "free reign" repertoire expansion.

The olfactory receptor (OR) subgenome harbors the largest known gene family in mammals, disposed in clusters on numerous chromosomes. We have carried out a comparative evolutionary analysis of the best characterized genomic OR gene cluster, on human chromosome 17p13. Fifteen orthologs from chimpanzee (localized to chromosome 19p15), as well as key OR counterparts from other primates, have been identified and sequenced. Comparison among orthologs and paralogs revealed a multiplicity of gene conversion events, which occurred exclusively within OR subfamilies. These appear to lead to segment shuffling in the odorant binding site, an evolutionary process reminiscent of somatic combinatorial diversification in the immune system. We also demonstrate that the functional mammalian OR repertoire has undergone a rapid decline in the past 10 million years: while for the common ancestor of all great apes an intact OR cluster is inferred, in present-day humans and great apes the cluster includes nearly 40% pseudogenes.

In order to elicit an olfactory response, a substance has to partition into the gas phase and diffuse into the nose. Such odorant molecules, usually low molecular-mass hydrophobic compounds, encounter the ciliated endings of sensory neuronal dendrites, which protrude into a mucus layer at the surface of the olfactory epithelium in the nasal cavity. Embedded in the membranes of such cilia are olfactory receptor (OR) proteins, which recognize odorants and elicit a transduction cascade that underlies the nerve cell response. The sensory axons project to the olfactory bulb in the brain, where they converge into synaptic structures called glomeruli. The specific convergence patterns of olfactory axons, which depend on OR expression, provide a model system for neuronal network development. Here, initial processing of odour information occurs, which is followed by additional analysis in higher olfactory brain centres.

A Graded Autocatalysis Replication Domain (GARD) model is proposed, which provides a rigorous kinetic analysis of simple chemical sets that manifest mutual catalysis. It is shown that catalytic closure can sustain self-replication up to a critical dilution rate, λ_c, related to the graded extent of mutual catalysis. We explore the behaviour of vesicles containing GARD species whose mutual catalysis is governed by a previously published statistical distribution. In the population thus generated, some GARD vesicles display a significantly higher replication efficiency than most others. GARD thus represents a simple model for primordial chemical selection of mutually catalytic sets.