Ludmila Prokunina-Olsson

A genome-wide association study of type 2 diabetes in Finns detects multiple susceptibility variants.

Identifying the genetic variants that increase the risk of type 2 diabetes (T2D) in humans has been a formidable challenge. Adopting a genome-wide association strategy, we genotyped 1161 Finnish T2D cases and 1174 Finnish normal glucose-tolerant (NGT) controls with >315,000 single-nucleotide polymorphisms (SNPs) and imputed genotypes for an additional >2 million autosomal SNPs. We carried out association analysis with these SNPs to identify genetic variants that predispose to T2D, compared our T2D association results with the results of two similar studies, and genotyped 80 SNPs in an additional 1215 Finnish T2D cases and 1258 Finnish NGT controls. We identify T2D-associated variants in an intergenic region of chromosome 11p12, contribute to the identification of T2D-associated variants near the genes IGF2BP2 and CDKAL1 and the region of CDKN2A and CDKN2B, and confirm that variants near TCF7L2, SLC30A8, HHEX, FTO, PPARG, and KCNJ11 are associated with T2D risk. This brings the number of T2D loci now confidently identified to at least 10.

From noncoding variant to phenotype via SORT1 at the 1p13 cholesterol locus

Recent genome-wide association studies (GWASs) have identified a locus on chromosome 1p13 strongly associated with both plasma low-density lipoprotein cholesterol (LDL-C) and myocardial infarction (MI) in humans. Here we show through a series of studies in human cohorts and human-derived hepatocytes that a common noncoding polymorphism at the 1p13 locus, rs12740374, creates a C/EBP (CCAAT/enhancer binding protein) transcription factor binding site and alters the hepatic expression of the SORT1 gene. With small interfering RNA (siRNA) knockdown and viral overexpression in mouse liver, we demonstrate that Sort1 alters plasma LDL-C and very low-density lipoprotein (VLDL) particle levels by modulating hepatic VLDL secretion. Thus, we provide functional evidence for a novel regulatory pathway for lipoprotein metabolism and suggest that modulation of this pathway may alter risk for MI in humans. We also demonstrate that common noncoding DNA variants identified by GWASs can directly contribute to clinical phenotypes.

A variant upstream of IFNL3 ( IL28B ) creating a new interferon gene IFNL4 is associated with impaired clearance of hepatitis C virus

Chronic infection with hepatitis C virus (HCV) is a common cause of liver cirrhosis and cancer. We performed RNA sequencing in primary human hepatocytes activated with synthetic double-stranded RNA to mimic HCV infection. Upstream of IFNL3 (IL28B) on chromosome 19q13.13, we discovered a new transiently induced region that harbors a dinucleotide variant ss469415590 (TT or ΔG), which is in high linkage disequilibrium with rs12979860, a genetic marker strongly associated with HCV clearance. ss469415590[ΔG] is a frameshift variant that creates a novel gene, designated IFNL4, encoding the interferon-λ4 protein (IFNL4), which is moderately similar to IFNL3. Compared to rs12979860, ss469415590 is more strongly associated with HCV clearance in individuals of African ancestry, although it provides comparable information in Europeans and Asians. Transient overexpression of IFNL4 in a hepatoma cell line induced STAT1 and STAT2 phosphorylation and the expression of interferon-stimulated genes. Our findings provide new insights into the genetic regulation of HCV clearance and its clinical management.

A multistage genome-wide association study in breast cancer identifies two new risk alleles at 1p11.2 and 14q24.1 (RAD51L1).

We conducted a three-stage genome-wide association study (GWAS) of breast cancer in 9,770 cases and 10,799 controls in the Cancer Genetic Markers of Susceptibility (CGEMS) initiative. In stage 1, we genotyped 528,173 SNPs in 1,145 cases of invasive breast cancer and 1,142 controls. In stage 2, we analyzed 24,909 top SNPs in 4,547 cases and 4,434 controls. In stage 3, we investigated 21 loci in 4,078 cases and 5,223 controls. Two new loci achieved genome-wide significance. A pericentromeric SNP on chromosome 1p11.2 (rs11249433; P = 6.74 × 10−10 adjusted genotype test, 2 degrees of freedom) resides in a large linkage disequilibrium block neighboring NOTCH2 and FCGR1B; this signal was stronger for estrogen-receptor–positive tumors. A second SNP on chromosome 14q24.1 (rs999737; P = 1.74 × 10−7) localizes to RAD51L1, a gene in the homologous recombination DNA repair pathway. We also confirmed associations with loci on chromosomes 2q35, 5p12, 5q11.2, 8q24, 10q26 and 16q12.1.

IL-29 Is the Dominant Type III Interferon Produced by Hepatocytes During Acute Hepatitis C Virus Infection

Early, vigorous intrahepatic induction of interferon (IFN)‐stimulated gene (ISG) induction is a feature of hepatitis C virus (HCV) infection, even though HCV inhibits the induction of type I IFNs in vitro. To identify the cytokines and cells that drive ISG induction and mediate antiviral activity during acute HCV infection, type I and III IFN responses were studied in (1) serial liver biopsies and plasma samples obtained from 6 chimpanzees throughout acute HCV infection and (2) primary human hepatocyte (PHH) cultures upon HCV infection. Type I IFNs were minimally induced at the messenger RNA (mRNA) level in the liver and were undetectable at the protein level in plasma during acute HCV infection of chimpanzees. In contrast, type III IFNs, in particular, interleukin (IL)‐29 mRNA and protein, were strongly induced and these levels correlated with ISG expression and viremia. However, there was no association between intrahepatic or peripheral type III IFN levels and the outcome of acute HCV infection. Infection of PHH with HCV recapitulated strong type III and weak type I IFN responses. Supernatants from HCV‐infected PHH cultures mediated antiviral activity upon transfer to HCV‐replicon–containing cells. This effect was significantly reduced by neutralization of type III IFNs and less by neutralization of type I IFNs. Furthermore, IL‐29 production by HCV‐infected PHH occurred independently from type I IFN signaling and was not enhanced by the presence of plasmacytoid dendritic cells. Conclusion: Hepatocyte‐derived type III IFNs contribute to ISG induction and antiviral activity, but are not the principal determinant of the outcome of HCV infection. (HEPATOLOGY 2012;56:2060–2070)

Endogenous intrahepatic IFNs and association with IFN-free HCV treatment outcome

BACKGROUND. Hepatitis C virus (HCV) infects approximately 170 million people worldwide and may lead to cirrhosis and hepatocellular carcinoma in chronically infected individuals. Treatment is rapidly evolving from IFN-α-based therapies to IFN-α-free regimens that consist of directly acting antiviral agents (DAAs), which demonstrate improved efficacy and tolerability in clinical trials. Virologic relapse after DAA therapy is a common cause of treatment failure; however, it is not clear why relapse occurs or whether certain individuals are more prone to recurrent viremia. METHODS. We conducted a clinical trial using the DAA sofosbuvir plus ribavirin (SOF/RBV) and performed detailed mRNA expression analysis in liver and peripheral blood from patients who achieved either a sustained virologic response (SVR) or relapsed. RESULTS. On-treatment viral clearance was accompanied by rapid downregulation of IFN-stimulated genes (ISGs) in liver and blood, regardless of treatment outcome. Analysis of paired pretreatment and end of treatment (EOT) liver biopsies from SVR patients showed that viral clearance was accompanied by decreased expression of type II and III IFNs, but unexpectedly increased expression of the type I IFN IFNA2. mRNA expression of ISGs was higher in EOT liver biopsies of patients who achieved SVR than in patients who later relapsed. CONCLUSION. These results suggest that restoration of type I intrahepatic IFN signaling by EOT may facilitate HCV eradication and prevention of relapse upon withdrawal of SOF/RBV. TRIAL REGISTRATION. ClinicalTrials.gov NCT01441180.

A functional variant at a prostate cancer predisposition locus at 8q24 is associated with PVT1 expression.

Genetic mapping studies have identified multiple cancer susceptibility regions at chromosome 8q24, upstream of the MYC oncogene. MYC has been widely presumed as the regulated target gene, but definitive evidence functionally linking these cancer regions with MYC has been difficult to obtain. Here we examined candidate functional variants of a haplotype block at 8q24 encompassing the two independent risk alleles for prostate and breast cancer, rs620861 and rs13281615. We used the mapping of DNase I hypersensitive sites as a tool to prioritise regions for further functional analysis. This approach identified rs378854, which is in complete linkage disequilibrium (LD) with rs620861, as a novel functional prostate cancer-specific genetic variant. We demonstrate that the risk allele (G) of rs378854 reduces binding of the transcription factor YY1 in vitro. This factor is known to repress global transcription in prostate cancer and is a candidate tumour suppressor. Additional experiments showed that the YY1 binding site is occupied in vivo in prostate cancer, but not breast cancer cells, consistent with the observed cancer-specific effects of this single nucleotide polymorphism (SNP). Using chromatin conformation capture (3C) experiments, we found that the region surrounding rs378854 interacts with the MYC and PVT1 promoters. Moreover, expression of the PVT1 oncogene in normal prostate tissue increased with the presence of the risk allele of rs378854, while expression of MYC was not affected. In conclusion, we identified a new functional prostate cancer risk variant at the 8q24 locus, rs378854 allele G, that reduces binding of the YY1 protein and is associated with increased expression of PVT1 located 0.5 Mb downstream.

Tissue-specific alternative splicing of TCF7L2.

Common variants in the transcription factor 7-like 2 (TCF7L2) gene have been identified as the strongest genetic risk factors for type 2 diabetes (T2D). However, the mechanisms by which these non-coding variants increase risk for T2D are not well-established. We used 13 expression assays to survey mRNA expression of multiple TCF7L2 splicing forms in up to 380 samples from eight types of human tissue (pancreas, pancreatic islets, colon, liver, monocytes, skeletal muscle, subcutaneous adipose tissue and lymphoblastoid cell lines) and observed a tissue-specific pattern of alternative splicing. We tested whether the expression of TCF7L2 splicing forms was associated with single nucleotide polymorphisms (SNPs), rs7903146 and rs12255372, located within introns 3 and 4 of the gene and most strongly associated with T2D. Expression of two splicing forms was lower in pancreatic islets with increasing counts of T2D-associated alleles of the SNPs: a ubiquitous splicing form (P = 0.018 for rs7903146 and P = 0.020 for rs12255372) and a splicing form found in pancreatic islets, pancreas and colon but not in other tissues tested here (P = 0.009 for rs12255372 and P = 0.053 for rs7903146). Expression of this form in glucose-stimulated pancreatic islets correlated with expression of proinsulin (r2 = 0.84–0.90, P < 0.00063). In summary, we identified a tissue-specific pattern of alternative splicing of TCF7L2. After adjustment for multiple tests, no association between expression of TCF7L2 in eight types of human tissue samples and T2D-associated genetic variants remained significant. Alternative splicing of TCF7L2 in pancreatic islets warrants future studies. GenBank Accession Numbers: FJ010164–FJ010174.

Association of PDCD1 genetic variation with risk and clinical manifestations of systemic lupus erythematosus in a multiethnic cohort

We evaluated the roles of five single-nucleotide polymorphisms (SNPs) within PDCD1, and haplotypes defined by these SNPs, for the development of systemic lupus erythematosus (SLE) and specific sub-phenotypes (nephritis, antiphospholipid antibody positive, arthritis and double-stranded DNA positive) within a multiethnic US cohort of 1036 patients. Family based analyses were performed using 844 simplex families from four ethnic groups (Caucasian, Asian, Hispanic and African American). Subjects were genotyped for five ‘tag’ SNPs (selected from 15) to provide complete genetic information in all main ethnic groups. We employed transmission disequilibrium testing to assess risk for SLE by allele or haplotype, and multiple logistic regression analysis of SLE cases to examine associations with specific sub-phenotypes. In family based analyses, a haplotype containing the PD1.3A allele was significantly associated with SLE susceptibility among Caucasian families (P=0.01). Among Hispanic families, two novel SNPs were associated with SLE risk (P=0.005 and 0.01). In multivariate logistic regression analyses, five haplotypes were associated with specific sub-phenotypes among the different ethnic groups. These results suggest that PDCD1 genetic variation influences the risk and expression of SLE and that these associations vary according to ethnic background.

A common variant associated with dyslexia reduces expression of the KIAA0319 gene

Numerous genetic association studies have implicated the KIAA0319 gene on human chromosome 6p22 in dyslexia susceptibility. The causative variant(s) remains unknown but may modulate gene expression, given that (1) a dyslexia-associated haplotype has been implicated in the reduced expression of KIAA0319, and (2) the strongest association has been found for the region spanning exon 1 of KIAA0319. Here, we test the hypothesis that variant(s) responsible for reduced KIAA0319 expression resides on the risk haplotype close to the genes transcription start site. We identified seven single-nucleotide polymorphisms on the risk haplotype immediately upstream of KIAA0319 and determined that three of these are strongly associated with multiple reading-related traits. Using luciferase-expressing constructs containing the KIAA0319 upstream region, we characterized the minimal promoter and additional putative transcriptional regulator regions. This revealed that the minor allele of rs9461045, which shows the strongest association with dyslexia in our sample (max p-value = 0.0001), confers reduced luciferase expression in both neuronal and non-neuronal cell lines. Additionally, we found that the presence of this rs9461045 dyslexia-associated allele creates a nuclear protein-binding site, likely for the transcriptional silencer OCT-1. Knocking down OCT-1 expression in the neuronal cell line SHSY5Y using an siRNA restores KIAA0319 expression from the risk haplotype to nearly that seen from the non-risk haplotype. Our study thus pinpoints a common variant as altering the function of a dyslexia candidate gene and provides an illustrative example of the strategic approach needed to dissect the molecular basis of complex genetic traits.