Sulfate sensitivity of early life stages of freshwater mussels Unio crassus and Margaritifera margaritifera.

Sulfate is increasingly found in elevated concentrations in freshwater ecosystems due to anthropogenic activities. Chronic exposure to sulfate has been reported to cause sublethal effects on freshwater invertebrates. Previous sulfate toxicity tests have mostly been conducted in hard or moderately hard waters, and research on species inhabiting soft water is needed, given that freshwater organisms face heightened sensitivity to toxicants in water of lower hardness. In the present study, we examined sulfate sensitivity of two endangered freshwater mussel species, Unio crassus, and Margaritifera margaritifera. Glochidia and juveniles of both species were subjected to acute and/or chronic sulfate exposures in soft water to compare sulfate sensitivity across age groups, and effective concentrations (EC)/lethal concentrations (LC) values were estimated. Mussels were individually exposed to allow relatively larger numbers of replicates per treatment. Chronic sulfate exposure significantly reduced growth, foot movement, and relative water content (RWC) in juvenile mussels of M. margaritifera. Mussels at younger stages were not necessarily more sensitive to sulfate. In the acute tests, LC50 of glochidia of M. margaritifera and U. crassus was 1301 and 857 mg/L, respectively. Chronic LC10 was 843 mg/L for 3-week-old U. crassus juveniles, 1051 mg/L for 7-week-old M. margaritifera juveniles, and 683 mg/L for 2-year-old M. margaritifera juveniles. True chronic Lowest Effective Concentration for 7-week-old M. margaritifera may be within the 95% interval of EC10 based on RWC (EC10 = 446 mg/L, 95%CI = 265-626 mg/L). Our study contributed to the understanding of sulfate toxicity to endangered freshwater mussel species in soft water.

Fish with slow life-history cope better with chronic manganese exposure than fish with fast life-history.

Animals with different life-history types vary in their stress-coping styles, which can affect their fitness and survival in changing environments. We studied how chronic exposure to manganese sulfate (MnSO4), a common aquatic pollutant, affects life-history traits, physiology, and behavior of zebrafish (Danio rerio) with two life-history types: fast (previously selected for fast juvenile growth, early maturation, and small adult body size) and slow life histories (selected for slow juvenile growth, late maturation, and large adult body size). We found that MnSO4 had negative effects on growth and condition factors, but the magnitude of these effects depended on the life-history type. Individuals with fast life histories were more susceptible to MnSO4 than fish with slow life histories as they had lower growth rate, condition factor and feeding probability in high MnSO4 concentrations. Our results demonstrate that MnSO4 can impair fish performance, and life-history variation can modulate the stress-coping ability of individuals.

High incidence and geographic distribution of cleft palate cases in Finland are associated with a regulatory variant in IRF6.

In Finland the frequency of isolated cleft palate (CP) is higher than that of isolated cleft lip with or without cleft palate (CL/P). This trend contrasts to that in other European countries but its genetic underpinnings are unknown. We performed a genome-wide association study for orofacial clefts, which include CL/P and CP, in the Finnish population. We identified rs570516915, a single nucleotide polymorphism that is highly enriched in Finns and Estonians, as being strongly associated with CP ( P = 5.25 × 10 -34 , OR = 8.65, 95% CI 6.11-12.25), but not with CL/P ( P = 7.2 × 10 -5 ), with genome-wide significance. The risk allele frequency of rs570516915 parallels the regional variation of CP prevalence in Finland, and the association was replicated in independent cohorts of CP cases from Finland ( P = 8.82 × 10 -28 ) and Estonia ( P = 1.25 × 10 -5 ). The risk allele of rs570516915 disrupts a conserved binding site for the transcription factor IRF6 within a previously characterized enhancer upstream of the IRF6 gene. Through reporter assay experiments we found that the risk allele of rs570516915 diminishes the enhancer activity. Oral epithelial cells derived from CRISPR-Cas9 edited induced pluripotent stem cells demonstrate that the CP-associated allele of rs570516915 concomitantly decreases the binding of IRF6 and the expression level of IRF6 , suggesting impaired IRF6 autoregulation as a molecular mechanism underlying the risk for CP.

Autoimmune hypothyroidism GWAS reveals independent autoimmune and thyroid-specific contributions and an inverse relation with cancer risk.

The high prevalence of autoimmune hypothyroidism (AIHT) - more than 5% in human populations - provides a unique opportunity to unlock the most complete picture to date of genetic loci that underlie systemic and organ-specific autoimmunity. Using a meta-analysis of 81,718 AIHT cases in FinnGen and the UK Biobank, we dissect associations along axes of thyroid dysfunction and autoimmunity. This largest-to-date scan of hypothyroidism identifies 418 independent associations (p < 5×10- 8), more than half of which have not previously been documented in thyroid disease. In 48 of these, a protein-coding variant is the lead SNP or is highly correlated (r2 > 0.95) with the lead SNP at the locus, including low-frequency coding variants at LAG3, ZAP70, TG, TNFSF11, IRF3, S1PR4, HABP2, ZNF429 as well as established variants at ADCY7, IFIH1 and TYK2. The variants at LAG3 (P67T), ZAP70 (T155M), and TG (Q655X) are highly enriched in Finland and functional experiments in T-cells demonstrate that the ZAP70:T155M allele reduces T-cell activation. By employing a large-scale scan of non-thyroid autoimmunity and a published meta-analysis of TSH levels, we use a Bayesian classifier to dissect the associated loci into distinct groupings and from this estimate, a significant proportion are involved in systemic (i.e., general to multiple autoimmune conditions) autoimmunity (34%) and another subset in thyroid-specific dysfunction (17%). By comparing these association results further to other common disease endpoints, we identify a noteworthy overlap with skin cancer, with 10% of AIHT loci showing a consistent but opposite pattern of association where alleles that increase the risk of hypothyroidism have protective effects for skin cancer. The association results, including genes encoding checkpoint inhibitors and other genes affecting protein levels of PD1, bolster the causal role of natural variation in autoimmunity influencing cancer outcomes.

Oral and non-oral lichen planus show genetic heterogeneity and differential risk for autoimmune disease and oral cancer.

Lichen planus (LP) is a T-cell-mediated inflammatory disease affecting squamous epithelia in many parts of the body, most often the skin and oral mucosa. Cutaneous LP is usually transient and oral LP (OLP) is most often chronic, so we performed a large-scale genetic and epidemiological study of LP to address whether the oral and non-oral subgroups have shared or distinct underlying pathologies and their overlap with autoimmune disease. Using lifelong records covering diagnoses, procedures, and clinic identity from 473,580 individuals in the FinnGen study, genome-wide association analyses were conducted on carefully constructed subcategories of OLP (n = 3,323) and non-oral LP (n = 4,356) and on the combined group. We identified 15 genome-wide significant associations in FinnGen and an additional 12 when meta-analyzed with UKBB (27 independent associations at 25 distinct genomic locations), most of which are shared between oral and non-oral LP. Many associations coincide with known autoimmune disease loci, consistent with the epidemiologic enrichment of LP with hypothyroidism and other autoimmune diseases. Notably, a third of the FinnGen associations demonstrate significant differences between OLP and non-OLP. We also observed a 13.6-fold risk for tongue cancer and an elevated risk for other oral cancers in OLP, in agreement with earlier reports that connect LP with higher cancer incidence. In addition to a large-scale dissection of LP genetics and comorbidities, our study demonstrates the use of comprehensive, multidimensional health registry data to address outstanding clinical questions and reveal underlying biological mechanisms in common but understudied diseases.

Cardiovascular Significance and Genetics of Epicardial and Pericardial Adiposity.

Importance: Epicardial and pericardial adipose tissue (EPAT) has been associated with cardiovascular diseases such as atrial fibrillation or flutter (AF) and coronary artery disease (CAD), but studies have been limited in sample size or drawn from selected populations. It has been suggested that the association between EPAT and cardiovascular disease could be mediated by local or paracrine effects. Objective: To evaluate the association of EPAT with prevalent and incident cardiovascular disease and to elucidate the genetic basis of EPAT in a large population cohort. Design, Setting, and Participants: A deep learning model was trained to quantify EPAT area from 4-chamber magnetic resonance images using semantic segmentation. Cross-sectional and prospective cardiovascular disease associations were evaluated, controlling for sex and age. Prospective associations were additionally controlled for abdominal visceral adipose tissue (VAT) volumes. A genome-wide association study was performed, and a polygenic score (PGS) for EPAT was examined in independent FinnGen cohort study participants. Data analyses were conducted from March 2022 to December 2023. Exposures: The primary exposures were magnetic resonance imaging-derived continuous measurements of epicardial and pericardial adipose tissue area and visceral adipose tissue volume. Main Outcomes and Measures: Prevalent and incident CAD, AF, heart failure (HF), stroke, and type 2 diabetes (T2D). Results: After exclusions, this study included 44 475 participants (mean [SD] age, 64.1 [7.7] years; 22 972 female [51.7%]) from the UK Biobank. Cross-sectional and prospective cardiovascular disease associations were evaluated for a mean (SD) of 3.2 (1.5) years of follow-up. Prospective associations were additionally controlled for abdominal VAT volumes for 38 527 participants. A PGS for EPAT was examined in 453 733 independent FinnGen cohort study participants. EPAT was positively associated with male sex (ß = +0.78 SD in EPAT; P < 3 × 10-324), age (Pearson r = 0.15; P = 9.3 × 10-229), body mass index (Pearson r = 0.47; P < 3 × 10-324), and VAT (Pearson r = 0.72; P < 3 × 10-324). EPAT was more elevated in prevalent HF (ß = +0.46 SD units) and T2D (ß = +0.56) than in CAD (ß = +0.23) or AF (ß = +0.18). EPAT was associated with incident HF (hazard ratio [HR], 1.29 per +1 SD in EPAT; 95% CI, 1.17-1.43), T2D (HR, 1.63; 95% CI, 1.51-1.76), and CAD (HR, 1.19; 95% CI, 1.11-1.28). However, the associations were no longer significant when controlling for VAT. Seven genetic loci were identified for EPAT, implicating transcriptional regulators of adipocyte morphology and brown adipogenesis (EBF1, EBF2, and CEBPA) and regulators of visceral adiposity (WARS2 and TRIB2). The EPAT PGS was associated with T2D (odds ratio [OR], 1.06; 95% CI, 1.05-1.07; P =3.6 × 10-44), HF (OR, 1.05; 95% CI, 1.04-1.06; P =4.8 × 10-15), CAD (OR, 1.04; 95% CI, 1.03-1.05; P =1.4 × 10-17), AF (OR, 1.04; 95% CI, 1.03-1.06; P =7.6 × 10-12), and stroke in FinnGen (OR, 1.02; 95% CI, 1.01-1.03; P =3.5 × 10-3) per 1 SD in PGS. Conclusions and Relevance: Results of this cohort study suggest that epicardial and pericardial adiposity was associated with incident cardiovascular diseases, but this may largely reflect a metabolically unhealthy adiposity phenotype similar to abdominal visceral adiposity.

Does size-selective harvesting erode adaptive potential to thermal stress?

Overharvesting is a serious threat to many fish populations. High mortality and directional selection on body size can cause evolutionary change in exploited populations via selection for a specific phenotype and a potential reduction in phenotypic diversity. Whether the loss of phenotypic diversity that accompanies directional selection impairs response to environmental stress is not known. To address this question, we exposed three zebrafish selection lines to thermal stress. Two lines had experienced directional selection for (1) large and (2) small body size, and one was (3) subject to random removal of individuals with respect to body size (i.e. line with no directional selection). Selection lines were exposed to three temperatures (elevated, 34°C; ambient, 28°C; low, 22°C) to determine the response to an environmental stressor (thermal stress). We assessed differences among selection lines in their life history (growth and reproduction), physiological traits (metabolic rate and critical thermal max) and behaviour (activity and feeding behaviour) when reared at different temperatures. Lines experiencing directional selection (i.e. size selected) showed reduced growth rate and a shift in average phenotype in response to lower or elevated thermal stress compared with fish from the random-selected line. Our data indicate that populations exposed to directional selection can have a more limited capacity to respond to thermal stress compared with fish that experience a comparable reduction in population size (but without directional selection). Future studies should aim to understand the impacts of environmental stressors on natural fish stocks.

The Cardiovascular Impact and Genetics of Pericardial Adiposity.

Background: While previous studies have reported associations of pericardial adipose tissue (PAT) with cardiovascular diseases such as atrial fibrillation and coronary artery disease, they have been limited in sample size or drawn from selected populations. Additionally, the genetic determinants of PAT remain largely unknown. We aimed to evaluate the association of PAT with prevalent and incident cardiovascular disease and to elucidate the genetic basis of PAT in a large population cohort. Methods: A deep learning model was trained to quantify PAT area from four-chamber magnetic resonance images in the UK Biobank using semantic segmentation. Cross-sectional and prospective cardiovascular disease associations were evaluated, controlling for sex and age. A genome-wide association study was performed, and a polygenic score (PGS) for PAT was examined in 453,733 independent FinnGen study participants. Results: A total of 44,725 UK Biobank participants (51.7% female, mean [SD] age 64.1 [7.7] years) were included. PAT was positively associated with male sex (ß = +0.76 SD in PAT), age (r = 0.15), body mass index (BMI; r = 0.47) and waist-to-hip ratio (r = 0.55) (P < 1×10-230). PAT was more elevated in prevalent heart failure (ß = +0.46 SD units) and type 2 diabetes (ß = +0.56) than in coronary artery disease (ß = +0.22) or AF (ß = +0.18). PAT was associated with incident heart failure (HR = 1.29 per +1 SD in PAT [95% CI 1.17-1.43]) and type 2 diabetes (HR = 1.63 [1.51-1.76]) during a mean 3.2 (±1.5) years of follow-up; the associations remained significant when controlling for BMI. We identified 5 novel genetic loci for PAT and implicated transcriptional regulators of adipocyte morphology and brown adipogenesis (EBF1, EBF2 and CEBPA) and regulators of visceral adiposity (WARS2 and TRIB2). The PAT PGS was associated with T2D, heart failure, coronary artery disease and atrial fibrillation in FinnGen (ORs 1.03-1.06 per +1 SD in PGS, P < 2×10-10). Conclusions: PAT shares genetic determinants with abdominal adiposity and is an independent predictor of incident type 2 diabetes and heart failure.

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.

Genetic Risk Factors Associated With Preeclampsia and Hypertensive Disorders of Pregnancy.

Importance: A genetic contribution to preeclampsia susceptibility has been established but is still incompletely understood. Objective: To disentangle the underlying genetic architecture of preeclampsia and preeclampsia or other maternal hypertension during pregnancy with a genome-wide association study (GWAS) of hypertensive disorders of pregnancy. Design, Setting, and Participants: This GWAS included meta-analyses in maternal preeclampsia and a combination phenotype encompassing maternal preeclampsia and preeclampsia or other maternal hypertensive disorders. Two overlapping phenotype groups were selected for examination, namely, preeclampsia and preeclampsia or other maternal hypertension during pregnancy. Data from the Finnish Genetics of Pre-eclampsia Consortium (FINNPEC, 1990-2011), Finnish FinnGen project (1964-2019), Estonian Biobank (1997-2019), and the previously published InterPregGen consortium GWAS were combined. Individuals with preeclampsia or other maternal hypertension during pregnancy and control individuals were selected from the cohorts based on relevant International Classification of Diseases codes. Data were analyzed from July 2020 to February 2023. Exposures: The association of a genome-wide set of genetic variants and clinical risk factors was analyzed for the 2 phenotypes. Results: A total of 16 743 women with prior preeclampsia and 15 200 with preeclampsia or other maternal hypertension during pregnancy were obtained from FINNPEC, FinnGen, Estonian Biobank, and the InterPregGen consortium study (respective mean [SD] ages at diagnosis: 30.3 [5.5], 28.7 [5.6], 29.7 [7.0], and 28 [not available] years). The analysis found 19 genome-wide significant associations, 13 of which were novel. Seven of the novel loci harbor genes previously associated with blood pressure traits (NPPA, NPR3, PLCE1, TNS2, FURIN, RGL3, and PREX1). In line with this, the 2 study phenotypes showed genetic correlation with blood pressure traits. In addition, novel risk loci were identified in the proximity of genes involved in the development of placenta (PGR, TRPC6, ACTN4, and PZP), remodeling of uterine spiral arteries (NPPA, NPPB, NPR3, and ACTN4), kidney function (PLCE1, TNS2, ACTN4, and TRPC6), and maintenance of proteostasis in pregnancy serum (PZP). Conclusions and Relevance: The findings indicate that genes related to blood pressure traits are associated with preeclampsia, but many of these genes have additional pleiotropic effects on cardiometabolic, endothelial, and placental function. Furthermore, several of the associated loci have no known connection with cardiovascular disease but instead harbor genes contributing to maintenance of successful pregnancy, with dysfunctions leading to preeclampsialike symptoms.

Deficit of homozygosity among 1.52 million individuals and genetic causes of recessive lethality.

Genotypes causing pregnancy loss and perinatal mortality are depleted among living individuals and are therefore difficult to find. To explore genetic causes of recessive lethality, we searched for sequence variants with deficit of homozygosity among 1.52 million individuals from six European populations. In this study, we identified 25 genes harboring protein-altering sequence variants with a strong deficit of homozygosity (10% or less of predicted homozygotes). Sequence variants in 12 of the genes cause Mendelian disease under a recessive mode of inheritance, two under a dominant mode, but variants in the remaining 11 have not been reported to cause disease. Sequence variants with a strong deficit of homozygosity are over-represented among genes essential for growth of human cell lines and genes orthologous to mouse genes known to affect viability. The function of these genes gives insight into the genetics of intrauterine lethality. We also identified 1077 genes with homozygous predicted loss-of-function genotypes not previously described, bringing the total set of genes completely knocked out in humans to 4785.

Sulfate sensitivity of aquatic organism in soft freshwaters explored by toxicity tests and species sensitivity distribution.

Elevated concentrations of sulfate in waterways are observed due to various anthropogenic activities. Elevated levels of sulfate can have harmful effects on aquatic life in freshwaters: sulfate can cause osmotic stress or specific ion toxicity in aquatic organisms, especially in soft waters where Ca2+ and Mg2+ concentrations are low. Formerly, chronic toxicity test data in soft water have been scarce. The chronic and acute sulfate toxicity tests conducted with aquatic organisms from 10 families across various trophic levels in this study multiplied the number of tests conducted in soft freshwater conditions and enabled derivation of the species sensitivity distribution (SSD) and sulfate hazardous concentrations for soft freshwaters. The cladoceran Daphnia longispina and freshwater snail Lymnaea stagnalis were the most sensitive to sulfate among the studied species. Harmful effects on the reproduction of D. longispina were observed at 49 mg SO4 /L while growth of L. stagnalis was inhibited at 217 mg SO4 /L. Most studied organisms tolerated high sulfate concentrations: the median of chronic effective concentrations (EC10 or LC10) was 1008 mg/L for all the species tested in this study. Based on the species sensitivity distribution of the studied species the hazardous concentration for 5 % of aquatic organism (HC5) in soft waters was 117-194 mg SO4/L. Different data set combinations were used to demonstrate the data variability in SSD-based HC5 estimates. The lowest values were produced from combining biotest results from the present study and earlier literature, while the highest values were calculated from the present study only. The derived chronic no-effect concentrations (PNEC) varied between 39 and 65 mg SO4/L.

Genome-wide screen of otosclerosis in population biobanks: 27 loci and shared associations with skeletal structure.

Otosclerosis is one of the most common causes of conductive hearing loss, affecting 0.3% of the population. It typically presents in adulthood and half of the patients have a positive family history. The pathophysiology of otosclerosis is poorly understood. A previous genome-wide association study (GWAS) identified a single association locus in an intronic region of RELN. Here, we report a meta-analysis of GWAS studies of otosclerosis in three population-based biobanks comprising 3504 cases and 861,198 controls. We identify 23 novel risk loci (p < 5 × 10-8) and report an association in RELN and three previously reported candidate gene or linkage regions (TGFB1, MEPE, and OTSC7). We demonstrate developmental stage-dependent immunostaining patterns of MEPE and RUNX2 in mouse otic capsules. In most association loci, the nearest protein-coding genes are implicated in bone remodelling, mineralization or severe skeletal disorders. We highlight multiple genes involved in transforming growth factor beta signalling for follow-up studies.

Inflammatory and infectious upper respiratory diseases associate with 41 genomic loci and type 2 inflammation.

Inflammatory and infectious upper respiratory diseases (ICD-10: J30-J39), such as diseases of the sinonasal tract, pharynx and larynx, are growing health problems yet their genomic similarity is not known. We analyze genome-wide association to eight upper respiratory diseases (61,195 cases) among 260,405 FinnGen participants, meta-analyzing diseases in four groups based on an underlying genetic correlation structure. Aiming to understand which genetic loci contribute to susceptibility to upper respiratory diseases in general and its subtypes, we detect 41 independent genome-wide significant loci, distinguishing impact on sinonasal or pharyngeal diseases, or both. Fine-mapping implicated non-synonymous variants in nine genes, including three linked to immune-related diseases. Phenome-wide analysis implicated asthma and atopic dermatitis at sinonasal disease loci, and inflammatory bowel diseases and other immune-mediated disorders at pharyngeal disease loci. Upper respiratory diseases also genetically correlated with autoimmune diseases such as rheumatoid arthritis, autoimmune hypothyroidism, and psoriasis. Finally, we associated separate gene pathways in sinonasal and pharyngeal diseases that both contribute to type 2 immunological reaction. We show shared heritability among upper respiratory diseases that extends to several immune-mediated diseases with diverse mechanisms, such as type 2 high inflammation.

Genome-wide association study of varicose veins identifies a protective missense variant in GJD3 enriched in the Finnish population.

Varicose veins is the most common manifestation of chronic venous disease that displays female-biased incidence. To identify protein-inactivating variants that could guide identification of drug target genes for varicose veins and genetic evidence for the disease prevalence difference between the sexes, we conducted a genome-wide association study of varicose veins in Finns using the FinnGen dataset with 17,027 cases and 190,028 controls. We identified 50 associated genetic loci (P < 5.0 × 10-8) of which 29 were novel including one near ERG with female-specificity (rs2836405-G, OR[95% CI] = 1.09[1.05-1.13], P = 3.1 × 10-8). These also include two X-chromosomal (ARHGAP6 and SRPX) and two autosomal novel loci (TGFB2 and GJD3) with protein-coding lead variants enriched above 56-fold in Finns over non-Finnish non-Estonian Europeans. A low-frequency missense variant in GJD3 (p.Pro59Thr) is exclusively associated with a lower risk for varicose veins (OR = 0.62 [0.55-0.70], P = 1.0 × 10-14) in a phenome-wide scan of the FinnGen data. The absence of observed pleiotropy and its membership of the connexin gene family underlines GJD3 as a potential connexin-modulating therapeutic strategy for varicose veins. Our results provide insights into varicose veins etiopathology and highlight the power of isolated populations, including Finns, to discover genetic variants that inform therapeutic development.

Genetic predictors of lifelong medication-use patterns in cardiometabolic diseases.

Little is known about the genetic determinants of medication use in preventing cardiometabolic diseases. Using the Finnish nationwide drug purchase registry with follow-up since 1995, we performed genome-wide association analyses of longitudinal patterns of medication use in hyperlipidemia, hypertension and type 2 diabetes in up to 193,933 individuals (55% women) in the FinnGen study. In meta-analyses of up to 567,671 individuals combining FinnGen with the Estonian Biobank and the UK Biobank, we discovered 333 independent loci (P < 5 × 10-9) associated with medication use. Fine-mapping revealed 494 95% credible sets associated with the total number of medication purchases, changes in medication combinations or treatment discontinuation, including 46 credible sets in 40 loci not associated with the underlying treatment targets. The polygenic risk scores (PRS) for cardiometabolic risk factors were strongly associated with the medication-use behavior. A medication-use enhanced multitrait PRS for coronary artery disease matched the performance of a risk factor-based multitrait coronary artery disease PRS in an independent sample (UK Biobank, n = 343,676). In summary, we demonstrate medication-based strategies for identifying cardiometabolic risk loci and provide genome-wide tools for preventing cardiovascular diseases.

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.

Efficient and accurate frailty model approach for genome-wide survival association analysis in large-scale biobanks.

With decades of electronic health records linked to genetic data, large biobanks provide unprecedented opportunities for systematically understanding the genetics of the natural history of complex diseases. Genome-wide survival association analysis can identify genetic variants associated with ages of onset, disease progression and lifespan. We propose an efficient and accurate frailty model approach for genome-wide survival association analysis of censored time-to-event (TTE) phenotypes by accounting for both population structure and relatedness. Our method utilizes state-of-the-art optimization strategies to reduce the computational cost. The saddlepoint approximation is used to allow for analysis of heavily censored phenotypes (>90%) and low frequency variants (down to minor allele count 20). We demonstrate the performance of our method through extensive simulation studies and analysis of five TTE phenotypes, including lifespan, with heavy censoring rates (90.9% to 99.8%) on ~400,000 UK Biobank participants with white British ancestry and ~180,000 individuals in FinnGen. We further analyzed 871 TTE phenotypes in the UK Biobank and presented the genome-wide scale phenome-wide association results with the PheWeb browser.