Highly-automated, high-throughput replication of yeast-based logic circuit design assessments.

We describe an experimental campaign that replicated the performance assessment of logic gates engineered into cells of Saccharomyces cerevisiae by Gander et al. Our experimental campaign used a novel high-throughput experimentation framework developed under Defense Advanced Research Projects Agency's Synergistic Discovery and Design program: a remote robotic lab at Strateos executed a parameterized experimental protocol. Using this protocol and robotic execution, we generated two orders of magnitude more flow cytometry data than the original experiments. We discuss our results, which largely, but not completely, agree with the original report and make some remarks about lessons learned. Graphical Abstract.

Prediction of whole-cell transcriptional response with machine learning.

MOTIVATION: Applications in synthetic and systems biology can benefit from measuring whole-cell response to biochemical perturbations. Execution of experiments to cover all possible combinations of perturbations is infeasible. In this paper, we present the host response model (HRM), a machine learning approach that maps response of single perturbations to transcriptional response of the combination of perturbations. RESULTS: The HRM combines high-throughput sequencing with machine learning to infer links between experimental context, prior knowledge of cell regulatory networks, and RNASeq data to predict a gene's dysregulation. We find that the HRM can predict the directionality of dysregulation to a combination of inducers with an accuracy of >90% using data from single inducers. We further find that the use of prior, known cell regulatory networks doubles the predictive performance of the HRM (an R2 from 0.3 to 0.65). The model was validated in two organisms, Escherichia coli and Bacillus subtilis, using new experiments conducted after training. Finally, while the HRM is trained with gene expression data, the direct prediction of differential expression makes it possible to also conduct enrichment analyses using its predictions. We show that the HRM can accurately classify >95% of the pathway regulations. The HRM reduces the number of RNASeq experiments needed as responses can be tested in silico prior to the experiment. AVAILABILITY AND IMPLEMENTATION: The HRM software and tutorial are available at https://github.com/sd2e/CDM and the configurable differential expression analysis tools and tutorials are available at https://github.com/SD2E/omics_tools. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The biochemical basis of mitochondrial dysfunction in Zellweger Spectrum Disorder.

Peroxisomal biogenesis disorders (PBDs) are genetic disorders of peroxisome biogenesis and metabolism that are characterized by profound developmental and neurological phenotypes. The most severe class of PBDs-Zellweger spectrum disorder (ZSD)-is caused by mutations in peroxin genes that result in both non-functional peroxisomes and mitochondrial dysfunction. It is unclear, however, how defective peroxisomes contribute to mitochondrial impairment. In order to understand the molecular basis of this inter-organellar relationship, we investigated the fate of peroxisomal mRNAs and proteins in ZSD model systems. We found that peroxins were still expressed and a subset of them accumulated on the mitochondrial membrane, which resulted in gross mitochondrial abnormalities and impaired mitochondrial metabolic function. We showed that overexpression of ATAD1, a mitochondrial quality control factor, was sufficient to rescue several aspects of mitochondrial function in human ZSD fibroblasts. Together, these data suggest that aberrant peroxisomal protein localization is necessary and sufficient for the devastating mitochondrial morphological and metabolic phenotypes in ZSDs.

The genetic basis of aneuploidy tolerance in wild yeast.

Aneuploidy is highly detrimental during development yet common in cancers and pathogenic fungi - what gives rise to differences in aneuploidy tolerance remains unclear. We previously showed that wild isolates of Saccharomyces cerevisiae tolerate chromosome amplification while laboratory strains used as a model for aneuploid syndromes do not. Here, we mapped the genetic basis to Ssd1, an RNA-binding translational regulator that is functional in wild aneuploids but defective in laboratory strain W303. Loss of SSD1 recapitulates myriad aneuploidy signatures previously taken as eukaryotic responses. We show that aneuploidy tolerance is enabled via a role for Ssd1 in mitochondrial physiology, including binding and regulating nuclear-encoded mitochondrial mRNAs, coupled with a role in mitigating proteostasis stress. Recapitulating ssd1Δ defects with combinatorial drug treatment selectively blocked proliferation of wild-type aneuploids compared to euploids. Our work adds to elegant studies in the sensitized laboratory strain to present a mechanistic understanding of eukaryotic aneuploidy tolerance.

The RNA Helicase DDX6 Controls Cellular Plasticity by Modulating P-Body Homeostasis.

Post-transcriptional mechanisms have the potential to influence complex changes in gene expression, yet their role in cell fate transitions remains largely unexplored. Here, we show that suppression of the RNA helicase DDX6 endows human and mouse primed embryonic stem cells (ESCs) with a differentiation-resistant, "hyper-pluripotent" state, which readily reprograms to a naive state resembling the preimplantation embryo. We further demonstrate that DDX6 plays a key role in adult progenitors where it controls the balance between self-renewal and differentiation in a context-dependent manner. Mechanistically, DDX6 mediates the translational suppression of target mRNAs in P-bodies. Upon loss of DDX6 activity, P-bodies dissolve and release mRNAs encoding fate-instructive transcription and chromatin factors that re-enter the ribosome pool. Increased translation of these targets impacts cell fate by rewiring the enhancer, heterochromatin, and DNA methylation landscapes of undifferentiated cell types. Collectively, our data establish a link between P-body homeostasis, chromatin organization, and stem cell potency.

Reduced MEK inhibition preserves genomic stability in naive human embryonic stem cells.

Human embryonic stem cells (hESCs) can be captured in a primed state in which they resemble the postimplantation epiblast, or in a naive state where they resemble the preimplantation epiblast. Naive-cell-specific culture conditions allow the study of preimplantation development ex vivo but reportedly lead to chromosomal abnormalities, which compromises their utility in research and potential therapeutic applications. Although MEK inhibition is essential for the naive state, here we show that reduced MEK inhibition facilitated the establishment and maintenance of naive hESCs that retained naive-cell-specific features, including global DNA hypomethylation, HERVK expression, and two active X chromosomes. We further show that hESCs cultured under these modified conditions proliferated more rapidly; accrued fewer chromosomal abnormalities; and displayed changes in the phosphorylation levels of MAPK components, regulators of DNA damage/repair, and cell cycle. We thus provide a simple modification to current methods that can enable robust growth and reduced genomic instability in naive hESCs.

ACP Acylation Is an Acetyl-CoA-Dependent Modification Required for Electron Transport Chain Assembly.

The electron transport chain (ETC) is an important participant in cellular energy conversion, but its biogenesis presents the cell with numerous challenges. To address these complexities, the cell utilizes ETC assembly factors, which include the LYR protein family. Each member of this family interacts with the mitochondrial acyl carrier protein (ACP), the scaffold protein upon which the mitochondrial fatty acid synthesis (mtFAS) pathway builds fatty acyl chains from acetyl-CoA. We demonstrate that the acylated form of ACP is an acetyl-CoA-dependent allosteric activator of the LYR protein family used to stimulate ETC biogenesis. By tuning ETC assembly to the abundance of acetyl-CoA, which is the major fuel of the TCA cycle and ETC, this system could provide an elegant mechanism for coordinating the assembly of ETC complexes with one another and with substrate availability.

Nudt21 Controls Cell Fate by Connecting Alternative Polyadenylation to Chromatin Signaling.

Cell fate transitions involve rapid gene expression changes and global chromatin remodeling, yet the underlying regulatory pathways remain incompletely understood. Here, we identified the RNA-processing factor Nudt21 as a novel regulator of cell fate change using transcription-factor-induced reprogramming as a screening assay. Suppression of Nudt21 enhanced the generation of induced pluripotent stem cells, facilitated transdifferentiation into trophoblast stem cells, and impaired differentiation of myeloid precursors and embryonic stem cells, suggesting a broader role for Nudt21 in cell fate change. We show that Nudt21 directs differential polyadenylation of over 1,500 transcripts in cells acquiring pluripotency, although only a fraction changed protein levels. Remarkably, these proteins were strongly enriched for chromatin regulators, and their suppression neutralized the effect of Nudt21 during reprogramming. Collectively, our data uncover Nudt21 as a novel post-transcriptional regulator of cell fate and establish a direct, previously unappreciated link between alternative polyadenylation and chromatin signaling.

Accelerating high-dimensional clustering with lossless data reduction.

MOTIVATION: For cluster analysis, high-dimensional data are associated with instability, decreased classification accuracy and high-computational burden. The latter challenge can be eliminated as a serious concern. For applications where dimension reduction techniques are not implemented, we propose a temporary transformation which accelerates computations with no loss of information. The algorithm can be applied for any statistical procedure depending only on Euclidean distances and can be implemented sequentially to enable analyses of data that would otherwise exceed memory limitations. RESULTS: The method is easily implemented in common statistical software as a standard pre-processing step. The benefit of our algorithm grows with the dimensionality of the problem and the complexity of the analysis. Consequently, our simple algorithm not only decreases the computation time for routine analyses, it opens the door to performing calculations that may have otherwise been too burdensome to attempt. AVAILABILITY AND IMPLEMENTATION: R, Matlab and SAS/IML code for implementing lossless data reduction is freely available in the Appendix. CONTACT: obrienj@hms.harvard.edu.

A unique ecological niche fosters hybridization of oak-tree and vineyard isolates of Saccharomyces cerevisiae.

Differential adaptation to distinct niches can restrict gene flow and promote population differentiation within a species. However, in some cases the distinction between niches can collapse, forming a hybrid niche with features of both environments. We previously reported that distinctions between vineyards and oak soil present an ecological barrier that restricts gene flow between lineages of Saccharomyces cerevisiae. Vineyard isolates are tolerant to stresses associated with grapes while North American oak strains are particularly tolerant to freeze-thaw cycles. Here, we report the isolation of S. cerevisiae strains from Wisconsin cherry trees, which display features common to vineyards (e.g. high sugar concentrations) and frequent freeze-thaw cycles. Genome sequencing revealed that the isolated strains are highly heterozygous and represent recent hybrids of the oak × vineyard lineages. We found that the hybrid strains are phenotypically similar to vineyard strains for some traits, but are more similar to oak strains for other traits. The cherry strains were exceptionally good at growing in cherry juice, raising the possibility that they have adapted to this niche. We performed transcriptome profiling in cherry, oak and vineyard strains and show that the cherry-tree hybrids display vineyard-like or oak-like expression, depending on the gene sets, and in some cases, the expression patterns linked back to shared stress tolerances. Allele-specific expression in these natural hybrids suggested concerted cis-regulatory evolution at sets of functionally regulated genes. Our results raise the possibility that hybridization of the two lineages provides a genetic solution to the thriving in this unique niche.

Ecological and Genetic Barriers Differentiate Natural Populations of Saccharomyces cerevisiae.

How populations that inhabit the same geographical area become genetically differentiated is not clear. To investigate this, we characterized phenotypic and genetic differences between two populations of Saccharomyces cerevisiae that in some cases inhabit the same environment but show relatively little gene flow. We profiled stress sensitivity in a group of vineyard isolates and a group of oak-soil strains and found several niche-related phenotypes that distinguish the populations. We performed bulk-segregant mapping on two of the distinguishing traits: The vineyard-specific ability to grow in grape juice and oak-specific tolerance to the cell wall damaging drug Congo red. To implicate causal genes, we also performed a chemical genomic screen in the lab-strain deletion collection and identified many important genes that fell under quantitative trait loci peaks. One gene important for growth in grape juice and identified by both the mapping and the screen was SSU1, a sulfite-nitrite pump implicated in wine fermentations. The beneficial allele is generated by a known translocation that we reasoned may also serve as a genetic barrier. We found that the translocation is prevalent in vineyard strains, but absent in oak strains, and presents a postzygotic barrier to spore viability. Furthermore, the translocation was associated with a fitness cost to the rapid growth rate seen in oak-soil strains. Our results reveal the translocation as a dual-function locus that enforces ecological differentiation while producing a genetic barrier to gene flow in these sympatric populations.

Genetic variation in senescence marker protein-30 is associated with natural variation in cold tolerance in Drosophila.

A comprehensive understanding of the genetic basis of phenotypic adaptation in nature requires the identification of the functional allelic variation underlying adaptive phenotypes. The manner in which organisms respond to temperature extremes is an adaptation in many species. In the current study, we investigate the role of molecular variation in senescence marker protein-30 (Smp-30) on natural phenotypic variation in cold tolerance in Drosophila melanogaster. Smp-30 encodes a product that is thought to be involved in the regulation of Ca2+ ion homeostasis and has been shown previously to be differentially expressed in response to cold stress. Thus, we sought to assess whether molecular variation in Smp-30 was associated with natural phenotypic variation in cold tolerance in a panel of naturally derived inbred lines from a population in Raleigh, North Carolina. We identified four non-coding polymorphisms that were strongly associated with natural phenotypic variation in cold tolerance. Interestingly, two polymorphisms that were in close proximity to one another (2 bp apart) exhibited opposite phenotypic effects. Consistent with the maintenance of a pair of antagonistically acting polymorphisms, tests of molecular evolution identified a significant excess of maintained variation in this region, suggesting balancing selection is acting to maintain this variation. These results suggest that multiple mutations in non-coding regions can have significant effects on phenotypic variation in adaptive traits within natural populations, and that balancing selection can maintain polymorphisms with opposite effects on phenotypic variation.