Engineering of a high-fidelity Cas12a nuclease variant capable of allele-specific editing.

CRISPR-Cas12a, often regarded as a precise genome editor, still requires improvements in specificity. In this study, we used a GFP-activation assay to screen 14 new Cas12a nucleases for mammalian genome editing, successfully identifying 9 active ones. Notably, these Cas12a nucleases prefer pyrimidine-rich PAMs. Among these nucleases, we extensively characterized Mb4Cas12a obtained from Moraxella bovis CCUG 2133, which recognizes a YYN PAM (Y = C or T). Our biochemical analysis demonstrates that Mb4Cas12a can cleave double-strand DNA across a wide temperature range. To improve specificity, we constructed a SWISS-MODEL of Mb4Cas12a based on the FnCas12a crystal structure and identified 8 amino acids potentially forming hydrogen bonds at the target DNA-crRNA interface. By replacing these amino acids with alanine to disrupt the hydrogen bond, we tested the influence of each mutation on Mb4Cas12a specificity. Interestingly, the F370A mutation improved specificity with minimal influence on activity. Further study showed that Mb4Cas12a-F370A is capable of discriminating single-nucleotide polymorphisms. These new Cas12a orthologs and high-fidelity variants hold substantial promise for therapeutic applications.

Molecular evolutionary process of advanced gastric cancer during sequential chemotherapy detected by circulating tumor DNA.

BACKGROUND: Efficacy of conventional sequential chemotherapy paradigm for advanced gastric cancer (AGC) patients has largely plateaued. Dynamic molecular changes during and after sequential chemotherapy have not been fully delineated. We aimed to profile the molecular evolutionary process of AGC patients during sequential chemotherapy by next generation sequencing (NGS) of plasma circulating tumor DNA (ctDNA). METHODS: A total of 30 chemo-naïve patients who were diagnosed with unresectable advanced or metastatic stomach adenocarcinoma were enrolled. All patients received sequential chemotherapy regimens following the clinical guideline. One hundred and eight serial peripheral blood samples were collected at baseline, radiographical assessment and disease progression. Plasma ctDNA was isolated and a customized NGS panel was used to detect the genomic features of ctDNA including single nucleotide variants (SNVs) and gene-level copy number variations (CNVs). KEGG pathway enrichment analysis was performed. RESULTS: Platinum-based combination chemotherapy was administrated as first-line regimen. Objective response rate was 50% (15/30). Patients with higher baseline values of copy number instability (CNI), CNVs and variant allel frequency (VAF) were more sensitive to platinum-based first-line regimens. Tumor mutation burden (TMB), CNI and CNV burden at partial response and stable disease were significantly lower than those at baseline, where at progressive disease they recovered to baseline levels. Dynamic change of TMB (ΔTMB) was correlated with progression-free survival of first-line treatment. Fluctuating changes of SNVs and gene-level CNVs could be observed during sequential chemotherapy. Under the pressure of conventional chemotherapy, the number of novel gene-level CNVs were found to be higher than that of novel SNVs. Such novel molecular alterations could be enriched into multiple common oncologic signaling pathways, including EGFR tyrosine kinase inhibitor resistance and platinum drug resistance pathways, where their distributions were found to be highly heterogenous among patients. The impact of subsequent regimens, including paclitaxel-based and irinotecan-based regimens, on the molecular changes driven by first-line therapy was subtle. CONCLUSION: Baseline and dynamic changes of genomic features of ctDNA could be biomarkers for predicting response of platinum-based first-line chemotherapy in AGC patients. After treatment with standard chemotherapy regimens, convergent oncologic pathway enrichment was identified, which is yet characterized by inter-patient heterogenous gene-level CNVs.

Mutational landscape of DNA damage response deficiency-related genes and its association with immune biomarkers in esophageal squamous cell carcinoma.

Esophageal squamous cell carcinoma (ESCC) has limited effective treatment strategies. DNA damage response (DDR) genes are of therapeutic interest in multiple cancer types. This study aimed to depict the landscape of DDR mutations in ESCC and evaluate the association between DDR mutations and known immunotherapy biomarkers. We recruited 250 Chinese patients with ESCC and performed next-generation sequencing. A total of 107 patients underwent a PD-L1 examination. Among the 250 patients, 73 (29.2%) harbored at least one DDR gene mutation and were defined as DDR-mut. Among the six functional DDR pathways, homologous recombination (HR) accounted for 12.4% (31/250). DDR-mut patients were significantly associated with higher tumor mutational burden than those in the DDR-wt group (p=7.4e-07). Patients with PDL1-H accounted for 21.2% (36/107) of the patients. PDL1-H was more prevalent in DDR-mut than DDR-wt, although the p-value did not reach a significant level (40.5% vs. 30%, p=0.29). Further analysis revealed that BRCA1, one of the most frequently mutated genes in the HR pathway, was significantly associated with PDL1-H (p=0.01). Our data revealed a subset of patients with ESCC harbored DDR gene mutations. Patients with these DDR gene mutations are significantly associated with immune biomarkers, implying the potential feasibility of combining DDR agents with immunotherapy in patients with DDR deficiency.

Germline ALK variations are associated with a poor prognosis in glioma and IDH-wildtype glioblastoma.

BACKGROUND: Glioma is the most common primary brain tumor. Clear classification is crucial for accurate diagnosis and individualized treatment. Histopathological characteristics and genetic alterations have shown to be related to prognosis and treatment response. Germline variants are important components of genetic alterations. However, the distribution of germline variations in glioma patients and their association with survival remain unknown. METHODS: We carried out whole-exome sequencing on 99 cases to explore germline variants in glioma. We also analyzed the association of germline variants with clinicopathological features and other prognostic indicators. RESULTS: All the glioma cases harbored rare germline variants. Germline ALK variants (gALK-Mut) were identified in 12/99 (12.12%) patients. The gALK-Mut patients had significantly shorter overall survival than germline ALK wildtype (gALK-WT) patients in the all glioma group (99 cases) and the subset of patients with IDH-wildtype glioblastoma (IDH-WT-GBM, 39 cases) (P = 0.013 and 0.027, respectively). The gALK-Mut patients also had higher frequency of BIRC5, PIK3CA and RPN1 somatic mutations than the gALK-WT patients in IDH-WT-GBM. Other confounding factors appeared to contribute to patient survival. The subgroup of patients in IDH-WT-GBM with gALK-Mut/TP53-Mut had worse prognosis than the gALK-WT/TP53-Mut subgroup (P = 0.031); The gALK-Mut/TERT-WT and gALK-Mut/TERT-Mut subgroups both had a worse prognosis than the gALK-WT/TERT-Mut subgroup (P = 0.031 and 0.018, respectively). CONCLUSIONS: Our study revealed ALK variation was an independent indicator of poor prognosis in glioma and IDH-WT-GBM. It could be a promising biomarker and tractable therapeutic target for this deadly disease.

Patient specific circulating tumor DNA fingerprints to monitor treatment response across multiple tumors.

BACKGROUND: Circulating tumor DNA (ctDNA) offers a convenient way to monitor tumor progression and treatment response. Because tumor mutational profiles are highly variable from person to person, a fixed content panel may be insufficient to track treatment response in all patients. METHODS: We design ctDNA fingerprint panels specific to individual patients which are based on whole exome sequencing and target to high frequency clonal population clusters in patients. We test the fingerprint panels in 313 patients who together have eight tumor types (colorectal, hepatocellular, gastric, breast, pancreatic, and esophageal carcinomas and lung cancer and cholangiocarcinoma) and exposed to multiple treatment methods (surgery, chemotherapy, radiotherapy, targeted-drug therapy, immunotherapy, and combinations of them). We also monitor drug-related mutations in the patients using a pre-designed panel with eight hotspot genes. RESULTS: 291 (93.0%) designed fingerprint panels harbor less than ten previously known tumor genes. We detected 7475 ctDNA mutations in 238 (76%) patients and 6196 (96.0%) of the mutations are detected in only one test. Both the level of ctDNA content fraction (CCF) and fold change of CCF (between the definitive and proceeding tests) are highly correlated with clinical outcomes (p-values 1.36e-6 for level and 5.64e-10 for fold change, Kruskal-Wallis test). The CCFs of PD patients are an order of magnitude higher than the CCFs of SD and OR patients (median/mean 2.22%/8.96% for SD, 0.18/0.21% for PD, and 0.31/0.54% for OR; pairwise p-values 7.8e-6 for SD ~ PD, 2.7e-4 for OR ~ PD, and 7.0e-3 for SD ~ OR, Wilcoxon rank sum test). The fold change of CCF distinguishes the patient groups even better, which increases for PD, remains stable for SD, and decreases for OR patients (p-values 0.002, ~ 1, and 0.0001 respectively, Wilcoxon signed-rank test). Eleven drug-related mutations are identified from nine out of the 313 patients. CONCLUSIONS: The ctDNA fingerprint method improves both specificity and sensitivity of monitoring treatment response across several tumor types. It can identify tumor relapse/recurrence potentially earlier than imaging-based diagnosis. When augmented with tumor hotspot genes, it can track acquired drug-related mutations in patients.

Integrated genomic and transcriptomic analysis revealed mutation patterns of de-differentiated liposarcoma and leiomyosarcoma.

BACKGROUND: Treating patients with advanced sarcomas is challenging due to great histologic diversity among its subtypes. Leiomyosarcoma (LMS) and de-differentiated liposarcoma (DDLPS) are two common and aggressive subtypes of soft tissue sarcoma (STS). They differ significantly in histology and clinical behaviors. However, the molecular driving force behind the difference is unclear. METHODS: We collected 20 LMS and 12 DDLPS samples and performed whole exome sequencing (WES) to obtain their somatic mutation profiles. We also performed RNA-Seq to analyze the transcriptomes of 8 each of the LMS and DDLPS samples and obtained information about differential gene expression, pathway enrichment, immune cell infiltration in tumor microenvironment, and chromosomal rearrangement including gene fusions. Selected gene fusion events from the RNA-seq prediction were checked by RT-PCR in tandem with Sanger sequencing. RESULTS: We detected loss of function mutation and deletion of tumor suppressors mostly in LMS, and oncogene amplification mostly in DDLPS. A focal amplification affecting chromosome 12q13-15 region which encodes MDM2, CDK4 and HMGA2 is notable in DDLPS. Mutations in TP53, ATRX, PTEN, and RB1 are identified in LMS but not DDLPS, while mutation of HERC2 is only identified in DDLPS but not LMS. RNA-seq revealed overexpression of MDM2, CDK4 and HMGA2 in DDLPS and down-regulation of TP53 and RB1 in LMS. It also detected more fusion events in DDLPS than LMS (4.5 vs. 1, p = 0.0195), and the ones involving chromosome 12 in DDLPS stand out. RT-PCR and Sanger sequencing verified the majority of the fusion events in DDLPS but only one event in LMS selected to be tested. The tumor microenvironmental signatures are highly correlated with histologic types. DDLPS has more endothelial cells and fibroblasts content than LMS. CONCLUSIONS: Our analysis revealed different recurrent genetic variations in LMS and DDLPS including simultaneous upregulation of gene expression and gene copy number amplification of MDM2 and CDK4. Up-regulation of tumor related genes is favored in DDLPS, while loss of suppressor function is favored in LMS. DDLPS harbors more frequent fusion events which can generate neoepitopes and potentially targeted by personalized immune treatment.

Analysis of circulating tumor DNA identifies distinct therapeutic response to intraperitoneal and intravenous paclitaxel plus S-1 in gastric cancer patients with peritoneal metastasis.

Background: Circulating tumor DNA (ctDNA) is a promising technique for predicting curative effects and monitoring tumor recurrence. The utility of ctDNA for gastric cancer with peritoneal dissemination remains elusive. Objectives: To assess the feasibility of ctDNA in predicting tumor response to chemotherapy in gastric cancer with peritoneal dissemination. Design: This was a prospective study. Methods: We enrolled 30 patients with gastric cancer peritoneal metastasis, treated with intraperitoneal and intravenous paclitaxel plus S-1. Peripheral blood samples of patients were prospectively collected at baseline, after treatment initiation accompanied by computed tomography scan and disease progression. Mutational profiles from ctDNA were analyzed to evaluate its association with chemotherapeutic response. Results: Tumor protein 53 (TP53) was the most frequently altered gene at baseline blood samples. Although baseline TP53 mutation was not related to therapeutic response, patients with TP53 mutation had worse progression-free survival (PFS) and overall survival (OS). Additionally, baseline ctDNA content fraction (CCF) was found to be significantly lower in responders than non-responders. Meanwhile, patients with high CCF had a trend of worse PFS and OS. Combining TP53 alteration and CCF, the prognosis of TP53-wt patients could be further stratified. Patients with CCF-low_TP53-wt had markedly longer survival than those with CCF-high_TP53-wt. Conclusion: Our study highlighted the significance of ctDNA in predicting potential clinical outcomes in gastric cancer patients during chemotherapy. Trial registration: ChiCTR-IIR-16009802 (Chinese Clinical Trial Registry).

A recurrence-predictive model based on eight genes and tumor mutational burden/microsatellite instability status in Stage II/III colorectal cancer.

BACKGROUND: Although adjuvant chemotherapy (ACT) is widely used to treat patients with Stage II/III colorectal cancer (CRC), administering ACT to specific patients remains a challenge. The decision to ACT requires an accurate assessment of recurrence risk and absolute treatment benefit. However, the traditional TNM staging system does not accurately assess a patient's individual risk of recurrence. METHODS: To identify recurrence risk-related genetic factors for Stage II/III CRC patients after radical surgery, we conducted an analysis of whole-exome sequencing of 47 patients with Stage II/III CRC who underwent radical surgery at five institutions. Patients were grouped into non-recurrence group (NR, n = 24, recurrence-free survival [RFS] > 5 years) and recurrence group (R, n = 23, RFS <2 years). The TCGA-COAD/READ cohort was employed as the validation dataset. RESULTS: A recurrence-predictive model (G8plus score) based on eight gene (CUL9, PCDHA12, HECTD3, DCX, SMARCA2, FAM193A, AATK, and SORCS2) mutations and tumor mutation burden/microsatellite instability (TMB/MSI) status was constructed, with 97.87% accuracy in our data and 100% negative predictive value in the TCGA-COAD/READ cohort. For the TCGA-COAD/READ cohort, the G8plus-high group had better RFS (HR = 0.22, p = 0.024); the G8plus-high tumors had significantly more infiltrated immune cell types, higher tertiary lymphoid structure signature scores, and higher immunological signature scores. The G8plus score was also a predict biomarker for immunotherapeutic in advanced CRC in the PUCH cohort. CONCLUSIONS: In conclusion, the G8plus score is a powerful biomarker for predicting the risk of recurrence in patients with stage II/III CRC. It can be used to stratify patients who benefit from ACT and immunotherapy.

Chronic obstructive pulmonary disease alters the genetic landscape and tumor immune microenvironment in lung cancer patients.

Background: Chronic obstructive pulmonary disease (COPD) and lung cancer are leading causes of morbidity and mortality worldwide. Studies have reported molecular alterations in patients with lung cancer and in patients with COPD. However, few investigation has been conducted on the molecular characteristics of lung cancer patients with COPD. Materials and methods: We performed a retrospective cohort study that included 435 patients with pathologically confirmed lung cancer at the Ruijin Hospital. For patients with documented spirometry, Global Initiative for Chronic Obstructive Lung Disease criteria were used to define COPD. For patients without documented spirometry, chest computed tomography and other clinical information were used to define COPD. Tumor tissue DNA was extracted from formalin-fixed paraffin-embedded samples. DNA mutation analysis, multiplex immunohistochemistry (mIHC), calculation of tumor mutational burden (TMB), mutant-allele tumor heterogeneity (MATH), and predication of neoantigens were performed. Results: Although SNV mutations in lung cancer patients with COPD (G1 group) were generally higher than those in lung cancer patients without COPD (G2 group), the difference in the number of mutations was insignificant between the two groups. Of the 35 mutated genes, the number of them was higher in G1 than in G2, except that of EGFR. PI3K-Akt signaling pathway was enriched from significantly different genes. While TMB and MATH levels were not significantly different, the tumor neoantigen burdenwas markedly higher in G1 than that in G2. The level of CD68+ macrophages was significant higher in the stroma and total areas in the G1 group than in G2 group. The level of CD8+ lymphocytes was markedly higher in the stroma and showed a clear tendency forhigher expression in the G1 group than inthe G2 group. No significant differences were observed for the level of programmed death-ligand 1+ (PD-L1+), programmed death 1+ (PD-1+), and CD68PD-L1 in the stroma, tumor and total areas. Conclusion: Our study revealed different genetic aberrations and pathways, higher neoantigen burden, and higher level of CD68+ macrophages and CD8+ T lymphocytes in lung cancer patients with COPD. Our investigation implies that the existence of COPD should be considered and immunotherapy is a potential choice when treating lung cancer patients with COPD.

Exploration of the immunogenetic landscape of hyperprogressive disease after combined immunotherapy in cancer patients.

The immune-genetic changes that occur in cancer patients experiencing hyperprogressive disease (HPD) during combined immunotherapy are unclear. In this study, HPD patients with pre- and post-HPD samples and non-HPD patients with solid tumors were molecularly characterized by genetic and tumor immune microenvironment (TiME) analyses of paired samples by whole-exome sequencing, RNA sequencing, and multiplex immunofluorescence. The genetic analysis of paired samples showed that almost all the tumor driver gene mutations were preserved between pre- and post-HPD tumors. HPD patients had higher frequencies of mutations in TP53 and CNN2, and a significantly higher mutant-allele tumor heterogeneity than non-HPD patients. Tumor IL-6 mRNA was upregulated in post-HPD samples vs. pre-HPD, accompanied by a potential immune suppressive TiME with an elevated M2/M1 ratio. Salvage treatment with irinotecan plus bevacizumab was effective in one HPD patient, who experienced prolonged survival. These genetic features and TiME characteristics might help identify the features of HPD after immunotherapy.

Orthogonal CRISPR-Cas tools for genome editing, inhibition, and CRISPR recording in zebrafish embryos.

The clustered regularly interspaced short palindromic repeats (CRISPR)-Cas universe continues to expand. The type II CRISPR-Cas system from Streptococcus pyogenes (SpyCas9) is the most widely used for genome editing due to its high efficiency in cells and organisms. However, concentrating on a single CRISPR-Cas system imposes limits on target selection and multiplexed genome engineering. We hypothesized that CRISPR-Cas systems originating from different bacterial species could operate simultaneously and independently due to their distinct single-guide RNAs (sgRNAs) or CRISPR-RNAs (crRNAs), and protospacer adjacent motifs (PAMs). Additionally, we hypothesized that CRISPR-Cas activity in zebrafish could be regulated through the expression of inhibitory anti-CRISPR (Acr) proteins. Here, we use a simple mutagenesis approach to demonstrate that CRISPR-Cas systems from S. pyogenes (SpyCas9), Streptococcus aureus (SauCas9), Lachnospiraceae bacterium (LbaCas12a, previously known as LbCpf1) are orthogonal systems capable of operating simultaneously in zebrafish. CRISPR systems from Acidaminococcus sp. (AspCas12a, previously known as AsCpf1) and Neisseria meningitidis (Nme2Cas9) were also active in embryos. We implemented multichannel CRISPR recording using three CRISPR systems and show that LbaCas12a may provide superior information density compared with previous methods. We also demonstrate that type II Acrs (anti-CRISPRs) are effective inhibitors of SpyCas9 in zebrafish. Our results indicate that at least five CRISPR-Cas systems and two anti-CRISPR proteins are functional in zebrafish embryos. These orthogonal CRISPR-Cas systems and Acr proteins will enable combinatorial and intersectional strategies for spatiotemporal control of genome editing and genetic recording in animals.

Gastrointestinal disturbance and effect of fecal microbiota transplantation in discharged COVID-19 patients.

BACKGROUND: To investigate the potential beneficial effect of fecal microbiota transplantation (FMT) on gastrointestinal symptoms, gut dysbiosis and immune status in discharged COVID-19 patients. CASE PRESENTATION: A total of 11 COVID-19 patients were recruited in April, 2020, about one month on average after they were discharged from the hospital. All subjects received FMT for 4 consecutive days by oral capsule administrations with 10 capsules for each day. In total, 5 out of 11 patients reported to be suffered from gastrointestinal symptoms, which were improved after FMT. After FMT, alterations of B cells were observed, which was characterized as decreased naive B cell (P = 0.012) and increased memory B cells (P = 0.001) and non-switched B cells (P = 0.012).The microbial community richness indicated by operational taxonomic units number, observed species and Chao1 estimator was marginally increased after FMT. Gut microbiome composition of discharged COVID-19 patients differed from that of the general population at both phylum and genera level, which was characterized with a lower proportion of Firmicutes (41.0%) and Actinobacteria (4.0%), higher proportion of Bacteroidetes (42.9%) and Proteobacteria (9.2%). FMT can partially restore the gut dysbiosis by increasing the relative abundance of Actinobacteria (15.0%) and reducing Proteobacteria (2.8%) at the phylum level. At the genera level, Bifidobacterium and Faecalibacterium had significantly increased after FMT. CONCLUSIONS: After FMT, altered peripheral lymphocyte subset, restored gut microbiota and alleviated gastrointestinal disorders were observe, suggesting that FMT may serve as a potential therapeutic and rehabilitative intervention for the COVID-19.

Whole-exome sequencing identifies prognostic mutational signatures in gastric cancer.

BACKGROUND: Gastric cancer (GC) is a heterogeneous disease, and is a leading cause of cancer deaths in Eastern Asia. Genomic analysis, such as whole-exome sequencing (WES), can help identify key genetic alterations leading to the malignancy and diversity of GC, and may help identify new drug targets. METHODS: We identified genomic alterations in a cohort of 38 GC patients, including 26 metastatic and 12 non-metastatic patients. We analyzed the association between novel gene mutations and copy number variations (CNVs) with tumor metastasis and patient survival. RESULTS: A number of significantly mutated genes in somatic and germline cells were identified. Among them, ATAD3B somatic mutation, a potential biomarker of immunotherapy in stomach cancers, was associated with better patient survival (P=0.0939) and metastasis (P=0.074). POLE germline variation was correlated with shorter overall survival (OS; P=0.0100). Novel CNVs were also identified and can potentially be used as biomarkers. These included 9p24.1 deletion (P=0.0376) and 16p11.2 amplification (P=0.0066), which were both associated with shorter OS. CNVs of several genes including MMP9, PTPN1, and SS18L1 were found to be significantly related to metastasis (P<0.05). CONCLUSIONS: We characterized the mutational landscape of 38 GC patients and discovered several potential new predictive markers of survival and metastasis in GC.

Establishment of a Gene Detection System for Hotspot Mutations of Hearing Loss.

Hearing loss is an etiologically heterogeneous trait with a high incidence in China. Though conventional newborn hearing screening program has been widely adopted, gene detection can significantly improve the means of early discovering genetic risk factors. Thus, simple and efficient methods with higher sensitivity and lower cost for detecting hotspot mutations of hearing loss are urgently requested. Here we established a mutation detection system based on multiple fluorescent probe technique, which can detect and genotype nine hotspot mutations of four prominent hearing loss-related genes in two reactions on a four-channel real-time PCR instrument, including GJB2 (rs750188782, rs80338943, rs1110333204, and rs80338939), GJB3 (rs74315319), SLC26A4 (rs111033313 and rs121908362), and mtDNA 12S rRNA (rs267606617 and rs267606619). This system is with high sensitivity that enables detecting as low as 10 DNA copies samples per reaction. A comparison study in 268 clinical samples showed that the detection system had 100% concordance to Sanger sequencing. Besides, blood and saliva samples can be directly detected without DNA extraction process, which greatly simplifies the manipulation. The new system with high sensitivity, accuracy, and specimen type compatibility can be expectedly a reliable tool in clinical application.