Robust genetic codes enhance protein evolvability.

The standard genetic code defines the rules of translation for nearly every life form on Earth. It also determines the amino acid changes accessible via single-nucleotide mutations, thus influencing protein evolvability-the ability of mutation to bring forth adaptive variation in protein function. One of the most striking features of the standard genetic code is its robustness to mutation, yet it remains an open question whether such robustness facilitates or frustrates protein evolvability. To answer this question, we use data from massively parallel sequence-to-function assays to construct and analyze 6 empirical adaptive landscapes under hundreds of thousands of rewired genetic codes, including those of codon compression schemes relevant to protein engineering and synthetic biology. We find that robust genetic codes tend to enhance protein evolvability by rendering smooth adaptive landscapes with few peaks, which are readily accessible from throughout sequence space. However, the standard genetic code is rarely exceptional in this regard, because many alternative codes render smoother landscapes than the standard code. By constructing low-dimensional visualizations of these landscapes, which each comprise more than 16 million mRNA sequences, we show that such alternative codes radically alter the topological features of the network of high-fitness genotypes. Whereas the genetic codes that optimize evolvability depend to some extent on the detailed relationship between amino acid sequence and protein function, we also uncover general design principles for engineering nonstandard genetic codes for enhanced and diminished evolvability, which may facilitate directed protein evolution experiments and the bio-containment of synthetic organisms, respectively.

Specificity, synergy, and mechanisms of splice-modifying drugs.

Drugs that target pre-mRNA splicing hold great therapeutic potential, but the quantitative understanding of how these drugs work is limited. Here we introduce mechanistically interpretable quantitative models for the sequence-specific and concentration-dependent behavior of splice-modifying drugs. Using massively parallel splicing assays, RNA-seq experiments, and precision dose-response curves, we obtain quantitative models for two small-molecule drugs, risdiplam and branaplam, developed for treating spinal muscular atrophy. The results quantitatively characterize the specificities of risdiplam and branaplam for 5' splice site sequences, suggest that branaplam recognizes 5' splice sites via two distinct interaction modes, and contradict the prevailing two-site hypothesis for risdiplam activity at SMN2 exon 7. The results also show that anomalous single-drug cooperativity, as well as multi-drug synergy, are widespread among small-molecule drugs and antisense-oligonucleotide drugs that promote exon inclusion. Our quantitative models thus clarify the mechanisms of existing treatments and provide a basis for the rational development of new therapies.

Designed active-site library reveals thousands of functional GFP variants.

Mutations in a protein active site can lead to dramatic and useful changes in protein activity. The active site, however, is sensitive to mutations due to a high density of molecular interactions, substantially reducing the likelihood of obtaining functional multipoint mutants. We introduce an atomistic and machine-learning-based approach, called high-throughput Functional Libraries (htFuncLib), that designs a sequence space in which mutations form low-energy combinations that mitigate the risk of incompatible interactions. We apply htFuncLib to the GFP chromophore-binding pocket, and, using fluorescence readout, recover >16,000 unique designs encoding as many as eight active-site mutations. Many designs exhibit substantial and useful diversity in functional thermostability (up to 96 °C), fluorescence lifetime, and quantum yield. By eliminating incompatible active-site mutations, htFuncLib generates a large diversity of functional sequences. We envision that htFuncLib will be used in one-shot optimization of activity in enzymes, binders, and other proteins.

Evolutionary paths that link orthogonal pairs of binding proteins.

Some protein binding pairs exhibit extreme specificities that functionally insulate them from homologs. Such pairs evolve mostly by accumulating single-point mutations, and mutants are selected if their affinity exceeds the threshold required for function1-4. Thus, homologous and high-specificity binding pairs bring to light an evolutionary conundrum: how does a new specificity evolve while maintaining the required affinity in each intermediate5,6? Until now, a fully functional single-mutation path that connects two orthogonal pairs has only been described where the pairs were mutationally close thus enabling experimental enumeration of all intermediates2. We present an atomistic and graph-theoretical framework for discovering low molecular strain single-mutation paths that connect two extant pairs, enabling enumeration beyond experimental capability. We apply it to two orthogonal bacterial colicin endonuclease-immunity pairs separated by 17 interface mutations7. We were not able to find a strain-free and functional path in the sequence space defined by the two extant pairs. But including mutations that bridge amino acids that cannot be exchanged through single-nucleotide mutations led us to a strain-free 19-mutation trajectory that is completely viable in vivo. Our experiments show that the specificity switch is remarkably abrupt, resulting from only one radical mutation on each partner. Furthermore, each of the critical specificity-switch mutations increases fitness, demonstrating that functional divergence could be driven by positive Darwinian selection. These results reveal how even radical functional changes in an epistatic fitness landscape may evolve.

MCRS1 modulates the heterogeneity of microtubule minus-end morphologies in mitotic spindles.

Faithful chromosome segregation requires the assembly of a bipolar spindle, consisting of two antiparallel microtubule (MT) arrays having most of their minus ends focused at the spindle poles and their plus ends overlapping in the spindle midzone. Spindle assembly, chromosome alignment, and segregation require highly dynamic MTs. The plus ends of MTs have been extensively investigated but their minus-end structure remains poorly characterized. Here, we used large-scale electron tomography to study the morphology of the MT minus ends in three dimensionally reconstructed metaphase spindles in HeLa cells. In contrast to the homogeneous open morphology of the MT plus ends at the kinetochores, we found that MT minus ends are heterogeneous, showing either open or closed morphologies. Silencing the minus end-specific stabilizer, MCRS1 increased the proportion of open MT minus ends. Altogether, these data suggest a correlation between the morphology and the dynamic state of the MT ends. Taking this heterogeneity of the MT minus-end morphologies into account, our work indicates an unsynchronized behavior of MTs at the spindle poles, thus laying the groundwork for further studies on the complexity of MT dynamics regulation.

Functional Impact and Regulation of Alternative Splicing in Mouse Heart Development and Disease.

Alternative splicing (AS) plays a major role in the generation of transcript diversity. In the heart, roles have been described for some AS variants, but the global impact and regulation of AS patterns are poorly understood. Here, we studied the AS profiles in heart disease, their relationship with heart development, and the regulatory mechanisms controlling AS dynamics in the mouse heart. We found that AS profiles characterized the different groups and that AS and gene expression changes affected independent genes and biological functions. Moreover, AS changes, specifically in heart disease, were associated with potential protein-protein interaction changes. While developmental transitions were mainly driven by the upregulation of MBNL1, AS changes in disease were driven by a complex regulatory network, where PTBP1 played a central role. Indeed, PTBP1 over-expression was sufficient to induce cardiac hypertrophy and diastolic dysfunction, potentially by perturbing AS patterns.

The SRSF4-GAS5-Glucocorticoid Receptor Axis Regulates Ventricular Hypertrophy.

[Figure: see text].

dSreg: a Bayesian model to integrate changes in splicing and RNA-binding protein activity.

MOTIVATION: Alternative splicing (AS) is an important mechanism in the generation of transcript diversity across mammals. AS patterns are dynamically regulated during development and in response to environmental changes. Defects or perturbations in its regulation may lead to cancer or neurological disorders, among other pathological conditions. The regulatory mechanisms controlling AS in a given biological context are typically inferred using a two-step framework: differential AS analysis followed by enrichment methods. These strategies require setting rather arbitrary thresholds and are prone to error propagation along the analysis. RESULTS: To overcome these limitations, we propose dSreg, a Bayesian model that integrates RNA-seq with data from regulatory features, e.g. binding sites of RNA-binding proteins. dSreg identifies the key underlying regulators controlling AS changes and quantifies their activity while simultaneously estimating the changes in exon inclusion rates. dSreg increased both the sensitivity and the specificity of the identified AS changes in simulated data, even at low read coverage. dSreg also showed improved performance when analyzing a collection of knock-down RNA-binding proteins' experiments from ENCODE, as opposed to traditional enrichment methods, such as over-representation analysis and gene set enrichment analysis. dSreg opens the possibility to integrate a large amount of readily available RNA-seq datasets at low coverage for AS analysis and allows more cost-effective RNA-seq experiments. AVAILABILITY AND IMPLEMENTATION: dSreg was implemented in python using stan and is freely available to the community at https://bitbucket.org/cmartiga/dsreg. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Severe Cardiac Dysfunction and Death Caused by Arrhythmogenic Right Ventricular Cardiomyopathy Type 5 Are Improved by Inhibition of Glycogen Synthase Kinase-3ß.

BACKGROUND: Arrhythmogenic cardiomyopathy/arrhythmogenic right ventricular cardiomyopathy (ARVC) is an inherited cardiac disease characterized by fibrofatty replacement of the myocardium, resulting in heart failure and sudden cardiac death. The most aggressive arrhythmogenic cardiomyopathy/ARVC subtype is ARVC type 5 (ARVC5), caused by a p.S358L mutation in TMEM43 (transmembrane protein 43). The function and localization of TMEM43 are unknown, as is the mechanism by which the p.S358L mutation causes the disease. Here, we report the characterization of the first transgenic mouse model of ARVC5. METHODS: We generated transgenic mice overexpressing TMEM43 in either its wild-type or p.S358L mutant (TMEM43-S358L) form in postnatal cardiomyocytes under the control of the α-myosin heavy chain promoter. RESULTS: We found that mice expressing TMEM43-S358L recapitulate the human disease and die at a young age. Mutant TMEM43 causes cardiomyocyte death and severe fibrofatty replacement. We also demonstrate that TMEM43 localizes at the nuclear membrane and interacts with emerin and ß-actin. TMEM43-S358L shows partial delocalization to the cytoplasm, reduced interaction with emerin and ß-actin, and activation of glycogen synthase kinase-3ß (GSK3ß). Furthermore, we show that targeting cardiac fibrosis has no beneficial effect, whereas overexpression of the calcineurin splice variant calcineurin Aß1 results in GSK3ß inhibition and improved cardiac function and survival. Similarly, treatment of TMEM43 mutant mice with a GSK3ß inhibitor improves cardiac function. Finally, human induced pluripotent stem cells bearing the p.S358L mutation also showed contractile dysfunction that was partially restored after GSK3ß inhibition. CONCLUSIONS: Our data provide evidence that TMEM43-S358L leads to sustained cardiomyocyte death and fibrofatty replacement. Overexpression of calcineurin Aß1 in TMEM43 mutant mice or chemical GSK3ß inhibition improves cardiac function and increases mice life span. Our results pave the way toward new therapeutic approaches for ARVC5.

Loss of SRSF3 in Cardiomyocytes Leads to Decapping of Contraction-Related mRNAs and Severe Systolic Dysfunction.

RATIONALE: RBPs (RNA binding proteins) play critical roles in the cell by regulating mRNA transport, splicing, editing, and stability. The RBP SRSF3 (serine/arginine-rich splicing factor 3) is essential for blastocyst formation and for proper liver development and function. However, its role in the heart has not been explored. OBJECTIVE: To investigate the role of SRSF3 in cardiac function. METHODS AND RESULTS: Cardiac SRSF3 expression was high at mid gestation and decreased during late embryonic development. Mice lacking SRSF3 in the embryonic heart showed impaired cardiomyocyte proliferation and died in utero. In the adult heart, SRSF3 expression was reduced after myocardial infarction, suggesting a possible role in cardiac homeostasis. To determine the role of this RBP in the adult heart, we used an inducible, cardiomyocyte-specific SRSF3 knockout mouse model. After SRSF3 depletion in cardiomyocytes, mice developed severe systolic dysfunction that resulted in death within 8 days. RNA-Seq analysis revealed downregulation of mRNAs encoding sarcomeric and calcium handling proteins. Cardiomyocyte-specific SRSF3 knockout mice also showed evidence of alternative splicing of mTOR (mammalian target of rapamycin) mRNA, generating a shorter protein isoform lacking catalytic activity. This was associated with decreased phosphorylation of 4E-BP1 (eIF4E-binding protein 1), a protein that binds to eIF4E (eukaryotic translation initiation factor 4E) and prevents mRNA decapping. Consequently, we found increased decapping of mRNAs encoding proteins involved in cardiac contraction. Decapping was partially reversed by mTOR activation. CONCLUSIONS: We show that cardiomyocyte-specific loss of SRSF3 expression results in decapping of critical mRNAs involved in cardiac contraction. The molecular mechanism underlying this effect likely involves the generation of a short mTOR isoform by alternative splicing, resulting in reduced 4E-BP1 phosphorylation. The identification of mRNA decapping as a mechanism of systolic heart failure may open the way to the development of urgently needed therapeutic tools.

De novo assembly of the zucchini genome reveals a whole-genome duplication associated with the origin of the Cucurbita genus.

The Cucurbita genus (squashes, pumpkins and gourds) includes important domesticated species such as C. pepo, C. maxima and C. moschata. In this study, we present a high-quality draft of the zucchini (C. pepo) genome. The assembly has a size of 263 Mb, a scaffold N50 of 1.8 Mb and 34 240 gene models. It includes 92% of the conserved BUSCO core gene set, and it is estimated to cover 93.0% of the genome. The genome is organized in 20 pseudomolecules that represent 81.4% of the assembly, and it is integrated with a genetic map of 7718 SNPs. Despite the small genome size, three independent lines of evidence support that the C. pepo genome is the result of a whole-genome duplication: the topology of the gene family phylogenies, the karyotype organization and the distribution of 4DTv distances. Additionally, 40 transcriptomes of 12 species of the genus were assembled and analysed together with all the other published genomes of the Cucurbitaceae family. The duplication was detected in all the Cucurbita species analysed, including C. maxima and C. moschata, but not in the more distant cucurbits belonging to the Cucumis and Citrullus genera, and it is likely to have occurred 30 ± 4 Mya in the ancestral species that gave rise to the genus.