Complex effects of sequence variants on lipid levels and coronary artery disease.

Many sequence variants have additive effects on blood lipid levels and, through that, on the risk of coronary artery disease (CAD). We show that variants also have non-additive effects and interact to affect lipid levels as well as affecting variance and correlations. Variance and correlation effects are often signatures of epistasis or gene-environmental interactions. These complex effects can translate into CAD risk. For example, Trp154Ter in FUT2 protects against CAD among subjects with the A1 blood group, whereas it associates with greater risk of CAD in others. His48Arg in ADH1B interacts with alcohol consumption to affect lipid levels and CAD. The effect of variants in TM6SF2 on blood lipids is greatest among those who never eat oily fish but absent from those who often do. This work demonstrates that variants that affect variance of quantitative traits can allow for the discovery of epistasis and interactions of variants with the environment.

Organoid Modeling of the Tumor Immune Microenvironment.

In vitro cancer cultures, including three-dimensional organoids, typically contain exclusively neoplastic epithelium but require artificial reconstitution to recapitulate the tumor microenvironment (TME). The co-culture of primary tumor epithelia with endogenous, syngeneic tumor-infiltrating lymphocytes (TILs) as a cohesive unit has been particularly elusive. Here, an air-liquid interface (ALI) method propagated patient-derived organoids (PDOs) from >100 human biopsies or mouse tumors in syngeneic immunocompetent hosts as tumor epithelia with native embedded immune cells (T, B, NK, macrophages). Robust droplet-based, single-cell simultaneous determination of gene expression and immune repertoire indicated that PDO TILs accurately preserved the original tumor T cell receptor (TCR) spectrum. Crucially, human and murine PDOs successfully modeled immune checkpoint blockade (ICB) with anti-PD-1- and/or anti-PD-L1 expanding and activating tumor antigen-specific TILs and eliciting tumor cytotoxicity. Organoid-based propagation of primary tumor epithelium en bloc with endogenous immune stroma should enable immuno-oncology investigations within the TME and facilitate personalized immunotherapy testing.

DNA demethylase activity maintains intestinal cells in an undifferentiated state following loss of APC.

Although genome-wide hypomethylation is a hallmark of many cancers, roles for active DNA demethylation during tumorigenesis are unknown. Here, loss of the APC tumor suppressor gene causes upregulation of a DNA demethylase system and the concomitant hypomethylation of key intestinal cell fating genes. Notably, this hypomethylation maintained zebrafish intestinal cells in an undifferentiated state that was released upon knockdown of demethylase components. Mechanistically, the demethylase genes are directly activated by Pou5f1 and Cebpß and are indirectly repressed by retinoic acid, which antagonizes Pou5f1 and Cebpß. Apc mutants lack retinoic acid as a result of the transcriptional repression of retinol dehydrogenase l1 via a complex that includes Lef1, Groucho2, Ctbp1, Lsd1, and Corest. Our findings imply a model wherein APC controls intestinal cell fating through a switch in DNA methylation dynamics. Wild-type APC and retinoic acid downregulate demethylase components, thereby promoting DNA methylation of key genes and helping progenitors commit to differentiation.

Blood-based tests for multicancer early detection (PATHFINDER): a prospective cohort study.

BACKGROUND: Multicancer early detection (MCED) blood tests can detect a cancer signal from circulating cell-free DNA (cfDNA). PATHFINDER was a prospective cohort study investigating the feasibility of MCED testing for cancer screening. METHODS: In this prospective cohort study done in oncology and primary care outpatient clinics at seven US health networks, a convenience sample of adults aged 50 years or older without signs or symptoms of cancer consented to MCED testing. We collected blood, analysed cfDNA, and returned results to participants' doctors. If a methylation signature indicative of cancer was detected, predicted cancer signal origin(s) informed diagnostic assessment. The primary outcome was time to, and extent of, diagnostic testing required to confirm the presence or absence of cancer. This trial is registered at ClinicalTrials.gov, NCT04241796, and is completed. FINDINGS: Between Dec 12, 2019, and Dec 4, 2020, we recruited 6662 participants. 4204 (63·5%) of 6621 participants with analysable results were women, 2417 (36·5%) were men, and 6071 (91·7%) were White. A cancer signal was detected in 92 (1·4%) of 6621 participants with analysable results. 35 (38%) participants were diagnosed with cancer (true positives) and 57 (62%) had no cancer diagnosis (false positives). Excluding two participants whose diagnostic assessments began before MCED test results were reported, median time to diagnostic resolution was 79 days (IQR 37-219): 57 days (33-143) in true-positive and 162 days (44-248) in false-positive participants. Most participants had both laboratory tests (26 [79%] of 33 with true-positive results and 50 [88%] of 57 with false-positive results) and imaging (30 [91%] of 33 with true-positive results and 53 [93%] of 57 with false-positive results). Fewer procedures were done in participants with false-positive results (17 [30%] of 57) than true-positive results (27 [82%] of 33) and few had surgery (one with a false-positive result and three with a true-positive result). INTERPRETATION: This study supports the feasibility of MCED screening for cancer and underscores the need for further research investigating the test's clinical utility. FUNDING: GRAIL.

Genetic variants associated with syncope implicate neural and autonomic processes.

AIMS: Syncope is a common and clinically challenging condition. In this study, the genetics of syncope were investigated to seek knowledge about its pathophysiology and prognostic implications. METHODS AND RESULTS: This genome-wide association meta-analysis included 56 071 syncope cases and 890 790 controls from deCODE genetics (Iceland), UK Biobank (United Kingdom), and Copenhagen Hospital Biobank Cardiovascular Study/Danish Blood Donor Study (Denmark), with a follow-up assessment of variants in 22 412 cases and 286 003 controls from Intermountain (Utah, USA) and FinnGen (Finland). The study yielded 18 independent syncope variants, 17 of which were novel. One of the variants, p.Ser140Thr in PTPRN2, affected syncope only when maternally inherited. Another variant associated with a vasovagal reaction during blood donation and five others with heart rate and/or blood pressure regulation, with variable directions of effects. None of the 18 associations could be attributed to cardiovascular or other disorders. Annotation with regard to regulatory elements indicated that the syncope variants were preferentially located in neural-specific regulatory regions. Mendelian randomization analysis supported a causal effect of coronary artery disease on syncope. A polygenic score (PGS) for syncope captured genetic correlation with cardiovascular disorders, diabetes, depression, and shortened lifespan. However, a score based solely on the 18 syncope variants performed similarly to the PGS in detecting syncope risk but did not associate with other disorders. CONCLUSION: The results demonstrate that syncope has a distinct genetic architecture that implicates neural regulatory processes and a complex relationship with heart rate and blood pressure regulation. A shared genetic background with poor cardiovascular health was observed, supporting the importance of a thorough assessment of individuals presenting with syncope.

Evaluating use of changing technologies for rapid next-generation sequencing in pediatrics.

BACKGROUND: Rapid next-generation sequencing (NGS) offers the potential to shorten the diagnostic process and improve the care of acutely ill children. The goal of this study was to report our findings, including benefits and limitations, of a targeted NGS panel and rapid genome sequencing (rGS) in neonatal and pediatric acute clinical care settings. METHODS: Retrospective analysis of patient characteristics, diagnostic yields, turnaround time, and changes in management for infants and children receiving either RapSeq, a targeted NGS panel for 4500+ genes, or rGS, at the University of Utah Hospital and Primary Children's Hospital, from 2015 to 2020. RESULTS: Over a 5-year period, 142 probands underwent rapid NGS: 66 received RapSeq and 76 rGS. Overall diagnostic yield was 39%. In the majority of diagnostic cases, there were one or more changes in clinical care management. Of note, 7% of diagnoses identified by rGS would not have been identified by RapSeq. CONCLUSIONS: Our results indicate that rapid NGS impacts acute pediatric care in real-life clinical settings. Although affected by patient selection criteria, diagnostic yields were similar to those from clinical trial settings. Future studies are needed to determine relative advantages, including cost, turnaround time, and benefits for patients, of each approach in specific clinical circumstances. IMPACT: The use of comprehensive Mendelian gene panels and genome sequencing in the clinical setting allows for early diagnosis of patients in neonatal, pediatric, and cardiac intensive care units and impactful change in management. Diagnoses led to significant changes in management for several patients in lower acuity inpatient units supporting further exploration of the utility of rapid sequencing in these settings. This study reviews the limitations of comparing sequencing platforms in the clinical setting and the variables that should be considered in evaluating diagnostic rates across studies.

Health State Utilities Associated with False-Positive Cancer Screening Results.

INTRODUCTION: Early cancer detection can significantly improve patient outcomes and reduce mortality rates. Novel cancer screening approaches, including multi-cancer early detection tests, have been developed. Cost-utility analyses will be needed to examine their value, and these models require health state utilities. The purpose of this study was to estimate the disutility (i.e., decrease in health state utility) associated with false-positive cancer screening results. METHODS: In composite time trade-off interviews using a 1-year time horizon, UK general population participants valued 10 health state vignettes describing cancer screening with true-negative or false-positive results. Each false-positive vignette described a common diagnostic pathway following a false-positive result suggesting lung, colorectal, breast, or pancreatic cancer. Every pathway ended with a negative result (no cancer detected). The disutility of each false positive was calculated as the difference between the true-negative and each false-positive health state, and because of the 1-year time horizon, each disutility can be interpreted as a quality-adjusted life-year decrement associated with each type of false-positive experience. RESULTS: A total of 203 participants completed interviews (49.8% male; mean age = 42.0 years). The mean (SD) utility for the health state describing a true-negative result was 0.958 (0.065). Utilities for false-positive health states ranged from 0.847 (0.145) to 0.932 (0.059). Disutilities for false positives ranged from - 0.031 to - 0.111 (- 0.041 to - 0.111 for lung cancer; - 0.079 for colorectal cancer; - 0.031 to - 0.067 for breast cancer; - 0.048 to - 0.088 for pancreatic cancer). CONCLUSION: All false-positive results were associated with a disutility. Greater disutility was associated with more invasive follow-up diagnostic procedures, longer duration of uncertainty regarding the eventual diagnosis, and perceived severity of the suspected cancer type. Utility values estimated in this study would be useful for economic modeling examining the value of cancer screening procedures.

Obesity Variants in the GIPR Gene Are not Associated With Risk of Fracture or Bone Mineral Density.

CONTEXT: It is not clear if antagonizing the GIP (glucose-dependent insulinotropic polypeptide) receptor (GIPR) for treatment of obesity is likely to increase the risk of fractures, or to lower bone mineral density (BMD) beyond what is expected with rapid weight loss. OBJECTIVE: The objective of this study was to investigate the risk of fracture and BMD of sequence variants in GIPR that reduce the activity of the GIP receptor and have been associated with reduced body mass index (BMI). METHODS: We analyzed the association of 3 missense variants in GIPR, a common variant, rs1800437 (p.Glu354Gln), and 2 rare variants, rs139215588 (p.Arg190Gln) and rs143430880 (p.Glu288Gly), as well as a burden of predicted loss-of-function (LoF) variants with risk of fracture and with BMD in a large meta-analysis of up to 1.2 million participants. We analyzed associations with fractures at different skeletal sites in the general population: any fractures, hip fractures, vertebral fractures and forearm fractures, and specifically nonvertebral and osteoporotic fractures in postmenopausal women. We also evaluated associations with BMD at the lumbar spine, femoral neck, and total body measured with dual-energy x-ray absorptiometry (DXA), and with BMD estimated from heel ultrasound (eBMD). RESULTS: None of the 3 missense variants in GIPR was significantly associated with increased risk of fractures or with lower BMD. Burden of LoF variants in GIPR was not associated with fractures or with BMD measured with clinically validated DXA, but was associated with eBMD. CONCLUSION: Missense variants in GIPR, or burden of LoF variants in the gene, are not associated with risk of fractures or with lower BMD.

Variant in the synaptonemal complex protein SYCE2 associates with pregnancy loss through effect on recombination.

Two-thirds of all human conceptions are lost, in most cases before clinical detection. The lack of detailed understanding of the causes of pregnancy losses constrains focused counseling for future pregnancies. We have previously shown that a missense variant in synaptonemal complex central element protein 2 (SYCE2), in a key residue for the assembly of the synaptonemal complex backbone, associates with recombination traits. Here we show that it also increases risk of pregnancy loss in a genome-wide association analysis on 114,761 women with reported pregnancy loss. We further show that the variant associates with more random placement of crossovers and lower recombination rate in longer chromosomes but higher in the shorter ones. These results support the hypothesis that some pregnancy losses are due to failures in recombination. They further demonstrate that variants with a substantial effect on the quality of recombination can be maintained in the population.

Start codon variant in LAG3 is associated with decreased LAG-3 expression and increased risk of autoimmune thyroid disease.

Autoimmune thyroid disease (AITD) is a common autoimmune disease. In a GWAS meta-analysis of 110,945 cases and 1,084,290 controls, 290 sequence variants at 225 loci are associated with AITD. Of these variants, 115 are previously unreported. Multiomics analysis yields 235 candidate genes outside the MHC-region and the findings highlight the importance of genes involved in T-cell regulation. A rare 5'-UTR variant (rs781745126-T, MAF = 0.13% in Iceland) in LAG3 has the largest effect (OR = 3.42, P = 2.2 × 10-16) and generates a novel start codon for an open reading frame upstream of the canonical protein translation initiation site. rs781745126-T reduces mRNA and surface expression of the inhibitory immune checkpoint LAG-3 co-receptor on activated lymphocyte subsets and halves LAG-3 levels in plasma among heterozygotes. All three homozygous carriers of rs781745126-T have AITD, of whom one also has two other T-cell mediated diseases, that is vitiligo and type 1 diabetes. rs781745126-T associates nominally with vitiligo (OR = 5.1, P = 6.5 × 10-3) but not with type 1 diabetes. Thus, the effect of rs781745126-T is akin to drugs that inhibit LAG-3, which unleash immune responses and can have thyroid dysfunction and vitiligo as adverse events. This illustrates how a multiomics approach can reveal potential drug targets and safety concerns.

GWAS meta-analysis reveals key risk loci in essential tremor pathogenesis.

Essential tremor (ET) is a prevalent neurological disorder with a largely unknown underlying biology. In this genome-wide association study meta-analysis, comprising 16,480 ET cases and 1,936,173 controls from seven datasets, we identify 12 sequence variants at 11 loci. Evaluating mRNA expression, splicing, plasma protein levels, and coding effects, we highlight seven putative causal genes at these loci, including CA3 and CPLX1. CA3 encodes Carbonic Anhydrase III and carbonic anhydrase inhibitors have been shown to decrease tremors. CPLX1, encoding Complexin-1, regulates neurotransmitter release. Through gene-set enrichment analysis, we identify a significant association with specific cell types, including dopaminergic and GABAergic neurons, as well as biological processes like Rho GTPase signaling. Genetic correlation analyses reveals a positive association between ET and Parkinson's disease, depression, and anxiety-related phenotypes. This research uncovers risk loci, enhancing our knowledge of the complex genetics of this common but poorly understood disorder, and highlights CA3 and CPLX1 as potential therapeutic targets.

Biallelic variants in POPDC2 cause a novel autosomal recessive syndrome presenting with cardiac conduction defects and variable hypertrophic cardiomyopathy.

POPDC2 encodes for the Popeye domain-containing protein 2 which has an important role in cardiac pacemaking and conduction, due in part to its cAMP-dependent binding and regulation of TREK-1 potassium channels. Loss of Popdc2 in mice results in sinus pauses and bradycardia and morpholino knockdown of popdc2 in zebrafish results in atrioventricular (AV) block. We identified bi-allelic variants in POPDC2 in 4 families that presented with a phenotypic spectrum consisting of sinus node dysfunction, AV conduction defects and hypertrophic cardiomyopathy. Using homology modelling we show that the identified POPDC2 variants are predicted to diminish the ability of POPDC2 to bind cAMP. In in vitro electrophysiological studies we demonstrated that, while co-expression of wild-type POPDC2 with TREK-1 increased TREK-1 current density, POPDC2 variants found in the patients failed to increase TREK-1 current density. While patient muscle biopsy did not show clear myopathic disease, it showed significant reduction of the expression of both POPDC1 and POPDC2, suggesting that stability and/or membrane trafficking of the POPDC1-POPDC2 complex is impaired by pathogenic variants in any of the two proteins. Single-cell RNA sequencing from human hearts demonstrated that co-expression of POPDC1 and 2 was most prevalent in AV node, AV node pacemaker and AV bundle cells. Sinoatrial node cells expressed POPDC2 abundantly, but expression of POPDC1 was sparse. Together, these results concur with predisposition to AV node disease in humans with loss-of-function variants in POPDC1 and POPDC2 and presence of sinus node disease in POPDC2, but not in POPDC1 related disease in human. Using population-level genetic data of more than 1 million individuals we showed that none of the familial variants were associated with clinical outcomes in heterozygous state, suggesting that heterozygous family members are unlikely to develop clinical manifestations and therefore might not necessitate clinical follow-up. Our findings provide evidence for POPDC2 as the cause of a novel Mendelian autosomal recessive cardiac syndrome, consistent with previous work showing that mice and zebrafish deficient in functional POPDC2 display sinus and AV node dysfunction.

Integrative common and rare variant analyses provide insights into the genetic architecture of liver cirrhosis.

We report a multi-ancestry genome-wide association study on liver cirrhosis and its associated endophenotypes, alanine aminotransferase (ALT) and γ-glutamyl transferase. Using data from 12 cohorts, including 18,265 cases with cirrhosis, 1,782,047 controls, up to 1 million individuals with liver function tests and a validation cohort of 21,689 cases and 617,729 controls, we identify and validate 14 risk associations for cirrhosis. Many variants are located near genes involved in hepatic lipid metabolism. One of these, PNPLA3 p.Ile148Met, interacts with alcohol intake, obesity and diabetes on the risk of cirrhosis and hepatocellular carcinoma (HCC). We develop a polygenic risk score that associates with the progression from cirrhosis to HCC. By focusing on prioritized genes from common variant analyses, we find that rare coding variants in GPAM associate with lower ALT, supporting GPAM as a potential target for therapeutic inhibition. In conclusion, this study provides insights into the genetic underpinnings of cirrhosis.

Variants at the Interleukin 1 Gene Locus and Pericarditis.

Importance: Recurrent pericarditis is a treatment challenge and often a debilitating condition. Drugs inhibiting interleukin 1 cytokines are a promising new treatment option, but their use is based on scarce biological evidence and clinical trials of modest sizes, and the contributions of innate and adaptive immune processes to the pathophysiology are incompletely understood. Objective: To use human genomics, transcriptomics, and proteomics to shed light on the pathogenesis of pericarditis. Design, Setting, and Participants: This was a meta-analysis of genome-wide association studies of pericarditis from 5 countries. Associations were examined between the pericarditis-associated variants and pericarditis subtypes (including recurrent pericarditis) and secondary phenotypes. To explore mechanisms, associations with messenger RNA expression (cis-eQTL), plasma protein levels (pQTL), and CpG methylation of DNA (ASM-QTL) were assessed. Data from Iceland (deCODE genetics, 1983-2020), Denmark (Copenhagen Hospital Biobank/Danish Blood Donor Study, 1977-2022), the UK (UK Biobank, 1953-2021), the US (Intermountain, 1996-2022), and Finland (FinnGen, 1970-2022) were included. Data were analyzed from September 2022 to August 2023. Exposure: Genotype. Main Outcomes and Measures: Pericarditis. Results: In this genome-wide association study of 4894 individuals with pericarditis (mean [SD] age at diagnosis, 51.4 [17.9] years, 2734 [67.6%] male, excluding the FinnGen cohort), associations were identified with 2 independent common intergenic variants at the interleukin 1 locus on chromosome 2q14. The lead variant was rs12992780 (T) (effect allele frequency [EAF], 31%-40%; odds ratio [OR], 0.83; 95% CI, 0.79-0.87; P = 6.67 × 10-16), downstream of IL1B and the secondary variant rs7575402 (A or T) (EAF, 45%-55%; adjusted OR, 0.89; 95% CI, 0.85-0.93; adjusted P = 9.6 × 10-8). The lead variant rs12992780 had a smaller odds ratio for recurrent pericarditis (0.76) than the acute form (0.86) (P for heterogeneity = .03) and rs7575402 was associated with CpG methylation overlapping binding sites of 4 transcription factors known to regulate interleukin 1 production: PU.1 (encoded by SPI1), STAT1, STAT3, and CCAAT/enhancer-binding protein ß (encoded by CEBPB). Conclusions and Relevance: This study found an association between pericarditis and 2 independent sequence variants at the interleukin 1 gene locus. This finding has the potential to contribute to development of more targeted and personalized therapy of pericarditis with interleukin 1-blocking drugs.

Rapid genome diagnosis of alveolar capillary dysplasia leading to treatment in a child with respiratory and cardiac failure.

Alveolar capillary dysplasia (ACD) is a fatal disorder that typically presents in the neonatal period with refractory hypoxemia and pulmonary hypertension. Lung biopsy is traditionally required to establish the diagnosis. We report a 22-mo-old male who presented with anemia, severe pulmonary hypertension, and right heart failure. He had a complicated hospital course resulting in cardiac arrest and requirement for extracorporeal membrane oxygenation. Computed tomography of the chest showed a heterogenous pattern of interlobular septal thickening and pulmonary edema. The etiology of his condition was unknown, lung biopsy was contraindicated because of his medical fragility, and discussions were held to move to palliative care. Rapid whole-genome sequencing (rWGS) was performed. In 2 d it resulted, revealing a novel FOXF1 gene pathogenic variant that led to the presumptive diagnosis of atypical ACD. Cases of atypical ACD have been reported with survival in patients using medical therapy or lung transplantation. Based on the rWGS diagnosis and more favorable potential of atypical ACD, aggressive medical treatment was pursued. The patient was discharged home after 67 d in the hospital; he is currently doing well more than 30 mo after his initial presentation with only one subsequent hospitalization and no requirement for lung transplantation. Our case reveals the potential for use of rWGS in a critically ill child in which the diagnosis is unknown. rWGS and other advanced genetic tests can guide clinical management and expand our understanding of atypical ACD and other conditions.

HerediGene Population Study IT infrastructure: A model to support genomic research recruitment and precision public health.

The HerediGene Population Study is a large research study focused on identifying new genetic biomarkers for disease prevention, diagnosis, prognosis, and development of new therapeutics. A substantial IT infrastructure evolved to reach enrollment targets and return results to participants. More than 170,000 participants have been enrolled in the study to date, with 5.87% of those whole genome sequenced and 0.46% of those genotyped harboring pathogenic variants. Among other purposes, this infrastructure supports: (1) identifying candidates from clinical criteria, (2) monitoring for qualifying clinical events (e.g., blood draw), (3) contacting candidates, (4) obtaining consent electronically, (5) initiating lab orders, (6) integrating consent and lab orders into clinical workflow, (7) de-identifying samples and clinical data, (8) shipping/transmitting samples and clinical data, (9) genotyping/sequencing samples, (10) and re-identifying and returning results for participants where applicable. This study may serve as a model for similar genomic research and precision public health initiatives.

Genome-wide meta-analysis identifies 93 risk loci and enables risk prediction equivalent to monogenic forms of venous thromboembolism.

We report a genome-wide association study of venous thromboembolism (VTE) incorporating 81,190 cases and 1,419,671 controls sampled from six cohorts. We identify 93 risk loci, of which 62 are previously unreported. Many of the identified risk loci are at genes encoding proteins with functions converging on the coagulation cascade or platelet function. A VTE polygenic risk score (PRS) enabled effective identification of both high- and low-risk individuals. Individuals within the top 0.1% of PRS distribution had a VTE risk similar to homozygous or compound heterozygous carriers of the variants G20210A (c.*97 G > A) in F2 and p.R534Q in F5. We also document that F2 and F5 mutation carriers in the bottom 10% of the PRS distribution had a risk similar to that of the general population. We further show that PRS improved individual risk prediction beyond that of genetic and clinical risk factors. We investigated the extent to which venous and arterial thrombosis share clinical risk factors using Mendelian randomization, finding that some risk factors for arterial thrombosis were directionally concordant with VTE risk (for example, body mass index and smoking) whereas others were discordant (for example, systolic blood pressure and triglyceride levels).

Multidimensional biomarker predicts disease control in response to immunotherapy in recurrent or metastatic head and neck squamous-cell carcinoma.

PURPOSE: Anti-PD-1 therapy provides clinical benefit in 40-50% of patients with relapsed and/or metastatic head and neck squamous cell carcinoma (RM-HNSCC). Selection of anti- PD-1 therapy is typically based on patient PD-L1 immunohistochemistry (IHC) which has low specificity for predicting disease control. Therefore, there is a critical need for a clinical biomarker that will predict clinical benefit to anti-PD-1 treatment with high specificity. METHODS: Clinical treatment and outcomes data for 103 RM-HNSCC patients were paired with RNA-sequencing data from formalin-fixed patient samples. Using logistic regression methods, we developed a novel biomarker classifier based on expression patterns in the tumor immune microenvironment to predict disease control with monotherapy PD-1 inhibitors (pembrolizumab and nivolumab). The performance of the biomarker was internally validated using out-of-bag methods. RESULTS: The biomarker significantly predicted disease control (65% in predicted non-progressors vs. 17% in predicted progressors, p < 0.001) and was significantly correlated with overall survival (OS; p = 0.004). In addition, the biomarker outperformed PD-L1 IHC across numerous metrics including sensitivity (0.79 vs 0.64, respectively; p = 0.005) and specificity (0.70 vs 0.61, respectively; p = 0.009). CONCLUSION: This novel assay uses tumor immune microenvironment expression data to predict disease control and OS with high sensitivity and specificity in patients with RM-HNSCC treated with anti-PD-1 monotherapy.

Rare variants with large effects provide functional insights into the pathology of migraine subtypes, with and without aura.

Migraine is a complex neurovascular disease with a range of severity and symptoms, yet mostly studied as one phenotype in genome-wide association studies (GWAS). Here we combine large GWAS datasets from six European populations to study the main migraine subtypes, migraine with aura (MA) and migraine without aura (MO). We identified four new MA-associated variants (in PRRT2, PALMD, ABO and LRRK2) and classified 13 MO-associated variants. Rare variants with large effects highlight three genes. A rare frameshift variant in brain-expressed PRRT2 confers large risk of MA and epilepsy, but not MO. A burden test of rare loss-of-function variants in SCN11A, encoding a neuron-expressed sodium channel with a key role in pain sensation, shows strong protection against migraine. Finally, a rare variant with cis-regulatory effects on KCNK5 confers large protection against migraine and brain aneurysms. Our findings offer new insights with therapeutic potential into the complex biology of migraine and its subtypes.