Amplified genes may be overexpressed, unchanged, or downregulated in cervical cancer cell lines

Several copy number-altered regions (CNAs) have been identified in the genome of cervical cancer, notably, amplifications of 3q and 5p. However, the contribution of copy-number alterations to cervical carcinogenesis is unresolved because genome-wide there exists a lack of correlation between copy-number alterations and gene expression. In this study, we investigated whether CNAs in the cell lines CaLo, CaSki, HeLa, and SiHa were associated with changes in gene expression. On average, 19.2% of the cell-line genomes had CNAs. However, only 2.4% comprised minimal recurrent regions (MRRs) common to all the cell lines. Whereas 3q had limited common gains (13%), 5p was entirely duplicated recurrently. Genome-wide, only 15.6% of genes located in CNAs changed gene expression; in contrast, the rate in MRRs was up to 3 times this. Chr 5p was confirmed entirely amplified by FISH; however, maximum 33.5% of the explored genes in 5p were deregulated. In 3q, this rate was 13.4%. Even in 3q26, which had 5 MRRs and 38.7% recurrently gained SNPs, the rate was only 15.1%. Interestingly, up to 19% of deregulated genes in 5p and 73% in 3q26 were downregulated, suggesting additional factors were involved in gene repression. The deregulated genes in 3q and 5p occurred in clusters, suggesting local chromatin factors may also influence gene expression. In regions amplified discontinuously, downregulated genes increased steadily as the number of amplified SNPs increased (p<0.01, Spearman's correlation). Therefore, partial gene amplification may function in silencing gene expression. Additional genes in 1q, 3q and 5p could be involved in cervical carcinogenesis, specifically in apoptosis. These include PARP1 in 1q, TNFSF10 and ECT2 in 3q and CLPTM1L, AHRR, PDCD6, and DAP in 5p. Overall, gene expression and copy-number profiles reveal factors other than gene dosage, like epigenetic or chromatin domains, may influence gene expression within the entirely amplified genome segments.     

Impact of Gene Dosage on Gene Expression,Biological Processes and Survival in Cervical Cancer: A Genome-Wide Follow-Up Study

We investigated the role of tumor copy number (CN)–altered genome (CN-AG) in the carcinogenesis of cervical cancer (CC), especially its effect on gene expression, biological processes, and patient survival. Fifty-nine human papillomavirus 16 (HPV16)-positive CCs were investigated with microarrays–31 for mapping CN-AG and 55 for global gene expression, with 27 CCs in common. Five-year survival was investigated in 55 patients. Deletions and amplifications >2.5 Mb were defined as CN alterations. The %CN-AG varied from 0 to 32.2% (mean = 8.1±8.9). Tumors were classified as low (mean = 0.5±0.6, n = 11), medium (mean = 5.4±2.4, n = 10), or high (mean = 19.2±6.6, n = 10) CN. The highest %CN-AG was found in 3q, which contributed an average of 55% of all CN alterations. Genome-wide, only 5.3% of CN-altered genes were deregulated directly by gene dosage. In contrast, the rate in fully duplicated 3q was twice as high. Amplification of 3q explained 23.2% of deregulated genes in whole tumors (r² = 0.232, p = 0.006; analysis of variance), including genes located in 3q and other chromosomes. A total of 862 genes were deregulated exclusively in high-CN tumors, but only 22.9% were CN altered. This suggests that the remaining genes are not deregulated directly by gene dosage, but by mechanisms induced in trans by CN-altered genes. Anaphase-promoting complex/cyclosome (APC/C)-dependent proteasome proteolysis, glycolysis, and apoptosis were upregulated, whereas cell adhesion and angiogenesis were downregulated exclusively in high-CN tumors. The high %CN-AG and upregulated gene expression profile of APC/C-dependent proteasome proteolysis were associated with poor patient survival (p<0.05, log-rank test). Along with glycolysis, they were linearly associated with FIGO stage (r>0.38, p<0.01, Spearman test). Therefore, inhibition of APC/C-dependent proteasome proteolysis and glycolysis could be useful for CC treatment. However, whether they are indispensable for tumor growth remains to be demonstrated.     

Integration of DNA copy number alterations and transcriptional expression analysis in human gastric cancer

Background

Genomic instability with frequent DNA copy number alterations is one of the key hallmarks of carcinogenesis. The chromosomal regions with frequent DNA copy number gain and loss in human gastric cancer are still poorly defined. It remains unknown how the DNA copy number variations contributes to the changes of gene expression profiles, especially on the global level.

Principal Findings

We analyzed DNA copy number alterations in 64 human gastric cancer samples and 8 gastric cancer cell lines using bacterial artificial chromosome (BAC) arrays based comparative genomic hybridization (aCGH). Statistical analysis was applied to correlate previously published gene expression data obtained from cDNA microarrays with corresponding DNA copy number variation data to identify candidate oncogenes and tumor suppressor genes. We found that gastric cancer samples showed recurrent DNA copy number variations, including gains at 5p, 8q, 20p, 20q, and losses at 4q, 9p, 18q, 21q. The most frequent regions of amplification were 20q12 (7/72), 20q12–20q13.1 (12/72), 20q13.1–20q13.2 (11/72) and 20q13.2–20q13.3 (6/72). The most frequent deleted region was 9p21 (8/72). Correlating gene expression array data with aCGH identified 321 candidate oncogenes, which were overexpressed and showed frequent DNA copy number gains; and 12 candidate tumor suppressor genes which were down-regulated and showed frequent DNA copy number losses in human gastric cancers. Three networks of significantly expressed genes in gastric cancer samples were identified by ingenuity pathway analysis.

Conclusions

This study provides insight into DNA copy number variations and their contribution to altered gene expression profiles during human gastric cancer development. It provides novel candidate driver oncogenes or tumor suppressor genes for human gastric cancer, useful pathway maps for the future understanding of the molecular pathogenesis of this malignancy, and the construction of new therapeutic targets.     

Copy Number Gains at 8q24 and 20q11-q13 in Gastric Cancer Are More Common in Intestinal-Type than Diffuse-Type

The present study was aimed at discovering DNA copy number alterations (CNAs) involved in the carcinogenesis of stomach and at understanding their clinicopathological significances in the Korean population. DNA copy numbers were analyzed using Agilent 244K or 400K array comparative genomic hybridization (aCGH) in fresh-frozen tumor and matched normal tissues from 40 gastric cancer patients. Some of the detected CNA regions were validated using multiplex ligation-dependent probe amplification (MLPA) in six of the 40 patients and customized Agilent 60K aCGH in an independent set of 48 gastric cancers. The mRNA levels of genes at common CNA regions were analyzed using quantitative real-time PCR. Copy number gains were more common than losses across the entire genome in tumor tissues compared to matched normal tissues. The mean number of alterations per case was 64 for gains and 40 for losses, and the median aberration length was 44016 bp for gains and 4732 bp for losses. Copy number gains were frequently detected at 7p22.1 (20%), 8q24.21 (27%–30%), 8q24.3 (22%–48%), 13q34 (20%–31%), and 20q11-q13 (25%–30%), and losses at 3p14.2 (43%), 4q35.2 (27%), 6q26 (23%), and 17p13.3 (20%–23%). CNAs at 7p22.1, 13q34, and 17p13.3 have not been reported in other populations. Most of the copy number losses were associated with down-regulation of mRNA levels, but the correlation between copy number gains and mRNA expression levels varied in a gene-dependent manner. In addition, copy number gains tended to occur more commonly in intestinal-type cancers than in diffuse-type cancers. In conclusion, the present study suggests that copy number gains at 8q24 and 20q11-q13 and losses at 3p14.2 may be common events in gastric cancer but CNAs at 7p22.1, 13q34, and 17p13.3 may be Korean-specific.     

Genomic profiling of advanced-stage oral cancers reveals chromosome 11q alterations as markers of poor clinical outcome

Identifying oral cancer lesions associated with high risk of relapse and predicting clinical outcome remain challenging questions in clinical practice. Genomic alterations may add prognostic information and indicate biological aggressiveness thereby emphasizing the need for genome-wide profiling of oral cancers. High-resolution array comparative genomic hybridization was performed to delineate the genomic alterations in clinically annotated primary gingivo-buccal complex and tongue cancers (n = 60). The specific genomic alterations so identified were evaluated for their potential clinical relevance. Copy-number changes were observed on chromosomal arms with most frequent gains on 3q (60%), 5p (50%), 7p (50%), 8q (73%), 11q13 (47%), 14q11.2 (47%), and 19p13.3 (58%) and losses on 3p14.2 (55%) and 8p (83%). Univariate statistical analysis with correction for multiple testing revealed chromosomal gain of region 11q22.1–q22.2 and losses of 17p13.3 and 11q23–q25 to be associated with loco-regional recurrence (P = 0.004, P = 0.003, and P = 0.0003) and shorter survival (P = 0.009, P = 0.003, and P 0.0001) respectively. The gain of 11q22 and loss of 11q23-q25 were validated by interphase fluorescent in situ hybridization (I-FISH). This study identifies a tractable number of genomic alterations with few underlying genes that may potentially be utilized as biological markers for prognosis and treatment decisions in oral cancers.     

Adipose tissue gene expression analysis reveals changes in inflammatory, mitochondrial respiratory and lipid metabolic pathways in obese insulin-resistant subjects

Background

To elucidate gene expression associated with copy number changes, we performed a genome-wide copy number and expression microarray analysis of 25 pairs of gastric tissues.

Methods

We applied laser capture microdissection (LCM) to obtain samples for microarray experiments and profiled DNA copy number and gene expression using 244K CGH Microarray and Human Exon 1.0 ST Microarray.

Results

Obviously, gain at 8q was detected at the highest frequency (70%) and 20q at the second (63%). We also identified molecular genetic divergences for different TNM-stages or histological subtypes of gastric cancers. Interestingly, the C20orf11 amplification and gain at 20q13.33 almost separated moderately differentiated (MD) gastric cancers from poorly differentiated (PD) type. A set of 163 genes showing the correlations between gene copy number and expression was selected and the identified genes were able to discriminate matched adjacent noncancerous samples from gastric cancer samples in an unsupervised two-way hierarchical clustering. Quantitative RT-PCR analysis for 4 genes (C20orf11, XPO5, PUF60, and PLOD3) of the 163 genes validated the microarray results. Notably, some candidate genes (MCM4 and YWHAZ) and its adjacent genes such as PRKDC, UBE2V2, ANKRD46, ZNF706, and GRHL2, were concordantly deregulated by genomic aberrations.

Conclusions

Taken together, our results reveal diverse chromosomal region alterations for different TNM-stages or histological subtypes of gastric cancers, which is helpful in researching clinicopathological classification, and highlight several interesting genes as potential biomarkers for gastric cancer.     

Heterogeneity of 11q13 region rearrangements in laryngeal squamous cell carcinoma analyzed by microarray platforms and fluorescence in situ hybridization

We reinvestigated rearrangements occurring in region q13 of chromosome 11 aiming to: (i) describe heterogeneity of the observed structural alterations, (ii) estimate amplicon size and (iii) identify of oncogenes involved in laryngeal cancer progression as potential targets for therapy. The study included 17 cell lines derived from laryngeal cancers and 34 specimens from primary laryngeal tumors. The region 11q13 was analyzed by fluorescence in situ hybridization (FISH), array comparative genomic hybridization (aCGH) and gene expression microarray. Next, quantitative real time PCR was used for chosen genes to confirm results from aCGH and gene expression microarray. The observed pattern of aberrations allows to distinguish three ways, in which gain and amplification involving 11q13 region may occur: formation of a homogeneously staining region; breakpoints in/near 11q13, which lead to the three to sevenfold increase of the copy number of 11q13 region; the presence of additional copies of the whole chromosome 11. The minimal altered region of gain and/or amplification was limited to ~1.8 Mb (chr.11:69,395,184–71,209,568) and comprised mostly 11q13.3 band which contain 12 genes. Five, out of these genes (CCND1, ORAOV1, FADD, PPFIA1, CTTN) had higher expression levels in comparison to healthy controls. Apart from CCND1 gene, which has an established role in pathogenesis of head and neck cancers, CTTN, ORAOV1 and FADD genes appear to be oncogene-candidates in laryngeal cancers, while a function of PPFIA1 requires further studies.     

Molecular characterization of acquired tolerance of tumor cells to picropodophyllin (PPP)

Background

Picropodophyllin (PPP) is a promising novel anti-neoplastic agent that efficiently kills tumor cells in vitro and causes tumor regression and increased survival in vivo. We have previously reported that PPP treatment induced moderate tolerance in two out of 10 cell lines only, and here report the acquired genomic and expression alterations associated with PPP selection over 1.5 years of treatment.

Methodology/Principal Findings

Copy number alterations monitored using metaphase and array-based comparative genomic hybridization analyses revealed largely overlapping alterations in parental and maximally tolerant cells. Gain/ amplification of the MYC and PVT1 loci in 8q24.21 were verified on the chromosome level. Abnormalities observed in connection to PPP treatment included regular gains and losses, as well as homozygous losses in 10q24.1-q24.2 and 12p12.3-p13.2 in one of the lines and amplification at 5q11.2 in the other. Abnormalities observed in both tolerant derivatives include amplification/gain of 5q11.2, gain of 11q12.1-q14.3 and gain of 13q33.3-qter. Using Nexus software analysis we combined the array-CGH data with data from gene expression profilings and identified genes that were altered in both inputs. A subset of genes identified as downregulated (ALDH1A3, ANXA1, TLR4 and RAB5A) or upregulated (COX6A1, NFIX, ME1, MAPK and TAP2) were validated by siRNA in the tolerant or parental cells to alter sensitivity to PPP and confirmed to alter sensitivity to PPP in further cell lines.

Conclusions

Long-term PPP selection lead to altered gene expression in PPP tolerant cells with increase as well as decrease of genes involved in cell death such as PTEN and BCL2. In addition, acquired genomic copy number alterations were observed that were often reflected by altered mRNA expression levels for genes in the same regions.     

Candidate Luminal B Breast Cancer Genes Identified by Genome,Gene Expression and DNA Methylation Profiling

Breast cancers (BCs) of the luminal B subtype are estrogen receptor-positive (ER+), highly proliferative, resistant to standard therapies and have a poor prognosis. To better understand this subtype we compared DNA copy number aberrations (CNAs), DNA promoter methylation, gene expression profiles, and somatic mutations in nine selected genes, in 32 luminal B tumors with those observed in 156 BCs of the other molecular subtypes. Frequent CNAs included 8p11-p12 and 11q13.1-q13.2 amplifications, 7q11.22-q34, 8q21.12-q24.23, 12p12.3-p13.1, 12q13.11-q24.11, 14q21.1-q23.1, 17q11.1-q25.1, 20q11.23-q13.33 gains and 6q14.1-q24.2, 9p21.3-p24,3, 9q21.2, 18p11.31-p11.32 losses. A total of 237 and 101 luminal B-specific candidate oncogenes and tumor suppressor genes (TSGs) presented a deregulated expression in relation with their CNAs, including 11 genes previously reported associated with endocrine resistance. Interestingly, 88% of the potential TSGs are located within chromosome arm 6q, and seven candidate oncogenes are potential therapeutic targets. A total of 100 candidate oncogenes were validated in a public series of 5,765 BCs and the overexpression of 67 of these was associated with poor survival in luminal tumors. Twenty-four genes presented a deregulated expression in relation with a high DNA methylation level. FOXO3, PIK3CA and TP53 were the most frequent mutated genes among the nine tested. In a meta-analysis of next-generation sequencing data in 875 BCs, KCNB2 mutations were associated with luminal B cases while candidate TSGs MDN1 (6q15) and UTRN (6q24), were mutated in this subtype. In conclusion, we have reported luminal B candidate genes that may play a role in the development and/or hormone resistance of this aggressive subtype.     

Synergistic Interactions between Drosophila Orthologues of Genes Spanned by De Novo Human CNVs Support Multiple-Hit Models of Autism

Autism spectrum disorders (ASDs) are highly heritable and characterised by deficits in social interaction and communication, as well as restricted and repetitive behaviours. Although a number of highly penetrant ASD gene variants have been identified, there is growing evidence to support a causal role for combinatorial effects arising from the contributions of multiple loci. By examining synaptic and circadian neurological phenotypes resulting from the dosage variants of unique human:fly orthologues in Drosophila, we observe numerous synergistic interactions between pairs of informatically-identified candidate genes whose orthologues are jointly affected by large de novo copy number variants (CNVs). These CNVs were found in the genomes of individuals with autism, including a patient carrying a 22q11.2 deletion. We first demonstrate that dosage alterations of the unique Drosophila orthologues of candidate genes from de novo CNVs that harbour only a single candidate gene display neurological defects similar to those previously reported in Drosophila models of ASD-associated variants. We then considered pairwise dosage changes within the set of orthologues of candidate genes that were affected by the same single human de novo CNV. For three of four CNVs with complete orthologous relationships, we observed significant synergistic effects following the simultaneous dosage change of gene pairs drawn from a single CNV. The phenotypic variation observed at the Drosophila synapse that results from these interacting genetic variants supports a concordant phenotypic outcome across all interacting gene pairs following the direction of human gene copy number change. We observe both specificity and transitivity between interactors, both within and between CNV candidate gene sets, supporting shared and distinct genetic aetiologies. We then show that different interactions affect divergent synaptic processes, demonstrating distinct molecular aetiologies. Our study illustrates mechanisms through which synergistic effects resulting from large structural variation can contribute to human disease.     

Expression of Alpha-type Platelet-derived Growth Factor Receptor–influenced Genes Predicts Clinical Outcome in Glioma

BACKGROUND: Alpha-type platelet-derived growth factor receptor (PDGFRα) is a cell surface tyrosine kinase receptor for members of the platelet-derived growth factor family. PDGFRα plays an important role in the regulation of several biological processes and contributes to the pathophysiology of a broad range of human cancers, including glioma. Here, we hypothesize that the genes directly or indirectly influenced by PDGFRα might be useful for prognosis in glioma. METHODS: By comparing the genome-wide gene expression pattern between PDGFRα⁺ and PDGFRα⁻ cells from human oligodendrocyte progenitor, we defined the genes potentially influenced by PDGFRα. RESULTS: The PDGFRα-influenced genes are strongly associated with cancer-related pathways. We subsequently developed a prognostic gene signature derived from the PDGFRα-influenced genes. This gene signature is able to predict clinical outcome of glioma. This signature is also independent of traditional prognostic factors of glioma. Resampling tests indicate that the prognostic power of this gene signature outperforms random gene sets selected from human genome. More importantly, this signature is superior to the random gene signatures selected from glioma related genes. CONCLUSIONS: Despite the absence of clear elucidation of molecular mechanisms, this study suggests the vital role of PDGFRα in carcinogenesis. Furthermore, the PDGFRα-based gene signature provides a promising prognostic tool for glioma and validates PDGFRα as a novel and effective therapeutic target in human cancers.     

The Landscape of Candidate Driver Genes Differs between Male and Female Breast Cancer

The rapidly growing collection of diverse genome-scale data from multiple tumor types sheds light on various aspects of the underlying tumor biology. With the objective to identify genes of importance for breast tumorigenesis in men and to enable comparisons with genes important for breast cancer development in women, we applied the computational framework COpy Number and EXpression In Cancer (CONEXIC) to detect candidate driver genes among all altered passenger genes. Unique to this approach is that each driver gene is associated with several gene modules that are believed to be altered by the driver. Thirty candidate drivers were found in the male breast cancers and 67 in the female breast cancers. We identified many known drivers of breast cancer and other types of cancer, in the female dataset (e.g. GATA3, CCNE1, GRB7, CDK4). In contrast, only three known cancer genes were found among male breast cancers; MAP2K4, LHP, and ZNF217. Many of the candidate drivers identified are known to be involved in processes associated with tumorigenesis, including proliferation, invasion and differentiation. One of the modules identified in male breast cancer was regulated by THY1, a gene involved in invasion and related to epithelial-mesenchymal transition. Furthermore, men with THY1 positive breast cancers had significantly inferior survival. THY1 may thus be a promising novel prognostic marker for male breast cancer. Another module identified among male breast cancers, regulated by SPAG5, was closely associated with proliferation. Our data indicate that male and female breast cancers display highly different landscapes of candidate driver genes, as only a few genes were found in common between the two. Consequently, the pathobiology of male breast cancer may differ from that of female breast cancer and can be associated with differences in prognosis; men diagnosed with breast cancer may consequently require different management and treatment strategies than women.     

A transcriptomic computational analysis of mastic oil-treated Lewis lung carcinomas reveals molecular mechanisms targeting tumor cell growth and survival

Background

DNA copy number alterations are frequently observed in ovarian cancer, but it remains a challenge to identify the most relevant alterations and the specific causal genes in those regions.

Methods

We obtained high-resolution 500K SNP array data for 52 ovarian tumors and identified the most statistically significant minimal genomic regions with the most prevalent and highest-level copy number alterations (recurrent CNAs). Within a region of recurrent CNA, comparison of expression levels in tumors with a given CNA to tumors lacking that CNA and to whole normal ovary samples was used to select genes with CNA-specific expression patterns. A public expression array data set of laser capture micro-dissected (LCM) non-malignant fallopian tube epithelia and LCM ovarian serous adenocarcinoma was used to evaluate the effect of cell-type mixture biases.

Results

Fourteen recurrent deletions were detected on chromosomes 4, 6, 9, 12, 13, 15, 16, 17, 18, 22 and most prevalently on X and 8. Copy number and expression data suggest several apoptosis mediators as candidate drivers of the 8p deletions. Sixteen recurrent gains were identified on chromosomes 1, 2, 3, 5, 8, 10, 12, 15, 17, 19, and 20, with the most prevalent gains localized to 8q and 3q. Within the 8q amplicon, PVT1, but not MYC, was strongly over-expressed relative to tumors lacking this CNA and showed over-expression relative to normal ovary. Likewise, the cell polarity regulators PRKCI and ECT2 were identified as putative drivers of two distinct amplicons on 3q. Co-occurrence analyses suggested potential synergistic or antagonistic relationships between recurrent CNAs. Genes within regions of recurrent CNA showed an enrichment of Cancer Census genes, particularly when filtered for CNA-specific expression.

Conclusion

These analyses provide detailed views of ovarian cancer genomic changes and highlight the benefits of using multiple reference sample types for the evaluation of CNA-specific expression changes.     

Transcriptome Analysis of the Hippocampus in Novel Rat Model of Febrile Seizures

Febrile seizures (FS) are the most common type of convulsive events in infants and young children, but the precise underlying genetic mechanism remains to be explored. To investigate the underlying pathogenic factors in FS and subsequent epilepsy, alterations in gene expression between the two new strains of rats (hyperthermia-prone [HP] vs hyperthermia-resistant [HR]), were investigated by using the Whole Rat Genome Oligo Microarray. This process identified 1,140 differentially expressed genes (DEGs; 602 upregulated and 538 downregulated), which were analyzed to determine significant Gene Ontology (GO) categories, signaling pathways and gene networks. Based on the GO analyses, the modified genes are closely related to various FS pathogenesis factors, including immune and inflammatory responses and ion transport. Certain DEGs identified have not been previously examined in relation to FS pathogenesis. Among these genes is dipeptidyl peptidase 4 (DPP4), a gene closely linked to interleukin 6 (IL-6), which played a key role in the gene network analysis. Furthermore, sitagliptin, a DPP4 inhibitor significantly decreased epileptic discharge in rats, observed via electroencephalogram, suggesting an important role for DPP4 in FS. The effectiveness of sitagliptin in reducing seizure activity may occur through a mechanism that stabilizes cellular Ca²⁺ homeostasis. In addition, DPP4 expression may be regulated by DNA methylation. The hippocampal gene expression profiles in novel rat models of FS provides a large database of candidate genes and pathways, which will be useful for researchers interested in disorders of neuronal excitability.     

Alterations in p16 and p53 genes and chromosomal findings in patients with lung cancer: Fluorescence in situ hybridization and cytogenetic studies

Background: Chromosomal aberrations and instability of gene(s) are two factors related to the genetic instability of cancer cells. A loss of the tumor-suppressor function of the genes p16 and p53 is the most common event leading to the development of human cancers. Carcinoma of the lung is the leading cause of cancer deaths in the world. Chromosomal abnormalities in lung cancer may provide a valuable clue to the identification of target loci and culminate in a successful search for the major genes. The aim of this study was to investigate (i) alterations of the p16 and p53 genes and (ii) chromosomal aberrations in patients with small cell and non-small cell lung cancer by fluorescence in situ hybridization (FISH) and cytogenetic studies. Methods: We carried out cytogenetic analysis by Giemsa-banding in 18 cases. FISH probes for the p16 and p53 genes were also used on interphase nuclei to screen the alterations in these genes in lung cancer (LC).Results: We observed a high frequency of losses of the p16 – in 8/18 (44%) – and p53 – in 7/18 (39%) – genes in the cases with LC. A total of 18 patients showed predominantly numerical and structural aberrations. Among these two types, structural aberrations predominated and usually consisted of deletions, breaks, and fragilities in various chromosomes. Both structural and numerical changes were observed in almost all patients. Chromosomes 3 and 1 were found to be most frequently involved in structural abnormalities, followed by chromosomes 6, 9, and 8. Autosomal aneuploidies were also observed to be the most frequent (chromosomes 22, 19, 18, 20, 9, and 17), followed by those of the X and Y chromosomes. The expression of fragile sites was also found to be significantly higher in seven chromosomal regions: 3p14, 1q21, 1q12, 6q26, 9q13, 8q22, and 8q24. Conclusion: Our data confirmed that DNA damage and genomic instability may be factors contributing to the mutation profile and development of lung cancer. The patients who developed lung cancer showed a high frequency of loss of both p16 and p53, in addition to chromosomal aberrations. Tobacco could be a major carcinogenic factor in lung-cancer progression. The loss of p16 and p53, and increased incidence of autosomal aneuploidy and chromatid breaks, along with other chromosomal alterations, can contribute to the progression of the disease.