E2F transcription factors promote tumorigenicity in pancreatic ductal adenocarcinoma.

BACKGROUND: Pancreatic ductal adenocarcinoma (PDAC) is one of the most lethal cancers with limited treatment options, illustrating an urgent need to identify new drugable targets in PDACs. OBJECTIVE: Using the similarities between tumor development and normal embryonic development, which is accompanied by rapid cell expansion, we aimed to identify and characterize embryonic signaling pathways that were reinitiated during tumor formation and expansion. METHODS AND RESULTS: Here, we report that the transcription factors E2F1 and E2F8 are potential key regulators in PDAC. E2F1 and E2F8 RNA expression is mainly localized in proliferating cells in the developing pancreas and in malignant ductal cells in PDAC. Silencing of E2F1 and E2F8 in PANC-1 pancreatic tumor cells inhibited cell proliferation and impaired cell spreading and migration. Moreover, loss of E2F1 also affected cell viability and apoptosis with E2F expression in PDAC tissues correlating with expression of apoptosis and mitosis pathway genes, suggesting that E2F factors promote cell cycle regulation and tumorigenesis in PDAC cells. CONCLUSION: Our findings illustrate that E2F1 and E2F8 transcription factors are expressed in pancreatic progenitor and PDAC cells, where they contribute to tumor cell expansion by regulation of cell proliferation, viability, and cell migration making these genes attractive therapeutic targets and potential prognostic markers for pancreatic cancer.

DNA methylation risk score for type 2 diabetes is associated with gestational diabetes.

BACKGROUND: Gestational diabetes mellitus (GDM) and type 2 diabetes mellitus (T2DM) share many pathophysiological factors including genetics, but whether epigenetic marks are shared is unknown. We aimed to test whether a DNA methylation risk score (MRS) for T2DM was associated with GDM across ancestry and GDM criteria. METHODS: In two independent pregnancy cohorts, EPIPREG (n = 480) and EPIDG (n = 32), DNA methylation in peripheral blood leukocytes was measured at a gestational age of 28 ± 2. We constructed an MRS in EPIPREG and EPIDG based on CpG hits from a published epigenome-wide association study (EWAS) of T2DM. RESULTS: With mixed models logistic regression of EPIPREG and EPIDG, MRS for T2DM was associated with GDM: odd ratio (OR)[95% CI]: 1.3 [1.1-1.8], P = 0.002 for the unadjusted model, and 1.4 [1.1-1.7], P = 0.00014 for a model adjusted by age, pre-pregnant BMI, family history of diabetes and smoking status. Also, we found 6 CpGs through a meta-analysis (cg14020176, cg22650271, cg14870271, cg27243685, cg06378491, cg25130381) associated with GDM, and some of their methylation quantitative loci (mQTLs) were related to T2DM and GDM. CONCLUSION: For the first time, we show that DNA methylation marks for T2DM are also associated with GDM, suggesting shared epigenetic mechanisms between GDM and T2DM.

Genome-wide mRNA profiling in urinary extracellular vesicles reveals stress gene signature for diabetic kidney disease.

Urinary extracellular vesicles (uEV) are a largely unexplored source of kidney-derived mRNAs with potential to serve as a liquid kidney biopsy. We assessed ∼200 uEV mRNA samples from clinical studies by genome-wide sequencing to discover mechanisms and candidate biomarkers of diabetic kidney disease (DKD) in Type 1 diabetes (T1D) with replication in Type 1 and 2 diabetes. Sequencing reproducibly showed >10,000 mRNAs with similarity to kidney transcriptome. T1D DKD groups showed 13 upregulated genes prevalently expressed in proximal tubules, correlated with hyperglycemia and involved in cellular/oxidative stress homeostasis. We used six of them (GPX3, NOX4, MSRB, MSRA, HRSP12, and CRYAB) to construct a transcriptional "stress score" that reflected long-term decline of kidney function and could even identify normoalbuminuric individuals showing early decline. We thus provide workflow and web resource for studying uEV transcriptomes in clinical urine samples and stress-linked DKD markers as potential early non-invasive biomarkers or drug targets.

Type 2 diabetes candidate genes, including PAX5, cause impaired insulin secretion in human pancreatic islets.

Type 2 diabetes (T2D) is caused by insufficient insulin secretion from pancreatic ß cells. To identify candidate genes contributing to T2D pathophysiology, we studied human pancreatic islets from approximately 300 individuals. We found 395 differentially expressed genes (DEGs) in islets from individuals with T2D, including, to our knowledge, novel (OPRD1, PAX5, TET1) and previously identified (CHL1, GLRA1, IAPP) candidates. A third of the identified expression changes in islets may predispose to diabetes, as expression of these genes associated with HbA1c in individuals not previously diagnosed with T2D. Most DEGs were expressed in human ß cells, based on single-cell RNA-Seq data. Additionally, DEGs displayed alterations in open chromatin and associated with T2D SNPs. Mouse KO strains demonstrated that the identified T2D-associated candidate genes regulate glucose homeostasis and body composition in vivo. Functional validation showed that mimicking T2D-associated changes for OPRD1, PAX5, and SLC2A2 impaired insulin secretion. Impairments in Pax5-overexpressing ß cells were due to severe mitochondrial dysfunction. Finally, we discovered PAX5 as a potential transcriptional regulator of many T2D-associated DEGs in human islets. Overall, we have identified molecular alterations in human pancreatic islets that contribute to ß cell dysfunction in T2D pathophysiology.

The type 1 diabetes gene TYK2 regulates ß-cell development and its responses to interferon-α.

Type 1 diabetes (T1D) is an autoimmune disease that results in the destruction of insulin producing pancreatic ß-cells. One of the genes associated with T1D is TYK2, which encodes a Janus kinase with critical roles in type-Ι interferon (IFN-Ι) mediated intracellular signalling. To study the role of TYK2 in ß-cell development and response to IFNα, we generated TYK2 knockout human iPSCs and directed them into the pancreatic endocrine lineage. Here we show that loss of TYK2 compromises the emergence of endocrine precursors by regulating KRAS expression, while mature stem cell-islets (SC-islets) function is not affected. In the SC-islets, the loss or inhibition of TYK2 prevents IFNα-induced antigen processing and presentation, including MHC Class Ι and Class ΙΙ expression, enhancing their survival against CD8+ T-cell cytotoxicity. These results identify an unsuspected role for TYK2 in ß-cell development and support TYK2 inhibition in adult ß-cells as a potent therapeutic target to halt T1D progression.

Robustness and lethality in multilayer biological molecular networks.

Robustness is a prominent feature of most biological systems. Most previous related studies have been focused on homogeneous molecular networks. Here we propose a comprehensive framework for understanding how the interactions between genes, proteins and metabolites contribute to the determinants of robustness in a heterogeneous biological network. We integrate heterogeneous sources of data to construct a multilayer interaction network composed of a gene regulatory layer, a protein-protein interaction layer, and a metabolic layer. We design a simulated perturbation process to characterize the contribution of each gene to the overall system's robustness, and find that influential genes are enriched in essential and cancer genes. We show that the proposed mechanism predicts a higher vulnerability of the metabolic layer to perturbations applied to genes associated with metabolic diseases. Furthermore, we find that the real network is comparably or more robust than expected in multiple random realizations. Finally, we analytically derive the expected robustness of multilayer biological networks starting from the degree distributions within and between layers. These results provide insights into the non-trivial dynamics occurring in the cell after a genetic perturbation is applied, confirming the importance of including the coupling between different layers of interaction in models of complex biological systems.

Genetic variant effects on gene expression in human pancreatic islets and their implications for T2D.

Most signals detected by genome-wide association studies map to non-coding sequence and their tissue-specific effects influence transcriptional regulation. However, key tissues and cell-types required for functional inference are absent from large-scale resources. Here we explore the relationship between genetic variants influencing predisposition to type 2 diabetes (T2D) and related glycemic traits, and human pancreatic islet transcription using data from 420 donors. We find: (a) 7741 cis-eQTLs in islets with a replication rate across 44 GTEx tissues between 40% and 73%; (b) marked overlap between islet cis-eQTL signals and active regulatory sequences in islets, with reduced eQTL effect size observed in the stretch enhancers most strongly implicated in GWAS signal location; (c) enrichment of islet cis-eQTL signals with T2D risk variants identified in genome-wide association studies; and (d) colocalization between 47 islet cis-eQTLs and variants influencing T2D or glycemic traits, including DGKB and TCF7L2. Our findings illustrate the advantages of performing functional and regulatory studies in disease relevant tissues.

Glucose-dependent insulinotropic peptide and risk of cardiovascular events and mortality: a prospective study.

AIMS/HYPOTHESIS: Evidence that glucose-dependent insulinotropic peptide (GIP) and/or the GIP receptor (GIPR) are involved in cardiovascular biology is emerging. We hypothesised that GIP has untoward effects on cardiovascular biology, in contrast to glucagon-like peptide 1 (GLP-1), and therefore investigated the effects of GIP and GLP-1 concentrations on cardiovascular disease (CVD) and mortality risk. METHODS: GIP concentrations were successfully measured during OGTTs in two independent populations (Malmö Diet Cancer-Cardiovascular Cohort [MDC-CC] and Prevalence, Prediction and Prevention of Diabetes in Botnia [PPP-Botnia]) in a total of 8044 subjects. GLP-1 (n = 3625) was measured in MDC-CC. The incidence of CVD and mortality was assessed via national/regional registers or questionnaires. Further, a two-sample Mendelian randomisation (2SMR) analysis between the GIP pathway and outcomes (coronary artery disease [CAD] and myocardial infarction) was carried out using a GIP-associated genetic variant, rs1800437, as instrumental variable. An additional reverse 2SMR was performed with CAD as exposure variable and GIP as outcome variable, with the instrumental variables constructed from 114 known genetic risk variants for CAD. RESULTS: In meta-analyses, higher fasting levels of GIP were associated with risk of higher total mortality (HR[95% CI] = 1.22 [1.11, 1.35]; p = 4.5 × 10-5) and death from CVD (HR[95% CI] 1.30 [1.11, 1.52]; p = 0.001). In accordance, 2SMR analysis revealed that increasing GIP concentrations were associated with CAD and myocardial infarction, and an additional reverse 2SMR revealed no significant effect of CAD on GIP levels, thus confirming a possible effect solely of GIP on CAD. CONCLUSIONS/INTERPRETATION: In two prospective, community-based studies, elevated levels of GIP were associated with greater risk of all-cause and cardiovascular mortality within 5-9 years of follow-up, whereas GLP-1 levels were not associated with excess risk. Further studies are warranted to determine the cardiovascular effects of GIP per se.

Genetic determinants of circulating GIP and GLP-1 concentrations.

The secretion of insulin and glucagon from the pancreas and the incretin hormones glucagon-like peptide-1 (GLP-1) and glucose-dependent insulinotropic peptide (GIP) from the gastrointestinal tract is essential for glucose homeostasis. Several novel treatment strategies for type 2 diabetes (T2D) mimic GLP-1 actions or inhibit incretin degradation (DPP4 inhibitors), but none is thus far aimed at increasing the secretion of endogenous incretins. In order to identify new potential therapeutic targets for treatment of T2D, we performed a meta-analysis of a GWAS and an exome-wide association study of circulating insulin, glucagon, GIP, and GLP-1 concentrations measured during an oral glucose tolerance test in up to 7,828 individuals. We identified 6 genome-wide significant functional loci associated with plasma incretin concentrations in or near the SLC5A1 (encoding SGLT1), GIPR, ABO, GLP2R, F13A1, and HOXD1 genes and studied the effect of these variants on mRNA expression in pancreatic islet and on metabolic phenotypes. Immunohistochemistry showed expression of GIPR, ABO, and HOXD1 in human enteroendocrine cells and expression of ABO in pancreatic islets, supporting a role in hormone secretion. This study thus provides candidate genes and insight into mechanisms by which secretion and breakdown of GIP and GLP-1 are regulated.

Excess maternal transmission of variants in the THADA gene to offspring with type 2 diabetes.

AIMS/HYPOTHESIS: Genome-wide association studies (GWAS) have identified more than 65 genetic loci associated with risk of type 2 diabetes. However, the contribution of distorted parental transmission of alleles to risk of type 2 diabetes has been mostly unexplored. Our goal was therefore to search for parent-of-origin effects (POE) among type 2 diabetes loci in families. METHODS: Families from the Botnia study (n = 4,211, 1,083 families) were genotyped for 72 single-nucleotide polymorphisms (SNPs) associated with type 2 diabetes and assessed for POE on type 2 diabetes. The family-based Hungarian Transdanubian Biobank (HTB) (n = 1,463, >135 families) was used to replicate SNPs showing POE. Association of type 2 diabetes loci within families was also tested. RESULTS: Three loci showed nominal POE, including the previously reported variants in KCNQ1, for type 2 diabetes in families from Botnia (rs2237895: p POE = 0.037), which can be considered positive controls. The strongest POE was seen for rs7578597 SNP in the THADA gene, showing excess transmission of the maternal risk allele T to diabetic offspring (Botnia: p POE = 0.01; HTB p POE = 0.045). These data are consistent with previous evidence of allelic imbalance for expression in islets, suggesting that the THADA gene can be imprinted in a POE-specific fashion. Five CpG sites, including those flanking rs7578597, showed differential methylation between diabetic and non-diabetic donor islets. CONCLUSIONS/INTERPRETATION: Taken together, the data emphasise the need for genetic studies to consider from which parent an offspring has inherited a susceptibility allele.

MTHFR and MTRR genotype and haplotype analysis and colorectal cancer susceptibility in a case-control study from the Czech Republic.

Polymorphic variants in genes involved in one-carbon metabolism, in particular of dietary folate, may modulate the risk for colorectal cancer through aberrant DNA-methylation and altered nucleotide synthesis and repair. In the present study, we have assessed the association of six polymorphisms and relative haplotypes in the MTHFR gene (rs1801133 and rs1801131) and in the MTRR gene (rs1801394, rs1532268, rs162036, and rs10380) with the risk for colorectal cancer in 666 patients and 1377 controls from the Czech Republic. We found that the 677 C>T polymorphism in the MTHFR gene significantly decreased the risk for colorectal cancer in homozygous carriers of the variant allele (OR, 0.58; 95% CI, 0.39-0.87). Also, we noted a significantly different distribution of genotypes between cases and controls for the 66A>G polymorphism in the MTRR gene. In particular, homozygous carriers of the G-containing allele of this polymorphism were at an increased risk for colorectal cancer (OR, 1.39; 95% CI, 1.04-1.85). Haplotype analysis of the two MTHFR polymorphisms showed a moderate difference in the distribution of the TA haplotype between cases and controls. In comparison to the most common haplotype (CA), the TA haplotype was associated with a decreased risk for colorectal cancer (OR, 0.84; 95% CI, 0.71-0.99). No difference in the distribution between cases and controls was observed for the haplotypes based on the four polymorphisms in the MTRR gene. The present study suggests that the 677TT genotype and the TA haplotype in the MTHFR gene may also have a role in colorectal cancer risk in the Czech population, indicating the importance of genes involved in folate metabolism with respect to cancer risk. For MTRR, additional studies on larger populations are needed to clarify the possible role of variation in this gene in colorectal carcinogenesis.

Variation in CDKN2A at 9p21.3 influences childhood acute lymphoblastic leukemia risk.

Using data from a genome-wide association study of 907 individuals with childhood acute lymphoblastic leukemia (cases) and 2,398 controls and with validation in samples totaling 2,386 cases and 2,419 controls, we have shown that common variation at 9p21.3 (rs3731217, intron 1 of CDKN2A) influences acute lymphoblastic leukemia risk (odds ratio = 0.71, P = 3.01 x 10(-11)), irrespective of cell lineage.

Verification of the susceptibility loci on 7p12.2, 10q21.2, and 14q11.2 in precursor B-cell acute lymphoblastic leukemia of childhood.

Recent genome-wide association data have implicated genetic variation at 7p12.2 (IKZF1), 10q21.2 (ARIDB5), and 14q11.2 (CEBPE) in the etiology of B-cell childhood acute lymphoblastic leukemia (ALL). To verify and further examine the relationship between these variants and ALL risk, we genotyped 1384 cases of precursor B-cell childhood ALL and 1877 controls from Germany and the United Kingdom. The combined data provided statistically significant support for an association between genotype at each of these loci and ALL risk; odds ratios (OR), 1.69 (P = 7.51 x10(-22)), 1.80 (P = 5.90 x 10(-28)), and 1.27 (P = 4.90 x 10(-6)), respectively. Furthermore, the risk of ALL increases with an increasing numbers of variant alleles for the 3 loci (OR(per-allele) = 1.53, 95% confidence interval, 1.44-1.62; P(trend) = 3.49 x 10(-42)), consistent with a polygenic model of disease susceptibility. These data provide unambiguous evidence for the role of these variants in defining ALL risk underscoring approximately 64% of cases.