FinnGen provides genetic insights from a well-phenotyped isolated population.

Population isolates such as those in Finland benefit genetic research because deleterious alleles are often concentrated on a small number of low-frequency variants (0.1% ≤ minor allele frequency < 5%). These variants survived the founding bottleneck rather than being distributed over a large number of ultrarare variants. Although this effect is well established in Mendelian genetics, its value in common disease genetics is less explored1,2. FinnGen aims to study the genome and national health register data of 500,000 Finnish individuals. Given the relatively high median age of participants (63 years) and the substantial fraction of hospital-based recruitment, FinnGen is enriched for disease end points. Here we analyse data from 224,737 participants from FinnGen and study 15 diseases that have previously been investigated in large genome-wide association studies (GWASs). We also include meta-analyses of biobank data from Estonia and the United Kingdom. We identified 30 new associations, primarily low-frequency variants, enriched in the Finnish population. A GWAS of 1,932 diseases also identified 2,733 genome-wide significant associations (893 phenome-wide significant (PWS), P < 2.6 × 10-11) at 2,496 (771 PWS) independent loci with 807 (247 PWS) end points. Among these, fine-mapping implicated 148 (73 PWS) coding variants associated with 83 (42 PWS) end points. Moreover, 91 (47 PWS) had an allele frequency of <5% in non-Finnish European individuals, of which 62 (32 PWS) were enriched by more than twofold in Finland. These findings demonstrate the power of bottlenecked populations to find entry points into the biology of common diseases through low-frequency, high impact variants.

Impact of the pre-examination phase on multicenter metabolomic studies.

The development of metabolomics in clinical applications has been limited by the lack of validation in large multicenter studies. Large population cohorts and their biobanks are a valuable resource for acquiring insights into molecular disease mechanisms. Nevertheless, most of their collections are not tailored for metabolomics and have been created without specific attention to the pre-analytical requirements for high-quality metabolome assessment. Thus, comparing samples obtained by different pre-analytical procedures remains a major challenge. Here, 1H NMR-based analyses are used to demonstrate how human serum and plasma samples collected with different operating procedures within several large European cohort studies from the Biobanking and Biomolecular Resources Infrastructure - Large Prospective Cohorts (BBMRI-LPC) consortium can be easily revealed by supervised multivariate statistical analyses at the initial stages of the process, to avoid biases in the downstream analysis. The inter-biobank differences are discussed in terms of deviations from the validated CEN/TS 16945:2016 / ISO 23118:2021 norms. It clearly emerges that biobanks must adhere to the evidence-based guidelines in order to support wider-scale application of metabolomics in biomedicine, and that NMR spectroscopy is informative in comparing the quality of different sample sources in multi cohort/center studies.

Observational study on variability between biobanks in the estimation of DNA concentration.

BACKGROUND: There is little confidence in the consistency of estimation of DNA concentrations when samples move between laboratories. Evidence on this consistency is largely anecdotal. Therefore there is a need first to measure this consistency among different laboratories and then identify and implement remedies. A pilot experiment to test logistics and provide initial data on consistency was therefore conceived. METHODS: DNA aliquots at nominal concentrations between 10 and 300 ng/mul were dispensed into the wells of 96-well plates by one participant - the coordinating centre. Participants estimated the concentration in each well and returned estimates to the coordinating centre. RESULTS: Considerable overall variability was observed among estimates. There were statistically significant differences between participants' measurements and between fluorescence emission and absorption spectroscopy. CONCLUSION: Anecdotal evidence of variability in DNA concentration estimation has been substantiated. Reduction in variability between participants will require the identification of major sources of variation, specification of effective remedies and their implementation.

Functional diversity of human protection of telomeres 1 isoforms in telomere protection and cellular senescence.

Protection of telomeres 1 (POT1) proteins in various organisms bind telomeres and regulate their structure and function. In contrast to mice carrying two distinct POT1 genes encoding two POT1 proteins (POT1a and POT1b), humans have the single POT1 gene. In addition to full-length POT1 protein (variant v1), the human POT1 gene encodes four other variants due to alternative RNA splicing (variants v2, v3, v4, and v5), whose functions are poorly understood. The functional analyses of the NH(2)-terminally and COOH-terminally truncated POT1 variants in this study showed that neither the single-stranded telomere-binding ability of the NH(2)-terminal oligonucleotide-binding (OB) folds nor the telomerase-dependent telomere elongation activity mediated by the COOH-terminal TPP1-interacting domain was telomere protective by itself. Importantly, a COOH-terminally truncated variant (v5), which consists of the NH(2)-terminal OB folds and the central region of unknown function, was found to protect telomeres and prevent cellular senescence as efficiently as v1. Our data revealed mechanistic and functional differences between v1 and v5: (a) v1, but not v5, functions through the maintenance of telomeric 3' overhangs; (b) p53 is indispensable to v5 knockdown-induced senescence; and (c) v5 functions at only a fraction of telomeres to prevent DNA damage signaling. Furthermore, v5 was preferentially expressed in mismatch repair (MMR)-deficient cells and tumor tissues, suggesting its role in chromosome stability associated with MMR deficiency. This study highlights a human-specific complexity in telomere protection and damage signaling conferred by functionally distinct isoforms from the single POT1 gene.

Increased HIF1 alpha in SDH and FH deficient tumors does not cause microsatellite instability.

Germline mutations in nuclear genes encoding mitochondrial enzymes fumarate hydratase (FH) and succinate dehydrogenase (subunits SDHB/C/D) have been implicated in the development of tumor syndromes referred to as hereditary leiomyomatosis and renal cell cancer (HLRCC) and hereditary paragangliomatosis (HPGL), respectively. FH and SDH are operating in the tricarboxylic acid cycle (the TCA cycle, the Krebs cycle). In the FH and SDH deficient tumors, accumulation of the substrates, fumarate and succinate, has been shown to cause stabilization of hypoxia inducible factor 1 alpha (HIF1 alpha). According to recent studies, HIF1 alpha could contribute to the hypoxia induced genomic instability seen in many cancers, through repression of mismatch repair (MMR) protein MSH2. In this study, in agreement with previous works, we found HIF1 alpha to be moderately or highly stabilized in 67% (16/24) and 77% (48/62) of HLRCC tumors and SDHB/C/D paragangliomas (PGL) and pheochromocytomas (PHEO), respectively. In addition, a set of 54 other familial and nonfamilial PGLs/PHEOs were studied. Moderately or highly stabilized HIF1 alpha was present in 68% (26/38) of the PGLs but in PHEOs (n = 16) no such pattern was observed. We then analyzed the suggested link between HIF1 alpha stabilization and MSH2 repression, in HLRCC and HPGL tumor material. No microsatellite instability (MSI) or lack of MSH2 expression was, however, observed. Thus we failed to provide in vivo evidence for the proposed link between HIF1 alpha stabilization and functional MMR deficiency, in TCAC deficient tumors.

EPHB2 germline variants in patients with colorectal cancer or hyperplastic polyposis.

BACKGROUND: Ephrin receptor B2 (EPHB2) has recently been proposed as a novel tumor suppressor gene in colorectal cancer (CRC). Inactivation of the gene has been shown to correlate with progression of colorectal tumorigenesis, and somatic mutations have been reported in both colorectal and prostate tumors. METHODS: Here we have analyzed the EPHB2 gene for germline alterations in 101 individuals either with 1) CRC and a personal or family history of prostate cancer (PC), or 2) intestinal hyperplastic polyposis (HPP), a condition associated with malignant degeneration such as serrated adenoma and CRC. RESULTS: Four previously unknown missense alterations were observed, which may be associated with the disease phenotype. Two of the changes, I361V and R568W, were identified in Finnish CRC patients, but not in over 300 Finnish familial CRC or PC patients or more than 200 population-matched healthy controls. The third change, D861N, was observed in a UK HPP patient, but not in additional 40 UK HPP patients or in 200 UK healthy controls. The fourth change R80H, originally identified in a Finnish CRC patient, was also found in 1/106 familial CRC patients and in 9/281 healthy controls and is likely to be a neutral polymorphism. CONCLUSION: We detected novel germline EPHB2 alterations in patients with colorectal tumors. The results suggest a limited role for these EPHB2 variants in colon tumor predisposition. Further studies including functional analyses are needed to confirm this.

Distinct expression profile in fumarate-hydratase-deficient uterine fibroids.

Defects in mitochondrial enzymes predispose to severe developmental defects as well as tumorigenesis. Heterozygous germline mutations in the nuclear gene encoding fumarate hydratase (FH), an enzyme catalyzing the hydration of fumarate in the Krebs tricarboxylic acid cycle, cause hereditary leiomyomatosis and renal cell cancer; yet the connection between disruption of mitochondrial metabolic pathways and neoplasia remains to be discovered. We have used an expression microarray approach for studying differences in global gene expression pattern caused by mutations in FH. Seven uterine fibroids carrying FH mutations were compared with 15 fibroids with wild-type FH. The two groups showed markedly different expression profiles, and multiple differentially expressed genes were detected. The most significant increase in FH mutants was seen in the expression of carbohydrate metabolism- and glycolysis-related genes. Other significantly up-regulated gene categories in FH mutants were, for example, iron ion homeostasis and oxidoreduction. Genes with lower expression in FH-mutant fibroids belonged to groups such as extracellular matrix, cell adhesion, muscle development and cell contraction. We show that FH mutations alter significantly the expression profiles of fibroids, most strikingly increasing the expression of genes involved in glycolysis.

Mechanisms of inactivation of the receptor tyrosine kinase EPHB2 in colorectal tumors.

The receptor tyrosine kinase EPHB2 has recently been shown to be a direct transcriptional target of TCF/beta-catenin. Premalignant lesions of the colon express high levels of EPHB2 but the expression of this kinase is reduced or lost in most colorectal carcinomas. In addition, inactivation of EPHB2 has been shown to accelerate tumorigenesis initiated by APC mutation in the colon and rectum. In this study, we investigated the molecular mechanisms responsible for the inactivation of EPHB2 in colorectal tumors. We show here the presence of mutations in repetitive sequences in exon 17 of EPHB2 in 6 of 29 adenomas with microsatellite instability (MSI), and 101 of 246 MSI carcinomas (21% and 41%, respectively). Moreover, we found EPHB2 promoter hypermethylation in 54 of the 101 colorectal tumors studied (53%). Importantly, EPHB2 expression was restored after treatment of EPHB2-methylated colon cancer cells with the DNA methyltransferase inhibitor 5-aza-2'-deoxycytidine. In conclusion, in this study, we elucidate the molecular mechanisms of inactivation of EPHB2 and show for the first time the high incidence of frameshift mutations in MSI colorectal tumors and aberrant methylation of the regulatory sequences of this important tumor suppressor gene.

Gene-expression profiling predicts recurrence in Dukes' C colorectal cancer.

BACKGROUND & AIMS: Although approximately 50% of Dukes' C colorectal cancer patients are surgically cured, it is currently not possible to distinguish these patients from those at high risk of recurrence. The recent advent of routine adjuvant chemotherapy for these patients has greatly complicated the identification of new markers predicting the response to surgery, which is now reliant on archived materials. Microarray analysis allows fine tumor classification but cannot be used with paraffin-embedded archival samples. METHODS: We used microarray analysis of a unique set of fresh-frozen tumor samples from Dukes' C patients who had surgery as the only form of treatment to identify molecular signatures that characterize tumors from patients with good and bad prognosis. RESULTS: Unsupervised hierarchical clustering and a K-nearest neighbors-based classifier identified groups of patients with significantly different survival (P = .019 and P = .0001). Expression profiling outperformed previously reported genetic markers of prognosis such as TP53 and K-RAS mutational status and allelic imbalance in chromosome 18q, which were of limited prognostic power in this study. Functional categories significantly enriched in gene-expression differences included protein transport and folding. The prognostic potential of the RAS homologue RHOA, one of the most differentially expressed genes, was further investigated using immunohistochemistry and a tissue microarray containing 137 independent Dukes' C tumor samples. Reduced RHOA expression was associated with significantly shorter survival (P = .01). CONCLUSIONS: This study shows that gene-expression profiling of surgical tumor samples can predict recurrence in Dukes' C patients. Therefore, this approach could be used to guide decisions concerning the clinical management of these patients.

7q deletion mapping and expression profiling in uterine fibroids.

Uterine fibroids are some of the most common tumours of females, but relatively little is known about their molecular basis. Several studies have suggested that deletions on chromosome 7q could have a role in fibroid formation. We analysed 165 sporadic uterine fibroids to define a small 3.2 megabase (Mb) commonly deleted region on 7q22.3-q31.1, flanked by clones AC005070 and AC007567. We also used oligonucleotide microarrays to compare the expression profiles of 10 samples of normal myometrium and 15 fibroids, nine of which displayed 7q-deletions. Activating transcription factor 3, patched homolog (Drosophila), homeo box A5, death-associated protein kinase 1, and retinoic acid receptor responder 3 were downregulated, and excision repair crosscomplementing 3, transcription factor AP-2 gamma and protein kinase C beta 1 were upregulated in fibroids. New pathways were discovered related to fibroid formation. The presence or absence of 7q-deletions did not dramatically affect the global expression pattern of the tumours; changes, however, were observed in genes related to vesicular transport and nucleic acid binding.

Distinct patterns of KRAS mutations in colorectal carcinomas according to germline mismatch repair defects and hMLH1 methylation status.

In sporadic colorectal tumours the BRAFV600E is associated with microsatellite instability (MSI-H) and inversely associated to KRAS mutations. Tumours from hereditary non-polyposis colorectal cancer (HNPCC) patients carrying germline mutations in hMSH2 or hMLH1 do not show BRAFV600E, however no consistent data exist regarding KRAS mutation frequency and spectrum in HNPCC tumours. We investigated KRAS in 158 HNPCC tumours from patients with germline hMLH1, hMSH2 or hMSH6 mutations, 166 MSI-H and 688 microsatellite stable (MSS) sporadic carcinomas. All tumours were characterized for MSI and 81 of 166 sporadic MSI-H colorectal cancer (CRCs) were analysed for hMLH1 promoter hypermethylation. KRAS mutations were observed in 40% of HNPCC tumours, and the mutation frequency varied upon the mismatch repair gene affected: 48% (29/61) in hMSH2, 32% (29/91) in hMLH1 and 83% (5/6) in hMSH6 (P = 0.01). KRAS mutation frequency was different between HNPCC, MSS and MSI-H CRCs (P = 0.002), and MSI-H with hMLH1 hypermethylation (P = 0.005). Furthermore, HNPCC CRCs had more G13D mutations than MSS (P < 0.0001), MSI-H (P = 0.02) or MSI-H tumours with hMLH1 hypermethylation (P = 0.03). HNPCC colorectal and sporadic MSI-H tumours without hMLH1 hypermethylation shared similar KRAS mutation frequency, in particular G13D. In conclusion, we show that depending on the genetic/epigenetic mechanism leading to MSI-H, the outcome in terms of oncogenic activation may be different, reinforcing the idea that HNPCC, sporadic MSI-H (depending on the hMLH1 status) and MSS CRCs, may target distinct kinases within the RAS/RAF/MAPK pathway.

Little evidence for involvement of MLH3 in colorectal cancer predisposition.

Mutations in the DNA MMR genes MSH2, MLH1, MSH6 and PMS2 underlie a large subset of HNPCC cases, and a hallmark of the tumors is MSI. In many HNPCC families, however, a causative mutation has not been found. Therefore, the involvement of additional, thus far unknown, genes in MSI as well as MSS colorectal tumor predisposition is possible. The role of a relatively recently cloned MMR gene, MLH3, in familial CRC has been studied; but the results appear somewhat conflicting. To further evaluate the role of MLH3 in CRC predisposition, we analyzed 30 Finnish CRC cases for germline mutations by sequencing. These cases were selected from a large series of Finnish CRC patients, to match features previously proposed to associate with MLH3 germline defects. We found 5 missense variants, 4 of which were also found in Finnish cancer-free controls. The only remaining variant does not appear to be an attractive candidate for a disease-associated mutation because the amino acid change is located outside the conserved residues. We also screened for the previously reported variants, including a frameshift change, the most likely pathogenic MLH3 mutation observed so far. The frameshift was not present in the 30 CRC cases or in 700 cancer-free controls. While it is a difficult task to exclude a role of MLH3 in HNPCC, our study could not confirm a role for MLH3 in CRC predisposition.

Genome-wide allelotyping of 104 Finnish colorectal cancers reveals an excess of allelic imbalance in chromosome 20q in familial cases.

We have allelotyped a series of 104 Finnish colorectal cancers (CRCs) using 372 polymorphic markers spaced, on average, at 10 cM intervals, and have made a comparison of the differences in the frequency of allelic imbalance (AI) between familial and sporadic cases. Differences in the frequency of allelic imbalance (loss of heterozygosity or amplification) at a number of loci were detected and these were evaluated through analysis of additional series of cancers using specific markers. The most consistent difference was observed at chromosome 20q13.1-13.3 characterized by a two fold difference between familial and nonfamilial disease in a total of 99 familial and 186 sporadic Finnish cases. This difference was not observed in a UK set of 67 familial and 96 sporadic CRCs. The genome-wide effort resulted in a large data set giving clues to the location of putative CRC predisposition genes in the genome. The approach provides an alternative strategy for detecting cancer predisposition genes solely reliant on the molecular analysis of single cases obviating the requirement to collect multiple samples from families.

Gene expression in colorectal cancer.

Understanding molecular alterations in colorectal cancer (CRC) is needed to define new biomarkers and treatment targets. We used oligonucleotide microarrays to monitor gene expression of about 6,800 known genes and 35,000 expressed sequence tags (ESTs) on five pools (four to six samples in each pool) of total RNA from left-sided sporadic colorectal carcinomas. We compared normal tissue to carcinoma tissue from Dukes' stages A-D (noninvasive to distant metastasis) and identified 908 known genes and 4,155 ESTs that changed remarkably from normal to tumor tissue. Based on intensive filtering 226 known genes and 157 ESTs were found to be highly relevant for CRC. The alteration of known genes was confirmed in >70% of the cases by array analysis of 25 single samples. Two-way hierarchical average linkage cluster analysis clustered normal tissue together with Dukes' A, clustered Dukes' B with Dukes' C, and clustered Dukes' D separately. Real-time PCR of 10 known genes and 5 ESTs demonstrated excellent reproducibility of the array-based findings. The most frequently altered genes belonged to functional categories of metabolism (22%), transcription and translation (11%), and cellular processes (9%). Fifteen nuclear encoded mitochondrial proteins were all down-regulated in CRC. We identified several chromosomal locations with clusters of either potential oncogenes or potential tumor suppressors. Some of these, such as aminopeptidase N/CD13 and sigma B3 protein on chromosome 15q25, coincided with a high frequency of loss of heterozygosity. The genes and ESTs presented in this study encode new potential tumor markers as well as potential novel therapeutic targets for prevention or therapy of CRC.

Germline mutations in FH predispose to dominantly inherited uterine fibroids, skin leiomyomata and papillary renal cell cancer.

Uterine leiomyomata (fibroids) are common and clinically important tumors, but little is known about their etiology and pathogenesis. We previously mapped a gene that predisposes to multiple fibroids, cutaneous leiomyomata and renal cell carcinoma to chromosome 1q42.3-q43 (refs 4-6). Here we show, through a combination of mapping critical recombinants, identifying individuals with germline mutations and screening known and predicted transcripts, that this gene encodes fumarate hydratase, an enzyme of the tricarboxylic acid cycle. Leiomyomatosis-associated mutations are predicted to result in absent or truncated protein, or substitutions or deletions of highly conserved amino acids. Activity of fumarate hydratase is reduced in lymphoblastoid cells from individuals with leiomyomatosis. This enzyme acts as a tumor suppressor in familial leiomyomata, and its measured activity is very low or absent in tumors from individuals with leiomyomatosis. Mutations in FH also occur in the recessive condition fumarate hydratase deficiency, and some parents of people with this condition are susceptible to leiomyomata. Thus, heterozygous and homozygous or compound heterozygous mutants have very different clinical phenotypes. Our results provide clues to the pathogenesis of fibroids and emphasize the importance of mutations of housekeeping and mitochondrial proteins in the pathogenesis of common types of tumor.

Low-level microsatellite instability in most colorectal carcinomas.

Twelve to 16% of colorectal cancers (CRCs) display a high degree of microsatellite instability (MSI-H), whereas most are believed to be microsatellite stable (MSS). The existence of a low degree of instability (MSI-L) group has also been proposed. By using the Bethesda panel of microsatellite markers, the microsatellite instability (MSI) status of CRCs can be determined. This set is recommended to distinguish between MSI-H and MSI-L/MSS. No definition for MSI-L has emerged. Most reports on MSI-L rely on the Bethesda panel, using 5-15markers. Tumors with more than 30% MSI are designated as MSI-H, but the lower limit for MSI-L is ambiguous. We hypothesized that if many markers are studied, almost all CRCs would show some MSI. It would be necessary to establish a cutoff level for MSI-L by showing that, above this cutoff level, tumors display molecular and/or clinical features different from those under the cutoff level. To perform this task, we analyzed 90 BAT26 stable CRC samples with 377 markers. MSI at 1-11 loci was observed in 71 (79%) of the 90 cases. K-RAS mutation, loss of heterozygosity, and MLH1 and MGMT hypermethylation analyses were performed, as well as clinical features being scrutinized, to examine possible differences between MSI-L and MSS tumors using all of the possible cutoff levels for MSI-L. Convincing differences between putative MSI-L and MSS groups were not observed. Our results show that the sensitivity of a typically used marker number to detect MSI-L is very low, and they suggest that MSS and MSI-L tumors have a common molecular background.