SPLASH: A statistical, reference-free genomic algorithm unifies biological discovery.

Today's genomics workflows typically require alignment to a reference sequence, which limits discovery. We introduce a unifying paradigm, SPLASH (Statistically Primary aLignment Agnostic Sequence Homing), which directly analyzes raw sequencing data, using a statistical test to detect a signature of regulation: sample-specific sequence variation. SPLASH detects many types of variation and can be efficiently run at scale. We show that SPLASH identifies complex mutation patterns in SARS-CoV-2, discovers regulated RNA isoforms at the single-cell level, detects the vast sequence diversity of adaptive immune receptors, and uncovers biology in non-model organisms undocumented in their reference genomes: geographic and seasonal variation and diatom association in eelgrass, an oceanic plant impacted by climate change, and tissue-specific transcripts in octopus. SPLASH is a unifying approach to genomic analysis that enables expansive discovery without metadata or references.

A Well-Mixed E. coli Genome: Widespread Contacts Revealed by Tracking Mu Transposition.

The three-dimensional structures of chromosomes are increasingly being recognized as playing a major role in cellular regulatory states. The efficiency and promiscuity of phage Mu transposition was exploited to directly measure in vivo interactions between genomic loci in E. coli. Two global organizing principles have emerged: first, the chromosome is well-mixed and uncompartmentalized, with transpositions occurring freely between all measured loci; second, several gene families/regions show "clustering": strong three-dimensional co-localization regardless of linear genomic distance. The activities of the SMC/condensin protein MukB and nucleoid-compacting protein subunit HU-α are essential for the well-mixed state; HU-α is also needed for clustering of 6/7 ribosomal RNA-encoding loci. The data are explained by a model in which the chromosomal structure is driven by dynamic competition between DNA replication and chromosomal relaxation, providing a foundation for determining how region-specific properties contribute to both chromosomal structure and gene regulation.

Crystal Structure of the COMPASS H3K4 Methyltransferase Catalytic Module.

The SET1/MLL family of histone methyltransferases is conserved in eukaryotes and regulates transcription by catalyzing histone H3K4 mono-, di-, and tri-methylation. These enzymes form a common five-subunit catalytic core whose assembly is critical for their basal and regulated enzymatic activities through unknown mechanisms. Here, we present the crystal structure of the intact yeast COMPASS histone methyltransferase catalytic module consisting of Swd1, Swd3, Bre2, Sdc1, and Set1. The complex is organized by Swd1, whose conserved C-terminal tail not only nucleates Swd3 and a Bre2-Sdc1 subcomplex, but also joins Set1 to construct a regulatory pocket next to the catalytic site. This inter-subunit pocket is targeted by a previously unrecognized enzyme-modulating motif in Swd3 and features a doorstop-style mechanism dictating substrate selectivity among SET1/MLL family members. By spatially mapping the functional components of COMPASS, our results provide a structural framework for understanding the multifaceted functions and regulation of the H3K4 methyltransferase family.

Sensory Neurons Co-opt Classical Immune Signaling Pathways to Mediate Chronic Itch.

Mammals have evolved neurophysiologic reflexes, such as coughing and scratching, to expel invading pathogens and noxious environmental stimuli. It is well established that these responses are also associated with chronic inflammatory diseases, including asthma and atopic dermatitis. However, the mechanisms by which inflammatory pathways promote sensations such as itch remain poorly understood. Here, we show that type 2 cytokines directly activate sensory neurons in both mice and humans. Further, we demonstrate that chronic itch is dependent on neuronal IL-4Rα and JAK1 signaling. We also observe that patients with recalcitrant chronic itch that failed other immunosuppressive therapies markedly improve when treated with JAK inhibitors. Thus, signaling mechanisms previously ascribed to the immune system may represent novel therapeutic targets within the nervous system. Collectively, this study reveals an evolutionarily conserved paradigm in which the sensory nervous system employs classical immune signaling pathways to influence mammalian behavior.

Targeted Apoptosis of Senescent Cells Restores Tissue Homeostasis in Response to Chemotoxicity and Aging.

The accumulation of irreparable cellular damage restricts healthspan after acute stress or natural aging. Senescent cells are thought to impair tissue function, and their genetic clearance can delay features of aging. Identifying how senescent cells avoid apoptosis allows for the prospective design of anti-senescence compounds to address whether homeostasis can also be restored. Here, we identify FOXO4 as a pivot in senescent cell viability. We designed a FOXO4 peptide that perturbs the FOXO4 interaction with p53. In senescent cells, this selectively causes p53 nuclear exclusion and cell-intrinsic apoptosis. Under conditions where it was well tolerated in vivo, this FOXO4 peptide neutralized doxorubicin-induced chemotoxicity. Moreover, it restored fitness, fur density, and renal function in both fast aging XpdTTD/TTD and naturally aged mice. Thus, therapeutic targeting of senescent cells is feasible under conditions where loss of health has already occurred, and in doing so tissue homeostasis can effectively be restored.

Ice-Binding Proteins and Their Function.

Ice-binding proteins (IBPs) are a diverse class of proteins that assist organism survival in the presence of ice in cold climates. They have different origins in many organisms, including bacteria, fungi, algae, diatoms, plants, insects, and fish. This review covers the gamut of IBP structures and functions and the common features they use to bind ice. We discuss mechanisms by which IBPs adsorb to ice and interfere with its growth, evidence for their irreversible association with ice, and methods for enhancing the activity of IBPs. The applications of IBPs in the food industry, in cryopreservation, and in other technologies are vast, and we chart out some possibilities.

Ileitis-associated tertiary lymphoid organs arise at lymphatic valves and impede mesenteric lymph flow in response to tumor necrosis factor.

Lymphangitis and the formation of tertiary lymphoid organs (TLOs) in the mesentery are features of Crohn's disease. Here, we examined the genesis of these TLOs and their impact on disease progression. Whole-mount and intravital imaging of the ileum and ileum-draining collecting lymphatic vessels (CLVs) draining to mesenteric lymph nodes from TNFΔARE mice, a model of ileitis, revealed TLO formation at valves of CLVs. TLOs obstructed cellular and molecular outflow from the gut and were sites of lymph leakage and backflow. Tumor necrosis factor (TNF) neutralization begun at early stages of TLO formation restored lymph transport. However, robustly developed, chronic TLOs resisted regression and restoration of flow after TNF neutralization. TNF stimulation of cultured lymphatic endothelial cells reprogrammed responses to oscillatory shear stress, preventing the induction of valve-associated genes. Disrupted transport of immune cells, driven by loss of valve integrity and TLO formation, may contribute to the pathology of Crohn's disease.

The Regulatory Factor ZFHX3 Modifies Circadian Function in SCN via an AT Motif-Driven Axis.

We identified a dominant missense mutation in the SCN transcription factor Zfhx3, termed short circuit (Zfhx3(Sci)), which accelerates circadian locomotor rhythms in mice. ZFHX3 regulates transcription via direct interaction with predicted AT motifs in target genes. The mutant protein has a decreased ability to activate consensus AT motifs in vitro. Using RNA sequencing, we found minimal effects on core clock genes in Zfhx3(Sci/+) SCN, whereas the expression of neuropeptides critical for SCN intercellular signaling was significantly disturbed. Moreover, mutant ZFHX3 had a decreased ability to activate AT motifs in the promoters of these neuropeptide genes. Lentiviral transduction of SCN slices showed that the ZFHX3-mediated activation of AT motifs is circadian, with decreased amplitude and robustness of these oscillations in Zfhx3(Sci/+) SCN slices. In conclusion, by cloning Zfhx3(Sci), we have uncovered a circadian transcriptional axis that determines the period and robustness of behavioral and SCN molecular rhythms.

The Cortical Motor Areas and the Emergence of Motor Skills: A Neuroanatomical Perspective.

What changes in neural architecture account for the emergence and expansion of dexterity in primates? Dexterity, or skill in performing motor tasks, depends on the ability to generate highly fractionated patterns of muscle activity. It also involves the spatiotemporal coordination of activity in proximal and distal muscles across multiple joints. Many motor skills require the generation of complex movement sequences that are only acquired and refined through extensive practice. Improvements in dexterity have enabled primates to manufacture and use tools and humans to engage in skilled motor behaviors such as typing, dance, musical performance, and sports. Our analysis leads to the following synthesis: The neural substrate that endows primates with their enhanced motor capabilities is due, in part, to (a) major organizational changes in the primary motor cortex and (b) the proliferation of output pathways from other areas of the cerebral cortex, especially from the motor areas on the medial wall of the hemisphere.

Axonally synthesized ATF4 transmits a neurodegenerative signal across brain regions.

In Alzheimer's disease (AD) brain, exposure of axons to Aß causes pathogenic changes that spread retrogradely by unknown mechanisms, affecting the entire neuron. We found that locally applied Aß1-42 initiates axonal synthesis of a defined set of proteins including the transcription factor ATF4. Inhibition of local translation and retrograde transport or knockdown of axonal Atf4 mRNA abolished Aß-induced ATF4 transcriptional activity and cell loss. Aß1-42 injection into the dentate gyrus (DG) of mice caused loss of forebrain neurons whose axons project to the DG. Protein synthesis and Atf4 mRNA were upregulated in these axons, and coinjection of Atf4 siRNA into the DG reduced the effects of Aß1-42 in the forebrain. ATF4 protein and transcripts were found with greater frequency in axons in the brain of AD patients. These results reveal an active role for intra-axonal translation in neurodegeneration and identify ATF4 as a mediator for the spread of AD pathology.

De novo evolution of macroscopic multicellularity.

While early multicellular lineages necessarily started out as relatively simple groups of cells, little is known about how they became Darwinian entities capable of sustained multicellular evolution1-3. Here we investigate this with a multicellularity long-term evolution experiment, selecting for larger group size in the snowflake yeast (Saccharomyces cerevisiae) model system. Given the historical importance of oxygen limitation4, our ongoing experiment consists of three metabolic treatments5-anaerobic, obligately aerobic and mixotrophic yeast. After 600 rounds of selection, snowflake yeast in the anaerobic treatment group evolved to be macroscopic, becoming around 2 × 104 times larger (approximately mm scale) and about 104-fold more biophysically tough, while retaining a clonal multicellular life cycle. This occurred through biophysical adaptation-evolution of increasingly elongate cells that initially reduced the strain of cellular packing and then facilitated branch entanglements that enabled groups of cells to stay together even after many cellular bonds fracture. By contrast, snowflake yeast competing for low oxygen5 remained microscopic, evolving to be only around sixfold larger, underscoring the critical role of oxygen levels in the evolution of multicellular size. Together, this research provides unique insights into an ongoing evolutionary transition in individuality, showing how simple groups of cells overcome fundamental biophysical limitations through gradual, yet sustained, multicellular evolution.

Shared and distinct mechanisms of UBA1 inactivation across different diseases.

Most cellular ubiquitin signaling is initiated by UBA1, which activates and transfers ubiquitin to tens of E2 enzymes. Clonally acquired UBA1 missense mutations cause an inflammatory-hematologic overlap disease called VEXAS (vacuoles, E1, X-linked, autoinflammatory, somatic) syndrome. Despite extensive clinical investigation into this lethal disease, little is known about the underlying molecular mechanisms. Here, by dissecting VEXAS-causing UBA1 mutations, we discovered that p.Met41 mutations alter cytoplasmic isoform expression, whereas other mutations reduce catalytic activity of nuclear and cytoplasmic isoforms by diverse mechanisms, including aberrant oxyester formation. Strikingly, non-p.Met41 mutations most prominently affect transthioesterification, revealing ubiquitin transfer to cytoplasmic E2 enzymes as a shared property of pathogenesis amongst different VEXAS syndrome genotypes. A similar E2 charging bottleneck exists in some lung cancer-associated UBA1 mutations, but not in spinal muscular atrophy-causing UBA1 mutations, which instead, render UBA1 thermolabile. Collectively, our results highlight the precision of conformational changes required for faithful ubiquitin transfer, define distinct and shared mechanisms of UBA1 inactivation in diverse diseases, and suggest that specific E1-E2 modules control different aspects of tissue differentiation and maintenance.

Deep physical neural networks trained with backpropagation.

Deep-learning models have become pervasive tools in science and engineering. However, their energy requirements now increasingly limit their scalability1. Deep-learning accelerators2-9 aim to perform deep learning energy-efficiently, usually targeting the inference phase and often by exploiting physical substrates beyond conventional electronics. Approaches so far10-22 have been unable to apply the backpropagation algorithm to train unconventional novel hardware in situ. The advantages of backpropagation have made it the de facto training method for large-scale neural networks, so this deficiency constitutes a major impediment. Here we introduce a hybrid in situ-in silico algorithm, called physics-aware training, that applies backpropagation to train controllable physical systems. Just as deep learning realizes computations with deep neural networks made from layers of mathematical functions, our approach allows us to train deep physical neural networks made from layers of controllable physical systems, even when the physical layers lack any mathematical isomorphism to conventional artificial neural network layers. To demonstrate the universality of our approach, we train diverse physical neural networks based on optics, mechanics and electronics to experimentally perform audio and image classification tasks. Physics-aware training combines the scalability of backpropagation with the automatic mitigation of imperfections and noise achievable with in situ algorithms. Physical neural networks have the potential to perform machine learning faster and more energy-efficiently than conventional electronic processors and, more broadly, can endow physical systems with automatically designed physical functionalities, for example, for robotics23-26, materials27-29 and smart sensors30-32.

The dawn of virtual cell biology.

To fulfill systems biology's promise of providing fundamental new insights will require the development of quantitative and predictive models of whole cells. In this issue, Karr et al. present the first integrated and dynamic computational model of a bacterium that accounts for all of its components and their interactions.

Personal omics profiling reveals dynamic molecular and medical phenotypes.

Personalized medicine is expected to benefit from combining genomic information with regular monitoring of physiological states by multiple high-throughput methods. Here, we present an integrative personal omics profile (iPOP), an analysis that combines genomic, transcriptomic, proteomic, metabolomic, and autoantibody profiles from a single individual over a 14 month period. Our iPOP analysis revealed various medical risks, including type 2 diabetes. It also uncovered extensive, dynamic changes in diverse molecular components and biological pathways across healthy and diseased conditions. Extremely high-coverage genomic and transcriptomic data, which provide the basis of our iPOP, revealed extensive heteroallelic changes during healthy and diseased states and an unexpected RNA editing mechanism. This study demonstrates that longitudinal iPOP can be used to interpret healthy and diseased states by connecting genomic information with additional dynamic omics activity.

Large-Scale, Quantitative Protein Assays on a High-Throughput DNA Sequencing Chip.

High-throughput DNA sequencing techniques have enabled diverse approaches for linking DNA sequence to biochemical function. In contrast, assays of protein function have substantial limitations in terms of throughput, automation, and widespread availability. We have adapted an Illumina high-throughput sequencing chip to display an immense diversity of ribosomally translated proteins and peptides and then carried out fluorescence-based functional assays directly on this flow cell, demonstrating that a single, widely available high-throughput platform can perform both sequencing-by-synthesis and protein assays. We quantified the binding of the M2 anti-FLAG antibody to a library of 1.3 × 104 variant FLAG peptides, exploring non-additive effects of combinations of mutations and discovering a "superFLAG" epitope variant. We also measured the enzymatic activity of 1.56 × 105 molecular variants of full-length human O6-alkylguanine-DNA alkyltransferase (SNAP-tag). This comprehensive corpus of catalytic rates revealed amino acid interaction networks and cooperativity, linked positive cooperativity to structural proximity, and revealed ubiquitous positively cooperative interactions with histidine residues.

Structural Basis of H2B Ubiquitination-Dependent H3K4 Methylation by COMPASS.

The COMPASS (complex of proteins associated with Set1) complex represents the prototype of the SET1/MLL family of methyltransferases that controls gene transcription by H3K4 methylation (H3K4me). Although H2B monoubiquitination (H2Bub) is well known as a prerequisite histone mark for COMPASS activity, how H2Bub activates COMPASS remains unclear. Here, we report the cryoelectron microscopy (cryo-EM) structures of an extended COMPASS catalytic module (CM) bound to the H2Bub and free nucleosome. The COMPASS CM clamps onto the nucleosome disk-face via an extensive interface to capture the flexible H3 N-terminal tail. The interface also sandwiches a critical Set1 arginine-rich motif (ARM) that autoinhibits COMPASS. Unexpectedly, without enhancing COMPASS-nucleosome interaction, H2Bub activates the enzymatic assembly by packing against Swd1 and alleviating the inhibitory effect of the Set1 ARM upon fastening it to the acidic patch. By delineating the spatial configuration of the COMPASS-H2Bub-nucleosome assembly, our studies establish the structural framework for understanding the long-studied H2Bub-H3K4me histone modification crosstalk.

Introducing a Spectrum of Long-Range Genomic Deletions in Human Embryonic Stem Cells Using Type I CRISPR-Cas.

CRISPR-Cas systems enable microbial adaptive immunity and provide eukaryotic genome editing tools. These tools employ a single effector enzyme of type II or V CRISPR to generate RNA-guided, precise genome breaks. Here we demonstrate the feasibility of using type I CRISPR-Cas to effectively introduce a spectrum of long-range chromosomal deletions with a single RNA guide in human embryonic stem cells and HAP1 cells. Type I CRISPR systems rely on the multi-subunit ribonucleoprotein (RNP) complex Cascade to identify DNA targets and on the helicase-nuclease enzyme Cas3 to degrade DNA processively. With RNP delivery of T. fusca Cascade and Cas3, we obtained 13%-60% editing efficiency. Long-range PCR-based and high-throughput-sequencing-based lesion analyses reveal that a variety of deletions, ranging from a few hundred base pairs to 100 kilobases, are created upstream of the target site. These results highlight the potential utility of type I CRISPR-Cas for long-range genome manipulations and deletion screens in eukaryotes.

Muscle eosinophilia is a hallmark of chronic disease in facioscapulohumeral muscular dystrophy.

Facioscapulohumeral muscular dystrophy (FSHD) is a progressive myopathy caused by the aberrant increased expression of the DUX4 retrogene in skeletal muscle cells. The DUX4 gene encodes a transcription factor that functions in zygotic genome activation and then is silenced in most adult somatic tissues. DUX4 expression in FSHD disrupts normal muscle cell function; however, the downstream pathogenic mechanisms are still unclear. Histologically, FSHD affected muscles show a characteristic dystrophic phenotype that is often accompanied by a pronounced immune cell infiltration, but the role of the immune system in FSHD is not understood. Previously, we used ACTA1;FLExDUX4 FSHD-like mouse models varying in severity as discovery tools to identify increased Interleukin 6 and microRNA-206 levels as serum biomarkers for FSHD disease severity. In this study, we use the ACTA1;FLExDUX4 chronic FSHD-like mouse model to provide insight into the immune response to DUX4 expression in skeletal muscles. We demonstrate that these FSHD-like muscles are enriched with the chemoattractant eotaxin and the cytotoxic eosinophil peroxidase, and exhibit muscle eosinophilia. We further identified muscle fibers with positive staining for eosinophil peroxidase in human FSHD muscle. Our data supports that skeletal muscle eosinophilia is a hallmark of FSHD pathology.

Regional and bilateral MRI and gene signatures in facioscapulohumeral dystrophy: implications for clinical trial design and mechanisms of disease progression.

Identifying the aberrant expression of DUX4 in skeletal muscle as the cause of facioscapulohumeral dystrophy (FSHD) has led to rational therapeutic development and clinical trials. Several studies support the use of MRI characteristics and the expression of DUX4-regulated genes in muscle biopsies as biomarkers of FSHD disease activity and progression. We performed lower-extremity MRI and muscle biopsies in the mid-portion of the tibialis anterior (TA) muscles bilaterally in FSHD subjects and validated our prior reports of the strong association between MRI characteristics and expression of genes regulated by DUX4 and other gene categories associated with FSHD disease activity. We further show that measurements of normalized fat content in the entire TA muscle strongly predict molecular signatures in the mid-portion of the TA, indicating that regional biopsies can accurately measure progression in the whole muscle and providing a strong basis for inclusion of MRI and molecular biomarkers in clinical trial design. An unanticipated finding was the strong correlations of molecular signatures in the bilateral comparisons, including markers of B-cells and other immune cell populations, suggesting that a systemic immune cell infiltration of skeletal muscle might have a role in disease progression.