Network proximity analysis as a theoretical model for identifying potential novel therapies in primary sclerosing cholangitis.

Primary Sclerosing Cholangitis (PSC) is a progressive cholestatic liver disease with no licensed therapies. Previous Genome Wide Association Studies (GWAS) have identified genes that correlate significantly with PSC, and these were identified by systematic review. Here we use novel Network Proximity Analysis (NPA) methods to identify already licensed candidate drugs that may have an effect on the genetically coded aspects of PSC pathophysiology.Over 2000 agents were identified as significantly linked to genes implicated in PSC by this method. The most significant results include previously researched agents such as metronidazole, as well as biological agents such as basiliximab, abatacept and belatacept. This in silico analysis could potentially serve as a basis for developing novel clinical trials in this rare disease.

A polygenic risk score for alcohol-associated cirrhosis among heavy drinkers with European ancestry.

BACKGROUND: Polygenic Risk Scores (PRS) based on results from genome-wide association studies offer the prospect of risk stratification for many common and complex diseases. We developed a PRS for alcohol-associated cirrhosis by comparing single-nucleotide polymorphisms among patients with alcohol-associated cirrhosis (ALC) versus drinkers who did not have evidence of liver fibrosis/cirrhosis. METHODS: Using a data-driven approach, a PRS for ALC was generated using a meta-genome-wide association study of ALC (N=4305) and an independent cohort of heavy drinkers with ALC and without significant liver disease (N=3037). It was validated in 2 additional independent cohorts from the UK Biobank with diagnosed ALC (N=467) and high-risk drinking controls (N=8981) and participants in the Indiana Biobank Liver cohort with alcohol-associated liver disease (N=121) and controls without liver disease (N=3239). RESULTS: A 20-single-nucleotide polymorphisms PRS for ALC (PRSALC) was generated that stratified risk for ALC comparing the top and bottom deciles of PRS in the 2 validation cohorts (ORs: 2.83 [95% CI: 1.82 -4.39] in UK Biobank; 4.40 [1.56 -12.44] in Indiana Biobank Liver cohort). Furthermore, PRSALC improved the prediction of ALC risk when added to the models of clinically known predictors of ALC risk. It also stratified the risk for metabolic dysfunction -associated steatotic liver disease -cirrhosis (3.94 [2.23 -6.95]) in the Indiana Biobank Liver cohort -based exploratory analysis. CONCLUSIONS: PRSALC incorporates 20 single-nucleotide polymorphisms, predicts increased risk for ALC, and improves risk stratification for ALC compared with the models that only include clinical risk factors. This new score has the potential for early detection of heavy drinking patients who are at high risk for ALC.

Genetic complexity of diagnostically unresolved Ehlers-Danlos syndrome.

BACKGROUND: The Ehlers-Danlos syndromes (EDS) are heritable disorders of connective tissue (HDCT), reclassified in the 2017 nosology into 13 subtypes. The genetic basis for hypermobile Ehlers-Danlos syndrome (hEDS) remains unknown. METHODS: Whole exome sequencing (WES) was undertaken on 174 EDS patients recruited from a national diagnostic service for complex EDS and a specialist clinic for hEDS. Patients had already undergone expert phenotyping, laboratory investigation and gene sequencing, but were without a genetic diagnosis. Filtered WES data were reviewed for genes underlying Mendelian disorders and loci reported in EDS linkage, transcriptome and genome-wide association studies (GWAS). A genetic burden analysis (Minor Allele Frequency (MAF) <0.05) incorporating 248 Avon Longitudinal Study of Parents and Children (ALSPAC) controls sequenced as part of the UK10K study was undertaken using TASER methodology. RESULTS: Heterozygous pathogenic (P) or likely pathogenic (LP) variants were identified in known EDS and Loeys-Dietz (LDS) genes. Multiple variants of uncertain significance where segregation and functional analysis may enable reclassification were found in genes associated with EDS, LDS, heritable thoracic aortic disease (HTAD), Mendelian disorders with EDS symptomatology and syndromes with EDS-like features. Genetic burden analysis revealed a number of novel loci, although none reached the threshold for genome-wide significance. Variants with biological plausibility were found in genes and pathways not currently associated with EDS or HTAD. CONCLUSIONS: We demonstrate the clinical utility of large panel-based sequencing and WES for patients with complex EDS in distinguishing rare EDS subtypes, LDS and related syndromes. Although many of the P and LP variants reported in this cohort would be identified with current panel testing, they were not at the time of this study, highlighting the use of extended panels and WES as a clinical tool for complex EDS. Our results are consistent with the complex genetic architecture of EDS and suggest a number of novel hEDS and HTAD candidate genes and pathways.

The Application of Genetic Risk Scores in Rheumatic Diseases: A Perspective.

Modest effect sizes have limited the clinical applicability of genetic associations with rheumatic diseases. Genetic risk scores (GRSs) have emerged as a promising solution to translate genetics into useful tools. In this review, we provide an overview of the recent literature on GRSs in rheumatic diseases. We describe six categories for which GRSs are used: (a) disease (outcome) prediction, (b) genetic commonalities between diseases, (c) disease differentiation, (d) interplay between genetics and environmental factors, (e) heritability and transferability, and (f) detecting causal relationships between traits. In our review of the literature, we identified current lacunas and opportunities for future work. First, the shortage of non-European genetic data restricts the application of many GRSs to European populations. Next, many GRSs are tested in settings enriched for cases that limit the transferability to real life. If intended for clinical application, GRSs are ideally tested in the relevant setting. Finally, there is much to elucidate regarding the co-occurrence of clinical traits to identify shared causal paths and elucidate relationships between the diseases. GRSs are useful instruments for this. Overall, the ever-continuing research on GRSs gives a hopeful outlook into the future of GRSs and indicates significant progress in their potential applications.

Comparison of regmed and BayesNetty for exploring causal models with many variables.

Here we compare a recently proposed method and software package, regmed, with our own previously developed package, BayesNetty, designed to allow exploratory analysis of complex causal relationships between biological variables. We find that regmed generally has poorer recall but much better precision than BayesNetty. This is perhaps not too surprising as regmed is specifically designed for use with high-dimensional data. BayesNetty is found to be more sensitive to the resulting multiple testing problem encountered in these circumstances. However, as regmed is not designed to handle missing data, its performance is severely affected when missing data is present, whereas the performance of BayesNetty is only slightly affected. The performance of regmed can be rescued in this situation by first using BayesNetty to impute the missing data, and then applying regmed to the resulting "filled-in" data set.

Oligogenic Architecture of Rare Noncoding Variants Distinguishes 4 Congenital Heart Disease Phenotypes.

BACKGROUND: Congenital heart disease (CHD) is highly heritable, but the power to identify inherited risk has been limited to analyses of common variants in small cohorts. METHODS: We performed reimputation of 4 CHD cohorts (n=55 342) to the TOPMed reference panel (freeze 5), permitting meta-analysis of 14 784 017 variants including 6 035 962 rare variants of high imputation quality as validated by whole genome sequencing. RESULTS: Meta-analysis identified 16 novel loci, including 12 rare variants, which displayed moderate or large effect sizes (median odds ratio, 3.02) for 4 separate CHD categories. Analyses of chromatin structure link 13 of the genome-wide significant loci to key genes in cardiac development; rs373447426 (minor allele frequency, 0.003 [odds ratio, 3.37 for Conotruncal heart disease]; P=1.49×10-8) is predicted to disrupt chromatin structure for 2 nearby genes BDH1 and DLG1 involved in Conotruncal development. A lead variant rs189203952 (minor allele frequency, 0.01 [odds ratio, 2.4 for left ventricular outflow tract obstruction]; P=1.46×10-8) is predicted to disrupt the binding sites of 4 transcription factors known to participate in cardiac development in the promoter of SPAG9. A tissue-specific model of chromatin conformation suggests that common variant rs78256848 (minor allele frequency, 0.11 [odds ratio, 1.4 for Conotruncal heart disease]; P=2.6×10-8) physically interacts with NCAM1 (PFDR=1.86×10-27), a neural adhesion molecule acting in cardiac development. Importantly, while each individual malformation displayed substantial heritability (observed h2 ranging from 0.26 for complex malformations to 0.37 for left ventricular outflow tract obstructive disease) the risk for different CHD malformations appeared to be separate, without genetic correlation measured by linkage disequilibrium score regression or regional colocalization. CONCLUSIONS: We describe a set of rare noncoding variants conferring significant risk for individual heart malformations which are linked to genes governing cardiac development. These results illustrate that the oligogenic basis of CHD and significant heritability may be linked to rare variants outside protein-coding regions conferring substantial risk for individual categories of cardiac malformation.

Regulation of immune responses in primary biliary cholangitis: a transcriptomic analysis of peripheral immune cells.

BACKGROUND AIMS: In patients with primary biliary cholangitis (PBC), the serum liver biochemistry measured during treatment with ursodeoxycholic acid-the UDCA response-accurately predicts long-term outcome. Molecular characterization of patients stratified by UDCA response can improve biological understanding of the high-risk disease, thereby helping to identify alternative approaches to disease-modifying therapy. In this study, we sought to characterize the immunobiology of the UDCA response using transcriptional profiling of peripheral blood mononuclear cell subsets. METHODS: We performed bulk RNA-sequencing of monocytes and TH1, TH17, TREG, and B cells isolated from the peripheral blood of 15 PBC patients with adequate UDCA response ("responders"), 16 PBC patients with inadequate UDCA response ("nonresponders"), and 15 matched controls. We used the Weighted Gene Co-expression Network Analysis to identify networks of co-expressed genes ("modules") associated with response status and the most highly connected genes ("hub genes") within them. Finally, we performed a Multi-Omics Factor Analysis of the Weighted Gene Co-expression Network Analysis modules to identify the principal axes of biological variation ("latent factors") across all peripheral blood mononuclear cell subsets. RESULTS: Using the Weighted Gene Co-expression Network Analysis, we identified modules associated with response and/or disease status (q<0.05) in each peripheral blood mononuclear cell subset. Hub genes and functional annotations suggested that monocytes are proinflammatory in nonresponders, but antiinflammatory in responders; TH1 and TH17 cells are activated in all PBC cases but better regulated in responders; and TREG cells are activated-but also kept in check-in responders. Using the Multi-Omics Factor Analysis, we found that antiinflammatory activity in monocytes, regulation of TH1 cells, and activation of TREG cells are interrelated and more prominent in responders. CONCLUSIONS: We provide evidence that adaptive immune responses are better regulated in patients with PBC with adequate UDCA response.

Multitrait genome-wide analyses identify new susceptibility loci and candidate drugs to primary sclerosing cholangitis.

Primary sclerosing cholangitis (PSC) is a rare autoimmune bile duct disease that is strongly associated with immune-mediated disorders. In this study, we implemented multitrait joint analyses to genome-wide association summary statistics of PSC and numerous clinical and epidemiological traits to estimate the genetic contribution of each trait and genetic correlations between traits and to identify new lead PSC risk-associated loci. We identified seven new loci that have not been previously reported and one new independent lead variant in the previously reported locus. Functional annotation and fine-mapping nominated several potential susceptibility genes such as MANBA and IRF5. Network-based in silico drug efficacy screening provided candidate agents for further study of pharmacological effect in PSC.

Biallelic variants in CEP164 cause a motile ciliopathy-like syndrome.

Ciliopathies may be classed as primary or motile depending on the underlying ciliary defect and are usually considered distinct clinical entities. Primary ciliopathies are associated with multisystem syndromes typically affecting the brain, kidney, and eye, as well as other organ systems such as the liver, skeleton, auditory system, and metabolism. Motile ciliopathies are a heterogenous group of disorders with defects in specialised motile ciliated tissues found within the lung, brain, and reproductive system, and are associated with primary ciliary dyskinesia, bronchiectasis, infertility and rarely hydrocephalus. Primary and motile cilia share defined core ultra-structures with an overlapping proteome, and human disease phenotypes can reflect both primary and motile ciliopathies. CEP164 encodes a centrosomal distal appendage protein vital for primary ciliogenesis. Human CEP164 mutations are typically described in patients with nephronophthisis-related primary ciliopathies but have also been implicated in motile ciliary dysfunction. Here we describe a patient with an atypical motile ciliopathy phenotype and biallelic CEP164 variants. This work provides further evidence that CEP164 mutations can contribute to both primary and motile ciliopathy syndromes, supporting their functional and clinical overlap, and informs the investigation and management of CEP164 ciliopathy patients.

Rationale, Design, and the Baseline Characteristics of the RHDGen (The Genetics of Rheumatic Heart Disease) Network Study†.

BACKGROUND: The genetics of rheumatic heart disease (RHDGen) Network was developed to assist the discovery and validation of genetic variations and biomarkers of risk for rheumatic heart disease (RHD) in continental Africans, as a part of the global fight to control and eradicate rheumatic fever/RHD. Thus, we describe the rationale and design of the RHDGen study, comprising participants from 8 African countries. METHODS: RHDGen screened potential participants using echocardiography, thereafter enrolling RHD cases and ethnically-matched controls for whom case characteristics were documented. Biological samples were collected for conducting genetic analyses, including a discovery case-control genome-wide association study (GWAS) and a replication trio family study. Additional biological samples were also collected, and processed, for the measurement of biomarker analytes and the biomarker analyses are underway. RESULTS: Participants were enrolled into RHDGen between December 2012 and March 2018. For GWAS, 2548 RHD cases and 2261 controls (3301 women [69%]; mean age [SD], 37 [16.3] years) were available. RHD cases were predominantly Black (66%), Admixed (24%), and other ethnicities (10%). Among RHD cases, 34% were asymptomatic, 26% had prior valve surgery, and 23% had atrial fibrillation. The trio family replication arm included 116 RHD trio probands and 232 parents. CONCLUSIONS: RHDGen presents a rare opportunity to identify relevant patterns of genetic factors and biomarkers in Africans that may be associated with differential RHD risk. Furthermore, the RHDGen Network provides a platform for further work on fully elucidating the causes and mechanisms associated with RHD susceptibility and development.

Comprehensive genetic and functional analyses of Fc gamma receptors influence on response to rituximab therapy for autoimmunity.

BACKGROUND: Rituximab is widely used to treat autoimmunity but clinical response varies. Efficacy is determined by the efficiency of B-cell depletion, which may depend on various Fc gamma receptor (FcγR)-dependent mechanisms. Study of FcγR is challenging due to the complexity of the FCGR genetic locus. We sought to assess the effect of FCGR variants on clinical response, B-cell depletion and NK-cell-mediated killing in rheumatoid arthritis (RA) and systemic lupus erythematosus (SLE). METHODS: A longitudinal cohort study was conducted in 835 patients [RA = 573; SLE = 262]. Clinical outcome measures were two-component disease activity score in 28-joints (2C-DAS28CRP) for RA and British Isles Lupus Assessment Group (BILAG)-2004 major clinical response (MCR) for SLE at 6 months. B-cells were evaluated by highly-sensitive flow cytometry. Single nucleotide polymorphism and copy number variation for genes encoding five FcγRs were measured using multiplex ligation-dependent probe amplification. Ex vivo studies assessed NK-cell antibody-dependent cellular cytotoxicity (ADCC) and FcγR expression. FINDINGS: In RA, carriage of FCGR3A-158V and increased FCGR3A-158V copies were associated with greater 2C-DAS28CRP response (adjusted for baseline 2C-DAS28CRP). In SLE, MCR was associated with increased FCGR3A-158V, OR 1.64 (95% CI 1.12-2.41) and FCGR2C-ORF OR 1.93 (95% CI 1.09-3.40) copies. 236/413 (57%) patients with B-cell data achieved complete depletion. Homozygosity for FCGR3A-158V and increased FCGR3A-158V copies were associated with complete depletion in combined analyses. FCGR3A genotype was associated with rituximab-induced ADCC, and increased NK-cell FcγRIIIa expression was associated with improved clinical response and depletion in vivo. Furthermore, disease status and concomitant therapies impacted both NK-cell FcγRIIIa expression and ADCC. INTERPRETATION: FcγRIIIa is the major low affinity FcγR associated with rituximab response. Increased copies of the FCGR3A-158V allele (higher affinity for IgG1), influences clinical and biological responses to rituximab in autoimmunity. Enhancing FcγR-effector functions could improve the next generation of CD20-depleting therapies and genotyping may stratify patients for optimal treatment protocols. FUNDING: Medical Research Council, National Institute for Health and Care Research, Versus Arthritis.

Confirmation of the superior performance of the causal Graphical Analysis Using Genetics (cGAUGE) pipeline in comparison to various competing alternatives.

Various methods exist that utilise information from genetic predictors to help identify potential causal relationships between measured biological or clinical traits. Here we conduct computer simulations to investigate the performance of a recently proposed causal Graphical Analysis Using Genetics (cGAUGE) pipeline, used as a precursor to Mendelian randomization analysis, in comparison to our previously proposed Bayesian Network approach for addressing this problem. We use the same simulation (and analysis) code as was used by the developers of cGAUGE, adding in a comparison with the Bayesian Network approach. Overall, we find the optimal method (in terms of giving high power and low false discovery rate) is the cGAUGE pipeline followed by subsequent analysis using the MR-PRESSO Mendelian randomization approach.

Investigating the prediction of CpG methylation levels from SNP genotype data to help elucidate relationships between methylation, gene expression and complex traits.

As popularised by PrediXcan (and related methods), transcriptome-wide association studies (TWAS), in which gene expression is imputed from single-nucleotide polymorphism (SNP) genotypes and tested for association with a phenotype, are a popular approach for investigating the role of gene expression in complex traits. Like gene expression, DNA methylation is an important biological process and, being under genetic regulation, may be imputable from SNP genotypes. Here, we investigate prediction of CpG methylation levels from SNP genotype data to help elucidate relationships between methylation, gene expression and complex traits. We start by examining how well CpG methylation can be predicted from SNP genotypes, comparing three penalised regression approaches and examining whether changing the window size improves prediction accuracy. Although methylation at most CpG sites cannot be accurately predicted from SNP genotypes, for a subset it can be predicted well. We next apply our methylation prediction models (trained using the optimal method and window size) to carry out a methylome-wide association study (MWAS) of primary biliary cholangitis. We intersect the regions identified via MWAS with those identified via TWAS, providing insight into the interplay between CpG methylation, gene expression and disease status. We conclude that MWAS has the potential to improve understanding of biological mechanisms in complex traits.

A multiancestry genome-wide association study of unexplained chronic ALT elevation as a proxy for nonalcoholic fatty liver disease with histological and radiological validation.

Nonalcoholic fatty liver disease (NAFLD) is a growing cause of chronic liver disease. Using a proxy NAFLD definition of chronic elevation of alanine aminotransferase (cALT) levels without other liver diseases, we performed a multiancestry genome-wide association study (GWAS) in the Million Veteran Program (MVP) including 90,408 cALT cases and 128,187 controls. Seventy-seven loci exceeded genome-wide significance, including 25 without prior NAFLD or alanine aminotransferase associations, with one additional locus identified in European American-only and two in African American-only analyses (P < 5 × 10-8). External replication in histology-defined NAFLD cohorts (7,397 cases and 56,785 controls) or radiologic imaging cohorts (n = 44,289) replicated 17 single-nucleotide polymorphisms (SNPs) (P < 6.5 × 10-4), of which 9 were new (TRIB1, PPARG, MTTP, SERPINA1, FTO, IL1RN, COBLL1, APOH and IFI30). Pleiotropy analysis showed that 61 of 77 multiancestry and all 17 replicated SNPs were jointly associated with metabolic and/or inflammatory traits, revealing a complex model of genetic architecture. Our approach integrating cALT, histology and imaging reveals new insights into genetic liability to NAFLD.

Rare complement factor I variants associated with reduced macular thickness and age-related macular degeneration in the UK Biobank.

To evaluate potential diagnostic and therapeutic biomarkers for age-related macular degeneration (AMD), we identified 8433 UK Biobank participants with rare complement Factor I gene (CFI) variants, 579 with optical coherence tomography-derived macular thickness data. We stratified these variants by predicted gene expression and measured their association with retinal pigment epithelium-Bruch's membrane (RPE-BM) complex and retinal thicknesses at nine macular subfields, as well as AMD risk, using multivariable regression models adjusted for the common complement Factor H gene (CFH) p.Y402H and age-related maculopathy susceptibility protein 2 gene (ARMS2) p.A69S risk genotypes. CFI variants associated with low Factor I levels predicted a thinner mean RPE-BM (95% confidence interval [CI] -1.66 to -0.37 µm, P = 0.002) and retina (95% CI -5.88 to -0.13 µm, P = 0.04) and a higher AMD risk (odds ratio [OR] = 2.26, 95% CI 1.56 to 3.27, P < 0.001). CFI variants associated with normal Factor I levels did not impact mean RPE-BM/retinal thickness (P = 0.28; P = 0.99) or AMD risk (P = 0.97). CFH p.Y402H was associated with a thinner RPE-BM (95% CI -0.31 to -0.18 µm, P < 0.001 heterozygous; 95% CI -0.62 to -0.42 µm, P < 0.001 homozygous) and retina (95% CI -0.73 to -0.12 µm, P = 0.007 heterozygous; 95% CI -1.08 to -0.21 µm, P = 0.004 homozygous). ARMS2 p.A69S did not influence RPE-BM (P = 0.80 heterozygous; P = 0.12 homozygous) or retinal thickness (P = 0.75 heterozygous; P = 0.07 homozygous). p.Y402H and p.A69S exhibited a significant allele-dose response with AMD risk. Thus, CFI rare variants associated with low Factor I levels are robust predictors of reduced macular thickness and AMD. The observed association between macular thickness and CFH p.Y402H, but not ARMS2 p.A69S, highlights the importance of complement dysregulation in early pathogenesis.

A genetic risk score and diabetes predict development of alcohol-related cirrhosis in drinkers.

BACKGROUND & AIMS: Only a minority of excess alcohol drinkers develop cirrhosis. We developed and evaluated risk stratification scores to identify those at highest risk. METHODS: Three cohorts (GenomALC-1: n = 1,690, GenomALC-2: n = 3,037, UK Biobank: relevant n = 6,898) with a history of heavy alcohol consumption (≥80 g/day (men), ≥50 g/day (women), for ≥10 years) were included. Cases were participants with alcohol-related cirrhosis. Controls had a history of similar alcohol consumption but no evidence of liver disease. Risk scores were computed from up to 8 genetic loci identified previously as associated with alcohol-related cirrhosis and 3 clinical risk factors. Score performance for the stratification of alcohol-related cirrhosis risk was assessed and compared across the alcohol-related liver disease spectrum, including hepatocellular carcinoma (HCC). RESULTS: A combination of 3 single nucleotide polymorphisms (SNPs) (PNPLA3:rs738409, SUGP1-TM6SF2:rs10401969, HSD17B13:rs6834314) and diabetes status best discriminated cirrhosis risk. The odds ratios (ORs) and (95% CIs) between the lowest (Q1) and highest (Q5) score quintiles of the 3-SNP score, based on independent allelic effect size estimates, were 5.99 (4.18-8.60) (GenomALC-1), 2.81 (2.03-3.89) (GenomALC-2), and 3.10 (2.32-4.14) (UK Biobank). Patients with diabetes and high risk scores had ORs of 14.7 (7.69-28.1) (GenomALC-1) and 17.1 (11.3-25.7) (UK Biobank) compared to those without diabetes and with low risk scores. Patients with cirrhosis and HCC had significantly higher mean risk scores than patients with cirrhosis alone (0.76 ± 0.06 vs. 0.61 ± 0.02, p = 0.007). Score performance was not significantly enhanced by information on additional genetic risk variants, body mass index or coffee consumption. CONCLUSIONS: A risk score based on 3 genetic risk variants and diabetes status enables the stratification of heavy drinkers based on their risk of cirrhosis, allowing for the provision of earlier preventative interventions. LAY SUMMARY: Excessive chronic drinking leads to cirrhosis in some people, but so far there is no way to identify those at high risk of developing this debilitating disease. We developed a genetic risk score that can identify patients at high risk. The risk of cirrhosis is increased >10-fold with just two risk factors - diabetes and a high genetic risk score. Risk assessment using this test could enable the early and personalised management of this disease in high-risk patients.