Machine learning-coupled combinatorial mutagenesis enables resource-efficient engineering of CRISPR-Cas9 genome editor activities.

The genome-editing Cas9 protein uses multiple amino-acid residues to bind the target DNA. Considering only the residues in proximity to the target DNA as potential sites to optimise Cas9's activity, the number of combinatorial variants to screen through is too massive for a wet-lab experiment. Here we generate and cross-validate ten in silico and experimental datasets of multi-domain combinatorial mutagenesis libraries for Cas9 engineering, and demonstrate that a machine learning-coupled engineering approach reduces the experimental screening burden by as high as 95% while enriching top-performing variants by ∼7.5-fold in comparison to the null model. Using this approach and followed by structure-guided engineering, we identify the N888R/A889Q variant conferring increased editing activity on the protospacer adjacent motif-relaxed KKH variant of Cas9 nuclease from Staphylococcus aureus (KKH-SaCas9) and its derived base editor in human cells. Our work validates a readily applicable workflow to enable resource-efficient high-throughput engineering of genome editor's activity.

High-fidelity KKH variant of Staphylococcus aureus Cas9 nucleases with improved base mismatch discrimination.

The Cas9 nuclease from Staphylococcus aureus (SaCas9) holds great potential for use in gene therapy, and variants with increased fidelity have been engineered. However, we find that existing variants have not reached the greatest accuracy to discriminate base mismatches and exhibited much reduced activity when their mutations were grafted onto the KKH mutant of SaCas9 for editing an expanded set of DNA targets. We performed structure-guided combinatorial mutagenesis to re-engineer KKH-SaCas9 with enhanced accuracy. We uncover that introducing a Y239H mutation on KKH-SaCas9's REC domain substantially reduces off-target edits while retaining high on-target activity when added to a set of mutations on REC and RuvC domains that lessen its interactions with the target DNA strand. The Y239H mutation is modelled to have removed an interaction from the REC domain with the guide RNA backbone in the guide RNA-DNA heteroduplex structure. We further confirmed the greatly improved genome-wide editing accuracy and single-base mismatch discrimination of our engineered variants, named KKH-SaCas9-SAV1 and SAV2, in human cells. In addition to generating broadly useful KKH-SaCas9 variants with unprecedented accuracy, our findings demonstrate the feasibility for multi-domain combinatorial mutagenesis on SaCas9's DNA- and guide RNA- interacting residues to optimize its editing fidelity.

Deep RNA Sequencing Revealed Fusion Junctional Heterogeneity May Predict Crizotinib Treatment Efficacy in ALK-Rearranged NSCLC.

INTRODUCTION: Gene fusion variants in ALK-rearranged NSCLC may predict patient outcomes, but previous results have been inconclusive. Fusion isoforms coexisting in the same tumor may affect the efficacy of targeted therapy, but they have not been investigated. METHODS: Patients with ALK-rearranged NSCLC who received crizotinib treatments were recruited. Precrizotinib tumor tissues were analyzed by the anchored multiplex polymerase chain reaction for targeted RNA sequencing. Kaplan-Meier and Cox regression were used to compare overall and progression-free survivals. RESULTS: Of the 51 studied subjects, EML4-ALK variant types v1, v2, v3, and others were detected in 23 (45.1%), five (9.8%), 19 (37.3%), and four patients (7.8%), respectively. Multiple EML4-ALK RNA isoforms were detected in 24 tumors (47.1%), and single isoform in 27 (52.9%). Most of the v3 tumors (16 of 19) harbored both v3a and v3b RNA isoforms. Multiple isoforms were also detected in eight non-v3 tumors (33.3% of all 24 multiple isoforms; five v1, two v5', and one v2). Compared with patients with single isoform, those with multiple isoforms had worse progression-free (hazard ratio and 95% confidence interval: 2.45 [1.06-5.69]) and overall (hazard ratio [95% confidence interval]: 3.74 [1.26-11.13]) survivals after adjusting for potential confounders including variant type. Using the patient-derived H2228 cells known to express v3a and v3b, our single-cell polymerase chain reaction detected either v3a or v3b in most single cells. Treatment of H2228 cells by three ALK inhibitors revealed increased ratios of v3a-to-v3b expression over time. CONCLUSIONS: Intratumoral EML4-ALK isoforms may predict the efficacy of targeted therapy in ALK-rearranged NSCLC. Temporal changes of intratumoral fusion isoforms may result from differential selection pressures that a drug might have on one isoform over another. Larger studies on fusion heterogeneity using RNA sequencing are warranted.

A Three-Way Combinatorial CRISPR Screen for Analyzing Interactions among Druggable Targets.

We present a CRISPR-based multi-gene knockout screening system and toolkits for extensible assembly of barcoded high-order combinatorial guide RNA libraries en masse. We apply this system for systematically identifying not only pairwise but also three-way synergistic therapeutic target combinations and successfully validate double- and triple-combination regimens for suppression of cancer cell growth and protection against Parkinson's disease-associated toxicity. This system overcomes the practical challenges of experimenting on a large number of high-order genetic and drug combinations and can be applied to uncover the rare synergistic interactions between druggable targets.