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.

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.

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.

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.

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.

Macrophage scavenger receptor 1 mediates lipid-induced inflammation in non-alcoholic fatty liver disease.

BACKGROUND & AIMS: Obesity-associated inflammation is a key player in the pathogenesis of non-alcoholic fatty liver disease (NAFLD). However, the role of macrophage scavenger receptor 1 (MSR1, CD204) remains incompletely understood. METHODS: A total of 170 NAFLD liver biopsies were processed for transcriptomic analysis and correlated with clinicopathological features. Msr1-/- and wild-type mice were subjected to a 16-week high-fat and high-cholesterol diet. Mice and ex vivo human liver slices were treated with a monoclonal antibody against MSR1. Genetic susceptibility was assessed using genome-wide association study data from 1,483 patients with NAFLD and 430,101 participants of the UK Biobank. RESULTS: MSR1 expression was associated with the occurrence of hepatic lipid-laden foamy macrophages and correlated with the degree of steatosis and steatohepatitis in patients with NAFLD. Mice lacking Msr1 were protected against diet-induced metabolic disorder, showing fewer hepatic foamy macrophages, less hepatic inflammation, improved dyslipidaemia and glucose tolerance, and altered hepatic lipid metabolism. Upon induction by saturated fatty acids, MSR1 induced a pro-inflammatory response via the JNK signalling pathway. In vitro blockade of the receptor prevented the accumulation of lipids in primary macrophages which inhibited the switch towards a pro-inflammatory phenotype and the release of cytokines such as TNF-ɑ. Targeting MSR1 using monoclonal antibody therapy in an obesity-associated NAFLD mouse model and human liver slices resulted in the prevention of foamy macrophage formation and inflammation. Moreover, we identified that rs41505344, a polymorphism in the upstream transcriptional region of MSR1, was associated with altered serum triglycerides and aspartate aminotransferase levels in a cohort of over 400,000 patients. CONCLUSIONS: Taken together, our data suggest that MSR1 plays a critical role in lipid-induced inflammation and could thus be a potential therapeutic target for the treatment of NAFLD. LAY SUMMARY: Non-alcoholic fatty liver disease (NAFLD) is a chronic disease primarily caused by excessive consumption of fat and sugar combined with a lack of exercise or a sedentary lifestyle. Herein, we show that the macrophage scavenger receptor MSR1, an innate immune receptor, mediates lipid uptake and accumulation in Kupffer cells, resulting in liver inflammation and thereby promoting the progression of NAFLD in humans and mice.

No Evidence That Genetic Variation at the Klotho Locus Is Associated With Longevity in Caucasians from the Newcastle 85+ Study and the UK Biobank.

The demographics of Western populations are changing, with an increase in the proportion of older adults. There is evidence to suggest that genetic factors may influence the aging process: studying these may lead to interventions to help individuals live a longer and healthier life. Evidence from several groups indicates that Klotho (KL), a gene encoding a single-pass transmembrane protein that acts as an FGF23 co-receptor, may be associated with longevity and healthy aging. We aimed to explore this area further by comparing the genotype counts in 642 long-lived individuals from the Newcastle 85+ Study with 18 295 middle-aged Newcastle-based controls from the UK Biobank to test whether variants at the KL gene locus are over- or under-represented in older individuals. If KL is associated with longevity, then we would expect the genotype counts to differ between the 2 cohorts. We found that the rs2283368 CC genotype and the rs9536338 C allele, but not the KL-VS haplotype, were associated with reaching very old age. However, these associations did not replicate in the remainder of the UK Biobank cohort. Thus, our results do not reliably support the role of KL as a longevity factor.