Digenic inheritance involving a muscle-specific protein kinase and the giant titin protein causes a skeletal muscle myopathy

In digenic inheritance, pathogenic variants in two genes must be inherited together to cause disease. Only very few examples of digenic inheritance have been described in the neuromuscular disease field. Here we show that predicted deleterious variants in SRPK3, encoding the X-linked serine/argenine protein kinase 3, lead to a progressive early onset skeletal muscle myopathy only when in combination with heterozygous variants in the TTN gene. The co-occurrence of predicted deleterious SRPK3/TTN variants was not seen among 76,702 healthy male individuals, and statistical modeling strongly supported digenic inheritance as the best-fitting model. Furthermore, double-mutant zebrafish (srpk3-/-; ttn.1+/-) replicated the myopathic phenotype and showed myofibrillar disorganization. Transcriptome data suggest that the interaction of srpk3 and ttn.1 in zebrafish occurs at a post-transcriptional level. We propose that digenic inheritance of deleterious changes impacting both the protein kinase SRPK3 and the giant muscle protein titin causes a skeletal myopathy and might serve as a model for other genetic diseases.

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.

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.

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.

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.

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.

All-cause and liver-related mortality risk factors in excessive drinkers: Analysis of data from the UK biobank

BACKGROUND

High alcohol intake is associated with increased mortality. We aimed to identify factors affecting mortality in people drinking extreme amounts of alcohol.

METHODS

We obtained information from the UK Biobank on approximately 500,000 participants aged 40-70 years at baseline assessment in 2006-2010. Habitual alcohol intake, lifestyle and physiological data, laboratory test results, and hospital diagnoses and death certificate data (to June 2020) for 5136 men (2.20% of male participants) and 1504 women (0.60%) who reported consuming ≥80 or ≥50 g/day, respectively, were used in survival analysis.

RESULTS

Mortality hazard ratios for these excessive drinkers, compared to all other participants, were 2.02 (95% CI 1.89-2.17) for all causes, 1.89 (1.69-2.12) for any cancer, 1.87 (1.61-2.17) for any circulatory disease, and 9.40 (7.00-12.64) for any liver disease. Liver disease diagnosis or abnormal liver function tests predicted not only deaths attributed to liver disease but also those from cancers or circulatory diseases. Mortality among excessive drinkers was also associated with quantitative alcohol intake; diagnosed alcohol dependence, harmful use, or withdrawal syndrome; and current smoking at assessment.

CONCLUSIONS

People with chronic excessive alcohol intake experience decreased average survival, but there is substantial variation in their mortality, with liver abnormality and alcohol dependence or other alcohol use disorders associated with a worse prognosis. Clinically, patients with these risk factors and high alcohol intake should be considered for early or intensive management. Research can usefully focus on the factors predisposing to dependence or liver abnormality.

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.

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.

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.

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.

Common Genetic Variants Contribute to Risk of Transposition of the Great Arteries

RATIONALE

Dextro-transposition of the great arteries (D-TGA) is a severe congenital heart defect which affects approximately 1 in 4,000 live births. While there are several reports of D-TGA patients with rare variants in individual genes, the majority of D-TGA cases remain genetically elusive. Familial recurrence patterns and the observation that most cases with D-TGA are sporadic suggest a polygenic inheritance for the disorder, yet this remains unexplored.

OBJECTIVE

We sought to study the role of common single nucleotide polymorphisms (SNPs) in risk for D-TGA.

METHODS AND RESULTS

We conducted a genome-wide association study in an international set of 1,237 patients with D-TGA and identified a genome-wide significant susceptibility locus on chromosome 3p14.3, which was subsequently replicated in an independent case-control set (rs56219800, meta-analysis P=8.6x10-10, OR=0.69 per C allele). SNP-based heritability analysis showed that 25% of variance in susceptibility to D-TGA may be explained by common variants. A genome-wide polygenic risk score derived from the discovery set was significantly associated to D-TGA in the replication set (P=4x10-5). The genome-wide significant locus (3p14.3) co-localizes with a putative regulatory element that interacts with the promoter of WNT5A, which encodes the Wnt Family Member 5A protein known for its role in cardiac development in mice. We show that this element drives reporter gene activity in the developing heart of mice and zebrafish and is bound by the developmental transcription factor TBX20. We further demonstrate that TBX20 attenuates Wnt5a expression levels in the developing mouse heart.

CONCLUSIONS

This work provides support for a polygenic architecture in D-TGA and identifies a susceptibility locus on chromosome 3p14.3 near WNT5A. Genomic and functional data support a causal role of WNT5A at the locus.

No Evidence That Genetic Variation At The Klotho Locus Is Associated With Longevity In Caucasians From The Newcastle 85 Plus 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.

Common and Rare Genetic Variants That Could Contribute to Severe Otitis Media in an Australian Aboriginal Population

BACKGROUND

Our goal was to identify genetic risk factors for severe otitis media (OM) in Aboriginal Australians.

METHODS

Illumina® Omni2.5 BeadChip and imputed data were compared between 21 children with severe OM (multiple episodes chronic suppurative OM and/or perforations or tympanic sclerosis) and 370 individuals without this phenotype, followed by FUnctional Mapping and Annotation (FUMA). Exome data filtered for common (EXaC_all ≥ 0.1) putative deleterious variants influencing protein coding (CADD-scaled scores ≥15] were used to compare 15 severe OM cases with 9 mild cases (single episode of acute OM recorded over ≥3 consecutive years). Rare (ExAC_all ≤ 0.01) such variants were filtered for those present only in severe OM. Enrichr was used to determine enrichment of genes contributing to pathways/processes relevant to OM.

RESULTS

FUMA analysis identified 2 plausible genetic risk loci for severe OM: NR3C1 (Pimputed_1000G = 3.62 × 10-6) encoding the glucocorticoid receptor, and NREP (Pimputed_1000G = 3.67 × 10-6) encoding neuronal regeneration-related protein. Exome analysis showed: (i) association of severe OM with variants influencing protein coding (CADD-scaled ≥ 15) in a gene-set (GRXCR1, CDH23, LRP2, FAT4, ARSA, EYA4) enriched for Mammalian Phenotype Level 4 abnormal hair cell stereociliary bundle morphology and related phenotypes; (ii) rare variants influencing protein coding only seen in severe OM provided gene-sets enriched for "abnormal ear" (LMNA, CDH23, LRP2, MYO7A, FGFR1), integrin interactions, transforming growth factor signaling, and cell projection phenotypes including hair cell stereociliary bundles and cilium assembly.

CONCLUSIONS

This study highlights interacting genes and pathways related to cilium structure and function that may contribute to extreme susceptibility to OM in Aboriginal Australian children.

Genome-wide association study identifies susceptibility loci for acute myeloid leukemia

Acute myeloid leukemia (AML) is a hematological malignancy with an undefined heritable risk. Here we perform a meta-analysis of three genome-wide association studies, with replication in a fourth study, incorporating a total of 4018 AML cases and 10488 controls. We identify a genome-wide significant risk locus for AML at 11q13.2 (rs4930561; P = 2.15 × 10-8; KMT5B). We also identify a genome-wide significant risk locus for the cytogenetically normal AML sub-group (N = 1287) at 6p21.32 (rs3916765; P = 1.51 × 10-10; HLA). Our results inform on AML etiology and identify putative functional genes operating in histone methylation (KMT5B) and immune function (HLA).

A Bayesian network approach incorporating imputation of missing data enables exploratory analysis of complex causal biological relationships

Bayesian networks can be used to identify possible causal relationships between variables based on their conditional dependencies and independencies, which can be particularly useful in complex biological scenarios with many measured variables. Here we propose two improvements to an existing method for Bayesian network analysis, designed to increase the power to detect potential causal relationships between variables (including potentially a mixture of both discrete and continuous variables). Our first improvement relates to the treatment of missing data. When there is missing data, the standard approach is to remove every individual with any missing data before performing analysis. This can be wasteful and undesirable when there are many individuals with missing data, perhaps with only one or a few variables missing. This motivates the use of imputation. We present a new imputation method that uses a version of nearest neighbour imputation, whereby missing data from one individual is replaced with data from another individual, their nearest neighbour. For each individual with missing data, the subsets of variables to be used to select the nearest neighbour are chosen by sampling without replacement the complete data and estimating a best fit Bayesian network. We show that this approach leads to marked improvements in the recall and precision of directed edges in the final network identified, and we illustrate the approach through application to data from a recent study investigating the causal relationship between methylation and gene expression in early inflammatory arthritis patients. We also describe a second improvement in the form of a pseudo-Bayesian approach for upweighting certain network edges, which can be useful when there is prior evidence concerning their directions.

Data analysis using Bayesian networks can help identify possible causal relationships between measured biological variables. Here we propose two improvements to an existing method for Bayesian network analysis. Our first improvement relates to the treatment of missing data. When there is missing data, the standard approach is to remove every individual with any missing data before performing analysis, even if only one or a few variables are missing. This is undesirable as it can reduce the ability of the approach to infer correct relationships. We propose a new method to instead fill in (impute) the missing data prior to analysis. We show through computer simulations that our method improves the reliability of the results obtained, and we illustrate the proposed approach by applying it to data from a recent study in early inflammatory arthritis. We also describe a second improvement involving the upweighting of certain network edges, which can be useful when there is prior evidence concerning their directions.

An international genome-wide meta-analysis of primary biliary cholangitis: Novel risk loci and candidate drugs

BACKGROUNDS & AIMS

Primary biliary cholangitis (PBC) is a chronic liver disease in which autoimmune destruction of the small intrahepatic bile ducts eventually leads to cirrhosis. Many patients have inadequate response to licensed medications, motivating the search for novel therapies. Previous genome-wide association studies (GWAS) and meta-analyses (GWMA) of PBC have identified numerous risk loci for this condition, providing insight into its aetiology. We undertook the largest GWMA of PBC to date, aiming to identify additional risk loci and prioritise candidate genes for in silico drug efficacy screening.

METHODS

We combined new and existing genotype data for 10,516 cases and 20,772 controls from 5 European and 2 East Asian cohorts.

RESULTS

We identified 56 genome-wide significant loci (20 novel) including 46 in European, 13 in Asian, and 41 in combined cohorts; and a 57th genome-wide significant locus (also novel) in conditional analysis of the European cohorts. Candidate genes at newly identified loci include FCRL3, INAVA, PRDM1, IRF7, CCR6, CD226, and IL12RB1, which each play key roles in immunity. Pathway analysis reiterated the likely importance of pattern recognition receptor and TNF signalling, JAK-STAT signalling, and differentiation of T helper (TH)1 and TH17 cells in the pathogenesis of this disease. Drug efficacy screening identified several medications predicted to be therapeutic in PBC, some of which are well-established in the treatment of other autoimmune disorders.

CONCLUSIONS

This study has identified additional risk loci for PBC, provided a hierarchy of agents that could be trialled in this condition, and emphasised the value of genetic and genomic approaches to drug discovery in complex disorders.

LAY SUMMARY

Primary biliary cholangitis (PBC) is a chronic liver disease that eventually leads to cirrhosis. In this study, we analysed genetic information from 10,516 people with PBC and 20,772 healthy individuals recruited in Canada, China, Italy, Japan, the UK, or the USA. We identified several genetic regions associated with PBC. Each of these regions contains several genes. For each region, we used diverse sources of evidence to help us choose the gene most likely to be involved in causing PBC. We used these 'candidate genes' to help us identify medications that are currently used for treatment of other conditions, which might also be useful for treatment of PBC.

Association of Novel Locus With Rheumatic Heart Disease in Black African Individuals: Findings From the RHDGen Study

Importance

Rheumatic heart disease (RHD), a sequela of rheumatic fever characterized by permanent heart valve damage, is the leading cause of cardiac surgery in Africa. However, its pathophysiologic characteristics and genetics are poorly understood. Understanding genetic susceptibility may aid in prevention, control, and interventions to eliminate RHD.

Objective

To identify common genetic loci associated with RHD susceptibility in Black African individuals.

Design, Setting, and Participants

This multicenter case-control genome-wide association study (GWAS), the Genetics of Rheumatic Heart Disease, examined more than 7 million genotyped and imputed single-nucleotide variations. The 4809 GWAS participants and 116 independent trio families were enrolled from 8 African countries between December 31, 2012, and March 31, 2018. All GWAS participants and trio probands were screened by use of echocardiography. Data analyses took place from May 15, 2017, until March 14, 2021.

Main Outcomes and Measures

Genetic associations with RHD.

Results

This study included 4809 African participants (2548 RHD cases and 2261 controls; 3301 women [69%]; mean [SD] age, 36.5 [16.3] years). The GWAS identified a single RHD risk locus, 11q24.1 (rs1219406 [odds ratio, 1.65; 95% CI, 1.48-1.82; P = 4.36 × 10-8]), which reached genome-wide significance in Black African individuals. Our meta-analysis of Black (n = 3179) and admixed (n = 1055) African individuals revealed several suggestive loci. The study also replicated a previously reported association in Pacific Islander individuals (rs11846409) at the immunoglobulin heavy chain locus, in the meta-analysis of Black and admixed African individuals (odds ratio, 1.16; 95% CI, 1.06-1.27; P = 1.19 × 10-3). The HLA (rs9272622) associations reported in Aboriginal Australian individuals could not be replicated. In support of the known polygenic architecture for RHD, overtransmission of a polygenic risk score from unaffected parents to affected probands was observed (polygenic transmission disequilibrium testing mean [SE], 0.27 [0.16] SDs; P = .04996), and the chip-based heritability was estimated to be high at 0.49 (SE = 0.12; P = 3.28 × 10-5) in Black African individuals.

Conclusions and Relevance

This study revealed a novel candidate susceptibility locus exclusive to Black African individuals and an important heritable component to RHD susceptibility in African individuals.

Rationale and design of the African Cardiomyopathy and Myocarditis Registry Program: The IMHOTEP study

BACKGROUND

Heart failure (HF), the dominant form of cardiovascular disease in Africans, is mainly due to hypertension, rheumatic heart disease and cardiomyopathy. Cardiomyopathies pose a great challenge because of poor prognosis and high prevalence in low- and middle-income countries (LMICs). Little is known about the etiology and outcome of cardiomyopathy in Africa. Specifically, the role of myocarditis and the genetic causes of cardiomyopathy are largely unidentified in Africans.

METHOD

The African Cardiomyopathy and Myocarditis Registry Program (the IMHOTEP study) is a pan-African multi-centre, hospital-based cohort study, designed with the primary aim of describing the clinical characteristics, genetic causes, prevalence, management and outcome of cardiomyopathy and myocarditis in children and adults. The secondary aim is to identify barriers to the implementation of evidence-based care and provide a platform for trials and other intervention studies to reduce morbidity and mortality in cardiomyopathy. The registry consists of a prospective cohort of newly diagnosed (i.e., incident) cases and a retrospective (i.e., prevalent) cohort of existing cases from participating centres. Patients with cardiomyopathy and myocarditis will be subjected to a standardized 3-stage diagnostic process. To date, 750 patients have been recruited into the multi-centre pilot phase of the study.

CONCLUSION

The IMHOTEP study will provide comprehensive and novel data on clinical features, genetic causes, prevalence and outcome of African children and adults with all forms of cardiomyopathy and myocarditis in Africa. Based on these findings, appropriate strategies for management and prevention of the cardiomyopathies in LMICs are likely to emerge.

X Chromosome Contribution to the Genetic Architecture of Primary Biliary Cholangitis

BACKGROUND & AIMS

Genome-wide association studies in primary biliary cholangitis (PBC) have failed to find X chromosome (chrX) variants associated with the disease. Here, we specifically explore the chrX contribution to PBC, a sexually dimorphic complex autoimmune disease.

METHODS

We performed a chrX-wide association study, including genotype data from 5 genome-wide association studies (from Italy, United Kingdom, Canada, China, and Japan; 5244 case patients and 11,875 control individuals).

RESULTS

Single-marker association analyses found approximately 100 loci displaying P < 5 × 10-4, with the most significant being a signal within the OTUD5 gene (rs3027490; P = 4.80 × 10-6; odds ratio [OR], 1.39; 95% confidence interval [CI], 1.028-1.88; Japanese cohort). Although the transethnic meta-analysis evidenced only a suggestive signal (rs2239452, mapping within the PIM2 gene; OR, 1.17; 95% CI, 1.09-1.26; P = 9.93 × 10-8), the population-specific meta-analysis showed a genome-wide significant locus in East Asian individuals pointing to the same region (rs7059064, mapping within the GRIPAP1 gene; P = 6.2 × 10-9; OR, 1.33; 95% CI, 1.21-1.46). Indeed, rs7059064 tags a unique linkage disequilibrium block including 7 genes: TIMM17B, PQBP1, PIM2, SLC35A2, OTUD5, KCND1, and GRIPAP1, as well as a superenhancer (GH0XJ048933 within OTUD5) targeting all these genes. GH0XJ048933 is also predicted to target FOXP3, the main T-regulatory cell lineage specification factor. Consistently, OTUD5 and FOXP3 RNA levels were up-regulated in PBC case patients (1.75- and 1.64-fold, respectively).

CONCLUSIONS

This work represents the first comprehensive study, to our knowledge, of the chrX contribution to the genetics of an autoimmune liver disease and shows a novel PBC-related genome-wide significant locus.

A Genome-wide Association Study Identifies SERPINB10, CRLF3, STX7, LAMP3, IFNG-AS1, and KRT80 As Risk Loci Contributing to Cutaneous Leishmaniasis in Brazil

BACKGROUND

Our goal was to identify genetic risk factors for cutaneous leishmaniasis (CL) caused by Leishmania braziliensis.

METHODS

Genotyping 2066 CL cases and 2046 controls using Illumina HumanCoreExomeBeadChips provided data for 4 498 586 imputed single-nucleotide variants (SNVs). A genome-wide association study (GWAS) using linear mixed models took account of genetic diversity/ethnicity/admixture. Post-GWAS positional, expression quantitative trait locus (eQTL) and chromatin interaction mapping was performed in Functional Mapping and Annotation (FUMA). Transcriptional data were compared between lesions and normal skin, and cytokines measured using flow cytometry and Bioplex assay.

RESULTS

Positional mapping identified 32 genomic loci associated with CL, none achieving genome-wide significance (P < 5 × 10-8). Lead SNVs at 23 loci occurred at protein coding or noncoding RNA genes, 15 with eQTLs for functionally relevant cells/tissues and/or showing differential expression in lesions. Of these, the 6 most plausible genetic risk loci were SERPINB10 (Pimputed_1000G = 2.67 × 10-6), CRLF3 (Pimputed_1000G = 5.12 × 10-6), STX7 (Pimputed_1000G = 6.06 × 10-6), KRT80 (Pimputed_1000G = 6.58 × 10-6), LAMP3 (Pimputed_1000G = 6.54 × 10-6), and IFNG-AS1 (Pimputed_1000G = 1.32 × 10-5). LAMP3 (Padjusted = 9.25 × 10-12; +6-fold), STX7 (Padjusted = 7.62 × 10-3; +1.3-fold), and CRLF3 (Padjusted = 9.19 × 10-9; +1.97-fold) were expressed more highly in CL biopsies compared to normal skin; KRT80 (Padjusted = 3.07 × 10-8; -3-fold) was lower. Multiple cis-eQTLs across SERPINB10 mapped to chromatin interaction regions of transcriptional/enhancer activity in neutrophils, monocytes, B cells, and hematopoietic stem cells. Those at IFNG-AS1 mapped to transcriptional/enhancer regions in T, natural killer, and B cells. The percentage of peripheral blood CD3+ T cells making antigen-specific interferon-γ differed significantly by IFNG-AS1 genotype.

CONCLUSIONS

This first GWAS for CL identified multiple genetic risk loci including a novel lead to understanding CL pathogenesis through regulation of interferon-γ by IFNG antisense RNA 1.

Genome-wide Association Study and Meta-analysis on Alcohol-Associated Liver Cirrhosis Identifies Genetic Risk Factors

BACKGROUND AND AIMS

Only a minority of heavy drinkers progress to alcohol-associated cirrhosis (ALC). The aim of this study was to identify common genetic variants that underlie risk for ALC.

APPROACH AND RESULTS

We analyzed data from 1,128 subjects of European ancestry with ALC and 614 heavy-drinking subjects without known liver disease from Australia, the United States, the United Kingdom, and three countries in Europe. A genome-wide association study (GWAS) was performed, adjusting for principal components and clinical covariates (alcohol use, age, sex, body mass index, and diabetes). We validated our GWAS findings using UK Biobank. We then performed a meta-analysis combining data from our study, the UK Biobank, and a previously published GWAS. Our GWAS found genome-wide significant risk association of rs738409 in patatin-like phospholipase domain containing 3 (PNPLA3) (odds ratio [OR] = 2.19 [G allele], P = 4.93 × 10^-17 ) and rs4607179 near HSD17B13 (OR = 0.57 [C allele], P = 1.09 × 10^-10 ) with ALC. Conditional analysis accounting for the PNPLA3 and HSD17B13 loci identified a protective association at rs374702773 in Fas-associated factor family member 2 (FAF2) (OR = 0.61 [del(T) allele], P = 2.56 × 10^-8 ) for ALC. This association was replicated in the UK Biobank using conditional analysis (OR = 0.79, P = 0.001). Meta-analysis (without conditioning) confirmed genome-wide significance for the identified FAF2 locus as well as PNPLA3 and HSD17B13. Two other previously known loci (SERPINA1 and SUGP1/TM6SF2) were also genome-wide significant in the meta-analysis. GeneOntology pathway analysis identified lipid droplets as the target for several identified genes. In conclusion, our GWAS identified a locus at FAF2 associated with reduced risk of ALC among heavy drinkers. Like the PNPLA3 and HSD17B13 gene products, the FAF2 product has been localized to fat droplets in hepatocytes.

CONCLUSIONS

Our genetic findings implicate lipid droplets in the biological pathway(s) underlying ALC.

Heritability and family-based GWAS analyses of the N-acyl ethanolamine and ceramide plasma lipidome

Signalling lipids of the N-acyl ethanolamine (NAE) and ceramide (CER) classes have emerged as potential biomarkers of cardiovascular disease (CVD). We sought to establish the heritability of plasma NAEs (including the endocannabinoid anandamide) and CERs, to identify common DNA variants influencing the circulating concentrations of the heritable lipids, and assess causality of these lipids in CVD using 2-sample Mendelian randomization (2SMR). Nine NAEs and 16 CERs were analyzed in plasma samples from 999 members of 196 British Caucasian families, using targeted ultra-performance liquid chromatography with tandem mass spectrometry. All lipids were significantly heritable (h2 = 36-62%). A missense variant (rs324420) in the gene encoding the enzyme fatty acid amide hydrolase (FAAH), which degrades NAEs, associated at genome-wide association study (GWAS) significance (P < 5 × 10-8) with four NAEs (DHEA, PEA, LEA and VEA). For CERs, rs680379 in the SPTLC3 gene, which encodes a subunit of the rate-limiting enzyme in CER biosynthesis, associated with a range of species (e.g. CER[N(24)S(19)]; P = 4.82 × 10-27). We observed three novel associations between SNPs at the CD83, SGPP1 and DEGS1 loci, and plasma CER traits (P < 5 × 10-8). 2SMR in the CARDIoGRAMplusC4D cohorts (60 801 cases; 123 504 controls) and in the DIAGRAM cohort (26 488 cases; 83 964 controls), using the genetic instruments from our family-based GWAS, did not reveal association between genetically determined differences in CER levels and CVD or diabetes. Two of the novel GWAS loci, SGPP1 and DEGS1, suggested a casual association between CERs and a range of haematological phenotypes, through 2SMR in the UK Biobank, INTERVAL and UKBiLEVE cohorts (n = 110 000-350 000).

Publisher Correction: A meta-analysis of genome-wide association studies identifies multiple longevity genes

A Correction to this paper has been published: https://doi.org/10.1038/s41467-021-22613-2

A PDCD1 Role in the Genetic Predisposition to NAFLD-HCC?

Obesity and non-alcoholic fatty liver disease (NAFLD) are contributing to the global rise in deaths from hepatocellular carcinoma (HCC). The pathogenesis of NAFLD-HCC is not well understood. The severity of hepatic steatosis, steatohepatitis and fibrosis are key pathogenic mechanisms, but animal studies suggest altered immune responses are also involved. Genetic studies have so far highlighted a major role of gene variants promoting fat deposition in the liver (PNPLA3 rs738409; TM6SF2 rs58542926). Here, we have considered single-nucleotide polymorphisms (SNPs) in candidate immunoregulatory genes (MICA rs2596542; CD44 rs187115; PDCD1 rs7421861 and rs10204525), in 594 patients with NAFLD and 391 with NAFLD-HCC, from three European centres. Associations between age, body mass index, diabetes, cirrhosis and SNPs with HCC development were explored. PNPLA3 and TM6SF2 SNPs were associated with both progression to cirrhosis and NAFLD-HCC development, while PDCD1 SNPs were specifically associated with NAFLD-HCC risk, regardless of cirrhosis. PDCD1 rs7421861 was independently associated with NAFLD-HCC development, while PDCD1 rs10204525 acquired significance after adjusting for other risks, being most notable in the smaller numbers of women with NAFLD-HCC. The study highlights the potential impact of inter individual variation in immune tolerance induction in patients with NAFLD, both in the presence and absence of cirrhosis.

Congenital heart disease risk loci identified by genome-wide association study in European patients

Genetic factors undoubtedly affect the development of congenital heart disease (CHD) but still remain ill defined. We sought to identify genetic risk factors associated with CHD and to accomplish a functional analysis of SNP-carrying genes. We performed a genome-wide association study (GWAS) of 4034 White patients with CHD and 8486 healthy controls. One SNP on chromosome 5q22.2 reached genome-wide significance across all CHD phenotypes and was also indicative for septal defects. One region on chromosome 20p12.1 pointing to the MACROD2 locus identified 4 highly significant SNPs in patients with transposition of the great arteries (TGA). Three highly significant risk variants on chromosome 17q21.32 within the GOSR2 locus were detected in patients with anomalies of thoracic arteries and veins (ATAV). Genetic variants associated with ATAV are suggested to influence the expression of WNT3, and the variant rs870142 related to septal defects is proposed to influence the expression of MSX1. We analyzed the expression of all 4 genes during cardiac differentiation of human and murine induced pluripotent stem cells in vitro and by single-cell RNA-Seq analyses of developing murine and human hearts. Our data show that MACROD2, GOSR2, WNT3, and MSX1 play an essential functional role in heart development at the embryonic and newborn stages.

Transcriptomic profiling across the nonalcoholic fatty liver disease spectrum reveals gene signatures for steatohepatitis and fibrosis

The mechanisms that drive nonalcoholic fatty liver disease (NAFLD) remain incompletely understood. This large multicenter study characterized the transcriptional changes that occur in liver tissue across the NAFLD spectrum as disease progresses to cirrhosis to identify potential circulating markers. We performed high-throughput RNA sequencing on a discovery cohort comprising histologically characterized NAFLD samples from 206 patients. Unsupervised clustering stratified NAFLD on the basis of disease activity and fibrosis stage with differences in age, aspartate aminotransferase (AST), type 2 diabetes mellitus, and carriage of PNPLA3 rs738409, a genetic variant associated with NAFLD. Relative to early disease, we consistently identified 25 differentially expressed genes as fibrosing steatohepatitis progressed through stages F2 to F4. This 25-gene signature was independently validated by logistic modeling in a separate replication cohort (n = 175), and an integrative analysis with publicly available single-cell RNA sequencing data elucidated the likely relative contribution of specific intrahepatic cell populations. Translating these findings to the protein level, SomaScan analysis in more than 300 NAFLD serum samples confirmed that circulating concentrations of proteins AKR1B10 and GDF15 were strongly associated with disease activity and fibrosis stage. Supporting the biological plausibility of these data, in vitro functional studies determined that endoplasmic reticulum stress up-regulated expression of AKR1B10, GDF15, and PDGFA, whereas GDF15 supplementation tempered the inflammatory response in macrophages upon lipid loading and lipopolysaccharide stimulation. This study provides insights into the pathophysiology of progressive fibrosing steatohepatitis, and proof of principle that transcriptomic changes represent potentially tractable and clinically relevant markers of disease progression.

Cardiac Imaging of Aortic Valve Area From 34 287 UK Biobank Participants Reveals Novel Genetic Associations and Shared Genetic Comorbidity With Multiple Disease Phenotypes

Background:

The aortic valve is an important determinant of cardiovascular physiology and anatomic location of common human diseases.

Methods:

From a sample of 34 287 white British ancestry participants, we estimated functional aortic valve area by planimetry from prospectively obtained cardiac magnetic resonance imaging sequences of the aortic valve. Aortic valve area measurements were submitted to genome-wide association testing, followed by polygenic risk scoring and phenome-wide screening, to identify genetic comorbidities.

Results:

A genome-wide association study of aortic valve area in these UK Biobank participants showed 3 significant associations, indexed by rs71190365 (chr13:50764607, DLEU1 , P =1.8×10 −9 ), rs35991305 (chr12:94191968, CRADD , P =3.4×10 −8 ), and chr17:45013271:C:T ( GOSR2 , P =5.6×10 −8 ). Replication on an independent set of 8145 unrelated European ancestry participants showed consistent effect sizes in all 3 loci, although rs35991305 did not meet nominal significance. We constructed a polygenic risk score for aortic valve area, which in a separate cohort of 311 728 individuals without imaging demonstrated that smaller aortic valve area is predictive of increased risk for aortic valve disease (odds ratio, 1.14; P =2.3×10 −6 ). After excluding subjects with a medical diagnosis of aortic valve stenosis (remaining n=308 683 individuals), phenome-wide association of >10 000 traits showed multiple links between the polygenic score for aortic valve disease and key health-related comorbidities involving the cardiovascular system and autoimmune disease. Genetic correlation analysis supports a shared genetic etiology with between aortic valve area and birth weight along with other cardiovascular conditions.

Conclusions:

These results illustrate the use of automated phenotyping of cardiac imaging data from the general population to investigate the genetic etiology of aortic valve disease, perform clinical prediction, and uncover new clinical and genetic correlates of cardiac anatomy.

Identifying genetic factors that contribute to the increased risk of congenital heart defects in infants with Down syndrome

Atrioventricular septal defects (AVSD) are a severe congenital heart defect present in individuals with Down syndrome (DS) at a > 2000-fold increased prevalence compared to the general population. This study aimed to identify risk-associated genes and pathways and to examine a potential polygenic contribution to AVSD in DS. We analyzed a total cohort of 702 individuals with DS with or without AVSD, with genomic data from whole exome sequencing, whole genome sequencing, and/or array-based imputation. We utilized sequence kernel association testing and polygenic risk score (PRS) methods to examine rare and common variants. Our findings suggest that the Notch pathway, particularly NOTCH4, as well as genes involved in the ciliome including CEP290 may play a role in AVSD in DS. These pathways have also been implicated in DS-associated AVSD in prior studies. A polygenic component for AVSD in DS has not been examined previously. Using weights based on the largest genome-wide association study of congenital heart defects available (2594 cases and 5159 controls; all general population samples), we found PRS to be associated with AVSD with odds ratios ranging from 1.2 to 1.3 per standard deviation increase in PRS and corresponding liability r2 values of approximately 1%, suggesting at least a small polygenic contribution to DS-associated AVSD. Future studies with larger sample sizes will improve identification and quantification of genetic contributions to AVSD in DS.

Heritability of haemodynamics in the ascending aorta

Blood flow in the vasculature can be characterised by dimensionless numbers commonly used to define the level of instabilities in the flow, for example the Reynolds number, Re. Haemodynamics play a key role in cardiovascular disease (CVD) progression. Genetic studies have identified mechanosensitive genes with causal roles in CVD. Given that CVD is highly heritable and abnormal blood flow may increase risk, we investigated the heritability of fluid metrics in the ascending aorta calculated using patient-specific data from cardiac magnetic resonance (CMR) imaging. 341 participants from 108 British Caucasian families were phenotyped by CMR and genotyped for 557,124 SNPs. Flow metrics were derived from the CMR images to provide some local information about blood flow in the ascending aorta, based on maximum values at systole at a single location, denoted max, and a 'peak mean' value averaged over the area of the cross section, denoted pm. Heritability was estimated using pedigree-based (QTDT) and SNP-based (GCTA-GREML) methods. Estimates of Reynolds number based on spatially averaged local flow during systole showed substantial heritability ([Formula: see text], [Formula: see text]), while the estimated heritability for Reynolds number calculated using the absolute local maximum velocity was not statistically significant (12-13%; [Formula: see text]). Heritability estimates of the geometric quantities alone; e.g. aortic diameter ([Formula: see text], [Formula: see text]), were also substantially heritable, as described previously. These findings indicate the potential for the discovery of genetic factors influencing haemodynamic traits in large-scale genotyped and phenotyped cohorts where local spatial averaging is used, rather than instantaneous values. Future Mendelian randomisation studies of aortic haemodynamic estimates, which are swift to derive in a clinical setting, will allow for the investigation of causality of abnormal blood flow in CVD.

Investigation of genetically regulated gene expression and response to treatment in rheumatoid arthritis highlights an association between IL18RAP expression and treatment response

Objectives

In this study, we sought to investigate whether there was any association between genetically regulated gene expression (as predicted using various reference panels) and anti-tumour necrosis factor (anti-TNF) treatment response (change in erythrocyte sedimentation rate (ESR)) using 3158 European ancestry patients with rheumatoid arthritis.

Methods

The genetically regulated portion of gene expression was estimated in the full cohort of 3158 subjects (as well as within a subcohort consisting of 1575 UK patients) using the PrediXcan software package with three different reference panels. Estimated expression was tested for association with anti-TNF treatment response. As a replication/validation experiment, we also investigated the correlation between change in ESR with measured gene expression at the Interleukin 18 Receptor Accessory Protein (IL18RAP) gene in whole blood and synovial tissue, using an independent replication data set of patients receiving conventional synthetic disease modifying anti-rheumatic drugs, with directly measured (via RNA sequencing) gene expression.

Results

We found that predicted expression of IL18RAP showed a consistent signal of association with treatment response across the reference panels. In our independent replication data set, IL18RAP expression in whole blood showed correlation with the change in ESR between baseline and follow-up (r=−0.35, p=0.0091). Change in ESR was also correlated with the expression of IL18RAP in synovial tissue (r=−0.28, p=0.02).

Conclusion

Our results suggest that IL18RAP expression is worthy of further investigation as a potential predictor of treatment response in rheumatoid arthritis that is not specific to a particular drug type.

Investigation of prediction accuracy and the impact of sample size, ancestry, and tissue in transcriptome-wide association studies

In transcriptome-wide association studies (TWAS), gene expression values are predicted using genotype data and tested for association with a phenotype. The power of this approach to detect associations relies, at least in part, on the accuracy of the prediction. Here we compare the prediction accuracy of six different methods-LASSO, Ridge regression, Elastic net, Best Linear Unbiased Predictor, Bayesian Sparse Linear Mixed Model, and Random Forests-by performing cross-validation using data from the Geuvadis Project. We also examine prediction accuracy (a) at different sample sizes, (b) when ancestry of the prediction model training and testing populations is different, and (c) when the tissue used to train the model is different from the tissue to be predicted. We find that, for most genes, the expression cannot be accurately predicted, but in general sparse statistical models tend to outperform polygenic models at prediction. Average prediction accuracy is reduced when the model training set size is reduced or when predicting across ancestries and is marginally reduced when predicting across tissues. We conclude that using sparse statistical models and the development of large reference panels across multiple ethnicities and tissues will lead to better prediction of gene expression, and thus may improve TWAS power.

Mutations in SPATA13/ASEF2 cause primary angle closure glaucoma

Current estimates suggest 50% of glaucoma blindness worldwide is caused by primary angle-closure glaucoma (PACG) but the causative gene is not known. We used genetic linkage and whole genome sequencing to identify Spermatogenesis Associated Protein 13, SPATA13 (NM_001166271; NP_001159743, SPATA13 isoform I), also known as ASEF2 (Adenomatous polyposis coli-stimulated guanine nucleotide exchange factor 2), as the causal gene for PACG in a large seven-generation white British family showing variable expression and incomplete penetrance. The 9 bp deletion, c.1432_1440del; p.478_480del was present in all affected individuals with angle-closure disease. We show ubiquitous expression of this transcript in cell lines derived from human tissues and in iris, retina, retinal pigment and ciliary epithelia, cornea and lens. We also identified eight additional mutations in SPATA13 in a cohort of 189 unrelated PACS/PAC/PACG samples. This gene encodes a 1277 residue protein which localises to the nucleus with partial co-localisation with nuclear speckles. In cells undergoing mitosis SPATA13 isoform I becomes part of the kinetochore complex co-localising with two kinetochore markers, polo like kinase 1 (PLK-1) and centrosome-associated protein E (CENP-E). The 9 bp deletion reported in this study increases the RAC1-dependent guanine nucleotide exchange factors (GEF) activity. The increase in GEF activity was also observed in three other variants identified in this study. Taken together, our data suggest that SPATA13 is involved in the regulation of mitosis and the mutations dysregulate GEF activity affecting homeostasis in tissues where it is highly expressed, influencing PACG pathogenesis.

Bayesian network analysis incorporating genetic anchors complements conventional Mendelian randomization approaches for exploratory analysis of causal relationships in complex data

Mendelian randomization (MR) implemented through instrumental variables analysis is an increasingly popular causal inference tool used in genetic epidemiology. But it can have limitations for evaluating simultaneous causal relationships in complex data sets that include, for example, multiple genetic predictors and multiple potential risk factors associated with the same genetic variant. Here we use real and simulated data to investigate Bayesian network analysis (BN) with the incorporation of directed arcs, representing genetic anchors, as an alternative approach. A Bayesian network describes the conditional dependencies/independencies of variables using a graphical model (a directed acyclic graph) with an accompanying joint probability. In real data, we found BN could be used to infer simultaneous causal relationships that confirmed the individual causal relationships suggested by bi-directional MR, while allowing for the existence of potential horizontal pleiotropy (that would violate MR assumptions). In simulated data, BN with two directional anchors (mimicking genetic instruments) had greater power for a fixed type 1 error than bi-directional MR, while BN with a single directional anchor performed better than or as well as bi-directional MR. Both BN and MR could be adversely affected by violations of their underlying assumptions (such as genetic confounding due to unmeasured horizontal pleiotropy). BN with no directional anchor generated inference that was no better than by chance, emphasizing the importance of directional anchors in BN (as in MR). Under highly pleiotropic simulated scenarios, BN outperformed both MR (and its recent extensions) and two recently-proposed alternative approaches: a multi-SNP mediation intersection-union test (SMUT) and a latent causal variable (LCV) test. We conclude that BN incorporating genetic anchors is a useful complementary method to conventional MR for exploring causal relationships in complex data sets such as those generated from modern "omics" technologies.

Whole Exome Sequencing Reveals the Major Genetic Contributors to Nonsyndromic Tetralogy of Fallot

RATIONALE

Familial recurrence studies provide strong evidence for a genetic component to the predisposition to sporadic, nonsyndromic Tetralogy of Fallot (TOF), the most common cyanotic congenital heart disease phenotype. Rare genetic variants have been identified as important contributors to the risk of congenital heart disease, but relatively small numbers of TOF cases have been studied to date.

OBJECTIVE

We used whole exome sequencing to assess the prevalence of unique, deleterious variants in the largest cohort of nonsyndromic TOF patients reported to date.

METHODS AND RESULTS

Eight hundred twenty-nine TOF patients underwent whole exome sequencing. The presence of unique, deleterious variants was determined; defined by their absence in the Genome Aggregation Database and a scaled combined annotation-dependent depletion score of ≥20. The clustering of variants in 2 genes, NOTCH1 and FLT4, surpassed thresholds for genome-wide significance (assigned as P<5×10^-8) after correction for multiple comparisons. NOTCH1 was most frequently found to harbor unique, deleterious variants. Thirty-one changes were observed in 37 probands (4.5%; 95% CI, 3.2%-6.1%) and included 7 loss-of-function variants 22 missense variants and 2 in-frame indels. Sanger sequencing of the unaffected parents of 7 cases identified 5 de novo variants. Three NOTCH1 variants (p.G200R, p.C607Y, and p.N1875S) were subjected to functional evaluation, and 2 showed a reduction in Jagged1-induced NOTCH signaling. FLT4 variants were found in 2.4% (95% CI, 1.6%-3.8%) of TOF patients, with 21 patients harboring 22 unique, deleterious variants. The variants identified were distinct to those that cause the congenital lymphoedema syndrome Milroy disease. In addition to NOTCH1, FLT4 and the well-established TOF gene, TBX1, we identified potential association with variants in several other candidates, including RYR1, ZFPM1, CAMTA2, DLX6, and PCM1.

CONCLUSIONS

The NOTCH1 locus is the most frequent site of genetic variants predisposing to nonsyndromic TOF, followed by FLT4. Together, variants in these genes are found in almost 7% of TOF patients.

Marked variation in heritability estimates of left ventricular mass depending on modality of measurement

Left ventricular (LV) hypertrophy is a strong risk factor for heart failure and cardiovascular death. ECG measures of LV mass are estimated as heritable in twin and family-based analyses and heritability estimates of LV mass measured by echocardiography are lower. We hypothesised that CMR-derived measurements, being more precise than echocardiographic measurements, would advance our understanding of heritable LV traits. We phenotyped 116 British families (427 individuals) by CMR and ECG, and undertook heritability analyses using variance-components (QTDT) and GWAS SNP-based (GCTA-GREML) methods. ECG-based traits such as LV mass and Sokolow-Lyon duration showed substantial estimates of heritability (60%), whereas CMR-derived LV mass was only modestly heritable (20%). However, the ECG LV mass was positively correlated with the lateral diameter of the chest (rho = 0.67), and adjustment for this attenuated the heritability estimate (42%). Finally, CMR-derived right ventricular mass showed considerable heritability (44%). Heritability estimates of LV phenotypes show substantial variation depending on the modality of measurement, being greater when measured by ECG than CMR. This may reflect the differences between electrophysiological as opposed to anatomical hypertrophy. However, ECG LV hypertrophy traits are likely to be influenced by genetic association with anthropometric measures, inflating their overall measured heritability.

Exploring Shared Susceptibility between Two Neural Crest Cells Originating Conditions: Neuroblastoma and Congenital Heart Disease

In the past years, genome wide association studies (GWAS) have provided evidence that inter-individual susceptibility to diverse pathological conditions can reveal a common genetic architecture. Through the analysis of congenital heart disease (CHD) and neuroblastoma (NB) GWAS data, we aimed to dissect the genetic susceptibility shared between these conditions, which are known to arise from neural crest cell (NCC) migration or development abnormalities, via identification and functional characterization of common regions of association. Two loci (2q35 and 3q25.32) harbor single nucleotide polymorphisms (SNPs) that are associated at a p-value < 10^-3 with conotruncal malformations and ventricular septal defect respectively, as well as with NB. In addition, the lead SNP in 4p16.2 for atrial septal defect and the lead SNP in 3q25.32 for tetralogy of Fallot are less than 250 Kb distant from the lead SNPs for NB at the same genomic regions. Some of these shared susceptibility loci regulate the expression of relevant genes involved in NCC formation and developmental processes (such as BARD1, MSX1, and SHOX2) and are enriched in several epigenetic markers from NB and fetal heart cell lines. Although the clinical correlation between NB and CHD is unclear, our exploration of a possible common genetic basis between NB and a subset of cardiac malformations can help shed light on their shared embryological origin and pathogenetic mechanisms.