Search tips
Search criteria 


Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
PLoS One. 2013; 8(2): e55596.
Published online 2013 February 6. doi:  10.1371/journal.pone.0055596
PMCID: PMC3566005

A High-Dimensional, Deep-Sequencing Study of Lung Adenocarcinoma in Female Never-Smokers

William C. S. Cho, Editor



Deep sequencing techniques provide a remarkable opportunity for comprehensive understanding of tumorigenesis at the molecular level. As omics studies become popular, integrative approaches need to be developed to move from a simple cataloguing of mutations and changes in gene expression to dissecting the molecular nature of carcinogenesis at the systemic level and understanding the complex networks that lead to cancer development.


Here, we describe a high-throughput, multi-dimensional sequencing study of primary lung adenocarcinoma tumors and adjacent normal tissues of six Korean female never-smoker patients. Our data encompass results from exome-seq, RNA-seq, small RNA-seq, and MeDIP-seq. We identified and validated novel genetic aberrations, including 47 somatic mutations and 19 fusion transcripts. One of the fusions involves the c-RET gene, which was recently reported to form fusion genes that may function as drivers of carcinogenesis in lung cancer patients. We also characterized gene expression profiles, which we integrated with genomic aberrations and gene regulations into functional networks. The most prominent gene network module that emerged indicates that disturbances in G2/M transition and mitotic progression are causally linked to tumorigenesis in these patients. Also, results from the analysis strongly suggest that several novel microRNA-target interactions represent key regulatory elements of the gene network.


Our study not only provides an overview of the alterations occurring in lung adenocarcinoma at multiple levels from genome to transcriptome and epigenome, but also offers a model for integrative genomics analysis and proposes potential target pathways for the control of lung adenocarcinoma.


Recent advances in DNA sequencing technology have revolutionized genomics and biomedical research, especially in the field of cancer research [1]. Various types of mutations as well as large scale chromosomal aberrations are being reported and cataloged, and the rate of data accumulation will likely accelerate for the foreseeable future. This should certainly apply to lung cancer which is currently the second most common cancer and the primary cause of mortality among cancer-related death in the United States [2].

The first complete sequence of a lung adenocarcinoma genome revealed about 50 000 single nucleotide variations in the tumor relative to normal lung [3]. This was followed by the sequencing study of a small-cell lung cancer genome which highlighted the role of tobacco carcinogens such as polycyclic aromatic hydrocarbons in shaping mutational patterns in lung cancers from smokers [4]. Transcriptome analysis of multiple lung adenocarcinoma patients using next-generation sequencing (NGS) recently showed the existence of a fusion gene containing the tyrosine kinase domain of the c-RET oncogene in 1%–2% of patients; this fusion leads to aberrant activation of RET kinase and is considered to be a new driver mutation of lung adenocarcinoma [5]. This finding was further confirmed through an independent study using a combination of targeted sequencing with an integrated molecular- and histopathology-based screening system [6]. Given that patients with c-RET fusions do not harbor mutations or fusions in EGFR, KRAS or ALK oncogenes, it is likely that c-RET fusion genes represent lung adenocarcinoma drivers and will lead to the definition of a new subclass of lung cancer [5].

Identifying mutations with high probabilities of being ‘drivers’, mutations that confer genes with oncogenic activity, is clearly a prototypical and certainly already a productive application of NGS, but the greater challenge is moving beyond the simple cataloging of mutations and establishing means for integrating diverse high-throughput data generated by NGS [7] to understand cancer at the multiple levels of gene networks and signaling pathways [8]. In this report, we describe a high-dimensional, high-throughput sequencing study of primary lung tumors and matched normal tissues isolated from 6 Korean female never-smoker patients with non-small cell lung cancer (NSCLC). This is the first multi-dimensional study of NSCLC that covers the exome-seq, RNA-seq, small RNA-seq, and methylated DNA immunoprecipitation-sequencing (MeDIP-seq). To complement the NGS data and obtain a full picture of sequence and structural variation, we also performed microarray-based gene expression profiling and array comparative genomic hybridization (array-CGH) study for DNA copy number variations (CNVs). Our study represents the simultaneous probing of the genome, transcriptome, and epigenome of biological samples revealing the full spectrum of cancer-associated alterations, including structural and genetic variations as well as changes in gene expression and epigenetic regulation. More importantly, we describe integrative analyses that entail the combination of the different types of omic data obtained in this study and that identify key regulators of NSCLC and elucidate relevant cellular processes at the systems level. Our findings show that the gene network modules that are highly relevant to the development of cancer, including those that govern progression through mitosis, are consistently disturbed in these NSCLC patients. We also report that multiple microRNAs are consistently inversely correlated with the predicted and validated target genes within these modules and throughout the whole network, indicating that microRNAs might be key regulatory agents of NSCLC development.


Multi-omic Data Description

To understand the genomic, transcriptomic and epigenomic changes in NSCLC, we performed high-throughput sequencing experiments for exome, transcriptome, and methylome on matched normal and tumor samples of 6 female non-smoker patients (see Figure S1 in File S8; data summary, experimental procedures are provided in the File S8; detailed sample/patient descriptions are provided in Table S1 in File S8 and File S1). CNV data were obtained from array-CGH assays. The genomic landscape of all NSCLC samples analyzed is visualized as a Circos plot of somatic mutations, transcriptome expression, CNVs, and structural variations_ENREF_6 (Figure 1; see Table S2 in File S8 for summary statistics of the exome data and Figure S2 in File S8 for Circos plots for individual patients) [9]. Experimental raw data and processed results are deposited in GEO (GSE37765) and SRA (SRA051952) databases. Raw additional data and information are also available at

Figure 1
Circos plot of somatic mutations, copy number variations, transcriptome expression, and structural variations.

NSCLC Somatic Mutations from Exome Sequencing Data

In our case, mutation calling by conventional programs such as Varscan (version 1.0) [10] did not show satisfactory performance, which was most likely due to the problem of normal cell contamination or heterogeneity of cancer cells. We therefore used the JointSNVmix program instead, to take advantage of the paired nature of samples (tumour and adjacent normal material) [11]. After validation by Sanger sequencing, we identified 47 somatic mutations that included 37 missense, 2 nonsense, and 7 silent mutations; there was also 1 mutation in the 3′ UTR (see Figure S3 in File S8). For several ambiguous cases, we subcloned PCR products and sequenced individual plasmid clones to confirm the mutation calls. Analyses of the validation process indicated that stringent criteria are required for the reliable prediction of somatic mutations if bulk clinical samples are used, as they were in our study. Cases with a predicted probability of over 0.999 often turned out to be false (45 positives and 55 negatives out of the 103 cases tested; PCR amplification failed in 3 cases). It should nevertheless be noted that some of these ‘false-positive’ somatic mutations may have occurred in a minority of tumor cells and are in fact positive, and newly discovered mutations in the future should be examined (e.g. for recurrence) with respect to all available raw data rather than just those confirmed by Sanger sequencing.

All confirmed mutations except for one were homozygous in the normal tissue. None of the somatic mutations identified in this study were identical to those reported in previous studies. In fact, none of the mutated genes isolated in this study except for CELF4 (G86C/S29T in the COSMIC database; G86A/S29N in patient 3) have been reported to be mutated in other studies or in the COSMIC database. A complete list of mutations is provided in Table S3 in File S8, and a summary of gene functions is provided in File S6. It should be noted that several genes identified are known to have functions that might be relevant to cancer development: BAZ1B regulates the DNA-damage response by phosphorylating the histone H2A.X; POLN is a DNA polymerase that performs translesion synthesis in response to DNA damage; and FBOX11, a component of the Skp1-Cullin1-F-box (SCF) complex, promotes neddylation of p53 and inhibits its transcriptional activity.

Differentially Expressed genes and Isoforms from RNA-Seq Data

We used the Bowtie and NEUMA applications for the mapping and quantification of RNA-Seq data, respectively [12], [13] (see Table S4 in File S8 for RNA-Seq data mapping summary). NEUMA, our in-house developed software (accessible at, provides a highly accurate estimation of transcript abundance both at the gene and individual splice variant (isoform) levels using an algorithm that mimics the real-time PCR process.

Determining differentially expressed genes (DEGs) and differentially expressed isoforms (DEIs) from RNA-Seq data was performed using the edgeR program, which supports the analysis of paired samples. A rigorous filtering procedure based on false discovery rates, minimum applicable patient numbers, and gene expression levels was devised to select reliable sets of DEGs and DEIs (see File S8 for details). For the final result, we obtained 1459 DEGs (543 upregulated and 916 downregulated) and 1320 DEIs (460 upregulated and 860 downregulated) in tumors when compared with normal tissues (see Table S5 in File S8). Imposing additional requirement of a minimum two-fold change yielded 387 DEGs (98 upregulated and 289 downregulated in tumors). The detailed procedure of the RNA-Seq analysis is described in the File S8, and the list of DEGs is provided in File S2.

Identification of Fusion Transcripts from RNA-Seq Data

We used FusionMap [14] and an in-house developed application, FusionScan, to predict fusion transcripts from RNA-Seq data. These two programs require the fusion boundary to be found inside one of the sequence reads, even in the case of paired-end data. The likelihood of missing fusion transcripts due to this requirement should be minimal since our RNA-Seq data have a high sequencing coverage (32.7X on average after mapping) and long read length (78 bp on average). Given that the two applications produced an overwhelmingly large number of candidates, we further filtered the initial output candidates by manual inspection of alignment against the hypothetical fusion transcripts. All candidate transcripts were examined for coherency of the 5′–3′ direction between the two fusion partner transcripts and strict adherence to the established wild-type exon-intron boundaries. Experimental validation was carried out by RT-PCR and Sanger sequencing. In total, we confirmed 19 fusion transcripts from 5 out of 6 patients as summarily presented in Table 1, including the MARK4-ERCC2 gene fusion shown in detail as an example (see Figure 2 and Table S6 in File S8).

Figure 2
MARK4-ERCC2 fusion transcript.
Table 1
List of experimentally confirmed fusion genes.

Most of the validated fusion transcripts were intra-chromosomal (15 out of 19) and presented tumor-specific events (16 out of 19; Table 1); the others were found in both tumor and normal cells. This finding implies that at least a subset of fusion transcripts likely arose from intra-chromosomal events such as copy number changes specific to cancer cells. Indeed, the fusion case of MARK4-ERCC2 showed the tumor-specific CNVs at the fusion points in the array-CGH data (Figure 2D). In-depth array-CGH analyses revealed 7 additional fusion events with strong association with DNA copy number changes (Table 1). Two fusions (RHPN2-PEPD and SIRT2-NPHS1) were observed in chromosome 19 of patient 1, and five cases were observed in chromosome 8 of patient 3 (PTK2-FAM84B, COX6C-LAPTM4B, STK3-PTK2, PKHD1L1-MATN2, NDUFB9-PGCP). Patient 3 harbored 7 fusion events in total two of which involved PTK2, also known as Focal Adhesion Kinase, a kinase with multiple functions including regulation of cell locomotion, survival and mitogen response [15]. In one patient (#6), all 3 fusions were detected in both normal and tumor samples, which suggests that these fusions are likely germline mutations. Of the 19 fusion events, 8 yielded in-frame gene fusions that potentially created proteins with novel functions.

Several of the gene fusion events could have an impact on cancer development and could potentially be driver mutations. For example, ERCC2 is involved in transcription-coupled DNA repair, and the tyrosine phosphatase PTPN12 is known to dephosphorylate and thereby inactivate the proto-oncogene c-ABL [16]. Perhaps most notably, one of the fusions involves CCDC6 and c-RET kinase, which is seen frequently in papillary thyroid carcinoma [17]. Preliminary structural analysis shows that the c-RET kinase domain is intact, raising the strong possibility that this fusion is a so-called ‘driver’ mutation of NSCLC. Recent publications have reported fusion events of c-RET kinase with KIF5B as well as with CCDC6 [5], [6], [18]. Finally, none of the identified fusions were found recurrently in our 6 patients, suggesting that a larger number of patients must be examined before the full significance of these fusion events can be evaluated.

Functional Interpretation of Somatic Mutations, DEGs, and Fusion Events

Although most of the somatic mutations and gene fusions are probably ‘passenger’ mutations, we cannot rule any of them out a priori as drivers of carcinogenesis. To facilitate the process of isolating functional DEGs and significant mutations, we performed a gene set analysis (GSA) and network analysis for 1536 genes, including the 47 genes that we found to have somatic mutations, the 37 genes involved in fusion, and the 1459 DEGs, some of which belong to more than one category.

A GSA, which tests the statistical enrichment or depletion of specific annotation terms, provides a comprehensive functional summary for genes in the list. We used the Ingenuity Pathway Analysis (IPA) software, which uses a database of evidence manually compiled from the literature. The most enriched term in the diseases and disorders category was cancer (p value = 2.13E-42), which supports the validity of our gene set. Other relevant terms in the molecular and cellular functions category included cellular growth and proliferation (p value = 1.71E-17) and cell death (p value = 1.97E-17). The IPA results are presented in Figure S9 in File S8. Gene ontology (GO) analysis produced similar results to IPA, albeit in a less comprehensive manner (data not shown).

We sought to gain further mechanistic and functional insights about the genes of interest using a network-based analysis. Protein-protein interaction data from the MIMI database was superimposed onto our gene list. The overall network thus obtained (Figure S10 in File S8) is complex but reveals a number of interesting interactions that may be connected to tumorigenesis. To identify network modules of coherent function, we used the MCODE application to find densely connected network components [19]. We found 8 network modules consisting of 66 genes in total (Figure S11 in File S8). Most genes in each network module showed expression changes in the same direction, suggesting a coherent and coordinated function in carcinogenesis as a gene network module. The largest network module features genes involved in mitotic cell-cycle regulation. In fact, this network module contains several key genes such as AURKB, PLK1, CCNE1, CCNB2, CHEK1, and PKMYT1, which are involved in the G2/M transition and/or M-phase progression [20], [21], [22]. All the genes within the cell-cycle regulation module were significantly upregulated (see below).

MicroRNA Analysis from Small RNA-Seq Data

Multiple studies have demonstrated that microRNAs could serve as viable tumor biomarkers and potential therapeutic targets or tools [23]. The computational pipeline for analyzing small RNA-Seq data (i.e., mapping, normalization, quantification, and identification of differentially expressed microRNAs) is illustrated in Figure S4 in File S8. On average, 70% and 76% of the total reads from small RNA-seq experiments were identified and mapped as microRNAs for normal and tumor tissues, respectively (see Table S7 in File S8).

Differentially expressed microRNAs (DEmiRs) were obtained using a process similar to that used to obtain DEGs. We limited the analysis to a subset of microRNAs belonging to the upper 25% in expression levels in at least one of the 12 samples. Other filtering conditions, such as fold change, were maintained. In total, we obtained 40 DEmiRs (23 up- and 17 downregulated in tumors compared to adjacent normal material; Table S8 in File S8 and File S3).

An inverse correlation in expression levels within a validated or predicted microRNA-mRNA pair provides strong evidence for an extant microRNA-target relationship in the biological context under examination. A search for inverse correlations between DEmiRs and DEGs yielded 151 relations (14 validated and 137 predicted) with a Pearson correlation cutoff of −0.5 and a p-value cutoff of 0.05 (see Table S9 and Figure S6 in File S8). Expanding the search scope beyond DEmiRs to identify other microRNAs of potential functional importance, we found 13 additional microRNAs with at least a two-fold change in expression between normal and tumor samples involved in 53 validated inverse correlations with DEGs. We did not use the predicted targets in this case, in order to avoid the inclusion of false positives. Overall, we identified 40 DEmiRs and 13 additional microRNAs that may play important roles in lung cancer development.

The MA-plot (logConcentration vs. logFoldChange) shows abundance and changes in expression and is thus an effective method for assessing the significance of potential biomarkers (see Figure 3 for the MA-plot of 40 DEmiRs and 13 additional microRNAs). Notably, the DEmiRs from a single genomic locus of chr7q32.2 (miR-96, miR-182, miR-183, and miR-183*) found upregulated in our study have also been previously reported as potential biomarkers for NSCLC [24]. By contrast, miR-144, miR-144* and miR-451 (all from chr17q11.2 locus) form a cluster of downregulated DEmiRs. Given that microRNAs derived from a polycistronic transcript often target the same set of genes, the change in expression of these microRNAs may have particularly strong effects. Among the 13 additional ‘non-DEmiR’ microRNAs, miR-21 (3.32 fold increase) and let-7b (2.20 fold decrease) deserve special attention. The miR-21 and let-7 families have previously been reported as oncomirs and tumor suppressors, respectively [25], [26]. Although neither family satisfied our stringent criteria for DEmiRs, their expression levels were among the highest, and their fold changes were statistically significant.

Figure 3
Differential expression of microRNAs.

Analysis of CNV and DNA Methylation Data

Copy number data from Agilent 1M CNV microarrays were analyzed using the within-slide lowess normalization and the circular binary segmentation (CBS) method [27]. We detected statistically significant somatic copy number alterations in tumor with the threshold of CNV change beyond log2 (tumor/normal) = ±0.3 (corresponding to the range outside 1.62–2.46 copies). We found multiple significant copy number gains in chromosomes 1, 5, 8, 16, 19, 20, and X, and copy number losses in chromosomes 1, 6, 9, 16, 17, 18, 19, and 22 (see Figure S7 and S8 in File S8 for CNV plots; Table S10 in File S8 for statistics). The genomic loci of copy number gains or losses in at least 3 patients included 8q22–24, 16p13, and 20q13 for gains and 1p12–13 and 9p21 for losses (File S4).

Comparison with previous studies (Table S11 in File S8) showed that our data agreed most closely with the Job et al. data that profiled 60 adenocarcinoma samples of never-smokers [28]. The genomic loci of 5p, 8q24, and 20q13 were repeatedly detected for copy number gains, and 9p21 was consistently detected for loss in copy numbers. Other markers may be never-smoker-female specific, which need to be validated with additional patient samples. Other chromosome-level events require further examination down the road. For example, the copy number gain in the large chr8q region, which has already been reported as a smoking-related biomarker for NSCLC patients [29], was detected in three of our six patients. Preliminary analysis indicated that many genes within this region were overexpressed in tumor samples (data not shown) and were involved in gene fusion in one of the patients (see Table 1).

DNA methylation patterns represent a potentially valuable biomarker in various types of cancer. We analyzed the DNA methylation data from MeDIP-seq using the Eland (version 2) and edgeR programs (see Table S14 in File S8 for mapping statistics; details in File S8). With conservative filtering options, we identified 558 differentially methylated regions (DMRs), almost 75% of which were in the promoter or 5′ UTR regions. The list of DMRs with genome annotation (from the UCSC genome browser database for hg19) is provided in File S5.

Copy number variations and DNA methylation are important factors in regulating gene expression. We investigated the correlation between CNV and mRNA expression in a similar manner to that performed in the microRNAs-target correlation study and obtained 107 positive correlations with a correlation cutoff of 0.5 (Table S12 in File S8). Highly correlated genes include several genes involved in gene fusion such as STK3, NDUFB9, COX6C, FAM84B, and PTK2. The inverse correlation analysis between mRNA expression and DNA methylation in promoter regions yielded 32 relations whose correlation coefficients were smaller than −0.3 (Table S13 in File S8).

Pathway Modeling with Network Modules and microRNAs

Several validated and predicted target genes of anti-correlated microRNAs were DEGs that were found in gene network modules using the MCODE clustering analysis described above. To obtain a more comprehensive picture of the regulatory networks, we integrated the 66 genes from MCODE clustering and microRNA-target relations using the IPA systems knowledgebase. The resulting gene network now incorporates the relevant microRNAs and additional genes interacting with the 66 DEGs (Figure 4). Genes were clustered initially into broad functional categories, and the inverse correlations with microRNAs were superimposed. Consistent with the initial clustering analysis, the cell-cycle cluster labeled as “cell death, cell cycle, cancer” formed the most comprehensive cluster. Furthermore, the sub-clustering analysis shows that the mitosis module, subsumed within the cell-cycle cluster, formed the most tightly organized gene network model (Figure 4).

Figure 4
NSCLC pathway modeling for female never-smokers.

The integrated network presented in this study enables the examination of gene- and microRNA-expression changes in combination with the interactions among them. Several microRNA-target relationships relevant to tumorigenesis are readily observed within this network. As described above, let-7b shows multiple functionally significant and inverse correlative interactions with genes within the network. The expression of two genes whose overexpression can drive tumorigenesis, AURKB and c-RET, is inversely correlated to the reduced expression of let-7b microRNA. AURKB (aurora kinase B) is the catalytic component of the chromosomal passenger complex that is responsible for chromatin condensation, bipolar spindle formation and attachment of chromatin to the bipolar spindle [21]. AURKB is also a DEG and a member of the network module described above whose members regulate the progression through the M phase. The c-RET receptor tyrosine kinase [17] has been found as a fusion protein in which the kinase domain remains intact. An overexpression of this fusion gene is therefore expected to have oncogenic activity, as recent studies have consistently suggested4. Two other cell cycle-related genes, NCAPG and UHRF1, are targeted by let-7b. NCAPG, the subunit G of non-SMC condensin I complex, which is required for conversion of interphase chromatin into condensed mitotic chromatin [30], is targeted by both let-7a and let-7b, exemplifying the regulation of a gene by multiple members of a microRNA cluster. UHRF1, a RING-finger type E3 ubiquitin ligase, which plays a major role in G1/S transition [31], is another target of let-7b. UHRF1 is also targeted by another microRNA, miR-1, which also regulates CENPF, a protein involved in kinetochore formation and chromosomal segregation. Whether the convergence of regulatory inputs into the genes from multiple microRNAs, as exemplified above, will lead to synergistic effects in the regulation of the mRNA level and/or in carcinogenesis remains an interesting question.

One of the most notable targets of the oncomir miR-21 is RECK, a membrane-anchored glycoprotein that inhibits matrix metalloproteinase-9. RECK, which is known to be strongly downregulated in multiple tumors and in cell lines transformed by oncogenes [32], is a predicted and inverse correlated target of several microRNAs (miR-96, miR-182, and miR-135b) within our network (Figure 4, Figure S6 in File S8). miR-96 and miR-182, which are DEmiRs from a single cluster, along with miR-183, have been proposed as RECK inhibitors, again exemplifying co-targeting by members of a single polycistronic cluster. Another microRNA-target relationship of significance involves MYBL2. This proto-oncogene produces the B-myb transcription factor, a well-known transcription factor with critical functions during G1/S and G2/M transitions. It is a validated target of miR-30a*, miR30b* and 30C-2*, all of which are downregulated in tumors consistent with upregulation of MYBL2.


In this study, we sought to define the nature of pathology at the molecular and systemic level within tumor cells of NSCLC female never-smoker patients, a group with few mutations known. Never-smokers account for 20% of men and 50% of all women with lung cancer and likely represent an etiologically distinct group from smokers [33]. Activating mutations of the EGFR tyrosine kinase are found at a far higher frequency among non-smokers, and subsequently a correspondingly higher efficacy is seen with EGFR tyrosine kinase inhibitors such as gefitinib and erlotinib in these patients [34], [35]. The Tumor Sequencing Project (TSP), which used a conventional Sanger sequencing technique to examine the coding exons of 623 candidate cancer genes in 188 lung adenocarcinomas, revealed that smokers suffer mutations at rates 5- to 10-times higher than never-smokers [36]. Given that most somatic mutations are expected to be “passengers,” the smaller number of mutations among never-smoker patients offers the possibility of more effectively isolating driver mutations.

We have identified and validated multiple point mutations and gene fusion events from 6 patients. It should be noted that these patients were negative mutations in KRAS and EGFR genes and did not sustain a fusion between EML4 and ALK which are the well-known transforming events of NSCLC. In fact, comparisons with other gene expression profiling studies of NSCLC cases indicate that the patient group tested in this study has acquired patterns in gene expression distinct from NSCLC cases with KRAS, EGFR or ALK mutations (File S7). A telling observation from our study was that none of the mutations or fusion genes was found in more than one patient. This strongly implies that for the large proportion of NSCLC cases negative for the established ‘driver’ mutations, the diversity or heterogeneity of mutations will likely remain a challenge for isolating biomarkers or drivers even in never-smokers. One of the most important goals for the future should be to expand the number of patient samples for NGS studies. This is critical for determining the recurrence and thereby the significance of each of the genetic aberrations.

Other than the CCDC6-RET fusion, which is expected to generate an oncogenic driver, it is difficult in most cases to predict the function of a given mutation (i.e. gain in oncogenic activity or a loss of tumor suppressor activity, or neither). This would be the case even when some of the information is available for the gene structure and function, unless an extensive characterization of the gene activity is performed in vitro and in vivo. Although such an analysis is obviously necessary, a different tactic based on integration of multiple omic data that investigates the interplay among genetic aberrations, transcripts, and regulatory factors in tumor cells can provide valuable information, as demonstrated in this study. Most importantly, we propose that the key cellular malfunction during tumorigenesis in our sampled population occurs in the control of the M-phase progression and that a set of specific microRNAs may be a source of viable biomarkers and functionally significant regulators of tumorigenesis. This will however require examination of individual genes in the context of lung cancer development as most of the genes promoting M-phase progression appear to up-regulated in diverse types of cancer as well according to our analysis on gene expression pattern using the GENT database (; data not shown).

In summary, integrative data analyses, such as those performed in this study, may be the only practically viable method for handling the anticipated volume of data from NGS studies. As more high-throughput data of multiple types from additional tumor samples are combined with reference datasets and gene networks, such as those reported here, the predictive power of integrative analyses should become more evident.

Supporting Information

File S1

Complete pathological reports of 6 patients.


File S2

List of differentially expressed genes and isoforms (DEGs and DEIs).


File S3

MicroRNA expression plot for tumor and normal tissues.


File S4

List of CNV genes in at least 3 out of 6 patients.


File S5

List of DMR genes.


File S6

Gene lists with brief summaries for somatic mutations and network modules.


File S7

Comparison with public expression data on NSCLC.


File S8

Supplemental information; Supplemental materials, methods, figures, tables and references.


Funding Statement

This work was funded by grants from the KRIBB Research Initiative Program (S.L.), National Research Foundation of Korea (grant numbers 20110019747 and 20110002321 to S.L. and grant number 20110014992 to WK), GIST Systems Biology Infrastructure Establishment Grant through ERCSB (SL and JK), and RP-Grant 2011 of Ewha Womans University (SC). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.


1. Meyerson M, Gabriel S, Getz G (2010) Advances in understanding cancer genomes through second-generation sequencing. Nature reviews Genetics 11: 685–696 [PubMed]
2. Jemal A, Siegel R, Ward E, Hao Y, Xu J, et al. (2009) Cancer statistics, 2009. CA: a cancer journal for clinicians 59: 225–249 [PubMed]
3. Lee W, Jiang Z, Liu J, Haverty PM, Guan Y, et al. (2010) The mutation spectrum revealed by paired genome sequences from a lung cancer patient. Nature 465: 473–477 [PubMed]
4. Pleasance ED, Stephens PJ, O'Meara S, McBride DJ, Meynert A, et al. (2010) A small-cell lung cancer genome with complex signatures of tobacco exposure. Nature 463: 184–190 [PMC free article] [PubMed]
5. Kohno T, Ichikawa H, Totoki Y, Yasuda K, Hiramoto M, et al. (2012) KIF5B-RET fusions in lung adenocarcinoma. Nature medicine 18: 375–377 [PubMed]
6. Takeuchi K, Soda M, Togashi Y, Suzuki R, Sakata S, et al. (2012) RET, ROS1 and ALK fusions in lung cancer. Nature medicine 18: 378–381 [PubMed]
7. Hawkins RD, Hon GC, Ren B (2010) Next-generation genomics: an integrative approach. Nature reviews Genetics 11: 476–486 [PMC free article] [PubMed]
8. Cho WC (2011) Science and technology: revealing genome networks for improving health Gastric Breast Cancer. 11: 58–64
9. Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, et al. (2009) Circos: an information aesthetic for comparative genomics. Genome research 19: 1639–1645 [PubMed]
10. Koboldt DC, Chen K, Wylie T, Larson DE, McLellan MD, et al. (2009) VarScan: variant detection in massively parallel sequencing of individual and pooled samples. Bioinformatics 25: 2283–2285 [PMC free article] [PubMed]
11. Roth A, Ding J, Morin R, Crisan A, Ha G, et al. (2012) JointSNVMix: a probabilistic model for accurate detection of somatic mutations in normal/tumour paired next-generation sequencing data. Bioinformatics 28: 907–913 [PMC free article] [PubMed]
12. Langmead B, Trapnell C, Pop M, Salzberg SL (2009) Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome biology 10: R25. [PMC free article] [PubMed]
13. Lee S, Seo CH, Lim B, Yang JO, Oh J, et al. (2011) Accurate quantification of transcriptome from RNA-Seq data by effective length normalization. Nucleic acids research 39: e9. [PMC free article] [PubMed]
14. Ge H, Liu K, Juan T, Fang F, Newman M, et al. (2011) FusionMap: detecting fusion genes from next-generation sequencing data at base-pair resolution. Bioinformatics 27: 1922–1928 [PubMed]
15. Gabarra-Niecko V, Schaller MD, Dunty JM (2003) FAK regulates biological processes important for the pathogenesis of cancer. Cancer metastasis reviews 22: 359–374 [PubMed]
16. Sun T, Aceto N, Meerbrey KL, Kessler JD, Zhou C, et al. (2011) Activation of multiple proto-oncogenic tyrosine kinases in breast cancer via loss of the PTPN12 phosphatase. Cell 144: 703–718 [PubMed]
17. Grieco M, Santoro M, Berlingieri MT, Melillo RM, Donghi R, et al. (1990) PTC is a novel rearranged form of the ret proto-oncogene and is frequently detected in vivo in human thyroid papillary carcinomas. Cell 60: 557–563 [PubMed]
18. Lipson D, Capelletti M, Yelensky R, Otto G, Parker A, et al. (2012) Identification of new ALK and RET gene fusions from colorectal and lung cancer biopsies. Nature medicine 18: 382–384 [PubMed]
19. Bader GD, Hogue CW (2003) An automated method for finding molecular complexes in large protein interaction networks. BMC bioinformatics 4: 2. [PMC free article] [PubMed]
20. Chan GK, Liu ST, Yen TJ (2005) Kinetochore structure and function. Trends in cell biology 15: 589–598 [PubMed]
21. Nigg EA (2001) Mitotic kinases as regulators of cell division and its checkpoints. Nature reviews Molecular cell biology 2: 21–32 [PubMed]
22. Smits VA, Klompmaker R, Arnaud L, Rijksen G, Nigg EA, et al. (2000) Polo-like kinase-1 is a target of the DNA damage checkpoint. Nature cell biology 2: 672–676 [PubMed]
23. Heneghan HM, Miller N, Kerin MJ (2010) MiRNAs as biomarkers and therapeutic targets in cancer. Current opinion in pharmacology 10: 543–550 [PubMed]
24. Zhu W, Liu X, He J, Chen D, Hunag Y, et al. (2011) Overexpression of members of the microRNA-183 family is a risk factor for lung cancer: a case control study. BMC cancer 11: 393. [PMC free article] [PubMed]
25. Pan X, Wang Z-X, Wang R (2010) MicroRNA-21: A novel therapeutic target in human cancer. Cancer Biology & Therapy 10: 1224–1232 [PubMed]
26. Yanaihara N, Caplen N, Bowman E, Seike M, Kumamoto K, et al. (2006) Unique microRNA molecular profiles in lung cancer diagnosis and prognosis. Cancer cell 9: 189–198 [PubMed]
27. Olshen AB, Venkatraman ES, Lucito R, Wigler M (2004) Circular binary segmentation for the analysis of array-based DNA copy number data. Biostatistics 5: 557–572 [PubMed]
28. Job B, Bernheim A, Beau-Faller M, Camilleri-Broet S, Girard P, et al. (2010) Genomic aberrations in lung adenocarcinoma in never smokers. PloS one 5: e15145. [PMC free article] [PubMed]
29. Huang YT, Lin X, Liu Y, Chirieac LR, McGovern R, et al. (2011) Cigarette smoking increases copy number alterations in nonsmall-cell lung cancer. Proceedings of the National Academy of Sciences of the United States of America 108: 16345–16350 [PubMed]
30. Geiman TM, Sankpal UT, Robertson AK, Chen Y, Mazumdar M, et al. (2004) Isolation and characterization of a novel DNA methyltransferase complex linking DNMT3B with components of the mitotic chromosome condensation machinery. Nucleic acids research 32: 2716–2729 [PMC free article] [PubMed]
31. Jeanblanc M, Mousli M, Hopfner R, Bathami K, Martinet N, et al. (2005) The retinoblastoma gene and its product are targeted by ICBP90: a key mechanism in the G1/S transition during the cell cycle. Oncogene 24: 7337–7345 [PubMed]
32. Loayza-Puch F, Yoshida Y, Matsuzaki T, Takahashi C, Kitayama H, et al. (2010) Hypoxia and RAS-signaling pathways converge on, and cooperatively downregulate, the RECK tumor-suppressor protein through microRNAs. Oncogene 29: 2638–2648 [PubMed]
33. Jemal A, Bray F (2011) Center MM, Ferlay J, Ward E, et al (2011) Global cancer statistics. CA: a cancer journal for clinicians 61: 69–90 [PubMed]
34. Marchetti A, Martella C, Felicioni L, Barassi F, Salvatore S, et al. (2005) EGFR mutations in non-small-cell lung cancer: analysis of a large series of cases and development of a rapid and sensitive method for diagnostic screening with potential implications on pharmacologic treatment. Journal of clinical oncology : official journal of the American Society of Clinical Oncology 23: 857–865 [PubMed]
35. Mok TS, Wu YL, Thongprasert S, Yang CH, Chu DT, et al. (2009) Gefitinib or carboplatin-paclitaxel in pulmonary adenocarcinoma. The New England journal of medicine 361: 947–957 [PubMed]
36. Ding L, Getz G, Wheeler DA, Mardis ER, McLellan MD, et al. (2008) Somatic mutations affect key pathways in lung adenocarcinoma. Nature 455: 1069–1075 [PMC free article] [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science