Receptor for advanced glycation end-products (RAGE) mediates phagocytosis in nonprofessional phagocytes.

In mammals, both professional phagocytes and nonprofessional phagocytes (NPPs) can perform phagocytosis. However, limited targets are phagocytosed by NPPs, and thus, the mechanism remains unclear. We find that spores of the yeast Saccharomyces cerevisiae are internalized efficiently by NPPs. Analyses of this phenomenon reveals that RNA fragments derived from cytosolic RNA species are attached to the spore wall, and these fragments serve as ligands to induce spore internalization. Furthermore, we show that a multiligand receptor, RAGE (receptor for advanced glycation end-products), mediates phagocytosis in NPPs. RAGE-mediated phagocytosis is not uniquely induced by spores but is an intrinsic mechanism by which NPPs internalize macromolecules containing RAGE ligands. In fact, artificial particles labeled with polynucleotides, HMGB1, or histone (but not bovine serum albumin) are internalized in NPPs. Our findings provide insight into the molecular basis of phagocytosis by NPPs, a process by which a variety of macromolecules are targeted for internalization.

A conserved expression signature predicts growth rate and reveals cell & lineage-specific differences.

Isogenic cells cultured together show heterogeneity in their proliferation rate. To determine the differences between fast and slow-proliferating cells, we developed a method to sort cells by proliferation rate, and performed RNA-seq on slow and fast proliferating subpopulations of pluripotent mouse embryonic stem cells (mESCs) and mouse fibroblasts. We found that slowly proliferating mESCs have a more naïve pluripotent character. We identified an evolutionarily conserved proliferation-correlated transcriptomic signature that is common to all eukaryotes: fast cells have higher expression of genes for protein synthesis and protein degradation. This signature accurately predicted growth rate in yeast and cancer cells, and identified lineage-specific proliferation dynamics during development, using C. elegans scRNA-seq data. In contrast, sorting by mitochondria membrane potential revealed a highly cell-type specific mitochondria-state related transcriptome. mESCs with hyperpolarized mitochondria are fast proliferating, while the opposite is true for fibroblasts. The mitochondrial electron transport chain inhibitor antimycin affected slow and fast subpopulations differently. While a major transcriptional-signature associated with cell-to-cell heterogeneity in proliferation is conserved, the metabolic and energetic dependency of cell proliferation is cell-type specific.

A rapidly reversible mutation generates subclonal genetic diversity and unstable drug resistance.

Most genetic changes have negligible reversion rates. As most mutations that confer resistance to an adverse condition (e.g., drug treatment) also confer a growth defect in its absence, it is challenging for cells to genetically adapt to transient environmental changes. Here, we identify a set of rapidly reversible drug-resistance mutations in Schizosaccharomyces pombe that are caused by microhomology-mediated tandem duplication (MTD) and reversion back to the wild-type sequence. Using 10,000× coverage whole-genome sequencing, we identify nearly 6,000 subclonal MTDs in a single clonal population and determine, using machine learning, how MTD frequency is encoded in the genome. We find that sequences with the highest-predicted MTD rates tend to generate insertions that maintain the correct reading frame, suggesting that MTD formation has shaped the evolution of coding sequences. Our study reveals a common mechanism of reversible genetic variation that is beneficial for adaptation to environmental fluctuations and facilitates evolutionary divergence.

Uncovering de novo gene birth in yeast using deep transcriptomics.

De novo gene origination has been recently established as an important mechanism for the formation of new genes. In organisms with a large genome, intergenic and intronic regions provide plenty of raw material for new transcriptional events to occur, but little is know about how de novo transcripts originate in more densely-packed genomes. Here, we identify 213 de novo originated transcripts in Saccharomyces cerevisiae using deep transcriptomics and genomic synteny information from multiple yeast species grown in two different conditions. We find that about half of the de novo transcripts are expressed from regions which already harbor other genes in the opposite orientation; these transcripts show similar expression changes in response to stress as their overlapping counterparts, and some appear to translate small proteins. Thus, a large fraction of de novo genes in yeast are likely to co-evolve with already existing genes.

Budding yeast complete DNA synthesis after chromosome segregation begins.

To faithfully transmit genetic information, cells must replicate their entire genome before division. This is thought to be ensured by the temporal separation of replication and chromosome segregation. Here we show that in 20-40% of unperturbed yeast cells, DNA synthesis continues during anaphase, late in mitosis. High cyclin-Cdk activity inhibits DNA synthesis in metaphase, and the decrease in cyclin-Cdk activity during mitotic exit allows DNA synthesis to finish at subtelomeric and some difficult-to-replicate regions. DNA synthesis during late mitosis correlates with elevated mutation rates at subtelomeric regions, including copy number variation. Thus, yeast cells temporally overlap DNA synthesis and chromosome segregation during normal growth, possibly allowing cells to maximize population-level growth rate while simultaneously exploring greater genetic space.

Impact of Chromosome Fusions on 3D Genome Organization and Gene Expression in Budding Yeast.

The three-dimensional (3D) organization of chromosomes can influence transcription. However, the frequency and magnitude of these effects remain debated. To determine how changes in chromosome positioning affect transcription across thousands of genes with minimal perturbation, we characterized nuclear organization and global gene expression in budding yeast containing chromosome fusions. We used computational modeling and single-cell imaging to determine chromosome positions, and integrated these data with genome-wide transcriptional profiles from RNA sequencing. We find that chromosome fusions dramatically alter 3D nuclear organization without leading to strong genome-wide changes in transcription. However, we observe a mild but significant and reproducible increase in the expression of genes displaced away from the periphery. The increase in transcription is inversely proportional to the propensity of a given locus to be at the nuclear periphery; for example, a 10% decrease in the propensity of a gene to reside at the nuclear envelope is accompanied by a 10% increase in gene expression. Modeling suggests that this is due to both deletion of telomeres and to displacement of genes relative to the nuclear periphery. These data suggest that basal transcriptional activity is sensitive to radial changes in gene position, and provide insight into the functional relevance of budding yeast chromosome-level 3D organization in gene expression.

Synchronized replication of genes encoding the same protein complex in fast-proliferating cells.

DNA replication perturbs the dosage balance among genes; at mid-S phase, early-replicating genes have doubled their copies while late-replicating ones have not. Dosage imbalance among genes, especially within members of a protein complex, is toxic to cells. However, the molecular mechanisms that cells use to deal with such imbalance remain not fully understood. Here, we validate at the genomic scale that the dosage between early- and late-replicating genes is imbalanced in HeLa cells. We propose the synchronized replication hypothesis that genes sensitive to stoichiometric relationships will be replicated simultaneously to maintain stoichiometry. In support of this hypothesis, we observe that genes encoding the same protein complex have similar replication timing but mainly in fast-proliferating cells such as embryonic stem cells and cancer cells. We find that the synchronized replication observed in cancer cells, but not in slow-proliferating differentiated cells, is due to convergent evolution during tumorigenesis that restores synchronized replication timing within protein complexes. Taken together, our study reveals that the demand for dosage balance during S phase plays an important role in the optimization of the replication-timing program; this selection is relaxed during differentiation as the cell cycle prolongs and is restored during tumorigenesis as the cell cycle shortens.

Overdosage of Balanced Protein Complexes Reduces Proliferation Rate in Aneuploid Cells.

Cells with complex aneuploidies display a wide range of phenotypic abnormalities. However, the molecular basis for this has been mainly studied in trisomic (2n + 1) and disomic (n + 1) cells. To determine how karyotype affects proliferation in cells with complex aneuploidies, we generated 92 2n + x yeast strains in which each diploid cell has between 3 and 12 extra chromosomes. Genome-wide and, for individual protein complexes, proliferation defects are caused by the presence of protein complexes in which all subunits are balanced at the 3-copy level. Proteomics revealed that over 50% of 3-copy members of imbalanced complexes were expressed at only 2n protein levels, whereas members of complexes in which all subunits are stoichiometrically balanced at 3 copies per cell had 3n protein levels. We validated this finding using orthogonal datasets from yeast and from human cancers. Taken together, our study provides an explanation of how aneuploidy affects phenotype.

Extensive post-transcriptional buffering of gene expression in the response to severe oxidative stress in baker's yeast.

Cells responds to diverse stimuli by changing the levels of specific effector proteins. These changes are usually examined using high throughput RNA sequencing data (RNA-Seq); transcriptional regulation is generally assumed to directly influence protein abundances. However, the correlation between RNA-Seq and proteomics data is in general quite limited owing to differences in protein stability and translational regulation. Here we perform RNA-Seq, ribosome profiling and proteomics analyses in baker's yeast cells grown in rich media and oxidative stress conditions to examine gene expression regulation at various levels. With the exception of a small set of genes involved in the maintenance of the redox state, which are regulated at the transcriptional level, modulation of protein expression is largely driven by changes in the relative ribosome density across conditions. The majority of shifts in mRNA abundance are compensated by changes in the opposite direction in the number of translating ribosomes and are predicted to result in no net change at the protein level. We also identify a subset of mRNAs which is likely to undergo specific translational repression during stress and which includes cell cycle control genes. The study suggests that post-transcriptional buffering of gene expression may be more common than previously anticipated.

Empirical mean-noise fitness landscapes reveal the fitness impact of gene expression noise.

The effects of cell-to-cell variation (noise) in gene expression have proven difficult to quantify because of the mechanistic coupling of noise to mean expression. To independently quantify the effects of changes in mean expression and noise we determine the fitness landscapes in mean-noise expression space for 33 genes in yeast. For most genes, short-lived (noise) deviations away from the expression optimum are nearly as detrimental as sustained (mean) deviations. Fitness landscapes can be classified by a combination of each gene's sensitivity to protein shortage or surplus. We use this classification to explore evolutionary scenarios for gene expression and find that certain landscape topologies can break the mechanistic coupling of mean and noise, thus promoting independent optimization of both properties. These results demonstrate that noise is detrimental for many genes and reveal non-trivial consequences of mean-noise-fitness topologies for the evolution of gene expression systems.

Transcriptomics data of 11 species of yeast identically grown in rich media and oxidative stress conditions.

OBJECTIVE: The objective of this experiment was to identify transcripts in baker's yeast (Saccharomyces cerevisiae) that could have originated from previously non-coding genomic regions, or de novo. We generated this data to be able to compare the transcriptomes of different species of Ascomycota. DATA DESCRIPTION: We generated high-depth RNA sequencing data for 11 species of yeast: Saccharomyces cerevisiae, Saccharomyces paradoxus, Saccharomyces mikatae, Saccharomyces kudriavzevii, Saccharomyces bayanus, Naumovia castelii, Kluyveromyces lactis, Lachancea waltii, Lachancea thermotolerans, Lachancea kluyveri, and Schizosaccharomyces pombe. Using RNA-Seq from yeast grown in rich and oxidative conditions we created genome-guided de novo assemblies of the transcriptomes for each species. We included synthetic spike-in transcripts in each sample to determine the lower limit of detection of the sequencing platform as well as the reliability of our de novo transcriptome assembly pipeline. We subsequently compared the de novo transcripts assemblies to the reference gene annotations and generated assemblies that comprised both annotated and novel transcripts.

An experimental assay of the interactions of amino acids from orthologous sequences shaping a complex fitness landscape.

Characterizing the fitness landscape, a representation of fitness for a large set of genotypes, is key to understanding how genetic information is interpreted to create functional organisms. Here we determined the evolutionarily-relevant segment of the fitness landscape of His3, a gene coding for an enzyme in the histidine synthesis pathway, focusing on combinations of amino acid states found at orthologous sites of extant species. Just 15% of amino acids found in yeast His3 orthologues were always neutral while the impact on fitness of the remaining 85% depended on the genetic background. Furthermore, at 67% of sites, amino acid replacements were under sign epistasis, having both strongly positive and negative effect in different genetic backgrounds. 46% of sites were under reciprocal sign epistasis. The fitness impact of amino acid replacements was influenced by only a few genetic backgrounds but involved interaction of multiple sites, shaping a rugged fitness landscape in which many of the shortest paths between highly fit genotypes are inaccessible.

Single cell functional genomics reveals the importance of mitochondria in cell-to-cell phenotypic variation.

Mutations frequently have outcomes that differ across individuals, even when these individuals are genetically identical and share a common environment. Moreover, individual microbial and mammalian cells can vary substantially in their proliferation rates, stress tolerance, and drug resistance, with important implications for the treatment of infections and cancer. To investigate the causes of cell-to-cell variation in proliferation, we used a high-throughput automated microscopy assay to quantify the impact of deleting >1500 genes in yeast. Mutations affecting mitochondria were particularly variable in their outcome. In both mutant and wild-type cells mitochondrial membrane potential - but not amount - varied substantially across individual cells and predicted cell-to-cell variation in proliferation, mutation outcome, stress tolerance, and resistance to a clinically used anti-fungal drug. These results suggest an important role for cell-to-cell variation in the state of an organelle in single cell phenotypic variation.

Recurrence-based information processing in gene regulatory networks.

Cellular information processing is generally attributed to the complex networks of genes and proteins that regulate cell behavior. It is still unclear, however, what are the main features of those networks that allow a cell to encode and interpret its ever changing environment. Here, we address this question by studying the computational capabilities of the transcriptional regulatory networks of five evolutionary distant organisms. We identify in all cases a cyclic recurrent structure, formed by a small core of genes, that is essential for dynamical encoding and information integration. The recent history of the cell is projected nonlinearly into this recurrent reservoir of nodes, where it is encoded by its transient dynamics, while the rest of the network forms a readout layer devoted to decode and interpret the high-dimensional dynamical state of the recurrent core. In that way, gene regulatory networks act as echo-state networks that perform optimally in standard memory-demanding tasks, with most of their memory residing in the recurrent reservoir. The biological significance of these results is analyzed in the particular case of the bacterium Escherichia coli. Our work thus suggests that recurrent nonlinear dynamics is a key element for the processing of complex time-dependent information by cells.

Reduced intrinsic DNA curvature leads to increased mutation rate.

BACKGROUND: Mutation rates vary across the genome. Many trans factors that influence mutation rates have been identified, as have specific sequence motifs at the 1-7-bp scale, but cis elements remain poorly characterized. The lack of understanding regarding why different sequences have different mutation rates hampers our ability to identify positive selection in evolution and to identify driver mutations in tumorigenesis. RESULTS: Here, we use a combination of synthetic genes and sequences of thousands of isolated yeast colonies to show that intrinsic DNA curvature is a major cis determinant of mutation rate. Mutation rate negatively correlates with DNA curvature within genes, and a 10% decrease in curvature results in a 70% increase in mutation rate. Consistently, both yeast and humans accumulate mutations in regions with small curvature. We further show that this effect is due to differences in the intrinsic mutation rate, likely due to differences in mutagen sensitivity and not due to differences in the local activity of DNA repair. CONCLUSIONS: Our study establishes a framework for understanding the cis properties of DNA sequence in modulating the local mutation rate and identifies a novel causal source of non-uniform mutation rates across the genome.

Poor codon optimality as a signal to degrade transcripts with frameshifts.

Frameshifting errors are common and mRNA quality control pathways, such as nonsense-mediated decay (NMD), exist to degrade these aberrant transcripts. Recent work has shown the existence of a genetic link between NMD and codon-usage mediated mRNA decay. Here we present computational evidence that these pathways are synergic for removing frameshifts.

Promoter Activity Buffering Reduces the Fitness Cost of Misregulation.

Organisms regulate gene expression through changes in the activity of transcription factors (TFs). In yeast, the response of genes to changes in TF activity is generally assumed to be encoded in the promoter. To directly test this assumption, we chose 42 genes and, for each, replaced the promoter with a synthetic inducible promoter and measured how protein expression changes as a function of TF activity. Most genes exhibited gene-specific TF dose-response curves not due to differences in mRNA stability, translation, or protein stability. Instead, most genes have an intrinsic ability to buffer the effects of promoter activity. This can be encoded in the open reading frame and the 3' end of genes and can be implemented by both autoregulatory feedback and by titration of limiting trans regulators. We show experimentally and computationally that, when misexpression of a gene is deleterious, this buffering insulates cells from fitness defects due to misregulation.

The L83L ORF of African swine fever virus strain Georgia encodes for a non-essential gene that interacts with the host protein IL-1ß.

African swine fever virus (ASFV) causes a contagious and frequently lethal disease of pigs causing significant economic consequences to the swine industry. The ASFV genome encodes for more than 150 genes, but only a few of them have been studied in detail. Here we report the characterization of open reading frame L83L which encodes a highly conserved protein across all ASFV isolates. A recombinant ASFV harboring a HA tagged L83L protein was developed (ASFV-G-L83L-HA) and used to demonstrate that L83L is a transiently expressed early virus protein. A recombinant ASFV lacking the L83L gene (ASFV-G-ΔL83L) was developed from the highly virulent field isolate Georgia2007 (ASFV-G) and was used to show that L83L is a non-essential gene. ASFV-G-ΔL83L had similar replication in primary swine macrophage cells when compared to its parental virus ASFV-G. Analysis of host-protein interactions for L83L identified IL-1ß as its host ligand. Experimental infection of domestic pigs showed that ASFV-G-ΔL83L is as virulent as the parental virus ASFV-G.

Promoter architecture determines cotranslational regulation of mRNA.

Information that regulates gene expression is encoded throughout each gene but if different regulatory regions can be understood in isolation, or if they interact, is unknown. Here we measure mRNA levels for 10,000 open reading frames (ORFs) transcribed from either an inducible or constitutive promoter. We find that the strength of cotranslational regulation on mRNA levels is determined by promoter architecture. By using a novel computational genetic screen of 6402 RNA-seq experiments, we identify the RNA helicase Dbp2 as the mechanism by which cotranslational regulation is reduced specifically for inducible promoters. Finally, we find that for constitutive genes, but not inducible genes, most of the information encoding regulation of mRNA levels in response to changes in growth rate is encoded in the ORF and not in the promoter. Thus, the ORF sequence is a major regulator of gene expression, and a nonlinear interaction between promoters and ORFs determines mRNA levels.