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.

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

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

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.

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).

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.

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.

The genome of common long-arm octopus Octopus minor.

Background: The common long-arm octopus (Octopus minor) is found in mudflats of subtidal zones and faces numerous environmental challenges. The ability to adapt its morphology and behavioral repertoire to diverse environmental conditions makes the species a promising model for understanding genomic adaptation and evolution in cephalopods. Findings: The final genome assembly of O. minor is 5.09 Gb, with a contig N50 size of 197 kb and longest size of 3.027 Mb, from a total of 419 Gb raw reads generated using the Pacific Biosciences RS II platform. We identified 30,010 genes; 44.43% of the genome is composed of repeat elements. The genome-wide phylogenetic tree indicated the divergence time between O. minor and Octopus bimaculoides was estimated to be 43 million years ago based on single-copy orthologous genes. In total, 178 gene families are expanded in O. minor in the 14 bilaterian species. Conclusions: We found that the O. minor genome was larger than that of closely related O. bimaculoides, and this difference could be explained by enlarged introns and recently diversified transposable elements. The high-quality O. minor genome assembly provides a valuable resource for understanding octopus genome evolution and the molecular basis of adaptations to mudflats.

Metabolic Syndrome Prediction Using Machine Learning Models with Genetic and Clinical Information from a Nonobese Healthy Population.

The prevalence of metabolic syndrome (MS) in the nonobese population is not low. However, the identification and risk mitigation of MS are not easy in this population. We aimed to develop an MS prediction model using genetic and clinical factors of nonobese Koreans through machine learning methods. A prediction model for MS was designed for a nonobese population using clinical and genetic polymorphism information with five machine learning algorithms, including naïve Bayes classification (NB). The analysis was performed in two stages (training and test sets). Model A was designed with only clinical information (age, sex, body mass index, smoking status, alcohol consumption status, and exercise status), and for model B, genetic information (for 10 polymorphisms) was added to model A. Of the 7,502 nonobese participants, 647 (8.6%) had MS. In the test set analysis, for the maximum sensitivity criterion, NB showed the highest sensitivity: 0.38 for model A and 0.42 for model B. The specificity of NB was 0.79 for model A and 0.80 for model B. In a comparison of the performances of models A and B by NB, model B (area under the receiver operating characteristic curve [AUC] = 0.69, clinical and genetic information input) showed better performance than model A (AUC = 0.65, clinical information only input). We designed a prediction model for MS in a nonobese population using clinical and genetic information. With this model, we might convince nonobese MS individuals to undergo health checks and adopt behaviors associated with a preventive lifestyle.

miR-204 downregulates EphB2 in aging mouse hippocampal neurons.

Hippocampal synaptic function and plasticity deteriorate with age, often resulting in learning and memory deficits. As MicroRNAs (miRNAs) are important regulators of neuronal protein expression, we examined whether miRNAs may contribute to this age-associated decline in hippocampal function. We first compared the small RNA transcriptome of hippocampal tissues from young and old mice. Among 269 hippocampal miRNAs, 80 were differentially expressed (≥ twofold) among the age groups. We focused on 36 miRNAs upregulated in the old mice compared with those in the young mice. The potential targets of these 36 miRNAs included 11 critical Eph/Ephrin synaptic signaling components. The expression levels of several genes in the Eph/Ephrin pathway, including EphB2, were significantly downregulated in the aged hippocampus. EphB2 is a known regulator of synaptic plasticity in hippocampal neurons, in part by regulating the surface expression of the NMDA receptor NR1 subunit. We found that EphB2 is a direct target of miR-204 among miRNAs that were upregulated with age. The transfection of primary hippocampal neurons with a miR-204 mimic suppressed both EphB2 mRNA and protein expression and reduced the surface expression of NR1. Transfection of miR-204 also decreased the total expression of NR1. miR-204 induces senescence-like phenotype in fully matured neurons as evidenced by an increase in p16-positive cells. We suggest that aging is accompanied by the upregulation of miR-204 in the hippocampus, which downregulates EphB2 and results in reduced surface and total NR1 expression. This mechanism may contribute to age-associated decline in hippocampal synaptic plasticity and the related cognitive functions.

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.

Brain somatic mutations in MTOR cause focal cortical dysplasia type II leading to intractable epilepsy.

Focal cortical dysplasia type II (FCDII) is a sporadic developmental malformation of the cerebral cortex characterized by dysmorphic neurons, dyslamination and medically refractory epilepsy. It has been hypothesized that FCD is caused by somatic mutations in affected regions. Here, we used deep whole-exome sequencing (read depth, 412-668×) validated by site-specific amplicon sequencing (100-347,499×) in paired brain-blood DNA from four subjects with FCDII and uncovered a de novo brain somatic mutation, mechanistic target of rapamycin (MTOR) c.7280T>C (p.Leu2427Pro) in two subjects. Deep sequencing of the MTOR gene in an additional 73 subjects with FCDII using hybrid capture and PCR amplicon sequencing identified eight different somatic missense mutations found in multiple brain tissue samples of ten subjects. The identified mutations accounted for 15.6% of all subjects with FCDII studied (12 of 77). The identified mutations induced the hyperactivation of mTOR kinase. Focal cortical expression of mutant MTOR by in utero electroporation in mice was sufficient to disrupt neuronal migration and cause spontaneous seizures and cytomegalic neurons. Inhibition of mTOR with rapamycin suppressed cytomegalic neurons and epileptic seizures. This study provides, to our knowledge, the first evidence that brain somatic activating mutations in MTOR cause FCD and identifies mTOR as a treatment target for intractable epilepsy in FCD.

Gene expression defines natural changes in mammalian lifespan.

Mammals differ more than 100-fold in maximum lifespan, which can be altered in either direction during evolution, but the molecular basis for natural changes in longevity is not understood. Divergent evolution of mammals also led to extensive changes in gene expression within and between lineages. To understand the relationship between lifespan and variation in gene expression, we carried out RNA-seq-based gene expression analyses of liver, kidney, and brain of 33 diverse species of mammals. Our analysis uncovered parallel evolution of gene expression and lifespan, as well as the associated life-history traits, and identified the processes and pathways involved. These findings provide direct insights into how nature reversibly adjusts lifespan and other traits during adaptive radiation of lineages.

Differential expression of microRNAs and their target genes in non-small-cell lung cancer.

MicroRNAs (miRNAs) are single­stranded RNA species that constitute a class of non­coding RNAs, and are emerging as key regulators of gene expression. Since each miRNA is capable of regulating multiple genes, miRNAs are attractive markers for studies of coordinated gene expression. In this study, we investigated miRNA expression profiling using a massively parallel sequencing technique to compare non­small­cell lung cancer (NSCLC) tissue and normal lung tissue. Lung cancer tissue and normal lung tissue were obtained from nine NSCLC patients. RNA isolated from these samples was processed using RNA sequencing (RNA Seq) and the HiSeq 2000 system. Differentially expressed miRNAs and mRNAs were analyzed using a t­test. We selected target pairs that showed a negative correlation among significantly differentially expressed miRNAs and their putative target mRNAs using miRBase Targets. The differences in the expression levels of 222 miRNAs and 1,597 genes were statistically significant, as indicated by an absolute fold change ≥1.5 and P<0.05. miR­577, miR­301b, miR­944, miR­891a and miR­615­3p were generally upregulated, and miR­338­3p was generally downregulated. miRNA­mRNA target pair analysis revealed that 49 miRNAs had 696 target mRNAs. There were significantly differentially expressed miRNAs and mRNAs between lung cancer and normal tissue. Further investigation of miRNAs and their target genes is warranted to better understand NSCLC.

An epigenomic roadmap to induced pluripotency reveals DNA methylation as a reprogramming modulator.

Reprogramming of somatic cells to induced pluripotent stem cells involves a dynamic rearrangement of the epigenetic landscape. To characterize this epigenomic roadmap, we have performed MethylC-seq, ChIP-seq (H3K4/K27/K36me3) and RNA-Seq on samples taken at several time points during murine secondary reprogramming as part of Project Grandiose. We find that DNA methylation gain during reprogramming occurs gradually, while loss is achieved only at the ESC-like state. Binding sites of activated factors exhibit focal demethylation during reprogramming, while ESC-like pluripotent cells are distinguished by extension of demethylation to the wider neighbourhood. We observed that genes with CpG-rich promoters demonstrate stable low methylation and strong engagement of histone marks, whereas genes with CpG-poor promoters are safeguarded by methylation. Such DNA methylation-driven control is the key to the regulation of ESC-pluripotency genes, including Dppa4, Dppa5a and Esrrb. These results reveal the crucial role that DNA methylation plays as an epigenetic switch driving somatic cells to pluripotency.

Targeted exon sequencing fails to identify rare coding variants with large effect in rheumatoid arthritis.

INTRODUCTION: Although it has been suggested that rare coding variants could explain the substantial missing heritability, very few sequencing studies have been performed in rheumatoid arthritis (RA). We aimed to identify novel functional variants with rare to low frequency using targeted exon sequencing of RA in Korea. METHODS: We analyzed targeted exon sequencing data of 398 genes selected from a multifaceted approach in Korean RA patients (n = 1,217) and controls (n = 717). We conducted a single-marker association test and a gene-based analysis of rare variants. For meta-analysis or enrichment tests, we also used ethnically matched independent samples of Korean genome-wide association studies (GWAS) (n = 4,799) or immunochip data (n = 4,722). RESULTS: After stringent quality control, we analyzed 10,588 variants of 398 genes from 1,934 Korean RA case controls. We identified 13 nonsynonymous variants with nominal association in single-variant association tests. In a meta-analysis, we did not find any novel variant with genome-wide significance for RA risk. Using a gene-based approach, we identified 17 genes with nominal burden signals. Among them, VSTM1 showed the greatest association with RA (P = 7.80 × 10-4). In the enrichment test using Korean GWAS, although the significant signal appeared to be driven by total genic variants, we found no evidence for enriched association of coding variants only with RA. CONCLUSIONS: We were unable to identify rare coding variants with large effect to explain the missing heritability for RA in the current targeted resequencing study. Our study raises skepticism about exon sequencing of targeted genes for complex diseases like RA.

Heterozygous mutations in cyclic AMP phosphodiesterase-4D (PDE4D) and protein kinase A (PKA) provide new insights into the molecular pathology of acrodysostosis.

Acrodysostosis without hormone resistance is a rare skeletal disorder characterized by brachydactyly, nasal hypoplasia, mental retardation and occasionally developmental delay. Recently, loss-of-function mutations in the gene encoding cAMP-hydrolyzing phosphodiesterase-4D (PDE4D) have been reported to cause this rare condition but the pathomechanism has not been fully elucidated. To understand the pathogenetic mechanism of PDE4D mutations, we conducted 3D modeling studies to predict changes in the binding efficacy of cAMP to the catalytic pocket in PDE4D mutants. Our results indicated diminished enzyme activity in the two mutants we analyzed (Gly673Asp and Ile678Thr; based on PDE4D4 residue numbering). Ectopic expression of PDE4D mutants in HEK293 cells demonstrated this reduction in activity, which was identified by increased cAMP levels. However, the cells from an acrodysostosis patient showed low cAMP accumulation, which resulted in a decrease in the phosphorylated cAMP Response Element-Binding Protein (pCREB)/CREB ratio. The reason for this discrepancy was due to a compensatory increase in expression levels of PDE4A and PDE4B isoforms, which accounted for the paradoxical decrease in cAMP levels in the patient cells expressing mutant isoforms with a lowered PDE4D activity. Skeletal radiographs of 10-week-old knockout (KO) rats showed that the distal part of the forelimb was shorter than in wild-type (WT) rats and that all the metacarpals and phalanges were also shorter in KO, as the name acrodysostosis implies. Like the G-protein α-stimulatory subunit and PRKAR1A, PDE4D critically regulates the cAMP signal transduction pathway and influences bone formation in a way that activity-compromising PDE4D mutations can result in skeletal dysplasia. We propose that specific inhibitory PDE4D mutations can lead to the molecular pathology of acrodysostosis without hormone resistance but that the pathological phenotype may well be dependent on an over-compensatory induction of other PDE4 isoforms that can be expected to be targeted to different signaling complexes and exert distinct effects on compartmentalized cAMP signaling.

Whole genome and global gene expression analyses of the model mushroom Flammulina velutipes reveal a high capacity for lignocellulose degradation.

Flammulina velutipes is a fungus with health and medicinal benefits that has been used for consumption and cultivation in East Asia. F. velutipes is also known to degrade lignocellulose and produce ethanol. The overlapping interests of mushroom production and wood bioconversion make F. velutipes an attractive new model for fungal wood related studies. Here, we present the complete sequence of the F. velutipes genome. This is the first sequenced genome for a commercially produced edible mushroom that also degrades wood. The 35.6-Mb genome contained 12,218 predicted protein-encoding genes and 287 tRNA genes assembled into 11 scaffolds corresponding with the 11 chromosomes of strain KACC42780. The 88.4-kb mitochondrial genome contained 35 genes. Well-developed wood degrading machinery with strong potential for lignin degradation (69 auxiliary activities, formerly FOLymes) and carbohydrate degradation (392 CAZymes), along with 58 alcohol dehydrogenase genes were highly expressed in the mycelium, demonstrating the potential application of this organism to bioethanol production. Thus, the newly uncovered wood degrading capacity and sequential nature of this process in F. velutipes, offer interesting possibilities for more detailed studies on either lignin or (hemi-) cellulose degradation in complex wood substrates. The mutual interest in wood degradation by the mushroom industry and (ligno-)cellulose biomass related industries further increase the significance of F. velutipes as a new model.

RNA sequencing identifies novel markers of non-small cell lung cancer.

INTRODUCTION: The development of reliable gene expression profiling technology increasingly impacts our understanding of lung cancer biology. Here, we used RNA sequencing (RNA-Seq) to compare the transcriptomes of non-small cell lung cancer (NSCLC) and normal lung tissues and to investigate expression in lung cancer tissues. METHODS: We enrolled 88 male patients (mean age, 61.2 years) with NSCLC. RNA-Seq was performed on 88 pairs of NSCLC tumor tissue and non-tumor tissue from 54 patients with adenocarcinoma and 34 patients with squamous cell carcinoma. Immunohistochemistry was performed to validate differential candidate gene expression in a different NSCLC group. RESULTS: RNA-Seq produced 25.41 × 10(6) (± 8.90 × 10(6)) reads in NSCLC tissues and 24.70×10(6) (± 4.70 × 10(6)) reads in normal lung tissues [mean (± standard deviation)]. Among the genes expressed in both tissues, 335 were upregulated and 728 were downregulated ≥ 2-fold (p < 0.001). Four upregulated genes - CBX3, GJB2, CRABP2, and DSP - not previously reported in lung cancer were studied further. Their altered expression was verified by immunohistochemistry in a different set of NSCLC tissues (n = 154). CBX3 was positive in 90.3% (139 cases) of the samples; GJB2, in 22.7% (35 cases); CRABP2, in 72.1% (111 cases); and DSP, in 17.5% (27 cases). The positive rate of CRABP2 was higher in adenocarcinoma than squamous cell carcinoma (p < 0.01). CONCLUSIONS: CBX3 and CRABP2 expression was markedly increased in lung cancer tissues and especially CRABP2 may be promising candidate genes in lung adenocarcinoma.