Citations

Abstract Background. Several recent studies have reported an inverse relationship between sun exposure and non-Hodgkin (NHL) risk. One hypothesis is that vitamin D mediates this effect. We evaluated the association of sun exposure at various ages and common variants in vitamin D related genes with risk of NHL and its major subtypes. Methods. Newly diagnosed NHL cases (N=1023) and frequency matched controls (N=1254) were enrolled in a clinic-based case-control study conducted at the Mayo Clinic. Thirty-eight single nucleotide polymorphisms (SNPs) from 5 candidate genes relevant to vitamin D metabolism and sun sensitivity (IRF4, RXR, VDR, CYP24A1, CYP27B1) were genotyped. Participants also completed a detailed questionnaire, including self-reported sun exposure at ages 13-21, 22-40, and 41+ years. Odds ratios (OR), 95% confidence intervals (CI), and tests for trend and effect modification were estimated using unconditional logistic regression adjusted for age and gender. ORs are reported for the highest vs. lowest exposure category or as the ordinal OR for SNPs. Results. The mean age of diagnosis/enrollment was 61 years for cases and 60 years for controls. The most common subtypes were chronic lymphocytic leukemia/small lymphocytic lymphoma (CLL/SLL; 34%), follicular (FL; 24%), and diffuse large B-cell (DLBCL; 18%) lymphomas. There was a significant decrease in NHL risk with increased sun exposure at ages 13-21 years (OR=0.66, 95% CI 0.44-1.00; ptrend=0.04), which attenuated with increasing age at exposure. Exploratory analysis by NHL subtype indicated a trend toward a protective effect of sun exposure at ages 13-21 years for each of the major subtypes, though these estimates did not reach statistical significance. Of the 38 SNPs evaluated, there was a significant association for 3 SNPs in VDR and 1 SNP in CYP24A1: rs3819545 (OR=1.24, 95% CI 1.10-1.40, p=0.0004), rs2239186 (OR=1.22, 95% CI 1.05-1.41, p=0.009), and rs886441 (OR=0.82, 95% CI 0.70-0.96, p=0.02) for VDR and rs2762939 (OR=0.85, 95% CI 0.75-0.98, p=0.02) for CYP24A1. Both the VDR and CYP24A1 associations were observed for CLL/SLL, FL, but were attenuated for DLBCL. There was no evidence of any significant effect modification of sun exposure at ages 13-21 years by genotype in these four SNPs for NHL overall, or for NHL subtypes. Conclusions. We observed an inverse association between the frequency of sun exposure at ages 13-21 years and NHL risk, which attenuated for sun exposure at older ages, suggesting that early life exposure may be more relevant. While genetic variation in VDR and CYP24A1 were associated with NHL risk, there was no evidence for an interaction of sun exposure with these genes, suggesting that a sun association may not act through the vitamin D metabolism pathway. This would be consistent with a lack of association of circulating vitamin D with NHL risk recently published from the Cohort Consortium pooling project. 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 1894. doi:10.1158/1538-7445.AM2011-1894

The 15q13.3 microdeletion syndrome (OMIM #612001) is characterized by a wide range of phenotypic features, including intellectual disability, seizures, autism, and psychiatric conditions. This deletion is inherited in approximately 75% of cases and has been found in mildly affected and normal parents, consistent with variable expressivity and incomplete penetrance. The common deletion is approximately 2 Mb and contains several genes; however, the gene(s) responsible for the resulting clinical features have not been clearly defined. Recently, four probands were reported with small deletions including only the CHRNA7 gene. These patients showed a wide range of phenotypic features similar to those associated with the larger 15q13.3 microdeletion. To further correlate genotype and phenotype, we queried our database of >15,000 patients tested in the Mayo Clinic Cytogenetics Laboratory from 2008 to 2011 and identified 19 individuals (10 probands and 9 family members) with isolated heterozygous CHRNA7 gene deletions. All but two infants displayed multiple features consistent with 15q13.3 microdeletion syndrome. We also identified the first de novo deletion confined to CHRNA7 as well as the second known case with homozygous deletion of CHRNA7 only. These results provide further evidence implicating CHRNA7 as the gene responsible for the clinical findings associated with 15q13.3 microdeletion.

Abstract Catabolism of fluoropyrimidine drugs such as 5-fluorouracil (5-FU) and its prodrug capecitabine depends on an exclusive enzymatic step, dihydropyrimidine dehydrogenase (DPD), which results in drug inactivation. DPD exhibits wide interindividual variation, with partial to complete DPD deficiency occurring in approximately 3 to 5% of the general population. Patients with less pronounced, but decreased DPD activity may also be unable to efficiently inactivate 5-FU, such that this may contribute to 5-FU toxicity. Though numerous tests have been developed to phenotypically assess the DPD enzyme, they remain labor intensive and time consuming with difficulty in defining threshold values for separating DPD deficient patients from those with normal DPD enzyme activity. To date, attempts to develop genotypic assays have been complicated by the size and complex structure of the gene encoding the DPD enzyme (DPYD). Numerous sequence variants have been reported with varying frequencies across different populations. Unfortunately many DPYD variants have not been examined to determine whether they are associated with loss of function in the DPD enzyme. Current genotypic methods are limited due to the inability to rapidly assess multiple variants in a large number of samples at one time. In the current study, we have developed a rapid, high-throughput genotypic method using 6 multiplexed PCR assays covering 115 DPYD variants analyzed with the Sequenom matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS). Selected variants included previously identified and functionally characterized SNPs as well as newly identified mutations within the coding region from a previous mutation discovery analysis in Coriell cell lines and patient samples. Due to the potential importance of variants in the untranslated regions (UTRs), we have included known variants from the 5′- and 3′- UTRs (9 and 16 respectively) within our plexes. A total of 95 DNA samples derived from 4 colorectal cancer cell lines, 76 Coriell cell lines (40 Caucasian and 36 African American), and 7 Caucasian patients with clinically observed 5-FU toxicity and/or decreased DPD activity were utilized for initial assay testing. To assess “calling” efficiency and presence of false positives, 8 Coriell DNA samples plated in duplicate and 1 negative control were included. Our results showed that identified variants were 100% in agreement across duplicated samples, and only four variants within the 6 different plexes were removed due to PCR/extension failure. In summary, this method provides a comprehensive and rapid analysis of 111 DPYD variants within 6 different plexes with both sensitive and accurate detection suitable for high-throughput patient screening. Functional studies in a large patient population are underway (CA116964). 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 5451. doi:10.1158/1538-7445.AM2011-5451

Abstract 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 <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

Abstract Pancreatic cancer is highly lethal with an estimated 5-year relative survival rate of less than 5%. A genome wide association study (GWAS) of pancreatic cancer (PanScan) conducted within the NCI Cohort Consortium and the Pancreatic Cancer Case-Control Consortium (PANC4) recently identified four common susceptibility loci for pancreatic cancer on chromosomes 1q32.1, 5p15.33, 9q34.2 and 13q22.1. To investigate how these germline variants influence pancreatic cancer risk, we are employing both genome wide and targeted approaches. For the former, we are characterizing the epigenome and transcriptome of pancreatic cell lines and tissue samples for genomic annotation and prioritization of methods to use in targeted approaches with the aim of directly linking pancreatic cancer risk variants to molecular phenotypes and mechanisms that underlie the association signals. As an example of our approach we will present results in a pancreatic cancer susceptibility locus on chromosome 1q32.2 marked by SNP (single nucleotide polymorphism) rs3790844. This SNP is located in the first intron of the NR5A2 gene that encodes a nuclear receptor of the Ftz-F1 family known to be involved in pancreatic development, cholesterol and bile acid homeostasis, steroidogenesis and cell proliferation. In ChIP-sequencing and MeDIP-chip experiments, we observed repressive histone modification marks (H3K27me3) and regions of DNA methylation in the NR5A2 gene locus in pancreatic cell lines whereas only a few weak positive histone modification marks were observed (H3K4me1 and H3K4me3). This pattern of epigenetic marks correlates well with low expression of NR5A2 in the same cell lines. Expression of the NR5A2 gene at the RNA level is high in normal pancreatic tissue samples as compared to stage I-III pancreatic adenocarcinoma, most islet cell and neuroendocrine tumors. Using immunohistochemistry, nuclear NR5A2 was detected in normal acinar cells and – at lower levels – in normal ductal cells. The protein was also expressed in low and high grade PanINs and in the majority of PDAC analyzed. Our analysis of the tissue specific expression of NR5A2 shows that it is mainly expressed in normal pancreas, colon, small intestine, liver and breast but detected at low levels or not at all in most other tissue types and cancer cell lines. Furthermore we have identified a novel NR5A2 transcript by rapid amplification of cDNA ends (RACE) and mapped a promoter in front of this transcript in the genomic region corresponding to first intron of the two annotated Refseq transcripts. Establishing a comprehensive catalog of epigenetic patterns coupled with targeted analysis to correlate risk variants to molecular phenotypes could provide functional plausibility for known pancreatic cancer association signals discovered by GWAS findings. 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 4965. doi:10.1158/1538-7445.AM2011-4965

Abstract Background: Estrogen receptor (ER) expression in breast caner is an important biomarker for targeted therapy and outcome prediction. Nevertheless, the molecular mechanism that controls the phenotypic difference is little known. We used deep sequence technology to profile the transcriptome, gene copy number, and CpG island methylation simultaneously in eight commonly used breast cell lines to develop a model for how these genomic features are integrated in ER+ and ER- breast cancer. Materials and methods: Total mRNA sequence (mRNA-seq), gene copy number (DNA-seq), and genomic CpG island methylation (Methyl-seq) were carried out using the Illumina Genome Analyzer. Sequences were mapped to the human genome (hg18) to obtain digitized gene expression data, DNA copy number in reference to the non-tumor cell line (MCF10A), and methylation status of 21,570 CpG islands; differentially expressed genes between ER+ and ER- cell lines were selected and then correlated with methylation status of these genes’ CpG islands within 5kb of transcript start or copy number changes. The genes that appeared controlled by methylation or copy number changes in cell lines were further evaluated for their expression in the dataset of 129 primary breast tumors. Results: ER+ breast cancer cells had very different gene expression pattern from ER- cancer cells. The two cell line types formed two distinct clusters in unsupervised cluster and 1,873 genes were differentially expressed by moderated t statistics. We identified 149 differentially expressed genes that exhibited differential methylation of one or more CpG islands within 5kb of the 5’ end of the gene and for which mRNA abundance was inversely correlated with CpG island methylation status. Eighty nine of these 149 genes were also differentially expressed in the primary tumor samples and 84 of them were consistent with cell line gene expression and methylation data. The set of 149 genes were significantly enriched in the ER+ tumors by Gene Set Enrichment Analysis. Part of chromosome 8 was deleted in all ER- cells and part of chromosome 17 amplified in all ER+ cells. These loci encoded 30 genes that were overexpressed in ER+ cells; however, only 9 of these genes were overexpressed in ER+ tumors. Many of the methylation or CNA affected genes were functionally significant and some are known to be associated with breast cancer outcome such as GATA3 and LYN. Conclusions: A global pattern of differential CpG island methylation influences expression of a cohort of genes that contribute to the transcriptome landscape of ER+ and ER- breast cancer cells and tumors. The role of gene amplification/deletion in defining the transcriptome landscape of ER+ and ER- cells appears to more modest, although several potentially significant genes appear to be globally regulated by copy number aberrations. 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 4908. doi:10.1158/1538-7445.AM2011-4908

Abstract Epithelial ovarian cancer (EOC) is the 6th most commonly diagnosed cancer among women in the world, accounting for 4% of all female cancers. On a worldwide basis, an estimated 204,000 new cases are diagnosed and 125,000 women die of EOC annually, with 5-year survival rate at 45%. Metastasis and its complications are the usual cause of EOC-related deaths. Thus, factors related to tumor aggressiveness, response to therapy, and the underlying health state of the patient often predict survival. Candidate gene and genome-wide association studies (GWAS) for survival from EOC have discovered only a few the genetic variants, accounting for only a small proportion of the heritability. One explanation for the missing heritability is that the common analysis approach _ assessing the effect of each single nucleotide polymorphism (SNP) individually _ is not well suited to the detection of small effects of multiple SNPs. Gene set analysis (GSA) overcomes this limitation by assessing the overall evidence of association of a phenotype with all measured variation in a set of genes or pathway. Therefore, we conducted a GSA using data from an EOC GWAS, using the Illumina 610-quad SNP array, to determine novel gene sets (GSs) associated with overall survival (OS) from EOC (N cases = 1057, N deaths = 495). All subjects were non-Hispanic, non-Jewish Caucasians, with sample and SNP call rates >95%. Subjects with ambiguous gender, unresolved identical genotypes, and less than 80% European ancestry were also excluded. SNPs were mapped to genes within 20 kb, with genes then mapped to GSs defined by the following: KEGG, Gene Ontology (biological, cellular, molecular) and PharmGKB. GSA was completed in which the association of each gene with OS was assessed using a principal component gene-level test (adjusting for study site, age and population sub-structure based on the first eigenvector from EigenStrat), followed by the aggregation of the gene-level p-values using the Gamma method with a soft truncation threshold of 0.15. Empirical p-values were determined using permutation methods. Of the 2,573 GSs assessed, 30 GSs had a p-value < 0.01. The top GSs were related to: mismatch repair (p = 3.0E-4), platelet formation (p=0.001), leukocyte adhesive activation (p=0.002), antigen processing and presentation of endogenous antigens (p = 0.002) or peptide antigens (p = 0.002), platinum pathway (p = 0.003), aging and lifespan (p = 0.004, p = 0.007), immune system development (p = 0.004) and somatic diversification (p=0.007). These results illustrate the usefulness of GSA for the development of novel hypotheses that can be followed-up with additional clinical or functional studies. Moreover, application of GSA to survival from EOC allowed the detection of significant GSs (and the corresponding genes/SNPs within these GSs) that impact OS that were not detected by assessing each SNP individually, where the smallest p-value from single SNP analysis of the GWAS for OS was 1.3E-05. 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 2738. doi:10.1158/1538-7445.AM2011-2738

Abstract The etiology of ovarian cancer is poorly understood but there is clearly a heritable component. Efforts to identify susceptibility alleles rarely consider interactions among alleles or their joint effects because it quickly becomes computationally intractable. In this study, we sought to identify multi-SNP effects jointly with a pathway-based analysis (GSEA-SNP) of 1,952 EOC cases and 2,042 frequency-matched controls genotyped with the Illumina 610K array. All subjects were self-reported non-Hispanic non-Jewish Caucasians, with SNP and sample call rates > 95%. Subjects with ambiguous gender, unresolved identical genotypes and < 80% European ancestry were excluded. SNPs with MAF < 1% were excluded. Missing genotypes were inferred using Mach based on the HapMap CEU population. SNPs not within introns were annotated to genes within 100 bp. We retrieved the following databases of gene sets (GSs) from a compiled database, MsigDB: human chromosome and cytogenetic band, chemical and genetic perturbations, canonical pathways, microRNA binding targets, transcription factor targets (TFT), cancer gene neighborhood (CGN), Cancer modules and Gene Ontology. The GSEA-SNP approach calculates and rank orders the trend statistic for association between each SNP and EOC risk. The Enrichment Score (ES) estimates the overrepresentation of top-ranked SNPs for each GS; the statistical significance was estimated using 10,000 permutations. The ES was normalized according to the size of GS to yield Normalized Enrichment Score (NES). The false discovery (FDR) rate was estimated using NES to adjust for multiple hypothesis testing. A total of 5181 gene-sets were included in the analysis. When controlling the FDR at 15%, 14 of the GSs were highly enriched with association signals, including two chromosomal regions, 8q13 (p = 0.0055, FDR = 0.15) and 6p24 (p = 0.0012, FDR = 0.08). BIOCARTA_LYM_PATHWAY, a pathway related to cell adhesion and diapedesis of lymphocytes, was also significant (p = 0.0002, FDR = 0.12). A TFT GS, composed of genes with promoter regions around transcription start sites containing the motif GGCNRNWCTTYS, was associated with risk (p = 0.0004, FDR = 0.08). Currently, no known transcription factors bind to this computationally predicted motif. Molecular function of double stranded RNA binding was also enriched (p = 0.005, FDR = 0.10). The remaining enriched sets were CGN sets, for which the neighborhoods around the cancer genes were originally defined using correlation of gene expression from 4 large data sets of various cancer types. The enriched CGN included: CD48 (p = 0.011), CD53 (p = 0.009), CD97 (p = 0.011), INPP5D (p = 0.010), PTPN6 (p = 0.008), VAV1 (p = 0.011), ITGAL (p = 0.012), PTPRC (p = 0.012), and STAT6 (p = 0.010). In summary, these analyses detected biologically plausible GSs related to etiology of EOC, highlighting SNPs in core enrichment groups that were not identified using individual SNP tests. 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 4729. doi:10.1158/1538-7445.AM2011-4729