Stepwise emergence of the neuronal gene expression program in early animal evolution.

The assembly of the neuronal and other major cell type programs occurred early in animal evolution. We can reconstruct this process by studying non-bilaterians like placozoans. These small disc-shaped animals not only have nine morphologically described cell types and no neurons but also show coordinated behaviors triggered by peptide-secreting cells. We investigated possible neuronal affinities of these peptidergic cells using phylogenetics, chromatin profiling, and comparative single-cell genomics in four placozoans. We found conserved cell type expression programs across placozoans, including populations of transdifferentiating and cycling cells, suggestive of active cell type homeostasis. We also uncovered fourteen peptidergic cell types expressing neuronal-associated components like the pre-synaptic scaffold that derive from progenitor cells with neurogenesis signatures. In contrast, earlier-branching animals like sponges and ctenophores lacked this conserved expression. Our findings indicate that key neuronal developmental and effector gene modules evolved before the advent of cnidarian/bilaterian neurons in the context of paracrine cell signaling.

Integrating physics in deep learning algorithms: a force field as a PyTorch module.

MOTIVATION: Deep learning algorithms applied to structural biology often struggle to converge to meaningful solutions when limited data is available, since they are required to learn complex physical rules from examples. State-of-the-art force-fields, however, cannot interface with deep learning algorithms due to their implementation. RESULTS: We present MadraX, a forcefield implemented as a differentiable PyTorch module, able to interact with deep learning algorithms in an end-to-end fashion. AVAILABILITY AND IMPLEMENTATION: MadraX documentation, together with tutorials and installation guide, is available at madrax.readthedocs.io.

Unique transcriptomes of sensory and non-sensory neurons: insights from Splicing Regulatory States.

Alternative Splicing (AS) programs serve as instructive signals of cell type specificity, particularly within the brain, which comprises dozens of molecularly and functionally distinct cell types. Among them, retinal photoreceptors stand out due to their unique transcriptome, making them a particularly well-suited system for studying how AS shapes cell type-specific molecular functions. Here, we use the Splicing Regulatory State (SRS) as a novel framework to discuss the splicing factors governing the unique AS pattern of photoreceptors, and how this pattern may aid in the specification of their highly specialized sensory cilia. In addition, we discuss how other sensory cells with ciliated structures, for which data is much scarcer, also rely on specific SRSs to implement a proteome specialized in the detection of sensory stimuli. By reviewing the general rules of cell type- and tissue-specific AS programs, firstly in the brain and subsequently in specialized sensory neurons, we propose a novel paradigm on how SRSs are established and how they can diversify. Finally, we illustrate how SRSs shape the outcome of mutations in splicing factors to produce cell type-specific phenotypes that can lead to various human diseases.

Engineering signal transduction pathways.

Cells respond to their environment by sensing signals and translating them into changes in gene expression. In recent years, synthetic networks have been designed in both prokaryotic and eukaryotic systems to create new functionalities and for specific applications. In this review, we discuss the challenges associated with engineering signal transduction pathways. Furthermore, we address advantages and disadvantages of engineering signaling pathways in prokaryotic and eukaryotic cells, highlighting recent examples, and discuss how progress in synthetic biology might impact biotechnology and biomedicine.

Single-read tRNA-seq analysis reveals coordination of tRNA modification and aminoacylation and fragmentation.

Transfer RNA (tRNA) utilizes multiple properties of abundance, modification, and aminoacylation in translational regulation. These properties were typically studied one-by-one; however, recent advance in high throughput tRNA sequencing enables their simultaneous assessment in the same sequencing data. How these properties are coordinated at the transcriptome level is an open question. Here, we develop a single-read tRNA analysis pipeline that takes advantage of the pseudo single-molecule nature of tRNA sequencing in NGS libraries. tRNAs are short enough that a single NGS read can represent one tRNA molecule, and can simultaneously report on the status of multiple modifications, aminoacylation, and fragmentation of each molecule. We find correlations among modification-modification, modification-aminoacylation and modification-fragmentation. We identify interdependencies among one of the most common tRNA modifications, m1A58, as coordinators of tissue-specific gene expression. Our method, SingLe-read Analysis of Crosstalks (SLAC), reveals tRNAome-wide networks of modifications, aminoacylation, and fragmentation. We observe changes of these networks under different stresses, and assign a function for tRNA modification in translational regulation and fragment biogenesis. SLAC leverages the richness of the tRNA-seq data and provides new insights on the coordination of tRNA properties.

An ancient founder mutation located between ROBO1 and ROBO2 is responsible for increased microtia risk in Amerindigenous populations.

Microtia is a congenital malformation that encompasses mild hypoplasia to complete loss of the external ear, or pinna. Although the contribution of genetic variation and environmental factors to microtia remains elusive, Amerindigenous populations have the highest reported incidence. Here, using both transmission disequilibrium tests and association studies in microtia trios (parents and affected child) and microtia cohorts enrolled in Latin America, we map an ∼10-kb microtia locus (odds ratio = 4.7; P = 6.78e-18) to the intergenic region between Roundabout 1 (ROBO1) and Roundabout 2 (ROBO2) (chr3: 78546526 to 78555137). While alleles at the microtia locus significantly increase the risk of microtia, their penetrance is low (<1%). We demonstrate that the microtia locus contains a polymorphic complex repeat element that is expanded in affected individuals. The locus is located near a chromatin loop region that regulates ROBO1 and ROBO2 expression in induced pluripotent stem cell­derived neural crest cells. Furthermore, we use single nuclear RNA sequencing to demonstrate ROBO1 and ROBO2 expression in both fibroblasts and chondrocytes of the mature human pinna. Because the microtia allele is enriched in Amerindigenous populations and is shared by some East Asian subjects with craniofacial malformations, we propose that both populations share a mutation that arose in a common ancestor prior to the ancient migration of Eurasian populations into the Americas and that the high incidence of microtia among Amerindigenous populations reflects the population bottleneck that occurred during the migration out of Eurasia.

Specialization of the photoreceptor transcriptome by Srrm3-dependent microexons is required for outer segment maintenance and vision.

Retinal photoreceptors have a distinct transcriptomic profile compared to other neuronal subtypes, likely reflecting their unique cellular morphology and function in the detection of light stimuli by way of the ciliary outer segment. We discovered a layer of this molecular specialization by revealing that the vertebrate retina expresses the largest number of tissue-enriched microexons of all tissue types. A subset of these microexons is included exclusively in photoreceptor transcripts, particularly in genes involved in cilia biogenesis and vesicle-mediated transport. This microexon program is regulated by Srrm3, a paralog of the neural microexon regulator Srrm4. Despite the fact that both proteins positively regulate retina microexons in vitro, only Srrm3 is highly expressed in mature photoreceptors. Its deletion in zebrafish results in widespread down-regulation of microexon inclusion from early developmental stages, followed by other transcriptomic alterations, severe photoreceptor defects, and blindness. These results shed light on the transcriptomic specialization and functionality of photoreceptors, uncovering unique cell type-specific roles for Srrm3 and microexons with implications for retinal diseases.

Comprehensive quantitative modeling of translation efficiency in a genome-reduced bacterium.

Translation efficiency has been mainly studied by ribosome profiling, which only provides an incomplete picture of translation kinetics. Here, we integrated the absolute quantifications of tRNAs, mRNAs, RNA half-lives, proteins, and protein half-lives with ribosome densities and derived the initiation and elongation rates for 475 genes (67% of all genes), 73 with high precision, in the bacterium Mycoplasma pneumoniae (Mpn). We found that, although the initiation rate varied over 160-fold among genes, most of the known factors had little impact on translation efficiency. Local codon elongation rates could not be fully explained by the adaptation to tRNA abundances, which varied over 100-fold among tRNA isoacceptors. We provide a comprehensive quantitative view of translation efficiency, which suggests the existence of unidentified mechanisms of translational regulation in Mpn.

Bacterial expression of a designed single-chain IL-10 prevents severe lung inflammation.

Interleukin-10 (IL-10) is an anti-inflammatory cytokine that is active as a swapped domain dimer and is used in bacterial therapy of gut inflammation. IL-10 can be used as treatment of a wide range of pulmonary diseases. Here we have developed a non-pathogenic chassis (CV8) of the human lung bacterium Mycoplasma pneumoniae (MPN) to treat lung diseases. We find that IL-10 expression by MPN has a limited impact on the lung inflammatory response in mice. To solve these issues, we rationally designed a single-chain IL-10 (SC-IL10) with or without surface mutations, using our protein design software (ModelX and FoldX). As compared to the IL-10 WT, the designed SC-IL10 molecules increase the effective expression in MPN four-fold, and the activity in mouse and human cell lines between 10 and 60 times, depending on the cell line. The SC-IL10 molecules expressed in the mouse lung by CV8 in vivo have a powerful anti-inflammatory effect on Pseudomonas aeruginosa lung infection. This rational design strategy could be used to other molecules with immunomodulatory properties used in bacterial therapy.

Silencing of SRRM4 suppresses microexon inclusion and promotes tumor growth across cancers.

RNA splicing is widely dysregulated in cancer, frequently due to altered expression or activity of splicing factors (SFs). Microexons are extremely small exons (3-27 nucleotides long) that are highly evolutionarily conserved and play critical roles in promoting neuronal differentiation and development. Inclusion of microexons in mRNA transcripts is mediated by the SF Serine/Arginine Repetitive Matrix 4 (SRRM4), whose expression is largely restricted to neural tissues. However, microexons have been largely overlooked in prior analyses of splicing in cancer, as their small size necessitates specialized computational approaches for their detection. Here, we demonstrate that despite having low expression in normal nonneural tissues, SRRM4 is further silenced in tumors, resulting in the suppression of normal microexon inclusion. Remarkably, SRRM4 is the most consistently silenced SF across all tumor types analyzed, implying a general advantage of microexon down-regulation in cancer independent of its tissue of origin. We show that this silencing is favorable for tumor growth, as decreased SRRM4 expression in tumors is correlated with an increase in mitotic gene expression, and up-regulation of SRRM4 in cancer cell lines dose-dependently inhibits proliferation in vitro and in a mouse xenograft model. Further, this proliferation inhibition is accompanied by induction of neural-like expression and splicing patterns in cancer cells, suggesting that SRRM4 expression shifts the cell state away from proliferation and toward differentiation. We therefore conclude that SRRM4 acts as a proliferation brake, and tumors gain a selective advantage by cutting off this brake.

Developing a platform for secretion of biomolecules in Mycoplasma feriruminatoris.

BACKGROUND: Having a simple and fast dividing organism capable of producing and exposing at its surface or secreting functional complex biomolecules with disulphide bridges is of great interest. The mycoplasma bacterial genus offers a set of relevant properties that make it an interesting chassis for such purposes, the main one being the absence of a cell wall. However, due to their slow growth, they have rarely been considered as a potential platform in this respect. This notion may be challenged with the recent discovery of Mycoplasma feriruminatoris, a species with a dividing time close to that of common microbial workhorses. So far, no tools for heterologous protein expression nor secretion have been described for it. RESULTS: The work presented here develops the fast-dividing M. feriruminatoris as a tool for secreting functional biomolecules of therapeutic interest that could be used for screening functional mutants as well as potentially for protein-protein interactions. Based on RNAseq, quantitative proteomics and promoter sequence comparison we have rationally designed optimal promoter sequences. Then, using in silico analysis, we have identified putative secretion signals that we validated using a luminescent reporter. The potential of the resulting secretion cassette has been shown with set of active clinically relevant proteins (interleukins and nanobodies). CONCLUSIONS: We have engineered Mycoplasma feriruminatoris for producing and secreting functional proteins of medical interest.

SURE editing: combining oligo-recombineering and programmable insertion/deletion of selection markers to efficiently edit the Mycoplasma pneumoniae genome.

The development of advanced genetic tools is boosting microbial engineering which can potentially tackle wide-ranging challenges currently faced by our society. Here we present SURE editing, a multi-recombinase engineering rationale combining oligonucleotide recombineering with the selective capacity of antibiotic resistance via transient insertion of selector plasmids. We test this method in Mycoplasma pneumoniae, a bacterium with a very inefficient native recombination machinery. Using SURE editing, we can seamlessly generate, in a single step, a wide variety of genome modifications at high efficiencies, including the largest possible deletion of this genome (30 Kb) and the targeted complementation of essential genes in the deletion of a region of interest. Additional steps can be taken to remove the selector plasmid from the edited area, to obtain markerless or even scarless edits. Of note, SURE editing is compatible with different site-specific recombinases for mediating transient plasmid integration. This battery of selector plasmids can be used to select different edits, regardless of the target sequence, which significantly reduces the cloning load associated to genome engineering projects. Given the proven functionality in several microorganisms of the machinery behind the SURE editing logic, this method is likely to represent a valuable advance for the synthetic biology field.

Geniposide and Harpagoside Functionalized Cerium Oxide Nanoparticles as a Potential Neuroprotective.

Alzheimer's disease is associated with protein aggregation, oxidative stress, and the role of acetylcholinesterase in the pathology of the disease. Previous investigations have demonstrated that geniposide and harpagoside protect the brain neurons, and cerium nanoparticles (CeO2 NPs) have potent redox and antioxidant properties. Thus, the effect of nanoparticles of Ce NPs and geniposide and harpagoside (GH/CeO2 NPs) on ameliorating AD pathogenesis was established on AlCl3-induced AD in mice and an aggregation proteins test in vitro. Findings of spectroscopy analysis have revealed that GH/CeO2 NPs are highly stable, nano-size, spherical in shape, amorphous nature, and a total encapsulation of GH in cerium. Treatments with CeO2 NPs, GH/CeO2 NPs, and donepezil used as positive control inhibit fibril formation and protein aggregation, protect structural modifications in the BSA-ribose system, have the ability to counteract Tau protein aggregation and amyloid-ß1-42 aggregation under fibrillation condition, and are able to inhibit AChE and BuChE. While the GH/CeO2 NPs, treatment in AD induced by AlCl3 inhibited amyloid-ß1-42, substantially enhanced the memory, the cognition coordination of movement in part AD pathogenesis may be alleviated through reducing amyloidogenic pathway and AChE and BuChE activities. The findings of this work provide important comprehension of the chemoprotective activities of iridoids combined with nanoparticles. This could be useful in the development of new therapeutic methods for the treatment of neurodegenerative diseases.

pyFoldX: enabling biomolecular analysis and engineering along structural ensembles.

SUMMARY: Recent years have seen an increase in the number of structures available, not only for new proteins but also for the same protein crystallized with different molecules and proteins. While protein design software has proven to be successful in designing and modifying proteins, they can also be overly sensitive to small conformational differences between structures of the same protein. To cope with this, we introduce here pyFoldX, a python library that allows the integrative analysis of structures of the same protein using FoldX, an established forcefield and modelling software. The library offers new functionalities for handling different structures of the same protein, an improved molecular parametrization module and an easy integration with the data analysis ecosystem of the python programming language. AVAILABILITY AND IMPLEMENTATION: pyFoldX rely on the FoldX software for energy calculations and modelling, which can be downloaded upon registration in http://foldxsuite.crg.eu/ and its licence is free of charge for academics. The pyFoldX library is open-source. Full details on installation, tutorials covering the library functionality and the scripts used to generate the data and figures presented in this paper are available at https://github.com/leandroradusky/pyFoldX. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Damaging variants in FOXI3 cause microtia and craniofacial microsomia.

PURPOSE: Craniofacial microsomia (CFM) represents a spectrum of craniofacial malformations, ranging from isolated microtia with or without aural atresia to underdevelopment of the mandible, maxilla, orbit, facial soft tissue, and/or facial nerve. The genetic causes of CFM remain largely unknown. METHODS: We performed genome sequencing and linkage analysis in patients and families with microtia and CFM of unknown genetic etiology. The functional consequences of damaging missense variants were evaluated through expression of wild-type and mutant proteins in vitro. RESULTS: We studied a 5-generation kindred with microtia, identifying a missense variant in FOXI3 (p.Arg236Trp) as the cause of disease (logarithm of the odds = 3.33). We subsequently identified 6 individuals from 3 additional kindreds with microtia-CFM spectrum phenotypes harboring damaging variants in FOXI3, a regulator of ectodermal and neural crest development. Missense variants in the nuclear localization sequence were identified in cases with isolated microtia with aural atresia and found to affect subcellular localization of FOXI3. Loss of function variants were found in patients with microtia and mandibular hypoplasia (CFM), suggesting dosage sensitivity of FOXI3. CONCLUSION: Damaging variants in FOXI3 are the second most frequent genetic cause of CFM, causing 1% of all cases, including 13% of familial cases in our cohort.

Therapeutic potential of fetal liver cell transplantation in hemophilia A mice.

Hemophilia A (HA) cell therapy approaches in pediatric individuals require suitable factor (F)VIII-producing cells for stable engraftment. Liver sinusoidal endothelial cells (LSEC) and hematopoietic stem cells (HSC) have been demonstrated to be suitable for the treatment of adult HA mice. However, after transplantation in busulfan (BU)-conditioned newborn mice, adult LSEC/HSC cannot efficiently engraft, while murine fetal liver (FL) hemato/vascular cells from embryonic day 11-13 of gestation (E11-E13), strongly engraft the hematopoietic and endothelial compartments while also secreting FVIII. Our aim was to investigate the engraftment of FL cells in newborn HA mice to obtain a suitable "proof of concept" for the development of a new HA treatment in neonates. Hence, we transplanted FL E11 or E13 cells and adult bone marrow (BM) cells into newborn HA mice with or without BU preconditioning. Engraftment levels and FVIII activity were assessed starting from 6 weeks after transplantation. FL E11-E13+ BU transplanted newborns reached up to 95% engraftment with stable FVIII activity levels observed for 16 months. FL E13 cells showed engraftment ability even in the absence of BU preconditioning, while FL E11 cells did not. BM BU transplanted newborn HA mice showed high levels of engraftment; nevertheless, in contrast to FL cells, BM cells cannot engraft HA newborns in BU non-conditioning regimen. Finally, none of the transplanted mice developed anti-FVIII antibodies. Overall, this study sheds some light on the therapeutic potential of healthy FL cells in the cure of HA neonatal/pediatric patients.

Mutation bias within oncogene families is related to proliferation-specific codon usage.

It is well known that in cancer gene families some members are more frequently mutated in tumor samples than their family counterparts. A paradigmatic case of this phenomenon is KRAS from the RAS family. Different explanations have been proposed ranging from differential interaction with other proteins to preferential expression or localization. Interestingly, it has been described that despite the high amino acid identity between RAS family members, KRAS employs an intriguing differential codon usage. Here, we found that this phenomenon is not exclusive to the RAS family. Indeed, in the RAS family and other oncogene families with two or three members, the most prevalently mutated gene in tumor samples employs a differential codon usage that is characteristic of genes involved in proliferation. Prompted by these observations, we chose the RAS family to experimentally demonstrate that the translation efficiency of oncogenes that are preferentially mutated in tumor samples is increased in proliferative cells compared to quiescent cells. These results were further validated by assessing the translation efficiency of KRAS in cell lines that differ in their tRNA expression profile. These differences are related to the cell division rate of the studied cells and thus suggest an important role in context-specific oncogene expression regulation. Altogether, our study demonstrates that dynamic translation programs contribute to shaping the expression profiles of oncogenes. Therefore, we propose this codon bias as a regulatory layer to control cell context-specific expression and explain the differential prevalence of mutations in certain members of oncogene families.

Role of the CB2 Cannabinoid Receptor in the Regulation of Food Intake: A Systematic Review.

The CB2 cannabinoid receptor has been found in brain areas that are part of the reward system and has been shown to play a role in food intake regulation. Herein, we conducted a systematic review of studies assessing the role of the CB2 receptor in food intake regulation. Records from the PubMed, Scopus, and EBSCO databases were screened, resulting in 13 studies that were used in the present systematic review, following the PRISMA guidelines. A risk of bias assessment was carried out using the tool of the Systematic Review Center for Laboratory Animal Experimentation (SYRCLE). The studies analyzed used two main strategies: (1) the intraperitoneal or intracerebroventricular administration of a CB2 agonist/antagonist; and (2) depletion of CB2 receptors via knockout in mice. Both strategies are useful in identifying the role of the CB2 receptor in food intake in standard and palatable diets. The conclusions derived from animal models showed that CB2 receptors are necessary for modulating food intake and mediating energy balance.

robustica: customizable robust independent component analysis.

BACKGROUND: Independent Component Analysis (ICA) allows the dissection of omic datasets into modules that help to interpret global molecular signatures. The inherent randomness of this algorithm can be overcome by clustering many iterations of ICA together to obtain robust components. Existing algorithms for robust ICA are dependent on the choice of clustering method and on computing a potentially biased and large Pearson distance matrix. RESULTS: We present robustica, a Python-based package to compute robust independent components with a fully customizable clustering algorithm and distance metric. Here, we exploited its customizability to revisit and optimize robust ICA systematically. Of the 6 popular clustering algorithms considered, DBSCAN performed the best at clustering independent components across ICA iterations. To enable using Euclidean distances, we created a subroutine that infers and corrects the components' signs across ICA iterations. Our subroutine increased the resolution, robustness, and computational efficiency of the algorithm. Finally, we show the applicability of robustica by dissecting over 500 tumor samples from low-grade glioma (LGG) patients, where we define two new gene expression modules with key modulators of tumor progression upon IDH1 and TP53 mutagenesis. CONCLUSION: robustica brings precise, efficient, and customizable robust ICA into the Python toolbox. Through its customizability, we explored how different clustering algorithms and distance metrics can further optimize robust ICA. Then, we showcased how robustica can be used to discover gene modules associated with combinations of features of biological interest. Taken together, given the broad applicability of ICA for omic data analysis, we envision robustica will facilitate the seamless computation and integration of robust independent components in large pipelines.

Protein haploinsufficiency drivers identify MYBPC3 variants that cause hypertrophic cardiomyopathy.

Hypertrophic cardiomyopathy (HCM) is the most common inherited cardiac disease. Variants in MYBPC3, the gene encoding cardiac myosin-binding protein C (cMyBP-C), are the leading cause of HCM. However, the pathogenicity status of hundreds of MYBPC3 variants found in patients remains unknown, as a consequence of our incomplete understanding of the pathomechanisms triggered by HCM-causing variants. Here, we examined 44 nontruncating MYBPC3 variants that we classified as HCM-linked or nonpathogenic according to cosegregation and population genetics criteria. We found that around half of the HCM-linked variants showed alterations in RNA splicing or protein stability, both of which can lead to cMyBP-C haploinsufficiency. These protein haploinsufficiency drivers associated with HCM pathogenicity with 100% and 94% specificity, respectively. Furthermore, we uncovered that 11% of nontruncating MYBPC3 variants currently classified as of uncertain significance in ClinVar induced one of these molecular phenotypes. Our strategy, which can be applied to other conditions induced by protein loss of function, supports the idea that cMyBP-C haploinsufficiency is a fundamental pathomechanism in HCM.