Precision targeting tumor cells using cancer-specific InDel mutations with CRISPR-Cas9.

An ideal cancer therapeutic strategy involves the selective killing of cancer cells without affecting the surrounding normal cells. However, researchers have failed to develop such methods for achieving selective cancer cell death because of shared features between cancerous and normal cells. In this study, we have developed a therapeutic strategy called the cancer-specific insertions-deletions (InDels) attacker (CINDELA) to selectively induce cancer cell death using the CRISPR-Cas system. CINDELA utilizes a previously unexplored idea of introducing CRISPR-mediated DNA double-strand breaks (DSBs) in a cancer-specific fashion to facilitate specific cell death. In particular, CINDELA targets multiple InDels with CRISPR-Cas9 to produce many DNA DSBs that result in cancer-specific cell death. As a proof of concept, we demonstrate here that CINDELA selectively kills human cancer cell lines, xenograft human tumors in mice, patient-derived glioblastoma, and lung patient-driven xenograft tumors without affecting healthy human cells or altering mouse growth.

PDX models of human lung squamous cell carcinoma: consideration of factors in preclinical and co-clinical applications.

BACKGROUND: Treatment of human lung squamous cell carcinoma (LUSC) using current targeted therapies is limited because of their diverse somatic mutations without any specific dominant driver mutations. These mutational diversities preventing the use of common targeted therapies or the combination of available therapeutic modalities would require a preclinical animal model of this tumor to acquire improved clinical responses. Patient-derived xenograft (PDX) models have been recognized as a potentially useful preclinical model for personalized precision medicine. However, whether the use of LUSC PDX models would be appropriate enough for clinical application is still controversial. METHODS: In the process of developing PDX models from Korean patients with LUSC, the authors investigated the factors influencing the successful initial engraftment of tumors in NOD scid gamma mice and the retainability of the pathological and genomic characteristics of the parental patient tumors in PDX tumors. CONCLUSIONS: The authors have developed 62 LUSC PDX models that retained the pathological and genomic features of parental patient tumors, which could be used in preclinical and co-clinical studies. Trial registration Tumor samples were obtained from 139 patients with LUSC between November 2014 and January 2019. All the patients provided signed informed consents. This study was approved by the institutional review board (IRB) of Samsung Medical Center (2018-03-110).

Recurrent fusion transcripts detected by whole-transcriptome sequencing of 120 primary breast cancer samples.

Relatively few recurrent gene fusion events have been associated with breast cancer to date. In an effort to uncover novel fusion transcripts, we performed whole-transcriptome sequencing of 120 fresh-frozen primary breast cancer samples and five adjacent normal breast tissues using the Illumina HiSeq2000 platform. Three different fusion-detecting tools (deFuse, Chimerascan, and TopHatFusion) were used, and the results were compared. These tools detected 3,831, 6,630 and 516 fusion transcripts (FTs) overall. We primarily focused on the results obtained using the deFuse software. More FTs were identified from HER2 subtype breast cancer samples than from the luminal or triple-negative subtypes (P < 0.05). Seventy fusion candidates were selected for validation, and 32 (45.7%) were confirmed by RT-PCR and Sanger sequencing. Of the validated fusions, six were recurrent (found in 2 or more samples), three were in-frame (PRDX1-AKR1A1, TACSTD2-OMA1, and C2CD2-TFF1) and three were off-frame (CEACAM7-CEACAM6, CYP4X1-CYP4Z2P, and EEF1DP3-FRY). Notably, the novel read-through fusion, EEF1DP3-FRY, was identified and validated in 6.7% (8/120) of the breast cancer samples. This off-frame fusion results in early truncation of the FRY gene, which plays a key role in the structural integrity during mitosis. Three previously reported fusions, PPP1R1B-STARD3, MFGE8-HAPL, and ETV6-NTRK3, were detected in 8.3, 3.3, and 0.8% of the 120 samples, respectively, by both deFuse and Chimerascan. The recently reported MAGI3-AKT3 fusion was not detected in our analysis. Although future work will be needed to examine the biological significance of our new findings, we identified a number of novel fusions and confirmed some previously reported fusions.

The genome of Diuraphis noxia, a global aphid pest of small grains.

BACKGROUND: The Russian wheat aphid, Diuraphis noxia Kurdjumov, is one of the most important pests of small grains throughout the temperate regions of the world. This phytotoxic aphid causes severe systemic damage symptoms in wheat, barley, and other small grains as a direct result of the salivary proteins it injects into the plant while feeding. RESULTS: We sequenced and de novo assembled the genome of D. noxia Biotype 2, the strain most virulent to resistance genes in wheat. The assembled genomic scaffolds span 393 MB, equivalent to 93% of its 421 MB genome, and contains 19,097 genes. D. noxia has the most AT-rich insect genome sequenced to date (70.9%), with a bimodal CpG(O/E) distribution and a complete set of methylation related genes. The D. noxia genome displays a widespread, extensive reduction in the number of genes per ortholog group, including defensive, detoxification, chemosensory, and sugar transporter groups in comparison to the Acyrthosiphon pisum genome, including a 65% reduction in chemoreceptor genes. Thirty of 34 known D. noxia salivary genes were found in this assembly. These genes exhibited less homology with those salivary genes commonly expressed in insect saliva, such as glucose dehydrogenase and trehalase, yet greater conservation among genes that are expressed in D. noxia saliva but not detected in the saliva of other insects. Genes involved in insecticide activity and endosymbiont-derived genes were also found, as well as genes involved in virus transmission, although D. noxia is not a viral vector. CONCLUSIONS: This genome is the second sequenced aphid genome, and the first of a phytotoxic insect. D. noxia's reduced gene content of may reflect the influence of phytotoxic feeding in shaping the D. noxia genome, and in turn in broadening its host range. The presence of methylation-related genes, including cytosine methylation, is consistent with other parthenogenetic and polyphenic insects. The D. noxia genome will provide an important contrast to the A. pisum genome and advance functional and comparative genomics of insects and other organisms.

The transcriptional landscape and mutational profile of lung adenocarcinoma.

All cancers harbor molecular alterations in their genomes. The transcriptional consequences of these somatic mutations have not yet been comprehensively explored in lung cancer. Here we present the first large scale RNA sequencing study of lung adenocarcinoma, demonstrating its power to identify somatic point mutations as well as transcriptional variants such as gene fusions, alternative splicing events, and expression outliers. Our results reveal the genetic basis of 200 lung adenocarcinomas in Koreans including deep characterization of 87 surgical specimens by transcriptome sequencing. We identified driver somatic mutations in cancer genes including EGFR, KRAS, NRAS, BRAF, PIK3CA, MET, and CTNNB1. Candidates for novel driver mutations were also identified in genes newly implicated in lung adenocarcinoma such as LMTK2, ARID1A, NOTCH2, and SMARCA4. We found 45 fusion genes, eight of which were chimeric tyrosine kinases involving ALK, RET, ROS1, FGFR2, AXL, and PDGFRA. Among 17 recurrent alternative splicing events, we identified exon 14 skipping in the proto-oncogene MET as highly likely to be a cancer driver. The number of somatic mutations and expression outliers varied markedly between individual cancers and was strongly correlated with smoking history of patients. We identified genomic blocks within which gene expression levels were consistently increased or decreased that could be explained by copy number alterations in samples. We also found an association between lymph node metastasis and somatic mutations in TP53. These findings broaden our understanding of lung adenocarcinoma and may also lead to new diagnostic and therapeutic approaches.

FX: an RNA-Seq analysis tool on the cloud.

UNLABELLED: FX is an RNA-Seq analysis tool, which runs in parallel on cloud computing infrastructure, for the estimation of gene expression levels and genomic variant calling. In the mapping of short RNA-Seq reads, FX uses a transcriptome-based reference primarily, generated from ~160 000 mRNA sequences from RefSeq, UCSC and Ensembl databases. This approach reduces the misalignment of reads originating from splicing junctions. Unmapped reads not aligned on known transcripts are then mapped on the human genome reference. FX allows analysis of RNA-Seq data on cloud computing infrastructures, supporting access through a user-friendly web interface. AVAILABILITY: FX is freely available on the web at (http://fx.gmi.ac.kr), and can be installed on local Hadoop clusters. Guidance for the installation and operation of FX can be found under the 'Documentation' menu on the website. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Exome sequencing is an efficient tool for genetic screening of Charcot-Marie-Tooth disease.

Charcot-Marie-Tooth disease (CMT) is one of the most common inherited neuropathies and is a genetically and clinically heterogeneous disorder with variable inheritance modes. As several molecules have been reported to have therapeutic effects on CMT, depending on the underlying genetic causes, exact genetic diagnostics have become very important for executing personalized therapy. Whole-exome sequencing has recently been introduced as an available method to identify rare or novel genetic defects from genetic disorders. Particularly, CMT is a model disease to apply exome sequencing because more than 50 genes (loci) are involved in its development with weak genotype-phenotype correlation. This study performed the exome sequencing in 25 unrelated CMT patients who revealed neither 17p12 duplication/deletion nor several major CMT genes. This study identified eight causative heterozygous mutations (32%). This detection rate seems rather high because each sample was tested before the study for major genetic causes. Therefore, this study suggests that the exome sequencing can be a highly exact, rapid, and economical molecular diagnostic tool for CMT patients who are tested for major genetic causes.

Exome sequencing and subsequent association studies identify five amino acid-altering variants influencing human height.

Height is a highly heritable trait that involves multiple genetic loci. To identify causal variants that influence stature, we sequenced whole exomes of four children with idiopathic short stature. Ninety-five nonsynonymous single-nucleotide polymorphisms (nsSNPs) were selected as potential candidate variants. We performed association analysis in 740 cohort individuals and identified 11 nsSNPs in 10 loci (DIS3L2, ZBTB38, FAM154A, PTCH1, TSSC4, KIF18A, GPR133, ACAN, FAM59A, and NINL) associated with adult height (P < 0.05), including five novel loci. Of these, two nsSNPs (TSSC4 and KIF18A loci) were significant at P < 0.05 in the replication study (n = 1,000) and five (ZBTB38, FAM154A, TSSC4, KIF18A, and FAM59A loci) were significant at P < 0.01 in the combined analysis (n = 1,740). Together, the five nsSNPs accounted for approximately 2.5% of the height variation. This study demonstrated the utility of next-generation sequencing in identifying genetic variants and loci associated with complex traits.

Erratum: Interleukin-34 Limits the Therapeutic Effects of Immune Checkpoint Blockade.

[This corrects the article DOI: 10.1016/j.isci.2020.101584.].

Identification of frequently mutated genes with relevance to nonsense mediated mRNA decay in the high microsatellite instability cancers.

Frameshift mutations at coding mononucleotide repeats (cMNR) are frequent in high-microsatellite instability (MSI-H) cancers. Frameshift mutations in cMNR result in the formation of a premature termination codon (PTC) in the transcribed mRNA, and these abnormal mRNAs are generally degraded by nonsense mediated mRNA decay (NMD). We have identified novel genes that are frequently mutated at their cMNR by blocking NMD in two MSI-H cancer cell lines. After blocking NMD, we screened for differentially expressed genes using DNA microarrays, and then used database analysis to select 28 candidate genes containing cMNR with more than 9 nucleotide repeats. cMNR mutations have not been previously reported in MSI-H cancers for 15 of the 28 genes. We analyzed the cMNR mutation of each of the 15 genes in 10 MSI-H cell lines and 21 MSI-H cancers, and found frequent mutations of 12 genes in MSI-H cell lines and cancers, but not in microsatellite stable (MSS) cancers. Among these genes, the most frequently mutated in MSI-H cell lines were MLL3 (70%), PHACTR4 (70%), RUFY2 (50%) and TBC1D23 (50%). MLL3, which has already been implicated in cancer, had the highest mutation frequency in MSI-H cancers (48%). Our combined approach of NMD block, database search, and mutation analysis has identified a large number of genes mutated in their cMNR in MSI-H cancers. The identified mutations are expected to contribute to MSI-H tumorigenesis by causing an absence of gene expression or low gene dosage effects.

Gene expression changes in patient-matched gastric normal mucosa, adenomas, and carcinomas.

A subset of gastric carcinomas shows histologic evidence of a multistep process, progressing from gastric adenoma to gastric carcinoma. We examined gene expression changes during the gastric adenoma-carcinoma sequence in 26 snap-frozen samples (normal mucosa, adenoma, and carcinoma samples from eight patients and two additional carcinomas) by oligonucleotide microarray. Unsupervised hierarchical clustering analysis demonstrated differential gene expression between gastric normal mucosa, adenomas and carcinomas. We identified 319 and 422 genes differentially regulated in adenoma and carcinoma, respectively, relative to normal mucosa, using a combination of Welch's t-test and fold-change analysis. Applying a combination of robust multi-category support vector machines to the data, reveal that 39 and 21 genes were gradually up- and down-regulated, respectively, in succession in normal mucosa, adenoma, and carcinoma samples. We validated gene expression levels of four genes: hydroxyprostaglandin dehydrogenase 15 (HPGD), follistatin-like 1, trefoil factor 1 (TTF1) and trefoil factor 2 (TFF2) by RT-PCR and found direct correlation with microarray results. The expressions of the TFF2 and HPGD genes were further evaluated by immunohistochemistry in 103 adenomas and 70 carcinomas; expression of both proteins was decreased in these tissues. The progressive alteration in gene expression in the transition from normal mucosa to carcinoma suggests that these changes may play critical roles in gastric carcinogenesis.

MicroRNA expression profile of gastrointestinal stromal tumors is distinguished by 14q loss and anatomic site.

MicroRNAs are known to regulate gene expression. Although unique microRNA expression profiles have been reported in several tumors, little is known about microRNA expression profiles in GISTs. To evaluate the relationship between microRNA expression and clinicopathologic findings of GISTs, we analyzed the microRNA expression profiles of GISTs. We used fresh frozen tissues from 20 GISTs and analyzed KIT and PDGFRA mutations and chromosomal loss status. MicroRNA expression was analyzed using a microRNA chip containing 470 microRNAs. Using unsupervised hierarchical clustering analysis, we found four distinct microRNA expression patterns in our 20 GISTs. Six GISTs that did not have 14q loss formed a separate cluster. In the 14 GISTs with 14q loss, 5 small bowel GISTs formed a separate cluster and the remaining 9 GISTs could be divided into two groups according to frequent chromosomal losses and tumor risk. We found 73 microRNAs that were significantly down-regulated in the GISTs with 14q loss; 38 of these microRNAs are encoded on 14q. We also found many microRNAs that were down-regulated in small bowel and high-risk group GISTs. Most of the microRNAs down-regulated in the high-risk group and small bowel GISTs are known to be involved in tumor progression, specifically by stimulating mitogen-activated protein kinase (MAPK) and the cell cycle. The microRNA expression patterns of GISTs are closely related to the status of 14q loss, anatomic site, and tumor risk. These findings suggest that microRNA expression patterns can differentiate several subsets of GISTs.

Validity of patient-derived xenograft mouse models for lung cancer based on exome sequencing data.

Patient-derived xenograft (PDX) mouse models are frequently used to test the drug efficacy in diverse types of cancer. They are known to recapitulate the patient characteristics faithfully, but a systematic survey with a large number of cases is yet missing in lung cancer. Here we report the comparison of genomic characters between mouse and patient tumor tissues in lung cancer based on exome sequencing data. We established PDX mouse models for 132 lung cancer patients and performed whole exome sequencing for trio samples of tumor-normal-xenograft tissues. Then we computed the somatic mutations and copy number variations, which were used to compare the PDX and patient tumor tissues. Genomic and histological conclusions for validity of PDX models agreed in most cases, but we observed eight (~7%) discordant cases. We further examined the changes in mutations and copy number alterations in PDX model production and passage processes, which highlighted the clonal evolution in PDX mouse models. Our study shows that the genomic characterization plays complementary roles to the histological examination in cancer studies utilizing PDX mouse models.

Genome-wide association study of metabolic syndrome in Korean populations.

Metabolic syndrome (MetS) which is caused by obesity and insulin resistance, is well known for its predictive capability for the risk of type 2 diabetes mellitus and cardiovascular disease. The development of MetS is associated with multiple genetic factors, environmental factors and lifestyle. We performed a genome-wide association study to identify single-nucleotide polymorphism (SNP) related to MetS in large Korean population based samples of 1,362 subjects with MetS and 6,061 controls using the Axiom® Korean Biobank Array 1.0. We replicated the data in another sample including 502 subjects with MetS and 1,751 controls. After adjusting for age and sex, rs662799 located in the APOA5 gene were significantly associated with MetS. 15 SNPs in GCKR, C2orf16, APOA5, ZPR1, and BUD13 were associated with high triglyceride (TG). 14 SNPs in APOA5, ALDH1A2, LIPC, HERPUD1, and CETP, and 2 SNPs in MTNR1B were associated with low high density lipoprotein cholesterol (HDL-C) and high fasting blood glucose respectively. Among these SNPs, 6 TG SNPs: rs1260326, rs1260333, rs1919127, rs964184, rs2075295 and rs1558861 and 11 HDL-C SNPs: rs4775041, rs10468017, rs1800588, rs72786786, rs173539, rs247616, rs247617, rs3764261, rs4783961, rs708272, and rs7499892 were first discovered in Koreans. Additional research is needed to confirm these 17 novel SNPs in Korean population.

Interleukin-34 Limits the Therapeutic Effects of Immune Checkpoint Blockade.

Interleukin-34 (IL-34) is an alternative ligand to colony-stimulating factor-1 (CSF-1) for the CSF-1 receptor that acts as a key regulator of monocyte/macrophage lineage. In this study, we show that tumor-derived IL-34 mediates resistance to immune checkpoint blockade regardless of CSF-1 existence in various murine cancer models. Consistent with its immunosuppressive characteristics, the expression of IL-34 in tumors correlates with decreased frequencies of cellular (such as CD8+ and CD4+ T cells and M1-biased macrophages) and molecular (including various cytokines and chemokines) effectors at the tumor microenvironment. Then, a neutralizing antibody against IL-34 improved the therapeutic effects of the immune checkpoint blockade in combinatorial therapeutic models, including a patient-derived xenograft model. Collectively, we revealed that tumor-derived IL-34 inhibits the efficacy of immune checkpoint blockade and proposed the utility of IL-34 blockade as a new strategy for cancer therapy.

MedRefSNP: a database of medically investigated SNPs.

Genetic association studies and linkage analyses using single nucleotide polymorphisms (SNPs) are rapidly increasing in number, and the results are important for evaluating the utility of SNPs in the biomedical sciences. Although many SNP databases have been established, there is no database focusing on published SNPs, where the research results of scientific investigations are available. To enhance the utilization of such SNP data, we have developed the MedRefSNP database (http://www.medclue.com/medrefsnp) to provide integrated information about SNPs collected from the PubMed and OMIM databases. The RefSNP identifiers are automatically identified and are linked to various information sources such as the dbSNP, the HapMap database, the Entrez Gene database, the UCSC genome browser, the CGAP Pathway Searcher, and genetic association databases. And, each SNP is checked to determine whether the PolyDoms, SNPs3D or PolyPhen databases predicts that the SNP affects the phenotype of the protein encoded by the gene carrying the SNP. Also, neighboring SNPs showing strong linkage disequilibrium (LD) with published SNPs are included, using HapMap data. Currently, 36199 unique SNPs (including 31368 neighboring SNPs) collected from 25906 PubMed abstracts and 590 OMIM records are stored along with 2491 human genes related to 466 molecular pathways. The MedRefSNP database will help researchers to review previously investigated results more efficiently, and will expand knowledge by using the genomic and functional contexts of the SNPs.

Impacts of Salt Stress on Locomotor and Transcriptomic Responses in the Intertidal Gastropod Batillaria attramentaria.

Salinity is one of the most crucial environmental factors that structures biogeographic boundaries of aquatic organisms, affecting distribution, abundance, and behavior. However, the association between behavior and gene regulation underlying acclimation to changes in salinity remains poorly understood. In this study, we investigated the effects of salinity stress on behavior (movement distance) and patterns of gene expression (using RNA sequencing) of the intertidal gastropod Batillaria attramentaria. We examined responses to short-term (1-hour) and long-term (30-day) acclimation to a range of salinities (43, 33 [control], 23, 13, and 3 psu). We found that the intertidal B. attramentaria is able to tolerate a broad range of salinity from 13 to 43 psu but not the acute low salinity of 3 psu. Behavioral experiments showed that salt stress significantly influenced snails' movement, with lower salinity resulting in shorter movement distance. Transcriptomic analyses revealed critical metabolic pathways and genes potentially involved in acclimation to salinity stress, including ionic and osmotic regulation, signal and hormonal transduction pathways, water exchange, cell protection, and gene regulation or epigenetic modification. In general, our study presents a robust, integrative laboratory-based approach to investigate the effects of salt stress on a nonmodel gastropod facing detrimental consequences of environmental change. The current genetic results provide a wealth of reference data for further research on mechanisms of ionic and osmotic regulation and adaptive evolution of this coastal gastropod.

A simple and accurate SNP scoring strategy based on typeIIS restriction endonuclease cleavage and matrix-assisted laser desorption/ionization mass spectrometry.

BACKGROUND: We describe the development of a novel matrix-assisted laser desorption ionization time-of-flight (MALDI-TOF)-based single nucleotide polymorphism (SNP) scoring strategy, termed Restriction Fragment Mass Polymorphism (RFMP) that is suitable for genotyping variations in a simple, accurate, and high-throughput manner. The assay is based on polymerase chain reaction (PCR) amplification and mass measurement of oligonucleotides containing a polymorphic base, to which a typeIIS restriction endonuclease recognition was introduced by PCR amplification. Enzymatic cleavage of the products leads to excision of oligonucleotide fragments representing base variation of the polymorphic site whose masses were determined by MALDI-TOF MS. RESULTS: The assay represents an improvement over previous methods because it relies on the direct mass determination of PCR products rather than on an indirect analysis, where a base-extended or fluorescent report tag is interpreted. The RFMP strategy is simple and straightforward, requiring one restriction digestion reaction following target amplification in a single vessel. With this technology, genotypes are generated with a high call rate (99.6%) and high accuracy (99.8%) as determined by independent sequencing. CONCLUSION: The simplicity, accuracy and amenability to high-throughput screening analysis should make the RFMP assay suitable for large-scale genotype association study as well as clinical genotyping in laboratories.

Differential gene expression profiles of metastases in paired primary and metastatic colorectal carcinomas.

BACKGROUND AND METHODS: Despite the overwhelming clinical significance of metastases, the cellular and molecular mechanisms involved are largely unknown. In order to define significant differences between primary colon carcinomas and their metastases, we analyzed gene expression profiles of 12 sets of triple-paired tissues using 19 K human oligonucleotide microarrays. A total of 36 microarray experiments were analyzed by unsupervised two-way hierarchical clustering and multi-dimensional scaling (MDS). RESULTS: Both methods completely distinguished normal mucosa from carcinoma, but failed to demonstrate a complete classification of primary and metastatic carcinomas. We found a separable tendency to be classified into the primary and metastatic colon carcinomas by MDS. In supervised hierarchical clustering, we identified 80 genes that were differentially expressed between paired primary and metastatic colon carcinomas. The 80 identified genes also successfully distinguished three validation sets of primary and lung-metastatic colon carcinomas. A specific set of genes was identified that distinguished the metastasis from the corresponding primary tumor in nearly half of the metastases analyzed. CONCLUSIONS: We suggest that a more accurate model of the metastatic potential is based on a global tumor expression pattern along with the appearance of distinct metastatic variants. This molecular profiling may be useful for the future study of colon cancer metastasis.

Macrophage migration inhibitory factor may be used as an early diagnostic marker in colorectal carcinomas.

Recent genetic studies have identified many differentially expressed genes in colorectal carcinomas. For validation of up-regulated genes in colorectal carcinomas, we performed an enzyme-linked immunosorbent assay. Candidate markers were selected from gene expression data for 40 colorectal cancers and 35 matched normal mucosal samples. Based on intensive filtering, 9 genes were selected for the further evaluations. Among them, macrophage migration inhibitory factor (MIF), inhibin betaA, and chemokine ligand 10 were screened, and the results were compared with carcinoembryonic antigen (CEA) in serum samples of 129 patients with colon cancer and 53 healthy control subjects. We found that the serum MIF level was significantly increased in patients with colorectal cancer. Compared with CEA, MIF was more sensitive in early cancer detection (47.3% vs 29.5%). However, the specificity was not as high as that of CEA (90.6% vs 100.0%). Our findings indicate that MIF may be used as a diagnostic marker in colorectal carcinomas.