PHEIGES: all-cell-free phage synthesis and selection from engineered genomes

Bacteriophages constitute an invaluable biological reservoir for biotechnology and medicine. The ability to exploit such vast resources is hampered by the lack of methods to rapidly engineer, assemble, package genomes, and select phages. Cell-free transcription-translation (TXTL) offers experimental settings to address such a limitation. Here, we describe PHage Engineering by In vitro Gene Expression and Selection (PHEIGES) using T7 phage genome and Escherichia coli TXTL. Phage genomes are assembled in vitro from PCR-amplified fragments and directly expressed in batch TXTL reactions to produce up to 1011 PFU/ml engineered phages within one day. We further demonstrate a significant genotype-phenotype linkage of phage assembly in bulk TXTL. This enables rapid selection of phages with altered rough lipopolysaccharides specificity from phage genomes incorporating tail fiber mutant libraries. We establish the scalability of PHEIGES by one pot assembly of such mutants with fluorescent gene integration and 10% length-reduced genome.

The mesencephalic locomotor region recruits V2a reticulospinal neurons to drive forward locomotion in larval zebrafish

The mesencephalic locomotor region (MLR) is a brain stem area whose stimulation triggers graded forward locomotion. How MLR neurons recruit downstream vsx2+ (V2a) reticulospinal neurons (RSNs) is poorly understood. Here, to overcome this challenge, we uncovered the locus of MLR in transparent larval zebrafish and show that the MLR locus is distinct from the nucleus of the medial longitudinal fasciculus. MLR stimulations reliably elicit forward locomotion of controlled duration and frequency. MLR neurons recruit V2a RSNs via projections onto somata in pontine and retropontine areas, and onto dendrites in the medulla. High-speed volumetric imaging of neuronal activity reveals that strongly MLR-coupled RSNs are active for steering or forward swimming, whereas weakly MLR-coupled medullary RSNs encode the duration and frequency of the forward component. Our study demonstrates how MLR neurons recruit specific V2a RSNs to control the kinematics of forward locomotion and suggests conservation of the motor functions of V2a RSNs across vertebrates.

Multi-feature clustering of CTCF binding creates robustness for loop extrusion blocking and Topologically Associating Domain boundaries

Topologically Associating Domains (TADs) separate vertebrate genomes into insulated regulatory neighborhoods that focus genome-associated processes. TADs are formed by Cohesin-mediated loop extrusion, with many TAD boundaries consisting of clustered binding sites of the CTCF insulator protein. Here we determine how this clustering of CTCF binding contributes to the blocking of loop extrusion and the insulation between TADs. We identify enrichment of three features of CTCF binding at strong TAD boundaries, consisting of strongly bound and closely spaced CTCF binding peaks, with a further enrichment of DNA-binding motifs within these peaks. Using multi-contact Nano-C analysis in cells with normal and perturbed CTCF binding, we establish that individual CTCF binding sites contribute to the blocking of loop extrusion, but in an incomplete manner. When clustered, individual CTCF binding sites thus create a stepwise insulation between neighboring TADs. Based on these results, we propose a model whereby multiple instances of temporal loop extrusion blocking create strong insulation between TADs.

Dimensionality and Ramping: Signatures of Sentence Integration in the Dynamics of Brains and Deep Language Models

A sentence is more than the sum of its words: its meaning depends on how they combine with one another. The brain mechanisms underlying such semantic composition remain poorly understood. To shed light on the neural vector code underlying semantic composition, we introduce two hypotheses: (1) the intrinsic dimensionality of the space of neural representations should increase as a sentence unfolds, paralleling the growing complexity of its semantic representation; and (2) this progressive integration should be reflected in ramping and sentence-final signals. To test these predictions, we designed a dataset of closely matched normal and jabberwocky sentences (composed of meaningless pseudo words) and displayed them to deep language models and to 11 human participants (5 men and 6 women) monitored with simultaneous MEG and intracranial EEG. In both deep language models and electrophysiological data, we found that representational dimensionality was higher for meaningful sentences than jabberwocky. Furthermore, multivariate decoding of normal versus jabberwocky confirmed three dynamic patterns: (1) a phasic pattern following each word, peaking in temporal and parietal areas; (2) a ramping pattern, characteristic of bilateral inferior and middle frontal gyri; and (3) a sentence-final pattern in left superior frontal gyrus and right orbitofrontal cortex. These results provide a first glimpse into the neural geometry of semantic integration and constrain the search for a neural code of linguistic composition.SIGNIFICANCE STATEMENT Starting from general linguistic concepts, we make two sets of predictions in neural signals evoked by reading multiword sentences. First, the intrinsic dimensionality of the representation should grow with additional meaningful words. Second, the neural dynamics should exhibit signatures of encoding, maintaining, and resolving semantic composition. We successfully validated these hypotheses in deep neural language models, artificial neural networks trained on text and performing very well on many natural language processing tasks. Then, using a unique combination of MEG and intracranial electrodes, we recorded high-resolution brain data from human participants while they read a controlled set of sentences. Time-resolved dimensionality analysis showed increasing dimensionality with meaning, and multivariate decoding allowed us to isolate the three dynamical patterns we had hypothesized.

Phenotypic effects of genetic variants associated with autism

While over 100 genes have been associated with autism, little is known about the prevalence of variants affecting them in individuals without a diagnosis of autism. Nor do we fully appreciate the phenotypic diversity beyond the formal autism diagnosis. Based on data from more than 13,000 individuals with autism and 210,000 undiagnosed individuals, we estimated the odds ratios for autism associated to rare loss-of-function (LoF) variants in 185 genes associated with autism, alongside 2,492 genes displaying intolerance to LoF variants. In contrast to autism-centric approaches, we investigated the correlates of these variants in individuals without a diagnosis of autism. We show that these variants are associated with a small but significant decrease in fluid intelligence, qualification level and income and an increase in metrics related to material deprivation. These effects were larger for autism-associated genes than in other LoF-intolerant genes. Using brain imaging data from 21,040 individuals from the UK Biobank, we could not detect significant differences in the overall brain anatomy between LoF carriers and non-carriers. Our results highlight the importance of studying the effect of the genetic variants beyond categorical diagnosis and the need for more research to understand the association between these variants and sociodemographic factors, to best support individuals carrying these variants.

Ultrafast light targeting for high-throughput precise control of neuronal networks

Two-photon, single-cell resolution optogenetics based on holographic light-targeting approaches enables the generation of precise spatiotemporal neuronal activity patterns and thus a broad range of experimental applications, such as high throughput connectivity mapping and probing neural codes for perception. Yet, current holographic approaches limit the resolution for tuning the relative spiking time of distinct cells to a few milliseconds, and the achievable number of targets to 100-200, depending on the working depth. To overcome these limitations and expand the capabilities of single-cell optogenetics, we introduce an ultra-fast sequential light targeting (FLiT) optical configuration based on the rapid switching of a temporally focused beam between holograms at kHz rates. We used FLiT to demonstrate two illumination protocols, termed hybrid- and cyclic-illumination, and achieve sub-millisecond control of sequential neuronal activation and high throughput multicell illumination in vitro (mouse organotypic and acute brain slices) and in vivo (zebrafish larvae and mice), while minimizing light-induced thermal rise. These approaches will be important for experiments that require rapid and precise cell stimulation with defined spatio-temporal activity patterns and optical control of large neuronal ensembles.

Structural basis and dynamics of Chikungunya alphavirus RNA capping by nsP1 capping pores

Alphaviruses are emerging positive-stranded RNA viruses which replicate and transcribe their genomes in membranous organelles formed in the cell cytoplasm. The nonstructural protein 1 (nsP1) is responsible for viral RNA capping and gates the replication organelles by assembling into monotopic membrane-associated dodecameric pores. The capping pathway is unique to Alphaviruses; beginning with the N⁷ methylation of a guanosine triphosphate (GTP) molecule, followed by the covalent linkage of an m⁷GMP group to a conserved histidine in nsP1 and the transfer of this cap structure to a diphosphate RNA. Here, we provide structural snapshots of different stages of the reaction pathway showing how nsP1 pores recognize the substrates of the methyl-transfer reaction, GTP and S-adenosyl methionine (SAM), how the enzyme reaches a metastable postmethylation state with SAH and m⁷GTP in the active site, and the subsequent covalent transfer of m⁷GMP to nsP1 triggered by the presence of RNA and postdecapping reaction conformational changes inducing the opening of the pore. In addition, we biochemically characterize the capping reaction, demonstrating specificity for the RNA substrate and the reversibility of the cap transfer resulting in decapping activity and the release of reaction intermediates. Our data identify the molecular determinants allowing each pathway transition, providing an explanation for the need for the SAM methyl donor all along the pathway and clues about the conformational rearrangements associated to the enzymatic activity of nsP1. Together, our results set ground for the structural and functional understanding of alphavirus RNA-capping and the design of antivirals.

Size- and position-dependent cytoplasm viscoelasticity through hydrodynamic interactions with the cell surface

Many studies of cytoplasm rheology have focused on small components in the submicrometer scale. However, the cytoplasm also baths large organelles like nuclei, microtubule asters, or spindles that often take significant portions of cells and move across the cytoplasm to regulate cell division or polarization. Here, we translated passive components of sizes ranging from few up to ~50 percents of the cell diameter, through the vast cytoplasm of live sea urchin eggs, with calibrated magnetic forces. Creep and relaxation responses indicate that for objects larger than the micron size, the cytoplasm behaves as a Jeffreys material, viscoelastic at short timescales, and fluidizing at longer times. However, as component size approached that of cells, cytoplasm viscoelastic resistance increased in a nonmonotonic manner. Flow analysis and simulations suggest that this size-dependent viscoelasticity emerges from hydrodynamic interactions between the moving object and the static cell surface. This effect also yields to position-dependent viscoelasticity with objects initially closer to the cell surface being harder to displace. These findings suggest that the cytoplasm hydrodynamically couples large organelles to the cell surface to restrain their motion, with important implications for cell shape sensing and cellular organization.

Base-editing-mediated dissection of a γ-globin cis-regulatory element for the therapeutic reactivation of fetal hemoglobin expression

Sickle cell disease and β-thalassemia affect the production of the adult β-hemoglobin chain. The clinical severity is lessened by mutations that cause fetal γ-globin expression in adult life (i.e., the hereditary persistence of fetal hemoglobin). Mutations clustering ~200 nucleotides upstream of the HBG transcriptional start sites either reduce binding of the LRF repressor or recruit the KLF1 activator. Here, we use base editing to generate a variety of mutations in the -200 region of the HBG promoters, including potent combinations of four to eight γ-globin-inducing mutations. Editing of patient hematopoietic stem/progenitor cells is safe, leads to fetal hemoglobin reactivation and rescues the pathological phenotype. Creation of a KLF1 activator binding site is the most potent strategy - even in long-term repopulating hematopoietic stem/progenitor cells. Compared with a Cas9-nuclease approach, base editing avoids the generation of insertions, deletions and large genomic rearrangements and results in higher γ-globin levels. Our results demonstrate that base editing of HBG promoters is a safe, universal strategy for treating β-hemoglobinopathies.

Systematic and quantitative view of the antiviral arsenal of prokaryotes

Bacteria and archaea have developed multiple antiviral mechanisms, and genomic evidence indicates that several of these antiviral systems co-occur in the same strain. Here, we introduce DefenseFinder, a tool that automatically detects known antiviral systems in prokaryotic genomes. We use DefenseFinder to analyse 21000 fully sequenced prokaryotic genomes, and find that antiviral strategies vary drastically between phyla, species and strains. Variations in composition of antiviral systems correlate with genome size, viral threat, and lifestyle traits. DefenseFinder will facilitate large-scale genomic analysis of antiviral defense systems and the study of host-virus interactions in prokaryotes.

Environmental cues from neural crest derivatives act as metastatic triggers in an embryonic neuroblastoma model

Embryonic malignant transformation is concomitant to organogenesis, often affecting multipotent and migratory progenitors. While lineage relationships between malignant cells and their physiological counterparts are extensively investigated, the contribution of exogenous embryonic signals is not fully known. Neuroblastoma (NB) is a childhood malignancy of the peripheral nervous system arising from the embryonic trunk neural crest (NC) and characterized by heterogeneous and interconvertible tumor cell identities. Here, using experimental models mimicking the embryonic context coupled to proteomic and transcriptomic analyses, we show that signals released by embryonic sympathetic ganglia, including Olfactomedin-1, induce NB cells to shift from a noradrenergic to mesenchymal identity, and to activate a gene program promoting NB metastatic onset and dissemination. From this gene program, we extract a core signature specifically shared by metastatic cancers with NC origin. This reveals non-cell autonomous embryonic contributions regulating the plasticity of NB identities and setting pro-dissemination gene programs common to NC-derived cancers.

Higher order genetic interactions switch cancer genes from two-hit to one-hit drivers

The classic two-hit model posits that both alleles of a tumor suppressor gene (TSG) must be inactivated to cause cancer. In contrast, for some oncogenes and haploinsufficient TSGs, a single genetic alteration can suffice to increase tumor fitness. Here, by quantifying the interactions between mutations and copy number alterations (CNAs) across 10,000 tumors, we show that many cancer genes actually switch between acting as one-hit or two-hit drivers. Third order genetic interactions identify the causes of some of these switches in dominance and dosage sensitivity as mutations in other genes in the same biological pathway. The correct genetic model for a gene thus depends on the other mutations in a genome, with a second hit in the same gene or an alteration in a different gene in the same pathway sometimes representing alternative evolutionary paths to cancer.

Myofibril and mitochondria morphogenesis are coordinated by a mechanical feedback mechanism in muscle

Complex animals build specialised muscles to match specific biomechanical and energetic needs. Hence, composition and architecture of sarcomeres and mitochondria are muscle type specific. However, mechanisms coordinating mitochondria with sarcomere morphogenesis are elusive. Here we use Drosophila muscles to demonstrate that myofibril and mitochondria morphogenesis are intimately linked. In flight muscles, the muscle selector spalt instructs mitochondria to intercalate between myofibrils, which in turn mechanically constrain mitochondria into elongated shapes. Conversely in cross-striated leg muscles, mitochondria networks surround myofibril bundles, contacting myofibrils only with thin extensions. To investigate the mechanism causing these differences, we manipulated mitochondrial dynamics and found that increased mitochondrial fusion during myofibril assembly prevents mitochondrial intercalation in flight muscles. Strikingly, this causes the expression of cross-striated muscle specific sarcomeric proteins. Consequently, flight muscle myofibrils convert towards a partially cross-striated architecture. Together, these data suggest a biomechanical feedback mechanism downstream of spalt synchronizing mitochondria with myofibril morphogenesis.

The desensitization pathway of GABAA receptors, one subunit at a time

GABAA receptors mediate most inhibitory synaptic transmission in the brain of vertebrates. Following GABA binding and fast activation, these receptors undergo a slower desensitization, the conformational pathway of which remains largely elusive. To explore the mechanism of desensitization, we used concatemeric α1β2γ2 GABAA receptors to selectively introduce gain-of-desensitization mutations one subunit at a time. A library of twenty-six mutant combinations was generated and their bi-exponential macroscopic desensitization rates measured. Introducing mutations at the different subunits shows a strongly asymmetric pattern with a key contribution of the γ2 subunit, and combining mutations results in marked synergistic effects indicating a non-concerted mechanism. Kinetic modelling indeed suggests a pathway where subunits move independently, the desensitization of two subunits being required to occlude the pore. Our work thus hints towards a very diverse and labile conformational landscape during desensitization, with potential implications in physiology and pharmacology.