ISLET: individual-specific reference panel recovery improves cell-type-specific inference

We propose a statistical framework ISLET to infer individual-specific and cell-type-specific transcriptome reference panels. ISLET models the repeatedly measured bulk gene expression data, to optimize the usage of shared information within each subject. ISLET is the first available method to achieve individual-specific reference estimation in repeated samples. Using simulation studies, we show outstanding performance of ISLET in the reference estimation and downstream cell-type-specific differentially expressed genes testing. We apply ISLET to longitudinal transcriptomes profiled from blood samples in a large observational study of young children and confirm the cell-type-specific gene signatures for pancreatic islet autoantibody. ISLET is available at https://bioconductor.org/packages/ISLET .

Transcriptional networks in at-risk individuals identify signatures of type 1 diabetes progression

Type 1 diabetes (T1D) is a disease of insulin deficiency that results from autoimmune destruction of pancreatic islet β cells. The exact cause of T1D remains unknown, although asymptomatic islet autoimmunity lasting from weeks to years before diagnosis raises the possibility of intervention before the onset of clinical disease. The number, type, and titer of islet autoantibodies are associated with long-term disease risk but do not cause disease, and robust early predictors of individual progression to T1D onset remain elusive. The Environmental Determinants of Diabetes in the Young (TEDDY) consortium is a prospective cohort study aiming to determine genetic and environmental interactions causing T1D. Here, we analyzed longitudinal blood transcriptomes of 2013 samples from 400 individuals in the TEDDY study before both T1D and islet autoimmunity. We identified and interpreted age-associated gene expression changes in healthy infancy and age-independent changes tracking with progression to both T1D and islet autoimmunity, beginning before other evidence of islet autoimmunity was present. We combined multivariate longitudinal data in a Bayesian joint model to predict individual risk of T1D onset and validated the association of a natural killer cell signature with progression and the model's predictive performance on an additional 356 samples from 56 individuals in the independent Type 1 Diabetes Prediction and Prevention study. Together, our results indicate that T1D is characterized by early and longitudinal changes in gene expression, informing the immunopathology of disease progression and facilitating prediction of its course.

Longitudinal Metabolome-Wide Signals Prior to the Appearance of a First Islet Autoantibody in Children Participating in the TEDDY Study

Children at increased genetic risk for type 1 diabetes (T1D) after environmental exposures may develop pancreatic islet autoantibodies (IA) at a very young age. Metabolic profile changes over time may imply responses to exposures and signal development of the first IA. Our present research in The Environmental Determinants of Diabetes in the Young (TEDDY) study aimed to identify metabolome-wide signals preceding the first IA against GAD (GADA-first) or against insulin (IAA-first). We profiled metabolomes by mass spectrometry from children's plasma at 3-month intervals after birth until appearance of the first IA. A trajectory analysis discovered each first IA preceded by reduced amino acid proline and branched-chain amino acids (BCAAs), respectively. With independent time point analysis following birth, we discovered dehydroascorbic acid (DHAA) contributing to the risk of each first IA, and γ-aminobutyric acid (GABAs) associated with the first autoantibody against insulin (IAA-first). Methionine and alanine, compounds produced in BCAA metabolism and fatty acids, also preceded IA at different time points. Unsaturated triglycerides and phosphatidylethanolamines decreased in abundance before appearance of either autoantibody. Our findings suggest that IAA-first and GADA-first are heralded by different patterns of DHAA, GABA, multiple amino acids, and fatty acids, which may be important to primary prevention of T1D.

Plasma ascorbic acid and the risk of islet autoimmunity and type 1 diabetes: the TEDDY study

AIMS/HYPOTHESIS

We studied the association of plasma ascorbic acid with the risk of developing islet autoimmunity and type 1 diabetes and examined whether SNPs in vitamin C transport genes modify these associations. Furthermore, we aimed to determine whether the SNPs themselves are associated with the risk of islet autoimmunity or type 1 diabetes.

METHODS

We used a risk set sampled nested case-control design within an ongoing international multicentre observational study: The Environmental Determinants of Diabetes in the Young (TEDDY). The TEDDY study followed children with increased genetic risk from birth to endpoints of islet autoantibodies (350 cases, 974 controls) and type 1 diabetes (102 cases, 282 controls) in six clinical centres. Control participants were matched for family history of type 1 diabetes, clinical centre and sex. Plasma ascorbic acid concentration was measured at ages 6 and 12 months and then annually up to age 6 years. SNPs in vitamin C transport genes were genotyped using the ImmunoChip custom microarray. Comparisons were adjusted for HLA genotypes and for background population stratification.

RESULTS

Childhood plasma ascorbic acid (mean ± SD 10.76 ± 3.54 mg/l in controls) was inversely associated with islet autoimmunity risk (adjusted OR 0.96 [95% CI 0.92, 0.99] per +1 mg/l), particularly islet autoimmunity, starting with insulin autoantibodies (OR 0.94 [95% CI 0.88, 0.99]), but not with type 1 diabetes risk (OR 0.93 [95% Cl 0.86, 1.02]). The SLC2A2 rs5400 SNP was associated with increased risk of type 1 diabetes (OR 1.77 [95% CI 1.12, 2.80]), independent of plasma ascorbic acid (OR 0.92 [95% CI 0.84, 1.00]).

CONCLUSIONS/INTERPRETATION

Higher plasma ascorbic acid levels may protect against islet autoimmunity in children genetically at risk for type 1 diabetes. Further studies are warranted to confirm these findings.

DATA AVAILABILITY

The datasets generated and analysed during the current study will be made available in the NIDDK Central Repository at https://www.niddkrepository.org/studies/teddy.

An open resource for accurately benchmarking small variant and reference calls

Benchmark small variant calls are required for developing, optimizing and assessing the performance of sequencing and bioinformatics methods. Here, as part of the Genome in a Bottle (GIAB) Consortium, we apply a reproducible, cloud-based pipeline to integrate multiple short- and linked-read sequencing datasets and provide benchmark calls for human genomes. We generate benchmark calls for one previously analyzed GIAB sample, as well as six genomes from the Personal Genome Project. These new genomes have broad, open consent, making this a 'first of its kind' resource that is available to the community for multiple downstream applications. We produce 17% more benchmark single nucleotide variations, 176% more indels and 12% larger benchmark regions than previously published GIAB benchmarks. We demonstrate that this benchmark reliably identifies errors in existing callsets and highlight challenges in interpreting performance metrics when using benchmarks that are not perfect or comprehensive. Finally, we identify strengths and weaknesses of callsets by stratifying performance according to variant type and genome context.

Serum biomarkers of glucocorticoid response and safety in anti-neutrophil cytoplasmic antibody-associated vasculitis and juvenile dermatomyositis

Glucocorticoids are standard of care for many chronic inflammatory conditions, including juvenile dermatomyositis (JDM) and anti-neutrophil cytoplasmic antibody (ANCA)-associated vasculitis (AAV). We sought to define pharmacodynamic biomarkers of therapeutic efficacy and safety concerns of glucocorticoid treatment for these two disorders. Previous proteomic profiling of patients with Duchenne muscular dystrophy (DMD) and inflammatory bowel disease (IBD) treated with glucocorticoids identified candidate biomarkers for efficacy and safety concerns of glucocorticoids. Serial serum samples from patients with AAV (n = 30) and JDM (n = 12) were obtained during active disease, and after treatment with glucocorticoids. For AAV, 8 of 11 biomarkers of the anti-inflammatory response to glucocorticoids were validated (P-value ≤0.05; CD23, macrophage-derived cytokine, interleukin-22 binding protein, matrix metalloproteinase-12, T lymphocyte surface antigen Ly9, fibrinogen gamma chain, angiopoietin-2 [all decreased], and protein C [increased]), as were 5 of 7 safety biomarkers (P-value ≤0.05; afamin, matrix metalloproteinase-3, insulin growth factor binding protein-5, angiotensinogen, leptin [all increased]). For JDM, 10 of 11 efficacy biomarkers were validated (P-value ≤0.05; all proteins except fibrinogen gamma chain) and 6 of 7 safety biomarkers (P-value ≤0.05; AAV proteins plus growth hormone binding protein). The identified efficacy biomarkers may be useful as objective outcome measures for early phase proof-of-concept studies when assessing novel anti-inflammatory drugs in JDM and AAV, and likely in other inflammatory disorders. Similarly, safety biomarkers may also be helpful assessing toxicity of alternatives to glucocorticoids.

Metabolic pathways and immunometabolism in rare kidney diseases

OBJECTIVES

To characterise renal tissue metabolic pathway gene expression in different forms of glomerulonephritis.

METHODS

Patients with nephrotic syndrome (NS), antineutrophil cytoplasmic antibody-associated vasculitis (AAV), systemic lupus erythematosus (SLE) and healthy living donors (LD) were studied. Clinically indicated renal biopsies were obtained at time of diagnosis and microdissected into glomerular and tubulointerstitial compartments. Microarray-derived differential gene expression of 88 genes representing critical enzymes of metabolic pathways and 25 genes related to immune cell markers was compared between disease groups. Correlation analyses measured relationships between metabolic pathways, kidney function and cytokine production.

RESULTS

Reduced steady state levels of mRNA species were enriched in pathways of oxidative phosphorylation and increased in the pentose phosphate pathway (PPP) with maximal perturbation in AAV and SLE followed by NS, and least in LD. Transcript regulation was isozymes specific with robust regulation in hexokinases, enolases and glucose transporters. Intercorrelation networks were observed between enzymes of the PPP (eg, transketolase) and macrophage markers (eg, CD68) (r=0.49, p<0.01). Increased PPP transcript levels were associated with reduced glomerular filtration rate in the glomerular (r=-0.49, p<0.01) and tubulointerstitial (r=-0.41, p<0.01) compartments. PPP expression and tumour necrosis factor activation were tightly co-expressed (r=0.70, p<0.01).

CONCLUSION

This study demonstrated concordant alterations of the renal transcriptome consistent with metabolic reprogramming across different forms of glomerulonephritis. Activation of the PPP was tightly linked with intrarenal macrophage marker expression, reduced kidney function and increased production of cytokines. Modulation of glucose metabolism may offer novel immune-modulatory therapeutic approaches in rare kidney diseases.

Characterising cis-regulatory variation in the transcriptome of histologically normal and tumour-derived pancreatic tissues

OBJECTIVE

To elucidate the genetic architecture of gene expression in pancreatic tissues.

DESIGN

We performed expression quantitative trait locus (eQTL) analysis in histologically normal pancreatic tissue samples (n=95) using RNA sequencing and the corresponding 1000 genomes imputed germline genotypes. Data from pancreatic tumour-derived tissue samples (n=115) from The Cancer Genome Atlas were included for comparison.

RESULTS

We identified 38 615 cis-eQTLs (in 484 genes) in histologically normal tissues and 39 713 cis-eQTL (in 237 genes) in tumour-derived tissues (false discovery rate <0.1), with the strongest effects seen near transcriptional start sites. Approximately 23% and 42% of genes with significant cis-eQTLs appeared to be specific for tumour-derived and normal-derived tissues, respectively. Significant enrichment of cis-eQTL variants was noted in non-coding regulatory regions, in particular for pancreatic tissues (1.53-fold to 3.12-fold, p≤0.0001), indicating tissue-specific functional relevance. A common pancreatic cancer risk locus on 9q34.2 (rs687289) was associated with ABO expression in histologically normal (p=5.8×10^-8) and tumour-derived (p=8.3×10^-5) tissues. The high linkage disequilibrium between this variant and the O blood group generating deletion variant in ABO (exon 6) suggested that nonsense-mediated decay (NMD) of the 'O' mRNA might explain this finding. However, knockdown of crucial NMD regulators did not influence decay of the ABO 'O' mRNA, indicating that a gene regulatory element influenced by pancreatic cancer risk alleles may underlie the eQTL.

CONCLUSIONS

We have identified cis-eQTLs representing potential functional regulatory variants in the pancreas and generated a rich data set for further studies on gene expression and its regulation in pancreatic tissues.

A common intronic variant of PARP1 confers melanoma risk and mediates melanocyte growth via regulation of MITF

Previous genome-wide association studies have identified a melanoma-associated locus at 1q42.1 that encompasses a ∼100-kb region spanning the PARP1 gene. Expression quantitative trait locus (eQTL) analysis in multiple cell types of the melanocytic lineage consistently demonstrated that the 1q42.1 melanoma risk allele (rs3219090[G]) is correlated with higher PARP1 levels. In silico fine-mapping and functional validation identified a common intronic indel, rs144361550 (-/GGGCCC; r2 = 0.947 with rs3219090), as displaying allele-specific transcriptional activity. A proteomic screen identified RECQL as binding to rs144361550 in an allele-preferential manner. In human primary melanocytes, PARP1 promoted cell proliferation and rescued BRAFV600E-induced senescence phenotypes in a PARylation-independent manner. PARP1 also transformed TERT-immortalized melanocytes expressing BRAFV600E. PARP1-mediated senescence rescue was accompanied by transcriptional activation of the melanocyte-lineage survival oncogene MITF, highlighting a new role for PARP1 in melanomagenesis.

Abstract 1442: Analysis of cis-eQTLs in normal and tumor-derived pancreatic tissues reveals functional insights, including for the 9q34.1 ABO pancreatic cancer risk locus

Objective: To elucidate the genetic architecture of gene expression in pancreatic tissues.

Design: We performed expression quantitative trait locus (eQTL) and allele specific expression (ASE) analyses using RNA-sequence data and 1000 Genomes (1000G) imputed GWAS genotypes from 95 fresh frozen histologically normal pancreatic tissue samples. Data from 115 pancreatic tumor-derived tissue samples from The Cancer Genome Atlas (TCGA) was included for comparison.

Results: We identified 38,615 cis-eQTLs (corresponding to 484 Genes) in histologically normal tissues and 39,713 cis-eQTL (corresponding to 237 Genes) in tumor tissues (FDR&lt;0.1), with the strongest effects seen near transcriptional start sites (TSS). Approximately 23% and 42% of genes with significant cis-eQTLs (eGenes) appeared to be specific for tumor and normal derived tissues, respectively. Significant enrichment of cis-eQTL variants was noted in noncoding regulatory regions marked by modified histones, DNAse hypersensitivity and bound transcription factors, in particular for pancreatic tissues (1.53-3.12 fold, P≤0.0001), indicating tissue-specific functional relevance. A common pancreatic cancer risk locus on 9q34.2 in the ABO gene (rs687289) was associated with ABO expression in histologically normal (P=5.8x10-8) and tumor-derived (P=8.3x10-5) pancreatic tissues. The high linkage disequilibrium (LD) between this variant and the O blood group generating deletion variant in exon 6 of ABO suggested that nonsense-mediated decay (NMD) of the “O” mRNA could explain the eQTL. However, knock-down of crucial NMD regulators did not influence decay of the ABO “O” mRNA, indicating that a gene regulatory element influenced by pancreatic cancer risk alleles may underlie the eQTL.

Conclusions: We have identified cis-eQTLs representing potential functional regulatory variants in the pancreas and generated a rich dataset for further studies on gene expression and regulation in pancreatic tissues.

Citation Format: Laufey T. Amundadottir, Soren Lykke Andresen, Wenming Xiao, Jason Hoskins, Ashley Jermusyk, Lauren Rost, Irene Collins, Jinping Jia, Michael Mobaraki, Bin Zhu, Robert Kurtz, Hemang Parikh, Lei Song, Meredith Yeager, Torben Jensen, William Bamlet, Nilanjan Chatterjee, Brian Wolpin, Jill Smith, Sara Olson, Gloria Petersen, Jianxin Shi, Mingfeng Zhang. Analysis of cis-eQTLs in normal and tumor-derived pancreatic tissues reveals functional insights, including for the 9q34.1 ABO pancreatic cancer risk locus [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 1442. doi:10.1158/1538-7445.AM2017-1442

Functional characterization of a chr13q22.1 pancreatic cancer risk locus reveals long-range interaction and allele-specific effects on DIS3 expression

Genome-wide association studies (GWAS) have identified multiple common susceptibility loci for pancreatic cancer. Here we report fine-mapping and functional analysis of one such locus residing in a 610 kb gene desert on chr13q22.1 (marked by rs9543325). The closest candidate genes, KLF5, KLF12, PIBF1, DIS3 and BORA, range in distance from 265-586 kb. Sequencing three sub-regions containing the top ranked SNPs by imputation P-value revealed a 30 bp insertion/deletion (indel) variant that was significantly associated with pancreatic cancer risk (rs386772267, P = 2.30 × 10-11, OR = 1.22, 95% CI 1.15-1.28) and highly correlated to rs9543325 (r2 = 0.97 in the 1000 Genomes EUR population). This indel was the most significant cis-eQTL variant in a set of 222 histologically normal pancreatic tissue samples (β = 0.26, P = 0.004), with the insertion (risk-increasing) allele associated with reduced DIS3 expression. DIS3 encodes a catalytic subunit of the nuclear RNA exosome complex that mediates RNA processing and decay, and is mutated in several cancers. Chromosome conformation capture revealed a long range (570 kb) physical interaction between a sub-region of the risk locus, containing rs386772267, and a region ∼6 kb upstream of DIS3 Finally, repressor regulatory activity and allele-specific protein binding by transcription factors of the TCF/LEF family were observed for the risk-increasing allele of rs386772267, indicating that expression regulation at this risk locus may be influenced by the Wnt signaling pathway. In conclusion, we have identified a putative functional indel variant at chr13q22.1 that associates with decreased DIS3 expression in carriers of pancreatic cancer risk-increasing alleles, and could therefore affect nuclear RNA processing and/or decay.

Abstract 4487: An INDEL variant confers melanoma risk through PARP1 expression regulation

Recent genome wide association studies identified several new loci for melanoma susceptibility including chr1q42.1 locus encompassing Poly [ADP-ribose] polymerase 1 (PARP1). As an effort to identify effecter genes and functional risk variants from this locus we performed expression quantitative trait loci (eQTL) analysis in 62 melanoma cell lines. Among the genes in 1Mb area eQTL was identified for PARP1 and subsequently validated by qPCR, where increased PARP1 levels are significantly correlated with the risk allele (p = 0.03, copy number adjusted). Further allele discrimination qPCR of PARP1 transcripts in 14 melanoma cell lines heterozygous and of neutral copy for the lead SNP indicated significantly higher proportion for the risk allele (p = .0001). Same allelic imbalance was also observed in 51 heterozygous melanomas with neutral copy numbers from The Cancer Genome Atlas (p = .028). To identify functional risk variants mediating these effects we annotated this locus using six melanoma relevant cell types available from ENCODE and Roadmap database. Based on recent fine mapping data suggesting single-SNP model for this locus we prioritized high LD variants for nomination. Among 65 SNPs of high LD with the lead SNP or imputed best SNP (r2&gt;0.6 using 1000 Genomes phase3), four exhibited strong evidence as potential transcriptional enhancers. Electro Mobility Shift Assays and luciferase assays on these four variants demonstrated that one of them, a six-base pair INDEL (-/GGGCCC) in the first intron, displayed strong allelic functionality. Consistent with the expression data melanoma-associated deletion allele results in higher luciferase activity while protective insertion allele binds a group of proteins with higher affinity. Subsequent comparative mass-spectrometry for these insertion-binding proteins identified a striking collection of Guanine-quadruplex binding proteins including RECQ1 as potential inhibitors of PARP1 expression. Consistent with this hypothesis over-expression of RECQL results in more pronounced allelic difference in luciferase activities indicating that RECQ1 contributes to allelic PARP1 expression. Further interrogation of RECQ1 and PARP1 expression correlation analysis suggested that RECQ1 levels are inversely correlated with PARP1 in melanomas carrying insertion alleles. These data demonstrate that increased PARP1 expression is correlated with melanoma risk and an INDEL variant mediates differential PARP1 expression possibly through secondary DNA structure binding proteins including RECQ1.

Citation Format: Jiyeon Choi, Matthew Makowski, Tongwu Zhang, Matthew Law, Wendy Kim, Michael Kovacs, Hemang Parikh, Lauren Aoude, Michael Gartside, Jeffrey Trent, Michiel Vermeulen, Stuart Macgregor, Nicholas Hayward, Mai Xu, Kevin Brown. An INDEL variant confers melanoma risk through PARP1 expression regulation. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 4487.

svclassify: a method to establish benchmark structural variant calls

Background

The human genome contains variants ranging in size from small single nucleotide polymorphisms (SNPs) to large structural variants (SVs). High-quality benchmark small variant calls for the pilot National Institute of Standards and Technology (NIST) Reference Material (NA12878) have been developed by the Genome in a Bottle Consortium, but no similar high-quality benchmark SV calls exist for this genome. Since SV callers output highly discordant results, we developed methods to combine multiple forms of evidence from multiple sequencing technologies to classify candidate SVs into likely true or false positives. Our method (svclassify) calculates annotations from one or more aligned bam files from many high-throughput sequencing technologies, and then builds a one-class model using these annotations to classify candidate SVs as likely true or false positives.

Results

We first used pedigree analysis to develop a set of high-confidence breakpoint-resolved large deletions. We then used svclassify to cluster and classify these deletions as well as a set of high-confidence deletions from the 1000 Genomes Project and a set of breakpoint-resolved complex insertions from Spiral Genetics. We find that likely SVs cluster separately from likely non-SVs based on our annotations, and that the SVs cluster into different types of deletions. We then developed a supervised one-class classification method that uses a training set of random non-SV regions to determine whether candidate SVs have abnormal annotations different from most of the genome. To test this classification method, we use our pedigree-based breakpoint-resolved SVs, SVs validated by the 1000 Genomes Project, and assembly-based breakpoint-resolved insertions, along with semi-automated visualization using svviz.

Conclusions

We find that candidate SVs with high scores from multiple technologies have high concordance with PCR validation and an orthogonal consensus method MetaSV (99.7 % concordant), and candidate SVs with low scores are questionable. We distribute a set of 2676 high-confidence deletions and 68 high-confidence insertions with high svclassify scores from these call sets for benchmarking SV callers. We expect these methods to be particularly useful for establishing high-confidence SV calls for benchmark samples that have been characterized by multiple technologies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2366-2) contains supplementary material, which is available to authorized users.

Abstract A1-09: Functional analysis of the chr13q22.1 pancreatic cancer risk locus suggests allele-specific effects on DIS3 expression with prognostic implications

A genome-wide association study (GWAS) of pancreatic cancer conducted within the NCI Cohort Consortium (PanScan I and II) identified pancreatic cancer susceptibility loci on chromosomes 1q32.1/NR5A2, 5p15.33/CLPTM1L/TERT, 9q34.2/ABO, and 13q22.1. The most significant single-nucleotide polymorphism (SNP) identified on 13q22.1, rs9543325, lies in a ~600 kb gene desert; the nearest genes are KLF5, KLF12, PIBF1, DIS3, and BORA (265–586 kb away). Imputation using the 1000 Genomes and DCEG reference datasets did not improve the GWAS signal, but produced a set of highly correlated SNPs for functional follow-up. We performed eQTL analysis to test for association between the genotypes of these functional candidate variants and expression of nearby genes. Among 64 normal derived pancreatic tissue samples, DIS3 expression showed the strongest association with a 30 bp indel variant in the risk locus (P = 4.8 × 10-4), indicating risk alleles associate with reduced DIS3 expression. Mutations in DIS3 have previously been identified in acute myeloid leukemia and multiple myeloma, and its expression has been correlated with metastatic potential in colorectal cancer, suggesting DIS3 is relevant to cancer biology. Chromosome conformation capture identified a physical interaction between the indel-containing locus and a region near the DIS3 promoter. Luciferase assay for regulatory function of this indel-containing locus revealed allele-specific silencer activity for the insertion allele. Supershift electromobility shift assay (EMSA) demonstrated binding of LEF1 specifically to the insertion allele of the indel, which contains two in silico predicted LEF1 binding elements. Finally, through immunohistochemical analysis, high DIS3 protein levels associated with better survival for pancreatic cancer patients (hazard ratio = 2.87, 95% CI = 1.49–5.53, P = 0.001). Our results suggest that at least one target gene for the pancreatic cancer risk variants on chr13q22.1 may be DIS3, and that the underlying biology may be mediated by the novel indel through a long-range repressive effect on DIS3 expression.

Citation Format: Jason W. Hoskins, Abdisamad Ibrahim, Mickey Emmanuel, Sarah Manmiller, Jinping Jia, Hemang Parikh, Irene Collins, Kris Ylaya, Sean F. Altekruse, Stephen M. Hewitt, Gloria M. Petersen, Laufey T. Amundadottir. Functional analysis of the chr13q22.1 pancreatic cancer risk locus suggests allele-specific effects on DIS3 expression with prognostic implications. [abstract]. In: Proceedings of the AACR Special Conference on Translation of the Cancer Genome; Feb 7-9, 2015; San Francisco, CA. Philadelphia (PA): AACR; Cancer Res 2015;75(22 Suppl 1):Abstract nr A1-09.

Less pronounced response to exercise in healthy relatives to type 2 diabetic subjects compared with controls

Healthy first-degree relatives with heredity of type 2 diabetes (FH+) are known to have metabolic inflexibility compared with subjects without heredity for diabetes (FH-). In this study, we aimed to test the hypothesis that FH+ individuals have an impaired response to exercise compared with FH-. Sixteen FH+ and 19 FH- insulin-sensitive men similar in age, peak oxygen consumption (V̇o2 peak), and body mass index completed an exercise intervention with heart rate monitored during exercise for 7 mo. Before and after the exercise intervention, the participants underwent a physical examination and tests for glucose tolerance and exercise capacity, and muscle biopsies were taken for expression analysis. The participants attended, on average, 39 training sessions during the intervention and spent 18.8 MJ on exercise. V̇o2 peak/kg increased by 14%, and the participants lost 1.2 kg of weight and 3 cm waist circumference. Given that the FH+ group expended 61% more energy during the intervention, we used regression analysis to analyze the response in the FH+ and FH- groups separately. Exercise volume had a significant effect on V̇o2 peak, weight, and waist circumference in the FH- group, but not in the FH+ group. After exercise, expression of genes involved in metabolism, oxidative phosphorylation, and cellular respiration increased more in the FH- compared with the FH+ group. This suggests that healthy, insulin-sensitive FH+ and FH- participants with similar age, V̇o2 peak, and body mass index may respond differently to an exercise intervention. The FH+ background might limit muscle adaptation to exercise, which may contribute to the increased susceptibility to type 2 diabetes in FH+ individuals.

Abstract B111: Functional analysis of the chr13q22.1 pancreatic cancer risk locus suggests allele-specific effects on DIS3 expression with prognostic implications

A genome-wide association study (GWAS) of pancreatic cancer conducted within the NCI Cohort Consortium (PanScan I and II) identified pancreatic cancer susceptibility loci on chromosomes 1q32.1/NR5A2, 5p15.33/CLPTM1L/TERT, 9q34.2/ABO, and 13q22.1. The most significant single-nucleotide polymorphism (SNP) identified on 13q22.1, rs9543325, lies in a 600 kb gene desert; the nearest genes are KLF5, KLF12, PIBF1, DIS3, and BORA (265–586 kb). Imputation using the 1000 Genomes and DCEG reference datasets did not improve the GWAS signal, but produced a set of highly correlated SNPs for functional follow-up. We performed eQTL analysis to test for association between the genotypes of these functional candidate variants and expression of nearby genes. Among 100 normal derived pancreatic tissue samples, DIS3 expression showed the strongest association with a novel 30 bp indel variant in the risk locus (P = 4.0 × 104), indicating risk alleles associate with reduced DIS3 expression. Mutations in DIS3 have previously been identified in acute myeloid leukemia and multiple myeloma, and its expression has been correlated with metastatic potential in colorectal cancer, suggesting DIS3 is relevant to cancer biology. Chromosome conformation capture identified a physical interaction between the indel-containing locus and a region near the DIS3 promoter. Luciferase assay for regulatory function of this indel-containing locus revealed allele-specific silencer activity for the insertion allele. Supershift electromobility shift assay (EMSA) demonstrated binding of LEF1 specifically to the insertion allele of the indel, which contains two in silico predicted LEF1 binding elements. Finally, through immunohistochemical analysis, high DIS3 protein levels associated with better survival for pancreatic cancer patients (hazard ratio = 2.87, 95% CI = 1.49–5.53, P = 0.001). Our results suggest that at least one target gene for the pancreatic cancer risk variants on chr13q22.1 may be DIS3, and that the underlying biology may be mediated by the novel indel through a long-range repressive effect on DIS3 expression.

Citation Format: Jason W. Hoskins, Abdisamad Ibrahim, Mickey Emmanuel, Hemang Parikh, Jinping Jia, Irene Collins, Kris Ylaya, Sean F. Altekruse, Stephen M. Hewitt, Gloria M. Petersen, Laufey T. Amundadottir. Functional analysis of the chr13q22.1 pancreatic cancer risk locus suggests allele-specific effects on DIS3 expression with prognostic implications. [abstract]. In: Proceedings of the AACR Special Conference on Pancreatic Cancer: Innovations in Research and Treatment; May 18-21, 2014; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2015;75(13 Suppl):Abstract nr B111.

Abstract B112: Transcriptome analysis in pancreatic cancer reveals a tumor suppressor function for HNF1A

Pancreatic ductal adenocarcinoma (PDAC), as with all cancer, is driven by dysregulation of genetic programs due to accumulation of somatic mutations, epigenetic modifications and changes in the micro-environment. Holistic views of perturbations in gene expression networks could illuminate key regulators and pathways for this highly lethal cancer. Toward this end, we performed massively-parallel mRNA-sequencing in normal (n = 10) and tumor (n = 8) derived pancreatic tissue samples as well as pancreatic cancer cell lines (n = 9), and determined differential gene expression (DE) patterns. Sub-network enrichment analyses of all expressed genes based on magnitude of DE identified HNF1A as the regulator of the most significantly and consistently dysregulated expression sub-network in our pancreatic tumor samples (median P = 7.56 x 10-7, median rank = 1, range = 1-25). To explore the effect of HNF1A expression in pancreatic tumor-derived cells we generated stable HNF1A-inducible clones in two pancreatic cancer cell lines. Induction of HNF1A overexpression caused severe growth inhibition (5.3-fold, P = 4.5x10-5 for MIA PaCa-2 clones; 7.2-fold, P = 2.2x10-5 for PANC-1 clones), and G0/G1 cell cycle arrest. This was accompanied by down-regulation of 49 out of 84 assayed cell cycle genes, while only 1 displayed increased expression. These data, combined with HNF1A’s direct regulation of pancreatic development and homeostasis genes (i.e. PDX1, PTF1A and NR5A2), suggest it may be an important tumor suppressor in pancreatic cells.

Citation Format: Jason W. Hoskins, Marta Flandez, Jinping Jia, Hemang Parikh, Irene Collins, Mickey Emmanuel, Abdisamad Ibrahim, Wenming Xiao, John Powell, Nuria Malats, Gloria M. Petersen, Fransisco X. Real, Laufey T. Amundadottir. Transcriptome analysis in pancreatic cancer reveals a tumor suppressor function for HNF1A. [abstract]. In: Proceedings of the AACR Special Conference on Pancreatic Cancer: Innovations in Research and Treatment; May 18-21, 2014; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2015;75(13 Suppl):Abstract nr B112.

Imputation and subset-based association analysis across different cancer types identifies multiple independent risk loci in the TERT-CLPTM1L region on chromosome 5p15.33

Genome-wide association studies (GWAS) have mapped risk alleles for at least 10 distinct cancers to a small region of 63 000 bp on chromosome 5p15.33. This region harbors the TERT and CLPTM1L genes; the former encodes the catalytic subunit of telomerase reverse transcriptase and the latter may play a role in apoptosis. To investigate further the genetic architecture of common susceptibility alleles in this region, we conducted an agnostic subset-based meta-analysis (association analysis based on subsets) across six distinct cancers in 34 248 cases and 45 036 controls. Based on sequential conditional analysis, we identified as many as six independent risk loci marked by common single-nucleotide polymorphisms: five in the TERT gene (Region 1: rs7726159, P = 2.10 × 10(-39); Region 3: rs2853677, P = 3.30 × 10(-36) and PConditional = 2.36 × 10(-8); Region 4: rs2736098, P = 3.87 × 10(-12) and PConditional = 5.19 × 10(-6), Region 5: rs13172201, P = 0.041 and PConditional = 2.04 × 10(-6); and Region 6: rs10069690, P = 7.49 × 10(-15) and PConditional = 5.35 × 10(-7)) and one in the neighboring CLPTM1L gene (Region 2: rs451360; P = 1.90 × 10(-18) and PConditional = 7.06 × 10(-16)). Between three and five cancers mapped to each independent locus with both risk-enhancing and protective effects. Allele-specific effects on DNA methylation were seen for a subset of risk loci, indicating that methylation and subsequent effects on gene expression may contribute to the biology of risk variants on 5p15.33. Our results provide strong support for extensive pleiotropy across this region of 5p15.33, to an extent not previously observed in other cancer susceptibility loci.

Smoking induces overexpression of immediate early genes in active Graves' ophthalmopathy

BACKGROUND

Cigarette smoking is a risk factor for the development of Graves' ophthalmopathy (GO). In a previous study of gene expression in intraorbital fat, adipocyte-related immediate early genes (IEGs) were overexpressed in patients with GO compared to controls. We investigated whether IEGs are upregulated by smoking, and examined other pathways that may be affected by smoking.

METHODS

Gene expression in intraorbital fat was studied in smokers (n=8) and nonsmokers (n=8) with severe active GO, as well as in subcutaneous fat in thyroid-healthy smokers (n=5) and nonsmokers (n=5) using microarray and real-time polymerase chain reaction (PCR).

RESULTS

With microarray, eight IEGs were upregulated more than 1.5-fold in smokers compared to nonsmokers with GO. Five were chosen for confirmation and were also overexpressed with real-time PCR. Interleukin-1 beta/IL-1B/(2.3-fold) and interleukin-6/IL-6/(2.4-fold) were upregulated both with microarray and with real-time PCR in smokers with GO compared to nonsmokers. Major histocompatibility complex, class II, DR beta 1/HLA-DRB1/was upregulated with microarray (2.1-fold) and with borderline significance with real-time PCR. None of these genes were upregulated in smokers compared to nonsmokers in subcutaneous fat.

CONCLUSIONS

IEGs, IL-1B, and IL-6 were overexpressed in smokers with severe active GO compared to nonsmokers, suggesting that smoking activates pathways associated with adipogenesis and inflammation. This study underlines the importance of IEGs in the pathogenesis of GO, and provides evidence for possible novel therapeutic interventions in GO. The mechanisms activated by smoking may be shared with other conditions such as rheumatoid arthritis.

Transcriptome analysis of pancreatic cancer reveals a tumor suppressor function for HNF1A

Pancreatic ductal adenocarcinoma (PDAC) is driven by the accumulation of somatic mutations, epigenetic modifications and changes in the micro-environment. New approaches to investigating disruptions of gene expression networks promise to uncover key regulators and pathways in carcinogenesis. We performed messenger RNA-sequencing in pancreatic normal (n = 10) and tumor (n = 8) derived tissue samples, as well as in pancreatic cancer cell lines (n = 9), to determine differential gene expression (DE) patterns. Sub-network enrichment analyses identified HNF1A as the regulator of the most significantly and consistently dysregulated expression sub-network in pancreatic tumor tissues and cells (median P = 7.56×10(-7), median rank = 1, range = 1-25). To explore the effects of HNF1A expression in pancreatic tumor-derived cells, we generated stable HNF1A-inducible clones in two pancreatic cancer cell lines (PANC-1 and MIA PaCa-2) and observed growth inhibition (5.3-fold, P = 4.5×10(-5) for MIA PaCa-2 clones; 7.2-fold, P = 2.2×10(-5) for PANC-1 clones), and a G0/G1 cell cycle arrest and apoptosis upon induction. These effects correlated with HNF1A-induced down-regulation of 51 of 84 cell cycle genes (e.g. E2F1, CDK2, CDK4, MCM2/3/4/5, SKP2 and CCND1), decreased expression of anti-apoptotic genes (e.g. BIRC2/5/6 and AKT) and increased expression of pro-apoptotic genes (e.g. CASP4/9/10 and APAF1). In light of the established role of HNF1A in the regulation of pancreatic development and homeostasis, our data suggest that it also functions as an important tumor suppressor in the pancreas.

Identification and validation of TGFBI as a promising prognosis marker of clear cell renal cell carcinoma

OBJECTIVE

To identify prognostic biomarkers in clear cell renal cell carcinoma (ccRCC) using a proteomic approach.

MATERIAL AND METHODS

We performed a comparative proteomic profiling of ccRCC and normal renal tissues from 9 different human specimens. We assessed differential protein expression by iTRAQ (isobaric tagging reagent for absolute quantify) labeling with regard to tumor aggressiveness according to the stage, size, grade, and necrosis (SSIGN) score and confirmed our results using Western blot (9 patients) and immunohistochemistry (135 patients) analysis.

RESULTS

After proteomic analysis, 928 constitutive proteins were identified. Among these proteins, 346 had a modified expression in tumor compared with that of normal tissue. Pathway and integrated analyses indicated the presence of an up-regulation of the pentose phosphate pathway in aggressive tumors. In total, 14 proteins were excreted and could potentially become biomarkers. Overexpression of transforming growth factor, beta-induced (TGFBI) in ccRCC was confirmed using Western blot and immunohistochemistry analysis. A significant association was found between the presence of TGFBI expression with tumor category T3-4 (P<0.0001), Fuhrman grades III and IV (P<0.0001), tumor size>4cm (P<0.0001), presence of tumor necrosis (P<0.0001), nodal involvement (n = 0.009), metastasis (P = 0.012), SSIGN score≥5 (P<0.0001), cancer progression (P<0.0001), and cancer-specific death (P<0.0001). Cancer-specific survival was significantly better for patients with no cytoplasmic TGFBI expression (1-, 3-, 5-y cancer-specific survival of 94.7%, 87.8%, and 73.4% vs. 92.9%, 71.2%, and 49.8%, respectively; P<0.0001).

CONCLUSION

We identified 346 proteins involved in renal carcinogenesis and confirmed the presence of a metabolic shift in aggressive tumors. TGFBI was overexpressed in tumors with high SSIGN scores and was significantly associated with oncologic outcomes.

CLPTM1L promotes growth and enhances aneuploidy in pancreatic cancer cells

Genome-wide association studies (GWAS) of 10 different cancers have identified pleiotropic cancer predisposition loci across a region of chromosome 5p15.33 that includes the TERT and CLPTM1L genes. Of these, susceptibility alleles for pancreatic cancer have mapped to the CLPTM1L gene, thus prompting an investigation of the function of CLPTM1L in the pancreas. Immunofluorescence analysis indicated that CLPTM1L localized to the endoplasmic reticulum where it is likely embedded in the membrane, in accord with multiple predicted transmembrane domains. Overexpression of CLPTM1L enhanced growth of pancreatic cancer cells in vitro (1.3-1.5-fold; PDAY7 < 0.003) and in vivo (3.46-fold; PDAY68 = 0.039), suggesting a role in tumor growth; this effect was abrogated by deletion of two hydrophilic domains. Affinity purification followed by mass spectrometry identified an interaction between CLPTM1L and non-muscle myosin II (NMM-II), a protein involved in maintaining cell shape, migration, and cytokinesis. The two proteins colocalized in the cytoplasm and, after treatment with a DNA-damaging agent, at the centrosomes. Overexpression of CLPTM1L and depletion of NMM-II induced aneuploidy, indicating that CLPTM1L may interfere with normal NMM-II function in regulating cytokinesis. Immunohistochemical analysis revealed enhanced staining of CLPTM1L in human pancreatic ductal adenocarcinoma (n = 378) as compared with normal pancreatic tissue samples (n = 17; P = 1.7 × 10(-4)). Our results suggest that CLPTM1L functions as a growth-promoting gene in the pancreas and that overexpression may lead to an abrogation of normal cytokinesis, indicating that it should be considered as a plausible candidate gene that could explain the effect of pancreatic cancer susceptibility alleles on chr5p15.33.

An integrated transcriptome and epigenome analysis identifies a novel candidate gene for pancreatic cancer

BACKGROUND

Pancreatic cancer is a highly lethal cancer with limited diagnostic and therapeutic modalities.

METHODS

To begin to explore the genomic landscape of pancreatic cancer, we used massively parallel sequencing to catalog and compare transcribed regions and potential regulatory elements in two human cell lines derived from normal and cancerous pancreas.

RESULTS

By RNA-sequencing, we identified 2,146 differentially expressed genes in these cell lines that were enriched in cancer related pathways and biological processes that include cell adhesion, growth factor and receptor activity, signaling, transcription and differentiation. Our high throughput Chromatin immunoprecipitation (ChIP) sequence analysis furthermore identified over 100,000 regions enriched in epigenetic marks, showing either positive (H3K4me1, H3K4me3, RNA Pol II) or negative (H3K27me3) correlation with gene expression. Notably, an overall enrichment of RNA Pol II binding and depletion of H3K27me3 binding were seen in the cancer derived cell line as compared to the normal derived cell line. By selecting genes for further assessment based on this difference, we confirmed enhanced expression of aldehyde dehydrogenase 1A3 (ALDH1A3) in two larger sets of pancreatic cancer cell lines and in tumor tissues as compared to normal derived tissues.

CONCLUSIONS

As aldehyde dehydrogenase (ALDH) activity is a key feature of cancer stem cells, our results indicate that a member of the ALDH superfamily, ALDH1A3, may be upregulated in pancreatic cancer, where it could mark pancreatic cancer stem cells.

Abstract 3137: Functional analysis of the chr13q22.1 pancreatic cancer risk locus suggests allele-specific effects on DIS3 expression

Pancreatic cancer is the 10th most common cancer and 4th most common cause of cancer mortality in the United States. A genome wide association study has revealed pancreatic cancer susceptibility loci on chromosomes 13q22.1, 1q32.1 and 5p15.33. The chr13q22.1 region contains the most significant of all risk SNPs, rs9543325, (P= 3.27 x 10−11) which is located in a gene desert. The nearest genes are KLF5, KLF12, PIBF1, DIS3 and BORA, which range in distance from 265kb to 586kb, respectively, from the most significant variants. Given the high linkage disequilibrium across this risk locus, there are many candidate functional variants to consider. Imputation in the region did not improve the signal but gave a set of highly correlated SNPs which likely includes the functional variant(s). In an effort to identify the functional variant(s) we performed eQTL analyses to test the association between the genotypes of these candidate SNPs and expression of nearby genes. Among 100 normal pancreatic tissue samples, DIS3 showed the strongest association with SNPs in our risk locus (P-values as low as 0.0004). Mutations in DIS3 have been identified in acute myeloid leukemia and multiple myeloma, and its expression has been correlated with metastatic potential in colorectal cancer, suggesting this gene could be important in pancreatic cancer biology. Chromosome Conformation Capture (3C) was performed to test for physical interactions between the risk locus and nearby genes. This assay confirmed the three dimensional proximity of the risk locus and DIS3 promoter. Sub-regions of the risk locus were then cloned upstream of a minimal promoter controlling luciferase expression to assay for potential allele specific enhancer/silencer activity. This assay revealed a sub-region downstream of rs9543325 that causes allele-specific silencing. These results suggest that a sub-region of the chr13q22.1 risk locus has allele-specific effects on DIS3 expression, which may have implications in the susceptibility to pancreatic cancer.

Citation Format: Jason Hoskins, Abdisamad Ibrahim, Jinping Jia, Irene Collins, Hemang Parikh, Gloria M. Petersen, Laufey Amundadottir. Functional analysis of the chr13q22.1 pancreatic cancer risk locus suggests allele-specific effects on DIS3 expression. [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 3137. doi:10.1158/1538-7445.AM2013-3137

Abstract 2556: Functional characterization of the pancreatic cancer TERT-CLPTM1L risk locus on chr5p15.33

Pancreatic cancer is a highly lethal cancer with few well established risk factors. PanScan, a genome wide association study (GWAS) of pancreatic cancer has identified four pancreatic cancer susceptibility loci in populations of European ancestry. One of these is in a multi cancer locus on chr5p15.33 where the most significant SNP from the GWAS (rs401681, P=3.7x10-7, ORAllele=1.19) lies in a region containing two genes, TERT and CLPTM1L. The TERT gene encodes the catalytic subunit of telomerase, well known for its essential role in maintaining telomere ends. The function of CLPTM1L is not as clear, although it has been proposed to play a role in apoptosis. It is predicted to encode a protein with 6 transmembrane (TM) domains and two large hydrophilic domains: a loop of 253 aa between the first and second TM domains, and a C-terminal tail of 89 aa.

We have performed imputation to fine-map the signal to a SNP three orders of magnitude more significant than the GWAS SNP (PImputed=1.4x10-10, ORAllele=1.30). As this SNP is located in the CLPTM1L gene we have performed a series of experiments to investigate the function of the CLPTM1L gene and its encoded protein. Immunofluorescence analysis in pancreatic cancer cells (PANC-1) indicates that it localizes to the endoplasmic reticulum. Affinity purification and mass spectrometry (HEK-293T, hTERT-HPNE and PANC-1 cells) identified MYH9, a non-muscle heavy chain myosin, as a potential interacting protein. The interaction has been validated by co-immunoprecipitation and co-localization experiments. To examine if CLPTM1L plays a role in growth control, we created stable PANC-1 cell lines overexpressing the full length CLPTM1L gene as well as two deletions, a C-terminal deletion and a loop deletion, and assayed growth in vitro and in vivo. Cell lines overexpressing full length CLPTM1L grow faster in vitro and in vivo as compared to cells containing empty vector. Interestingly, the two CLPTM1L mutants abolish this effect. Furthermore, we have shown by RNA-seq that the CLPTM1L gene is overexpressed in pancreatic tumors as compared to normal pancreatic tissues. Our results indicate that CLPTM1L may play a role in the control of cell growth and oncogenesis in the pancreas. Our current efforts aim at further characterizing the function of CLPTM1L and to correlate pancreatic cancer risk variants on 5p15.33 to molecular phenotypes to attempt to explain the underlying biology of the risk.

Citation Format: Jinping Jia, Irene Collins, Marta Dzyadyk, Abbey Thompson, Adam Cheuk, Hemang Parikh, Zhaoming Wang, Chris Westlake, Allen Bosley, Gloria Petersen, Thorkell Andresson, Laufey Amundadottir. Functional characterization of the pancreatic cancer TERT-CLPTM1L risk locus on chr5p15.33. [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 2556. doi:10.1158/1538-7445.AM2013-2556

Note: This abstract was not presented at the AACR Annual Meeting 2013 because the presenter was unable to attend.

Characterization of SNPs associated with prostate cancer in men of Ashkenazic descent from the set of GWAS identified SNPs: impact of cancer family history and cumulative SNP risk prediction

BACKGROUND

Genome-wide association studies (GWAS) have identified multiple SNPs associated with prostate cancer (PrCa). Population isolates may have different sets of risk alleles for PrCa constituting unique population and individual risk profiles.

METHODS

To test this hypothesis, associations between 31 GWAS SNPs of PrCa were examined among 979 PrCa cases and 1,251 controls of Ashkenazic descent using logistic regression. We also investigated risks by age at diagnosis, pathological features of PrCa, and family history of cancer. Moreover, we examined associations between cumulative number of risk alleles and PrCa and assessed the utility of risk alleles in PrCa risk prediction by comparing the area under the curve (AUC) for different logistic models.

RESULTS

Of the 31 genotyped SNPs, 8 were associated with PrCa at p ≤ 0.002 (corrected p-value threshold) with odds ratios (ORs) ranging from 1.22 to 1.42 per risk allele. Four SNPs were associated with aggressive PrCa, while three other SNPs showed potential interactions for PrCa by family history of PrCa (rs8102476; 19q13), lung cancer (rs17021918; 4q22), and breast cancer (rs10896449; 11q13). Men in the highest vs. lowest quartile of cumulative number of risk alleles had ORs of 3.70 (95% CI 2.76-4.97); 3.76 (95% CI 2.57-5.50), and 5.20 (95% CI 2.94-9.19) for overall PrCa, aggressive cancer and younger age at diagnosis, respectively. The addition of cumulative risk alleles to the model containing age at diagnosis and family history of PrCa yielded a slightly higher AUC (0.69 vs. 0.64).

CONCLUSION

These data define a set of risk alleles associated with PrCa in men of Ashkenazic descent and indicate possible genetic differences for PrCa between populations of European and Ashkenazic ancestry. Use of genetic markers might provide an opportunity to identify men at highest risk for younger age of onset PrCa; however, their clinical utility in identifying men at highest risk for aggressive cancer remains limited.

Impact of an exercise intervention on DNA methylation in skeletal muscle from first-degree relatives of patients with type 2 diabetes

To identify epigenetic patterns, which may predispose to type 2 diabetes (T2D) due to a family history (FH) of the disease, we analyzed DNA methylation genome-wide in skeletal muscle from individuals with (FH(+)) or without (FH(-)) an FH of T2D. We found differential DNA methylation of genes in biological pathways including mitogen-activated protein kinase (MAPK), insulin, and calcium signaling (P ≤ 0.007) and of individual genes with known function in muscle, including MAPK1, MYO18B, HOXC6, and the AMP-activated protein kinase subunit PRKAB1 in skeletal muscle of FH(+) compared with FH(-) men. We further validated our findings from FH(+) men in monozygotic twin pairs discordant for T2D, and 40% of 65 analyzed genes exhibited differential DNA methylation in muscle of both FH(+) men and diabetic twins. We further examined if a 6-month exercise intervention modifies the genome-wide DNA methylation pattern in skeletal muscle of the FH(+) and FH(-) individuals. DNA methylation of genes in retinol metabolism and calcium signaling pathways (P < 3 × 10(-6)) and with known functions in muscle and T2D including MEF2A, RUNX1, NDUFC2, and THADA decreased after exercise. Methylation of these human promoter regions suppressed reporter gene expression in vitro. In addition, both expression and methylation of several genes, i.e., ADIPOR1, BDKRB2, and TRIB1, changed after exercise. These findings provide new insights into how genetic background and environment can alter the human epigenome.

Abstract 28: Characterization of SNPs associated with prostate cancer in men of Ashkenazic descent from the set of GWAS identified SNPs: Impact of cancer family history

Background: Genome-wide association studies (GWAS) have identified several SNPs that are independently associated with prostate cancer (PrCa) in Caucasian men. Population isolates with unique allele frequencies are likely to have different sets of risk alleles for PrCa constituting personalized population risk stratification in addition to individual risk profiles.

Methods: To test this hypothesis, we examined associations between 31 GWAS SNPs and PrCa among 979 cases and 1,251 controls of Ashkenazic descent. Odds ratios (OR) and 95% confidence intervals (CIs) were estimated for each SNP using logistic regression models adjusted for age. We also investigated risk in strata by age at diagnosis, pathological features of PrCa, as well as family history of PrCa and other cancers.

Results: Of 31 SNPs, 19 (62%) were associated with PrCa at p&lt;0.05 in one or more genetic risk models. The strongest associations with PrCa were for rs6983267 at 8q24 (OR=1.78, p=5.7 x 10−7;comparison of GG vs. TT genotype), and rs1465618 at 2p21 (OR=2.19; p=0.003; AA vs. AG/GG). Six SNPs were associated with aggressive PrCa, while six others with less aggressive cancer. Four SNPs showed statistically significant interactions between PrCa and family history of PrCa (rs8102476 at 19q13), lung cancer (rs17021918 at 4q22), colorectal cancer (rs6983267 at 8q24), and breast cancer (rs10896449 at 11q13).

Conclusion: These data define a specific set of risk alleles associated with PrCa in Ashkenazic men. Results demonstrate differences between men of European and Ashkenazic ancestry in relation to genetic susceptibility of PrCa. Interactions of genetic variants with family history of specific cancers indicate a complex network of inherited cancer risk syndromes still to be defined. Use of genetic markers to identify men at higher risk of PrCa will provide an opportunity for personalized care of PrCa prevention through screening of high-risk individuals, independent of family history information.

Citation Format: Ilir Agalliu, Zhaoming Wang, Anne Dunn, Hemang Parikh, Timothy Myers, Robert Burk, Laufey Amundadottir. Characterization of SNPs associated with prostate cancer in men of Ashkenazic descent from the set of GWAS identified SNPs: Impact of cancer family history. [abstract]. In: Proceedings of the AACR Special Conference on Post-GWAS Horizons in Molecular Epidemiology: Digging Deeper into the Environment; 2012 Nov 11-14; Hollywood, FL. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2012;21(11 Suppl):Abstract nr 28.

A systems genetics approach identifies genes and pathways for type 2 diabetes in human islets

Close to 50 genetic loci have been associated with type 2 diabetes (T2D), but they explain only 15% of the heritability. In an attempt to identify additional T2D genes, we analyzed global gene expression in human islets from 63 donors. Using 48 genes located near T2D risk variants, we identified gene coexpression and protein-protein interaction networks that were strongly associated with islet insulin secretion and HbA(1c). We integrated our data to form a rank list of putative T2D genes, of which CHL1, LRFN2, RASGRP1, and PPM1K were validated in INS-1 cells to influence insulin secretion, whereas GPR120 affected apoptosis in islets. Expression variation of the top 20 genes explained 24% of the variance in HbA(1c) with no claim of the direction. The data present a global map of genes associated with islet dysfunction and demonstrate the value of systems genetics for the identification of genes potentially involved in T2D.

First-degree relatives of type 2 diabetic patients have reduced expression of genes involved in fatty acid metabolism in skeletal muscle

CONTEXT

First-degree relatives of patients with type 2 diabetes (FH+) have been shown to have decreased energy expenditure and decreased expression of mitochondrial genes in skeletal muscle. In previous studies, it has been difficult to distinguish whether mitochondrial dysfunction and differential regulation of genes are primary (genetic) or due to reduced physical activity, obesity, or other correlated factors.

OBJECTIVE

The aim of this study was to investigate whether mitochondrial dysfunction is a primary defect or results from an altered metabolic state.

DESIGN

We compared gene expression in skeletal muscle from 24 male subjects with FH and 26 without FH matched for age, glucose tolerance, VO(2peak) (peak oxygen uptake), and body mass index using microarrays. Additionally, type fiber composition, mitochondrial DNA content, and citrate synthase activity were measured. The results were followed up in an additional cohort with measurements of in vivo metabolism.

RESULTS

FH+ vs. FH- subjects showed reduced expression of mitochondrial genes (P = 2.75 × 10(-6)), particularly genes involved in fatty acid metabolism (P = 4.08 × 10(-7)), despite similar mitochondrial DNA content. Strikingly, a 70% reduced expression of the monoamine oxidase A (MAOA) gene was found in FH+ vs. FH- individuals (P = 0.0009). Down-regulation of the genes involved in fat metabolism was associated with decreased in vivo fat oxidation and increased glucose oxidation examined in an additional cohort of elderly men.

CONCLUSIONS

These results suggest that genetically altered fatty acid metabolism predisposes to type 2 diabetes and propose a role for catecholamine-metabolizing enzymes like MAOA in the regulation of energy metabolism.

Differential gene expression in adipose tissue from obese human subjects during weight loss and weight maintenance

BACKGROUND

Differential gene expression in adipose tissue during diet-induced weight loss followed by a weight stability period is poorly characterized. Markers of these processes may provide a deeper understanding of underlying mechanisms.

OBJECTIVE

We aimed to identify differentially expressed genes in human adipose tissue during weight loss and weight maintenance after weight loss.

DESIGN

RNA from subcutaneous abdominal adipose tissue from 9 obese subjects was analyzed by using a complementary DNA microarray at baseline after weight loss on a low-calorie diet and after weight maintenance.

RESULTS

Subjects lost 18.8 ± 1.8% of weight and maintained this loss during weight maintenance (1.1 ± 2.1%; range: -9.3 to 10.6%). Most differentially expressed genes exhibited a reciprocal regulation and returned to baseline after weight loss (2163 genes) and weight maintenance (3175 genes). CETP and ABCG1, both of which participate in the HDL-mediated reverse cholesterol transport (RCT), were among the most upregulated of the 750 genes that were differentially expressed after both processes. Several genes involved in inflammation were downregulated. The use of real-time polymerase chain reaction confirmed or partially confirmed the previously implicated genes TNMD and MMP9 (both downregulated), PNPLA3 (upregulated), and CIDEA and SCD (both reciprocally regulated).

CONCLUSIONS

The beneficial effects of weight loss should be investigated after long-term weight maintenance. The processes of weight loss and weight maintenance should be viewed as biologically distinct. CETP and ABCG1 may be important mediators of these effects through HDL-mediated RCT.

Y chromosome haplogroups and prostate cancer in populations of European and Ashkenazi Jewish ancestry

Genetic variation on the Y chromosome has not been convincingly implicated in prostate cancer risk. To comprehensively analyze the role of inherited Y chromosome variation in prostate cancer risk in individuals of European ancestry, we genotyped 34 binary Y chromosome markers in 3,995 prostate cancer cases and 3,815 control subjects drawn from four studies. In this set, we identified nominally significant association between a rare haplogroup, E1b1b1c, and prostate cancer in stage I (P = 0.012, OR = 0.51; 95% confidence interval 0.30–0.87). Population substructure of E1b1b1c carriers suggested Ashkenazi Jewish ancestry, prompting a replication phase in individuals of both European and Ashkenazi Jewish ancestry. The association was not significant for prostate cancer overall in studies of either Ashkenazi Jewish (1,686 cases and 1,597 control subjects) or European (686 cases and 734 control subjects) ancestry (P_meta = 0.078), but a meta-analysis of stage I and II studies revealed a nominally significant association with prostate cancer risk (P_meta = 0.010, OR = 0.77; 95% confidence interval 0.62–0.94). Comparing haplogroup frequencies between studies, we noted strong similarities between those conducted in the US and France, in which the majority of men carried R1 haplogroups, resembling Northwestern European populations. On the other hand, Finns had a remarkably different haplogroup distribution with a preponderance of N1c and I1 haplogroups. In summary, our results suggest that inherited Y chromosome variation plays a limited role in prostate cancer etiology in European populations but warrant follow-up in additional large and well characterized studies of multiple ethnic backgrounds.

Estimation of electrode location in a rat motor cortex by laminar analysis of electrophysiology and intracortical electrical stimulation

While the development of microelectrode arrays has enabled access to disparate regions of a cortex for neurorehabilitation, neuroprosthetic and basic neuroscience research, accurate interpretation of the signals and manipulation of the cortical neurons depend upon the anatomical placement of the electrode arrays in a layered cortex. Toward this end, this report compares two in vivo methods for identifying the placement of electrodes in a linear array spaced 100 µm apart based on in situ laminar analysis of (1) ketamine-xylazine-induced field potential oscillations in a rat motor cortex and (2) an intracortical electrical stimulation-induced movement threshold. The first method is based on finding the polarity reversal in laminar oscillations which is reported to appear at the transition between layers IV and V in laminar 'high voltage spindles' of the rat cortical column. Analysis of histological images in our dataset indicates that polarity reversal is detected 150.1 ± 104.2 µm below the start of layer V. The second method compares the intracortical microstimulation currents that elicit a physical movement for anodic versus cathodic stimulation. It is based on the hypothesis that neural elements perpendicular to the electrode surface are preferentially excited by anodic stimulation while cathodic stimulation excites those with a direction component parallel to its surface. With this method, we expect to see a change in the stimulation currents that elicits a movement at the beginning of layer V when comparing anodic versus cathodic stimulation as the upper cortical layers contain neuronal structures that are primarily parallel to the cortical surface and lower layers contain structures that are primarily perpendicular. Using this method, there was a 78.7 ± 68 µm offset in the estimate of the depth of the start of layer V. The polarity reversal method estimates the beginning of layer V within ±90 µm with 95% confidence and the intracortical stimulation method estimates it within ±69.3 µm. We propose that these methods can be used to estimate the in situ location of laminar electrodes implanted in the rat motor cortex.

Gene expression in Graves' ophthalmopathy and arm lymphedema: similarities and differences

BACKGROUND

Graves' ophthalmopathy (GO) and lymphedema share some pathogenetic mechanisms, such as edema, inflammation, and adipogenesis. The aim of this study was to examine similarities and differences between chronic GO and chronic lymphedema.

METHODS

Intraorbital adipose tissue was collected from patients with active (n = 10) or chronic GO (n = 10) and thyroid-healthy controls (n = 10). Arm subcutaneous adipose tissue was obtained from patients with chronic arm lymphedema (n = 10), where the unaffected arm served as a control. Gene expression was studied using microarray and real-time polymerase chain reaction.

RESULTS

The following genes were significantly upregulated (p < 0.05) in lymphedema but not in GO and have functions in wound healing, fibrosis, fat metabolism, inflammation, differentiation, development, adhesion, and the cytoskeleton: ATP-binding cassette, sub-family G (WHITE), member 1 (ABCG1), actin, alpha 2, smooth muscle, aorta (ACTA2), secreted frizzled-related protein 2 (SFRP2), tenascin C (TNC), pentraxin-related gene, rapidly induced by IL-1 beta (PTX3), and carboxypeptidase X (M14 family), member 1 (CPMX1). In chronic GO, but not in lymphedema, adipocyte-related immediate early genes known to be overexpressed in patients with active GO were upregulated but at a lower level than previously shown for the active phase. Genes of the Wnt pathway, such as secreted frizzled-related protein 1, 2, and 3, were up- and downregulated in both chronic GO and lymphedema. Parathyroid hormone-like hormone (PTHLH) was downregulated (p = 0.01) and apolipoprotein L domain containing 1 (APOLD1) was upregulated (p = 0.05) in both active and chronic GO.

CONCLUSIONS

There are more differences than similarities between chronic ophthalmopathy and chronic lymphedema, but both conditions exhibit less inflammation and adipogenesis compared to the active phases. In lymphedema, fibrosis dominates. PTHLH, which can inhibit adipogenesis, is downregulated both in active and chronic ophthalmopathy, indicating the possibility of an increased risk of adipogenesis.

Two common genetic variants near nuclear-encoded OXPHOS genes are associated with insulin secretion in vivo

CONTEXT

Mitochondrial ATP production is important in the regulation of glucose-stimulated insulin secretion. Genetic factors may modulate the capacity of the β-cells to secrete insulin and thereby contribute to the risk of type 2 diabetes.

OBJECTIVE

The aim of this study was to identify genetic loci in or adjacent to nuclear-encoded genes of the oxidative phosphorylation (OXPHOS) pathway that are associated with insulin secretion in vivo.

DESIGN AND METHODS

To find polymorphisms associated with glucose-stimulated insulin secretion, data from a genome-wide association study (GWAS) of 1467 non-diabetic individuals, including the Diabetes Genetic Initiative (DGI), was examined. A total of 413 single nucleotide polymorphisms with a minor allele frequency ≥0.05 located in or adjacent to 76 OXPHOS genes were included in the DGI GWAS. A more extensive population-based study of 4323 non-diabetics, the PPP-Botnia, was used as a replication cohort. Insulinogenic index during an oral glucose tolerance test was used as a surrogate marker of glucose-stimulated insulin secretion. Multivariate linear regression analyses were used to test genotype-phenotype associations.

RESULTS

Two common variants were identified in the DGI, where the major C-allele of rs606164, adjacent to NADH dehydrogenase (ubiquinone) 1 subunit C2 (NDUFC2), and the minor G-allele of rs1323070, adjacent to cytochrome c oxidase subunit VIIa polypeptide 2 (COX7A2), showed nominal associations with decreased glucose-stimulated insulin secretion (P=0.0009, respective P=0.003). These associations were replicated in PPP-Botnia (P=0.002 and P=0.05).

CONCLUSION

Our study shows that genetic variation near genes involved in OXPHOS may influence glucose-stimulated insulin secretion in vivo.

Abstract 4965: Genome analysis of a pancreatic cancer susceptibility locus in the NR5A2 gene on chromosome 1q32.1

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

A common variant in TFB1M is associated with reduced insulin secretion and increased future risk of type 2 diabetes

Type 2 diabetes (T2D) evolves when insulin secretion fails. Insulin release from the pancreatic β cell is controlled by mitochondrial metabolism, which translates fluctuations in blood glucose into metabolic coupling signals. We identified a common variant (rs950994) in the human transcription factor B1 mitochondrial (TFB1M) gene associated with reduced insulin secretion, elevated postprandial glucose levels, and future risk of T2D. Because islet TFB1M mRNA levels were lower in carriers of the risk allele and correlated with insulin secretion, we examined mice heterozygous for Tfb1m deficiency. These mice displayed lower expression of TFB1M in islets and impaired mitochondrial function and released less insulin in response to glucose in vivo and in vitro. Reducing TFB1M mRNA and protein in clonal β cells by RNA interference impaired complexes of the mitochondrial oxidative phosphorylation system. Consequently, nutrient-stimulated ATP generation was reduced, leading to perturbed insulin secretion. We conclude that a deficiency in TFB1M and impaired mitochondrial function contribute to the pathogenesis of T2D.

Abstract 164: High-throughput epigenetic analysis of susceptibility loci identified by GWAS in pancreatic cancer

Pancreatic cancer is a highly lethal cancer with few well established risk factors. A genome wide association study (GWAS) of pancreatic cancer (PanScan) is being conducted within the framework of the NCI-sponsored Cohort Consortium and the Pancreatic Cancer Case-Control Consortium (PANC4). Susceptibility loci discovered to date in PanScan are non-coding variants that lie in intronic or intergenic regions, suggesting that the underlying signals may function through regulatory mechanisms that influence gene expression or splicing. Alternatively, these may lie in unannotated transcripts and directly affect their function.

Epigenetic mechanisms, such as DNA methylation, can affect the regulation of gene expression and plays a critical role in the development of many human diseases including cancer. Powerful methods exist to analyze DNA methylation patterns in higher eukaryotes including methylated DNA immunoprecipitation (MeDIP), an affinity based approach to enrich methylated DNA regions from genomic DNA, which can be combined with microarrays to profile genomic DNA methylation patterns. We created a new semi-custom MeDIP-optimized array design based on our Human DNA Methylation 3×720K CpG Island Plus RefSeq Promoter array by tiling additional regions associated with pancreatic cancer that were identified in the PanScan study. To complement DNA methylation data a genome-wide transcriptome analysis was performed with RNA-sequencing (RNA-seq). Here we describe a comprehensive genome wide analysis of 8 pancreatic cancer cell lines to examine methylation patterns of Refseq promoters and annotated CpG islands as well as transcribed sequences.

One of the susceptibility loci from PanScan is in the vicinity of the ABO gene on Chr9q34 where four SNPs (rs505922, rs495828, rs657152 and rs630014) are associated with a significantly increased risk of pancreatic cancer. As a pilot study, DNA methylation patterns and expressed sequences in this locus were investigated in cell lines derived from pancreatic tumors and normal pancreatic tissues. Our genome wide DNA methylation and RNA-seq analysis aims at establishing a comprehensive catalog of epigenetic patterns in pancreatic cell lines that could provide plausibility for the association signal in the ABO gene and other GWAS regions and thereby, initiate the characterization of the molecular phenotype of the susceptibility alleles for pancreatic cancer risk.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 101st Annual Meeting of the American Association for Cancer Research; 2010 Apr 17-21; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2010;70(8 Suppl):Abstract nr 164.

Prioritizing genes for follow-up from genome wide association studies using information on gene expression in tissues relevant for type 2 diabetes mellitus

Background

Genome-wide association studies (GWAS) have emerged as a powerful approach for identifying susceptibility loci associated with polygenetic diseases such as type 2 diabetes mellitus (T2DM). However, it is still a daunting task to prioritize single nucleotide polymorphisms (SNPs) from GWAS for further replication in different population. Several recent studies have shown that genetic variation often affects gene-expression at proximal (cis) as well as distal (trans) genomic locations by different mechanisms such as altering rate of transcription or splicing or transcript stability.

Methods

To prioritize SNPs from GWAS, we combined results from two GWAS related to T2DM, the Diabetes Genetics Initiative (DGI) and the Wellcome Trust Case Control Consortium (WTCCC), with genome-wide expression data from pancreas, adipose tissue, liver and skeletal muscle of individuals with or without T2DM or animal models thereof to identify T2DM susceptibility loci.

Results

We identified 1,170 SNPs associated with T2DM with P < 0.05 in both GWAS and 243 genes that were located in the vicinity of these SNPs. Out of these 243 genes, we identified 115 differentially expressed in publicly available gene expression profiling data. Notably five of them, IGF2BP2, KCNJ11, NOTCH2, TCF7L2 and TSPAN8, have subsequently been shown to be associated with T2DM in different populations. To provide further validation of our approach, we reversed the approach and started with 26 known SNPs associated with T2DM and related traits. We could show that 12 (57%) (HHEX, HNF1B, IGF2BP2, IRS1, KCNJ11, KCNQ1, NOTCH2, PPARG, TCF7L2, THADA, TSPAN8 and WFS1) out of 21 genes located in vicinity of these SNPs were showing aberrant expression in T2DM from the gene expression profiling studies.

Conclusions

Utilizing of gene expression profiling data from different tissues of individuals with or without T2DM or animal models thereof is a powerful tool for prioritizing SNPs from WGAS for further replication studies.

Molecular correlates for maximal oxygen uptake and type 1 fibers

Maximal oxygen uptake (Vo(2max)) and the amount of type 1 fibers are interrelated, but the underlying unifying molecular mechanisms are poorly understood. To explore these mechanisms, we related gene expression profiles in skeletal muscle biopsies of 43 age-matched men from published datasets with Vo(2max) and the amount of type 1 fibers and replicated some of the findings in muscle biopsies from 154 young and elderly individuals using real-time PCR. We identified 66 probe sets (genes or expressed sequence tags) positively and 83 probe sets inversely correlated with Vo(2max) and 171 probe sets positively and 217 probe sets inversely correlated with percentage of type 1 fibers in human skeletal muscle. Genes involved in oxidative phosphorylation (OXPHOS) showed high expression in individuals with high Vo(2max), whereas the opposite was not the case in individuals with low Vo(2max). Instead, genes such as AHNAK and BCL6 were associated with low Vo(2max). Also, expression of the OXPHOS genes NDUFB5 and ATP5C1 increased with exercise training and decreased with aging. In contrast, expression of AHNAK in skeletal muscle decreased with exercise training and increased with aging. Eleven genes (NDUFB4, COX5A, UQCRB, ATP5C1, ATP5G3, ETHE1, FABP3, ISCA1, MYST4, C9orf3, and PKIA) were positively correlated with both Vo(2max) and the percentage of type 1 fibers. Vo(2max) closely reflects expression of OXPHOS genes, particularly that of NDUFB5 and ATP5C1, in skeletal muscle, suggesting good muscle fitness. In contrast, a high expression of AHNAK was associated with a low Vo(2max) and poor muscle fitness.

A phase 1 study of the selective cyclin dependent kinase inhibitor P276–00 in patients with advanced refractory neoplasms

14117

Background: In human cancers, genetic and epigenetic events result in over-expression of cyclins or absence or diminished levels of Cdk inhibitors, providing tumor cells with selective growth advantage. This has prompted the development of pharmacological Cdk inhibitors that could potentially produce anti-tumor effect. P276–00 is a selective Cdk4-D1 and Cdk1-B inhibitor. This study was designed to determine the maximum tolerated dose (MTD), toxicity profile, pharmacokinetics, and antitumour activity of P276–00 given intravenously to patients with advanced refractory solid tumours. Methods: P276–00 was administered in escalating doses to cohorts of eligible patients (pts), starting with a dose of 9 mg/m2 as a 30 minute iv infusion day 1 to 5, and day 8 to 12, q 3 weekly. To date 22 pts have been entered on the study (cohort 1 - 4 pts at 9 mg/m2, cohort 2 - 4 pts at 12.6 mg/m2, cohort 3 - 6 pts at 17.6 mg/m2, cohort 4 - 8 pts at 24.6 mg/m2) with PS 0–2, and mean age of 56 years. Pharmacokinetic profiles were obtained on cycle 1 days 1 and 5. Skin biopsies were obtained immediately prior to starting study treatment and on day 21 of cycle 2 and will be analyzed for Ki67, cleaved caspase 3, phospho-Rb, cyclin D1 and cdk4, and microarray. Results: To date dose limiting toxicity has occurred in one pt. Grade 3 fatigue occurred in 1 pt at 17.6 mg/m2. The most common drug-related adverse events, which were all grade 1 or 2, were fatigue, nausea, hypotension, sweating, and dry mouth. No Grade 3 biochemical toxicities have been reported so far. There have been no responses noted to date. 4 pts have stable disease after 2 cycles. Pharmacokinetic results: The Cmax, t1/2, and AUC0–8 on day 1 were as follows: 9 mg/m2- 315 ng/mL, 6.6 hr, 883 ng.h/mL; 12.6 mg/m2- 402 ng/mL, 5.5 hr, 848 ng.h/mL; 17.6 mg/m2- 589 ng/mL, 5.3 hr, 1289 ng.h/mL; 24.6 mg/m2- 621 ng/mL, 5.6 hr, 1286 ng.h/mL. Conclusions: P276–00 is well tolerated, but grade 3 fatigue has been noted in 1 pt at 17.6 mg/m2 dose level. We have observed confirmed stable disease in one patient. PK results indicate that at 9 mg/m2,12.6 mg/m2, 17.6 mg/m2 and 24.6 mg/m2 we are able to cross the cdk4 enzyme IC50 approximately 10, 13, 19 and 20 times and cross the anti-proliferative IC50 1.1, 1.4, 2.1 and 2.2 times respectively. Accrual continues at the 34.4 mg/m2 dose level.

No significant financial relationships to disclose.

COX-2 and SCD, markers of inflammation and adipogenesis, are related to disease activity in Graves' ophthalmopathy

CONTEXT

Inflammation and adipogenesis are two parallel processes with increased activity in severe Graves' ophthalmopathy.

OBJECTIVE

The aim of this work was to define target genes for therapeutic intervention in adipogenesis and inflammation in Graves' ophthalmopathy.

DESIGN

Orbital tissue was obtained from patients with ophthalmopathy in acute or chronic phase undergoing orbital surgery to study gene expression followed by the study of potential intervention mechanisms in preadipocytes.

SETTING

Clinic of Endocrinology, University Hospital, Malmö, Sweden.

PARTICIPANTS

Patients in acute severe or in chronic phase of ophthalmopathy.

INTERVENTIONS

Lateral orbital decompression in acute phase and restorative surgery in chronic phase. In vitro treatment of preadipocytes with rosiglitazone and diclofenac.

MAIN OUTCOME MEASURE

Gene expression in intraorbital tissue or preadipocytes and differentiation of preadipocytes.

RESULTS

A marker of adipose tissue, stearoyl-coenzyme A desaturase (SCD), and the proinflammatory gene, cyclooxygenase-2 (COX-2), were overexpressed in patients in active phase compared to the chronic phase of ophthalmopathy. In growth-arrested preadipocytes stimulated with rosiglitazone, COX-2 expression increased temporarily within 1 hour and decreased to undetectable levels after 48 hours. In contrast, SCD and peroxisome proliferator-activated receptor-gamma (PPAR-gamma) expression increased continuously from day 2 to day 7 during adipogenesis. Diclofenac, an inhibitor of cyclooxygenases with antagonistic effects on PPAR-gamma, reduced the number of mature adipocytes by approximately 50%.

CONCLUSION

We conclude that inflammation and adipogenesis decrease with a decrease in activity of ophthalmopathy and that the nonsteroidal antiinflammatory drug diclofenac inhibits adipogenesis. This may represent a putative future treatment of endocrine ophthalmopathy.

TXNIP regulates peripheral glucose metabolism in humans

BACKGROUND

Type 2 diabetes mellitus (T2DM) is characterized by defects in insulin secretion and action. Impaired glucose uptake in skeletal muscle is believed to be one of the earliest features in the natural history of T2DM, although underlying mechanisms remain obscure.

METHODS AND FINDINGS

We combined human insulin/glucose clamp physiological studies with genome-wide expression profiling to identify thioredoxin interacting protein (TXNIP) as a gene whose expression is powerfully suppressed by insulin yet stimulated by glucose. In healthy individuals, its expression was inversely correlated to total body measures of glucose uptake. Forced expression of TXNIP in cultured adipocytes significantly reduced glucose uptake, while silencing with RNA interference in adipocytes and in skeletal muscle enhanced glucose uptake, confirming that the gene product is also a regulator of glucose uptake. TXNIP expression is consistently elevated in the muscle of prediabetics and diabetics, although in a panel of 4,450 Scandinavian individuals, we found no evidence for association between common genetic variation in the TXNIP gene and T2DM.

CONCLUSIONS

TXNIP regulates both insulin-dependent and insulin-independent pathways of glucose uptake in human skeletal muscle. Combined with recent studies that have implicated TXNIP in pancreatic beta-cell glucose toxicity, our data suggest that TXNIP might play a key role in defective glucose homeostasis preceding overt T2DM.

A variant in the transcription factor 7-like 2 (TCF7L2) gene is associated with an increased risk of gestational diabetes mellitus

AIMS/HYPOTHESIS

Genetic and epidemiological studies suggest an association between gestational diabetes mellitus and type 2 diabetes. Both are polygenic multifactorial disorders characterised by beta cell dysfunction and insulin resistance. Our aim was to investigate whether common genetic variants that have previously been associated with type 2 diabetes or related phenotypes would also confer risk for gestational diabetes mellitus.

MATERIALS AND METHODS

In 1,881 unrelated pregnant Scandinavian women (649 women with gestational diabetes mellitus, 1,232 non-diabetic control subjects) we genotyped the transcription factor 7-like 2 (TCF7L2 rs7903146), adiponectin (ADIPOQ +276G > T), peroxisome-proliferator activated receptor, gamma 2 (PPARG Pro12Ala), PPARG-coactivator, 1 alpha (PPARGC1A Gly482Ser), forkhead box C2 (FOXC2 -512C > T) and beta3-adrenergic receptor (ADRB3 Trp64Arg) polymorphisms using TaqMan allelic discrimination assay or RFLP.

RESULTS

The CC, CT and TT genotype frequencies of the TCF7L2 rs7903146 variant differed significantly between women with gestational diabetes mellitus and control women (46.3, 43.6 and 10.1% vs 58.5, 35.3 and 6.2%, p = 3.7 x 10(-6), corrected p value [Pc] for multiple testing Pc = 2.2 x 10(-5)). The T-allele was associated with an increased risk of gestational diabetes mellitus (odds ratio 1.49 [95% CI 1.28-1.75], p = 4.9 x 10(-7) [Pc = 2.8 x 10(-6)]). Compared with wild-type CC-genotype carriers, heterozygous (CT-genotype) and homozygous (TT-genotype) carriers had a 1.6-fold (95% CI 1.26-1.93, p = 3.7 x 10(-5) [Pc = 0.0002]) and a 2.1-fold (95% CI 1.41-2.99, p = 0.0001 [Pc = 0.0008]) increased risk of gestational diabetes mellitus, respectively. The other polymorphisms studied were not significantly associated with gestational diabetes mellitus (ADIPOQ +276G > T: 1.17 [1.01-1.36], p = 0.039 [Pc = 0.23]; PPARG Pro12Ala: 1.06 [0.87-1.29], p = 0.53; PPARGC1A Gly482Ser: 0.96 [0.83-1.10], p = 0.54; FOXC2 -512C > T: 1.01 [0.87-1.16], p = 0.94; and ADRB3 Trp64Arg: 1.22 [0.95-1.56], p = 0.12).

CONCLUSIONS/INTERPRETATION

The TCF7L2 rs7903146 variant is associated with an increased risk of gestational diabetes mellitus in Scandinavian women.