Genome-Wide Association Study and Gene Expression Analysis Identifies CD84 as a Predictor of Response to Etanercept Therapy in Rheumatoid Arthritis

Anti-tumor necrosis factor alpha (anti-TNF) biologic therapy is a widely used treatment for rheumatoid arthritis (RA). It is unknown why some RA patients fail to respond adequately to anti-TNF therapy, which limits the development of clinical biomarkers to predict response or new drugs to target refractory cases. To understand the biological basis of response to anti-TNF therapy, we conducted a genome-wide association study (GWAS) meta-analysis of more than 2 million common variants in 2,706 RA patients from 13 different collections. Patients were treated with one of three anti-TNF medications: etanercept (n = 733), infliximab (n = 894), or adalimumab (n = 1,071). We identified a SNP (rs6427528) at the 1q23 locus that was associated with change in disease activity score (ΔDAS) in the etanercept subset of patients (P = 8×10-8), but not in the infliximab or adalimumab subsets (P>0.05). The SNP is predicted to disrupt transcription factor binding site motifs in the 3′ UTR of an immune-related gene, CD84, and the allele associated with better response to etanercept was associated with higher CD84 gene expression in peripheral blood mononuclear cells (P = 1×10-11 in 228 non-RA patients and P = 0.004 in 132 RA patients). Consistent with the genetic findings, higher CD84 gene expression correlated with lower cross-sectional DAS (P = 0.02, n = 210) and showed a non-significant trend for better ΔDAS in a subset of RA patients with gene expression data (n = 31, etanercept-treated). A small, multi-ethnic replication showed a non-significant trend towards an association among etanercept-treated RA patients of Portuguese ancestry (n = 139, P = 0.4), but no association among patients of Japanese ancestry (n = 151, P = 0.8). Our study demonstrates that an allele associated with response to etanercept therapy is also associated with CD84 gene expression, and further that CD84 expression correlates with disease activity. These findings support a model in which CD84 genotypes and/or expression may serve as a useful biomarker for response to etanercept treatment in RA patients of European ancestry. © 2013 Cui et al

Integration of sequence data from a consanguineous family with genetic data from an outbred population identifies PLB1 as a candidate rheumatoid arthritis risk gene

Integrating genetic data from families with highly penetrant forms of disease together with genetic data from outbred populations represents a promising strategy to uncover the complete frequency spectrum of risk alleles for complex traits such as rheumatoid arthritis (RA). Here, we demonstrate that rare, low-frequency and common alleles at one gene locus, phospholipase B1 (PLB1), might contribute to risk of RA in a 4-generation consanguineous pedigree (Middle Eastern ancestry) and also in unrelated individuals from the general population (European ancestry). Through identity-by-descent (IBD) mapping and whole-exome sequencing, we identified a non-synonymous c.2263G>C (p.G755R) mutation at the PLB1 gene on 2q23, which significantly co-segregated with RA in family members with a dominant mode of inheritance (P = 0.009). We further evaluated PLB1 variants and risk of RA using a GWAS meta-analysis of 8,875 RA cases and 29,367 controls of European ancestry. We identified significant contributions of two independent non-coding variants near PLB1 with risk of RA (rs116018341 [MAF = 0.042] and rs116541814 [MAF = 0.021], combined P = 3.2×10-6). Finally, we performed deep exon sequencing of PLB1 in 1,088 RA cases and 1,088 controls (European ancestry), and identified suggestive dispersion of rare protein-coding variant frequencies between cases and controls (P = 0.049 for C-alpha test and P = 0.055 for SKAT). Together, these data suggest that PLB1 is a candidate risk gene for RA. Future studies to characterize the full spectrum of genetic risk in the PLB1 genetic locus are warranted. © 2014 Plenge et al

Integration of Sequence Data from a Consanguineous Family with Genetic Data from an Outbred Population Identifies PLB1 as a Candidate Rheumatoid Arthritis Risk Gene

Integrating genetic data from families with highly penetrant forms of disease together with genetic data from outbred populations represents a promising strategy to uncover the complete frequency spectrum of risk alleles for complex traits such as rheumatoid arthritis (RA). Here, we demonstrate that rare, low-frequency and common alleles at one gene locus, phospholipase B1 (PLB1), might contribute to risk of RA in a 4-generation consanguineous pedigree (Middle Eastern ancestry) and also in unrelated individuals from the general population (European ancestry). Through identity-by-descent (IBD) mapping and whole-exome sequencing, we identified a non-synonymous c.2263G>C (p.G755R) mutation at the PLB1 gene on 2q23, which significantly co-segregated with RA in family members with a dominant mode of inheritance (P = 0.009). We further evaluated PLB1 variants and risk of RA using a GWAS meta-analysis of 8,875 RA cases and 29,367 controls of European ancestry. We identified significant contributions of two independent non-coding variants near PLB1 with risk of RA (rs116018341 [MAF = 0.042] and rs116541814 [MAF = 0.021], combined P = 3.2×10−6). Finally, we performed deep exon sequencing of PLB1 in 1,088 RA cases and 1,088 controls (European ancestry), and identified suggestive dispersion of rare protein-coding variant frequencies between cases and controls (P = 0.049 for C-alpha test and P = 0.055 for SKAT). Together, these data suggest that PLB1 is a candidate risk gene for RA. Future studies to characterize the full spectrum of genetic risk in the PLB1 genetic locus are warranted

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

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

Chromosome Xq23 is associated with lower atherogenic lipid concentrations and favorable cardiometabolic indices

AbstractAutosomal genetic analyses of blood lipids have yielded key insights for coronary heart disease (CHD). However, X chromosome genetic variation is understudied for blood lipids in large sample sizes. We now analyze genetic and blood lipid data in a high-coverage whole X chromosome sequencing study of 65,322 multi-ancestry participants and perform replication among 456,893 European participants. Common alleles on chromosome Xq23 are strongly associated with reduced total cholesterol, LDL cholesterol, and triglycerides (min P = 8.5 × 10−72), with similar effects for males and females. Chromosome Xq23 lipid-lowering alleles are associated with reduced odds for CHD among 42,545 cases and 591,247 controls (P = 1.7 × 10−4), and reduced odds for diabetes mellitus type 2 among 54,095 cases and 573,885 controls (P = 1.4 × 10−5). Although we observe an association with increased BMI, waist-to-hip ratio adjusted for BMI is reduced, bioimpedance analyses indicate increased gluteofemoral fat, and abdominal MRI analyses indicate reduced visceral adiposity. Co-localization analyses strongly correlate increased CHRDL1 gene expression, particularly in adipose tissue, with reduced concentrations of blood lipids.Abstract Autosomal genetic analyses of blood lipids have yielded key insights for coronary heart disease (CHD). However, X chromosome genetic variation is understudied for blood lipids in large sample sizes. We now analyze genetic and blood lipid data in a high-coverage whole X chromosome sequencing study of 65,322 multi-ancestry participants and perform replication among 456,893 European participants. Common alleles on chromosome Xq23 are strongly associated with reduced total cholesterol, LDL cholesterol, and triglycerides (min P = 8.5 × 10−72), with similar effects for males and females. Chromosome Xq23 lipid-lowering alleles are associated with reduced odds for CHD among 42,545 cases and 591,247 controls (P = 1.7 × 10−4), and reduced odds for diabetes mellitus type 2 among 54,095 cases and 573,885 controls (P = 1.4 × 10−5). Although we observe an association with increased BMI, waist-to-hip ratio adjusted for BMI is reduced, bioimpedance analyses indicate increased gluteofemoral fat, and abdominal MRI analyses indicate reduced visceral adiposity. Co-localization analyses strongly correlate increased CHRDL1 gene expression, particularly in adipose tissue, with reduced concentrations of blood lipids

The genomic basis of parasitism in the Strongyloides clade of nematodes.

Soil-transmitted nematodes, including the Strongyloides genus, cause one of the most prevalent neglected tropical diseases. Here we compare the genomes of four Strongyloides species, including the human pathogen Strongyloides stercoralis, and their close relatives that are facultatively parasitic (Parastrongyloides trichosuri) and free-living (Rhabditophanes sp. KR3021). A significant paralogous expansion of key gene families--families encoding astacin-like and SCP/TAPS proteins--is associated with the evolution of parasitism in this clade. Exploiting the unique Strongyloides life cycle, we compare the transcriptomes of the parasitic and free-living stages and find that these same gene families are upregulated in the parasitic stages, underscoring their role in nematode parasitism

A phase 1/2 safety and efficacy study of TAK-754 gene therapy: The challenge of achieving durable factor VIII expression in haemophilia A clinical trials

Takeda does not plan to share data supporting the results reported in this article because there is a reasonable likelihood that individual patients could be re-identified due to the very limited (small) study population. Data from UK biobank cannot be shared openly and are subject to UK Biobank ethical approval. Further information about applying for UK Biobank data access can be obtained from the UK Biobank website (https://www.ukbiobank.ac.uk) or by emailing UK Biobank ([email protected]).[Introduction] Haemophilia A is an X-linked bleeding disorder resulting from a deficiency of factor VIII (FVIII). To date, multiple gene therapies have entered clinical trials with the goal of providing durable haemostatic protection from a single dose. TAK 754 (BAX 888) is an investigational AAV8-based gene therapy containing a FVIII transgene. Reduction in CpG motifs was performed to reduce immunogenicity based on prior observations. Here, we describe the results of the first two cohorts treated with TAK 754.[Aim] To report clinical and translational results of the TAK-754 phase 1/2 AAV gene therapy study for the treatment of haemophilia A.[Methods] A phase 1/2 single arm open-label dose escalation study of TAK-754 was performed in participants with severe haemophilia A (NCT03370172). Participants were monitored for safety events, endogenous FVIII activity and bleeding rates. Glucocorticoids were implemented to preserve transgene expression. A transcriptomics analysis was performed to evaluate immunogenicity along with additional post-hoc analyses.[Results] Four participants were dosed in two cohorts. Infusion of TAK 754 was well-tolerated. All participants developed mild transient transaminase elevation and subsequent loss of FVIII expression within the first 12 months of treatment despite use of glucocorticoids. Transcriptomic analysis did not demonstrate significant changes in immunogenicity signals in peripheral blood. One serious adverse event of hypophosphatemia occurred in the second cohort without obvious risk factors.[Conclusions] Sustained FVIII expression remains a challenge in haemophilia A AAV gene therapy trials. Mechanisms of transgene expression loss require further study as clinical studies enter long term follow-up periods.We acknowledge the study participants, their families, and the investigators and study teams who participated in this clinical trial. We express our gratitude to Erin Smith for her internal curation of UKBB data. This study used the UK Biobank Resource under Application Number 26041. We acknowledge the Shire Vienna team for their contributions to the development of BAX-888/TAK-754. This research is funded by Takeda Pharmaceuticals, Ltd.Peer reviewe

Using genotype data to distinguish pleiotropy from heterogeneity: deciphering coheritability in autoimmune and neuropsychiatric diseases

Shared genetic architecture between phenotypes may be driven by a common genetic basis (pleiotropy) or a subset of genetically similar individuals (heterogeneity). We developed BUHMBOX, a well-powered statistical method to distinguish pleiotropy from heterogeneity using genotype data. We observed a shared genetic basis between 11 of 17 tested autoimmune diseases and type I diabetes (T1D, p&lt;10 12) and 11 of 17 tested autoimmune diseases and rheumatoid arthritis (RA, p&lt;10-7). This sharing could not be explained by heterogeneity (corrected pBUHMBOX&gt;0.2 using 6,670 T1D cases and 7,279 RA cases), suggesting that shared genetic features in autoimmunity are due to pleiotropy. We observed a shared genetic basis between seronegative and seropostive RA (p&lt;10-22), explained by heterogeneity (pBUHMBOX=0.008 in 2,406 seronegative RA cases). Consistent with previous observations, we observed genetic sharing between major depressive disorder (MDD) and schizophrenia (p&lt;10 9). This sharing is not explained by heterogeneity (pBUHMBOX=0.28 in 9,238 MDD cases).</jats:p

Brief Report:The Role of Rare Protein-Coding Variants in Anti–Tumor Necrosis Factor Treatment Response in Rheumatoid Arthritis

In many rheumatoid arthritis (RA) patients, disease is controlled with anti-tumor necrosis factor (anti-TNF) biologic therapies. However, in a significant number of patients, the disease fails to respond to anti-TNF therapy. We undertook the present study to examine the hypothesis that rare and low-frequency genetic variants might influence response to anti-TNF treatment. We sequenced the coding region of 750 genes in 1,094 RA patients of European ancestry who were treated with anti-TNF. After quality control, 690 genes were included in the analysis. We applied single-variant association and gene-based association tests to identify variants associated with anti-TNF treatment response. In addition, given the key mechanistic role of TNF, we performed gene set analyses of 27 TNF pathway genes. We identified 14,420 functional variants, of which 6,934 were predicted as nonsynonymous 2,136 of which were further predicted to be "damaging." Despite the fact that the study was well powered, no single variant or gene showed study-wide significant association with change in the outcome measures disease activity or European League Against Rheumatism response. Intriguingly, we observed 3 genes, of 27 with nominal signals of association (P  < 0.05), that were involved in the TNF signaling pathway. However, when we performed a rigorous gene set enrichment analysis based on association P value ranking, we observed no evidence of enrichment of association at genes involved in the TNF pathway (Penrichment  = 0.15, based on phenotype permutations). Our findings suggest that rare and low-frequency protein-coding variants in TNF signaling pathway genes or other genes do not contribute substantially to anti-TNF treatment response in patients with R