Prediction of efficiencies for diverse prime editing systems in multiple cell types.

Applications of prime editing are often limited due to insufficient efficiencies, and it can require substantial time and resources to determine the most efficient pegRNAs and prime editors (PEs) to generate a desired edit under various experimental conditions. Here, we evaluated prime editing efficiencies for a total of 338,996 pairs of pegRNAs including 3,979 epegRNAs and target sequences in an error-free manner. These datasets enabled a systematic determination of factors affecting prime editing efficiencies. Then, we developed computational models, named DeepPrime and DeepPrime-FT, that can predict prime editing efficiencies for eight prime editing systems in seven cell types for all possible types of editing of up to 3 base pairs. We also extensively profiled the prime editing efficiencies at mismatched targets and developed a computational model predicting editing efficiencies at such targets. These computational models, together with our improved knowledge about prime editing efficiency determinants, will greatly facilitate prime editing applications.

Local Axonal Conduction Shapes the Spatiotemporal Properties of Neural Sequences.

Sequential activation of neurons has been observed during various behavioral and cognitive processes, but the underlying circuit mechanisms remain poorly understood. Here, we investigate premotor sequences in HVC (proper name) of the adult zebra finch forebrain that are central to the performance of the temporally precise courtship song. We use high-density silicon probes to measure song-related population activity, and we compare these observations with predictions from a range of network models. Our results support a circuit architecture in which heterogeneous delays between sequentially active neurons shape the spatiotemporal patterns of HVC premotor neuron activity. We gauge the impact of several delay sources, and we find the primary contributor to be slow conduction through axonal collaterals within HVC, which typically adds between 1 and 7.5 ms for each link within the sequence. Thus, local axonal "delay lines" can play an important role in determining the dynamical repertoire of neural circuits.

A Species-Wide Inventory of NLR Genes and Alleles in Arabidopsis thaliana.

Infectious disease is both a major force of selection in nature and a prime cause of yield loss in agriculture. In plants, disease resistance is often conferred by nucleotide-binding leucine-rich repeat (NLR) proteins, intracellular immune receptors that recognize pathogen proteins and their effects on the host. Consistent with extensive balancing and positive selection, NLRs are encoded by one of the most variable gene families in plants, but the true extent of intraspecific NLR diversity has been unclear. Here, we define a nearly complete species-wide pan-NLRome in Arabidopsis thaliana based on sequence enrichment and long-read sequencing. The pan-NLRome largely saturates with approximately 40 well-chosen wild strains, with half of the pan-NLRome being present in most accessions. We chart NLR architectural diversity, identify new architectures, and quantify selective forces that act on specific NLRs and NLR domains. Our study provides a blueprint for defining pan-NLRomes.

Yeast Ataxin-2 Forms an Intracellular Condensate Required for the Inhibition of TORC1 Signaling during Respiratory Growth.

Yeast ataxin-2, also known as Pbp1 (polyA binding protein-binding protein 1), is an intrinsically disordered protein implicated in stress granule formation, RNA biology, and neurodegenerative disease. To understand the endogenous function of this protein, we identify Pbp1 as a dedicated regulator of TORC1 signaling and autophagy under conditions that require mitochondrial respiration. Pbp1 binds to TORC1 specifically during respiratory growth, but utilizes an additional methionine-rich, low complexity (LC) region to inhibit TORC1. This LC region causes phase separation, forms reversible fibrils, and enables self-association into assemblies required for TORC1 inhibition. Mutants that weaken phase separation in vitro exhibit reduced capacity to inhibit TORC1 and induce autophagy. Loss of Pbp1 leads to mitochondrial dysfunction and reduced fitness during nutritional stress. Thus, Pbp1 forms a condensate in response to respiratory status to regulate TORC1 signaling.

Redox State Controls Phase Separation of the Yeast Ataxin-2 Protein via Reversible Oxidation of Its Methionine-Rich Low-Complexity Domain.

Yeast ataxin-2, also known as Pbp1, senses the activity state of mitochondria in order to regulate TORC1. A domain of Pbp1 required to adapt cells to mitochondrial activity is of low sequence complexity. The low-complexity (LC) domain of Pbp1 forms labile, cross-ß polymers that facilitate phase transition of the protein into liquid-like or gel-like states. Phase transition for other LC domains is reliant upon widely distributed aromatic amino acids. In place of tyrosine or phenylalanine residues prototypically used for phase separation, Pbp1 contains 24 similarly disposed methionine residues. Here, we show that the Pbp1 methionine residues are sensitive to hydrogen peroxide (H2O2)-mediated oxidation in vitro and in living cells. Methionine oxidation melts Pbp1 liquid-like droplets in a manner reversed by methionine sulfoxide reductase enzymes. These observations explain how reversible formation of labile polymers by the Pbp1 LC domain enables the protein to function as a sensor of cellular redox state.

Combinatorial Genetics Reveals a Scaling Law for the Effects of Mutations on Splicing.

Despite a wealth of molecular knowledge, quantitative laws for accurate prediction of biological phenomena remain rare. Alternative pre-mRNA splicing is an important regulated step in gene expression frequently perturbed in human disease. To understand the combined effects of mutations during evolution, we quantified the effects of all possible combinations of exonic mutations accumulated during the emergence of an alternatively spliced human exon. This revealed that mutation effects scale non-monotonically with the inclusion level of an exon, with each mutation having maximum effect at a predictable intermediate inclusion level. This scaling is observed genome-wide for cis and trans perturbations of splicing, including for natural and disease-associated variants. Mathematical modeling suggests that competition between alternative splice sites is sufficient to cause this non-linearity in the genotype-phenotype map. Combining the global scaling law with specific pairwise interactions between neighboring mutations allows accurate prediction of the effects of complex genotype changes involving >10 mutations.

Structures of the Catalytically Activated Yeast Spliceosome Reveal the Mechanism of Branching.

Pre-mRNA splicing is executed by the spliceosome. Structural characterization of the catalytically activated complex (B∗) is pivotal for understanding the branching reaction. In this study, we assembled the B∗ complexes on two different pre-mRNAs from Saccharomyces cerevisiae and determined the cryo-EM structures of four distinct B∗ complexes at overall resolutions of 2.9-3.8 Å. The duplex between U2 small nuclear RNA (snRNA) and the branch point sequence (BPS) is discretely away from the 5'-splice site (5'SS) in the three B∗ complexes that are devoid of the step I splicing factors Yju2 and Cwc25. Recruitment of Yju2 into the active site brings the U2/BPS duplex into the vicinity of 5'SS, with the BPS nucleophile positioned 4 Å away from the catalytic metal M2. This analysis reveals the functional mechanism of Yju2 and Cwc25 in branching. These structures on different pre-mRNAs reveal substrate-specific conformations of the spliceosome in a major functional state.

Characterizing the Major Structural Variant Alleles of the Human Genome.

In order to provide a comprehensive resource for human structural variants (SVs), we generated long-read sequence data and analyzed SVs for fifteen human genomes. We sequence resolved 99,604 insertions, deletions, and inversions including 2,238 (1.6 Mbp) that are shared among all discovery genomes with an additional 13,053 (6.9 Mbp) present in the majority, indicating minor alleles or errors in the reference. Genotyping in 440 additional genomes confirms the most common SVs in unique euchromatin are now sequence resolved. We report a ninefold SV bias toward the last 5 Mbp of human chromosomes with nearly 55% of all VNTRs (variable number of tandem repeats) mapping to this portion of the genome. We identify SVs affecting coding and noncoding regulatory loci improving annotation and interpretation of functional variation. These data provide the framework to construct a canonical human reference and a resource for developing advanced representations capable of capturing allelic diversity.

Rational Engineering of XCaMPs, a Multicolor GECI Suite for In Vivo Imaging of Complex Brain Circuit Dynamics.

To decipher dynamic brain information processing, current genetically encoded calcium indicators (GECIs) are limited in single action potential (AP) detection speed, combinatorial spectral compatibility, and two-photon imaging depth. To address this, here, we rationally engineered a next-generation quadricolor GECI suite, XCaMPs. Single AP detection was achieved within 3-10 ms of spike onset, enabling measurements of fast-spike trains in parvalbumin (PV)-positive interneurons in the barrel cortex in vivo and recording three distinct (two inhibitory and one excitatory) ensembles during pre-motion activity in freely moving mice. In vivo paired recording of pre- and postsynaptic firing revealed spatiotemporal constraints of dendritic inhibition in layer 1 in vivo, between axons of somatostatin (SST)-positive interneurons and apical tufts dendrites of excitatory pyramidal neurons. Finally, non-invasive, subcortical imaging using red XCaMP-R uncovered somatosensation-evoked persistent activity in hippocampal CA1 neurons. Thus, the XCaMPs offer a critical enhancement of solution space in studies of complex neuronal circuit dynamics. VIDEO ABSTRACT.

Protein Sequence Editing of SKN-1A/Nrf1 by Peptide:N-Glycanase Controls Proteasome Gene Expression.

The proteasome mediates selective protein degradation and is dynamically regulated in response to proteotoxic challenges. SKN-1A/Nrf1, an endoplasmic reticulum (ER)-associated transcription factor that undergoes N-linked glycosylation, serves as a sensor of proteasome dysfunction and triggers compensatory upregulation of proteasome subunit genes. Here, we show that the PNG-1/NGLY1 peptide:N-glycanase edits the sequence of SKN-1A protein by converting particular N-glycosylated asparagine residues to aspartic acid. Genetically introducing aspartates at these N-glycosylation sites bypasses the requirement for PNG-1/NGLY1, showing that protein sequence editing rather than deglycosylation is key to SKN-1A function. This pathway is required to maintain sufficient proteasome expression and activity, and SKN-1A hyperactivation confers resistance to the proteotoxicity of human amyloid beta peptide. Deglycosylation-dependent protein sequence editing explains how ER-associated and cytosolic isoforms of SKN-1 perform distinct cytoprotective functions corresponding to those of mammalian Nrf1 and Nrf2. Thus, we uncover an unexpected mechanism by which N-linked glycosylation regulates protein function and proteostasis.

Optogenetic Editing Reveals the Hierarchical Organization of Learned Action Sequences.

The organization of action into sequences underlies complex behaviors that are essential for organismal survival and reproduction. Despite extensive studies of innate sequences in relation to central pattern generators, how learned action sequences are controlled and whether they are organized as a chain or a hierarchy remain largely unknown. By training mice to perform heterogeneous action sequences, we demonstrate that striatal direct and indirect pathways preferentially encode different behavioral levels of sequence structure. State-dependent closed-loop optogenetic stimulation of the striatal direct pathway can selectively insert a single action element into the sequence without disrupting the overall sequence length. Optogenetic manipulation of the striatal indirect pathway completely removes the ongoing subsequence while leaving the following subsequence to be executed with the appropriate timing and length. These results suggest that learned action sequences are not organized in a serial but rather a hierarchical structure that is distinctly controlled by basal ganglia pathways.

Dicer-like Enzymes with Sequence Cleavage Preferences.

Dicer proteins are known to produce small RNAs (sRNAs) from long double-stranded RNA (dsRNA) templates. These sRNAs are bound by Argonaute proteins, which select the guide strand, often with a 5' end sequence bias. However, Dicer proteins have never been shown to have sequence cleavage preferences. In Paramecium development, two classes of sRNAs that are required for DNA elimination are produced by three Dicer-like enzymes: Dcl2, Dcl3, and Dcl5. Through in vitro cleavage assays, we demonstrate that Dcl2 has a strict size preference for 25 nt and a sequence preference for 5' U and 5' AGA, while Dcl3 has a sequence preference for 5' UNG. Dcl5, however, has cleavage preferences for 5' UAG and 3' CUAC/UN, which leads to the production of RNAs precisely matching short excised DNA elements with corresponding end base preferences. Thus, we characterize three Dicer-like enzymes that are involved in Paramecium development and propose a biological role for their sequence-biased cleavage products.

Site-Specific Self-Catalyzed DNA Depurination: A Biological Mechanism That Leads to Mutations and Creates Sequence Diversity.

Self-catalyzed DNA depurination is a sequence-specific physiological mechanism mediated by spontaneous extrusion of a stem-loop catalytic intermediate. Hydrolysis of the 5'G residue of the 5'GA/TGG loop and of the first 5'A residue of the 5'GAGA loop, together with particular first stem base pairs, specifies their hydrolysis without involving protein, cofactor, or cation. As such, this mechanism is the only known DNA catalytic activity exploited by nature. The consensus sequences for self-depurination of such G- and A-loop residues occur in all genomes examined across the phyla, averaging one site every 2,000-4,000 base pairs. Because apurinic sites are subject to error-prone repair, leading to substitution and short frameshift mutations, they are both a source of genome damage and a means for creating sequence diversity. Their marked overrepresentation in genomes, and largely unchanging density from the lowest to the highest organisms, indicate their selection over the course of evolution. The mutagenicity at such sites in many human genes is associated with loss of function of key proteins responsible for diverse diseases.

HLA-II immunopeptidome profiling and deep learning reveal features of antigenicity to inform antigen discovery.

CD4+ T cell responses are exquisitely antigen specific and directed toward peptide epitopes displayed by human leukocyte antigen class II (HLA-II) on antigen-presenting cells. Underrepresentation of diverse alleles in ligand databases and an incomplete understanding of factors affecting antigen presentation in vivo have limited progress in defining principles of peptide immunogenicity. Here, we employed monoallelic immunopeptidomics to identify 358,024 HLA-II binders, with a particular focus on HLA-DQ and HLA-DP. We uncovered peptide-binding patterns across a spectrum of binding affinities and enrichment of structural antigen features. These aspects underpinned the development of context-aware predictor of T cell antigens (CAPTAn), a deep learning model that predicts peptide antigens based on their affinity to HLA-II and full sequence of their source proteins. CAPTAn was instrumental in discovering prevalent T cell epitopes from bacteria in the human microbiome and a pan-variant epitope from SARS-CoV-2. Together CAPTAn and associated datasets present a resource for antigen discovery and the unraveling genetic associations of HLA alleles with immunopathologies.

Structure of FUS Protein Fibrils and Its Relevance to Self-Assembly and Phase Separation of Low-Complexity Domains.

Polymerization and phase separation of proteins containing low-complexity (LC) domains are important factors in gene expression, mRNA processing and trafficking, and localization of translation. We have used solid-state nuclear magnetic resonance methods to characterize the molecular structure of self-assembling fibrils formed by the LC domain of the fused in sarcoma (FUS) RNA-binding protein. From the 214-residue LC domain of FUS (FUS-LC), a segment of only 57 residues forms the fibril core, while other segments remain dynamically disordered. Unlike pathogenic amyloid fibrils, FUS-LC fibrils lack hydrophobic interactions within the core and are not polymorphic at the molecular structural level. Phosphorylation of core-forming residues by DNA-dependent protein kinase blocks binding of soluble FUS-LC to FUS-LC hydrogels and dissolves phase-separated, liquid-like FUS-LC droplets. These studies offer a structural basis for understanding LC domain self-assembly, phase separation, and regulation by post-translational modification.

Interleukin-33-activated neuropeptide CGRP-producing memory Th2 cells cooperate with somatosensory neurons to induce conjunctival itch.

Allergic conjunctivitis is a chronic inflammatory disease that is characterized by severe itch in the conjunctiva, but how neuro-immune interactions shape the pathogenesis of severe itch remains unclear. We identified a subset of memory-type pathogenic Th2 cells that preferentially expressed Il1rl1-encoding ST2 and Calca-encoding calcitonin-gene-related peptide (CGRP) in the inflammatory conjunctiva using a single-cell analysis. The IL-33-ST2 axis in memory Th2 cells controlled the axonal elongation of the peripheral sensory C-fiber and the induction of severe itch. Pharmacological blockade and genetic deletion of CGRP signaling in vivo attenuated scratching behavior. The analysis of giant papillae from patients with severe allergic conjunctivitis revealed ectopic lymphoid structure formation with the accumulation of IL-33-producing epithelial cells and CGRP-producing pathogenic CD4+ T cells accompanied by peripheral nerve elongation. Thus, the IL-33-ST2-CGRP axis directs severe itch with neuro-reconstruction in the inflammatory conjunctiva and is a potential therapeutic target for severe itch in allergic conjunctivitis.

A large-scale systematic survey reveals recurring molecular features of public antibody responses to SARS-CoV-2.

Global research to combat the COVID-19 pandemic has led to the isolation and characterization of thousands of human antibodies to the SARS-CoV-2 spike protein, providing an unprecedented opportunity to study the antibody response to a single antigen. Using the information derived from 88 research publications and 13 patents, we assembled a dataset of ∼8,000 human antibodies to the SARS-CoV-2 spike protein from >200 donors. By analyzing immunoglobulin V and D gene usages, complementarity-determining region H3 sequences, and somatic hypermutations, we demonstrated that the common (public) responses to different domains of the spike protein were quite different. We further used these sequences to train a deep-learning model to accurately distinguish between the human antibodies to SARS-CoV-2 spike protein and those to influenza hemagglutinin protein. Overall, this study provides an informative resource for antibody research and enhances our molecular understanding of public antibody responses.

PAM-Dependent Target DNA Recognition and Cleavage by C2c1 CRISPR-Cas Endonuclease.

C2c1 is a newly identified guide RNA-mediated type V-B CRISPR-Cas endonuclease that site-specifically targets and cleaves both strands of target DNA. We have determined crystal structures of Alicyclobacillus acidoterrestris C2c1 (AacC2c1) bound to sgRNA as a binary complex and to target DNAs as ternary complexes, thereby capturing catalytically competent conformations of AacC2c1 with both target and non-target DNA strands independently positioned within a single RuvC catalytic pocket. Moreover, C2c1-mediated cleavage results in a staggered seven-nucleotide break of target DNA. crRNA adopts a pre-ordered five-nucleotide A-form seed sequence in the binary complex, with release of an inserted tryptophan, facilitating zippering up of 20-bp guide RNA:target DNA heteroduplex on ternary complex formation. Notably, the PAM-interacting cleft adopts a "locked" conformation on ternary complex formation. Structural comparison of C2c1 ternary complexes with their Cas9 and Cpf1 counterparts highlights the diverse mechanisms adopted by these distinct CRISPR-Cas systems, thereby broadening and enhancing their applicability as genome editing tools.

Toxic PR Poly-Dipeptides Encoded by the C9orf72 Repeat Expansion Target LC Domain Polymers.

Two complementary approaches were used in search of the intracellular targets of the toxic PR poly-dipeptide encoded by the repeat sequences expanded in the C9orf72 form of amyotrophic lateral sclerosis. The top categories of PRn-bound proteins include constituents of non-membrane invested cellular organelles and intermediate filaments. PRn targets are enriched for the inclusion of low complexity (LC) sequences. Evidence is presented indicating that LC sequences represent the direct target of PRn binding and that interaction between the PRn poly-dipeptide and LC domains is polymer-dependent. These studies indicate that PRn-mediated toxicity may result from broad impediments to the dynamics of cell structure and information flow from gene to message to protein.

U1 snRNP increases RNA Pol II elongation rate to enable synthesis of long genes.

The expansion of introns within mammalian genomes poses a challenge for the production of full-length messenger RNAs (mRNAs), with increasing evidence that these long AT-rich sequences present obstacles to transcription. Here, we investigate RNA polymerase II (RNAPII) elongation at high resolution in mammalian cells and demonstrate that RNAPII transcribes faster across introns. Moreover, we find that this acceleration requires the association of U1 snRNP (U1) with the elongation complex at 5' splice sites. The role of U1 to stimulate elongation rate through introns reduces the frequency of both premature termination and transcriptional arrest, thereby dramatically increasing RNA production. We further show that changes in RNAPII elongation rate due to AT content and U1 binding explain previous reports of pausing or termination at splice junctions and the edge of CpG islands. We propose that U1-mediated acceleration of elongation has evolved to mitigate the risks that long AT-rich introns pose to transcript completion.