KIF15 missense variant is associated with the early onset of idiopathic pulmonary fibrosis.

BACKGROUND: Idiopathic pulmonary fibrosis (IPF) has an unknown aetiology and limited treatment options. A recent meta-analysis identified three novel causal variants in the TERT, SPDL1, and KIF15 genes. This observational study aimed to investigate whether the aforementioned variants cause clinical phenotypes in a well-characterised IPF cohort. METHODS: The study consisted of 138 patients with IPF who were diagnosed and treated at the Helsinki University Hospital and genotyped in the FinnGen FinnIPF study. Data on > 25 clinical parameters were collected by two pulmonologists who were blinded to the genetic data for patients with TERT loss of function and missense variants, SPDL1 and KIF15 missense variants, and a MUC5B variant commonly present in patients with IPF, or no variants were separately analysed. RESULTS: The KIF15 missense variant is associated with the early onset of the disease, leading to progression to early-age transplantation or death. In patients with the KIF15 variant, the median age at diagnosis was 54.0 years (36.5-69.5 years) compared with 72.0 years (65.8-75.3 years) in the other patients (P = 0.023). The proportion of KIF15 variant carriers was 9- or 3.6-fold higher in patients aged < 55 or 65 years, respectively. The variants for TERT and MUC5B had similar effects on the patient's clinical course, as previously described. No distinct phenotypes were observed in patients with the SPDL1 variant. CONCLUSIONS: Our study indicated the potential of KIF15 to be used in the genetic diagnostics of IPF. Further studies are needed to elucidate the biological mechanisms of KIF15 in IPF.

Geographic Variation and Bias in the Polygenic Scores of Complex Diseases and Traits in Finland.

Polygenic scores (PSs) are becoming a useful tool to identify individuals with high genetic risk for complex diseases, and several projects are currently testing their utility for translational applications. It is also tempting to use PSs to assess whether genetic variation can explain a part of the geographic distribution of a phenotype. However, it is not well known how the population genetic properties of the training and target samples affect the geographic distribution of PSs. Here, we evaluate geographic differences, and related biases, of PSs in Finland in a geographically well-defined sample of 2,376 individuals from the National FINRISK study. First, we detect geographic differences in PSs for coronary artery disease (CAD), rheumatoid arthritis, schizophrenia, waist-hip ratio (WHR), body-mass index (BMI), and height, but not for Crohn disease or ulcerative colitis. Second, we use height as a model trait to thoroughly assess the possible population genetic biases in PSs and apply similar approaches to the other phenotypes. Most importantly, we detect suspiciously large accumulations of geographic differences for CAD, WHR, BMI, and height, suggesting bias arising from the population's genetic structure rather than from a direct genotype-phenotype association. This work demonstrates how sensitive the geographic patterns of current PSs are for small biases even within relatively homogeneous populations and provides simple tools to identify such biases. A thorough understanding of the effects of population genetic structure on PSs is essential for translational applications of PSs.

Correction: Insights into the genetic epidemiology of Crohn's and rare diseases in the Ashkenazi Jewish population.

[This corrects the article DOI: 10.1371/journal.pgen.1007329.].

Lifetime risk of rheumatoid arthritis-associated interstitial lung disease in MUC5B mutation carriers.

OBJECTIVES: To estimate lifetime risk of developing rheumatoid arthritis-associated interstitial lung disease (RA-ILD) with respect to the strongest known risk factor for pulmonary fibrosis, a MUC5B promoter variant. METHODS: FinnGen is a collection of epidemiological cohorts and hospital biobank samples, integrating genetic data with up to 50 years of follow-up within nationwide registries in Finland. Patients with RA and ILD were identified from the Finnish national hospital discharge, medication reimbursement and cause-of-death registries. We estimated lifetime risks of ILD by age 80 with respect to the common variant rs35705950, a MUC5B promoter variant. RESULTS: Out of 293 972 individuals, 1965 (0.7%) developed ILD by age 80. Among all individuals in the dataset, MUC5B increased the risk of ILD with a HR of 2.44 (95% CI: 2.22 to 2.68). Out of 6869 patients diagnosed with RA, 247 (3.6%) developed ILD. In patients with RA, MUC5B was a strong risk factor of ILD with a HR similar to the full dataset (HR: 2.27, 95% CI: 1.75 to 2.95). In patients with RA, lifetime risks of ILD were 16.8% (95% CI: 13.1% to 20.2%) for MUC5B carriers and 6.1% (95% CI: 5.0% to 7.2%) for MUC5B non-carriers. The difference between risks started to emerge at age 65, with a higher risk among men. CONCLUSION: Our findings provide estimates of lifetime risk of RA-ILD based on MUC5B mutation carrier status, demonstrating the potential of genomics for risk stratification of RA-ILD.

Exome sequencing in patient-parent trios suggests new candidate genes for early-onset primary sclerosing cholangitis.

BACKGROUND & AIMS: Primary sclerosing cholangitis (PSC) is a rare bile duct disease strongly associated with inflammatory bowel disease (IBD). Whole-exome sequencing (WES) has contributed to understanding the molecular basis of very early-onset IBD, but rare protein-altering genetic variants have not been identified for early-onset PSC. We performed WES in patients diagnosed with PSC ≤ 12 years to investigate the contribution of rare genetic variants to early-onset PSC. METHODS: In this multicentre study, WES was performed on 87 DNA samples from 29 patient-parent trios with early-onset PSC. We selected rare (minor allele frequency < 2%) coding and splice-site variants that matched recessive (homozygous and compound heterozygous variants) and dominant (de novo) inheritance in the index patients. Variant pathogenicity was predicted by an in-house developed algorithm (GAVIN), and PSC-relevant variants were selected using gene expression data and gene function. RESULTS: In 22 of 29 trios we identified at least 1 possibly pathogenic variant. We prioritized 36 genes, harbouring a total of 54 variants with predicted pathogenic effects. In 18 genes, we identified 36 compound heterozygous variants, whereas in the other 18 genes we identified 18 de novo variants. Twelve of 36 candidate risk genes are known to play a role in transmembrane transport, adaptive and innate immunity, and epithelial barrier function. CONCLUSIONS: The 36 candidate genes for early-onset PSC need further verification in other patient cohorts and evaluation of gene function before a causal role can be attributed to its variants.

Insights into the genetic epidemiology of Crohn's and rare diseases in the Ashkenazi Jewish population.

As part of a broader collaborative network of exome sequencing studies, we developed a jointly called data set of 5,685 Ashkenazi Jewish exomes. We make publicly available a resource of site and allele frequencies, which should serve as a reference for medical genetics in the Ashkenazim (hosted in part at https://ibd.broadinstitute.org, also available in gnomAD at http://gnomad.broadinstitute.org). We estimate that 34% of protein-coding alleles present in the Ashkenazi Jewish population at frequencies greater than 0.2% are significantly more frequent (mean 15-fold) than their maximum frequency observed in other reference populations. Arising via a well-described founder effect approximately 30 generations ago, this catalog of enriched alleles can contribute to differences in genetic risk and overall prevalence of diseases between populations. As validation we document 148 AJ enriched protein-altering alleles that overlap with "pathogenic" ClinVar alleles (table available at https://github.com/macarthur-lab/clinvar/blob/master/output/clinvar.tsv), including those that account for 10-100 fold differences in prevalence between AJ and non-AJ populations of some rare diseases, especially recessive conditions, including Gaucher disease (GBA, p.Asn409Ser, 8-fold enrichment); Canavan disease (ASPA, p.Glu285Ala, 12-fold enrichment); and Tay-Sachs disease (HEXA, c.1421+1G>C, 27-fold enrichment; p.Tyr427IlefsTer5, 12-fold enrichment). We next sought to use this catalog, of well-established relevance to Mendelian disease, to explore Crohn's disease, a common disease with an estimated two to four-fold excess prevalence in AJ. We specifically attempt to evaluate whether strong acting rare alleles, particularly protein-truncating or otherwise large effect-size alleles, enriched by the same founder-effect, contribute excess genetic risk to Crohn's disease in AJ, and find that ten rare genetic risk factors in NOD2 and LRRK2 are enriched in AJ (p < 0.005), including several novel contributing alleles, show evidence of association to CD. Independently, we find that genomewide common variant risk defined by GWAS shows a strong difference between AJ and non-AJ European control population samples (0.97 s.d. higher, p<10-16). Taken together, the results suggest coordinated selection in AJ population for higher CD risk alleles in general. The results and approach illustrate the value of exome sequencing data in case-control studies along with reference data sets like ExAC (sites VCF available via FTP at ftp.broadinstitute.org/pub/ExAC_release/release0.3/) to pinpoint genetic variation that contributes to variable disease predisposition across populations.

Association of Genetic Variants in NUDT15 With Thiopurine-Induced Myelosuppression in Patients With Inflammatory Bowel Disease.

Importance: Use of thiopurines may be limited by myelosuppression. TPMT pharmacogenetic testing identifies only 25% of at-risk patients of European ancestry. Among patients of East Asian ancestry, NUDT15 variants are associated with thiopurine-induced myelosuppression (TIM). Objective: To identify genetic variants associated with TIM among patients of European ancestry with inflammatory bowel disease (IBD). Design, Setting, and Participants: Case-control study of 491 patients affected by TIM and 679 thiopurine-tolerant unaffected patients who were recruited from 89 international sites between March 2012 and November 2015. Genome-wide association studies (GWAS) and exome-wide association studies (EWAS) were conducted in patients of European ancestry. The replication cohort comprised 73 patients affected by TIM and 840 thiopurine-tolerant unaffected patients. Exposures: Genetic variants associated with TIM. Main Outcomes and Measures: Thiopurine-induced myelosuppression, defined as a decline in absolute white blood cell count to 2.5 × 109/L or less or a decline in absolute neutrophil cell count to 1.0 × 109/L or less leading to a dose reduction or drug withdrawal. Results: Among 1077 patients (398 affected and 679 unaffected; median age at IBD diagnosis, 31.0 years [interquartile range, 21.2 to 44.1 years]; 540 [50%] women; 602 [56%] diagnosed as having Crohn disease), 919 (311 affected and 608 unaffected) were included in the GWAS analysis and 961 (328 affected and 633 unaffected) in the EWAS analysis. The GWAS analysis confirmed association of TPMT (chromosome 6, rs11969064) with TIM (30.5% [95/311] affected vs 16.4% [100/608] unaffected patients; odds ratio [OR], 2.3 [95% CI, 1.7 to 3.1], P = 5.2 × 10-9). The EWAS analysis demonstrated an association with an in-frame deletion in NUDT15 (chromosome 13, rs746071566) and TIM (5.8% [19/328] affected vs 0.2% [1/633] unaffected patients; OR, 38.2 [95% CI, 5.1 to 286.1], P = 1.3 × 10-8), which was replicated in a different cohort (2.7% [2/73] affected vs 0.2% [2/840] unaffected patients; OR, 11.8 [95% CI, 1.6 to 85.0], P = .03). Carriage of any of 3 coding NUDT15 variants was associated with an increased risk (OR, 27.3 [95% CI, 9.3 to 116.7], P = 1.1 × 10-7) of TIM, independent of TPMT genotype and thiopurine dose. Conclusions and Relevance: Among patients of European ancestry with IBD, variants in NUDT15 were associated with increased risk of TIM. These findings suggest that NUDT15 genotyping may be considered prior to initiation of thiopurine therapy; however, further study including additional validation in independent cohorts is required.

Individual FEV1 Trajectories Can Be Identified from a COPD Cohort.

OBJECTIVE: We aim to make use of clinical spirometry data in order to identify individual COPD-patients with divergent trajectories of lung function over time. STUDY DESIGN AND SETTING: Hospital-based COPD cohort (N = 607) was followed on average 4.6 years. Each patient had a mean of 8.4 spirometries available. We used a Hierarchical Bayesian Model (HBM) to identify the individuals presenting constant trends in lung function. RESULTS: At a probability level of 95%, one third of the patients (180/607) presented rapidly declining FEV1 (mean -78 ml/year, 95% CI -73 to -83 ml) compared to that in the rest of the patients (mean -26 ml/year, 95% CI -23 to -29 ml, p ≤ 2.2 × 10(-16)). Constant improvement of FEV1 was very rare. The rapid decliners more frequently suffered from exacerbations measured by various outcome markers. CONCLUSION: Clinical data of unique patients can be utilized to identify diverging trajectories of FEV1 with a high probability. Frequent exacerbations were more prevalent in FEV1-decliners than in the rest of the patients. The result confirmed previously reported association between FEV1 decline and exacerbation rate and further suggested that in clinical practice HBM could improve the identification of high-risk individuals at early stages of the disease.

Differences in plasma and sputum biomarkers between COPD and COPD-asthma overlap.

The pathophysiological features of chronic obstructive pulmonary disease (COPD)-asthma overlap are poorly understood and there has been no study of plasma or sputum biomarkers in overlap patients. In order to clarify the similarity and differences between overlap and COPD or asthma, we have investigated four potential biomarkers of COPD: surfactant protein A (SP-A), soluble receptor for advanced glycation end-products (sRAGE), myeloperoxidase (MPO) and neutrophil gelatinase-associated lipocalin (NGAL). SP-A and sRAGE are pneumocyte-derived markers. MPO and NGAL are neutrophil-derived molecules, but NGAL can also be expressed by respiratory epithelial cells. Plasma levels of SP-A and sRAGE and induced sputum levels of MPO and NGAL were measured by enzyme immunoassay/ELISA in 134 subjects: nonsmokers (n=26), smokers (n=23), asthma (n=32), COPD (n=39) and COPD-asthma overlap patients (n=14). In patients with COPD-asthma overlap, sputum MPO and plasma SP-A were significantly elevated whereas plasma sRAGE levels were reduced compared with asthma patients. Only sputum NGAL was significantly elevated in COPD-asthma overlap compared with COPD (p=0.00016) and could be used to differentiate patients with overlap from those with COPD. Increased induced sputum levels of NGAL might be a characteristic feature of overlap, suggesting enhanced neutrophilic airway inflammation and/or airway epithelial injury in COPD-asthma overlap.

Co-morbidities are the key nominators of the health related quality of life in mild and moderate COPD.

BACKGROUND: Co-morbidities are common in chronic obstructive pulmonary disease (COPD). We assessed the contribution of common co-morbidities on health related quality of life (HRQoL) among COPD patients. METHODS: Using both generic (15D) and respiratory-specific (AQ20) instruments, HRQoL was assessed in a hospital based COPD population (N = 739, 64% males, mean age 64 years, SD 7 years) in this observational study with inferential analysis. The prevalence of their co-morbidities was compared with those of 5000 population controls. The patients represented all severity stages of COPD and the patterns of common concomitant disorders differed between patients. RESULTS: Co-morbidities such as psychiatric conditions, alcohol abuse, cardiovascular diseases, and diabetes were more common among COPD patients than in age and gender matched controls. Psychiatric conditions and alcohol abuse were the strongest determinants of HRQoL in COPD and could be detected by both 15D (Odds Ratio 4.7 and 2.3 respectively) and AQ20 (OR 2.0 and 3.0) instruments. Compared to respiratory specific AQ20, generic 15D was more sensitive to the effects of comorbidities while AQ20 was slightly more sensitive for the low FEV1. FEV1 was a strong determinant of HRQoL only at more severe stages of disease (FEV1 < 40% of predicted). Poor HRQoL also predicted death during the next five years. CONCLUSIONS: The results suggest that co-morbidities may impair HRQoL at an early stage of the disease, while bronchial obstruction becomes a significant determinant of HRQoL only in severe COPD.

Genetic Analysis of Obstructive Sleep Apnea and Its Relationship with Severe COVID-19.

Rationale: Although patients with obstructive sleep apnea (OSA) have a higher risk for coronavirus disease (COVID-19) hospitalization, the causal relationship has remained unexplored. Objectives: To understand the causal relationship between OSA and COVID-19 by leveraging data from vaccination and electronic health records, genetic risk factors from genome-wide association studies, and Mendelian randomization. Methods: We elucidated genetic risk factors for OSA using FinnGen (total N = 377,277), performing genome-wide association. We used the associated variants as instruments for univariate and multivariate Mendelian randomization (MR) analyses and computed absolute risk reduction against COVID-19 hospitalization with or without vaccination. Results: We identified nine novel loci for OSA and replicated our findings in the Million Veteran Program. Furthermore, MR analysis showed that OSA was a causal risk factor for severe COVID-19 (P = 9.41 × 10-4). Probabilistic modeling showed that the strongest genetic risk factor for OSA at the FTO locus reflected a signal of higher body mass index (BMI), whereas BMI-independent association was seen with the earlier reported SLC9A4 locus and a MECOM locus, which is a transcriptional regulator with 210-fold enrichment in the Finnish population. Similarly, multivariate MR analysis showed that the causality for severe COVID-19 was driven by BMI (multivariate MR P = 5.97 × 10-6, ß = 0.47). Finally, vaccination reduced the risk for COVID-19 hospitalization more in the patients with OSA than in the non-OSA controls, with respective absolute risk reductions of 13.3% versus 6.3%. Conclusions: Our analysis identified novel genetic risk factors for OSA and showed that OSA is a causal risk factor for severe COVID-19. The effect is predominantly explained by higher BMI and suggests BMI-dependent effects at the level of individual variants and at the level of comorbid causality.

Face masks to prevent transmission of respiratory infections: Systematic review and meta-analysis of randomized controlled trials on face mask use.

OBJECTIVES: To examine the use of face mask intervention in mitigating the risk of spreading respiratory infections and whether the effect of face mask intervention differs in different exposure settings and age groups. DESIGN: Systematic review and meta-analysis. We evaluated the risk of bias using the Cochrane Risk of Bias 2 tool (ROB2). DATA SOURCES: We searched PubMed, Embase, Cochrane Central Register of Controlled Trials, and Web of Science were searched for randomized controlled trials investigating the effect of face masks on respiratory infections published between 1981 and February 9, 2022. We followed the PRISMA 2020 guidelines. ELIGIBILITY CRITERIA FOR SELECTING STUDIES: We included randomized controlled trials investigating the use of face mask intervention in mitigating the risk of spreading respiratory infections across different exposure settings. RESULTS: We identified 2,400 articles for screening. 18 articles passed the inclusion criteria for both evidence synthesis and meta-analysis. There were N = 189,145 individuals in the face mask intervention arm and N = 173,536 in the control arm, and the follow-up times ranged from 4 days to 19 months. Our results showed between-study heterogeneity (p < 0.0001). While there was no statistically significant association over all studies when the covariate unadjusted intervention effect estimates were used (RR = 0.977 [0.858-1.113], p = 0.728), our subgroup analyses revealed that a face mask intervention reduced respiratory infections in the adult subgroup (RR = 0.8795 [0.7861-0.9839], p = 0.0249) and in a community setting (RR = 0.890 [0.812-0.975], p = 0.0125). Furthermore, our leave-one-out analysis found that one study biased the results towards a null effect. Consequently, when using covariate adjusted odds ratio estimates to have a more precise effect estimates of the intervention effect to account for differences at the baseline, the results showed that a face mask intervention did reduce respiratory infections when the biasing study was excluded from the analysis (OR = 0.8892 [0.8061-0.9810], p = 0.0192). CONCLUSION: Our findings support the use of face masks particularly in a community setting and for adults. We also observed substantial between-study heterogeneity and varying adherence to protocol. Notably, many studies were subject to contamination bias thus affecting the efficacy of the intervention, that is when also some controls used masks or when the intervention group did not comply with mask use leading to a downward biased effect of treatment receipt and efficacy. TRIAL REGISTRATION: PROSPERO registration number CRD42020205523.

Leveraging global multi-ancestry meta-analysis in the study of idiopathic pulmonary fibrosis genetics.

The research of rare and devastating orphan diseases, such as idiopathic pulmonary fibrosis (IPF) has been limited by the rarity of the disease itself. The prognosis is poor-the prevalence of IPF is only approximately four times the incidence, limiting the recruitment of patients to trials and studies of the underlying biology. Global biobanking efforts can dramatically alter the future of IPF research. We describe a large-scale meta-analysis of IPF, with 8,492 patients and 1,355,819 population controls from 13 biobanks around the globe. Finally, we combine this meta-analysis with the largest available meta-analysis of IPF, reaching 11,160 patients and 1,364,410 population controls. We identify seven novel genome-wide significant loci, only one of which would have been identified if the analysis had been limited to European ancestry individuals. We observe notable pleiotropy across IPF susceptibility and severe COVID-19 infection and note an unexplained sex-heterogeneity effect at the strongest IPF locus MUC5B.

Genetic Associations and Architecture of Asthma-COPD Overlap.

BACKGROUND: Some people have characteristics of both asthma and COPD (asthma-COPD overlap), and evidence suggests they experience worse outcomes than those with either condition alone. RESEARCH QUESTION: What is the genetic architecture of asthma-COPD overlap, and do the determinants of risk for asthma-COPD overlap differ from those for COPD or asthma? STUDY DESIGN AND METHODS: We conducted a genome-wide association study in 8,068 asthma-COPD overlap case subjects and 40,360 control subjects without asthma or COPD of European ancestry in UK Biobank (stage 1). We followed up promising signals (P < 5 × 10-6) that remained associated in analyses comparing (1) asthma-COPD overlap vs asthma-only control subjects, and (2) asthma-COPD overlap vs COPD-only control subjects. These variants were analyzed in 12 independent cohorts (stage 2). RESULTS: We selected 31 independent variants for further investigation in stage 2, and discovered eight novel signals (P < 5 × 10-8) for asthma-COPD overlap (meta-analysis of stage 1 and 2 studies). These signals suggest a spectrum of shared genetic influences, some predominantly influencing asthma (FAM105A, GLB1, PHB, TSLP), others predominantly influencing fixed airflow obstruction (IL17RD, C5orf56, HLA-DQB1). One intergenic signal on chromosome 5 had not been previously associated with asthma, COPD, or lung function. Subgroup analyses suggested that associations at these eight signals were not driven by smoking or age at asthma diagnosis, and in phenome-wide scans, eosinophil counts, atopy, and asthma traits were prominent. INTERPRETATION: We identified eight signals for asthma-COPD overlap, which may represent loci that predispose to type 2 inflammation, and serious long-term consequences of asthma.

Sleep apnoea is a risk factor for severe COVID-19.

BACKGROUND: Obstructive sleep apnoea (OSA) is associated with higher body mass index (BMI), diabetes, older age and male gender, which are all risk factors for severe COVID-19.We aimed to study if OSA is an independent risk factor for COVID-19 infection or for severe COVID-19. METHODS: OSA diagnosis and COVID-19 infection were extracted from the hospital discharge, causes of death and infectious diseases registries in individuals who participated in the FinnGen study (n=260 405). Severe COVID-19 was defined as COVID-19 requiring hospitalisation. Multivariate logistic regression model was used to examine association. Comorbidities for either COVID-19 or OSA were selected as covariates. We performed a meta-analysis with previous studies. RESULTS: We identified 445 individuals with COVID-19, and 38 (8.5%) of them with OSA of whom 19 out of 91 (20.9%) were hospitalised. OSA associated with COVID-19 hospitalisation independent from age, sex, BMI and comorbidities (p-unadjusted=5.13×10-5, OR-adjusted=2.93 (95% CI 1.02 to 8.39), p-adjusted=0.045). OSA was not associated with the risk of contracting COVID-19 (p=0.25). A meta-analysis of OSA and severe COVID-19 showed association across 15 835 COVID-19 positive controls, and n=1294 patients with OSA with severe COVID-19 (OR=2.37 (95% 1.14 to 4.95), p=0.021). CONCLUSION: Risk for contracting COVID-19 was the same for patients with OSA and those without OSA. In contrast, among COVID-19 positive patients, OSA was associated with higher risk for hospitalisation. Our findings are in line with earlier works and suggest OSA as an independent risk factor for severe COVID-19.

An expanded analysis framework for multivariate GWAS connects inflammatory biomarkers to functional variants and disease.

Multivariate methods are known to increase the statistical power to detect associations in the case of shared genetic basis between phenotypes. They have, however, lacked essential analytic tools to follow-up and understand the biology underlying these associations. We developed a novel computational workflow for multivariate GWAS follow-up analyses, including fine-mapping and identification of the subset of traits driving associations (driver traits). Many follow-up tools require univariate regression coefficients which are lacking from multivariate results. Our method overcomes this problem by using Canonical Correlation Analysis to turn each multivariate association into its optimal univariate Linear Combination Phenotype (LCP). This enables an LCP-GWAS, which in turn generates the statistics required for follow-up analyses. We implemented our method on 12 highly correlated inflammatory biomarkers in a Finnish population-based study. Altogether, we identified 11 associations, four of which (F5, ABO, C1orf140 and PDGFRB) were not detected by biomarker-specific analyses. Fine-mapping identified 19 signals within the 11 loci and driver trait analysis determined the traits contributing to the associations. A phenome-wide association study on the 19 representative variants from the signals in 176,899 individuals from the FinnGen study revealed 53 disease associations (p < 1 × 10-4). Several reported pQTLs in the 11 loci provided orthogonal evidence for the biologically relevant functions of the representative variants. Our novel multivariate analysis workflow provides a powerful addition to standard univariate GWAS analyses by enabling multivariate GWAS follow-up and thus promoting the advancement of powerful multivariate methods in genomics.

Polygenic and clinical risk scores and their impact on age at onset and prediction of cardiometabolic diseases and common cancers.

Polygenic risk scores (PRSs) have shown promise in predicting susceptibility to common diseases1-3. We estimated their added value in clinical risk prediction of five common diseases, using large-scale biobank data (FinnGen; n = 135,300) and the FINRISK study with clinical risk factors to test genome-wide PRSs for coronary heart disease, type 2 diabetes, atrial fibrillation, breast cancer and prostate cancer. We evaluated the lifetime risk at different PRS levels, and the impact on disease onset and on prediction together with clinical risk scores. Compared to having an average PRS, having a high PRS contributed 21% to 38% higher lifetime risk, and 4 to 9 years earlier disease onset. PRSs improved model discrimination over age and sex in type 2 diabetes, atrial fibrillation, breast cancer and prostate cancer, and over clinical risk in type 2 diabetes, breast cancer and prostate cancer. In all diseases, PRSs improved reclassification over clinical thresholds, with the largest net reclassification improvements for early-onset coronary heart disease, atrial fibrillation and prostate cancer. This study provides evidence for the additional value of PRSs in clinical disease prediction. The practical applications of polygenic risk information for stratified screening or for guiding lifestyle and medical interventions in the clinical setting remain to be defined in further studies.

Genomic prediction of alcohol-related morbidity and mortality.

While polygenic risk scores (PRS) have been shown to predict many diseases and risk factors, the potential of genomic prediction in harm caused by alcohol use has not yet been extensively studied. Here, we built a novel polygenic risk score of 1.1 million variants for alcohol consumption and studied its predictive capacity in 96,499 participants from the FinnGen study and 39,695 participants from prospective cohorts with detailed baseline data and up to 25 years of follow-up time. A 1 SD increase in the PRS was associated with 11.2 g (=0.93 drinks) higher weekly alcohol consumption (CI = 9.85-12.58 g, p = 2.3 × 10-58). The PRS was associated with alcohol-related morbidity (4785 incident events) and the risk estimate between the highest and lowest quintiles of the PRS was 1.83 (95% CI = 1.66-2.01, p = 1.6 × 10-36). When adjusted for self-reported alcohol consumption, education, marital status, and gamma-glutamyl transferase blood levels in 28,639 participants with comprehensive baseline data from prospective cohorts, the risk estimate between the highest and lowest quintiles of the PRS was 1.58 (CI = 1.26-1.99, p = 8.2 × 10-5). The PRS was also associated with all-cause mortality with a risk estimate of 1.33 between the highest and lowest quintiles (CI = 1.20-1.47, p = 4.5 × 10-8) in the adjusted model. In conclusion, the PRS for alcohol consumption independently associates for both alcohol-related morbidity and all-cause mortality. Together, these findings underline the importance of heritable factors in alcohol-related health burden while highlighting how measured genetic risk for an important behavioral risk factor can be used to predict related health outcomes.

Genetic architecture of human plasma lipidome and its link to cardiovascular disease.

Understanding genetic architecture of plasma lipidome could provide better insights into lipid metabolism and its link to cardiovascular diseases (CVDs). Here, we perform genome-wide association analyses of 141 lipid species (n = 2,181 individuals), followed by phenome-wide scans with 25 CVD related phenotypes (n = 511,700 individuals). We identify 35 lipid-species-associated loci (P <5 ×10-8), 10 of which associate with CVD risk including five new loci-COL5A1, GLTPD2, SPTLC3, MBOAT7 and GALNT16 (false discovery rate<0.05). We identify loci for lipid species that are shown to predict CVD e.g., SPTLC3 for CER(d18:1/24:1). We show that lipoprotein lipase (LPL) may more efficiently hydrolyze medium length triacylglycerides (TAGs) than others. Polyunsaturated lipids have highest heritability and genetic correlations, suggesting considerable genetic regulation at fatty acids levels. We find low genetic correlations between traditional lipids and lipid species. Our results show that lipidomic profiles capture information beyond traditional lipids and identify genetic variants modifying lipid levels and risk of CVD.

Characterization of sputum biomarkers for asthma-COPD overlap syndrome.

Asthma-COPD overlap syndrome (ACOS) is a commonly encountered chronic airway disease. However, ACOS is still a consensus-based clinical phenotype and the underlying inflammatory mechanisms are inadequately characterized. To clarify the inflammatory mediatypical for ACOS, five biomarkers, namely interleukin (IL)-13, myeloperoxidase (MPO), neutrophil gelatinase-associated lipocalin (NGAL), chitinase-like protein (YKL-40), and IL-6, were selected. This study hypothesized that sputum biomarkers relevant for airway inflammation in asthma (IL-13), COPD (MPO, NGAL), or in both asthma and COPD (YKL-40, IL-6) could be used to differentiate ACOS from COPD and asthma. The aim of this study was to characterize the inflammatory profile and improve the recognition of ACOS. Induced sputum levels of IL-13, MPO, NGAL, YKL-40, and IL-6 were measured by enzyme-linked immunosorbent assay/Luminex assay in a Finnish discovery cohort (n=90) of nonsmokers, smokers, and patients with asthma, COPD, and ACOS and validated in a Japanese cohort (n=135). The classification accuracy of potential biomarkers was compared with area under the receiver operating characteristic curves. Only sputum NGAL levels could differentiate ACOS from asthma (P<0.001 and P<0.001) and COPD (P<0.05 and P=0.002) in the discovery and replication cohorts, respectively. Sputum NGAL levels were independently correlated with the percentage of pre-bronchodilator forced expiratory volume in 1 second predicted in multivariate analysis in the discovery and replication cohorts (P=0.001 and P=0.002, respectively). In conclusion, sputum biomarkers reflecting both airway inflammation and remodeling of the tissue show potential in differentiation between asthma, COPD, and ACOS.