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.

GWAS meta-analyses clarify the genetics of cervical phenotypes and inform risk stratification for cervical cancer.

Genome-wide association studies (GWAS) have successfully identified associations for cervical cancer, but the underlying mechanisms of cervical biology and pathology remain uncharacterised. Our GWAS meta-analyses fill this gap, as we characterise the genetic architecture of cervical phenotypes, including cervical ectropion, cervicitis, cervical dysplasia, as well as up to 9229 cases and 490 304 controls for cervical cancer from diverse ancestries. Leveraging the latest computational methods and gene expression data, we refine the association signals for cervical cancer and propose potential causal variants and genes at each locus. We prioritise PAX8/PAX8-AS1, LINC00339, CDC42, CLPTM1L, HLA-DRB1 and GSDMB as the most likely candidate genes for cervical cancer signals, providing insights into cervical cancer pathogenesis and supporting the involvement of reproductive tract development, immune response and cellular proliferation/apoptosis. We construct a genetic risk score (GRS) that is associated with cervical cancer [hazard ratios (HR) = 3.1 (1.7-5.6) for the top 15% vs lowest 15% of individuals], and with other HPV- and immune-system-related diagnoses in a phenome-wide association study analysis. Our results propose valuable leads for further functional studies and present a GRS for cervical cancer that allows additional risk stratification and could potentially be used to personalise the conventional screening strategies for groups more susceptible to cervical cancer.

Genome-Wide Association Study Identifies 4 Novel Risk Loci for Small Intestinal Neuroendocrine Tumors Including a Missense Mutation in LGR5.

BACKGROUND & AIMS: Small intestinal neuroendocrine tumor (SI-NET) is a rare disease, but its incidence has increased over the past 4 decades. Understanding the genetic risk factors underlying SI-NETs can help in disease prevention and may provide clinically beneficial markers for diagnosis. Here the results of the largest genome-wide association study of SI-NETs performed to date with 405 cases and 614,666 controls are reported. METHODS: Samples from 307 patients with SI-NETs and 287,137 controls in the FinnGen study were used for the identification of SI-NET risk-associated genetic variants. The results were also meta-analyzed with summary statistics from the UK Biobank (n = 98 patients with SI-NET and n = 327,529 controls). RESULTS: We identified 6 genome-wide significant (P < 5 × 10-8) loci associated with SI-NET risk, of which 4 (near SEMA6A, LGR5, CDKAL1, and FERMT2) are novel and 2 (near LTA4H-ELK and in KIF16B) have been reported previously. Interestingly, the top hit (rs200138614; P = 1.80 × 10-19) was a missense variant (p.Cys712Phe) in the LGR5 gene, a bona-fide marker of adult intestinal stem cells and a potentiator of canonical WNT signaling. The association was validated in an independent Finnish collection of 70 patients with SI-NETs, as well as in the UK Biobank exome sequence data (n = 92 cases and n = 392,814 controls). Overexpression of LGR5 p.Cys712Phe in intestinal organoids abolished the ability of R-Spondin1 to support organoid growth, indicating that the mutation perturbed R-Spondin-LGR5 signaling. CONCLUSIONS: Our study is the largest genome-wide association study to date on SI-NETs and reported 4 new associated genome-wide association study loci, including a novel missense mutation (rs200138614, p.Cys712Phe) in LGR5, a canonical marker of adult intestinal stem cells.

Primary age-related tauopathy in a Finnish population-based study of the oldest old (Vantaa 85+).

AIMS: Few studies have investigated primary age-related tauopathy (PART) in a population-based setting. Here, we assessed its prevalence, genetic background, comorbidities and features of cognitive decline in an unselected elderly population. METHODS: The population-based Vantaa 85+ study includes all 601 inhabitants of Vantaa aged ≥ 85 years in 1991. Neuropathological assessment was possible in 301. Dementia (DSM IIIR criteria) and Mini-Mental State Examination (MMSE) scores were assessed at the baseline of the study and follow-ups. PART subjects were identified according to the criteria by Crary et al and were compared with subjects with mild and severe Alzheimer's disease (AD) neuropathological changes. The effects of other neuropathologies were taken into account using multivariate and sensitivity assays. Genetic analyses included APOE genotypes and 29 polymorphisms of the MAPT 3' untranslated region (3'UTR region). RESULTS: The frequency of PART was 20% (n = 61/301, definite PART 5%). When PART subjects were compared with those with severe AD pathology, dementia was less common, its age at onset was higher and duration shorter. No such differences were seen when compared with those with milder AD pathology. However, both AD groups showed a steeper decline in MMSE scores in follow-ups compared with PART. APOE ε4 frequency was lower, and APOE ε2 frequency higher in the PART group compared with each AD group. The detected nominally significant associations between PART and two MAPT 3'UTR polymorphisms and haplotypes did not survive Bonferroni correction. CONCLUSIONS: PART is common among very elderly. PART subjects differ from individuals with AD-type changes in the pattern of cognitive decline, associated genetic and neuropathological features.

Systematic evaluation of NIPT aneuploidy detection software tools with clinically validated NIPT samples.

Non-invasive prenatal testing (NIPT) is a powerful screening method for fetal aneuploidy detection, relying on laboratory and computational analysis of cell-free DNA. Although several published computational NIPT analysis tools are available, no prior comprehensive, head-to-head accuracy comparison of the various tools has been published. Here, we compared the outcome accuracies obtained for clinically validated samples with five commonly used computational NIPT aneuploidy analysis tools (WisecondorX, NIPTeR, NIPTmer, RAPIDR, and GIPseq) across various sequencing depths (coverage) and fetal DNA fractions. The sample set included cases of fetal trisomy 21 (Down syndrome), trisomy 18 (Edwards syndrome), and trisomy 13 (Patau syndrome). We determined that all of the compared tools were considerably affected by lower sequencing depths, such that increasing proportions of undetected trisomy cases (false negatives) were observed as the sequencing depth decreased. We summarised our benchmarking results and highlighted the advantages and disadvantages of each computational NIPT software. To conclude, trisomy detection for lower coverage NIPT samples (e.g. 2.5M reads per sample) is technically possible but can, with some NIPT tools, produce troubling rates of inaccurate trisomy detection, especially in low-FF samples.

Genetic analysis of obstructive sleep apnoea discovers a strong association with cardiometabolic health.

There is currently limited understanding of the genetic aetiology of obstructive sleep apnoea (OSA). We aimed to identify genetic loci associated with OSA risk, and to test if OSA and its comorbidities share a common genetic background.We conducted the first large-scale genome-wide association study of OSA using the FinnGen study (217 955 individuals) with 16 761 OSA patients identified using nationwide health registries.We estimated 0.08 (95% CI 0.06-0.11) heritability and identified five loci associated with OSA (p<5.0×10-8): rs4837016 near GAPVD1 (GTPase activating protein and VPS9 domains 1), rs10928560 near CXCR4 (C-X-C motif chemokine receptor type 4), rs185932673 near CAMK1D (calcium/calmodulin-dependent protein kinase ID) and rs9937053 near FTO (fat mass and obesity-associated protein; a variant previously associated with body mass index (BMI)). In a BMI-adjusted analysis, an association was observed for rs10507084 near RMST/NEDD1 (rhabdomyosarcoma 2 associated transcript/NEDD1 γ-tubulin ring complex targeting factor). We found high genetic correlations between OSA and BMI (rg=0.72 (95% CI 0.62-0.83)), and with comorbidities including hypertension, type 2 diabetes, coronary heart disease, stroke, depression, hypothyroidism, asthma and inflammatory rheumatic disease (rg>0.30). The polygenic risk score for BMI showed 1.98-fold increased OSA risk between the highest and the lowest quintile, and Mendelian randomisation supported a causal relationship between BMI and OSA.Our findings support the causal link between obesity and OSA, and the joint genetic basis between OSA and comorbidities.

De novo mutations in schizophrenia implicate synaptic networks.

Inherited alleles account for most of the genetic risk for schizophrenia. However, new (de novo) mutations, in the form of large chromosomal copy number changes, occur in a small fraction of cases and disproportionally disrupt genes encoding postsynaptic proteins. Here we show that small de novo mutations, affecting one or a few nucleotides, are overrepresented among glutamatergic postsynaptic proteins comprising activity-regulated cytoskeleton-associated protein (ARC) and N-methyl-d-aspartate receptor (NMDAR) complexes. Mutations are additionally enriched in proteins that interact with these complexes to modulate synaptic strength, namely proteins regulating actin filament dynamics and those whose messenger RNAs are targets of fragile X mental retardation protein (FMRP). Genes affected by mutations in schizophrenia overlap those mutated in autism and intellectual disability, as do mutation-enriched synaptic pathways. Aligning our findings with a parallel case-control study, we demonstrate reproducible insights into aetiological mechanisms for schizophrenia and reveal pathophysiology shared with other neurodevelopmental disorders.

Comprehensive population-based genome sequencing provides insight into hematopoietic regulatory mechanisms.

Genetic variants affecting hematopoiesis can influence commonly measured blood cell traits. To identify factors that affect hematopoiesis, we performed association studies for blood cell traits in the population-based Estonian Biobank using high-coverage whole-genome sequencing (WGS) in 2,284 samples and SNP genotyping in an additional 14,904 samples. Using up to 7,134 samples with available phenotype data, our analyses identified 17 associations across 14 blood cell traits. Integration of WGS-based fine-mapping and complementary epigenomic datasets provided evidence for causal mechanisms at several loci, including at a previously undiscovered basophil count-associated locus near the master hematopoietic transcription factor CEBPA The fine-mapped variant at this basophil count association near CEBPA overlapped an enhancer active in common myeloid progenitors and influenced its activity. In situ perturbation of this enhancer by CRISPR/Cas9 mutagenesis in hematopoietic stem and progenitor cells demonstrated that it is necessary for and specifically regulates CEBPA expression during basophil differentiation. We additionally identified basophil count-associated variation at another more pleiotropic myeloid enhancer near GATA2, highlighting regulatory mechanisms for ordered expression of master hematopoietic regulators during lineage specification. Our study illustrates how population-based genetic studies can provide key insights into poorly understood cell differentiation processes of considerable physiologic relevance.

Creating basis for introducing non-invasive prenatal testing in the Estonian public health setting.

OBJECTIVE: The study aimed to validate a whole-genome sequencing-based NIPT laboratory method and our recently developed NIPTmer aneuploidy detection software with the potential to integrate the pipeline into prenatal clinical care in Estonia. METHOD: In total, 424 maternal blood samples were included. Analysis pipeline involved cell-free DNA extraction, library preparation and massively parallel sequencing on Illumina platform. Aneuploidies were determined with NIPTmer software, which is based on counting pre-defined per-chromosome sets of unique k-mers from sequencing raw data. SeqFF was implemented to estimate cell-free fetal DNA (cffDNA) fraction. RESULTS: NIPTmer identified correctly all samples of non-mosaic trisomy 21 (T21, 15/15), T18 (9/9), T13 (4/4) and monosomy X (4/4) cases, with the 100% sensitivity. However, one mosaic T18 remained undetected. Six false-positive (FP) results were observed (FP rate of 1.5%, 6/398), including three for T18 (specificity 99.3%) and three for T13 (specificity 99.3%). The level of cffDNA of <4% was estimated in eight samples, including one sample with T13 and T18. Despite low cffDNA level, these two samples were determined as aneuploid. CONCLUSION: We believe that the developed NIPT method can successfully be used as a universal primary screening test in combination with ultrasound scan for the first trimester fetal examination.

Epigenetic DNA methylation changes associated with headache chronification: A retrospective case-control study.

Background The biological mechanisms of headache chronification are poorly understood. We aimed to identify changes in DNA methylation associated with the transformation from episodic to chronic headache. Methods Participants were recruited from the population-based Norwegian HUNT Study. Thirty-six female headache patients who transformed from episodic to chronic headache between baseline and follow-up 11 years later were matched against 35 controls with episodic headache. DNA methylation was quantified at 485,000 CpG sites, and changes in methylation level at these sites were compared between cases and controls by linear regression analysis. Data were analyzed in two stages (Stages 1 and 2) and in a combined meta-analysis. Results None of the top 20 CpG sites identified in Stage 1 replicated in Stage 2 after multiple testing correction. In the combined meta-analysis the strongest associated CpG sites were related to SH2D5 and NPTX2, two brain-expressed genes involved in the regulation of synaptic plasticity. Functional enrichment analysis pointed to processes including calcium ion binding and estrogen receptor pathways. Conclusion In this first genome-wide study of DNA methylation in headache chronification several potentially implicated loci and processes were identified. The study exemplifies the use of prospectively collected population cohorts to search for epigenetic mechanisms of disease.

Distribution and medical impact of loss-of-function variants in the Finnish founder population.

Exome sequencing studies in complex diseases are challenged by the allelic heterogeneity, large number and modest effect sizes of associated variants on disease risk and the presence of large numbers of neutral variants, even in phenotypically relevant genes. Isolated populations with recent bottlenecks offer advantages for studying rare variants in complex diseases as they have deleterious variants that are present at higher frequencies as well as a substantial reduction in rare neutral variation. To explore the potential of the Finnish founder population for studying low-frequency (0.5-5%) variants in complex diseases, we compared exome sequence data on 3,000 Finns to the same number of non-Finnish Europeans and discovered that, despite having fewer variable sites overall, the average Finn has more low-frequency loss-of-function variants and complete gene knockouts. We then used several well-characterized Finnish population cohorts to study the phenotypic effects of 83 enriched loss-of-function variants across 60 phenotypes in 36,262 Finns. Using a deep set of quantitative traits collected on these cohorts, we show 5 associations (p<5×10⁻8) including splice variants in LPA that lowered plasma lipoprotein(a) levels (P = 1.5×10⁻¹¹7). Through accessing the national medical records of these participants, we evaluate the LPA finding via Mendelian randomization and confirm that these splice variants confer protection from cardiovascular disease (OR = 0.84, P = 3×10⁻4), demonstrating for the first time the correlation between very low levels of LPA in humans with potential therapeutic implications for cardiovascular diseases. More generally, this study articulates substantial advantages for studying the role of rare variation in complex phenotypes in founder populations like the Finns and by combining a unique population genetic history with data from large population cohorts and centralized research access to National Health Registers.

Genome-wide methylation analyses of primary human leukocyte subsets identifies functionally important cell-type-specific hypomethylated regions.

DNA methylation is an important mechanism by which gene transcription and hence cellular function are regulated. Here, we provide detailed functional genome-wide methylome maps of 5 primary peripheral blood leukocyte subsets including T cells, B cells, monocytes/macrophages, and neutrophils obtained from healthy individuals. A comparison of these methylomes revealed highly specific cell-lineage and cell-subset methylation profiles. DNA hypomethylation is known to be permissive for gene expression and we identified cell-subset-specific hypomethylated regions (HMRs) that strongly correlate with gene transcription levels suggesting these HMRs may regulate corresponding cell functions. Single-nucleotide polymorphisms associated with immune-mediated disease in genome-wide association studies preferentially localized to these cell-specific regulatory HMRs, offering insight into pathogenesis by highlighting cell subsets in which specific epigenetic changes may drive disease. Our data provide a valuable reference tool for researchers aiming to investigate the role of DNA methylation in regulating primary leukocyte function in health and immune-mediated disease.

Structural genomic variation as risk factor for idiopathic recurrent miscarriage.

Recurrent miscarriage (RM) is a multifactorial disorder with acknowledged genetic heritability that affects ∼3% of couples aiming at childbirth. As copy number variants (CNVs) have been shown to contribute to reproductive disease susceptibility, we aimed to describe genome-wide profile of CNVs and identify common rearrangements modulating risk to RM. Genome-wide screening of Estonian RM patients and fertile controls identified excessive cumulative burden of CNVs (5.4 and 6.1 Mb per genome) in two RM cases possibly increasing their individual disease risk. Functional profiling of all rearranged genes within RM study group revealed significant enrichment of loci related to innate immunity and immunoregulatory pathways essential for immune tolerance at fetomaternal interface. As a major finding, we report a multicopy duplication (61.6 kb) at 5p13.3 conferring increased maternal risk to RM in Estonia and Denmark (meta-analysis, n = 309/205, odds ratio = 4.82, P = 0.012). Comparison to Estonian population-based cohort (total, n = 1000) confirmed the risk for Estonian female cases (P = 7.9 × 10(-4) ). Datasets of four cohorts from the Database of Genomic Variants (total, n = 5,846 subjects) exhibited similar low duplication prevalence worldwide (0.7%-1.2%) compared to RM cases of this study (6.6%-7.5%). The CNV disrupts PDZD2 and GOLPH3 genes predominantly expressed in placenta and it may represent a novel risk factor for pregnancy complications.

Genetic Study of Psoriasis Highlights its Close Link with Socioeconomic Status and Affective Symptoms.

Psoriasis is an inflammatory skin disease with an estimated heritability of around 70%. Previous GWASs have detected several risk loci for psoriasis. To further improve the understanding of the genetic risk factors impacting the disease, we conducted a discovery GWAS in FinnGen and a subsequent replication and meta-analysis with data from the Estonian Biobank and the UK Biobank; the study sample included 925,649 individuals (22,659 cases and 902,990 controls), the largest sample for psoriasis yet. In addition, we conducted downstream analyses to find out more about psoriasis' cross-trait genetic correlations and causal relationships. We report 6 risk loci, which, to our knowledge, are previously unreported, most of which harbor genes related to NF-κB signaling pathway and overall immunity. Genetic correlations highlight the relationship between psoriasis and smoking, higher body weight, and lower education level. In addition, we report causal relationships between psoriasis and mood symptoms as well as 2-directioned causal relationship between psoriasis and lower education level. Our results provide further knowledge on psoriasis risk factors, which may be useful in the development of future treatment strategies.

Job-related exhaustion risk variant in UST is associated with dementia and DNA methylation.

Previous genome-wide association and replication study for job-related exhaustion indicated a risk variant, rs13219957 in the UST gene. Epidemiological studies suggest connection of stress-related conditions and dementia risk. Therefore, we first studied association of rs13219957 and register-based incident dementia using survival models in the Finnish National FINRISK study surveys (N = 26,693). The AA genotype of rs13219957 was significantly associated with 40% increased risk of all-cause dementia. Then we analysed the UST locus association with brain pathology in the Vantaa 85+ cohort and found association with tau pathology (Braak stage) but not with amyloid pathology. Finally, in the functional analyses, rs13219957 showed a highly significant association with two DNA methylation sites of UST, and UST expression. Thus, the results suggest a common risk variant for a stress-related condition and dementia. Mechanisms to mediate the connection may include differential DNA methylation and transcriptional regulation of UST.

Distinct and shared genetic architectures of gestational diabetes mellitus and type 2 diabetes.

Gestational diabetes mellitus (GDM) is a common metabolic disorder affecting more than 16 million pregnancies annually worldwide1,2. GDM is related to an increased lifetime risk of type 2 diabetes (T2D)1-3, with over a third of women developing T2D within 15 years of their GDM diagnosis. The diseases are hypothesized to share a genetic predisposition1-7, but few studies have sought to uncover the genetic underpinnings of GDM. Most studies have evaluated the impact of T2D loci only8-10, and the three prior genome-wide association studies of GDM11-13 have identified only five loci, limiting the power to assess to what extent variants or biological pathways are specific to GDM. We conducted the largest genome-wide association study of GDM to date in 12,332 cases and 131,109 parous female controls in the FinnGen study and identified 13 GDM-associated loci, including nine new loci. Genetic features distinct from T2D were identified both at the locus and genomic scale. Our results suggest that the genetics of GDM risk falls into the following two distinct categories: one part conventional T2D polygenic risk and one part predominantly influencing mechanisms disrupted in pregnancy. Loci with GDM-predominant effects map to genes related to islet cells, central glucose homeostasis, steroidogenesis and placental expression.

Targeted gene expression profiling for accurate endometrial receptivity testing.

Expressional profiling of the endometrium enables the personalised timing of the window of implantation (WOI). This study presents and evaluates a novel analytical pipeline based on a TAC-seq (Targeted Allele Counting by sequencing) method for endometrial dating. The expressional profiles were clustered, and differential expression analysis was performed on the model development group, using 63 endometrial biopsies spanning over proliferative (PE, n = 18), early-secretory (ESE, n = 18), mid-secretory (MSE, n = 17) and late-secretory (LSE, n = 10) endometrial phases of the natural cycle. A quantitative predictor model was trained on the development group and validated on sequenced samples from healthy women, consisting of 52 paired samples taken from ESE and MSE phases and five LSE phase samples from 31 individuals. Finally, the developed test was applied to 44 MSE phase samples from a study group of patients diagnosed with recurrent implantation failure (RIF). In validation samples (n = 57), we detected displaced WOI in 1.8% of the samples from fertile women. In the RIF study group, we detected a significantly higher proportion of the samples with shifted WOI than in the validation set of samples from fertile women, 15.9% and 1.8% (p = 0.012), respectively. The developed model was evaluated with an average cross-validation accuracy of 98.8% and an accuracy of 98.2% in the validation group. The developed beREADY screening model enables sensitive and dynamic detection of selected transcriptome biomarkers, providing a quantitative and accurate prediction of endometrial receptivity status.