The Transpeptidase Sortase A Binds Nucleic Acids and Mediates Mammalian Cell Labeling.

Wild-type sortase A is an important virulence factor displaying a diverse array of proteins on the surface of bacteria. This protein display relies on the transpeptidase activity of sortase A, which is widely engineered to allow protein ligation and protein engineering based on the interaction between sortase A and peptides. Here an unknown interaction is found between sortase A from Staphylococcus aureus and nucleic acids, in which exogenously expressed engineered sortase A binds oligonucleotides in vitro and is independent of its canonical transpeptidase activity. When incubated with mammalian cells, engineered sortase A further mediates oligonucleotide labeling to the cell surface, where sortase A attaches itself and is part of the labeled moiety. The labeling reaction can also be mediated by many classes of wild-type sortases as well. Cell surface GAG appears involved in sortase-mediated oligonucleotide cell labeling, as demonstrated by CRISPR screening. This interaction property is utilized to develop a technique called CellID to facilitate sample multiplexing for scRNA-seq and shows the potential of using sortases to label cells with diverse oligonucleotides. Together, the binding between sortase A and nucleic acids opens a new avenue to understanding the virulence of wild-type sortases and exploring the application of sortases in biotechnology.

Machine learning-based prediction models to guide the selection of Cas9 variants for efficient gene editing.

The increasing emergence of Cas9 variants has attracted broad interest, as these variants were designed to expand CRISPR applications. New Cas9 variants typically feature higher editing efficiency, improved editing specificity, or alternative PAM sequences. To select Cas9 variants and gRNAs for high-fidelity and efficient genome editing, it is crucial to systematically quantify the editing performances of gRNAs and develop prediction models based on high-quality datasets. Using synthetic gRNA-target paired libraries and next-generation sequencing, we compared the activity and specificity of gRNAs of four SpCas9 variants. The nucleotide composition in the PAM-distal region had more influence on the editing efficiency of HiFi Cas9 and LZ3 Cas9. We further developed machine learning models to predict the gRNA efficiency and specificity for the four Cas9 variants. To aid users from broad research areas, the machine learning models for the predictions of gRNA editing efficiency within human genome sites are available on our website.

A top variant identification pipeline for protein engineering.

Understanding the fitness of protein variants with combinatorial mutations is critical for effective protein engineering. In this issue of Cell Systems, Chu et al. present TopVIP, a top variant identification pipeline that enables accurate picking of the greatest number of best-performing protein variants with high-fitness leveraging zero-shot predictor and low-N iterative sampling.

The RNA m6A demethylase ALKBH5 drives emergency granulopoiesis and neutrophil mobilization by upregulating G-CSFR expression.

Emergency granulopoiesis and neutrophil mobilization that can be triggered by granulocyte colony-stimulating factor (G-CSF) through its receptor G-CSFR are essential for antibacterial innate defense. However, the epigenetic modifiers crucial for intrinsically regulating G-CSFR expression and the antibacterial response of neutrophils remain largely unclear. N6-methyladenosine (m6A) RNA modification and the related demethylase alkB homolog 5 (ALKBH5) are key epigenetic regulators of immunity and inflammation, but their roles in neutrophil production and mobilization are still unknown. We used cecal ligation and puncture (CLP)-induced polymicrobial sepsis to model systemic bacterial infection, and we report that ALKBH5 is required for emergency granulopoiesis and neutrophil mobilization. ALKBH5 depletion significantly impaired the production of immature neutrophils in the bone marrow of septic mice. In addition, Alkbh5-deficient septic mice exhibited higher retention of mature neutrophils in the bone marrow and defective neutrophil release into the circulation, which led to fewer neutrophils at the infection site than in their wild-type littermates. During bacterial infection, ALKBH5 imprinted production- and mobilization-promoting transcriptome signatures in both mouse and human neutrophils. Mechanistically, ALKBH5 erased m6A methylation on the CSF3R mRNA to increase the mRNA stability and protein expression of G-CSFR, consequently upregulating cell surface G-CSFR expression and downstream STAT3 signaling in neutrophils. The RIP-qPCR results confirmed the direct binding of ALKBH5 to the CSF3R mRNA, and the binding strength declined upon bacterial infection, accounting for the decrease in G-CSFR expression on bacteria-infected neutrophils. Considering these results collectively, we define a new role of ALKBH5 in intrinsically driving neutrophil production and mobilization through m6A demethylation-dependent posttranscriptional regulation, indicating that m6A RNA modification in neutrophils is a potential target for treating bacterial infections and neutropenia.

Optimized minimal genome-wide human sgRNA library.

Genome-wide clustered regularly interspaced short palindromic repeats (CRISPR)-based knockout screening is revolting the genetic analysis of a cellular or molecular phenotype in question but is challenged by the large size of single-guide RNA (sgRNA) library. Here we designed a minimal genome-wide human sgRNA library, H-mLib, which is composed of 21,159 sgRNA pairs assembled based on a dedicated selection strategy from all potential SpCas9/sgRNAs in the human genome. These sgRNA pairs were cloned into a dual-gRNA vector each targeting one gene, resulting in a compact library size nearly identical to the number of human protein-coding genes. The performance of the H-mLib was benchmarked to other CRISPR libraries in a proliferation screening conducted in K562 cells. We also identified groups of core essential genes and cell-type specific essential genes by comparing the screening results from the K562 and Jurkat cells. Together, the H-mLib exemplified high specificity and sensitivity in identifying essential genes while containing minimal library complexity, emphasizing its advantages and applications in CRISPR screening with limited cell numbers.

Deep sampling of gRNA in the human genome and deep-learning-informed prediction of gRNA activities.

Life science studies involving clustered regularly interspaced short palindromic repeat (CRISPR) editing generally apply the best-performing guide RNA (gRNA) for a gene of interest. Computational models are combined with massive experimental quantification on synthetic gRNA-target libraries to accurately predict gRNA activity and mutational patterns. However, the measurements are inconsistent between studies due to differences in the designs of the gRNA-target pair constructs, and there has not yet been an integrated investigation that concurrently focuses on multiple facets of gRNA capacity. In this study, we analyzed the DNA double-strand break (DSB)-induced repair outcomes and measured SpCas9/gRNA activities at both matched and mismatched locations using 926,476 gRNAs covering 19,111 protein-coding genes and 20,268 non-coding genes. We developed machine learning models to forecast the on-target cleavage efficiency (AIdit_ON), off-target cleavage specificity (AIdit_OFF), and mutational profiles (AIdit_DSB) of SpCas9/gRNA from a uniformly collected and processed dataset by deep sampling and massively quantifying gRNA capabilities in K562 cells. Each of these models exhibited superlative performance in predicting SpCas9/gRNA activities on independent datasets when benchmarked with previous models. A previous unknown parameter was also empirically determined regarding the "sweet spot" in the size of datasets used to establish an effective model to predict gRNA capabilities at a manageable experimental scale. In addition, we observed cell type-specific mutational profiles and were able to link nucleotidylexotransferase as the key factor driving these outcomes. These massive datasets and deep learning algorithms have been implemented into the user-friendly web service http://crispr-aidit.com to evaluate and rank gRNAs for life science studies.

PEAC-seq adopts Prime Editor to detect CRISPR off-target and DNA translocation.

CRISPR technology holds significant promise for biological studies and gene therapies because of its high flexibility and efficiency when applied in mammalian cells. But endonuclease (e.g., Cas9) potentially generates undesired edits; thus, there is an urgent need to comprehensively identify off-target sites so that the genotoxicities can be accurately assessed. To date, it is still challenging to streamline the entire process to specifically label and efficiently enrich the cleavage sites from unknown genomic locations. Here we develop PEAC-seq, in which we adopt the Prime Editor to insert a sequence-optimized tag to the editing sites and enrich the tagged regions with site-specific primers for high throughput sequencing. Moreover, we demonstrate that PEAC-seq could identify DNA translocations, which are more genotoxic but usually overlooked by other off-target detection methods. As PEAC-seq does not rely on exogenous oligodeoxynucleotides to label the editing site, we also conduct in vivo off-target identification as proof of concept. In summary, PEAC-seq provides a comprehensive and streamlined strategy to identify CRISPR off-targeting sites in vitro and in vivo, as well as DNA translocation events. This technique further diversified the toolkit to evaluate the genotoxicity of CRISPR applications in research and clinics.

An in-library ligation strategy and its application in CRISPR/Cas9 screening of high-order gRNA combinations.

Simultaneous targeting multiple genes is a big advantage of CRISPR (clustered regularly interspaced short palindromic repeats) genome editing but challenging to achieve in CRISPR screening. The crosstalk among genes or gene products is a common and fundamental mechanism to ensure cellular stability and functional diversity. However, the screening approach to map high-order gene combinations to the interesting phenotype is still lacking. Here, we developed a universal in-library ligation strategy and applied it to generate multiplexed CRISPR library, which could perturb four pre-designed targets in a cell. We conducted in vivo CRISPR screening for potential guide RNA (gRNA) combinations inducing anti-tumor immune responses. Simultaneously disturbing a combination of three checkpoints in CD8+ T cells was demonstrated to be more effective than disturbing Pdcd1 only for T cell activation in the tumor environment. This study developed a novel in-library ligation strategy to facilitate the multiplexed CRISPR screening, which could extend our ability to explore the combinatorial outcomes from coordinated gene behaviors.

m6A demethylase ALKBH5 is required for antibacterial innate defense by intrinsic motivation of neutrophil migration.

Neutrophil migration into the site of infection is necessary for antibacterial innate defense, whereas impaired neutrophil migration may result in excessive inflammation and even sepsis. The neutrophil migration directed by extracellular signals such as chemokines has been extensively studied, yet the intrinsic mechanism for determining neutrophil ability to migrate needs further investigation. N6-methyladenosine (m6A) RNA modification is important in immunity and inflammation, and our preliminary data indicate downregulation of RNA m6A demethylase alkB homolog 5 (ALKBH5) in neutrophils during bacterial infection. Whether m6A modification and ALKBH5 might intrinsically modulate neutrophil innate response remain unknown. Here we report that ALKBH5 is required for antibacterial innate defense by enhancing intrinsic ability of neutrophil migration. We found that deficiency of ALKBH5 increased mortality of mice with polymicrobial sepsis induced by cecal ligation and puncture (CLP), and Alkbh5-deficient CLP mice exhibited higher bacterial burden and massive proinflammatory cytokine production in the peritoneal cavity and blood because of less neutrophil migration. Alkbh5-deficient neutrophils had lower CXCR2 expression, thus exhibiting impaired migration toward chemokine CXCL2. Mechanistically, ALKBH5-mediated m6A demethylation empowered neutrophils with high migration capability through altering the RNA decay, consequently regulating protein expression of its targets, neutrophil migration-related molecules, including increased expression of neutrophil migration-promoting CXCR2 and NLRP12, but decreased expression of neutrophil migration-suppressive PTGER4, TNC, and WNK1. Our findings reveal a previously unknown role of ALKBH5 in imprinting migration-promoting transcriptome signatures in neutrophils and intrinsically promoting neutrophil migration for antibacterial defense, highlighting the potential application of targeting neutrophil m6A modification in controlling bacterial infections.

New Chromatin Run-On Reaction Enables Global Mapping of Active RNA Polymerase Locations in an Enrichment-free Manner.

The development of a simple and cost-effective method to map the distribution of RNA polymerase II (RNPII) genome-wide at a high resolution is highly beneficial to study cellular transcriptional activity. Here we report a mutation-based and enrichment-free global chromatin run-on sequencing (mGRO-seq) technique to locate active RNPII sites genome-wide at near-base resolution. An adenosine triphosphate (ATP) analog named N6-allyladenosine triphosphate (a6ATP) was designed and could be incorporated into nascent RNAs at RNPII-located positions during a chromatin run-on reaction. By treatment of the run-on RNAs with a mild iodination reaction and subjection of the products to reverse transcription into complementary DNA (cDNA), base mismatch occurs at the original a6A incorporation sites, thus making the RNPII locations detected in the high-throughput cDNA sequencing. The mGRO-seq yields both the map of RNPII sites and the chromatin RNA abundance and holds great promise for the study of single-cell transcriptional activity.

Correction to: YTHDF2 reduction fuels inflammation and vascular abnormalization in hepatocellular carcinoma.

An amendment to this paper has been published and can be accessed via the original article.

Direct-seq: programmed gRNA scaffold for streamlined scRNA-seq in CRISPR screen.

CRISPR-based genome perturbation provides a new avenue to conveniently change DNA sequences, transcription, and epigenetic modifications in genetic screens. However, it remains challenging to assay the complex molecular readouts after perturbation at high resolution and at scale. By introducing an A/G mixed capture sequence into the gRNA scaffold, we demonstrate that gRNA transcripts could be directly reverse transcribed by poly (dT) primer together with the endogenous mRNA, followed by high-content molecular phenotyping in scRNA-seq (Direct-seq). With this method, the CRISPR perturbation and its transcriptional readouts can be profiled together in a streamlined workflow.

A metabolic labeling method detects m6A transcriptome-wide at single base resolution.

Transcriptome-wide mapping of N6-methyladenosine (m6A) at base resolution remains an issue, impeding our understanding of m6A roles at the nucleotide level. Here, we report a metabolic labeling method to detect mRNA m6A transcriptome-wide at base resolution, called 'm6A-label-seq'. Human and mouse cells could be fed with a methionine analog, Se-allyl-L-selenohomocysteine, which substitutes the methyl group on the enzyme cofactor SAM with the allyl. Cellular RNAs could therefore be metabolically modified with N6-allyladenosine (a6A) at supposed m6A-generating adenosine sites. We pinpointed the mRNA a6A locations based on iodination-induced misincorporation at the opposite site in complementary DNA during reverse transcription. We identified a few thousand mRNA m6A sites in human HeLa, HEK293T and mouse H2.35 cells, carried out a parallel comparison of m6A-label-seq with available m6A sequencing methods, and validated selected sites by an orthogonal method. This method offers advantages in detecting clustered m6A sites and holds promise to locate nuclear nascent RNA m6A modifications.

The RNA-binding protein FMRP facilitates the nuclear export of N6-methyladenosine-containing mRNAs.

N6-Methyladenosine (m6A) is the most abundant post-transcriptional mRNA modification in eukaryotes and exerts many of its effects on gene expression through reader proteins that bind specifically to m6A-containing transcripts. Fragile X mental retardation protein (FMRP), an RNA-binding protein, has previously been shown to affect the translation of target mRNAs and trafficking of mRNA granules. Loss of function of FMRP causes fragile X syndrome, the most common form of inherited intellectual disability in humans. Using HEK293T cells, siRNA-mediated gene knockdown, cytoplasmic and nuclear fractions, RNA-Seq, and LC-MS/MS analyses, we demonstrate here that FMRP binds directly to a collection of m6A sites on mRNAs. FMRP depletion increased mRNA m6A levels in the nucleus. Moreover, the abundance of FMRP targets in the cytoplasm relative to the nucleus was decreased in Fmr1-KO mice, an effect also observed in highly methylated genes. We conclude that FMRP may affect the nuclear export of m6A-modified RNA targets.

YTHDF2 reduction fuels inflammation and vascular abnormalization in hepatocellular carcinoma.

BACKGROUND: Dynamic N6-methyladenosine (m6A) modification was previously identified as a ubiquitous post-transcriptional regulation that affected mRNA homeostasis. However, the m6A-related epitranscriptomic alterations and functions remain elusive in human cancer. Here we aim to identify the profile and outcome of m6A-methylation in hepatocellular carcinoma (HCC). RESULTS: Using liquid chromatography-tandem mass spectrometry and m6A-immunoprecipitation in combination with high-throughput sequencing, we determined the m6A-mRNA levels in human HCC. Human HCC exhibited a characteristic gain of m6A modification in tandem with an increase of mRNA expression, owing to YTH domain family 2 (YTHDF2) reduction. The latter predicted poor classification and prognosis of HCC patients, and highly correlated with HCC m6A landscape. YTHDF2 silenced in human HCC cells or ablated in mouse hepatocytes provoked inflammation, vascular reconstruction and metastatic progression. Mechanistically, YTHDF2 processed the decay of m6A-containing interleukin 11 (IL11) and serpin family E member 2 (SERPINE2) mRNAs, which were responsible for the inflammation-mediated malignancy and disruption of vascular normalization. Reciprocally, YTHDF2 transcription succumbed to hypoxia-inducible factor-2α (HIF-2α). Administration of a HIF-2α antagonist (PT2385) restored YTHDF2-programed epigenetic machinery and repressed liver cancer. CONCLUSION: Our results have characterized the m6A-mRNA landscape in human HCC and revealed YTHDF2 as a molecular 'rheostat' in epitranscriptome and cancer progression.

Viral N6-methyladenosine upregulates replication and pathogenesis of human respiratory syncytial virus.

N6-methyladenosine (m6A) is the most prevalent internal modification of mRNAs in most eukaryotes. Here we show that RNAs of human respiratory syncytial virus (RSV) are modified by m6A within discreet regions and that these modifications enhance viral replication and pathogenesis. Knockdown of m6A methyltransferases decreases RSV replication and gene expression whereas knockdown of m6A demethylases has the opposite effect. The G gene transcript contains the most m6A modifications. Recombinant RSV variants expressing G transcripts that lack particular clusters of m6A display reduced replication in A549 cells, primary well differentiated human airway epithelial cultures, and respiratory tracts of cotton rats. One of the m6A-deficient variants is highly attenuated yet retains high immunogenicity in cotton rats. Collectively, our results demonstrate that viral m6A methylation upregulates RSV replication and pathogenesis and identify viral m6A methylation as a target for rational design of live attenuated vaccine candidates for RSV and perhaps other pneumoviruses.

N 6-methyladenosine RNA modification-mediated cellular metabolism rewiring inhibits viral replication.

Host cell metabolism can be modulated by viral infection, affecting viral survival or clearance. Yet the cellular metabolism rewiring mediated by the N 6-methyladenosine (m6A) modification in interactions between virus and host remains largely unknown. Here we report that in response to viral infection, host cells impair the enzymatic activity of the RNA m6A demethylase ALKBH5. This behavior increases the m6A methylation on α-ketoglutarate dehydrogenase (OGDH) messenger RNA (mRNA) to reduce its mRNA stability and protein expression. Reduced OGDH decreases the production of the metabolite itaconate that is required for viral replication. With reduced OGDH and itaconate production in vivo, Alkbh5-deficient mice display innate immune response-independent resistance to viral exposure. Our findings reveal that m6A RNA modification-mediated down-regulation of the OGDH-itaconate pathway reprograms cellular metabolism to inhibit viral replication, proposing potential targets for controlling viral infection.

m6A mRNA demethylase FTO regulates melanoma tumorigenicity and response to anti-PD-1 blockade.

Melanoma is one of the most deadly and therapy-resistant cancers. Here we show that N6-methyladenosine (m6A) mRNA demethylation by fat mass and obesity-associated protein (FTO) increases melanoma growth and decreases response to anti-PD-1 blockade immunotherapy. FTO level is increased in human melanoma and enhances melanoma tumorigenesis in mice. FTO is induced by metabolic starvation stress through the autophagy and NF-κB pathway. Knockdown of FTO increases m6A methylation in the critical protumorigenic melanoma cell-intrinsic genes including PD-1 (PDCD1), CXCR4, and SOX10, leading to increased RNA decay through the m6A reader YTHDF2. Knockdown of FTO sensitizes melanoma cells to interferon gamma (IFNγ) and sensitizes melanoma to anti-PD-1 treatment in mice, depending on adaptive immunity. Our findings demonstrate a crucial role of FTO as an m6A demethylase in promoting melanoma tumorigenesis and anti-PD-1 resistance, and suggest that the combination of FTO inhibition with anti-PD-1 blockade may reduce the resistance to immunotherapy in melanoma.

Histone H3 trimethylation at lysine 36 guides m6A RNA modification co-transcriptionally.

DNA and histone modifications have notable effects on gene expression1. Being the most prevalent internal modification in mRNA, the N6-methyladenosine (m6A) mRNA modification is as an important post-transcriptional mechanism of gene regulation2-4 and has crucial roles in various normal and pathological processes5-12. However, it is unclear how m6A is specifically and dynamically deposited in the transcriptome. Here we report that histone H3 trimethylation at Lys36 (H3K36me3), a marker for transcription elongation, guides m6A deposition globally. We show that m6A modifications are enriched in the vicinity of H3K36me3 peaks, and are reduced globally when cellular H3K36me3 is depleted. Mechanistically, H3K36me3 is recognized and bound directly by METTL14, a crucial component of the m6A methyltransferase complex (MTC), which in turn facilitates the binding of the m6A MTC to adjacent RNA polymerase II, thereby delivering the m6A MTC to actively transcribed nascent RNAs to deposit m6A co-transcriptionally. In mouse embryonic stem cells, phenocopying METTL14 knockdown, H3K36me3 depletion also markedly reduces m6A abundance transcriptome-wide and in pluripotency transcripts, resulting in increased cell stemness. Collectively, our studies reveal the important roles of H3K36me3 and METTL14 in determining specific and dynamic deposition of m6A in mRNA, and uncover another layer of gene expression regulation that involves crosstalk between histone modification and RNA methylation.