Oncoprotein ZNF322A transcriptionally deregulates alpha-adducin, cyclin D1 and p53 to promote tumor growth and metastasis in lung cancer

ZNF322A encoding a classical Cys2His2 zinc finger transcription factor was previously revealed as a potential oncogene in lung cancer patients. However, the oncogenic role of ZNF322A and its underlying mechanism in lung tumorigenesis remain elusive. Here we show ZNF322A protein overexpression in 123 Asian and 74 Caucasian lung cancer patients. Multivariate Cox regression analysis indicated that ZNF322A was an independent risk factor for a poor outcome in lung cancer, corroborating the Kaplan–Meier results that patients with ZNF322A protein overexpression had significantly poorer overall survival than other patients. Overexpression of ZNF322A promoted cell proliferation and soft agar growth by prolonging cell cycle in S phase in multiple lung cell lines, including the immortalized lung cell BEAS-2B. In addition, ZNF322A overexpression enhanced cell migration and invasion, whereas knockdown of ZNF322A reduced cell growth, invasion and metastasis abilities in vitro and in vivo. Quantitative proteomic analysis revealed potential ZNF322A-regulated downstream targets, including alpha-adducin (ADD1), cyclin D1 (CCND1), and p53. Using luciferase promoter activity assay combined with site-directed mutagenesis and sequential chromatin immunoprecipitation-PCR assay, we found that ZNF322A could form a complex with c-Jun and cooperatively activate ADD1 and CCND1 but repress p53 gene transcription by recruiting differential chromatin modifiers, such as histone deacetylase 3, in an AP-1 element dependent manner. Reconstitution experiments indicated that CCND1 and p53 were important to ZNF322A-mediated promotion of cell proliferation, whereas ADD1 was necessary for ZNF322A-mediated cell migration and invasion. Our results provide compelling evidence that ZNF322A overexpression transcriptionally dysregulates genes involved in cell growth and motility therefore contributes to lung tumorigenesis and poor prognosis.

Gene expression patterns in CD4+ peripheral blood cells in healthy subjects and stage IV melanoma patients

Melanoma patients exhibit changes in immune responsiveness in the local tumor environment, draining lymph nodes, and peripheral blood. Immune-targeting therapies are revolutionizing melanoma patient care increasingly, and studies show that patients derive clinical benefit from these newer agents. Nonetheless, predicting which patients will benefit from these costly therapies remains a challenge. In an effort to capture individual differences in immune responsiveness, we are analyzing patterns of gene expression in human peripheral blood cells using RNAseq. Focusing on CD4+ peripheral blood cells, we describe multiple categories of immune regulating genes, which are expressed in highly ordered patterns shared by cohorts of healthy subjects and stage IV melanoma patients. Despite displaying conservation in overall transcriptome structure, CD4+ peripheral blood cells from melanoma patients differ quantitatively from healthy subjects in the expression of more than 2000 genes. Moreover, 1300 differentially expressed genes are found in transcript response patterns following activation of CD4+ cells ex vivo, suggesting that widespread functional discrepancies differentiate the immune systems of healthy subjects and melanoma patients. While our analysis reveals that the transcriptome architecture characteristic of healthy subjects is maintained in cancer patients, the genes expressed differentially among individuals and across cohorts provide opportunities for understanding variable immune states as well as response potentials, thus establishing a foundation for predicting individual responses to stimuli such as immunotherapeutic agents.

Abstract 3895: Genomic alterations modify ubiquitination pathways in pulmonary carcinoid tumors

Introduction: Pulmonary carcinoid tumors account for up to 5% of all lung malignancies in adults, and comprise 30% of all carcinoid malignancies. They are defined histologically as typical carcinoid (TC) and atypical carcinoid (AC) tumors, and are characterized by neuroendocrine differentiation and the potential to metastasize.

Method: We characterized genomic alterations in pulmonary carcinoid tumors using whole genome and exome sequencing in addition to mRNA expression and SNP genotyping from specimens of normal lung, typical and atypical carcinoid, and SCLC.

Results: Analysis of the sequencing data identified novel mutations in the pulmonary neuroendocrine tumor spectrum, including KRAS, ABCA9, GALNT10, HSFX1, PTHLH, RANBP2, BRIP1, FXR1, BACE2, NCOR2, DGKI, HERC1, HERC5 and GLI3. Pathway analysis found 45% of the cases to harbor a mutation in genes that are involved in protein-protein interaction with ubiquitin c (UBC) and may therefore play a role in ubiquitin protein degradation pathway. We also identified genes differentially expressed in corresponding regions of copy number variation, 52 of which were linked to ubiquitination pathways in 88% of the samples genotyped. Genes amplified in 15 to 23% of the cases included ITM2A, MAOA, EGR2, HAGH, TSC22D3, NBEA, TCEAL2, and TSPYL2 whereas deleted regions in 15% of the cases implicated CFD, GDI1, IFITM1 and IFITM2.

Conclusion: These results suggest that in addition to previously identified pathways, dysregulation of the ubiquitin pathway may contribute to lung carcinoid tumorigenesis and could be exploited for therapeutic purposes.

Citation Format: Michael K. Asiedu, Charles L. Thomas, Jie Dong, Prasidda Khadka, Zhifu Sun, Farhad Kosari, Jin Jen, Julian Molina, George Vasmatzis, Marie Christine Aubry, Ping Yang, Dennis Wigle. Genomic alterations modify ubiquitination pathways in pulmonary carcinoid tumors. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 3895. doi:10.1158/1538-7445.AM2015-3895

CHFR methylation strongly correlates with methylation of DNA damage repair and apoptotic pathway genes in non-small cell lung cancer

DNA methylation occurs commonly in non-small cell lung cancer (NSCLC). We sought to determine the frequency and relationship of methylation of key genes involved in the pathways of mitotic checkpoint control, DNA damage repair, apoptosis, and growth factor signaling in these patients. We analyzed the DNA methylation status of eight genes (CHFR, FANCF, MGMT, p16, DAPK, ASC or TMS-1, RAR-B, and CRBP1) using nested methylation-specific PCR (MSP) on over 314 paraffin-embedded, human non-small cell lung cancer samples. We determined the methylation frequency of each gene in addition to the association of the methylation of each gene with other members of the panel. Methylation was a common event in these samples. Our methylation analysis showed frequencies of methylation of 10% for CHFR, 14% for FANCF, 30% for MGMT, 29% for p16, 17% for DAPK, 33% for ASC, 38% for RAR-B, and 7% for CRBP1. There was a strong correlation between methylation of the mitotic G2-M checkpoint gene, CHFR, and methylation of other genes in our panel involved in DNA damage repair (FANCF and MGMT) and apoptosis (DAPK and ASC) but not with other genes in our panel, including p16 (the G1-S checkpoint gene), CRBP1, or RAR-B. In addition, MGMT methylation strongly correlated with the pro-apoptotic gene, ASC. There are distinct associations of methylated genes in non-small cell lung cancer involving DNA damage repair, apoptosis, and the G2-M mitotic checkpoint control. Further studies are warranted to determine whether these methylation patterns have implications for prognosis in addition to prediction of response to chemotherapeutic agents commonly used in the treatment of non-small cell lung cancer, such as radiotherapy and platinum- or taxane-based chemotherapy.

Genomic analysis reveals that immune function genes are strongly linked to clinical outcome in the North Central Cancer Treatment Group n9831 Adjuvant Trastuzumab Trial

PURPOSE

To develop a genomic signature that predicts benefit from trastuzumab in human epidermal growth factor receptor 2-positive breast cancer.

PATIENTS AND METHODS

DASL technology was used to quantify mRNA in samples from 1,282 patients enrolled onto the Combination Chemotherapy With or Without Trastuzumab in Treating Women With Breast Cancer (North Central Cancer Treatment Group N9831 [NCCTG-N9831]) adjuvant trastuzumab trial. Cox proportional hazard ratios (HRs), adjusted for significant clinicopathologic risk factors, were used to determine the association of each gene with relapse-free survival (RFS) for 433 patients who received chemotherapy alone (arm A) and 849 patients who received chemotherapy plus trastuzumab (arms B and C). Network and pathway analyses were used to identify key biologic processes linked to RFS. The signature was built by using a voting scheme.

RESULTS

Network and functional ontology analyses suggested that increased RFS was linked to a subset of immune function genes. A voting scheme model was used to define immune gene enrichment based on the expression of any nine or more of 14 immune function genes at or above the 0.40 quantile for the population. This model was used to identify immune gene-enriched tumors in arm A and arms B and C. Immune gene enrichment was linked to increased RFS in arms B and C (HR, 0.35; 95% CI, 0.22 to 0.55; P < .001), whereas arm B and C patients who did not exhibit immune gene enrichment did not benefit from trastuzumab (HR, 0.89; 95% CI, 0.62 to 1.28; P = .53). Enriched immune function gene expression as defined by our predictive signature was not associated with increased RFS in arm A (HR, 0.90; 95% CI, 0.60 to 1.37; P = .64).

CONCLUSION

Increased expression of a subset of immune function genes may provide a means of predicting benefit from adjuvant trastuzumab.

DNA methylation and RNA expression profiles in lung adenocarcinomas of never-smokers

Lung cancer occurs in never-smokers. Epigenetic changes in lung cancer potentially represent important diagnostic, prognostic, and therapeutic targets. We compared DNA methylation profiles of 28 adenocarcinomas of the lungs of never-smokers with paired adjacent nonmalignant lung tissue. We correlated differential methylation changes with gene expression changes from the same 28 sample pairs. Using principal component analysis, we observed a distinct separation in methylation profiles between tumor and adjacent nonmalignant lung tissue. Tumors were generally hypomethylated compared with adjacent nonmalignant tissue. Of 1,906 CpG sites differentially methylated between tumor and nonmalignant tissue, 1,198 were within classically defined CpG islands where tumors were hypermethylated compared with nonmalignant tissue. A total of 708 sites were outside CpG islands where tumors were hypomethylated compared with nonmalignant tissue. There were significant differences in expression of 351 genes (23%) of the 1,522 genes matched to the differentially methylated CpG sites. Genes that were not significantly differentially expressed and were hypermethylated within CpG sites were enriched for homeobox genes. These results suggest that the methylation profiles of lung adenocarcinomas of never-smokers and adjacent nonmalignant lung tissue are significantly different. Despite the differential methylation of homeobox genes, no significant changes in expression of these genes were detected.

A DRD1 polymorphism predisposes to lung cancer among those exposed to secondhand smoke during childhood

Lung cancer has a familial component which suggests a genetic contribution to its etiology. Given the strong evidence linking smoking with lung cancer, we studied miRNA-related loci in genes associated with smoking behavior. CHRNA, CHRNB gene families, CYP2A6, and DRD1 (dopamine receptor D1) were mined for SNPs that fell within the seed region of miRNA binding sites and then tested for associations with risk in a three-stage validation approach. A 3'UTR (untranslated region) SNP in DRD1 was associated with a lower risk of lung cancer among individuals exposed to secondhand smoke during childhood [OR, 0.69; 95% confidence interval (CI), 0.60-0.79; P < 0.0001]. This relationship was evident in both ever (OR, 0.74; 95% CI, 0.62-0.88; P = 0.001) and never smokers (OR, 0.61; 95% CI, 0.47-0.79; P < 0.0001), European American (OR, 0.65; 95% CI, 0.53-0.80; P < 0.0001), and African American (OR, 0.73; 95% CI, 0.62-0.88; P = 0.001) populations. Although much remains undefined about the long-term risks associated with exposure to secondhand smoke and heterogeneity between individuals in regard to their susceptibility to the effects of secondhand smoke, our data show an interaction between an SNP in the 3'UTR of DRD1 and exposure to secondhand smoke during childhood. Further work is needed to explore the mechanistic underpinnings of this SNP and the nature of the interaction between DRD1 and exposure to secondhand smoke during childhood.

Abstract 2156: Childhood exposure to secondhand smoke, nicotine dependence, and DRD1 are associated with lung cancer risk

Children who are exposed to secondhand smoke (SHS) are faced with an increased risk of nicotine dependence later in life. We have previously shown that children exposed to SHS have an increased risk of lung cancer as adults. Furthermore, genetic variation plays a role in the susceptibility of individuals to lung cancer. The gene DRD1 mediates dopamine signaling following nicotine exposure, thus we reasoned that DRD1 may be related to lung cancer via modulation of the nicotine reward pathway. Our previous studies identified 4 SNPs (rs4809294, rs2292975, rs4861327, rs686) that were predicted to modulate a microRNA binding site and we assessed the association of these SNPs with lung cancer risk. Our current study investigated whether these SNPs were associated with nicotine dependence and whether this had a subsequent effect on lung cancer risk. We used the National Cancer Institute-Maryland lung cancer case-control study. Cases and controls were recruited from the greater Baltimore metropolitan area. We used ordinal logistic regression to examine the association between SNPs and nicotine dependence. Logistic regression adjusted for age, gender and smoking was used to assess the association between SNPs and lung cancer risk. No significant associations were found between genotypes in the 4 SNPs and nicotine dependence levels (n=1,450). Individuals with the heterozygous allele (AG) exposed to SHS during childhood had increased odds of lung cancer compared to those with the homozygous wild type allele (AA) for rs686 (OR=1.52, CI= 1.12, 2.07). Individuals with the homozygous variant allele (GG) exposed to childhood SHS had increased odds of lung cancer compared to individuals with the homozygous wild type allele (AA) (OR=2.11, CI= 1.45, 3.10). Individuals with the heterozygous allele (GA) and the highest levels of nicotine dependence who were also exposed to childhood SHS, had higher odds of lung cancer than those with the homozygous recessive allele (AA) (OR=3.50, CI=1.37, 8.96). Consistent with earlier findings, heterozygotes (AG) had increased odds of lung cancer but this was not statistically significant (OR=1.90, CI= 0.82, 4.40). To further explore this association, we examined the effect of race on this relationship. Caucasians with the heterozygous allele (AG), the highest levels of nicotine dependence and childhood SHS, had increased odds of lung cancer compared to those with the homozygous recessive allele (AA) for rs686 (OR=1.47, CI=1.06, 2.05). Caucasians with the homozygous allele (GG) who were exposed to childhood SHS, had higher odds of lung cancer compared to individuals with the homozygous recessive allele (AA) for rs686 (OR=1.84, CI=1.15, 2.98). No association was found among African Americans. Although none of the SNPs examined were found to be associated with nicotine dependence, it appears that variations in rs686 may modify an individual's susceptibility to lung cancer due to SHS exposure.

Citation Format: Victoria A. Zigmont, Brid M. Ryan, Jin Jen, Ana I. Robles, Cain McClary, Kara Calhoun, Elise D. Bowman, Kirsi Vähäkangas, K. Leigh Greathouse, Wang Yi, Angela S. Wenzlaff, Bo Deng, Ping Yang, Ann G. Schwartz, Curtis C. Harris, Susan Olivo Marston. Childhood exposure to secondhand smoke, nicotine dependence, and DRD1 are associated with lung cancer risk. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 2156. doi:10.1158/1538-7445.AM2014-2156

Gene expression profiling in human lung development: an abundant resource for lung adenocarcinoma prognosis

A tumor can be viewed as a special “organ” that undergoes aberrant and poorly regulated organogenesis. Progress in cancer prognosis and therapy might be facilitated by re-examining distinctive processes that operate during normal development, to elucidate the intrinsic features of cancer that are significantly obscured by its heterogeneity. The global gene expression signatures of 44 human lung tissues at four development stages from Asian descent and 69 lung adenocarcinoma (ADC) tissue samples from ethnic Chinese patients were profiled using microarrays. All of the genes were classified into 27 distinct groups based on their expression patterns (named as PTN1 to PTN27) during the developmental process. In lung ADC, genes whose expression levels decreased steadily during lung development (genes in PTN1) generally had their expression reactivated, while those with uniformly increasing expression levels (genes in PTN27) had their expression suppressed. The genes in PTN1 contain many n-gene signatures that are of prognostic value for lung ADC. The prognostic relevance of a 12-gene demonstrator for patient survival was characterized in five cohorts of healthy and ADC patients [ADC_CICAMS (n = 69, p = 0.007), ADC_PNAS (n = 125, p = 0.0063), ADC_GSE13213 (n = 117, p = 0.0027), ADC_GSE8894 (n = 62, p = 0.01), and ADC_NCI (n = 282, p = 0.045)] and in four groups of stage I patients [ADC_CICAMS (n = 22, p = 0.017), ADC_PNAS (n = 76, p = 0.018), ADC_GSE13213 (n = 79, p = 0.02), and ADC_qPCR (n = 62, p = 0.006)]. In conclusion, by comparison of gene expression profiles during human lung developmental process and lung ADC progression, we revealed that the genes with a uniformly decreasing expression pattern during lung development are of enormous prognostic value for lung ADC.

Clinical biomarkers of pulmonary carcinoid tumors in never smokers via profiling miRNA and target mRNA

Background

miRNAs play key regulatory roles in cellular pathological processes. We aimed to identify clinically meaningful biomarkers in pulmonary carcinoid tumors (PCTs), a member of neuroendocrine neoplasms, via profiling miRNAs and mRNAs.

Results

From the total of 1145 miRNAs, we obtained 16 and 17 miRNAs that showed positive and negative fold changes (FCs, tumors vs. normal tissues) in the top 1% differentially expressed miRNAs, respectively. We uncovered the target genes that were predicted by at least two prediction tools and overlapped by at least one-half of the top miRNAs, which yielded 44 genes (FC<-2) and 56 genes (FC>2), respectively. Higher expressions of CREB5, PTPRB and COL4A3 predicted favorable disease free survival (Hazard ratio: 0.03, 0.19 and 0.36; P value: 0.03, 0.03 and 0.08). Additionally, 79 mutated genes have been found in nine PCTs where TP53 was the only repeated mutation.

Conclusion

We identified that the expressions of three genes have clinical implications in PCTs. The biological functions of these biomarkers warrant further studies.

Association of MAPT haplotypes with Alzheimer's disease risk and MAPT brain gene expression levels

Introduction

MAPT encodes for tau, the predominant component of neurofibrillary tangles that are neuropathological hallmarks of Alzheimer’s disease (AD). Genetic association of MAPT variants with late-onset AD (LOAD) risk has been inconsistent, although insufficient power and incomplete assessment of MAPT haplotypes may account for this.

Methods

We examined the association of MAPT haplotypes with LOAD risk in more than 20,000 subjects (n-cases = 9,814, n-controls = 11,550) from Mayo Clinic (n-cases = 2,052, n-controls = 3,406) and the Alzheimer’s Disease Genetics Consortium (ADGC, n-cases = 7,762, n-controls = 8,144). We also assessed associations with brain MAPT gene expression levels measured in the cerebellum (n = 197) and temporal cortex (n = 202) of LOAD subjects. Six single nucleotide polymorphisms (SNPs) which tag MAPT haplotypes with frequencies greater than 1% were evaluated.

Results

H2-haplotype tagging rs8070723-G allele associated with reduced risk of LOAD (odds ratio, OR = 0.90, 95% confidence interval, CI = 0.85-0.95, p = 5.2E-05) with consistent results in the Mayo (OR = 0.81, p = 7.0E-04) and ADGC (OR = 0.89, p = 1.26E-04) cohorts. rs3785883-A allele was also nominally significantly associated with LOAD risk (OR = 1.06, 95% CI = 1.01-1.13, p = 0.034). Haplotype analysis revealed significant global association with LOAD risk in the combined cohort (p = 0.033), with significant association of the H2 haplotype with reduced risk of LOAD as expected (p = 1.53E-04) and suggestive association with additional haplotypes. MAPT SNPs and haplotypes also associated with brain MAPT levels in the cerebellum and temporal cortex of AD subjects with the strongest associations observed for the H2 haplotype and reduced brain MAPT levels (β = -0.16 to -0.20, p = 1.0E-03 to 3.0E-03).

Conclusions

These results confirm the previously reported MAPT H2 associations with LOAD risk in two large series, that this haplotype has the strongest effect on brain MAPT expression amongst those tested and identify additional haplotypes with suggestive associations, which require replication in independent series. These biologically congruent results provide compelling evidence to screen the MAPT region for regulatory variants which confer LOAD risk by influencing its brain gene expression.

Conserved recurrent gene mutations correlate with pathway deregulation and clinical outcomes of lung adenocarcinoma in never-smokers

BACKGROUND

Novel and targetable mutations are needed for improved understanding and treatment of lung cancer in never-smokers.

METHODS

Twenty-seven lung adenocarcinomas from never-smokers were sequenced by both exome and mRNA-seq with respective normal tissues. Somatic mutations were detected and compared with pathway deregulation, tumor phenotypes and clinical outcomes.

RESULTS

Although somatic mutations in DNA or mRNA ranged from hundreds to thousands in each tumor, the overlap mutations between the two were only a few to a couple of hundreds. The number of somatic mutations from either DNA or mRNA was not significantly associated with clinical variables; however, the number of overlap mutations was associated with cancer subtype. These overlap mutants were preferentially expressed in mRNA with consistently higher allele frequency in mRNA than in DNA. Ten genes (EGFR, TP53, KRAS, RPS6KB2, ATXN2, DHX9, PTPN13, SP1, SPTAN1 and MYOF) had recurrent mutations and these mutations were highly correlated with pathway deregulation and patient survival.

CONCLUSIONS

The recurrent mutations present in both DNA and RNA are likely the driver for tumor biology, pathway deregulation and clinical outcomes. The information may be used for patient stratification and therapeutic target development.

Evidence that the lung Adenocarcinoma EML4-ALK fusion gene is not caused by exposure to secondhand tobacco smoke during childhood

BACKGROUND

The EML4-ALK fusion gene is more frequently found in younger, never smoking patients with lung cancer. Meanwhile, never smokers exposed to secondhand tobacco smoke (SHS) during childhood are diagnosed at a younger age compared with never smoking patients with lung cancer who are not exposed. We, therefore, hypothesized that SHS, which can induce DNA damage, is associated with the EML4-ALK fusion gene.

METHODS

We compared the frequency of the EML4-ALK fusion gene among 197 never smoker patients with lung cancer with and without a history of exposure to SHS during childhood at Mayo Clinic.

RESULTS

The EML4-ALK fusion gene was detected in 33% of cases from never smokers with a history of SHS exposure during childhood, whereas 47% of never smoking lung cancer cases without a history of childhood SHS exposure tested positive for the fusion gene.

CONCLUSIONS

The EML4-ALK fusion gene is not enriched in tumors from individuals exposed to SHS during childhood.

IMPACT

These data suggest that childhood exposure to SHS is not a significant etiologic cause of the EML4-ALK fusion gene in lung cancer.

Gene expression, single nucleotide variant and fusion transcript discovery in archival material from breast tumors

Advantages of RNA-Seq over array based platforms are quantitative gene expression and discovery of expressed single nucleotide variants (eSNVs) and fusion transcripts from a single platform, but the sensitivity for each of these characteristics is unknown. We measured gene expression in a set of manually degraded RNAs, nine pairs of matched fresh-frozen, and FFPE RNA isolated from breast tumor with the hybridization based, NanoString nCounter (226 gene panel) and with whole transcriptome RNA-Seq using RiboZeroGold ScriptSeq V2 library preparation kits. We performed correlation analyses of gene expression between samples and across platforms. We then specifically assessed whole transcriptome expression of lincRNA and discovery of eSNVs and fusion transcripts in the FFPE RNA-Seq data. For gene expression in the manually degraded samples, we observed Pearson correlations of >0.94 and >0.80 with NanoString and ScriptSeq protocols, respectively. Gene expression data for matched fresh-frozen and FFPE samples yielded mean Pearson correlations of 0.874 and 0.783 for NanoString (226 genes) and ScriptSeq whole transcriptome protocols respectively, p<2x10(-16). Specifically for lincRNAs, we observed superb Pearson correlation (0.988) between matched fresh-frozen and FFPE pairs. FFPE samples across NanoString and RNA-Seq platforms gave a mean Pearson correlation of 0.838. In FFPE libraries, we detected 53.4% of high confidence SNVs and 24% of high confidence fusion transcripts. Sensitivity of fusion transcript detection was not overcome by an increase in depth of sequencing up to 3-fold (increase from ~56 to ~159 million reads). Both NanoString and ScriptSeq RNA-Seq technologies yield reliable gene expression data for degraded and FFPE material. The high degree of correlation between NanoString and RNA-Seq platforms suggests discovery based whole transcriptome studies from FFPE material will produce reliable expression data. The RiboZeroGold ScriptSeq protocol performed particularly well for lincRNA expression from FFPE libraries, but detection of eSNV and fusion transcripts was less sensitive.

High expression of SNIP1 correlates with poor prognosis in non-small cell lung cancer and SNIP1 interferes with the recruitment of HDAC1 to RB in vitro

The Rb tumor suppressor gene performs a critical role in controlling cell proliferation and tumorigenesis; it recruits HDAC1 protein into the E2F complexes to repress transcription. In this study, we demonstrate that SNIP1, RB and HDAC1 were significantly expressed in same lung cancer tissues in a tissue microarray (TMA) containing 300 non-small cell lung cancers (NSCLC). High expression level of SNIP1 in tumor patients was significantly correlated with poor prognosis in NSCLC (log-rank P for OS = 0.01, log-rank P for DFS = 0.001). Functionally, SNIP1 competes with HDAC1 for binding to RB and reduces HDAC activity in vitro. Knockdown of SNIP1 reduced colony formation ability of lung cancer cells. These findings may indicate the involvement of SNIP1 in progression of lung cancer by regulating the RB/HDAC1 interaction.

Review of 200 consecutive patients with mutation profiling in a lung cancer individualized medicine clinic

e19013

Background: Mutation profiling to assess for potentially druggable mutations in NSCLC is being offered at an increasing number of cancer centers throughout North America and internationally. Although data continue to accumulate for the potential value of mutation testing in designing chemotherapeutic regimens, the treatment impact of obtaining information beyond assessment of EGFR and ALK status remains unclear. How best to obtain and clinically utilize these data, including information from exome and whole genome sequencing, also remains unclear. Methods: Patients were reviewed electronically in a multidisciplinary conference regarding indications for testing and results of mutation profiling from various methods, including the mass-spec based LungCarta test, targeted NexGen sequencing, exome, and whole genome sequencing. Outcomes of the multidisciplinary review were communicated back to treating physicians. Results: Mutation testing was performed on 200 patients using a variety of approaches. The majority (>150) were surgically resected stage I and II tumors. Mutations in at least 1 major cancer driver gene, including EGFR, KRAS, MET, BRAF and PIK3CA, were found in 47% of all patients tested. EGFR mutations were present in 14.8% of patients tested, KRAS 21.3%, BRAF 2.6%, PIK3CA 3.2%, and MET 4.5%. A total of 8 patients underwent either exome or whole genome sequencing. A limited number of patients (<10) had mutation results that impacted treatment decisions from this cohort. Conclusions: Mutation profiling can influence treatment decisions in NSCLC, but at a low frequency. The role of exome and or whole genome sequencing for patients with NSCLC is evolving and remains undefined.

Abstract 4581: Interaction between DRD1 and childhood exposure to environmental tobacco smoke modulates lung cancer risk in smokers and never smokers

In searching for a role for genetic variation in susceptibility to lung cancer, we postulated that miRNA modulating SNPs within smoking behavior associated genes, such as the nicotinic acetylcholine and dopamine receptor families, could modulate susceptibility to tobacco carcinogens and lung cancer risk. Bioinformatic analyses of these gene families identified 5 SNPs predicted to modulate a miRNA binding site. We assessed these SNPs in the National Cancer Institute-MD case-control study (European American [EA] & African Americans [AA] n=1,438). rs686 in DRD1 (dopamine receptor D1) was associated with a reduced risk of lung cancer after adjustment for age, gender and race (OR 0.63, 95% C.I. 0.47-0.85; P=0.002). As DRD1 mediates dopamine signaling following nicotine exposure, we reasoned that the rs686 variant associates with lung cancer risk via modulation of the nicotine reward pathway. However, adjustment for smoking both in terms of smoking status and pack-years did not significantly modulate the association (OR 0.57, 95% C.I. 0.40-0.81, P=0002). In addition, age at smoking initiation, pack-years of smoking and levels of urinary nicotine metabolites (cotinine and hydroxycotinine) did not vary across rs686 genotypes. Moreover, in an analysis of a never smoking cohort from Mayo Clinic, rs686 was also associated with a reduced risk of lung cancer (OR 0.77, 95% C.I. 0.62-0.97; P=0.027; n=629) (pooled analysis OR 0.61, 95% C.I. 0.47-0.78; P&lt;0.0001; n=2,067). Intriguingly, the association between rs686 and lung cancer risk was only observed in subjects exposed to tobacco smoke during childhood (EA smokers: [OR 0.63, 95% C.I. 0.45-0.88; P=0.007] EA never smokers [OR 0.64, 95% C.I. 0.49-0.85; P=0.002] AA smokers OR 0.31, 95% C.I. 0.11-0.89; P=0.029) (model adjusted as above and for education, adult-second-hand smoke exposure). To further explore a relationship between tobacco and DRD1, we exposed lung cancer cell lines to increasing doses of nicotine and observed a significant increase in DRD1 expression (P=0.009). As DRD1 lies within a CpG island, we are currently testing if this regulation occurs via epigenetic modulation of the DRD1 promoter. While a direct role for DRD1 in lung cancer has not previously been reported, a recent report linked DRD to the cancer stem cell hypothesis. Indeed, our work demonstrates that DRD1 is significantly upregulated in lung cells with cancer stem-like properties. Furthermore, miR-296, predicted to bind to variant rs686 allele, is increased in paired lung cancer tissues (P=0.0003), thus linking a decrease in DRD1 with reduced lung cancer risk. Collectively, our work suggests a novel lung cancer pathway whereby rs686 in DRD1 counteracts the increase in DRD1 imposed by exposure to tobacco smoke in both current and never smokers. Furthermore, it hints at a novel role for the dopamine receptor/signaling pathway in lung carcinogenesis.

Citation Format: Brid M. Ryan, Ana I. Robles, Andrew C. McClary, Elise Bowman, Kirsi Vahakangas, Susan Olivo-Marston, Ping Yang, Jin Jen, Curtis C. Harris. Interaction between DRD1 and childhood exposure to environmental tobacco smoke modulates lung cancer risk in smokers and never smokers. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 4581. doi:10.1158/1538-7445.AM2013-4581

Abstract 2005: Analysis of gene expression and copy number variation in breast tumors using both sequencing and hybridization-based platforms

Next Generation Sequencing (NGS) technologies provide rapid genomic analyses of single nucleotide variants, RNA expression and DNA copy number. Application of these technologies to material isolated from formalin fixed paraffin embedded (FFPE) tissue and even degraded frozen material could provide powerful replication samples but remains challenging. We tested the nanoString platform to validate deep sequence analysis of gene expression and DNA copy number in degraded and FFPE material. Firstly, RNA from the Universal Human Reference RNA and a breast cancer cell line (MDA-MB-436) was artificially degraded to different degrees (RIN 1.2-6.8). We used the nanoString platform to simultaneously measure RNA expression across 226 genes in each degraded sample and the corresponding undegraded RNA. Secondly we isolated RNA and DNA from matched fresh frozen and FFPE tissues from nine breast cancer patients (3 HER2+/ER+/PR+, 2 HER2+/ER-/PR-, 2 HER2+/ER+/PR-, 2 HER2-/ER+/PR+) using the nanoString platform to compare expression and copy number across 226 and 86 genes respectively. Finally, we correlated expression and copy number data generated by nanoString with Illumina transcriptome and whole genome sequencing (WGS). NanoString log2 expression fold-change between all artificially degraded samples and their undegraded counterpart showed extremely high correlation (r2&gt;0.91). NanoString DNA copy number between matched fresh-frozen and FFPE showed a high degree of correlation (r2=0.71). All gene amplifications with copy number ≥ 5 in DNA from fresh-frozen material (N=9) were successfully identified in DNA from FFPE material. We also observed good correlation of gene expression between whole transcriptome sequencing and the nanoString platform (r2 0.59 - 0.72) in FFPE and artificially degraded material and for DNA copy number between WGS and nanoString in DNA isolated from cancer cell lines (r2=0.96). The nanoString platform provides reliable data from highly degraded and FFPE material and correlates with sequence analysis of both expression and copy number from NGS platforms demonstrating potential for large-scale replication studies in FFPE material.

Citation Format: Nadine Norton, Edith A. Perez, Yan W. Asmann, Jennifer M. Carr, Brian M. Necela, Jennifer M. Kachergus, Jin Jen, Bruce W. Eckloff, E Aubrey Thompson. Analysis of gene expression and copy number variation in breast tumors using both sequencing and hybridization-based platforms. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 2005. doi:10.1158/1538-7445.AM2013-2005

Regulation of the SK3 channel by microRNA-499--potential role in atrial fibrillation

BACKGROUND

MicroRNAs are important regulators of gene expression, including those involving electrical remodeling in atrial fibrillation (AF). Recently, KCNN3, the gene that encodes the small-conductance calcium-activated potassium channel 3 (SK3), was found to be strongly associated with AF.

OBJECTIVES

To evaluate the changes in atrial myocardial microRNAs in patients with permanent AF and to determine the role of microRNA on the regulation of cardiac SK3 expression.

METHODS

Atrial tissue obtained during cardiac surgery from patients (4 sinus rhythm and 4 permanent AF) was analyzed by using microRNA arrays. Potential targets of microRNAs were predicted by using software programs. The effects of specific microRNAs on target gene expression were evaluated in HL-1 cells from a continuously proliferating mouse hyperplastic atrial cardiomyocyte cell line. Interactions between microRNAs and targets were further evaluated by using luciferase reporter assay and by using Argonaute pull-down assay.

RESULTS

Twenty-one microRNAs showed significant (>2-fold) changes in AF. MicroRNA 499 (miR-499) was upregulated by 2.33-fold (P < .01) in AF atria, whereas SK3 protein expression was downregulated by 46% (P < .05). Transfection of miR-499 mimic in HL-1 cells resulted in the downregulation of SK3 protein expression, while that of miR-499 inhibitor upregulated SK3 expression. Binding of miR-499 to the 3' untranslated region of KCNN3 was confirmed by luciferase reporter assay and by the increased presence of SK3 mRNA in Argonaute pulled-down microRNA-induced silencing complexes after transfection with miR-499.

CONCLUSION

Atrial miR-499 is significantly upregulated in AF, leading to SK3 downregulation and possibly contributing to the electrical remodeling in AF.

Quality assessment metrics for whole genome gene expression profiling of paraffin embedded samples

BACKGROUND

Formalin fixed, paraffin embedded tissues are most commonly used for routine pathology analysis and for long term tissue preservation in the clinical setting. Many institutions have large archives of Formalin fixed, paraffin embedded tissues that provide a unique opportunity for understanding genomic signatures of disease. However, genome-wide expression profiling of Formalin fixed, paraffin embedded samples have been challenging due to RNA degradation. Because of the significant heterogeneity in tissue quality, normalization and analysis of these data presents particular challenges. The distribution of intensity values from archival tissues are inherently noisy and skewed due to differential sample degradation raising two primary concerns; whether a highly skewed array will unduly influence initial normalization of the data and whether outlier arrays can be reliably identified.

FINDINGS

Two simple extensions of common regression diagnostic measures are introduced that measure the stress an array undergoes during normalization and how much a given array deviates from the remaining arrays post-normalization. These metrics are applied to a study involving 1618 formalin-fixed, paraffin-embedded HER2-positive breast cancer samples from the N9831 adjuvant trial processed with Illumina's cDNA-mediated Annealing Selection extension and Ligation assay.

CONCLUSION

Proper assessment of array quality within a research study is crucial for controlling unwanted variability in the data. The metrics proposed in this paper have direct biological interpretations and can be used to identify arrays that should either be removed from analysis all together or down-weighted to reduce their influence in downstream analyses.

Skeletal muscle and bone marrow derived stromal cells: a comparison of tenocyte differentiation capabilities

This study investigated the comparative ability of bone marrow and skeletal muscle derived stromal cells (BMSCs and SMSCs) to express a tenocyte phenotype, and whether this expression could be augmented by growth and differentiation factor-5 (GDF-5). Tissue harvest was performed on the hind limbs of seven dogs. Stromal cells were isolated via serial expansion in culture. After four passages, tenogenesis was induced using either ascorbic acid alone or in conjunction with GDF-5. CD44, tenomodulin, collagen I, and collagen III expression levels were compared for each culture condition at 7 and 14 days following induction. Immunohistochemistry (IHC) was performed to evaluate cell morphology and production of tenomodulin and collagen I. SMSCs and BMSCs were successfully isolated in culture. Following tenocytic induction, SMSCs demonstrated an increased mean relative expression of tenomodulin, collagen I, and collagen III at 14 days. BMSCs only showed increased mean relative expression of collagen I, and collagen III at 14 days. IHC revealed positive staining for tenomodulin and collagen I at 14 days for both cell types. The morphology of skeletal muscle derived stromal cells at 14 days had an organized appearance in contrast to the haphazard arrangement of the bone marrow derived cells. GDF-5 did not affect gene expression, cell staining, or cell morphology significantly. Stromal cells from either bone marrow or skeletal muscle can be induced to increase expression of matrix genes; however, based on expression of tenomodulin and cell culture morphology SMSCs may be a more ideal candidate for tenocytic differentiation.

Novel late-onset Alzheimer disease loci variants associate with brain gene expression

OBJECTIVE

Recent genome-wide association studies (GWAS) of late-onset Alzheimer disease (LOAD) identified 9 novel risk loci. Discovery of functional variants within genes at these loci is required to confirm their role in Alzheimer disease (AD). Single nucleotide polymorphisms that influence gene expression (eSNPs) constitute an important class of functional variants. We therefore investigated the influence of the novel LOAD risk loci on human brain gene expression.

METHODS

We measured gene expression levels in the cerebellum and temporal cortex of autopsied AD subjects and those with other brain pathologies (∼400 total subjects). To determine whether any of the novel LOAD risk variants are eSNPs, we tested their cis-association with expression of 6 nearby LOAD candidate genes detectable in human brain (ABCA7, BIN1, CLU, MS4A4A, MS4A6A, PICALM) and an additional 13 genes ±100 kb of these SNPs. To identify additional eSNPs that influence brain gene expression levels of the novel candidate LOAD genes, we identified SNPs ±100 kb of their location and tested for cis-associations.

RESULTS

CLU rs11136000 (p = 7.81 × 10(-4)) and MS4A4A rs2304933/rs2304935 (p = 1.48 × 10(-4)-1.86 × 10(-4)) significantly influence temporal cortex expression levels of these genes. The LOAD-protective CLU and risky MS4A4A locus alleles associate with higher brain levels of these genes. There are other cis-variants that significantly influence brain expression of CLU and ABCA7 (p = 4.01 × 10(-5)-9.09 × 10(-9)), some of which also associate with AD risk (p = 2.64 × 10(-2)-6.25 × 10(-5)).

CONCLUSIONS

CLU and MS4A4A eSNPs may at least partly explain the LOAD risk association at these loci. CLU and ABCA7 may harbor additional strong eSNPs. These results have implications in the search for functional variants at the novel LOAD risk loci.

Brain expression genome-wide association study (eGWAS) identifies human disease-associated variants

Genetic variants that modify brain gene expression may also influence risk for human diseases. We measured expression levels of 24,526 transcripts in brain samples from the cerebellum and temporal cortex of autopsied subjects with Alzheimer's disease (AD, cerebellar n = 197, temporal cortex n = 202) and with other brain pathologies (non–AD, cerebellar n = 177, temporal cortex n = 197). We conducted an expression genome-wide association study (eGWAS) using 213,528 cisSNPs within ±100 kb of the tested transcripts. We identified 2,980 cerebellar cisSNP/transcript level associations (2,596 unique cisSNPs) significant in both ADs and non–ADs (q<0.05, p = 7.70×10⁻⁵–1.67×10⁻⁸²). Of these, 2,089 were also significant in the temporal cortex (p = 1.85×10⁻⁵–1.70×10⁻¹⁴¹). The top cerebellar cisSNPs had 2.4-fold enrichment for human disease-associated variants (p<10⁻⁶). We identified novel cisSNP/transcript associations for human disease-associated variants, including progressive supranuclear palsy SLCO1A2/rs11568563, Parkinson's disease (PD) MMRN1/rs6532197, Paget's disease OPTN/rs1561570; and we confirmed others, including PD MAPT/rs242557, systemic lupus erythematosus and ulcerative colitis IRF5/rs4728142, and type 1 diabetes mellitus RPS26/rs1701704. In our eGWAS, there was 2.9–3.3 fold enrichment (p<10⁻⁶) of significant cisSNPs with suggestive AD–risk association (p<10⁻³) in the Alzheimer's Disease Genetics Consortium GWAS. These results demonstrate the significant contributions of genetic factors to human brain gene expression, which are reliably detected across different brain regions and pathologies. The significant enrichment of brain cisSNPs among disease-associated variants advocates gene expression changes as a mechanism for many central nervous system (CNS) and non–CNS diseases. Combined assessment of expression and disease GWAS may provide complementary information in discovery of human disease variants with functional implications. Our findings have implications for the design and interpretation of eGWAS in general and the use of brain expression quantitative trait loci in the study of human disease genetics.

Genetic variants that regulate gene expression levels can also influence human disease risk. Discovery of genomic loci that alter brain gene expression levels (brain expression quantitative trait loci = eQTLs) can be instrumental in the identification of genetic risk underlying both central nervous system (CNS) and non–CNS diseases. To systematically assess the role of brain eQTLs in human disease and to evaluate the influence of brain region and pathology in eQTL mapping, we performed an expression genome-wide association study (eGWAS) in 773 brain samples from the cerebellum and temporal cortex of ∼200 autopsied subjects with Alzheimer's disease (AD) and ∼200 with other brain pathologies (non–AD). We identified ∼3,000 significant associations between cisSNPs near ∼700 genes and their cerebellar transcript levels, which replicate in ADs and non–ADs. More than 2,000 of these associations were reproducible in the temporal cortex. The top cisSNPs are enriched for both CNS and non–CNS disease-associated variants. We identified novel and confirmed previous cisSNP/transcript associations for many disease loci, suggesting gene expression regulation as their mechanism of action. These findings demonstrate the reproducibility of the eQTL approach across different brain regions and pathologies, and advocate the combined use of gene expression and disease GWAS for identification and functional characterization of human disease-associated variants.

Increased miR-708 expression in NSCLC and its association with poor survival in lung adenocarcinoma from never smokers

PURPOSE

miRNA plays an important role in human disease and cancer. We seek to investigate the expression status, clinical relevance, and functional role of miRNA in non-small cell lung cancer.

EXPERIMENTAL DESIGN

We conducted miRNA expression profiling in matched lung adenocarcinoma and uninvolved lung using 56 pairs of fresh-frozen (FF) and 47 pairs of formalin-fixed, paraffin-embedded (FFPE) samples from never smokers. The most differentially expressed miRNA genes were evaluated by Cox analysis and log-rank test. Among the best candidate, miR-708 was further examined for differential expression in two independent cohorts. Functional significance of miR-708 expression in lung cancer was examined by identifying its candidate mRNA target and through manipulating its expression levels in cultured cells.

RESULTS

Among the 20 miRNAs most differentially expressed between tested tumor and normal samples, high expression level of miR-708 in the tumors was most strongly associated with an increased risk of death after adjustments for all clinically significant factors including age, sex, and tumor stage (FF cohort: HR, 1.90; 95% CI, 1.08-3.35; P = 0.025 and FFPE cohort: HR, 1.93; 95% CI, 1.02-3.63; P = 0.042). The transcript for TMEM88 gene has a miR-708 binding site in its 3' UTR and was significantly reduced in tumors high of miR-708. Forced miR-708 expression reduced TMEM88 transcript levels and increased the rate of cell proliferation, invasion, and migration in culture.

CONCLUSIONS

miRNA-708 acts as an oncogene contributing to tumor growth and disease progression by directly downregulating TMEM88, a negative regulator of the Wnt signaling pathway in lung cancer.

EGFR somatic mutations in lung tumors: radon exposure and passive smoking in former- and never-smoking U.S. women

BACKGROUND

Patients with lung cancer with mutations in EGF receptor (EGFR) tyrosine kinase have improved prognosis when treated with EGFR inhibitors. We hypothesized that EGFR mutations may be related to residential radon or passive tobacco smoke.

METHODS

This hypothesis was investigated by analyzing EGFR mutations in 70 lung tumors from a population of never and long-term former female smokers from Missouri with detailed exposure assessments. The relationship with passive smoking was also examined in never-smoking female lung cancer cases from the Mayo clinic.

RESULTS

Overall, the frequency of EGFR mutation was 41% [95% confidence interval (CI), 32%-49%]. Neither radon nor passive-smoking exposure was consistently associated with EGFR mutations in lung tumors.

CONCLUSIONS

The results suggest that EGFR mutations are common in female, never-smoking lung cancer cases from the United States, and EGFR mutations are unlikely due to exposure to radon or passive smoking.

Glutathione S-transferase omega genes in Alzheimer and Parkinson disease risk, age-at-diagnosis and brain gene expression: an association study with mechanistic implications

BACKGROUND

Glutathione S-transferase omega-1 and 2 genes (GSTO1, GSTO2), residing within an Alzheimer and Parkinson disease (AD and PD) linkage region, have diverse functions including mitigation of oxidative stress and may underlie the pathophysiology of both diseases. GSTO polymorphisms were previously reported to associate with risk and age-at-onset of these diseases, although inconsistent follow-up study designs make interpretation of results difficult. We assessed two previously reported SNPs, GSTO1 rs4925 and GSTO2 rs156697, in AD (3,493 ADs vs. 4,617 controls) and PD (678 PDs vs. 712 controls) for association with disease risk (case-controls), age-at-diagnosis (cases) and brain gene expression levels (autopsied subjects).

RESULTS

We found that rs156697 minor allele associates with significantly increased risk (odds ratio = 1.14, p = 0.038) in the older ADs with age-at-diagnosis > 80 years. The minor allele of GSTO1 rs4925 associates with decreased risk in familial PD (odds ratio = 0.78, p = 0.034). There was no other association with disease risk or age-at-diagnosis. The minor alleles of both GSTO SNPs associate with lower brain levels of GSTO2 (p = 4.7 × 10-11-1.9 × 10-27), but not GSTO1. Pathway analysis of significant genes in our brain expression GWAS, identified significant enrichment for glutathione metabolism genes (p = 0.003).

CONCLUSION

These results suggest that GSTO locus variants may lower brain GSTO2 levels and consequently confer AD risk in older age. Other glutathione metabolism genes should be assessed for their effects on AD and other chronic, neurologic diseases.

Loss of cytoplasmic CDK1 predicts poor survival in human lung cancer and confers chemotherapeutic resistance

The dismal lethality of lung cancer is due to late stage at diagnosis and inherent therapeutic resistance. The incorporation of targeted therapies has modestly improved clinical outcomes, but the identification of new targets could further improve clinical outcomes by guiding stratification of poor-risk early stage patients and individualizing therapeutic choices. We hypothesized that a sequential, combined microarray approach would be valuable to identify and validate new targets in lung cancer. We profiled gene expression signatures during lung epithelial cell immortalization and transformation, and showed that genes involved in mitosis were progressively enhanced in carcinogenesis. 28 genes were validated by immunoblotting and 4 genes were further evaluated in non-small cell lung cancer tissue microarrays. Although CDK1 was highly expressed in tumor tissues, its loss from the cytoplasm unexpectedly predicted poor survival and conferred resistance to chemotherapy in multiple cell lines, especially microtubule-directed agents. An analysis of expression of CDK1 and CDK1-associated genes in the NCI60 cell line database confirmed the broad association of these genes with chemotherapeutic responsiveness. These results have implications for personalizing lung cancer therapy and highlight the potential of combined approaches for biomarker discovery.

Frequent homozygous deletion of the LKB1/STK11 gene in non-small cell lung cancer

LKB1/STK11 is a tumor suppressor and a negative regulator of mammalian target of rapamycin signaling. It is inactivated in 30% of lung cancer cell lines but only 5-15% of primary lung adenocarcinomas. There is evidence that homozygous deletion (HD) of chromosome 19p at the LKB locus contributes to the inactivation of the gene in primary human lung cancers. Here, we used several complementary genetic approaches to assess the LKB1 locus in primary non-small cell lung cancers (NSCLCs). We first analyzed 124 NSCLC cases for allelic imbalance using eight microsatellite markers on chromosome 19p, which revealed an overall rate of 65% (80 of 124) loss of heterozygosity (LOH). We next used chromogenic in situ hybridization (CISH) to directly examine the chromosomal status of the LKB1 locus. In all, 65 of 124 LOH tested samples were available for CISH and 58 of those (89%) showed either loss of one copy of chromosome 19p (LOH, 40 of 65 cases, 62%) or both copies (HD 18 of 65 cases, 28%). The occurrence of HD was significantly more frequent in Caucasian (35%) than in African-American patients (6%) (P=0.04). A total of 62 of 124 samples with LOH at one or both markers immediately flanking the LKB1 gene were further analyzed by directly sequencing the complete coding region, which identified 7 of 62 (11%) tumors with somatic mutations in the gene. Jointly, our data identified total inactivation of the LKB1 gene by either HD or LOH with somatic mutation in 39% of tested samples, whereas loss of chromosome 19p region by HD or LOH at the LKB1 region occured in 90% of NSCLC.

Abstract 5060: Gene network predicts overall survival in patients with primary lung adenocarcinoma

Lung adenocarcinoma is the most common type of primary lung cancer. Previous studies have shown that gene expression signatures may predict clinical outcome. In this study, we used a systems biology approach, weighted gene co-expression network analysis (WGCNA), to identify survival related networks and gene signatures. We first analyzed whole genome gene expression profiles in 82 lung adenocarcinoma tumor tissues from never smokers (smoked &lt;100 cigarettes lifetime) and found a significant correlation of an expression module with tumor grade of differentiation (r=0.38, p=4×10−4) and overall survival (OS, p =0.027, HR, 1.71; 95% CI, 1.06-2.74) after adjusting for potential confounders (age and tumor grade). When this module signature was incorporated and compared to the conventional model consisting of clinical variables only, the survival prediction accuracy was increased from 0.70 to 0.75. Gene Ontology (GO) enrichment analysis indicated this module was significantly enriched in biological process GO term “cell cycle” (p = 2.1×10−53). To explore the effect of this network module at germline level, we performed the WGCNA on normal lung tissues derived from 78 of the 82 patients and observed significant association of an expression module with OS (adjusted p=3×10−4, HR=1.98; 95% CI, 1.36-2.87). Interestingly, this module was also significantly enriched in cell cycle-related genes (p=4.1×10−41). To further validate this result, we analyzed an independent microarray gene expression dataset including 442 lung adenocarcinoma patients from NCI Director's Challenge Consortium. We found the same association: the cell cycle-enriched module was correlated with tumor grade of cell differentiation (r=0.55, p-value=5×10−37) and OS (adjusted p=7.4×10−3; HR, 1.20; 95% CI, 1.05-1.37). We further identified 30 genes in the OS-related module that overlapped between the discovery and validation data sets. Furthermore, all these 30 genes were found to be up-regulated in adenocarcinoma compared to adjacent normal lung tissues in the discovery data set. Of the 30 genes, three key genes, UBE2C, TPX2 and MELK, were identified with highest connectivity in cell cycle-enriched module; their expression levels were significantly associated with OS in both discovery and validation sets. UBE2C and TPX2 have been reported as gene signatures for human lung cancer prognosis in vitro lung carcinogenesis system; MELK was identified as a key regulator of the proliferation of malignant brain tumors and malignant grade in human astrocytomas, and was also associated with breast cancer prognosis. Our results suggest that genes involved in cell-cycle are likely related with tumor aggressiveness and therefore can predict survival in lung adenocarcinoma. Further studies are needed to validate the key genes that could improve survival prediction and provide potential new targets to intervene progression of lung adenocarcinoma.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 102nd Annual Meeting of the American Association for Cancer Research; 2011 Apr 2-6; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2011;71(8 Suppl):Abstract nr 5060. doi:10.1158/1538-7445.AM2011-5060

Quantitative miRNA expression analysis using fluidigm microfluidics dynamic arrays

Background

MicroRNAs (miRNAs) represent a growing class of small non-coding RNAs that are important regulators of gene expression in both plants and animals. Studies have shown that miRNAs play a critical role in human cancer and they can influence the level of cell proliferation and apoptosis by modulating gene expression. Currently, methods for the detection and measurement of miRNA expression include small and moderate-throughput technologies, such as standard quantitative PCR and microarray based analysis. However, these methods have several limitations when used in large clinical studies where a high-throughput and highly quantitative technology needed for the efficient characterization of a large number of miRNA transcripts in clinical samples. Furthermore, archival formalin fixed, paraffin embedded (FFPE) samples are increasingly becoming the primary resource for gene expression studies because fresh frozen (FF) samples are often difficult to obtain and requires special storage conditions. In this study, we evaluated the miRNA expression levels in FFPE and FF samples as well as several lung cancer cell lines employing a high throughput qPCR-based microfluidic technology. The results were compared to standard qPCR and hybridization-based microarray platforms using the same samples.

Results

We demonstrated highly correlated Ct values between multiplex and singleplex RT reactions in standard qPCR assays for miRNA expression using total RNA from A549 (R = 0.98; p < 0.0001) and H1299 (R = 0.95; p < 0.0001) lung cancer cell lines. The Ct values generated by the microfluidic technology (Fluidigm 48.48 dynamic array systems) resulted in a left-shift toward lower Ct values compared to those observed by ABI 7900 HT (mean difference, 3.79), suggesting that the microfluidic technology exhibited a greater sensitivity. In addition, we show that as little as 10 ng total RNA can be used to reliably detect all 48 or 96 tested miRNAs using a 96-multiplexing RT reaction in both FFPE and FF samples. Finally, we compared miRNA expression measurements in both FFPE and FF samples by qPCR using the 96.96 dynamic array and Affymetrix microarrays. Fold change comparisons for comparable genes between the two platforms indicated that the overall correlation was R = 0.60. The maximum fold change detected by the Affymetrix microarray was 3.5 compared to 13 by the 96.96 dynamic array.

Conclusion

The qPCR-array based microfluidic dynamic array platform can be used in conjunction with multiplexed RT reactions for miRNA gene expression profiling. We showed that this approach is highly reproducible and the results correlate closely with the existing singleplex qPCR platform at a throughput that is 5 to 20 times higher and a sample and reagent usage that was approximately 50-100 times lower than conventional assays. We established optimal conditions for using the Fluidigm microfluidic technology for rapid, cost effective, and customizable arrays for miRNA expression profiling and validation.

Expression profiling of formalin-fixed paraffin-embedded primary breast tumors using cancer-specific and whole genome gene panels on the DASL® platform

Background

The cDNA-mediated Annealing, extension, Selection and Ligation (DASL) assay has become a suitable gene expression profiling system for degraded RNA from paraffin-embedded tissue. We examined assay characteristics and the performance of the DASL 502-gene Cancer Panel^v1 (1.5K) and 24,526-gene panel (24K) platforms at differentiating nine human epidermal growth factor receptor 2- positive (HER2+) and 11 HER2-negative (HER2-) paraffin-embedded breast tumors.

Methods

Bland-Altman plots and Spearman correlations evaluated intra/inter-panel agreement of normalized expression values. Unequal-variance t-statistics tested for differences in expression levels between HER2 + and HER2 - tumors. Regulatory network analysis was performed using Metacore (GeneGo Inc., St. Joseph, MI).

Results

Technical replicate correlations ranged between 0.815-0.956 and 0.986-0.997 for the 1.5K and 24K panels, respectively. Inter-panel correlations of expression values for the common 498 genes across the two panels ranged between 0.485-0.573. Inter-panel correlations of expression values of 17 probes with base-pair sequence matches between the 1.5K and 24K panels ranged between 0.652-0.899. In both panels, erythroblastic leukemia viral oncogene homolog 2 (ERBB2) was the most differentially expressed gene between the HER2 + and HER2 - tumors and seven additional genes had p-values < 0.05 and log2 -fold changes > |0.5| in expression between HER2 + and HER2 - tumors: topoisomerase II alpha (TOP2A), cyclin a2 (CCNA2), v-fos fbj murine osteosarcoma viral oncogene homolog (FOS), wingless-type mmtv integration site family, member 5a (WNT5A), growth factor receptor-bound protein 7 (GRB7), cell division cycle 2 (CDC2), and baculoviral iap repeat-containing protein 5 (BIRC5). The top 52 discriminating probes from the 24K panel are enriched with genes belonging to the regulatory networks centered around v-myc avian myelocytomatosis viral oncogene homolog (MYC), tumor protein p53 (TP53), and estrogen receptor α (ESR1). Network analysis with a two-step extension also showed that the eight discriminating genes common to the 1.5K and 24K panels are functionally linked together through MYC, TP53, and ESR1.

Conclusions

The relative RNA abundance obtained from two highly differing density gene panels are correlated with eight common genes differentiating HER2 + and HER2 - breast tumors. Network analyses demonstrated biological consistency between the 1.5K and 24K gene panels.

Genetic variants and risk of lung cancer in never smokers: a genome-wide association study

BACKGROUND

Lung cancer in individuals who have never smoked tobacco products is an increasing medical and public-health issue. We aimed to unravel the genetic basis of lung cancer in never smokers.

METHODS

We did a four-stage investigation. First, a genome-wide association study of single nucleotide polymorphisms (SNPs) was done with 754 never smokers (377 matched case-control pairs at Mayo Clinic, Rochester, MN, USA). Second, the top candidate SNPs from the first study were validated in two independent studies among 735 (MD Anderson Cancer Center, Houston, TX, USA) and 253 (Harvard University, Boston, MA, USA) never smokers. Third, further replication of the top SNP was done in 530 never smokers (UCLA, Los Angeles, CA, USA). Fourth, expression quantitative trait loci (eQTL) and gene-expression differences were analysed to further elucidate the causal relation between the validated SNPs and the risk of lung cancer in never smokers.

FINDINGS

44 top candidate SNPs were identified that might alter the risk of lung cancer in never smokers. rs2352028 at chromosome 13q31.3 was subsequently replicated with an additive genetic model in the four independent studies, with a combined odds ratio of 1.46 (95% CI 1.26-1.70, p=5.94x10(-6)). A cis eQTL analysis showed there was a strong correlation between genotypes of the replicated SNPs and the transcription level of the gene GPC5 in normal lung tissues (p=1.96x10(-4)), with the high-risk allele linked with lower expression. Additionally, the transcription level of GPC5 in normal lung tissue was twice that detected in matched lung adenocarcinoma tissue (p=6.75x10(-11)).

INTERPRETATION

Genetic variants at 13q31.3 alter the expression of GPC5, and are associated with susceptibility to lung cancer in never smokers. Downregulation of GPC5 might contribute to the development of lung cancer in never smokers.

FUNDING

US National Institutes of Health; Mayo Foundation.

Childhood exposure to secondhand smoke and functional mannose binding lectin polymorphisms are associated with increased lung cancer risk

BACKGROUND

Exposure to secondhand smoke during adulthood has detrimental health effects, including increased lung cancer risk. Compared with adults, children may be more susceptible to secondhand smoke. This susceptibility may be exacerbated by alterations in inherited genetic variants of innate immunity genes. We hypothesized a positive association between childhood secondhand smoke exposure and lung cancer risk that would be modified by genetic polymorphisms in the mannose binding lectin-2 (MBL2) gene resulting in well-known functional changes in innate immunity.

METHODS

Childhood secondhand smoke exposure and lung cancer risk was assessed among men and women in the ongoing National Cancer Institute-Maryland Lung Cancer (NCI-MD) study, which included 624 cases and 348 controls. Secondhand smoke history was collected via in-person interviews. DNA was used for genotyping the MBL2 gene. To replicate, we used an independent case-control study from Mayo Clinic consisting of 461 never smokers, made up of 172 cases and 289 controls. All statistical tests were two-sided.

RESULTS

In the NCI-MD study, secondhand smoke exposure during childhood was associated with increased lung cancer risk among never smokers [odds ratio (OR), 2.25; 95% confidence interval (95% CI), 1.04-4.90]. This was confirmed in the Mayo study (OR, 1.47; 95% CI, 1.00-2.15). A functional MBL2 haplotype associated with high circulating levels of MBL and increased MBL2 activity was associated with increased lung cancer risk among those exposed to childhood secondhand smoke in both the NCI-MD and Mayo studies (OR, 2.52; 95% CI, 1.13-5.60, and OR, 2.78; 95% CI, 1.18-3.85, respectively).

CONCLUSIONS

Secondhand smoke exposure during childhood is associated with increased lung cancer risk among never smokers, particularly among those possessing a haplotype corresponding to a known overactive complement pathway of the innate immune system.

Overexpression of phospho-eIF4E is associated with survival through AKT pathway in non-small cell lung cancer

PURPOSE

The eukaryotic translation initiation factor complex 4E (eIF4E) is downstream in the mammalian target of rapamycin (mTOR) pathway. This study explored expression of eIF4E and its relationship with the PTEN/AKT and RAS/MEK/ERK pathways in non-small cell lung carcinoma (NSCLC).

EXPERIMENTAL DESIGN

The status of phosphorylated eIF4E (p-eIF4E), phosphorylated AKT (p-AKT), PTEN, phosphorylated tuberin (p-TSC2), phosphorylated mTOR (p-mTOR), phosphorylated S6 (p-S6), and phosphorylated Erk1/2 (p-Erk1/2) was studied using immunohistochemical analysis applied to a tissue microarray containing 300 NSCLCs. Staining results for each antibody were compared with clinical and pathologic features, and the relationship between staining results was explored.

RESULTS

Overexpression of p-eIF4E, p-AKT, p-TSC2, p-mTOR, p-S6, and p-Erk1/2 in NSCLC was found in 39.9%, 78.8%, 5.1%, 46.7%, 27.1%, and 16.6% of tumors, respectively. The phenotype of p-eIF4E correlated positively with that of p-AKT, p-TSC2, and p-S6 (P < 0.001). Overall survival in NSCLC patients was significantly shorter in cases with overexpression of p-eIF4E and p-AKT alone and in combination (log-rank P < 0.001, each). Cases with underexpression of PTEN were limited (6.4%), and this phenotype did not correlate with any clinical variable. In cluster analysis, the p-AKT/p-mTOR/p-eIF4E/p-S6-positive group had significantly shorter survival compared with the survival of all cases (P < 0.001). Multivariate analysis showed that p-eIF4E overexpression is an independent prognostic factor for NSCLC (P = 0.004).

CONCLUSIONS

This study shows that p-eIF4E expression in addition to p-AKT predicts poor prognosis in NSCLC. Moreover, the correlation between expression of p-eIF4E with p-AKT, as well as p-TSC2 and p-S6, indicates that eIF4E activation through the AKT pathway plays an important role in the progression of NSCLC.

SMAD6 contributes to patient survival in non-small cell lung cancer and its knockdown reestablishes TGF-beta homeostasis in lung cancer cells

The malignant transformation in several types of cancer, including lung cancer, results in a loss of growth inhibition by transforming growth factor-beta (TGF-beta). Here, we show that SMAD6 expression is associated with a reduced survival in lung cancer patients. Short hairpin RNA (shRNA)-mediated knockdown of SMAD6 in lung cancer cell lines resulted in reduced cell viability and increased apoptosis as well as inhibition of cell cycle progression. However, these results were not seen in Beas2B, a normal bronchial epithelial cell line. To better understand the mechanism underlying the association of SMAD6 with poor patient survival, we used a lentivirus construct carrying shRNA for SMAD6 to knock down expression of the targeted gene. Through gene expression analysis, we observed that knockdown of SMAD6 led to the activation of TGF-beta signaling through up-regulation of plasminogen activator inhibitor-1 and phosphorylation of SMAD2/3. Furthermore, SMAD6 knockdown activated the c-Jun NH2-terminal kinase pathway and reduced phosphorylation of Rb-1, resulting in increased G0-G1 cell arrest and apoptosis in the lung cancer cell line H1299. These results jointly suggest that SMAD6 plays a critical role in supporting lung cancer cell growth and survival. Targeted inactivation of SMAD6 may provide a novel therapeutic strategy for lung cancers expressing this gene.