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 AT1 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-51), DCLRE1A (p = 8.4 × 10-13), and HELB (p = 5.7 × 10-7) with later menopause and TOP3A (p = 7.6 × 10-8) and CLPB (p = 8.1 × 10-7) with earlier menopause. Two additional genes are suggestive: RAD54L (p = 2.4 × 10-6) with later menopause and HROB (p = 2.9 × 10-6) 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.

Advancing human genetics research and drug discovery through exome sequencing of the UK Biobank.

The UK Biobank Exome Sequencing Consortium (UKB-ESC) is a private-public partnership between the UK Biobank (UKB) and eight biopharmaceutical companies that will complete the sequencing of exomes for all ~500,000 UKB participants. Here, we describe the early results from ~200,000 UKB participants and the features of this project that enabled its success. The biopharmaceutical industry has increasingly used human genetics to improve success in drug discovery. Recognizing the need for large-scale human genetics data, as well as the unique value of the data access and contribution terms of the UKB, the UKB-ESC was formed. As a result, exome data from 200,643 UKB enrollees are now available. These data include ~10 million exonic variants-a rich resource of rare coding variation that is particularly valuable for drug discovery. The UKB-ESC precompetitive collaboration has further strengthened academic and industry ties and has provided teams with an opportunity to interact with and learn from the wider research community.

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.

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.

Author Correction: Phenotypes associated with genes encoding drug targets are predictive of clinical trial side effects.

In the original version of this article, there were errors in the labelling of the colours in the key of Figure 2, whereby the labeling of the third and fourth of the four colours was reversed. This has been corrected in both the PDF and HTML versions of the article.

Phenotypes associated with genes encoding drug targets are predictive of clinical trial side effects.

Only a small fraction of early drug programs progress to the market, due to safety and efficacy failures, despite extensive efforts to predict safety. Characterizing the effect of natural variation in the genes encoding drug targets should present a powerful approach to predict side effects arising from drugging particular proteins. In this retrospective analysis, we report a correlation between the organ systems affected by genetic variation in drug targets and the organ systems in which side effects are observed. Across 1819 drugs and 21 phenotype categories analyzed, drug side effects are more likely to occur in organ systems where there is genetic evidence of a link between the drug target and a phenotype involving that organ system, compared to when there is no such genetic evidence (30.0 vs 19.2%; OR = 1.80). This result suggests that human genetic data should be used to predict safety issues associated with drug targets.

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.

A rare missense variant in NR1H4 associates with lower cholesterol levels.

Searching for novel sequence variants associated with cholesterol levels is of particular interest due to the causative role of non-HDL cholesterol levels in cardiovascular disease. Through whole-genome sequencing of 15,220 Icelanders and imputation of the variants identified, we discovered a rare missense variant in NR1H4 (R436H) associating with lower levels of total cholesterol (effect = -0.47 standard deviations or -0.55 mmol L-1, p = 4.21 × 10-10, N = 150,211). Importantly, NR1H4 R436H also associates with lower levels of non-HDL cholesterol and, consistent with this, protects against coronary artery disease. NR1H4 encodes FXR that regulates bile acid homeostasis, however, we do not detect a significant association between R436H and biological markers of liver function. Transcriptional profiling of hepatocytes carrying R436H shows that it is not a loss-of-function variant. Rather, we observe changes in gene expression compatible with effects on lipids. These findings highlight the role of FXR in regulation of cholesterol levels in humans.