Treatment response and neurofilament light chain levels with long-term patisiran in hereditary transthyretin-mediated amyloidosis with polyneuropathy: 24-month results of an open-label extension study

BACKGROUND

Longitudinal changes in neurofilament light chain (NfL) levels were evaluated alongside prespecified clinical assessments 24 months into the patisiran Global open-label extension (OLE) study in patients with ATTRv amyloidosis with polyneuropathy.

METHODS

All patients enrolled in the Global OLE, from phase III APOLLO and phase II OLE parent studies, received patisiran. Assessments included measures of polyneuropathy (modified Neuropathy Impairment Score+7 (mNIS+7)), quality of life (QOL; Norfolk QOL-Diabetic Neuropathy questionnaire (Norfolk QOL-DN)), and plasma NfL.

RESULTS

Patients receiving patisiran in the parent study (APOLLO-patisiran, n = 137; phase II OLE-patisiran, n = 25) demonstrated sustained improvements in mNIS+7 (mean change from parent study baseline (95% confidence interval): APOLLO-patisiran -4.8 (-8.9, -0.6); phase II OLE-patisiran -5.8 (-10.5, -1.2)) and Norfolk QOL-DN (APOLLO-patisiran -2.4 (-7.2, 2.3)), and maintained reduced NfL levels at Global OLE 24 months. After initiating patisiran in the Global OLE, APOLLO-placebo patients (n = 49) demonstrated stabilized mNIS+7, improved Norfolk QOL-DN, and significantly reduced NfL levels. Patisiran continued to demonstrate an acceptable safety profile. Earlier patisiran initiation was associated with a lower exposure-adjusted mortality rate.

CONCLUSIONS

Long-term patisiran treatment led to sustained improvements in neuropathy and QOL, with NfL demonstrating potential as a biomarker for disease progression and treatment response in ATTRv amyloidosis with polyneuropathy.

Transcript-aware analysis of rare predicted loss-of-function variants in the UK Biobank elucidate new isoform-trait associations

A single gene can produce multiple transcripts with distinct molecular functions. Rare-variant association tests often aggregate all coding variants across individual genes, without accounting for the variants' presence or consequence in resulting transcript isoforms. To evaluate the utility of transcript-aware variant sets, rare predicted loss-of-function (pLOF) variants were aggregated for 17,035 protein-coding genes using 55,558 distinct transcript-specific variant sets. These sets were tested for their association with 728 circulating proteins and 188 quantitative phenotypes across 406,921 individuals in the UK Biobank. The transcript-specific approach resulted in larger estimated effects of pLOF variants decreasing serum cis-protein levels compared to the gene-based approach (pbinom ≤ 2x10-16). Additionally, 251 quantitative trait associations were identified as being significant using the transcript-specific approach but not the gene-based approach, including PCSK5 transcript ENST00000376752 and standing height (transcript-specific statistic, P = 1.3x10-16, effect = 0.7 SD decrease; gene-based statistic, P = 0.02, effect = 0.05 SD decrease) and LDLR transcript ENST00000252444 and apolipoprotein B (transcript-specific statistic, P = 5.7x10-20, effect = 1.0 SD increase; gene-based statistic, P = 3.0x10-4, effect = 0.2 SD increase). This approach demonstrates the importance of considering the effect of pLOFs on specific transcript isoforms when performing rare-variant association studies.

Plasma proteomic associations with genetics and health in the UK Biobank

The Pharma Proteomics Project is a precompetitive biopharmaceutical consortium characterizing the plasma proteomic profiles of 54,219 UK Biobank participants. Here we provide a detailed summary of this initiative, including technical and biological validations, insights into proteomic disease signatures, and prediction modelling for various demographic and health indicators. We present comprehensive protein quantitative trait locus (pQTL) mapping of 2,923 proteins that identifies 14,287 primary genetic associations, of which 81% are previously undescribed, alongside ancestry-specific pQTL mapping in non-European individuals. The study provides an updated characterization of the genetic architecture of the plasma proteome, contextualized with projected pQTL discovery rates as sample sizes and proteomic assay coverages increase over time. We offer extensive insights into trans pQTLs across multiple biological domains, highlight genetic influences on ligand-receptor interactions and pathway perturbations across a diverse collection of cytokines and complement networks, and illustrate long-range epistatic effects of ABO blood group and FUT2 secretor status on proteins with gastrointestinal tissue-enriched expression. We demonstrate the utility of these data for drug discovery by extending the genetic proxied effects of protein targets, such as PCSK9, on additional endpoints, and disentangle specific genes and proteins perturbed at loci associated with COVID-19 susceptibility. This public-private partnership provides the scientific community with an open-access proteomics resource of considerable breadth and depth to help to elucidate the biological mechanisms underlying proteo-genomic discoveries and accelerate the development of biomarkers, predictive models and therapeutics1.

Blood pressure-independent renoprotective effects of small interference RNA targeting liver angiotensinogen in experimental diabetes

BACKGROUND AND PURPOSE

Small interfering RNA (siRNA) targeting liver angiotensinogen lowers blood pressure, but its effects in hypertensive diabetes are unknown.

EXPERIMENTAL APPROACH

To address this, TGR (mRen2)27 rats (angiotensin II-dependent hypertension model) were made diabetic with streptozotocin over 18 weeks and treated with either vehicle, angiotensinogen siRNA, the AT₁ antagonist valsartan, the ACE inhibitor captopril, valsartan + siRNA or valsartan + captopril for the final 3 weeks. Mean arterial pressure (MAP) was measured via radiotelemetry.

KEY RESULTS

MAP before treatment was 153 ± 2 mmHg. Diabetes resulted in albuminuria, accompanied by glomerulosclerosis and podocyte effacement, without a change in glomerular filtration rate. All treatments lowered MAP and cardiac hypertrophy, and the largest drop in MAP was observed with siRNA + valsartan. Treatment with siRNA lowered circulating angiotensinogen by >99%, and the lowest circulating angiotensin II and aldosterone levels occurred in the dual treatment groups. Angiotensinogen siRNA did not affect renal angiotensinogen mRNA expression, confirming its liver-specificity. Furthermore, only siRNA with or without valsartan lowered renal angiotensin I. All treatments lowered renal angiotensin II and the reduction was largest (>95%) in the siRNA + valsartan group. All treatments identically lowered albuminuria, whereas only siRNA with or without valsartan restored podocyte foot processes and reduced glomerulosclerosis.

CONCLUSION AND IMPLICATIONS

Angiotensinogen siRNA exerts renoprotection in diabetic TGR (mRen2)27 rats and this relies, at least in part, on the suppression of renal angiotensin II formation from liver-derived angiotensinogen. Clinical trials should now address whether this is also beneficial in human diabetic kidney disease.

Conventional Vasopressor and Vasopressor-Sparing Strategies to Counteract the Blood Pressure-Lowering Effect of Small Interfering RNA Targeting Angiotensinogen

Background

A single dose of small interfering RNA (siRNA) targeting liver angiotensinogen eliminates hepatic angiotensinogen and lowers blood pressure. Angiotensinogen elimination raises concerns for clinical application because an angiotensin rise is needed to maintain perfusion pressure during hypovolemia. Here, we investigated whether conventional vasopressors can raise arterial pressure after angiotensinogen depletion.

Methods and Results

Spontaneously hypertensive rats on a low‐salt diet were treated with siRNA (10 mg/kg fortnightly) for 4 weeks, supplemented during the final 2 weeks with fludrocortisone (6 mg/kg per day), the α‐adrenergic agonist midodrine (4 mg/kg per day), or a high‐salt diet (all groups n=6–7). Pressor responsiveness to angiotensin II and norepinephrine was assessed before and after siRNA administration. Blood pressure was measured via radiotelemetry. Depletion of liver angiotensinogen by siRNA lowered plasma angiotensinogen concentrations by 99.2±0.1% and mean arterial pressure by 19 mm Hg. siRNA‐mediated blood pressure lowering was rapidly reversed by intravenous angiotensin II or norepinephrine, or gradually reversed by fludrocortisone or high salt intake. Midodrine had no effect. Unexpectedly, fludrocortisone partially restored plasma angiotensinogen concentrations in siRNA‐treated rats, and nearly abolished plasma renin concentrations. To investigate whether this angiotensinogen originated from nonhepatic sources, fludrocortisone was administered to mice lacking hepatic angiotensinogen. Fludrocortisone did not increase angiotensinogen in these mice, implying that the rise in angiotensinogen in the siRNA‐treated rats must have depended on the liver, most likely reflecting diminished cleavage by renin.

Conclusions

Intact pressor responsiveness to conventional vasopressors provides pharmacological means to regulate the blood pressure–lowering effect of angiotensinogen siRNA and may support future therapeutic implementation of siRNA.

Rare loss of function variants in the hepatokine gene INHBE protect from abdominal obesity

Identifying genetic variants associated with lower waist-to-hip ratio can reveal new therapeutic targets for abdominal obesity. We use exome sequences from 362,679 individuals to identify genes associated with waist-to-hip ratio adjusted for BMI (WHRadjBMI), a surrogate for abdominal fat that is causally linked to type 2 diabetes and coronary heart disease. Predicted loss of function (pLOF) variants in INHBE associate with lower WHRadjBMI and this association replicates in data from AMP-T2D-GENES. INHBE encodes a secreted protein, the hepatokine activin E. In vitro characterization of the most common INHBE pLOF variant in our study, indicates an in-frame deletion resulting in a 90% reduction in secreted protein levels. We detect associations with lower WHRadjBMI for variants in ACVR1C, encoding an activin receptor, further highlighting the involvement of activins in regulating fat distribution. These findings highlight activin E as a potential therapeutic target for abdominal obesity, a phenotype linked to cardiometabolic disease.

Rare coding variants in DNA damage repair genes associated with timing of natural menopause

The age of menopause is associated with fertility and disease risk, and its genetic control is of great interest. We use whole-exome sequences from 132,370 women in the UK Biobank to test for associations between rare damaging variants and age at natural menopause. Rare damaging variants in five genes are significantly associated with menopause: CHEK2 (p = 3.3 × 10⁻⁵¹), DCLRE1A (p = 8.4 × 10⁻¹³), and HELB (p = 5.7 × 10⁻⁷) with later menopause and TOP3A (p = 7.6 × 10⁻⁸) and CLPB (p = 8.1 × 10⁻⁷) with earlier menopause. Two additional genes are suggestive: RAD54L (p = 2.4 × 10⁻⁶) with later menopause and HROB (p = 2.9 × 10⁻⁶) with earlier menopause. In a follow-up analysis of repeated questionnaires in women who were initially premenopausal, CHEK2, TOP3A, and RAD54L genotypes are associated with subsequent menopause. Consistent with previous genome-wide association studies (GWASs), six of the seven genes are involved in the DNA damage repair pathway. Phenome-wide scans across 398,569 men and women revealed that in addition to known associations with cancers and blood cell counts, rare variants in CHEK2 are also associated with increased risk for uterine fibroids, polycystic ovary syndrome, and prostate hypertrophy; these associations are not shared with higher-penetrance breast cancer genes. Causal mediation analysis suggests that approximately 8% of the breast cancer risk conferred by CHEK2 pathogenic variants after menopause is mediated through delayed menopause.

Gene-level analysis of rare variants in 379,066 whole exome sequences identifies an association of GIGYF1 loss of function with type 2 diabetes

Sequencing of large cohorts offers an unprecedented opportunity to identify rare genetic variants and to find novel contributors to human disease. We used gene-based collapsing tests to identify genes associated with glucose, HbA1c and type 2 diabetes (T2D) diagnosis in 379,066 exome-sequenced participants in the UK Biobank. We identified associations for variants in GCK, HNF1A and PDX1, which are known to be involved in Mendelian forms of diabetes. Notably, we uncovered novel associations for GIGYF1, a gene not previously implicated by human genetics in diabetes. GIGYF1 predicted loss of function (pLOF) variants associated with increased levels of glucose (0.77 mmol/L increase, p = 4.42 × 10^–12) and HbA1c (4.33 mmol/mol, p = 1.28 × 10^–14) as well as T2D diagnosis (OR = 4.15, p = 6.14 × 10^–11). Multiple rare variants contributed to these associations, including singleton variants. GIGYF1 pLOF also associated with decreased cholesterol levels as well as an increased risk of hypothyroidism. The association of GIGYF1 pLOF with T2D diagnosis replicated in an independent cohort from the Geisinger Health System. In addition, a common variant association for glucose and T2D was identified at the GIGYF1 locus. Our results highlight the role of GIGYF1 in regulating insulin signaling and protecting from diabetes.

GWAS of serum ALT and AST reveals an association of SLC30A10 Thr95Ile with hypermanganesemia symptoms

Understanding mechanisms of hepatocellular damage may lead to new treatments for liver disease, and genome-wide association studies (GWAS) of alanine aminotransferase (ALT) and aspartate aminotransferase (AST) serum activities have proven useful for investigating liver biology. Here we report 100 loci associating with both enzymes, using GWAS across 411,048 subjects in the UK Biobank. The rare missense variant SLC30A10 Thr95Ile (rs188273166) associates with the largest elevation of both enzymes, and this association replicates in the DiscovEHR study. SLC30A10 excretes manganese from the liver to the bile duct, and rare homozygous loss of function causes the syndrome hypermanganesemia with dystonia-1 (HMNDYT1) which involves cirrhosis. Consistent with hematological symptoms of hypermanganesemia, SLC30A10 Thr95Ile carriers have increased hematocrit and risk of iron deficiency anemia. Carriers also have increased risk of extrahepatic bile duct cancer. These results suggest that genetic variation in SLC30A10 adversely affects more individuals than patients with diagnosed HMNDYT1.

Pan-ancestry exome-wide association analyses of COVID-19 outcomes in 586,157 individuals

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) causes coronavirus disease 2019 (COVID-19), a respiratory illness that can result in hospitalization or death. We used exome sequence data to investigate associations between rare genetic variants and seven COVID-19 outcomes in 586,157 individuals, including 20,952 with COVID-19. After accounting for multiple testing, we did not identify any clear associations with rare variants either exome wide or when specifically focusing on (1) 13 interferon pathway genes in which rare deleterious variants have been reported in individuals with severe COVID-19, (2) 281 genes located in susceptibility loci identified by the COVID-19 Host Genetics Initiative, or (3) 32 additional genes of immunologic relevance and/or therapeutic potential. Our analyses indicate there are no significant associations with rare protein-coding variants with detectable effect sizes at our current sample sizes. Analyses will be updated as additional data become available, and results are publicly available through the Regeneron Genetics Center COVID-19 Results Browser.

Association of the transthyretin variant V122I with polyneuropathy among individuals of African ancestry

Hereditary transthyretin-mediated (hATTR) amyloidosis is an underdiagnosed, progressively debilitating disease caused by mutations in the transthyretin (TTR) gene. V122I, a common pathogenic TTR mutation, is found in 3-4% of individuals of African ancestry in the United States and has been associated with cardiomyopathy and heart failure. To better understand the phenotypic consequences of carrying V122I, we conducted a phenome-wide association study scanning 427 ICD diagnosis codes in UK Biobank participants of African ancestry (n = 6062). Significant associations were tested for replication in the Penn Medicine Biobank (n = 5737) and the Million Veteran Program (n = 82,382). V122I was significantly associated with polyneuropathy in the UK Biobank (odds ratio [OR] = 6.4, 95% confidence interval [CI] 2.6-15.6, p = 4.2 × 10-5), which was replicated in the Penn Medicine Biobank (OR = 1.6, 95% CI 1.2-2.4, p = 6.0 × 10-3) and Million Veteran Program (OR = 1.5, 95% CI 1.2-1.8, p = 1.8 × 10-4). Polyneuropathy prevalence among V122I carriers was 2.1%, 9.0%, and 4.8% in the UK Biobank, Penn Medicine Biobank, and Million Veteran Program, respectively. The cumulative incidence of common hATTR amyloidosis manifestations (carpal tunnel syndrome, polyneuropathy, cardiomyopathy, heart failure) was significantly enriched in V122I carriers compared with non-carriers (HR = 2.8, 95% CI 1.7-4.5, p = 2.6 × 10-5) in the UK Biobank, with 37.4% of V122I carriers having at least one of these manifestations by age 75. Our findings show that V122I carriers are at increased risk of polyneuropathy. These results also emphasize the underdiagnosis of disease in V122I carriers with a significant proportion of subjects showing phenotypic changes consistent with hATTR amyloidosis. Greater understanding of the manifestations associated with V122I is critical for earlier diagnosis and treatment.

A catalog of associations between rare coding variants and COVID-19 outcomes

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) causes coronavirus disease-19 (COVID-19), a respiratory illness that can result in hospitalization or death. We investigated associations between rare genetic variants and seven COVID-19 outcomes in 543,213 individuals, including 8,248 with COVID-19. After accounting for multiple testing, we did not identify any clear associations with rare variants either exome-wide or when specifically focusing on (i) 14 interferon pathway genes in which rare deleterious variants have been reported in severe COVID-19 patients; (ii) 167 genes located in COVID-19 GWAS risk loci; or (iii) 32 additional genes of immunologic relevance and/or therapeutic potential. Our analyses indicate there are no significant associations with rare protein-coding variants with detectable effect sizes at our current sample sizes. Analyses will be updated as additional data become available, with results publicly browsable at https://rgc-covid19.regeneron.com.

The V122I mutation in hereditary transthyretin-mediated amyloidosis is significantly associated with an increased incidence of polyneuropathy

Background

Hereditary transthyretin-mediated (hATTR) amyloidosis is a progressive, life-threatening disease caused by mutations in the transthyretin (TTR) gene. hATTR amyloidosis phenotypes can vary by patient and mutation. The V122I (Val122Ile; p.V142I) variant is one of the most common pathogenic TTR mutations, is primarily found in people of West African descent and has historically been associated with cardiomyopathy (CM).

Purpose

To characterize the cumulative incidence of diagnoses frequently seen with hATTR amyloidosis in V122I carriers and non-carriers in the UK Biobank (UKBB) and the Penn Medicine BioBank (PMBB).

Methods

UKBB and PMBB are prospective studies with ∼500,000 and ∼60,000 subjects, respectively. Clinical presentations frequently seen with hATTR amyloidosis were assessed using ICD10 diagnosis codes: G62–polyneuropathy (PN), I50 or I098–heart failure (HF), G560–carpal tunnel syndrome (CTS), I42–CM, and E85–amyloidosis. The cumulative incidence of diagnoses was estimated using Kaplan-Meier curves. Time to first hATTR amyloidosis-related diagnosis was compared between V122I carriers and non-carriers using Cox proportional hazards regression, controlling for age, sex, smoking, and genetic ancestry.

Results

Of the 6,062 unrelated black participants in the UKBB, 243 were V122I carriers. Only 0.8% of V122I carriers had a formal diagnosis of hATTR amyloidosis. V122I carriers were significantly more likely to have a PN diagnosis than non-carriers (p=6.35x10–5), a finding which was replicated in the PMBB. Of the ICD10 codes assessed, Cox proportional hazards regression revealed a significant association between V122I genotype and time to first diagnosis (p=2.6x10–5), with 11.1% of V122I carriers having at least one diagnosis during follow-up versus 4.9% of non-carriers. The calculated population attributable risk showed an excess risk of 16.7% for a PN diagnosis, 6.5% for HF, 4.1% for CTS, and 2.4% for CM in V122I carriers. The cumulative incidence of any hATTR amyloidosis-related diagnosis among V122I carriers by age 65 was 11.9% (95% CI=3.1–19.8%). The incidence increased to 37.4% (95% CI=20.5–50.7%) by age 75, which was significantly higher than non-carriers (13.8%, 95% CI=11.6–16%). Additionally, an assessment of the cumulative incidence of each diagnosis separately revealed that PN became more prevalent at younger ages, while HF and CM became more prevalent at older ages.

Conclusions

V122I carriers were significantly more likely to receive diagnoses frequently seen with hATTR amyloidosis than non-carriers. Although these diagnoses may not all be attributed to amyloidosis, the small fraction of patients diagnosed with hATTR amyloidosis suggests a potential underdiagnosis of disease. The V122I mutation has historically been associated with a cardiac phenotype, yet these data also suggest an increased incidence of PN. Increased vigilance for the mixed phenotype may lead to earlier diagnoses and treatment of V122I carriers.

Funding Acknowledgement

Type of funding source: Private company. Main funding source(s): Alnylam Pharmaceuticals

Neurofilament Light Chain as a Biomarker of Hereditary Transthyretin-Mediated Amyloidosis

OBJECTIVE

To identify changes in the proteome associated with onset and progression of hereditary transthyretin-mediated (hATTR) amyloidosis, also known as ATTRv amyloidosis, we performed an observational, case-controlled study that compared proteomes of patients with ATTRv amyloidosis and healthy controls.

METHODS

Plasma levels of >1,000 proteins were measured in patients with ATTRv amyloidosis with polyneuropathy who received either placebo or patisiran in a Phase 3 study of patisiran (APOLLO), and in healthy controls. The effect of patisiran on the time profile of each protein was determined by linear mixed model at 0, 9, and 18 months. Neurofilament light chain (NfL) was further assessed with an orthogonal quantitative approach.

RESULTS

Levels of 66 proteins were significantly changed with patisiran vs placebo, with NfL change most significant (p < 10^-20). Analysis of changes in protein levels demonstrated that the proteome of patients treated with patisiran trended toward that of healthy controls at 18 months. Healthy controls' NfL levels were 4-fold lower than in patients with ATTRv amyloidosis with polyneuropathy (16.3 pg/mL vs 69.4 pg/mL, effect -53.1 pg/mL [95% confidence interval -60.5 to -45.9]). NfL levels at 18 months increased with placebo (99.5 pg/mL vs 63.2 pg/mL, effect 36.3 pg/mL [16.5-56.1]) and decreased with patisiran treatment (48.8 pg/mL vs 72.1 pg/mL, effect -23.3 pg/mL [-33.4 to -13.1]) from baseline. At 18 months, improvement in modified Neuropathy Impairment Score +7 score after patisiran treatment significantly correlated with reduced NfL (R = 0.43 [0.29-0.55]).

CONCLUSIONS

Findings suggest that NfL may serve as a biomarker of nerve damage and polyneuropathy in ATTRv amyloidosis, enable earlier diagnosis of patients with ATTRv amyloidosis, and facilitate monitoring of disease progression.

CLASSIFICATION OF EVIDENCE

This study provides Class III evidence that NfL levels may enable earlier diagnosis of polyneuropathy in patients with ATTRv amyloidosis and facilitate monitoring of disease progression.

Transthyretin-stabilising mutation T119M is not associated with protection against vascular disease or death in the UK Biobank

Background: Destabilised transthyretin (TTR) can result in the progressive, fatal disease transthyretin-mediated (ATTR) amyloidosis. A stabilising TTR mutation, T119M, is the basis for a therapeutic strategy to reduce destabilised TTR. Recently, T119M was associated with extended lifespan and lower risk of cerebrovascular disease in a Danish cohort. We aimed to determine whether this finding could be replicated in the UK Biobank.Methods: TTR T119M carriers were identified in the UK Biobank, a large prospective cohort of ∼500,000 individuals. Association between T119M genotype and inpatient diagnosis of vascular disease, cardiovascular disease, cerebrovascular disease, and mortality was analysed.Results: Frequency of T119M within the white UK Biobank population (n = 337,148) was 0.4%. Logistic regression comparing T119M carriers to non-carriers found no association between T119M and vascular disease (odds ratio [OR] = 1.08; p = .27), cardiovascular disease (OR = 1.08; p = .31), cerebrovascular disease (OR = 1.1; p = .42), or death (OR = 1.2; p = .06). Cox proportional hazards regression showed similar results (hazard ratio >1, p>.05). Age at death and vascular disease diagnosis were similar between T119M carriers and non-carriers (p = .12 and p = .38, respectively).Conclusions: There was no association between the TTR T119M genotype and risk of vascular disease or death in a large prospective cohort study, indicating that TTR tetramer stabilisation through T119M is not protective in this setting.

Characterising a healthy adult with a rare HAO1 knockout to support a therapeutic strategy for primary hyperoxaluria

By sequencing autozygous human populations, we identified a healthy adult woman with lifelong complete knockout of HAO1 (expected ~1 in 30 million outbred people). HAO1 (glycolate oxidase) silencing is the mechanism of lumasiran, an investigational RNA interference therapeutic for primary hyperoxaluria type 1. Her plasma glycolate levels were 12 times, and urinary glycolate 6 times, the upper limit of normal observed in healthy reference individuals (n = 67). Plasma metabolomics and lipidomics (1871 biochemicals) revealed 18 markedly elevated biochemicals (>5 sd outliers versus n = 25 controls) suggesting additional HAO1 effects. Comparison with lumasiran preclinical and clinical trial data suggested she has <2% residual glycolate oxidase activity. Cell line p.Leu333SerfsTer4 expression showed markedly reduced HAO1 protein levels and cellular protein mis-localisation. In this woman, lifelong HAO1 knockout is safe and without clinical phenotype, de-risking a therapeutic approach and informing therapeutic mechanisms. Unlocking evidence from the diversity of human genetic variation can facilitate drug development.

Rationalizing Secondary Pharmacology Screening Using Human Genetic and Pharmacological Evidence

Safety-related drug failures remain a major challenge for the pharmaceutical industry. One approach to ensuring drug safety involves assessing small molecule drug specificity by examining the ability of a drug candidate to interact with a panel of “off-target” proteins, referred to as secondary pharmacology screening. Information from human genetics and pharmacology can be used to select proteins associated with adverse effects for such screening. In an analysis of marketed drugs, we found a clear relationship between the genetic and pharmacological phenotypes of a drug’s off-target proteins and the observed drug side effects. In addition to using this phenotypic information for the selection of secondary pharmacology screens, we also show that it can be used to help identify drug off-target protein interactions responsible for drug-related adverse events. We anticipate that this phenotype-driven approach to secondary pharmacology screening will help to reduce safety-related drug failures due to drug off-target protein interactions.

Combining hiPSCs and Human Genetics: Major Applications in Drug Development

Merging iPSC models and human genetic research has opened up new avenues in understanding disease mechanisms and target biology, which facilitate exciting translation of this research to many areas of drug development. We highlight recent applications of these combined disciplines and discuss remaining challenges and potential solutions.

Variant ASGR1 Associated with a Reduced Risk of Coronary Artery Disease

BACKGROUND

Several sequence variants are known to have effects on serum levels of non-high-density lipoprotein (HDL) cholesterol that alter the risk of coronary artery disease.

METHODS

We sequenced the genomes of 2636 Icelanders and found variants that we then imputed into the genomes of approximately 398,000 Icelanders. We tested for association between these imputed variants and non-HDL cholesterol levels in 119,146 samples. We then performed replication testing in two populations of European descent. We assessed the effects of an implicated loss-of-function variant on the risk of coronary artery disease in 42,524 case patients and 249,414 controls from five European ancestry populations. An augmented set of genomes was screened for additional loss-of-function variants in a target gene. We evaluated the effect of an implicated variant on protein stability.

RESULTS

We found a rare noncoding 12-base-pair (bp) deletion (del12) in intron 4 of ASGR1, which encodes a subunit of the asialoglycoprotein receptor, a lectin that plays a role in the homeostasis of circulating glycoproteins. The del12 mutation activates a cryptic splice site, leading to a frameshift mutation and a premature stop codon that renders a truncated protein prone to degradation. Heterozygous carriers of the mutation (1 in 120 persons in our study population) had a lower level of non-HDL cholesterol than noncarriers, a difference of 15.3 mg per deciliter (0.40 mmol per liter) (P=1.0×10(-16)), and a lower risk of coronary artery disease (by 34%; 95% confidence interval, 21 to 45; P=4.0×10(-6)). In a larger set of sequenced samples from Icelanders, we found another loss-of-function ASGR1 variant (p.W158X, carried by 1 in 1850 persons) that was also associated with lower levels of non-HDL cholesterol (P=1.8×10(-3)).

CONCLUSIONS

ASGR1 haploinsufficiency was associated with reduced levels of non-HDL cholesterol and a reduced risk of coronary artery disease. (Funded by the National Institutes of Health and others.).

Time-dependent cellular and transcriptional changes in the osteoblast lineage associated with sclerostin antibody treatment in ovariectomized rats

Inhibition of sclerostin with sclerostin antibody (Scl-Ab) has been shown to stimulate bone formation, decrease bone resorption, and increase bone mass in both animals and humans. To obtain insight into the temporal cellular and transcriptional changes in the osteoblast (OB) lineage associated with long-term Scl-Ab treatment, stereological and transcriptional analyses of the OB lineage were performed on lumbar vertebrae from aged ovariectomized rats. Animals were administered Scl-Ab 3 or 50mg/kg/wk or vehicle (VEH) for up to 26weeks (d183), followed by a treatment-free period (TFP). At 50mg/kg/wk, bone volume (BV/total volume [TV]) increased through d183 and declined during the TFP. Bone formation rate (BFR/bone surface [BS]) and total OB number increased through d29, then progressively declined, coincident with a decrease in total osteoprogenitor (OP) numbers from d29 through d183. Analysis of differentially expressed genes (DEGs) from microarray analysis of mRNA isolated from laser capture microdissection samples enriched for OB, lining cells, and osteocytes (OCy) revealed modules of genes that correlated with BFR/BS, BV/TV, and osteoblastic surface (Ob.S)/BS. Expression change of canonical Wnt target genes was similar in all three cell types at d8, including upregulation of Twist1 and Wisp1. At d29, the pattern of Wnt target gene expression changed in the OCy, with Twist1 returning to VEH level, sustained upregulation of Wisp1, and upregulation of several other Wnt targets that continued into the TFP. Predicted activation of pathways recognized to integrate with and regulate canonical Wnt signaling were also activated at d29 in the OCy. The most significantly affected pathways represented transcription factor signaling known to inhibit cell cycle progression (notably p53) and mitogenesis (notably c-Myc). These changes occurred at the time of peak BFR/BS and continued as BFR/BS declined during treatment, then trended toward VEH level in the TFP. Concurrent with this transcriptional switch was a reduction in OP numbers, an effect that would ultimately limit bone formation. This study confirms that the initial transcriptional response in response to Scl-Ab is activation of canonical Wnt signaling and the data demonstrate that there is induction of additional regulatory pathways in OCy with long-term treatment. The interactions between Wnt and p53/c-Myc signaling may be key in limiting OP populations, thus contributing to self-regulation of bone formation with continued Scl-Ab administration.

Transcriptional Profiling of Laser Capture Microdissected Subpopulations of the Osteoblast Lineage Provides Insight Into the Early Response to Sclerostin Antibody in Rats

Sclerostin antibody (Scl-Ab) increases bone formation through a process dependent on the activation of canonical Wnt signaling, although the specific signaling in the osteoblast lineage in vivo is largely unknown. To gain insight into the signaling pathways acutely modulated by Scl-Ab, the transcriptional response of subpopulations of the osteoblast lineage was assessed by TaqMan and microarray analyses of mRNA isolated from laser capture microdissection (LCM)-enriched samples from the vertebrae of ovariectomized rats during the first week after Scl-Ab administration. Briefly, 6-month-old Sprague-Dawley rats were ovariectomized and, after 2 months, received a single dose of vehicle (VEH) or 100 mg/kg Scl-Ab (n = 20/group). Lumbar vertebrae were collected at 6, 24, 72, and 168 hours postdose and cryosectioned for LCM. Osteocytes were captured from bone matrix, and osteoblasts and lining cells were captured from bone surfaces based on fluorochrome labeling. mRNA was isolated, amplified, and profiled by TaqMan and microarray. Expression analysis revealed that Scl-Ab caused strikingly similar transcriptional profiles across all three cell types. Only 13 known canonical Wnt target genes, the majority with known functions in bone, showed a significant change in expression by microarray in response to Scl-Ab, with Wisp1 and Twist1 being the most responsive. Coincident with increased expression of Wnt target genes was the upregulation of numerous extracellular matrix (ECM) genes. The acute and progressive upregulation of ECM genes in lining cells supports their activation into matrix-producing osteoblasts, consistent with modeling-based bone formation. A similar transcriptional profile in osteocytes may indicate that Scl-Ab stimulates perilacunar/pericanalicular matrix deposition. Pathway analyses indicated that Scl-Ab regulated a limited number of genes related to cell cycle arrest and B-cell development. These data describe the acute downstream signaling in response to Scl-Ab in vivo and demonstrate selected canonical Wnt target gene activation associated with increased bone formation in all mature osteoblast subpopulations.

Utilization of human nuclear receptors as an early counter screen for off-target activity: a case study with a compendium of 615 known drugs

Off-target effects of drugs on nuclear hormone receptors (NHRs) may result in adverse effects in multiple organs/physiological processes. Reliable assessments of the NHR activities for drug candidates are therefore crucial for drug development. However, the highly permissive structures of NHRs for vastly different ligands make it challenging to predict interactions by examining the chemical structures of the ligands. Here, we report a detailed investigation on the agonistic and antagonistic activities of 615 known drugs or drug candidates against a panel of 6 NHRs: androgen, progesterone, estrogen α/β, and thyroid hormone α/β receptors. Our study revealed that 4.7 and 12.4% compounds have agonistic and antagonistic activities, respectively, against this panel of NHRs. Nonetheless, potent, unintended NHR hits are relatively rare among the known drugs, indicating that such interactions are perhaps not tolerated during drug development. However, we uncovered examples of compounds that unintentionally agonize or antagonize NHRs. In addition, a number of compounds showed multi-NHR activities, suggesting that the cross-talk between multiple NHRs co-operate to elicit in vivo effects. These data highlight the merits of counter screening drug candidate against NHRs during drug discovery/development.

The Nrf2-antioxidant response element signaling pathway and its activation by oxidative stress

A major mechanism in the cellular defense against oxidative or electrophilic stress is activation of the Nrf2-antioxidant response element signaling pathway, which controls the expression of genes whose protein products are involved in the detoxication and elimination of reactive oxidants and electrophilic agents through conjugative reactions and by enhancing cellular antioxidant capacity. At the molecular level, however, the regulatory mechanisms involved in mediating Nrf2 activation are not fully understood. It is well established that Nrf2 activity is controlled, in part, by the cytosolic protein Keap1, but the nature of this pathway and the mechanisms by which Keap1 acts to repress Nrf2 activity remain to be fully characterized and are the topics of discussion in this minireview. In addition, a possible role of the Nrf2-antioxidant response element transcriptional pathway in neuroprotection will also be discussed.

Monitoring the accumulation of fluorescently labeled phospholipids in cell cultures provides an accurate screen for drugs that induce phospholipidosis

A large number of cationic amphiphilic drugs (CADs) are known to cause phospholipidosis (PLD) in vivo. In the present study, we have built upon our previous findings to further qualify the use of a fluorescently labeled phospholipid-based cell-culture assay to detect PLD-inducing drugs. In this paper, we demonstrate that 12 PLD-negative compounds and 11 drugs known to cause PLD in vivo are all correctly identified by using this assay. Interestingly, we found that in cells treated with certain CADs, the fluorescent phospholipid was sequestered in a very specific punctate pattern, which overlapped strongly with the staining pattern seen with a lysosomal marker protein. Our data also show that false positives can be generated with the fluorescence assay when compounds are used at concentrations that cause a >30% decrease in cell number in this assay. Confocal microscopy demonstrated that the staining pattern of fluorescent phospholipids in these cases may be differentiated from those of true positives by the fact that diffuse, rather than punctuate, fluorescence is observed. These studies confirm and expand our previous results showing that the fluorescent phospholipid assay is a highly sensitive, specific tool for detecting PLD-inducing drugs, if care is taken to rule out cytotoxicity-related artifact.

Erratum to "Prediction of non-genotoxic carcinogenesis in rats using changes in gene expression following acute dosing"

Carcinogenicity of chemicals can currently only be evaluated in 2-year rodent bioassays. Therefore, the development of early biomarkers for carcinogenesis would result in substantial savings in time and expense. The current study investigates whether early changes in gene expression may be developed as markers for cancer. Animals were treated for 1 or 5 days with either non-genotoxic carcinogens or non-carcinogens and gene expression was analyzed by quantitative PCR (qPCR).We tested two gene signatures previously reported to detect non-genotoxic carcinogens. Using one gene signature it was confirmed that 3/3 nongenotoxic carcinogens and 2/2 non-carcinogens are correctly identified with data from 1 or 5 days of dosing. In contrast an alternative signature correctly identified 0/3 and 2/3 nongenotoxic carcinogens at 1 and 5 days of treatment, respectively and 2/2 non-carcinogens at both time-points. Additionally, we evaluated a novel panel of putative biomarker genes, from the literature, many of which have roles in cell growth and division, including myc, cdc2 and mcm6. These genes were significantly induced by non-genotoxic carcinogens and not by non-carcinogens. Using the average fold-induction across this panel, 2/3 non-genotoxic carcinogens were detected on both day 1 and day 5. These data support the idea that acute changes in gene expression may provide biomarkers for non-genotoxic carcinogenesis but also highlight interesting differences in the sensitivities of distinct gene signatures.

Interlaboratory evaluation of genomic signatures for predicting carcinogenicity in the rat

The Critical Path Institute recently established the Predictive Safety Testing Consortium, a collaboration between several companies and the U.S. Food and Drug Administration, aimed at evaluating and qualifying biomarkers for a variety of toxicological endpoints. The Carcinogenicity Working Group of the Predictive Safety Testing Consortium has concentrated on sharing data to test the predictivity of two published hepatic gene expression signatures, including the signature by Fielden et al. (2007, Toxicol. Sci. 99, 90-100) for predicting nongenotoxic hepatocarcinogens, and the signature by Nie et al. (2006, Mol. Carcinog. 45, 914-933) for predicting nongenotoxic carcinogens. Although not a rigorous prospective validation exercise, the consortium approach created an opportunity to perform a meta-analysis to evaluate microarray data from short-term rat studies on over 150 compounds. Despite significant differences in study designs and microarray platforms between laboratories, the signatures proved to be relatively robust and more accurate than expected by chance. The accuracy of the Fielden et al. signature was between 63 and 69%, whereas the accuracy of the Nie et al. signature was between 55 and 64%. As expected, the predictivity was reduced relative to internal validation estimates reported under identical test conditions. Although the signatures were not deemed suitable for use in regulatory decision making, they were deemed worthwhile in the early assessment of drugs to aid decision making in drug development. These results have prompted additional efforts to rederive and evaluate a QPCR-based signature using these samples. When combined with a standardized test procedure and prospective interlaboratory validation, the accuracy and potential utility in preclinical applications can be ascertained.

A mutation of Keap1 found in breast cancer impairs its ability to repress Nrf2 activity

Keap1 is the substrate recognition module of a Cullin 3-based E3 ubiquitin ligase. Its primary role is to catalyze the ubiquitylation of the Nrf2 transcription factor. Oxidative stress blocks the E3 ligase activity of Keap1 which stabilizes Nrf2 allowing it to drive the expression of certain antioxidant and drug metabolizing enzymes. A recent study identified a mutation in the Keap1 gene (Keap1C23Y) that is present in breast cancer. Using reporter gene assays we show that Keap1C23Y is impaired in its ability to repress Nrf2 dependent transcription. Unlike wild-type Keap1, we found that Keap1C23Y failed to stimulate the degradation of Nrf2. Co-immunopreciptation experiments showed that Keap1C23Y retains its ability to interact with Nrf2 and Cullin 3. In contrast, we found that Keap1C23Y could not efficiently promote the ubiquitylation of Nrf2, suggesting that its intrinsic biological activity might have been compromised. These results revealed an unexpected role for the N-terminal region of Keap1 in regulating its E3 ligase activity. Importantly, our findings suggest that a paradox exists whereby Nrf2 activity is beneficial in non-malignant cells but in cancer cells it may provide a selective advantage for clonal expansion.

In Vitro Detection of Drug-Induced Phospholipidosis Using Gene Expression and Fluorescent Phospholipid Based Methodologies

Phospholipidosis (PLD) is characterized by the excessive intracellular accumulation of phospholipids. It is well established that a large number of cationic amphiphilic drugs have the potential to induce PLD. In the present study, we describe two facile in vitro methods to determine the PLD-inducing potential of a molecule. The first approach is based on a recent study by (Sawada et al., 2005, Toxicol. Sci. 83, 282-292) in which 17 genes were identified as potential biomarkers of PLD in HepG2 cells. To confirm the utility of this gene panel, we treated HepG2 cells with PLD-positive and -negative compounds and then analyzed gene expression using real-time PCR. Our initial analysis, which used a single dose of each drug, correctly identified five of eight positive compounds and four of four negative compounds. We then increased the doses of the three false negatives (amiodarone, tamoxifen, and loratadine) and found that the changes in gene expression became large enough to correctly identify them as PLD-inducing drugs. Our results suggest that a range of concentrations should be used to increase the accuracy of prediction in this assay. Our second approach utilized a fluorescently labeled phospholipid (LipidTox) which was added to the media of growing HepG2 cells along with compounds positive and negative for PLD. Phospholipid accumulation was determined using confocal microscopy and, more quantitatively, using a 96-well plate assay and a fluorescent plate reader. Using an expanded set of compounds, we show that this assay correctly identified 100% of PLD-positive and -negative compounds. Dose-dependent increases in intracellular fluorescent phospholipid accumulation were observed. We found that this assay was less time consuming, more sensitive, and higher throughput than gene expression analysis. To our knowledge, this study represents the first validation of the use of LipidTox in identifying drugs that can induce PLD.

The carboxy-terminal Neh3 domain of Nrf2 is required for transcriptional activation

Nrf2 is a transcription factor critical for the maintenance of cellular redox homeostasis. We have previously found that Nrf2 is a labile protein, and its activation in cells under stress involves mechanisms leading to its stabilization. As a modular protein, Nrf2 possesses distinct transactivation and DNA binding domains essential for its transcriptional activity. In this study, we found that the C-terminal "Neh3" domain of Nrf2 is also important for its activity. Deletion of the last 16 amino acids of the protein completely abolishes its ability to activate both reporter and endogenous gene expression. Using site-directed mutagenesis, we have identified a stretch of amino acids within this region that are essential for its activity and that are found to be conserved across species and among other members of the CNC-bZIP family. Importantly, deletion of the final 16 amino acids of Nrf2 does not influence its dimerizing capability, DNA binding activity, or subcellular localization, although it does increase the half-life of the protein. In addition, this region was found to be important for interaction with CHD6 (a chromo-ATPase/helicase DNA binding protein) in a yeast two-hybrid screen. RNA interference-mediated knockdown of CHD6 reduced both the basal and tert-butylhydroquinone-inducible expression of NQO1, a prototypical Nrf2 target gene. These data suggest that the Neh3 domain may act as a transactivation domain and that it is possibly involved in interaction with components of the transcriptional apparatus to affect its transcriptional activity.

Nrf2 controls constitutive and inducible expression of ARE-driven genes through a dynamic pathway involving nucleocytoplasmic shuttling by Keap1

Nrf2 regulates the expression of genes encoding antioxidant proteins involved in cellular redox homeostasis. Previous studies have suggested that activation of Nrf2 is mediated by mechanisms promoting its dissociation from Keap1, a cytosolic repressor that acts to sequester the transcription factor in the cytoplasm. As a short-lived protein, Nrf2 is also activated by mechanisms leading to its stabilization in cells under stress, and recent evidence indicates that Keap1 has an active role in the control of its stability. In this report, using immunocytochemistry, cell fractionation, and chromatin immunoprecipitation analyses, we found that Nrf2 is primarily a nuclear protein and that it is expressed and recruited to the chromatin constitutively to drive basal gene expression. Furthermore, we found evidence indicating that Keap1 may repress Nrf2 activity by transiently shuttling into the nucleus to promote its ubiquitylation. The data suggested that the steady-state level of Nrf2 is maintained by a dynamic pathway that balances its constitutive expression with a Keap1-regulated degradation process downstream of its role as a transcriptional activator. We suggest that the stabilization of Nrf2 in cells under stress represents the central regulatory response mediated by mechanisms that interfere with its interaction with Keap1, leading to the induction of antioxidant enzymes important to maintain cellular redox homeostasis.

Contribution of NAD(P)H:quinone oxidoreductase 1 to protection against carcinogenesis, and regulation of its gene by the Nrf2 basic-region leucine zipper and the arylhydrocarbon receptor basic helix-loop-helix transcription factors

NAD(P)H:quinone oxidoreductase 1 (NQO1) is a key enzyme involved in defence against reactive forms of oxygen and inhibition of neoplasia. Under conditions of oxidative stress, expression of NQO1 is induced, and the resulting increase in oxidoreductase protein provides the cell with multiple layers of protection against environmental insults. Firstly, the catalytic activity of NQO1 is directed towards the complete reduction and detoxication of highly reactive quinones. Secondly, the oxidoreductase maintains the endogenous lipid-soluble antioxidants, alpha-tocopherol-hydroquinone and ubiquinol in their reduced and active forms. Thirdly, NQO1 is required for the stabilisation of p53 protein in response to DNA-damaging stimuli, and it thereby influences cell fate decisions. In view of the anticarcinogenic actions of NQO1, an understanding of the mechanisms that govern its expression is desirable. The redox sensitivity of NQO1 transcription occurs through a cis-acting antioxidant response element (ARE) located within the regulatory region of the mouse, rat and human genes. This element recruits the positively acting basic leucine zipper (bZip) transcription factor NF-E2 p45-related factor 2 (Nrf2). Under normal constitutive conditions, Nrf2 associates with the cytoskeletal-binding protein Keap1, which regulates the subcellular distribution of the bZip factor and also targets it for proteasome-dependent degradation. Oxidative stress inhibits the Nrf2-Keap1 interaction, thus promoting nuclear accumulation of the transcription factor and transactivation of NQO1 and other ARE-driven genes. Mouse, rat and human NQO1 can also be induced by planar aromatic hydrocarbons through a cis-acting xenobiotic response element (XRE) located in their gene promoters. The XRE recruits the arylhydrocarbon receptor (AhR) and AhR nuclear translocator. Cross-talk may occur between Nrf2 and AhR, but the details of this process remain to be elucidated.

Identification of a novel Nrf2-regulated antioxidant response element (ARE) in the mouse NAD(P)H:quinone oxidoreductase 1 gene: reassessment of the ARE consensus sequence

NQO1 [NAD(P)H:quinone oxidoreductase 1] has an integral role in cellular responses to oxidative stress. The expression of NQO1 is up-regulated in the mouse following challenge with electrophilic chemicals, in an Nrf2 (NF-E2 p45-related factor 2)-dependent fashion, but the molecular basis for this observation remains unexplained. Through characterization of the murine nqo1 5'-upstream region, we now show that Nrf2 regulates this gene directly via an ARE (antioxidant response element) that lies within a 24 bp region spanning nt -444 to -421. A comprehensive mutation study of this ARE revealed that it does not conform to the currently accepted ARE consensus sequence [(5'-TMAnnRTGAYnnnGCRwwww-3', with essential nucleotides shown in capitals); two cytosine residues (shown in bold in the following sequence) that have been designated 'n' previously because they were thought to be redundant (5'-gagTcA C aGTgAGt C ggCAaaatt-3') have now been found to be essential for enhancer activity; two guanines (also shown in bold) previously regarded as essential for ARE function (5'-gagTcACaGT g AGtCg g CAaaatt-3') have proven to be dispensable]. Examination of wild-type and nrf2 (-/-) mouse embryonic fibroblasts demonstrated that Nrf2 is essential for both constitutive expression of NQO1 and its induction by sulphoraphane. Electrophoretic mobility-shift and chromatin immunoprecipitation assays revealed that Nrf2 associates, in low amounts, with the nqo1 ARE under constitutive conditions, and following sulphoraphane challenge of cells, Nrf2 is recruited to the ARE in substantially greater quantities, as a heterodimer with the small Maf (musculoaponeurotic fibrosarcoma virus) protein, MafK. Also, MafK was found to bind the nqo1 ARE in an Nrf2-independent fashion, and may contribute to transcriptional repression of the oxidoreductase gene. These findings allow a model for transcriptional control of nqo1 through the ARE to be proposed. Furthermore, our results indicate that distinct AREs have differential sequence requirements, and a universally applicable consensus sequence cannot be derived.