TemStaPro: protein thermostability prediction using sequence representations from protein language models.

MOTIVATION: Reliable prediction of protein thermostability from its sequence is valuable for both academic and industrial research. This prediction problem can be tackled using machine learning and by taking advantage of the recent blossoming of deep learning methods for sequence analysis. These methods can facilitate training on more data and, possibly, enable the development of more versatile thermostability predictors for multiple ranges of temperatures. RESULTS: We applied the principle of transfer learning to predict protein thermostability using embeddings generated by protein language models (pLMs) from an input protein sequence. We used large pLMs that were pre-trained on hundreds of millions of known sequences. The embeddings from such models allowed us to efficiently train and validate a high-performing prediction method using over one million sequences that we collected from organisms with annotated growth temperatures. Our method, TemStaPro (Temperatures of Stability for Proteins), was used to predict thermostability of CRISPR-Cas Class II effector proteins (C2EPs). Predictions indicated sharp differences among groups of C2EPs in terms of thermostability and were largely in tune with previously published and our newly obtained experimental data. AVAILABILITY AND IMPLEMENTATION: TemStaPro software and the related data are freely available from https://github.com/ievapudz/TemStaPro and https://doi.org/10.5281/zenodo.7743637.

Progress and harmonization of gene editing to treat human diseases: Proceeding of COST Action CA21113 GenE-HumDi.

The European Cooperation in Science and Technology (COST) is an intergovernmental organization dedicated to funding and coordinating scientific and technological research in Europe, fostering collaboration among researchers and institutions across countries. Recently, COST Action funded the "Genome Editing to treat Human Diseases" (GenE-HumDi) network, uniting various stakeholders such as pharmaceutical companies, academic institutions, regulatory agencies, biotech firms, and patient advocacy groups. GenE-HumDi's primary objective is to expedite the application of genome editing for therapeutic purposes in treating human diseases. To achieve this goal, GenE-HumDi is organized in several working groups, each focusing on specific aspects. These groups aim to enhance genome editing technologies, assess delivery systems, address safety concerns, promote clinical translation, and develop regulatory guidelines. The network seeks to establish standard procedures and guidelines for these areas to standardize scientific practices and facilitate knowledge sharing. Furthermore, GenE-HumDi aims to communicate its findings to the public in accessible yet rigorous language, emphasizing genome editing's potential to revolutionize the treatment of many human diseases. The inaugural GenE-HumDi meeting, held in Granada, Spain, in March 2023, featured presentations from experts in the field, discussing recent breakthroughs in delivery methods, safety measures, clinical translation, and regulatory aspects related to gene editing.

A new family of CRISPR-type V nucleases with C-rich PAM recognition.

Most CRISPR-type V nucleases are stimulated to cleave double-stranded (ds) DNA targets by a T-rich PAM, which restricts their targeting range. Here, we identify and characterize a new family of type V RNA-guided nuclease, Cas12l, that exclusively recognizes a C-rich (5'-CCY-3') PAM. The organization of genes within its CRISPR locus is similar to type II-B CRISPR-Cas9 systems, but both sequence analysis and functional studies establish it as a new family of type V effector. Biochemical experiments show that Cas12l nucleases function optimally between 37 and 52°C, depending on the ortholog, and preferentially cut supercoiled DNA. Like other type V nucleases, it exhibits collateral nonspecific ssDNA and ssRNA cleavage activity that is triggered by ssDNA or dsDNA target recognition. Finally, we show that one family member, Asp2Cas12l, functions in a heterologous cellular environment, altogether, suggesting that this new group of CRISPR-associated nucleases may be harnessed as genome editing reagents.

Diversification of the CRISPR Toolbox: Applications of CRISPR-Cas Systems Beyond Genome Editing.

The discovery of CRISPR has revolutionized the field of genome engineering, but the potential of this technology is far from reaching its limits. In this review, we explore the broad range of applications of CRISPR technology to highlight the rapid expansion of the field beyond gene editing alone. It has been demonstrated that CRISPR technology can control gene expression, spatiotemporally image the genome in vivo, and detect specific nucleic acid sequences for diagnostics. In addition, new technologies are under development to improve CRISPR quality controls for gene editing, thereby improving the reliability of these technologies for therapeutics and beyond. These are just some of the many CRISPR tools that have been developed in recent years, and the toolbox continues to diversify.

Discrimination of single-point mutations in unamplified genomic DNA via Cas9 immobilized on a graphene field-effect transistor.

Simple and fast methods for the detection of target genes with single-nucleotide specificity could open up genetic research and diagnostics beyond laboratory settings. We recently reported a biosensor for the electronic detection of unamplified target genes using liquid-gated graphene field-effect transistors employing an RNA-guided catalytically deactivated CRISPR-associated protein 9 (Cas9) anchored to a graphene monolayer. Here, using unamplified genomic samples from patients and by measuring multiple types of electrical response, we show that the biosensors can discriminate within one hour between wild-type and homozygous mutant alleles differing by a single nucleotide. We also show that biosensors using a guide RNA-Cas9 orthologue complex targeting genes within the protospacer-adjacent motif discriminated between homozygous and heterozygous DNA samples from patients with sickle cell disease, and that the biosensors can also be used to rapidly screen for guide RNA-Cas9 complexes that maximize gene-targeting efficiency.

A catalogue of biochemically diverse CRISPR-Cas9 orthologs.

Bacterial Cas9 nucleases from type II CRISPR-Cas antiviral defence systems have been repurposed as genome editing tools. Although these proteins are found in many microbes, only a handful of variants are used for these applications. Here, we use bioinformatic and biochemical analyses to explore this largely uncharacterized diversity. We apply cell-free biochemical screens to assess the protospacer adjacent motif (PAM) and guide RNA (gRNA) requirements of 79 Cas9 proteins, thus identifying at least 7 distinct gRNA classes and 50 different PAM sequence requirements. PAM recognition spans the entire spectrum of T-, A-, C-, and G-rich nucleotides, from single nucleotide recognition to sequence strings longer than 4 nucleotides. Characterization of a subset of Cas9 orthologs using purified components reveals additional biochemical diversity, including both narrow and broad ranges of temperature dependence, staggered-end DNA target cleavage, and a requirement for long stretches of homology between gRNA and DNA target. Our results expand the available toolset of RNA-programmable CRISPR-associated nucleases.

5' modifications to CRISPR-Cas9 gRNA can change the dynamics and size of R-loops and inhibit DNA cleavage.

A key aim in exploiting CRISPR-Cas is gRNA engineering to introduce additional functionalities, ranging from individual nucleotide changes that increase efficiency of on-target binding to the inclusion of larger functional RNA aptamers or ribonucleoproteins (RNPs). Cas9-gRNA interactions are crucial for complex assembly, but several distinct regions of the gRNA are amenable to modification. We used in vitro ensemble and single-molecule assays to assess the impact of gRNA structural alterations on RNP complex formation, R-loop dynamics, and endonuclease activity. Our results indicate that RNP formation was unaffected by any of our modifications. R-loop formation and DNA cleavage activity were also essentially unaffected by modification of the Upper Stem, first Hairpin and 3' end. In contrast, we found that 5' additions of only two or three nucleotides could reduce R-loop formation and cleavage activity of the RuvC domain relative to a single nucleotide addition. Such modifications are a common by-product of in vitro transcribed gRNA. We also observed that addition of a 20 nt RNA hairpin to the 5' end of a gRNA still supported RNP formation but produced a stable ∼9 bp R-loop that could not activate DNA cleavage. Consideration of these observations will assist in successful gRNA design.

Structure of the DNA-Bound Spacer Capture Complex of a Type II CRISPR-Cas System.

CRISPR and associated Cas proteins function as an adaptive immune system in prokaryotes to combat bacteriophage infection. During the immunization step, new spacers are acquired by the CRISPR machinery, but the molecular mechanism of spacer capture remains enigmatic. We show that the Cas9, Cas1, Cas2, and Csn2 proteins of a Streptococcus thermophilus type II-A CRISPR-Cas system form a complex and provide cryoelectron microscopy (cryo-EM) structures of three different assemblies. The predominant form, with the stoichiometry Cas18-Cas24-Csn28, referred to as monomer, contains ∼30 bp duplex DNA bound along a central channel. A minor species, termed a dimer, comprises two monomers that sandwich a further eight Cas1 and four Cas2 subunits and contains two DNA ∼30-bp duplexes within the channel. A filamentous form also comprises Cas18-Cas24-Csn28 units (typically 2-6) but with a different Cas1-Cas2 interface between them and a continuous DNA duplex running along a central channel.

Genomic Characterization of Cyanophage vB_AphaS-CL131 Infecting Filamentous Diazotrophic Cyanobacterium Aphanizomenon flos-aquae Reveals Novel Insights into Virus-Bacterium Interactions.

While filamentous cyanobacteria play a crucial role in food web dynamics and biogeochemical cycling of many aquatic ecosystems around the globe, the knowledge regarding the phages infecting them is limited. Here, we describe the complete genome of the virulent cyanophage vB_AphaS-CL131 (here, CL 131), a Siphoviridae phage that infects the filamentous diazotrophic bloom-forming cyanobacterium Aphanizomenon flos-aquae in the brackish Baltic Sea. CL 131 features a 112,793-bp double-stranded DNA (dsDNA) genome encompassing 149 putative open reading frames (ORFs), of which the majority (86%) lack sequence homology to genes with known functions in other bacteriophages or bacteria. Phylogenetic analysis revealed that CL 131 possibly represents a new evolutionary lineage within the group of cyanophages infecting filamentous cyanobacteria, which form a separate cluster from phages infecting unicellular cyanobacteria. CL 131 encodes a putative type V-U2 CRISPR-Cas system with one spacer (out of 10) targeting a DNA primase pseudogene in a cyanobacterium and a putative type II toxin-antitoxin system, consisting of a GNAT family N-acetyltransferase and a protein of unknown function containing the PRK09726 domain (characteristic of HipB antitoxins). Comparison of CL 131 proteins to reads from Baltic Sea and other available fresh- and brackish-water metagenomes and analysis of CRISPR-Cas arrays in publicly available A. flos-aquae genomes demonstrated that phages similar to CL 131 are present and dynamic in the Baltic Sea and share a common history with their hosts dating back at least several decades. In addition, different CRISPR-Cas systems within individual A. flos-aquae genomes targeted several sequences in the CL 131 genome, including genes related to virion structure and morphogenesis. Altogether, these findings revealed new genomic information for exploring viral diversity and provide a model system for investigation of virus-host interactions in filamentous cyanobacteria.IMPORTANCE The genomic characterization of novel cyanophage vB_AphaS-CL131 and the analysis of its genomic features in the context of other viruses, metagenomic data, and host CRISPR-Cas systems contribute toward a better understanding of aquatic viral diversity and distribution in general and of brackish-water cyanophages infecting filamentous diazotrophic cyanobacteria in the Baltic Sea in particular. The results of this study revealed previously undescribed features of cyanophage genomes (e.g., self-excising intein-containing putative dCTP deaminase and putative cyanophage-encoded CRISPR-Cas and toxin-antitoxin systems) and can therefore be used to predict potential interactions between bloom-forming cyanobacteria and their cyanophages.

A closely-related clade of globally distributed bloom-forming cyanobacteria within the Nostocales.

In order to better understand the relationships among current Nostocales cyanobacterial blooms, eight genomes were sequenced from cultured isolates or from environmental metagenomes of recent planktonic Nostocales blooms. Phylogenomic analysis of publicly available sequences placed the new genomes among a group of 15 genomes from four continents in a distinct ADA clade (Anabaena/Dolichospermum/Aphanizomenon) within the Nostocales. This clade contains four species-level groups, two of which include members with both Anabaena-like and Aphanizomenon flos-aquae-like morphology. The genomes contain many repetitive genetic elements and a sizable pangenome, in which ABC-type transporters are highly represented. Alongside common core genes for photosynthesis, the differentiation of N2-fixing heterocysts, and the uptake and incorporation of the major nutrients P, N and S, we identified several gene pathways in the pangenome that may contribute to niche partitioning. Genes for problematic secondary metabolites-cyanotoxins and taste-and-odor compounds-were sporadically present, as were other polyketide synthase (PKS) and nonribosomal peptide synthetase (NRPS) gene clusters. By contrast, genes predicted to encode the ribosomally generated bacteriocin peptides were found in all genomes.

DnaQ exonuclease-like domain of Cas2 promotes spacer integration in a type I-E CRISPR-Cas system.

CRISPR-Cas systems constitute an adaptive immune system that provides acquired resistance against phages and plasmids in prokaryotes. Upon invasion of foreign nucleic acids, some cells integrate short fragments of foreign DNA as spacers into the CRISPR locus to memorize the invaders and acquire resistance in the subsequent round of infection. This immunization step called adaptation is the least understood part of the CRISPR-Cas immunity. We have focused here on the adaptation stage of Streptococcus thermophilus DGCC7710 type I-E CRISPR4-Cas (St4) system. Cas1 and Cas2 proteins conserved in nearly all CRISPR-Cas systems are required for spacer acquisition. The St4 CRISPR-Cas system is unique because the Cas2 protein is fused to an additional DnaQ exonuclease domain. Here, we demonstrate that St4 Cas1 and Cas2-DnaQ form a multimeric complex, which is capable of integrating DNA duplexes with 3'-overhangs (protospacers) in vitro We further show that the DnaQ domain of Cas2 functions as a 3'-5'-exonuclease that processes 3'-overhangs of the protospacer to promote integration.

Harnessing the natural diversity and in vitro evolution of Cas9 to expand the genome editing toolbox.

In the past few years, the Cas9 endonuclease from the type II CRISPR-Cas bacterial antiviral defense system has revolutionized the genome editing field. Guided by an RNA molecule, Cas9 can be reprogrammed to target almost any DNA sequence: the only limitation being the short nucleotide sequence in the vicinity of the target, termed the PAM, which is characteristic for each Cas9 protein. Streptococcus pyogenes Cas9 which recognizes the NGG PAM is currently most widely used for genome manipulation. However, Cas9 orthologues and engineered Cas9 variants offer expanded genome targeting capabilities, improved specificity and biochemical properties.

Methods for decoding Cas9 protospacer adjacent motif (PAM) sequences: A brief overview.

Recently the Cas9, an RNA guided DNA endonuclease, emerged as a powerful tool for targeted genome manipulations. Cas9 protein can be reprogrammed to cleave, bind or nick any DNA target by simply changing crRNA sequence, however a short nucleotide sequence, termed PAM, is required to initiate crRNA hybridization to the DNA target. PAM sequence is recognized by Cas9 protein and must be determined experimentally for each Cas9 variant. Exploration of Cas9 orthologs could offer a diversity of PAM sequences and novel biochemical properties that may be beneficial for genome editing applications. Here we briefly review and compare Cas9 PAM identification assays that can be adopted for other PAM-dependent CRISPR-Cas systems.

Design of a CRISPR-Cas system to increase resistance of Bacillus subtilis to bacteriophage SPP1.

Clustered regularly interspaced short palindromic repeats (CRISPR) together with CRISPR-associated (cas) genes form an adaptive prokaryotic immune system which provides acquired resistance against viruses and plasmids. Bacillus subtilis presently is the best-characterized laboratory model for Gram-positive bacteria and also widely used for industrial production of enzymes, vitamins and antibiotics. In this study, we show that type II-A CRISPR-Cas system from Streptococcus thermophilus can be transferred into B. subtilis and provides heterologous protection against phage infection. We engineered a heterologous host by cloning S. thermophilus Cas9 and a spacer targeting bacteriophage SPP1 into the chromosome of B. subtilis, which does not harbor its own CRISPR-Cas systems. We found that the heterologous CRISPR-Cas system is functionally active in B. subtilis and provides resistance against bacteriophage SPP1 infection. The high efficiency of the acquired immunity against phage could be useful in generation of biotechnologically important B. subtilis strains with engineered chromosomes.

Lactase nonpersistence is directed by DNA-variation-dependent epigenetic aging.

The inability to digest lactose, due to lactase nonpersistence, is a common trait in adult mammals, except in certain human populations that exhibit lactase persistence. It is not known how the lactase gene is dramatically downregulated with age in most individuals but remains active in some individuals. We performed a comprehensive epigenetic study of human and mouse small intestines, by using chromosome-wide DNA-modification profiling and targeted bisulfite sequencing. Epigenetically controlled regulatory elements accounted for the differences in lactase mRNA levels among individuals, intestinal cell types and species. We confirmed the importance of these regulatory elements in modulating lactase mRNA levels by using CRISPR-Cas9-induced deletions. Genetic factors contribute to epigenetic changes occurring with age at the regulatory elements, because lactase-persistence and lactase-nonpersistence DNA haplotypes demonstrated markedly different epigenetic aging. Thus, genetic factors enable a gradual accumulation of epigenetic changes with age, thereby influencing phenotypic outcome.

Rewiring Cas9 to Target New PAM Sequences.

Papers by Anders et al. (2016) and Hirano et al. (2016b), published in this issue of Molecular Cell, show that SpCas9 uses an induced fit mechanism to recognize altered protospacer adjacent motif (PAM) sequences.

Streptococcus thermophilus CRISPR-Cas9 Systems Enable Specific Editing of the Human Genome.

RNA-guided nucleases (RGNs) based on the type II CRISPR-Cas9 system of Streptococcus pyogenes (Sp) have been widely used for genome editing in experimental models. However, the nontrivial level of off-target activity reported in several human cells may hamper clinical translation. RGN specificity depends on both the guide RNA (gRNA) and the protospacer adjacent motif (PAM) recognized by the Cas9 protein. We hypothesized that more stringent PAM requirements reduce the occurrence of off-target mutagenesis. To test this postulation, we generated RGNs based on two Streptococcus thermophilus (St) Cas9 proteins, which recognize longer PAMs, and performed a side-by-side comparison of the three RGN systems targeted to matching sites in two endogenous human loci, PRKDC and CARD11. Our results demonstrate that in samples with comparable on-target cleavage activities, significantly lower off-target mutagenesis was detected using St-based RGNs as compared to the standard Sp-RGNs. Moreover, similarly to SpCas9, the StCas9 proteins accepted truncated gRNAs, suggesting that the specificities of St-based RGNs can be further improved. In conclusion, our results show that Cas9 proteins with longer or more restrictive PAM requirements provide a safe alternative to SpCas9-based RGNs and hence a valuable option for future human gene therapy applications.

Rapid characterization of CRISPR-Cas9 protospacer adjacent motif sequence elements.

To expand the repertoire of Cas9s available for genome targeting, we present a new in vitro method for the simultaneous examination of guide RNA and protospacer adjacent motif (PAM) requirements. The method relies on the in vitro cleavage of plasmid libraries containing a randomized PAM as a function of Cas9-guide RNA complex concentration. Using this method, we accurately reproduce the canonical PAM preferences for Streptococcus pyogenes, Streptococcus thermophilus CRISPR3 (Sth3), and CRISPR1 (Sth1). Additionally, PAM and sgRNA solutions for a novel Cas9 protein from Brevibacillus laterosporus are provided by the assay and are demonstrated to support functional activity in vitro and in plants.

Cas3 nuclease-helicase activity assays.

Cas3 is a signature protein of the type I CRISPR-Cas systems and typically contains HD phosphohydrolase and Superfamily 2 (SF2) helicase domains. In the type I CRISPR-Cas systems Cas3 functions as a slicer that provides foreign DNA degradation. Biochemical analysis indicate that Cas3 of the Streptococcus thermophilus DGCC7710 (St-Cas3) CRISPR4 system is a single-stranded DNA nuclease which also possesses a single-stranded DNA-stimulated ATPase activity, which is coupled to unwinding of DNA/DNA and RNA/DNA duplexes in 3' to 5' direction. The interplay between the nuclease and ATPase/helicase activities of St-Cas3 results in DNA degradation. Here, we describe assays for monitoring of St-Cas3 nuclease, ATPase and helicase activities in a stand-alone form and in the presence of the Cascade ribonucleoprotein complex. These assays can be easily adapted for biochemical analysis of Cas3 proteins from different microorganisms.