Rare variant associations with plasma protein levels in the UK Biobank.

Integrating human genomics and proteomics can help elucidate disease mechanisms, identify clinical biomarkers and discover drug targets1-4. Because previous proteogenomic studies have focused on common variation via genome-wide association studies, the contribution of rare variants to the plasma proteome remains largely unknown. Here we identify associations between rare protein-coding variants and 2,923 plasma protein abundances measured in 49,736 UK Biobank individuals. Our variant-level exome-wide association study identified 5,433 rare genotype-protein associations, of which 81% were undetected in a previous genome-wide association study of the same cohort5. We then looked at aggregate signals using gene-level collapsing analysis, which revealed 1,962 gene-protein associations. Of the 691 gene-level signals from protein-truncating variants, 99.4% were associated with decreased protein levels. STAB1 and STAB2, encoding scavenger receptors involved in plasma protein clearance, emerged as pleiotropic loci, with 77 and 41 protein associations, respectively. We demonstrate the utility of our publicly accessible resource through several applications. These include detailing an allelic series in NLRC4, identifying potential biomarkers for a fatty liver disease-associated variant in HSD17B13 and bolstering phenome-wide association studies by integrating protein quantitative trait loci with protein-truncating variants in collapsing analyses. Finally, we uncover distinct proteomic consequences of clonal haematopoiesis (CH), including an association between TET2-CH and increased FLT3 levels. Our results highlight a considerable role for rare variation in plasma protein abundance and the value of proteogenomics in therapeutic discovery.

Misexpression of inactive genes in whole blood is associated with nearby rare structural variants.

Gene misexpression is the aberrant transcription of a gene in a context where it is usually inactive. Despite its known pathological consequences in specific rare diseases, we have a limited understanding of its wider prevalence and mechanisms in humans. To address this, we analyzed gene misexpression in 4,568 whole-blood bulk RNA sequencing samples from INTERVAL study blood donors. We found that while individual misexpression events occur rarely, in aggregate they were found in almost all samples and a third of inactive protein-coding genes. Using 2,821 paired whole-genome and RNA sequencing samples, we identified that misexpression events are enriched in cis for rare structural variants. We established putative mechanisms through which a subset of SVs lead to gene misexpression, including transcriptional readthrough, transcript fusions, and gene inversion. Overall, we develop misexpression as a type of transcriptomic outlier analysis and extend our understanding of the variety of mechanisms by which genetic variants can influence gene expression.

Rare variant contribution to human disease in 281,104 UK Biobank exomes.

Genome-wide association studies have uncovered thousands of common variants associated with human disease, but the contribution of rare variants to common disease remains relatively unexplored. The UK Biobank contains detailed phenotypic data linked to medical records for approximately 500,000 participants, offering an unprecedented opportunity to evaluate the effect of rare variation on a broad collection of traits1,2. Here we study the relationships between rare protein-coding variants and 17,361 binary and 1,419 quantitative phenotypes using exome sequencing data from 269,171 UK Biobank participants of European ancestry. Gene-based collapsing analyses revealed 1,703 statistically significant gene-phenotype associations for binary traits, with a median odds ratio of 12.4. Furthermore, 83% of these associations were undetectable via single-variant association tests, emphasizing the power of gene-based collapsing analysis in the setting of high allelic heterogeneity. Gene-phenotype associations were also significantly enriched for loss-of-function-mediated traits and approved drug targets. Finally, we performed ancestry-specific and pan-ancestry collapsing analyses using exome sequencing data from 11,933 UK Biobank participants of African, East Asian or South Asian ancestry. Our results highlight a significant contribution of rare variants to common disease. Summary statistics are publicly available through an interactive portal ( http://azphewas.com/ ).

Deep Learning-based Segmentation of Computed Tomography Scans Predicts Disease Progression and Mortality in Idiopathic Pulmonary Fibrosis.

Rationale: Despite evidence demonstrating a prognostic role for computed tomography (CT) scans in idiopathic pulmonary fibrosis (IPF), image-based biomarkers are not routinely used in clinical practice or trials. Objectives: To develop automated imaging biomarkers using deep learning-based segmentation of CT scans. Methods: We developed segmentation processes for four anatomical biomarkers, which were applied to a unique cohort of treatment-naive patients with IPF enrolled in the PROFILE (Prospective Observation of Fibrosis in the Lung Clinical Endpoints) study and tested against a further United Kingdom cohort. The relationships among CT biomarkers, lung function, disease progression, and mortality were assessed. Measurements and Main Results: Data from 446 PROFILE patients were analyzed. Median follow-up duration was 39.1 months (interquartile range, 18.1-66.4 mo), with a cumulative incidence of death of 277 (62.1%) over 5 years. Segmentation was successful on 97.8% of all scans, across multiple imaging vendors, at slice thicknesses of 0.5-5 mm. Of four segmentations, lung volume showed the strongest correlation with FVC (r = 0.82; P < 0.001). Lung, vascular, and fibrosis volumes were consistently associated across cohorts with differential 5-year survival, which persisted after adjustment for baseline gender, age, and physiology score. Lower lung volume (hazard ratio [HR], 0.98 [95% confidence interval (CI), 0.96-0.99]; P = 0.001), increased vascular volume (HR, 1.30 [95% CI, 1.12-1.51]; P = 0.001), and increased fibrosis volume (HR, 1.17 [95% CI, 1.12-1.22]; P < 0.001) were associated with reduced 2-year progression-free survival in the pooled PROFILE cohort. Longitudinally, decreasing lung volume (HR, 3.41 [95% CI, 1.36-8.54]; P = 0.009) and increasing fibrosis volume (HR, 2.23 [95% CI, 1.22-4.08]; P = 0.009) were associated with differential survival. Conclusions: Automated models can rapidly segment IPF CT scans, providing prognostic near and long-term information, which could be used in routine clinical practice or as key trial endpoints.

MiCellAnnGELo: annotate microscopy time series of complex cell surfaces with 3D virtual reality.

SUMMARY: Advances in 3D live cell microscopy are enabling high-resolution capture of previously unobserved processes. Unleashing the power of modern machine learning methods to fully benefit from these technologies is, however, frustrated by the difficulty of manually annotating 3D training data. MiCellAnnGELo virtual reality software offers an immersive environment for viewing and interacting with 4D microscopy data, including efficient tools for annotation. We present tools for labelling cell surfaces with a wide range of applications, including cell motility, endocytosis and transmembrane signalling. AVAILABILITY AND IMPLEMENTATION: MiCellAnnGELo employs the cross-platform (Mac/Unix/Windows) Unity game engine and is available under the MIT licence at https://github.com/CellDynamics/MiCellAnnGELo.git, together with sample data. MiCellAnnGELo can be run in desktop mode on a 2D screen or in 3D using a standard VR headset with a compatible GPU. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Type I interferon blockade with anifrolumab in patients with systemic lupus erythematosus modulates key immunopathological pathways in a gene expression and proteomic analysis of two phase 3 trials.

INTRODUCTION: Anifrolumab is a type I interferon (IFN) receptor 1 (IFNAR1) blocking antibody approved for treating patients with systemic lupus erythematosus (SLE). Here, we investigated the immunomodulatory mechanisms of anifrolumab using longitudinal transcriptomic and proteomic analyses of the 52-week, randomised, phase 3 TULIP-1 and TULIP-2 trials. METHODS: Patients with moderate to severe SLE were enrolled in TULIP-1 and TULIP-2 and received intravenous anifrolumab or placebo alongside standard therapy. Whole-blood expression of 18 017 genes using genome-wide RNA sequencing (RNA-seq) (pooled TULIP; anifrolumab, n=244; placebo, n=258) and 184 plasma proteins using Olink and Simoa panels (TULIP-1; anifrolumab, n=124; placebo, n=132) were analysed. We compared treatment groups via gene set enrichment analysis using MetaBase pathway analysis, blood transcriptome modules, in silico deconvolution of RNA-seq and longitudinal linear mixed effect models for gene counts and protein levels. RESULTS: Compared with placebo, anifrolumab modulated >2000 genes by week 24, with overlapping results at week 52, and 41 proteins by week 52. IFNAR1 blockade with anifrolumab downregulated multiple type I and II IFN-induced gene modules/pathways and type III IFN-λ protein levels, and impacted apoptosis-associated and neutrophil extracellular traps-(NET)osis-associated transcriptional pathways, innate cell activating chemokines and receptors, proinflammatory cytokines and B-cell activating cytokines. In silico deconvolution of RNA-seq data indicated an increase from baseline of mucosal-associated invariant and γδT cells and a decrease of monocytes following anifrolumab treatment. DISCUSSION: Type I IFN blockade with anifrolumab modulated multiple inflammatory pathways downstream of type I IFN signalling, including apoptotic, innate and adaptive mechanisms that play key roles in SLE immunopathogenesis.

Rare and Common Variants in KIF15 Contribute to Genetic Risk of Idiopathic Pulmonary Fibrosis.

Rationale: Genetic studies of idiopathic pulmonary fibrosis (IPF) have improved our understanding of this disease, but not all causal loci have been identified. Objectives: To identify genes enriched with rare deleterious variants in IPF and familial pulmonary fibrosis. Methods: We performed gene burden analysis of whole-exome data, tested single variants for disease association, conducted KIF15 (kinesin family member 15) functional studies, and examined human lung single-cell RNA sequencing data. Measurements and Main Results: Gene burden analysis of 1,725 cases and 23,509 control subjects identified heterozygous rare deleterious variants in KIF15, a kinesin involved in spindle separation during mitosis, and three telomere-related genes (TERT [telomerase reverse transcriptase], RTEL1 [regulator of telomere elongation helicase 1], and PARN [poly(A)-specific ribonuclease]). KIF15 was implicated in autosomal-dominant models of rare deleterious variants (odds ratio [OR], 4.9; 95% confidence interval [CI], 2.7-8.8; P = 2.55 × 10-7) and rare protein-truncating variants (OR, 7.6; 95% CI, 3.3-17.1; P = 8.12 × 10-7). Meta-analyses of the discovery and replication cohorts, including 2,966 cases and 29,817 control subjects, confirm the involvement of KIF15 plus the three telomere-related genes. A common variant within a KIF15 intron (rs74341405; OR, 1.6; 95% CI, 1.4-1.9; P = 5.63 × 10-10) is associated with IPF risk, confirming a prior report. Lymphoblastoid cells from individuals heterozygous for the common variant have decreased KIF15 and reduced rates of cell growth. Cell proliferation is dependent on KIF15 in the presence of an inhibitor of Eg5/KIF11, which has partially redundant function. KIF15 is expressed specifically in replicating human lung cells and shows diminished expression in replicating epithelial cells of patients with IPF. Conclusions: Both rare deleterious variants and common variants in KIF15 link a nontelomerase pathway of cell proliferation with IPF susceptibility.

Diagnostic Utility of Exome Sequencing for Kidney Disease.

BACKGROUND: Exome sequencing is emerging as a first-line diagnostic method in some clinical disciplines, but its usefulness has yet to be examined for most constitutional disorders in adults, including chronic kidney disease, which affects more than 1 in 10 persons globally. METHODS: We conducted exome sequencing and diagnostic analysis in two cohorts totaling 3315 patients with chronic kidney disease. We assessed the diagnostic yield and, among the patients for whom detailed clinical data were available, the clinical implications of diagnostic and other medically relevant findings. RESULTS: In all, 3037 patients (91.6%) were over 21 years of age, and 1179 (35.6%) were of self-identified non-European ancestry. We detected diagnostic variants in 307 of the 3315 patients (9.3%), encompassing 66 different monogenic disorders. Of the disorders detected, 39 (59%) were found in only a single patient. Diagnostic variants were detected across all clinically defined categories, including congenital or cystic renal disease (127 of 531 patients [23.9%]) and nephropathy of unknown origin (48 of 281 patients [17.1%]). Of the 2187 patients assessed, 34 (1.6%) had genetic findings for medically actionable disorders that, although unrelated to their nephropathy, would also lead to subspecialty referral and inform renal management. CONCLUSIONS: Exome sequencing in a combined cohort of more than 3000 patients with chronic kidney disease yielded a genetic diagnosis in just under 10% of cases. (Funded by the National Institutes of Health and others.).

Exome-Based Rare-Variant Analyses in CKD.

BACKGROUND: Studies have identified many common genetic associations that influence renal function and all-cause CKD, but these explain only a small fraction of variance in these traits. The contribution of rare variants has not been systematically examined. METHODS: We performed exome sequencing of 3150 individuals, who collectively encompassed diverse CKD subtypes, and 9563 controls. To detect causal genes and evaluate the contribution of rare variants we used collapsing analysis, in which we compared the proportion of cases and controls carrying rare variants per gene. RESULTS: The analyses captured five established monogenic causes of CKD: variants in PKD1, PKD2, and COL4A5 achieved study-wide significance, and we observed suggestive case enrichment for COL4A4 and COL4A3. Beyond known disease-associated genes, collapsing analyses incorporating regional variant intolerance identified suggestive dominant signals in CPT2 and several other candidate genes. Biallelic mutations in CPT2 cause carnitine palmitoyltransferase II deficiency, sometimes associated with rhabdomyolysis and acute renal injury. Genetic modifier analysis among cases with APOL1 risk genotypes identified a suggestive signal in AHDC1, implicated in Xia-Gibbs syndrome, which involves intellectual disability and other features. On the basis of the observed distribution of rare variants, we estimate that a two- to three-fold larger cohort would provide 80% power to implicate new genes for all-cause CKD. CONCLUSIONS: This study demonstrates that rare-variant collapsing analyses can validate known genes and identify candidate genes and modifiers for kidney disease. In so doing, these findings provide a motivation for larger-scale investigation of rare-variant risk contributions across major clinical CKD categories.

A retrospective analysis of RET translocation, gene copy number gain and expression in NSCLC patients treated with vandetanib in four randomized Phase III studies.

BACKGROUND: To determine the prevalence of RET rearrangement genes, RET copy number gains and expression in tumor samples from four Phase III non-small-cell lung cancer (NSCLC) trials of vandetanib, a selective inhibitor of VEGFR, RET and EGFR signaling, and to determine any association with outcome to vandetanib treatment. METHODS: Archival tumor samples from the ZODIAC ( NCT00312377 , vandetanib ± docetaxel), ZEAL ( NCT00418886 , vandetanib ± pemetrexed), ZEPHYR ( NCT00404924 , vandetanib vs placebo) and ZEST ( NCT00364351 , vandetanib vs erlotinib) studies were evaluated by fluorescence in situ hybridization (FISH) and immunohistochemistry (IHC) in 944 and 1102 patients. RESULTS: The prevalence of RET rearrangements by FISH was 0.7% (95% CI 0.3-1.5%) among patients with a known result. Seven tumor samples were positive for RET rearrangements (vandetanib, n = 3; comparator, n = 4). 2.8% (n = 26) of samples had RET amplification (innumerable RET clusters, or ≥7 copies in > 10% of tumor cells), 8.1% (n = 76) had low RET gene copy number gain (4-6 copies in ≥40% of tumor cells) and 8.3% (n = 92) were RET expression positive (signal intensity ++ or +++ in >10% of tumor cells). Of RET-rearrangement-positive patients, none had an objective response in the vandetanib arm and one patient responded in the comparator arm. Radiologic evidence of tumor shrinkage was observed in two patients treated with vandetanib and one treated with comparator drug. The objective response rate was similar in the vandetanib and comparator arms for patients positive for RET copy number gains or RET protein expression. CONCLUSIONS: We have identified prevalence for three RET biomarkers in a population predominated by non-Asians and smokers. RET rearrangement prevalence was lower than previously reported. We found no evidence of a differential benefit for efficacy by IHC and RET gene copy number gains. The low prevalence of RET rearrangements (0.7%) prevents firm conclusions regarding association of vandetanib treatment with efficacy in the RET rearrangement NSCLC subpopulation. TRIAL REGISTRATION: Randomized Phase III clinical trials ( NCT00312377 , ZODIAC; NCT00418886 , ZEAL; NCT00364351 , ZEST; NCT00404924 , ZEPHYR).

Direct RNA sequencing.

Our understanding of human biology and disease is ultimately dependent on a complete understanding of the genome and its functions. The recent application of microarray and sequencing technologies to transcriptomics has changed the simplistic view of transcriptomes to a more complicated view of genome-wide transcription where a large fraction of transcripts emanates from unannotated parts of genomes, and underlined our limited knowledge of the dynamic state of transcription. Most of this broad body of knowledge was obtained indirectly because current transcriptome analysis methods typically require RNA to be converted to complementary DNA (cDNA) before measurements, even though the cDNA synthesis step introduces multiple biases and artefacts that interfere with both the proper characterization and quantification of transcripts. Furthermore, cDNA synthesis is not particularly suitable for the analysis of short, degraded and/or small quantity RNA samples. Here we report direct single molecule RNA sequencing without prior conversion of RNA to cDNA. We applied this technology to sequence femtomole quantities of poly(A)(+) Saccharomyces cerevisiae RNA using a surface coated with poly(dT) oligonucleotides to capture the RNAs at their natural poly(A) tails and initiate sequencing by synthesis. We observed transcript 3' end heterogeneity and polyadenylated small nucleolar RNAs. This study provides a path to high-throughput and low-cost direct RNA sequencing and achieving the ultimate goal of a comprehensive and bias-free understanding of transcriptomes.

Genetic architecture of telomere length in 462,666 UK Biobank whole-genome sequences.

Telomeres protect chromosome ends from damage and their length is linked with human disease and aging. We developed a joint telomere length metric, combining quantitative PCR and whole-genome sequencing measurements from 462,666 UK Biobank participants. This metric increased SNP heritability, suggesting that it better captures genetic regulation of telomere length. Exome-wide rare-variant and gene-level collapsing association studies identified 64 variants and 30 genes significantly associated with telomere length, including allelic series in ACD and RTEL1. Notably, 16% of these genes are known drivers of clonal hematopoiesis-an age-related somatic mosaicism associated with myeloid cancers and several nonmalignant diseases. Somatic variant analyses revealed gene-specific associations with telomere length, including lengthened telomeres in individuals with large SRSF2-mutant clones, compared with shortened telomeres in individuals with clonal expansions driven by other genes. Collectively, our findings demonstrate the impact of rare variants on telomere length, with larger effects observed among genes also associated with clonal hematopoiesis.

Genome-wide association analysis implicates the involvement of eight loci with response to tocilizumab for the treatment of rheumatoid arthritis.

Rheumatoid arthritis (RA) is an immune-mediated inflammatory disease affecting the joints. A heterogeneous response to available therapies demonstrates the need to identify those patients likely to benefit from a particular therapy. Our objective was to identify genetic factors associated with response to tocilizumab, a humanized monoclonal antibody targeting the interleukin (IL)-6 receptor, recently approved for treating RA. We report the first genome-wide association study on the response to tocilizumab in 1683 subjects with RA from six clinical studies. Putative associations were identified with eight loci, previously unrecognized as linked to the IL-6 pathway or associated with RA risk. This study suggests that it is unlikely that a major genetic determinant of response exists, and it illustrates the complexity of performing genome-wide association scans in clinical trials.

AZD9567 versus prednisolone in patients with active rheumatoid arthritis: A phase IIa, randomized, double-blind, efficacy, and safety study.

Oral corticosteroid use is limited by side effects, some caused by off-target actions on the mineralocorticoid receptor that disrupt electrolyte balance. AZD9567 is a selective, nonsteroidal glucocorticoid receptor modulator. The efficacy, safety, and tolerability of AZD9567 and prednisolone were assessed in a phase IIa study. Anti-inflammatory mechanism of action was also evaluated in vitro in monocytes from healthy donors. In this randomized, double-blind, parallel-group, multicenter study, patients with active rheumatoid arthritis were randomized 1:1 to AZD9567 40 mg or prednisolone 20 mg once daily orally for 14 days. The primary end point was change from baseline in DAS28-CRP at day 15. Secondary end points included components of DAS28-CRP, American College of Rheumatology (ACR) response criteria (ACR20, ACR50, and ACR70), and safety end points, including serum electrolytes. Overall, 21 patients were randomized to AZD9567 (n = 11) or prednisolone (n = 10), and all completed the study. As anticipated, AZD9567 had a similar efficacy profile to prednisolone, with no clinically meaningful (i.e., >1.0) difference in change from baseline to day 15 in DAS28-CRP between AZD9567 and prednisolone (least-squares mean difference: 0.47, 95% confidence interval: -0.49 to 1.43). Similar results were observed for the secondary efficacy end points. In vitro transcriptomic analysis showed that anti-inflammatory responses were similar for AZD9567, prednisolone, and dexamethasone. Unlike prednisolone, AZD9567 had no effect on the serum sodium:potassium ratio. The safety profile was not different from that of prednisolone. Larger studies of longer duration are required to determine whether AZD9567 40 mg may in the future be an alternative to prednisolone in patients with inflammatory disease.

Multiomics links global surfactant dysregulation with airflow obstruction and emphysema in COPD.

Rationale: Pulmonary surfactant is vital for lung homeostasis as it reduces surface tension to prevent alveolar collapse and provides essential immune-regulatory and antipathogenic functions. Previous studies demonstrated dysregulation of some individual surfactant components in COPD. We investigated relationships between COPD disease measures and dysregulation of surfactant components to gain new insights into potential disease mechanisms. Methods: Bronchoalveolar lavage proteome and lipidome were characterised in ex-smoking mild/moderate COPD subjects (n=26) and healthy ex-smoking (n=20) and never-smoking (n=16) controls using mass spectrometry. Serum surfactant protein analysis was performed. Results: Total phosphatidylcholine, phosphatidylglycerol, phosphatidylinositol, surfactant protein (SP)-B, SP-A and SP-D concentrations were lower in COPD versus controls (log2 fold change (log2FC) -2.0, -2.2, -1.5, -0.5, -0.7 and -0.5 (adjusted p<0.02), respectively) and correlated with lung function. Total phosphatidylcholine, phosphatidylglycerol, phosphatidylinositol, SP-A, SP-B, SP-D, napsin A and CD44 inversely correlated with computed tomography small airways disease measures (expiratory to inspiratory mean lung density) (r= -0.56, r= -0.58, r= -0.45, r= -0.36, r= -0.44, r= -0.37, r= -0.40 and r= -0.39 (adjusted p<0.05)). Total phosphatidylcholine, phosphatidylglycerol, phosphatidylinositol, SP-A, SP-B, SP-D and NAPSA inversely correlated with emphysema (% low-attenuation areas): r= -0.55, r= -0.61, r= -0.48, r= -0.51, r= -0.41, r= -0.31 and r= -0.34, respectively (adjusted p<0.05). Neutrophil elastase, known to degrade SP-A and SP-D, was elevated in COPD versus controls (log2FC 0.40, adjusted p=0.0390), and inversely correlated with SP-A and SP-D. Serum SP-D was increased in COPD versus healthy ex-smoking volunteers, and predicted COPD status (area under the curve 0.85). Conclusions: Using a multiomics approach, we demonstrate, for the first time, global surfactant dysregulation in COPD that was associated with emphysema, giving new insights into potential mechanisms underlying the cause or consequence of disease.

Virtual terminator nucleotides for next-generation DNA sequencing.

We synthesized reversible terminators with tethered inhibitors for next-generation sequencing. These were efficiently incorporated with high fidelity while preventing incorporation of additional nucleotides, and we used them to sequence canine bacterial artificial chromosomes in a single-molecule system that provided even coverage for over 99% of the region sequenced. This single-molecule approach generated high-quality sequence data without the need for target amplification and thus avoided concomitant biases.

The commutability of enzyme linked immunosorbent assays for the quantification of serum eosinophil-derived neurotoxin (EDN).

Eosinophil-derived neurotoxin (EDN) is a surrogate biomarker of eosinophil activation and has considerable potential as a precision medicine biomarker in diseases where eosinophils may play a causative role. Clinical data for EDN have been generated using different quantitative immunoassays, but comparisons between these individual data sets are challenging as no internationally recognised EDN standards or orthogonal methods exist. In this study we aimed to compare commercial EDN assays from ALPCO, MBL, LSBio and CUSABIO for sample commutability. Firstly, we analytically validated the ALPCO enzyme linked immunosorbent assay (ELISA) and demonstrated appropriate analytical characteristics, including an intra/inter-assay precision coefficient-of-variation of between 1.9 and 6.8%. EDN purified from blood proved to be a good quality control material, whereas recombinant EDN, expressed in E.coli, did not react in the ALPCO immunoassay. Using healthy and asthma patient serum samples we confirmed that the ALPCO assay correlated well with the MBL assay, with a coefficient of determination (R2) of 0.92. However, the results from LSBio and CUSABIO assays were not commutable to the other assays.

Genetic variation in UGT1A1 typical of Gilbert syndrome is associated with unconjugated hyperbilirubinemia in patients receiving tocilizumab.

OBJECTIVE: Tocilizumab, a monoclonal antibody to interleukin-6 receptor, was recently approved for the treatment of moderate-to-severe rheumatoid arthritis. Two patients during clinical development met laboratory, but not clinical, criteria for Hy's law with bilirubin elevations suspected as a result of genetic variation in uridine diphosphoglucose glucuronosyltransferase (UGT1A1) typical of Gilbert syndrome. METHODS: Genotyping of the two cases potentially meeting with Hy's law was performed using commercially available procedures. UGT1A1 single nucleotide polymorphism data were extracted from a genome-wide array database for 1187 patients from tocilizumab trials, and associations of UGT1A1 genotypes with bilirubin elevations were analyzed using logistic regression for associations with baseline and change from baseline in bilirubin levels as continuous variables. RESULTS: Bilirubin elevations were not associated with clinical adverse events. Both patients potentially meeting Hy's law carry homozygous UGT1A1*28 alleles and UGT1A1*60 alleles. UGT1A1*28 and three additional single nucleotide polymorphisms showed odds ratios greater than 25 for associations with elevated bilirubin. The presence of rs6742078 accounted for 32% of the total variance in bilirubin (P=2.2×10). CONCLUSION: Bilirubin increases occurring with tocilizumab appear to be related to anti-inflammatory effects extending to the liver. Thus, in the absence of other signs of hepatic dysfunction, bilirubin elevations after treatment with tocilizumab have a high probability of association with UGT1A1 polymorphism, which should alleviate concerns of serious hepatotoxicity. Our results underscore the value of genotyping in the clinical trial setting to avoid misinterpretations that could lead to terminating development of a promising new agent.

c-Jun N-terminal kinase primes endothelial cells at atheroprone sites for apoptosis.

OBJECTIVE: Atherosclerosis is a focal disease that occurs predominantly at branches and bends of the arterial tree. Endothelial cells (EC) at atherosusceptible sites are prone to injury, which can contribute to lesion formation, whereas EC at atheroprotected sites are resistant. The c-Jun N-terminal kinase (JNK) is activated constitutively in EC at atherosusceptible sites but is inactivated at atheroprotected sites by mitogen-activated protein kinase phosphatase-1 (MKP-1). Here, we examined the effects of JNK activation on EC physiology at atherosusceptible sites. METHODS AND RESULTS: We identified transcriptional programs regulated by JNK by applying a specific pharmacological inhibitor to cultured EC and assessing the transcriptome using microarrays. This approach and subsequent validation by gene silencing revealed that JNK positively regulates the expression of numerous proapoptotic molecules. Analysis of aortae of wild-type, JNK1(-/-), and MKP-1(-/-) mice revealed that EC at an atherosusceptible site express proapoptotic proteins and are primed for apoptosis and proliferation in response to lipopolysaccharide through a JNK1-dependent mechanism, whereas EC at a protected site expressed lower levels of proapoptotic molecules and were protected from injury by MKP-1. CONCLUSIONS: Spatial variation of JNK1 activity delineates the spatial distribution of apoptosis and turnover of EC in arteries.