3D genomic analysis reveals novel enhancer-hijacking caused by complex structural alterations that drive oncogene overexpression

Enhancer hijacking, caused by structural alterations on chromosomes as well as extrachromosomal DNA (ecDNA), is a common cancer driver event. The complexity and ubiquity of structural alterations in cancer genomes make it difficult to identify enhancer hijacking using genome sequencing alone. Here we describe a 3D genomics-based analysis called HAPI (Highly Active Promoter Interactions) to characterize enhancer hijacking caused by structural alterations. HAPI analysis of HiChIP data from 34 cancer cell lines identified novel enhancer hijacking events that involve chromosomal rearrangements and activate both known and potentially novel oncogenes such as MYC, CCND1, ETV1, CRKL, and ID4, which we validated using CRISPRi assays and RNA-seq analysis. Furthermore, we found that ecDNAs often contain multiple oncogenes from different chromosomes, which causes nested enhancer hijacking among them. We found that ecDNAs containing MYC often harbor additional oncogenes from other chromosomes such as CDX2, ERBB2, or CD44 that co-opt MYC's enhancers for their overexpression, which we validated using dual-color DNA FISH and CRISPRi assays. These enhancer hijacking events involving multiple oncogenes on ecDNAs have important implications for therapeutic strategies that either target the co-opting oncogenes or the hijacked enhancers. Our publicly available HAPI analysis tool provides a robust strategy to detect enhancer hijacking and reveals novel insights into oncogene activation caused by chromosomal and extrachromosomal structural alterations.

Soluble factors in malignant ascites promote the metastatic adhesion of gastric adenocarcinoma cells

BACKGROUND

Adenocarcinoma of the proximal stomach is the fastest rising malignancy in North America. It is commonly associated with peritoneal accumulation of malignant ascites (MA), a fluid containing cancer and inflammatory cells and soluble proteins. Peritoneal metastasis (PM) is the most common site of gastric cancer (GC) progression after curative-intent surgery and is the leading cause of death among GC patients.

METHODS/RESULTS

Using a panel of gastric adenocarcinoma cell lines (human: MKN 45, SNU-5; murine: NCC-S1M), we demonstrate that prior incubation of GC cells with MA results in a significant (> 1.7-fold) increase in the number of cells capable of adhering to human peritoneal mesothelial cells (HPMC) (p < 0.05). We then corroborate these findings using an ex vivo PM model and show that MA also significantly enhances the ability of GC cells to adhere to strips of human peritoneum (p < 0.05). Using a multiplex ELISA, we identify MIF and VEGF as consistently elevated across MA samples from GC patients (p < 0.05). We demonstrate that agents that block the effects of MIF or VEGF abrogate the ability of MA to stimulate the adhesion of GC cells to adhere to human peritoneum and promote both ex vivo and in vivo metastases.

CONCLUSION

Agents targeting MIF or VEGF may be relevant to the treatment or prevention of PM in GC patients.

To bind or not to bind: Cistromic reprogramming in prostate cancer

The term “cistrome” refers to the genome-wide location of regulatory elements associated with transcription factor binding-sites. The cistrome of key regulatory factors in prostate cancer etiology are substantially reprogrammed and altered during prostatic transformation and disease progression. For instance, the cistrome of the androgen receptor (AR), a ligand-inducible transcription factor central in normal prostate epithelium biology, is directly impacted and substantially reprogrammed during malignant transformation. Accumulating evidence demonstrates that additional transcription factors that are frequently mutated, or aberrantly expressed in prostate cancer, such as the pioneer transcription factors Forkhead Box A1 (FOXA1), the homeobox protein HOXB13, and the GATA binding protein 2 (GATA2), and the ETS-related gene (ERG), and the MYC proto-oncogene, contribute to the reprogramming of the AR cistrome. In addition, recent findings have highlighted key roles for the SWI/SNF complex and the chromatin-modifying helicase CHD1 in remodeling the epigenome and altering the AR cistrome during disease progression. In this review, we will cover the role of cistromic reprogramming in prostate cancer initiation and progression. Specifically, we will discuss the impact of key prostate cancer regulators, as well as the role of epigenetic and chromatin regulators in relation to the AR cistrome and the transformation of normal prostate epithelium. Given the importance of chromatin-transcription factor dynamics in normal cellular differentiation and cancer, an in-depth assessment of the factors involved in producing these altered cistromes is of great relevance and provides insight into new therapeutic strategies for prostate cancer.

Failure of human rhombic lip differentiation underlies medulloblastoma formation

Medulloblastoma (MB) comprises a group of heterogeneous paediatric embryonal neoplasms of the hindbrain with strong links to early development of the hindbrain^1-4. Mutations that activate Sonic hedgehog signalling lead to Sonic hedgehog MB in the upper rhombic lip (RL) granule cell lineage^5-8. By contrast, mutations that activate WNT signalling lead to WNT MB in the lower RL^9,10. However, little is known about the more commonly occurring group 4 (G4) MB, which is thought to arise in the unipolar brush cell lineage^3,4. Here we demonstrate that somatic mutations that cause G4 MB converge on the core binding factor alpha (CBFA) complex and mutually exclusive alterations that affect CBFA2T2, CBFA2T3, PRDM6, UTX and OTX2. CBFA2T2 is expressed early in the progenitor cells of the cerebellar RL subventricular zone in Homo sapiens, and G4 MB transcriptionally resembles these progenitors but are stalled in developmental time. Knockdown of OTX2 in model systems relieves this differentiation blockade, which allows MB cells to spontaneously proceed along normal developmental differentiation trajectories. The specific nature of the split human RL, which is destined to generate most of the neurons in the human brain, and its high level of susceptible EOMES⁺KI67⁺ unipolar brush cell progenitor cells probably predisposes our species to the development of G4 MB.

A predominant enhancer co-amplified with the SOX2 oncogene is necessary and sufficient for its expression in squamous cancer

Amplification and overexpression of the SOX2 oncogene represent a hallmark of squamous cancers originating from diverse tissue types. Here, we find that squamous cancers selectively amplify a 3' noncoding region together with SOX2, which harbors squamous cancer-specific chromatin accessible regions. We identify a single enhancer e1 that predominantly drives SOX2 expression. Repression of e1 in SOX2-high cells causes collapse of the surrounding enhancers, remarkable reduction in SOX2 expression, and a global transcriptional change reminiscent of SOX2 knockout. The e1 enhancer is driven by a combination of transcription factors including SOX2 itself and the AP-1 complex, which facilitates recruitment of the co-activator BRD4. CRISPR-mediated activation of e1 in SOX2-low cells is sufficient to rebuild the e1-SOX2 loop and activate SOX2 expression. Our study shows that squamous cancers selectively amplify a predominant enhancer to drive SOX2 overexpression, uncovering functional links among enhancer activation, chromatin looping, and lineage-specific copy number amplifications of oncogenes.

Chromatin Looping Shapes KLF5-Dependent Transcriptional Programs in Human Epithelial Cancers

Activation of transcription factors is a key driver event in cancer. We and others have recently reported that the Krüppel-like transcription factor KLF5 is activated in multiple epithelial cancer types including squamous cancer and gastrointestinal adenocarcinoma, yet the functional consequences and the underlying mechanisms of this activation remain largely unknown. Here we demonstrate that activation of KLF5 results in strongly selective KLF5 dependency for these cancer types. KLF5 bound lineage-specific regulatory elements and activated gene expression programs essential to cancer cells. HiChIP analysis revealed that multiple distal KLF5 binding events cluster and synergize to activate individual target genes. Immunoprecipitation-mass spectrometry assays showed that KLF5 interacts with other transcription factors such as TP63 and YAP1, as well as the CBP/EP300 acetyltransferase complex. Furthermore, KLF5 guided the CBP/EP300 complex to increase acetylation of H3K27, which in turn enhanced recruitment of the bromodomain protein BRD4 to chromatin. The 3D chromatin architecture aggregated KLF5-dependent BRD4 binding to activate polymerase II elongation at KLF5 target genes, which conferred a transcriptional vulnerability to proteolysis-targeting chimera-induced degradation of BRD4. Our study demonstrates that KLF5 plays an essential role in multiple epithelial cancers by activating cancer-related genes through 3D chromatin loops, providing an evidence-based rationale for targeting the KLF5 pathway. SIGNIFICANCE: An integrative 3D genomics methodology delineates mechanisms underlying the function of KLF5 in multiple epithelial cancers and suggests potential strategies to target cancers with aberrantly activated KLF5.

Author Correction: Human somatic cell mutagenesis creates genetically tractable sarcomas

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Recurrent noncoding U1 snRNA mutations drive cryptic splicing in SHH medulloblastoma

Recurrent somatic single nucleotide variants (SNVs) in cancer are largely confined to protein coding genes, and are rare in most pediatric cancers^1–3. We report highly recurrent hotspot mutations of U1 spliceosomal small nuclear RNAs (snRNAs) in ~50% of Sonic Hedgehog medulloblastomas (Shh-MB), which were not present across other medulloblastoma subgroups. This U1-snRNA hotspot mutation (r.3a>g), was identified in <0.1% of 2,442 cancers across 36 other tumor types. Largely absent from infant Shh-MB, the mutation occurs in 97% of adults (Shhδ), and 25% of adolescents (Shhα). The U1-snRNA mutation occurs in the 5′ splice site binding region, and snRNA mutant tumors have significantly disrupted RNA splicing with an excess of 5′ cryptic splicing events. Mutant U1-snRNA mediated alternative splicing inactivates tumor suppressor genes (PTCH1), and activates oncogenes (GLI2, CCND2), represents a novel target for therapy, and constitutes a highly recurrent and tissue-specific mutation of a non-protein coding gene in cancer.

Pre-neoplastic epigenetic disruption of transcriptional enhancers in chronic inflammation

Chronic periodontitis (CP) is a chronic inflammatory disease independently associated with higher incidence of oral cavity squamous cell carcinoma (OSCC). However, the molecular mechanism responsible for this increased incidence is unknown. Here we profiled the DNA methylome of CP patients and healthy controls and compared to a large set of OSCC samples from TCGA. We observed a significant overlap between the altered DNA methylation patterns in CP and in OSCC, suggesting an emergence of a pre-neoplastic epigenome in CP. Remarkably, the hypermethylated CpGs in CP were significantly enriched for enhancer elements. This aberrant enhancer methylation is functional and able to disrupt enhancer activity by preventing the binding of chromatin looping factors. This study provides new insights on the molecular mechanisms linking chronic inflammation and tumor predisposition, highlighting the role of epigenetic disruption of transcriptional enhancers.

Noncoding somatic and inherited single-nucleotide variants converge to promote ESR1 expression in breast cancer

Sustained expression of the oestrogen receptor alpha (ESR1) drives two-thirds of breast cancer and defines the ESR1-positive subtype. ESR1 engages enhancers upon oestrogen stimulation to establish an oncogenic expression program¹. Somatic copy number alterations involving the ESR1 gene occur in approximately 1% of ESR1-positive breast cancers^2–5, implying that other mechanisms underlie the persistent expression of ESR1. We report the significant enrichment of somatic mutations within the set of regulatory elements (SRE) regulating ESR1 in 7% of ESR1-positive breast cancers. These mutations regulate ESR1 expression by modulating transcription factor binding to the DNA. The SRE includes a recurrently mutated enhancer whose activity is also affected by a functional inherited single nucleotide variant (SNV) rs9383590 that accounts for several breast cancer risk-loci. Our work highlights the importance of considering the combinatorial activity of regulatory elements as a single unit to delineate the impact of noncoding genetic alterations on single genes in cancer.

Abstract B04: In vivo BioID identifies novel Myc interacting partners

Myc oncoprotein is a major driver of cancer initiation and progression, and thus targeting its activity would mark a key therapeutic advance. In a genetic preclinical mouse model, systemic Myc inhibition using the dominant-negative Myc mutant, termed Omomyc, showed that Ras-driven lung cancer could be eradicated without any harmful long-term effects to the animal. However, developing an anti-cancer agent that directly binds and inhibits Myc has not been possible, to date. Therefore, new strategies are required to inhibit Myc in cancer. Understanding the Myc interactome may unravel novel approaches to target Myc in cancer. The BioID proximity-based biotin labeling technique was recently developed for the characterization of protein-protein interaction networks. In BioID, the protein of interest is expressed as a fusion partner biotin ligase (BirA*), which activates biotin. The active biotin reacts with lysine residues on nearby polypeptides. Following a stringent cell lysis and streptavidin-sepharose pulldown, biotinylated proteins can be identified using MS. To date, this method has been applied to a number of different polypeptides expressed in cultured cells. Here we report the adaptation of BioID to the identification of protein-protein interactions surrounding the Myc oncoprotein in human cells grown both under standard culture conditions and in mice as tumor xenografts. Notably, in vivo BioID yielded &gt;100 high confidence Myc interacting proteins, including &gt;30 known binding partners such as MAX (Myc-associated factor X), TRRAP (transformation/transcription domain-associated protein), the enhancer of polycomb homologs 1 and 2 (EPC1, EPC2), lysine acetyltransferase 5 (KAT5). Putative novel Myc interactors include components of the STAGA/KAT5 and SWI/SNF chromatin remodelling complexes (see Penn lab abstract Tu et al), DNA repair and replication factors, general transcription and elongation factors, and transcriptional co-regulators such as the DNA helicase chromodomain 8 (CHD8). Providing additional confidence in these findings, ENCODE ChIP-seq datasets highlight significant coincident binding throughout the genome for the Myc interactors identified here, and we validate the previously unreported CHD8 (an ATP-dependent helicase)-Myc interaction using both a yeast two hybrid analysis and the proximity-based ligation assay (PLA). Additionally, we also validate Myc-BRD4 and Myc-TRIM24 interaction by PLA. In sum, here we identify bona fide interacting partners of Myc in vivo by use of BioID. Our study shows for the first time Myc interactome in vivo, understanding these interactors will shed more light on Myc oncogenesis, which can be used to therapeutically target Myc in cancer.

Citation Format: Dharmendra Dingar, Manpreet Kalkat, Pak-Kei Chan, Swneke D. Bailey, Tharan Srikumar, William B. Tu, Etienne Coyaud, Romina Ponzielli, Max Kolyar, Igor Jurisica, Annie Huang, Mathieu Lupien, Brian Raught, Linda Z. Penn. In vivo BioID identifies novel Myc interacting partners. [abstract]. In: Proceedings of the AACR Special Conference on Myc: From Biology to Therapy; Jan 7-10, 2015; La Jolla, CA. Philadelphia (PA): AACR; Mol Cancer Res 2015;13(10 Suppl):Abstract nr B04.

ABC: a tool to identify SNVs causing allele-specific transcription factor binding from ChIP-Seq experiments

MOTIVATION

Detection of allelic imbalances in ChIP-Seq reads is a powerful approach to identify functional non-coding single nucleotide variants (SNVs), either polymorphisms or mutations, which modulate the affinity of transcription factors for chromatin. We present ABC, a computational tool that identifies allele-specific binding of transcription factors from aligned ChIP-Seq reads at heterozygous SNVs. ABC controls for potential false positives resulting from biases introduced by the use of short sequencing reads in ChIP-Seq and can efficiently process a large number of heterozygous SNVs.

RESULTS

ABC successfully identifies previously characterized functional SNVs, such as the rs4784227 breast cancer risk associated SNP that modulates the affinity of FOXA1 for the chromatin.

AVAILABILITY AND IMPLEMENTATION

The code is open-source under an Artistic-2.0 license and versioned on GitHub (https://github.com/mlupien/ABC/). ABC is written in PERL and can be run on any platform with both PERL (≥5.18.1) and R (≥3.1.1) installed. The script requires the PERL Statistics::R module.

CONTACT

mlupien@uhnres.utoronto.ca

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

ZNF143 provides sequence specificity to secure chromatin interactions at gene promoters

Chromatin interactions connect distal regulatory elements to target gene promoters guiding stimulus- and lineage-specific transcription. Few factors securing chromatin interactions have so far been identified. Here, by integrating chromatin interaction maps with the large collection of transcription factor-binding profiles provided by the ENCODE project, we demonstrate that the zinc-finger protein ZNF143 preferentially occupies anchors of chromatin interactions connecting promoters with distal regulatory elements. It binds directly to promoters and associates with lineage-specific chromatin interactions and gene expression. Silencing ZNF143 or modulating its DNA-binding affinity using single-nucleotide polymorphisms (SNPs) as a surrogate of site-directed mutagenesis reveals the sequence dependency of chromatin interactions at gene promoters. We also find that chromatin interactions alone do not regulate gene expression. Together, our results identify ZNF143 as a novel chromatin-looping factor that contributes to the architectural foundation of the genome by providing sequence specificity at promoters connected with distal regulatory elements.

Chromatin interactions can connect distal regulatory elements to promoters via protein factors, but few such factors have been identified. Here, the authors show that zinc-finger protein ZNF143 is a sequence-specific chromatin-looping factor that connects promoters with distal regulatory elements.

BioID identifies novel c-MYC interacting partners in cultured cells and xenograft tumors

The BioID proximity-based biotin labeling technique was recently developed for the characterization of protein-protein interaction networks [1]. To date, this method has been applied to a number of different polypeptides expressed in cultured cells. Here we report the adaptation of BioID to the identification of protein-protein interactions surrounding the c-MYC oncoprotein in human cells grown both under standard culture conditions and in mice as tumor xenografts. Notably, in vivo BioID yielded >100 high confidence MYC interacting proteins, including >30 known binding partners. Putative novel MYC interactors include components of the STAGA/KAT5 and SWI/SNF chromatin remodeling complexes, DNA repair and replication factors, general transcription and elongation factors, and transcriptional co-regulators such as the DNA helicase protein chromodomain 8 (CHD8). Providing additional confidence in these findings, ENCODE ChIP-seq datasets highlight significant coincident binding throughout the genome for the MYC interactors identified here, and we validate the previously unreported MYC-CHD8 interaction using both a yeast two hybrid analysis and the proximity-based ligation assay. In sum, we demonstrate that BioID can be utilized to identify bona fide interacting partners for a chromatin-associated protein in vivo. This technique will allow for a much improved understanding of protein-protein interactions in a previously inaccessible biological setting.

BIOLOGICAL SIGNIFICANCE

The c-MYC (MYC) oncogene is a transcription factor that plays important roles in cancer initiation and progression. MYC expression is deregulated in more than 50% of human cancers, but the role of this protein in normal cell biology and tumor progression is still not well understood, in part because identifying MYC-interacting proteins has been technically challenging: MYC-containing chromatin-associated complexes are difficult to isolate using traditional affinity purification methods, and the MYC protein is exceptionally labile, with a half-life of only ~30 min. Developing a new strategy to gain insight into MYC-containing protein complexes would thus mark a key advance in cancer research. The recently described BioID proximity-based labeling technique represents a promising new complementary approach for the characterization of protein-protein interactions (PPIs) in cultured cells. Here we report that BioID can also be used to characterize protein-protein interactions for a chromatin-associated protein in tumor xenografts, and present a comprehensive, high confidence in vivo MYC interactome. This article is part of a Special Issue entitled: Protein dynamics in health and disease. Guest Editors: Pierre Thibault and Anne-Claude Gingras.

Enhancer alterations in cancer: a source for a cell identity crisis

Enhancers are selectively utilized to orchestrate gene expression programs that first govern pluripotency and then proceed to highly specialized programs required for the process of cellular differentiation. Whereas gene-proximal promoters are typically active across numerous cell types, distal enhancer activation is cell-type-specific and central to cell fate determination, thereby accounting for cell identity. Recent studies have highlighted the diversity of enhancer usage, cataloguing millions of such elements in the human genome. The disruption of enhancer activity, through genetic or epigenetic alterations, can impact cell-type-specific functions, resulting in a wide range of pathologies. In cancer, these alterations can promote a 'cell identity crisis', in which enhancers associated with oncogenes and multipotentiality are activated, while those promoting cell fate commitment are inactivated. Overall, these alterations favor an undifferentiated cellular phenotype. Here, we review the current knowledge regarding the role of enhancers in normal cell function, and discuss how genetic and epigenetic changes in enhancer elements potentiate oncogenesis. In addition, we discuss how understanding the mechanisms regulating enhancer activity can inform therapeutic opportunities in cancer cells and highlight key challenges that remain in understanding enhancer biology as it relates to oncology.

Variation at the DPP4 locus influences apolipoprotein B levels in South Asians and exhibits heterogeneity in Europeans related to BMI

AIMS/HYPOTHESIS

Dyslipidaemia, a common feature of type 2 diabetes, is characterised by an increase in atherogenic particles, quantifiable through apolipoprotein B (ApoB) levels. Genetic studies of lipid levels have focused on Europeans; a study in South Asians could identify novel genes.

METHODS

We tested 31,739 single nucleotide polymorphisms (SNPs) from ∼ 2,000 genes in 2,573 South Asians from the epidemiological arm of the Diabetes Reduction Assessment with Ramipril and Rosiglitazone Medication (DREAM) study (EpiDREAM) for association with ApoB and we tested two novel associations for replication in 1,181 South Asians from the INTERHEART case-control study.

RESULTS

The SNP, rs4664443, within DPP4 was associated with ApoB (p = 7.98 × 10(-5)) in EpiDREAM. The observed association was replicated in the INTERHEART South Asians (one-sided p = 9.65 × 10(-3); combined two-sided p = 4.68 × 10(-6)). The rs4664443 SNP was not associated with ApoB among five other EpiDREAM ethnicities. However, because South Asians had a significantly lower mean BMI compared with other EpiDREAM ethnicities, we tested for and found an interaction between rs4664443 and BMI for ApoB among the Europeans, the largest subgroup in EpiDREAM (p = 4.14 × 10(-3) for interaction), observing an association with ApoB in Europeans with a BMI <25 kg/m(2) (p = 2.35 × 10(-3)), but not with a BMI ≥ 25 kg/m(2) (p = 0.21). The association between rs4664443 and ApoB among all EpiDREAM individuals with BMI <25 kg/m(2) was significant (n = 2,972; p = 1.44 × 10(-5)) compared with those with a BMI ≥ 25 kg/m(2) (n = 11,559; p = 0.81), and there was evidence of association among all genotyped individuals with a BMI <25 kg/m(2), including the INTERHEART South Asians (n = 3,601; p = 9.52 × 10(-7)).

CONCLUSION/INTERPRETATION

Variation at the DPP4 locus is associated with ApoB in South Asians and displays heterogeneity related to BMI in other ethnicities.

Laying a solid foundation for Manhattan--'setting the functional basis for the post-GWAS era'

Genome-wide association studies (GWAS) have identified more than 8900 genetic variants, mainly single-nucleotide polymorphisms (SNPs), associated with hundreds of human traits and diseases, which define risk-associated loci. Variants that map to coding regions can affect protein sequence, translation rate, and alternative splicing, all of which influence protein function. However, the vast majority of sequence variants map to non-coding intergenic and intronic regions, and it has been much more challenging to assess the functional nature of these variants. Recent work annotating the non-coding regions of the genome has contributed to post-GWAS studies by facilitating the identification of the functional targets of risk-associated loci. Many non-coding genetic variants within risk-associated loci alter gene expression by modulating the activity of cis-regulatory elements. We review here these recent findings, discuss their implication for the post-GWAS era, and relate their importance to the interpretation of disease-associated mutations identified through whole-genome sequencing.

Genetic information and the prediction of incident type 2 diabetes in a high-risk multiethnic population: the EpiDREAM genetic study

OBJECTIVE

To determine if 16 single nucleotide polymorphisms (SNPs) associated with type 2 diabetes (T2DM) in Europeans are also associated with T2DM in South Asians and Latinos and if they can add to the prediction of incident T2DM in a high-risk population.

RESEARCH DESIGN AND METHODS

In the EpiDREAM prospective cohort study, physical measures, questionnaires, and blood samples were collected from 25,063 individuals at risk for dysglycemia. Sixteen SNPs that have been robustly associated with T2DM in Europeans were genotyped. Among 15,466 European, South Asian, and Latino subjects, we examined the association of these 16 SNPs alone and combined in a gene score with incident cases of T2DM (n = 1,016) that developed during 3.3 years of follow-up.

RESULTS

Nine of the 16 SNPs were significantly associated with T2DM, and their direction of effect was consistent across the three ethnic groups. The gene score was significantly higher among subjects who developed incident T2DM (cases vs. noncases: 16.47 [2.50] vs. 15.99 [2.56]; P = 0.00001). The gene score remained an independent predictor of incident T2DM, with an odds ratio of 1.08 (95% CI 1.05-1.11) per additional risk allele after adjustment for T2DM risk factors. The gene score in those with no family history of T2DM was 16.02, whereas it was 16.19 in those with one parent with T2DM and it was 16.32 in those with two parents with T2DM (P trend = 0.0004). The C statistic of T2DM risk factors was 0.708 (0.691-0.725) and increased only marginally to 0.714 (0.698-0.731) with the addition of the gene score (P for C statistic change = 0.0052).

CONCLUSIONS

T2DM genetic associations are generally consistent across ethnic groups, and a gene score only adds marginal information to clinical factors for T2DM prediction.

Breast cancer risk-associated SNPs modulate the affinity of chromatin for FOXA1 and alter gene expression

Genome-wide association studies (GWAS) have identified thousands of SNPs that are associated with human traits and diseases. But, because the vast majority of these SNPs are located in non-coding regions of the genome, the mechanisms by which they promote disease risk have remained elusive. Employing a new methodology that combines cistromics, epigenomics and genotype imputation, we annotate the non-coding regions of the genome in breast cancer cells and systematically identify the functional nature of SNPs associated with breast cancer risk. Our results show that breast cancer risk-associated SNPs are enriched in the cistromes of FOXA1 and ESR1 and the epigenome of histone H3 lysine 4 monomethylation (H3K4me1) in a cancer- and cell type-specific manner. Furthermore, the majority of the risk-associated SNPs modulate the affinity of chromatin for FOXA1 at distal regulatory elements, thereby resulting in allele-specific gene expression, which is exemplified by the effect of the rs4784227 SNP on the TOX3 gene within the 16q12.1 risk locus.

The effect of chromosome 9p21 variants on cardiovascular disease may be modified by dietary intake: evidence from a case/control and a prospective study

BACKGROUND

One of the most robust genetic associations for cardiovascular disease (CVD) is the Chromosome 9p21 region. However, the interaction of this locus with environmental factors has not been extensively explored. We investigated the association of 9p21 with myocardial infarction (MI) in individuals of different ethnicities, and tested for an interaction with environmental factors.

METHODS AND FINDINGS

We genotyped four 9p21 SNPs in 8,114 individuals from the global INTERHEART study. All four variants were associated with MI, with odds ratios (ORs) of 1.18 to 1.20 (1.85×10(-8)≤p≤5.21×10(-7)). A significant interaction (p = 4.0×10(-4)) was observed between rs2383206 and a factor-analysis-derived "prudent" diet pattern score, for which a major component was raw vegetables. An effect of 9p21 on MI was observed in the group with a low prudent diet score (OR = 1.32, p = 6.82×10(-7)), but the effect was diminished in a step-wise fashion in the medium (OR = 1.17, p = 4.9×10(-3)) and high prudent diet scoring groups (OR = 1.02, p = 0.68) (p = 0.014 for difference). We also analyzed data from 19,129 individuals (including 1,014 incident cases of CVD) from the prospective FINRISK study, which used a closely related dietary variable. In this analysis, the 9p21 risk allele demonstrated a larger effect on CVD risk in the groups with diets low or average for fresh vegetables, fruits, and berries (hazard ratio [HR] = 1.22, p = 3.0×10(-4), and HR = 1.35, p = 4.1×10(-3), respectively) compared to the group with high consumption of these foods (HR = 0.96, p = 0.73) (p = 0.0011 for difference). The combination of the least prudent diet and two copies of the risk allele was associated with a 2-fold increase in risk for MI (OR = 1.98, p = 2.11×10(-9)) in the INTERHEART study and a 1.66-fold increase in risk for CVD in the FINRISK study (HR = 1.66, p = 0.0026).

CONCLUSIONS

The risk of MI and CVD conferred by Chromosome 9p21 SNPs appears to be modified by a prudent diet high in raw vegetables and fruits. Please see later in the article for the Editors' Summary.

The HDL proteome in acute coronary syndromes shifts to an inflammatory profile

Inflammation is a major factor underlying acute coronary syndromes (ACS). HDL particles may be remodeled, becoming functionally defective, under the inflammatory conditions seen in ACS. Shotgun proteomics was used to monitor changes in the HDL proteome between male age-matched control, stable CAD, and ACS subjects (n=10/group). HDL was isolated by ultracentrifugation and separated by 1D-gel followed by LC-MS/MS. We identified 67 HDL-associated proteins, 20 of which validated recently identified proteins including vitronectin and complement C4B, and 5 of which were novel. Using gene ontology analysis, we found that the HDL-proteome consisted of proteins involved in cholesterol homeostasis (~50%), with significant contributions by proteins involved in lipid binding, antioxidant, acute-phase response, immune response, and endopeptidase/protease inhibition. Importantly, levels of apoA-IV were significantly reduced in ACS patients, whereas levels of serum amyloid A (SAA) and complement C3 (C3) were significantly increased (spectral counting; t-test p≤0.05), as confirmed by immunoblot or ELISA. Despite differences in protein composition, ABCA1, ABCG1, and SR-BI mediated cholesterol efflux assays did not indicate that HDL from ACS patients is functionally deficient as compared to controls, when corrected for apoA-I mass. Our results support that the HDL proteome differs between control, CAD and ACS patients. Increased abundance of SAA, C3, and other inflammatory proteins in HDL from ACS patients suggests that HDL reflects a shift to an inflammatory profile which, in turn, might alter the protective effects of HDL on the atherosclerotic plaque. This article is part of a Special Issue entitled Advances in High Density Lipoprotein Formation and Metabolism: A Tribute to John F. Oram (1945-2010).

Fine Mapping of the Insulin-Induced Gene 2 Identifies a Variant Associated With LDL Cholesterol and Total Apolipoprotein B Levels

Background—

In a whole-genome scan, a single nucleotide polymorphism (SNP) (rs7566605) upstream of the insulin-induced gene 2 ( INSIG2 ) was shown to influence body mass index and obesity in the Framingham Heart Study, with replication of these results in an additional 4 of 5 studies. However, other studies could not replicate the association. Because INSIG2 plays an important role in cholesterol biosynthesis, we hypothesized that human INSIG2 variants might play a role in the regulation of plasma lipid and lipoprotein levels.

Methods and Results—

We selected tagging SNPs spanning >100 kb of INSIG2 locus and sequenced 18 434 base pairs to discover novel SNPs. Thirty-two SNPs were genotyped in 645 individuals from the Quebec Family Study. Two SNPs (rs10490626 and rs12464355) were associated with plasma low-density lipoprotein cholesterol (LDL-C) ( P <0.0015) and total apolipoprotein B (apoB) levels ( P <0.014), whereas no association was found between any SNP and body mass index. We replicated the finding of rs10490626 for both LDL-C and total apoB in additional study samples, including 758 individuals from Saguenay–Lac St. Jean, Quebec ( P =0.040 for LDL-C, P =0.044 for apoB), 3247 Europeans ( P =0.028 for LDL-C, P =0.030 for apoB), and 1695 South Asians ( P =0.0036 for LDL-C, P =0.034 for apoB) from the INTERHEART study (for LDL-C, the combined 2-sided P =6.2×10 −5 and for total apoB, P =0.0011). Furthermore, we identified a variant in the human sorbin and SH 3 -domain–containing-1 gene that was associated with INSIG2 mRNA levels, and this SNP was shown to act in combination with rs10490626 to affect LDL-C ( P =0.022) in the Quebec Family Study and in INTERHEART South Asians ( P =0.019) and Europeans ( P =0.052).

Conclusion—

These results suggest that INSIG2 genetic variants may have a more direct role in lipid and lipoprotein metabolism than in obesity.

Variation at the NFATC2 locus increases the risk of thiazolidinedione-induced edema in the Diabetes REduction Assessment with ramipril and rosiglitazone Medication (DREAM) study

OBJECTIVE

Thiazolidinediones are used to treat type 2 diabetes. Their use has been associated with peripheral edema and congestive heart failure-outcomes that may have a genetic etiology.

RESEARCH DESIGN AND METHODS

We genotyped 4,197 participants of the multiethnic DREAM (Diabetes REduction Assessment with ramipril and rosiglitazone Medication) trial with a 50k single nucleotide polymorphisms (SNP) array, which captures ∼2000 cardiovascular, inflammatory, and metabolic genes. We tested 32,088 SNPs for an association with edema among Europeans who received rosiglitazone (n = 965).

RESULTS

One SNP, rs6123045, in NFATC2 was significantly associated with edema (odds ratio 1.89 [95% CI 1.47-2.42]; P = 5.32 × 10(-7), corrected P = 0.017). Homozygous individuals had the highest edema rate (hazard ratio 2.89, P = 4.22 × 10(-4)) when compared with individuals homozygous for the protective allele, with heterozygous individuals having an intermediate risk. The interaction between the SNP and rosiglitazone for edema was significant (P = 7.68 × 10(-3)). Six SNPs in NFATC2 were significant in both Europeans and Latin Americans (P < 0.05).

CONCLUSIONS

Genetic variation at the NFATC2 locus contributes to edema among individuals who receive rosiglitazone.

Concept, design and implementation of a cardiovascular gene-centric 50 k SNP array for large-scale genomic association studies

A wealth of genetic associations for cardiovascular and metabolic phenotypes in humans has been accumulating over the last decade, in particular a large number of loci derived from recent genome wide association studies (GWAS). True complex disease-associated loci often exert modest effects, so their delineation currently requires integration of diverse phenotypic data from large studies to ensure robust meta-analyses. We have designed a gene-centric 50 K single nucleotide polymorphism (SNP) array to assess potentially relevant loci across a range of cardiovascular, metabolic and inflammatory syndromes. The array utilizes a "cosmopolitan" tagging approach to capture the genetic diversity across approximately 2,000 loci in populations represented in the HapMap and SeattleSNPs projects. The array content is informed by GWAS of vascular and inflammatory disease, expression quantitative trait loci implicated in atherosclerosis, pathway based approaches and comprehensive literature searching. The custom flexibility of the array platform facilitated interrogation of loci at differing stringencies, according to a gene prioritization strategy that allows saturation of high priority loci with a greater density of markers than the existing GWAS tools, particularly in African HapMap samples. We also demonstrate that the IBC array can be used to complement GWAS, increasing coverage in high priority CVD-related loci across all major HapMap populations. DNA from over 200,000 extensively phenotyped individuals will be genotyped with this array with a significant portion of the generated data being released into the academic domain facilitating in silico replication attempts, analyses of rare variants and cross-cohort meta-analyses in diverse populations. These datasets will also facilitate more robust secondary analyses, such as explorations with alternative genetic models, epistasis and gene-environment interactions.

Genetic variants of FTO influence adiposity, insulin sensitivity, leptin levels, and resting metabolic rate in the Quebec Family Study

OBJECTIVE

A genome-wide association study conducted by the Wellcome Trust Case Control Consortium recently associated single nucleotide polymorphisms (SNPs) in the FTO (fatso/fat mass and obesity associated) gene with type 2 diabetes. These associations were shown to be mediated by obesity. Other research groups found similar results in Europeans and Hispanics but not African Americans. The mechanism by which FTO influences obesity and type 2 diabetes is currently unknown. The present study investigated the role of two FTO SNPs (rs17817449 and rs1421085) in adiposity, insulin sensitivity, and body weight regulation, including energy intake and expenditure.

RESEARCH DESIGN AND METHODS

We genotyped 908 individuals from the Quebec City metropolitan area that participated in the Quebec Family Study, a long-term study of extensively phenotyped individuals designed to investigate factors involved in adiposity.

RESULTS

We found significant associations for both SNPs with several obesity-related phenotypes. In particular, rs17817449 was associated with BMI (P = 0.0014), weight (P = 0.0059), and waist circumference (P = 0.0021) under an additive model. In addition, this FTO SNP influenced fasting insulin (P = 0.011), homeostasis model assessment of insulin resistance (P = 0.038), and an insulin sensitivity index derived from an oral glucose tolerance test (P = 0.0091). Associations were also found with resting metabolic rate (RMR) (P = 0.042) and plasma leptin levels (P = 0.036). Adjustment for BMI abolished the associations with insulin sensitivity, RMR, and plasma leptin levels.

CONCLUSIONS

These results confirm that genetic variation at the FTO locus contributes to the etiology of obesity, insulin resistance, and increased plasma leptin levels.

Familial combined hyperlipidaemia: how can genetic disorders be common, complex and comprehensible?

FCHL (familial combined hyperlipidaemia) is characterized by multiple phenotypes that are shaped by genes, the environment and time. A longitudinal study by Brouwers and co-workers, which appears in this issue of Clinical Science, points to the central role of the liver in defining the FCHL phenotypes and demonstrates how they vary over time in relation to energy excess. On the basis of their work and that of others, we propose that FCHL is a multiple gene/multiple pathway/multiple phenotype disease. The key feature of this model of common complex disease is that it posits testable faults in definable metabolic pathways, which supply the genetic underpinning of the disorder.

Identification of a chromosome 8p locus for early-onset coronary heart disease in a French Canadian population

Susceptibility to coronary heart disease (CHD) has long been known to exhibit familial aggregation, with heritability estimated to be greater than 50%. The French Canadian population of the Saguenay-Lac Saint-Jean region of Quebec, Canada is descended from a founder population that settled this region 300-400 years ago and this may provide increased power to detect genes contributing to complex traits such as CHD. Probands with early-onset CHD, defined by angiographically determined coronary stenosis, and their relatives were recruited from this population (average sibship size of 6.4). Linkage analysis was performed following a genome-wide microsatellite marker scan on 42 families with 284 individuals. Nonparametric linkage (NPL) analysis provided suggestive evidence for a CHD susceptibility locus on chromosome 8 with an NPL score of 3.14 (P=0.001) at D8S1106. Linkage to this locus was verified by fine mapping in an enlarged sample of 50 families with 320 individuals. This analysis provided evidence of linkage at D8S552 (NPL score=3.53, P=0.0003), a marker that maps to the same location as D8S1106. Candidate genes in this region, including macrophage scavenger receptor 1, farnesyl-diphosphate farnesyltransferase 1, fibrinogen-like 1, and GATA-binding protein 4, were resequenced in all coding exons in both affected and unaffected individuals. Association studies with variants in these and five other genes did not identify a disease-associated mutation. In conclusion, a genome-wide scan and additional fine mapping provide evidence for a locus on chromosome 8 that contributes to CHD in a French Canadian population.

Common polymorphisms in the promoter of the visfatin gene (PBEF1) influence plasma insulin levels in a French-Canadian population

The adipokine visfatin (PBEF1) exhibits insulin-mimetic effects and correlates strongly with visceral adiposity. We sequenced visfatin gene exons and 1,480 bp of the promoter in 23 individuals, including 18 individuals from the Quebec Family Study (QFS) with varying degrees of abdominal visceral fat, assessed by computed tomography, and 5 individuals from the Saguenay-Lac-Saint-Jean region of Québec. We identified a synonymous polymorphism in exon 7 (SER301SER) but no nonsynonymous mutations. We observed an additional 10 polymorphisms, including 5 intronic, 4 within the promoter, and 1 within the 3' untranslated region. Further promoter sequencing (816 bp) identified five additional single nucleotide polymorphisms (SNPs) in the QFS population. To investigate the role of visfatin gene variants in obesity-related phenotypes, we genotyped a total of 13 SNPs in the promoter region of the gene. From these, we analyzed the seven common SNPs in the QFS sample (918 participants from 208 families). A significant association was found between two SNPs (rs9770242 and rs1319501), in perfect linkage disequilibrium, and fasting insulin levels (P = 0.002). These SNPs were also associated with fasting glucose (P <or= 0.02). In addition, a more distal SNP (rs7789066) was significantly associated with the apolipoprotein B component of VLDL (P = 0.012).

Effect of hip position on gravity effect torque

Hip position has been hypothesized to influence gravity effect torque (GET) at the knee during isokinetic testing; however, no data exist to support or refute this hypothesis. Therefore, the purposes of this study were 1) to determine if a significant difference exists between GET in seated and supine positions, 2) to determine the effect of the supine and seated GET on isokinetic peak torque values, and 3) to determine the relationship between hamstring flexibility and GET. Gravity effect torque was recorded in supine and seated positions. Peak torque values in flexion and extension were obtained on a isokinetic dynamometer at 1.047 and 5.235 rads.s-1 (60 and 300 degrees.s-1, respectively). Hamstring flexibility was assessed by the active knee extension test (AKET). The mean seated GET value was 5.64 Nm higher than the mean supine GET value (F(1,82) = 97.85, P = 0.0001). Significant correlations existed between hamstring flexibility and GET values measured in the seated and supine positions (r = 0.45, P = 0.0001, and r = 0.30, P = 0.0058, respectively). Significant differences in peak torque values occurred for three of the four isokinetic conditions when using different GET values (P-value 0.0002-0.0049). Although mean differences in peak torque values were only 2.43-4.23 Nm, these differences may translate to significant errors in the isokinetic measurement of the injured population undergoing rehabilitation. Furthermore, every attempt should be made to improve the validity of isokinetic testing. Therefore, we recommend the supine position for GET determination.