RNA quality control factors nucleate Clr4/SUV39H and trigger constitutive heterochromatin assembly.

In eukaryotes, the Suv39 family of proteins tri-methylate lysine 9 of histone H3 (H3K9me) to form constitutive heterochromatin. However, how Suv39 proteins are nucleated at heterochromatin is not fully described. In the fission yeast, current models posit that Argonaute1-associated small RNAs (sRNAs) nucleate the sole H3K9 methyltransferase, Clr4/SUV39H, to centromeres. Here, we show that in the absence of all sRNAs and H3K9me, the Mtl1 and Red1 core (MTREC)/PAXT complex nucleates Clr4/SUV39H at a heterochromatic long noncoding RNA (lncRNA) at which the two H3K9 deacetylases, Sir2 and Clr3, also accumulate by distinct mechanisms. Iterative cycles of H3K9 deacetylation and methylation spread Clr4/SUV39H from the nucleation center in an sRNA-independent manner, generating a basal H3K9me state. This is acted upon by the RNAi machinery to augment and amplify the Clr4/H3K9me signal at centromeres to establish heterochromatin. Overall, our data reveal that lncRNAs and RNA quality control factors can nucleate heterochromatin and function as epigenetic silencers in eukaryotes.

Mathematical model for the role of multiple pericentromeric repeats on heterochromatin assembly.

Although the length and constituting sequences for pericentromeric repeats are highly variable across eukaryotes, the presence of multiple pericentromeric repeats is one of the conserved features of the eukaryotic chromosomes. Pericentromeric heterochromatin is often misregulated in human diseases, with the expansion of pericentromeric repeats in human solid cancers. In this article, we have developed a mathematical model of the RNAi-dependent methylation of H3K9 in the pericentromeric region of fission yeast. Our model, which takes copy number as an explicit parameter, predicts that the pericentromere is silenced only if there are many copies of repeats. It becomes bistable or desilenced if the copy number of repeats is reduced. This suggests that the copy number of pericentromeric repeats alone can determine the fate of heterochromatin silencing in fission yeast. Through sensitivity analysis, we identified parameters that favor bistability and desilencing. Stochastic simulation shows that faster cell division and noise favor the desilenced state. These results show the unexpected role of pericentromeric repeat copy number in gene silencing and provide a quantitative basis for how the copy number allows or protects repetitive and unique parts of the genome from heterochromatin silencing, respectively.

YAP1 maintains active chromatin state in head and neck squamous cell carcinomas that promotes tumorigenesis through cooperation with BRD4.

Analysis of The Cancer Genome Atlas and other published data of head and neck squamous cell carcinoma (HNSCC) reveals somatic alterations of the Hippo-YAP pathway in approximately 50% of HNSCC. Better strategies to target the YAP1 transcriptional complex are sought. Here, we show that FAT1, an upstream inhibitor of YAP1, is mutated either by missense or by truncating mutation in 29% of HNSCC. Comprehensive proteomic and drug-screening studies across pan-cancer models confirm that FAT1-mutant HNSCC exhibits selective and higher sensitivity to BRD4 inhibition by JQ1. Epigenomic analysis reveals an active chromatin state in FAT1-mutant HNSCC cells that is driven by the YAP/TAZ transcriptional complex through recruitment of BRD4 to deposit active histone marks, thereby maintaining an oncogenic transcriptional state. This study reveals a detailed cooperative mechanism between YAP1 and BRD4 in HNSCC and suggests a specific therapeutic opportunity for the treatment of this subset of head and neck cancer patients.

Survival in Quiescence Requires the Euchromatic Deployment of Clr4/SUV39H by Argonaute-Associated Small RNAs.

Quiescence (G0) is a ubiquitous stress response through which cells enter reversible dormancy, acquiring distinct properties including reduced metabolism, resistance to stress, and long life. G0 entry involves dramatic changes to chromatin and transcription of cells, but the mechanisms coordinating these processes remain poorly understood. Using the fission yeast, here, we track G0-associated chromatin and transcriptional changes temporally and show that as cells enter G0, their survival and global gene expression programs become increasingly dependent on Clr4/SUV39H, the sole histone H3 lysine 9 (H3K9) methyltransferase, and RNAi proteins. Notably, G0 entry results in RNAi-dependent H3K9 methylation of several euchromatic pockets, prior to which Argonaute1-associated small RNAs from these regions emerge. Overall, our data reveal another function for constitutive heterochromatin proteins (the establishment of the global G0 transcriptional program) and suggest that stress-induced alterations in Argonaute-associated sRNAs can target the deployment of transcriptional regulatory proteins to specific sequences.

Nuclear Noncoding RNAs and Genome Stability.

In modern molecular biology, RNA has emerged as a versatile macromolecule capable of mediating an astonishing number of biological functions beyond its role as a transient messenger of genetic information. The recent discovery and functional analyses of new classes of noncoding RNAs (ncRNAs) have revealed their widespread use in many pathways, including several in the nucleus. This Review focuses on the mechanisms by which nuclear ncRNAs directly contribute to the maintenance of genome stability. We discuss how ncRNAs inhibit spurious recombination among repetitive DNA elements, repress mobilization of transposable elements (TEs), template or bridge DNA double-strand breaks (DSBs) during repair, and direct developmentally regulated genome rearrangements in some ciliates. These studies reveal an unexpected repertoire of mechanisms by which ncRNAs contribute to genome stability and even potentially fuel evolution by acting as templates for genome modification.

PRIMED: PRIMEr database for deleting and tagging all fission and budding yeast genes developed using the open-source genome retrieval script (GRS).

The fission (Schizosaccharomyces pombe) and budding (Saccharomyces cerevisiae) yeasts have served as excellent models for many seminal discoveries in eukaryotic biology. In these organisms, genes are deleted or tagged easily by transforming cells with PCR-generated DNA inserts, flanked by short (50-100 bp) regions of gene homology. These PCR reactions use especially designed long primers, which, in addition to the priming sites, carry homology for gene targeting. Primer design follows a fixed method but is tedious and time-consuming especially when done for a large number of genes. To automate this process, we developed the Python-based Genome Retrieval Script (GRS), an easily customizable open-source script for genome analysis. Using GRS, we created PRIMED, the complete PRIMEr D atabase for deleting and C-terminal tagging genes in the main S. pombe and five of the most commonly used S. cerevisiae strains. Because of the importance of noncoding RNAs (ncRNAs) in many biological processes, we also included the deletion primer set for these features in each genome. PRIMED are accurate and comprehensive and are provided as downloadable Excel files, removing the need for future primer design, especially for large-scale functional analyses. Furthermore, the open-source GRS can be used broadly to retrieve genome information from custom or other annotated genomes, thus providing a suitable platform for building other genomic tools by the yeast or other research communities.

Regulation of histone methylation by noncoding RNAs.

Cells can adapt to their environment and develop distinct identities by rewiring their transcriptional networks to regulate the output of key biological pathways without concomitant mutations to the underlying genes. These alterations, called epigenetic changes, persist stably through mitotic or, in some instances, meiotic cell divisions. In eukaryotes, heritable changes to chromatin structure are a prominent, but not exclusive, mechanism by which epigenetic changes are mediated. These changes are initiated by sequence-specific events, which trigger a cascade of molecular interactions resulting in feedback mechanisms, alterations in chromatin structure, histone posttranslational modifications (PTMs), and ultimately establishment of distinct transcriptional states. In recent years, advances in next generation sequencing have led to the discovery of several novel classes of noncoding RNAs (ncRNAs). In addition to their well-established cytoplasmic roles in posttranscriptional regulation of gene expression, ncRNAs have emerged as key regulators of epigenetic changes via chromatin-dependent mechanisms in organisms ranging from yeast to man. They function by affecting chromatin structure, histone PTMs, and the recruitment of transcriptional activating or repressing complexes. Among histone PTMs, lysine methylation serves as the binding substrate for the recruitment of key protein complexes involved in the regulation of genome architecture, stability, and gene expression. In this review, we will outline the known mechanisms by which ncRNAs of different origins regulate histone methylation, and in doing so contribute to a variety of genome regulatory functions in eukaryotes.

Dynamics of shigellosis epidemics: estimating individual-level transmission and reporting rates from national epidemiologic data sets.

Shigellosis, a diarrheal disease, is endemic worldwide and is responsible for approximately 15,000 laboratory-confirmed cases in the United States every year. However, patients with shigellosis often do not seek medical care. To estimate the burden of shigellosis, we extended time-series susceptible-infected-recovered models to infer epidemiologic parameters from underreported case data. We applied the time-series susceptible-infected-recovered-based inference schemes to analyze the largest surveillance data set of Shigella sonnei in the United States from 1967 to 2007 with county-level resolution. The dynamics of shigellosis transmission show strong annual and multiyear cycles, as well as seasonality. By using the schemes, we inferred individual-level parameters of shigellosis infection, including seasonal transmissibilities and basic reproductive number (R0). In addition, this study provides quantitative estimates of the reporting rate, suggesting that the shigellosis burden in the United States may be more than 10 times the number of laboratory-confirmed cases. Although the estimated reporting rate is generally under 20%, and R0 is generally under 1.5, there is a strong negative correlation between estimates of the reporting rate and R0. Such negative correlations are likely to pose identifiability problems in underreported diseases. We discuss complementary approaches that might further disentangle the true reporting rate and R0.

To lyse or not to lyse: transient-mediated stochastic fate determination in cells infected by bacteriophages.

Cell fate determination is usually described as the result of the stochastic dynamics of gene regulatory networks (GRNs) reaching one of multiple steady-states each of which corresponds to a specific decision. However, the fate of a cell is determined in finite time suggesting the importance of transient dynamics in cellular decision making. Here we consider cellular decision making as resulting from first passage processes of regulatory proteins and examine the effect of transient dynamics within the initial lysis-lysogeny switch of phage λ. Importantly, the fate of an infected cell depends, in part, on the number of coinfecting phages. Using a quantitative model of the phage λ GRN, we find that changes in the likelihood of lysis and lysogeny can be driven by changes in phage co-infection number regardless of whether or not there exists steady-state bistability within the GRN. Furthermore, two GRNs which yield qualitatively distinct steady state behaviors as a function of phage infection number can show similar transient responses, sufficient for alternative cell fate determination. We compare our model results to a recent experimental study of cell fate determination in single cell assays of multiply infected bacteria. Whereas the experimental study proposed a "quasi-independent" hypothesis for cell fate determination consistent with an observed data collapse, we demonstrate that observed cell fate results are compatible with an alternative form of data collapse consistent with a partial gene dosage compensation mechanism. We show that including partial gene dosage compensation at the mRNA level in our stochastic model of fate determination leads to the same data collapse observed in the single cell study. Our findings elucidate the importance of transient gene regulatory dynamics in fate determination, and present a novel alternative hypothesis to explain single-cell level heterogeneity within the phage λ lysis-lysogeny decision switch.

Dynamics of indirectly transmitted infectious diseases with immunological threshold.

There are numerous examples of human pathogens which persist in environmental reservoirs while infectious outbreaks remain rare. In this manuscript, we consider the dynamics of infectious diseases for which the primary mode of transmission is indirect and mediated by contact with a contaminated reservoir. We evaluate the realistic scenario in which the number of ingested pathogens must be above a critical threshold to cause infection in susceptible individuals. This minimal infectious dose is a consequence of the clearance effect of the innate immune system. Infected individuals shed pathogens back into the aquatic reservoir, indirectly increasing the transmittability of the pathogen to the susceptible. Building upon prior works in the study of cholera dynamics, we introduce and analyze a family of reservoir mediated SIR models with a threshold pathogen density for infection. Analyzing this family of models, we show that an outbreak can result from noninfinitesimal introductions of either infected individuals or additional pathogens in the reservoir. We devise two new measures of how likely it is that an environmentally persistent pathogen will cause an outbreak: (i) the minimum fraction of infected individuals; and (ii) the minimum fluctuation size of in-reservoir pathogens. We find an additional control parameter involving the shedding rate of infected individuals, which we term the pathogen enhancement ratio, which determines whether outbreaks lead to epidemics or endemic disease states. Thus, the ultimate outcome of disease is controlled by the strength of fluctuations and the global stability of a nonlinear dynamical system, as opposed to conventional analysis in which disease reflects the linear destabilization of a disease free equilibrium. Our model predicts that in the case of waterborne diseases, suppressing the pathogen density in aquatic reservoirs may be more effective than minimizing the number of infected individuals.

Small-scale copy number variation and large-scale changes in gene expression.

The expression dynamics of interacting genes depends, in part, on the structure of regulatory networks. Genetic regulatory networks include an overrepresentation of subgraphs commonly known as network motifs. In this article, we demonstrate that gene copy number is an omnipresent parameter that can dramatically modify the dynamical function of network motifs. We consider positive feedback, bistable feedback, and toggle switch motifs and show that variation in gene copy number, on the order of a single or few copies, can lead to multiple orders of magnitude change in gene expression and, in some cases, switches in deterministic control. Further, small changes in gene copy number for a 3-gene motif with successive inhibition (the "repressilator") can lead to a qualitative switch in system behavior among oscillatory and equilibrium dynamics. In all cases, the qualitative change in expression is due to the nonlinear nature of transcriptional feedback in which duplicated motifs interact via common pools of transcription factors. We are able to implicitly determine the critical values of copy number which lead to qualitative shifts in system behavior. In some cases, we are able to solve for the sufficient condition for the existence of a bifurcation in terms of kinetic rates of transcription, translation, binding, and degradation. We discuss the relevance of our findings to ongoing efforts to link copy number variation with cell fate determination by viruses, dynamics of synthetic gene circuits, and constraints on evolutionary adaptation.

Collective decision making in bacterial viruses.

For many bacterial viruses, the choice of whether to kill host cells or enter a latent state depends on the multiplicity of coinfection. Here, we present a mathematical theory of how bacterial viruses can make collective decisions concerning the fate of infected cells. We base our theory on mechanistic models of gene regulatory dynamics. Unlike most previous work, we treat the copy number of viral genes as variable. Increasing the viral copy number increases the rate of transcription of viral mRNAs. When viral regulation of cell fate includes nonlinear feedback loops, very small changes in transcriptional rates can lead to dramatic changes in steady-state gene expression. Hence, we prove that deterministic decisions can be reached, e.g., lysis or latency, depending on the cellular multiplicity of infection within a broad class of gene regulatory models of viral decision-making. Comparisons of a parameterized version of the model with molecular studies of the decision structure in the temperate bacteriophage lambda are consistent with our conclusions. Because the model is general, it suggests that bacterial viruses can respond adaptively to changes in population dynamics, and that features of collective decision-making in viruses are evolvable life history traits.