Genome-wide association studies (GWAS) have identified a large number of disease-associated SNPs, but in few cases the functional variant and the gene it controls have been identified. To systematically identify candidate regulatory variants, we sequenced ENCODE cell lines and used public ChIP-seq data to look for transcription factors binding preferentially to one allele. We found 9962 candidate regulatory SNPs, of which 16 % were rare and showed evidence of larger functional effect than common ones. Functionally rare variants may explain divergent GWAS results between populations and are candidates for a partial explanation of the missing heritability. The majority of allele-specific variants (96 %) were specific to a cell type. Furthermore, by examining GWAS loci we found >400 allele-specific candidate SNPs, 141 of which were highly relevant in our cell types. Functionally validated SNPs support identification of an SNP in SYNGR1 which may expose to the risk of rheumatoid arthritis and primary biliary cirrhosis, as well as an SNP in the last intron of COG6 exposing to the risk of psoriasis. We propose that by repeating the ChIP-seq experiments of 20 selected transcription factors in three to ten people, the most common polymorphisms can be interrogated for allele-specific binding. Our strategy may help to remove the current bottleneck in functional annotation of the genome.Electronic supplementary materialThe online version of this article (doi:10.1007/s00439-016-1654-x) contains supplementary material, which is available to authorized users.

To further explore the extent of structural large-scale variation in the human genome, we assessed copy number variations (CNVs) in a series of 71 healthy subjects from three ethnic groups. CNVs were analyzed using comparative genomic hybridization (CGH) to a BAC array covering the human genome, using DNA extracted from peripheral blood, thus avoiding any culture-induced rearrangements. By applying a newly developed computational algorithm based on Hidden Markov modeling, we identified 1,078 autosomal CNVs, including at least two neighboring/overlapping BACs, which represent 315 distinct regions. The average size of the sequence polymorphisms was approximately 350 kb and involved in total approximately 117 Mb or approximately 3.5% of the genome. Gains were about four times more common than deletions, and segmental duplications (SDs) were overrepresented, especially in larger deletion variants. This strengthens the notion that SDs often define hotspots of chromosomal rearrangements. Over 60% of the identified autosomal rearrangements match previously reported CNVs, recognized with various platforms. However, results from chromosome X do not agree well with the previously annotated CNVs. Furthermore, data from single BACs deviating in copy number suggest that our above estimate of total variation is conservative. This report contributes to the establishment of the common baseline for CNV, which is an important resource in human genetics.

Pheochromocytomas and abdominal paragangliomas are adrenal and extra-adrenal catecholamine-producing tumours. They arise due to heritable cancer syndromes, or more frequently occur sporadically due to an unknown genetic cause. The majority of cases are benign, but malignant tumours are observed. Previous comparative genomic hybridization (CGH) and loss of heterozygosity studies have shown frequent deletions of chromosome arms 1p, 3q and 22q in pheochromocytomas. We applied high-resolution whole-genome array CGH on 53 benign and malignant pheochromocytomas and paragangliomas to narrow down candidate regions as well as to identify chromosomal alterations more specific to malignant tumours. Minimal overlapping regions (MORs) were identified on 16 chromosomes, with the most frequent MORs of deletion (R32%) occurring on chromosome arms 1p, 3q, 11p/q, 17p and 22q, while the chromosome arms 1q, 7p, 12q and 19p harboured the most common MORs of gain (R14%). The most frequent MORs (61-75%) in the pheochromocytomas were identified at 1p, and the four regions of common losses encompassed 1p36, 1p32-31, 1p22-21 and 1p13. Tumours that did not show 1p loss generally demonstrated aberrations on chromosome 11. Gain of chromosomal material was significantly more frequent among the malignant cases. Moreover, gain at 19q, trisomy 12 and loss at 11q were positively associated with malignant pheochromocytomas, while 1q gain was commonly observed in the malignant paragangliomas. Our study revealed novel and narrow recurrent chromosomal regions of loss and gain at several autosomes, a prerequisite for identifying candidate tumour suppressor genes and oncogenes involved in the development of adrenal and extra-adrenal catecholamine-producing tumours.

Urinary bladder cancer is a heterogeneous disease with tumors ranging from papillary noninvasive (stage Ta) to solid muscle infiltrating tumors (stage T21). The risk of progression and death for the most frequent diagnosed type, Ta, is low, but the high incidence of recurrences has a significant effect on the patients' quality of life and poses substantial costs for health care systems. Consequently, the purpose of this study was to search for predictive factors of recurrence on the basis of genetic profiling. A clinically well characterized cohort of Ta bladder carcinomas, selected by the presence or absence of recurrences, was evaluated by an integrated analysis of DNA copy number changes and gene expression (clone-based 32K, respectively, U133Plus2.0 arrays). Only a few chromosomal aberrations have previously been defined in superficial bladder cancer. Surprisingly, the profiling of Ta tumors with a high-resolution array showed that DNA copy alterations are relatively common in this tumor type. Furthermore, we observed an overrepresentation of focal amplifications within high-grade and recurrent cases. Known (FGFR3, CCND1, MYC, MDM2) and novel candidate genes were identified within the loci. For example, MYBL2, a nuclear transcription factor involved in cell-cycle progression; YWHAB, an antiapoptotic protein; and SDC4, an important component of focal adhesions represent interesting candidates detected within two amplicons on chromosome 20, for which DNA amplification correlated with transcript up-regulation. The observed overrepresentation of amplicons within high-grade and recurrent cases may be clinically useful for the identification of patients who will benefit from a more aggressive therapy.Urinary bladder cancer is a common malignant disease that ranks fourth among all cancers in Europe with 91,000 new cases and 37,000 deaths annually. 1 The prevalence is 3 to 8 times higher than the incidence, making bladder cancer one of the most prevalent neoplasms, and hence, a major burden for all health care systems. The biology of this disease is heterogeneous with tumors ranging from papillary noninvasive (stage Ta) to solid muscle infiltrating high-grade tumors (stage T2þ). The histologic grade usually parallels the stage, and thus, low-grade is common in lower, whereas high-grade predominates in higher stages. The most frequent form of all newly diagnosed cancers is of the stage Ta category, usually of low grade, which constitute more than half of all newly diagnosed cases. These tumors are considered to evolve from urothelial hyperplasia and are often multifocal. After initial transurethral resection (TUR), approximately 70% recur in the bladder, which makes this tumor type responsible for the high prevalence rate. The risk of progression and death is small, but the frequency of recurrences has a significant effect on the patients' quality of life and poses a substantial cost for the health care systems. Consequently, an attempt to prevent recurrences is frequently made by intravesical instillations either by...

Glioblastomas (GBs) are malignant CNS tumors often associated with devastating symptoms. Patients with GB have a very poor prognosis, and despite treatment, most of them die within 12 months from diagnosis. Several pathways, such as the RAS, tumor protein 53 (TP53), and phosphoinositide kinase 3 (PIK3) pathways, as well as the cell cycle control pathway, have been identified to be disrupted in this tumor. However, emerging data suggest that these aberrations represent only a fraction of the genetic changes involved in gliomagenesis. In this study, we have applied a 32K clone-based genomic array, covering 99% of the current assembly of the human genome, to the detailed genetic profiling of a set of 78 GBs. Complex patterns of aberrations, including high and narrow copy number amplicons, as well as a number of homozygously deleted loci, were identified. Amplicons that varied both in number (three on average) and in size (1.4 Mb on average) were frequently detected (81% of the samples). The loci encompassed not only previously reported oncogenes (EGFR, PDGFRA, MDM2, and CDK4) but also numerous novel oncogenes as GRB10, MKLN1, PPARGC1A, HGF, NAV3, CNTN1, SYT1, and ADAMTSL3. BNC2, PTPLAD2, and PTPRE, on the other hand, represent novel candidate tumor suppressor genes encompassed within homozygously deleted loci. Many of these genes are already linked to several forms of cancer; others represent new candidate genes that may serve as prognostic markers or even as therapeutic targets in the future. The large individual variation observed between the samples demonstrates the underlying complexity of the disease and strengthens the demand for an individualized therapy based on the genetic profile of the patient.

The FADS1 and FADS2 genes in the FADS cluster encode the rate-limiting enzymes in the synthesis of long-chain polyunsaturated fatty acids (LC-PUFAs). Genetic variation in this region has been associated with a large number of diseases and traits many of them correlated to differences in metabolism of PUFAs. However, the causative variants leading to these associations have not been identified. Here we find that the multiallelic rs174557 located in an AluYe5 element in intron 1 of FADS1 is functional and lies within a PATZ1 binding site. The derived allele of rs174557, which is the common variant in most populations, diminishes binding of PATZ1, a transcription factor conferring allele-specific downregulation of FADS1. The PATZ1 binding site overlaps with a SP1 site. The competitive binding between the suppressive PATZ1 and the activating complex of SP1 and SREBP1c determines the enhancer activity of this region, which regulates expression of FADS1.

Genome-wide association studies (GWAS) point to regions with associated genetic variants but rarely to a specific gene and therefore detailed knowledge regarding the genes contributing to complex traits and diseases remains elusive. The functional role of GWAS-SNPs is also affected by linkage disequilibrium with many variants on the same haplotype and sometimes in the same regulatory element almost equally likely to mediate the effect. Using ChIP-seq data on many transcription factors, we pinpointed genetic variants in HepG2 and HeLa-S3 cell lines which show a genome-wide significant difference in binding between alleles. We identified a collection of 3713 candidate functional regulatory variants many of which are likely drivers of GWAS signals or genetic difference in expression. A recent study investigated many variants before finding the functional ones at the GALNT2 locus, which we found in our genome-wide screen in HepG2. This illustrates the efficiency of our approach.

Background: Dissemination of nonesmall-cell lung cancer (NSCLC) in the central nervous system is a frequent and challenging clinical problem. Systemic or local therapies rarely prolong survival and have modest activity regarding local control. Alterations in gene expression in brain metastasis versus primary tumour may increase aggressiveness and impair therapeutic efforts. Methods: We identified 25 patients with surgically removed NSCLC brain metastases in two different patient cohorts. For 13 of these patients, primary tumour samples were available. Gene expression analysis using the nCounter â PanCancer Immune Profiling gene expression panel (nanoString technologies Inc.) was performed in brain metastases and primary tumour samples. Identification of differentially expressed genes was conducted on normalized data using the nSolver analysis software.