ABSTRACT
SUMMARY: The genetic architecture of complex traits can be influenced by both many common regulatory variants with small effect sizes and rare deleterious variants in coding regions with larger effect sizes. However, the two kinds of genetic contributions are typically analyzed independently. Here, we present GenRisk, a python package for the computation and the integration of gene scores based on the burden of rare deleterious variants and common-variants-based polygenic risk scores. The derived scores can be analyzed within GenRisk to perform association tests or to derive phenotype prediction models by testing multiple classification and regression approaches. GenRisk is compatible with VCF input file formats. AVAILABILITY AND IMPLEMENTATION: GenRisk is an open source publicly available python package that can be downloaded or installed from Github (https://github.com/AldisiRana/GenRisk). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
Subject(s)
Multifactorial Inheritance , Software , Phenotype , Open Reading Frames , Risk FactorsABSTRACT
PURPOSE: We aimed to investigate to what extent polygenic risk scores (PRS), rare pathogenic germline variants (PVs), and family history jointly influence breast cancer and prostate cancer risk. METHODS: A total of 200,643 individuals from the UK Biobank were categorized as follows: (1) heterozygotes or nonheterozygotes for PVs in moderate to high-risk cancer genes, (2) PRS strata, and (3) with or without a family history of cancer. Multivariable logistic regression and Cox proportional hazards models were used to compute the odds ratio across groups and the cumulative incidence through life. RESULTS: Cumulative incidence by age 70 years among the nonheterozygotes across PRS strata ranged from 9% to 32% and from 9% to 35% for breast cancer and prostate cancer, respectively. Among the PV heterozygotes it ranged from 20% to 48% in moderate-risk genes and from 51% to 74% in high-risk genes for breast cancer, and it ranged from 30% to 59% in prostate cancer risk genes. Family history was always associated with an increased cancer odds ratio. CONCLUSION: PRS alone provides a meaningful risk gradient leading to a cancer risk stratification comparable to PVs in moderate risk genes, whereas acts as a risk modifier when considering high-risk genes. Including family history along with PV and PRS further improves cancer risk stratification.
Subject(s)
Breast Neoplasms , Prostatic Neoplasms , Aged , Breast Neoplasms/epidemiology , Breast Neoplasms/genetics , Genetic Predisposition to Disease , Germ Cells , Humans , Male , Multifactorial Inheritance/genetics , Prostatic Neoplasms/epidemiology , Prostatic Neoplasms/genetics , Risk FactorsABSTRACT
Male-pattern hair loss (MPHL) is common and highly heritable. While genome-wide association studies (GWAS) have generated insights into the contribution of common variants to MPHL etiology, the relevance of rare variants remains unclear. To determine the contribution of rare variants to MPHL etiology, we perform gene-based and single-variant analyses in exome-sequencing data from 72,469 male UK Biobank participants. While our population-level risk prediction suggests that rare variants make only a minor contribution to general MPHL risk, our rare variant collapsing tests identified a total of five significant gene associations. These findings provide additional evidence for previously implicated genes (EDA2R, WNT10A) and highlight novel risk genes at and beyond GWAS loci (HEPH, CEPT1, EIF3F). Furthermore, MPHL-associated genes are enriched for genes considered causal for monogenic trichoses. Together, our findings broaden the MPHL-associated allelic spectrum and provide insights into MPHL pathobiology and a shared basis with monogenic hair loss disorders.
Subject(s)
Biological Specimen Banks , Exome , Humans , Male , Exome/genetics , Genome-Wide Association Study , Alopecia/genetics , United KingdomABSTRACT
BACKGROUND: A relevant part of the genetic architecture of complex traits is still unknown; despite the discovery of many disease-associated common variants. Polygenic risk score (PRS) models are based on the evaluation of the additive effects attributable to common variants and have been successfully implemented to assess the genetic susceptibility for many phenotypes. In contrast, burden tests are often used to identify an enrichment of rare deleterious variants in specific genes. Both kinds of genetic contributions are typically analyzed independently. Many studies suggest that complex phenotypes are influenced by both low effect common variants and high effect rare deleterious variants. The aim of this paper is to integrate the effect of both common and rare functional variants for a more comprehensive genetic risk modeling. METHODS: We developed a framework combining gene-based scores based on the enrichment of rare functionally relevant variants with genome-wide PRS based on common variants for association analysis and prediction models. We applied our framework on UK Biobank dataset with genotyping and exome data and considered 28 blood biomarkers levels as target phenotypes. For each biomarker, an association analysis was performed on full cohort using gene-based scores (GBS). The cohort was then split into 3 subsets for PRS construction and feature selection, predictive model training, and independent evaluation, respectively. Prediction models were generated including either PRS, GBS or both (combined). RESULTS: Association analyses of the cohort were able to detect significant genes that were previously known to be associated with different biomarkers. Interestingly, the analyses also revealed heterogeneous effect sizes and directionality highlighting the complexity of the blood biomarkers regulation. However, the combined models for many biomarkers show little or no improvement in prediction accuracy compared to the PRS models. CONCLUSION: This study shows that rare variants play an important role in the genetic architecture of complex multifactorial traits such as blood biomarkers. However, while rare deleterious variants play a strong role at an individual level, our results indicate that classical common variant based PRS might be more informative to predict the genetic susceptibility at the population level.
Subject(s)
Exome , Genetic Predisposition to Disease , Humans , Genetic Predisposition to Disease/genetics , Biomarkers , Phenotype , Multifactorial Inheritance/geneticsABSTRACT
Polygenic risk score (PRS) predictions often show bias toward the population of available genome-wide association studies (GWASs), which is typically of European ancestry. This study aimed to assess the performance differences of ancestry-specific PRS and test the implementation of multi-ancestry PRS to enhance the generalizability of low-density lipoprotein (LDL) cholesterol predictions in the East Asian (EAS) population. In this study, we computed ancestry-specific and multi-ancestry PRSs for LDL using data obtained from the Global Lipid Genetics Consortium, while accounting for population-specific linkage disequilibrium patterns using the PRS-CSx method in the United Kingdom Biobank dataset (UKB, n = 423,596) and Taiwan Biobank dataset (TWB, n = 68,978). Population-specific PRSs were able to predict LDL levels better within the target population, whereas multi-ancestry PRSs were more generalizable. In the TWB dataset, covariate-adjusted R 2 values were 9.3% for ancestry-specific PRS, 6.7% for multi-ancestry PRS, and 4.5% for European-specific PRS. Similar trends (8.6%, 7.8%, and 6.2%) were observed in the smaller EAS population of the UKB (n = 1,480). Consistent with R 2 values, PRS stratification in EAS regions (TWB) effectively captured a heterogenous variability in LDL blood cholesterol levels across PRS strata. The mean difference in LDL levels between the lowest and highest EAS-specific PRS (EAS_PRS) deciles was 0.82, compared to 0.59 for European-specific PRS (EUR_PRS) and 0.76 for multi-ancestry PRS. Notably, the mean LDL values in the top decile of multi-ancestry PRS were comparable to those of EAS_PRS (3.543 vs. 3.541, p = 0.86). Our analysis of the PRS prediction model for LDL cholesterol further supports the issue of PRS generalizability across populations. Our targeted analysis of the EAS population revealed that integrating non-European genotyping data with a powerful European-based GWAS can enhance the generalizability of LDL PRS.
ABSTRACT
BACKGROUND AND AIMS: Summarised in polygenic risk scores (PRS), the effect of common, low penetrant genetic variants associated with colorectal cancer (CRC), can be used for risk stratification. METHODS: To assess the combined impact of the PRS and other main factors on CRC risk, 163,516 individuals from the UK Biobank were stratified as follows: 1. carriers status for germline pathogenic variants (PV) in CRC susceptibility genes (APC, MLH1, MSH2, MSH6, PMS2), 2. low (< 20%), intermediate (20-80%), or high PRS (> 80%), and 3. family history (FH) of CRC. Multivariable logistic regression and Cox proportional hazards models were applied to compare odds ratios and to compute the lifetime incidence, respectively. RESULTS: Depending on the PRS, the CRC lifetime incidence for non-carriers ranges between 6 and 22%, compared to 40% and 74% for carriers. A suspicious FH is associated with a further increase of the cumulative incidence reaching 26% for non-carriers and 98% for carriers. In non-carriers without FH, but high PRS, the CRC risk is doubled, whereas a low PRS even in the context of a FH results in a decreased risk. The full model including PRS, carrier status, and FH improved the area under the curve in risk prediction (0.704). CONCLUSION: The findings demonstrate that CRC risks are strongly influenced by the PRS for both a sporadic and monogenic background. FH, PV, and common variants complementary contribute to CRC risk. The implementation of PRS in routine care will likely improve personalized risk stratification, which will in turn guide tailored preventive surveillance strategies in high, intermediate, and low risk groups.
Subject(s)
Colorectal Neoplasms , Germ-Line Mutation , Humans , Incidence , Risk Factors , Germ CellsABSTRACT
The rapid accumulation of new biomedical literature not only causes curated knowledge graphs (KGs) to become outdated and incomplete, but also makes manual curation an impractical and unsustainable solution. Automated or semi-automated workflows are necessary to assist in prioritizing and curating the literature to update and enrich KGs. We have developed two workflows: one for re-curating a given KG to assure its syntactic and semantic quality and another for rationally enriching it by manually revising automatically extracted relations for nodes with low information density. We applied these workflows to the KGs encoded in Biological Expression Language from the NeuroMMSig database using content that was pre-extracted from MEDLINE abstracts and PubMed Central full-text articles using text mining output integrated by INDRA. We have made this workflow freely available at https://github.com/bel-enrichment/bel-enrichment.
Subject(s)
Data Mining , MEDLINE , Pattern Recognition, Automated , SemanticsABSTRACT
BACKGROUND: The number of diagnostic assays for the detection of herpes simplex virus type 1 (HSV-1) antibodies has increased over the years. However, their performance characteristics could vary among global populations. OBJECTIVE: To investigate performance of two commercial ELISA kits, HerpeSelect® 1 ELISA and Euroimmun Anti-HSV-1 (gC1) ELISA (IgG); and two commercial immunoblot (IB)/Western blot (WB) assays, HerpeSelect® 1 and 2 Immunoblot IgG, and Euroimmun Anti-HSV-1/HSV-2 gG2 Euroline-WB (IgG/IgM); in detecting HSV-1 antibodies in a Middle East and North Africa (MENA) population. STUDY DESIGN: Blood specimens were collected from blood donors in Doha, Qatar, June 2013-2016. Twenty specimens were randomly selected from 10 MENA nationalities (Egypt, Iran, Jordan, Lebanon, Pakistan, Palestine, Qatar, Sudan, Syria, and Yemen; totalâ¯=â¯200), and tested for HSV-1 antibodies. RESULTS: Across all six comparisons between assays, positive percent agreement ranged between 95.7% (95% CI: 91.4-98.3%) and 100.0% (95% CI: 97.8-100.0%). Negative percent agreement ranged between 86.2% (95% CI: 68.3-96.1%) and 96.2% (95% CI: 80.4-99.9%). Overall percent agreement ranged between 95.7% (95% CI: 91.7-97.8%) and 99.4% (95% CI: 96.7-99.9%). Cohen's kappa statistic ranged between 0.84 (95% CI: 0.73-0.95) and 0.98 (95% CI: 0.93-1.00). Compared against IB/WB, HerpeSelect® and Euroimmun had sensitivities and specificities >96% and >86%, respectively. Positive and negative predictive values were >97% and >83%, respectively. CONCLUSION: The assays showed excellent concordance with one another, and with a high kappa statistic. The ELISA kits demonstrated robust diagnostic performance compared to the IB/WB assays. These findings support the assays' utility in clinical diagnosis and research in MENA populations.