Molecular and clinical analyses of PHF6 mutant myeloid neoplasia provide their pathogenesis and therapeutic targeting.

PHF6 mutations (PHF6MT) are identified in various myeloid neoplasms (MN). However, little is known about the precise function and consequences of PHF6 in MN. Here we show three main findings in our comprehensive genomic and proteomic study. Firstly, we show a different pattern of genes correlating with PHF6MT in male and female cases. When analyzing male and female cases separately, in only male cases, RUNX1 and U2AF1 are co-mutated with PHF6. In contrast, female cases reveal co-occurrence of ASXL1 mutations and X-chromosome deletions with PHF6MT. Next, proteomics analysis reveals a direct interaction between PHF6 and RUNX1. Both proteins co-localize in active enhancer regions that define the context of lineage differentiation. Finally, we demonstrate a negative prognostic role of PHF6MT, especially in association with RUNX1. The negative effects on survival are additive as PHF6MT cases with RUNX1 mutations have worse outcomes when compared to cases carrying single mutation or wild-type.

Exome-wide association study to identify rare variants influencing COVID-19 outcomes: Results from the Host Genetics Initiative.

Host genetics is a key determinant of COVID-19 outcomes. Previously, the COVID-19 Host Genetics Initiative genome-wide association study used common variants to identify multiple loci associated with COVID-19 outcomes. However, variants with the largest impact on COVID-19 outcomes are expected to be rare in the population. Hence, studying rare variants may provide additional insights into disease susceptibility and pathogenesis, thereby informing therapeutics development. Here, we combined whole-exome and whole-genome sequencing from 21 cohorts across 12 countries and performed rare variant exome-wide burden analyses for COVID-19 outcomes. In an analysis of 5,085 severe disease cases and 571,737 controls, we observed that carrying a rare deleterious variant in the SARS-CoV-2 sensor toll-like receptor TLR7 (on chromosome X) was associated with a 5.3-fold increase in severe disease (95% CI: 2.75-10.05, p = 5.41x10-7). This association was consistent across sexes. These results further support TLR7 as a genetic determinant of severe disease and suggest that larger studies on rare variants influencing COVID-19 outcomes could provide additional insights.

Identification of potential COVID-19 treatment compounds which inhibit SARS Cov2 prototypic, Delta and Omicron variant infection.

Recurrent waves of COVID19 remain a major global health concern. Repurposing either FDA-approved or clinically advanced drug candidates can save time and effort required for validating the safety profile and FDA approval. However, the selection of appropriate screening approaches is key to identifying novel candidate drugs with a higher probability of clinical success. Here, we report a rapid, stratified two-step screening approach using pseudovirus entry inhibition assay followed by an infectious prototypic SARS CoV2 cytotoxic effect inhibition assay in multiple cell lines. Using this approach, we screened a library of FDA-approved and clinical-stage drugs and identified four compounds, apilimod, berbamine, cepharanthine and (S)-crizotinib which potently inhibited SARS CoV2-induced cell death. Importantly, these drugs exerted similar inhibitory effect on the delta and omicron variants although they replicated less efficiently than the prototypic strain. Apilimod is currently under clinical trial (NCT04446377) for COVID19 supporting the validity and robustness of our screening approach.

Elevated plasma levels of CXCL16 in severe COVID-19 patients.

Genome-wide association studies have recently identified 3p21.31, with lead variant pointing to the CXCR6 gene, as the strongest thus far reported susceptibility risk locus for severe manifestation of COVID-19. In order the determine its role, we measured plasma levels of Chemokine (C-X-C motif) ligand 16 (CXCL16) in the plasma of COVID-19 hospitalized patients. CXCL16 interacts with CXCR6 promoting chemotaxis or cell adhesion. The CXCR6/CXCL16 axis mediates homing of T cells to the lungs in disease and hyper-expression is associated with localised cellular injury. To characterize the CXCR6/CXCL16 axis in the pathogenesis of severe COVID-19, plasma concentrations of CXCL16 collected at baseline from 115 hospitalized COVID-19 patients participating in ODYSSEY COVID-19 clinical trial were assessed together with a set of controls. We report elevated levels of CXCL16 in a cohort of COVID-19 hospitalized patients. Specifically, we report significant elevation of CXCL16 plasma levels in association with severity of COVID-19 (as defined by WHO scale) (P-value < 0.02). Our current study is the largest thus far study reporting CXCL16 levels in COVID-19 hospitalized patients (with whole-genome sequencing data available). The results further support the significant role of the CXCR6/CXCL16 axis in the immunopathogenesis of severe COVID-19 and warrants further studies to understand which patients would benefit most from targeted treatments.

PBRM1 loss in kidney cancer unbalances the proximal tubule master transcription factor hub to repress proximal tubule differentiation.

PBRM1, a subunit of the PBAF coactivator complex that transcription factors use to activate target genes, is genetically inactivated in almost all clear cell renal cell cancers (RCCs). Using unbiased proteomic analyses, we find that PAX8, a master transcription factor driver of proximal tubule epithelial fates, recruits PBRM1/PBAF. Reverse analyses of the PAX8 interactome confirm recruitment specifically of PBRM1/PBAF and not functionally similar BAF. More conspicuous in the PAX8 hub in RCC cells, however, are corepressors, which functionally oppose coactivators. Accordingly, key PAX8 target genes are repressed in RCC versus normal kidneys, with the loss of histone lysine-27 acetylation, but intact lysine-4 trimethylation, activation marks. Re-introduction of PBRM1, or depletion of opposing corepressors using siRNA or drugs, redress coregulator imbalance and release RCC cells to terminal epithelial fates. These mechanisms thus explain RCC resemblance to the proximal tubule lineage but with suppression of the late-epithelial program that normally terminates lineage-precursor proliferation.

Assessing the potential correlation of polymorphisms in the IL6R with relative IL6 elevation in severely ill COVID-19 patients'.

BACKGROUND: Elevated Interleukin-6 (IL-6) may play an important role in the pathophysiology of COVID-19 yet attenuated response is not seen across all severe patients. We aimed to determine the effect of IL-6 baseline levels and other clinical variables on mortality and outcomes in hospitalized COVID-19 patients as well as to explore genetic variants associated with attenuated IL-6 response. METHODS: Baseline IL-6 cytokine levels were measured in hospitalized patients participating in ongoing ODYSSEY phase 3 randomized study of tradipitant and placebo in hospitalized patients with severe COVID-19 who are receiving supplemental oxygen support. Furthermore blood samples for whole genome sequencing analysis were collected from 150 participants. RESULTS: We report significantly elevated IL-6 in COVID-19 infected hospitalized patients, n = 100 (p-value < 0.0001) when compared to controls n = 324. We also report a significantly increased level of IL-6 (p-value < 0.01) between the severe and mild COVID-19 patients with severity defined on a WHO scale. Excessive IL-6 plasma levels correlate with higher mortality (p-value 0.001). Additionally, based on our classification analysis, combination of IL-6 elevation and high levels of serum glucose can identify highest risk-group of COVID19 patients. Furthermore, we explore the role of genetic regulatory variants affecting baseline IL-6 levels specifically in COVID-19 patients. We have directly tested the association between variants in the IL6 and IL6R gene region and IL6 plasma levels. We provide results for a common IL-6 variant previously associated with pneumonia, rs1800795, and rs2228145 that was previously shown to affect IL-6 plasma levels, as well as report on novel variants associated with IL-6 plasma levels detected in our study patients. CONCLUSIONS: While it is unlikely that "cytokine storm" is the norm in severe COVID19, baseline elevations above 150 pg/ml may be associated with worst outcomes and as such may warrant treatment considerations. So far no clinical studies used IL-6 baseline assessment to stratify the patient population participating in clinical studies. We believe that careful examination and interpretation of the IL-6 levels and genetic variants can help to determine a patient population with a potentially very robust clinical response to IL-6 inhibition. TRIAL REGISTRATION: Clinicaltrials.gov: NCT04326426.

Distinctive and common features of moderate aplastic anaemia.

The therapy algorithm for severe aplastic anaemia (sAA) is established but moderate AA (mAA), which likely reflects a more diverse pathogenic mechanism, often represents a treatment/management conundrum. A cohort of AA patients (n = 325) was queried for those with non-severe disease using stringent criteria including bone marrow hypocellularity and chronic persistence of moderately depressed blood counts. As a result, we have identified and analyzed pathological and clinical features in 85 mAA patients. Progression to sAA and direct clonal evolution (paroxysmal nocturnal haemoglobinuria/acute myeloid leukaemia; PNH/AML) occurred in 16%, 11% and 1% of mAA cases respectively. Of the mAA patients who received immunosuppressive therapy, 67% responded irrespective of time of initiation of therapy while conservatively managed patients showed no spontaneous remissions. Genomic analysis of mAA identified evidence of clonal haematopoiesis with both persisting and remitting patterns at low allelic frequencies; with more pronounced mutational burden in sAA. Most of the mAA patients have autoimmune pathogenesis similar to those with sAA, but mAA contains a mix of patients with diverse aetiologies. Although progression rates differed between mAA and sAA (P = 0·003), cumulative incidences of mortalities were only marginally different (P = 0·095). Our results provide guidance for diagnosis/management of mAA, a condition for which no current standard of care is established.

Therapy-related acute lymphoblastic leukemia is a distinct entity with adverse genetic features and clinical outcomes.

Patients with therapy-related acute lymphoblastic leukemia (t-ALL) represent a small subset of acute lymphoblastic leukemia (ALL) patients who received genotoxic therapy (ie, chemotherapy or radiation) for a prior malignancy. These patients should be distinguished from patients with de novo ALL (dn-ALL) and ALL patients who have a history of prior malignancy but have not received cytotoxic therapies in the past (acute lymphoblastic leukemia with prior malignancy [pm-ALL]). We report a retrospective multi-institutional study of patients with t-ALL (n = 116), dn-ALL (n = 100), and pm-ALL (n = 20) to investigate the impact of prior cytotoxic therapies on clinical outcomes. Compared with patients with pm-ALL, t-ALL patients had a significantly shorter interval between the first malignancy and ALL diagnosis and a higher frequency of poor-risk cytogenetic features, including KMT2A rearrangements and myelodysplastic syndrome-like abnormalities (eg, monosomal karyotype). We observed a variety of mutations among t-ALL patients, with the majority of patients exhibiting mutations that were more common with myeloid malignancies (eg, DNMT3A, RUNX1, ASXL1), whereas others had ALL-type mutations (eg, CDKN2A, IKZF1). Median overall survival was significantly shorter in the t-ALL cohort compared with patients with dn-ALL or pm-ALL. Patients who were eligible for hematopoietic cell transplantation had improved long-term survival. Collectively, our results support t-ALL as a distinct entity based on its biologic and clinical features.

Non-del(5q) myelodysplastic syndromes-associated loci detected by SNP-array genome-wide association meta-analysis.

Myelodysplastic syndromes (MDS) are hematopoietic stem cell malignancies. Known predisposing factors to adult MDS include rare germline mutations, cytotoxic therapy, age-related clonal hematopoiesis, and autoimmune or chronic inflammatory disorders. To date, no published studies characterizing MDS-associated germline susceptibility polymorphisms exist. We performed a genome-wide association study of 2 sample sets (555 MDS cases vs 2964 control subjects; 352 MDS cases vs 2640 control subjects) in non-del(5q) MDS cases of European genomic ancestry. Meta-analysis identified 8 MDS-associated loci at 1q31.1 (PLA2G4A), 3p14.1 (FAM19A4), 5q21.3 (EFNA5), 6p21.33, 10q23.1 (GRID1), 12q24.32, 15q26.1, and 20q13.12 (EYA2) that approached genome-wide significance. Gene expression for 5 loci that mapped within or near genes was significantly upregulated in MDS bone marrow cells compared with those of control subjects (P < .01). Higher PLA2G4A expression and lower EYA2 expression were associated with poorer overall survival (P = .039 and P = .037, respectively). Higher PLA2G4A expression is associated with mutations in NRAS (P < .001), RUNX1 (P = .012), ASXL1 (P = .007), and EZH2 (P = .038), all of which are known to contribute to MDS development. EYA2 expression was an independently favorable risk factor irrespective of age, sex, and Revised International Scoring System score (relative risk, 0.67; P = .048). Notably, these genes have regulatory roles in innate immunity, a critical driver of MDS pathogenesis. EYA2 overexpression induced innate immune activation, whereas EYA2 inhibition restored colony-forming potential in primary MDS cells indicative of hematopoietic restoration and possible clinical relevance. In conclusion, among 8 suggestive MDS-associated loci, 5 map to genes upregulated in MDS with functional roles in innate immunity and potential biological relevance to MDS.

Invariant patterns of clonal succession determine specific clinical features of myelodysplastic syndromes.

Myelodysplastic syndromes (MDS) arise in older adults through stepwise acquisitions of multiple somatic mutations. Here, analyzing 1809 MDS patients, we infer clonal architecture by using a stringent, the single-cell sequencing validated PyClone bioanalytic pipeline, and assess the position of the mutations within the clonal architecture. All 3,971 mutations are grouped based on their rank in the deduced clonal hierarchy (dominant and secondary). We evaluated how they affect the resultant morphology, progression, survival and response to therapies. Mutations of SF3B1, U2AF1, and TP53 are more likely to be dominant, those of ASXL1, CBL, and KRAS are secondary. Among distinct combinations of dominant/secondary mutations we identified 37 significant relationships, of which 12 affect clinical phenotypes, 5 cooperatively associate with poor prognosis. They also predict response to hypomethylating therapies. The clonal hierarchy has distinct ranking and the resultant invariant combinations of dominant/secondary mutations yield novel insights into the specific clinical phenotype of MDS.

Distinct clinical and biological implications of CUX1 in myeloid neoplasms.

Somatic mutations of the CUT-like homeobox 1 (CUX1) gene (CUX1 MT) can be found in myeloid neoplasms (MNs), in particular, in myelodysplastic syndromes (MDSs). The CUX1 locus is also deleted in 3 of 4 MN cases with -7/del(7q). A cohort of 1480 MN patients was used to characterize clinical features and clonal hierarchy associated with CUX1 MT and CUX1 deletions (CUX1 DEL) and to analyze their functional consequences in vitro. CUX1 MT were present in 4% of chronic MNs. CUX1 DEL were preferentially found in advanced cases (6%). Most MDS and acute myeloid leukemia (AML) patients with -7/del(7q) and up to 15% of MDS patients and 5% of AML patients diploid for the CUX1 locus exhibited downmodulated CUX1 expression. In 75% of mutant cases, CUX1 MT were heterozygous, whereas microdeletions and homozygous and compound-heterozygous mutations were less common. CUX MT/DEL were associated with worse survival compared with CUX1 WT Within the clonal hierarchy, 1 of 3 CUX1 MT served as founder events often followed by secondary BCOR and ASXL1 subclonal hits, whereas TET2 was the most common ancestral lesion, followed by subclonal CUX1 MT Comet assay of patients' bone marrow progenitor cells and leukemic cell lines performed in various experimental conditions revealed that frameshift mutations, hemizygous deletions, or experimental CUX1 knockdown decrease the repair of oxidized bases. These functional findings may explain why samples with either CUX1 MT or low CUX1 expression coincided with significantly higher numbers of somatic hits by whole-exome sequencing. Our findings implicate the DNA repair dysfunction resulting from CUX1 lesions in the pathogenesis of MNs, in which they lead to a mutator phenotype.

Misidentification of MLL3 and other mutations in cancer due to highly homologous genomic regions.

The MLL3 gene has been shown to be recurrently mutated in many malignancies including in families with acute myeloid leukemia. We demonstrate that many MLL3 variant calls made by exome sequencing are false positives due to misalignment to homologous regions, including a region on chr21, and can only be validated by long-range PCR. Numerous other recurrently mutated genes reported in COSMIC and TCGA databases have pseudogenes and cannot also be validated by conventional short read-based sequencing approaches. Genome-wide identification of pseudogene regions demonstrates that frequency of these homologous regions is increased with sequencing read lengths below 200 bps. To enable identification of poor quality sequencing variants in prospective studies, we generated novel genome-wide maps of regions with poor mappability that can be used in variant calling algorithms. Taken together, our findings reveal that pseudogene regions are a source of false-positive mutations in cancers.

Gene-centric functional dissection of human genetic variation uncovers regulators of hematopoiesis.

Genome-wide association studies (GWAS) have identified thousands of variants associated with human diseases and traits. However, the majority of GWAS-implicated variants are in non-coding regions of the genome and require in depth follow-up to identify target genes and decipher biological mechanisms. Here, rather than focusing on causal variants, we have undertaken a pooled loss-of-function screen in primary hematopoietic cells to interrogate 389 candidate genes contained in 75 loci associated with red blood cell traits. Using this approach, we identify 77 genes at 38 GWAS loci, with most loci harboring 1-2 candidate genes. Importantly, the hit set was strongly enriched for genes validated through orthogonal genetic approaches. Genes identified by this approach are enriched in specific and relevant biological pathways, allowing regulators of human erythropoiesis and modifiers of blood diseases to be defined. More generally, this functional screen provides a paradigm for gene-centric follow up of GWAS for a variety of human diseases and traits.

Development of a Cx46 Targeting Strategy for Cancer Stem Cells.

Gap-junction-mediated cell-cell communication enables tumor cells to synchronize complex processes. We previously found that glioblastoma cancer stem cells (CSCs) express higher levels of the gap junction protein Cx46 compared to non-stem tumor cells (non-CSCs) and that this was necessary and sufficient for CSC maintenance. To understand the mechanism underlying this requirement, we use point mutants to disrupt specific functions of Cx46 and find that Cx46-mediated gap-junction coupling is critical for CSCs. To develop a Cx46 targeting strategy, we screen a clinically relevant small molecule library and identify clofazimine as an inhibitor of Cx46-specific cell-cell communication. Clofazimine attenuates proliferation, self-renewal, and tumor growth and synergizes with temozolomide to induce apoptosis. Although clofazimine does not cross the blood-brain barrier, the combination of clofazimine derivatives optimized for brain penetrance with standard-of-care therapies may target glioblastoma CSCs. Furthermore, these results demonstrate the importance of targeting cell-cell communication as an anti-cancer therapy.

Subclonal STAT3 mutations solidify clonal dominance.

T large granular lymphocyte leukemia (T-LGLL) is a clonal lymphoproliferative disorder that can arise in the context of pathologic or physiologic cytotoxic T-cell (CTL) responses. STAT3 mutations are often absent in typical T-LGLL, suggesting that in a significant fraction of patients, antigen-driven expansion alone can maintain LGL clone persistence. We set out to determine the relationship between activating STAT3 hits and CTL clonal selection at presentation and in response to therapy. Thus, a group of patients with T-LGLL were serially subjected to deep next-generation sequencing (NGS) of the T-cell receptor (TCR) Vß complementarity-determining region 3 (CDR3) and STAT3 to recapitulate clonal hierarchy and dynamics. The results of this complex analysis demonstrate that STAT3 mutations produce either a sweeping or linear subclone within a monoclonal CTL population either early or during the course of disease. Therapy can extinguish a LGL clone, silence it, or adapt mechanisms to escape elimination. LGL clones can persist on elimination of STAT3 subclones, and alternate STAT3-negative CTL clones can replace therapy-sensitive CTL clones. LGL clones can evolve and are fueled by a nonextinguished antigenic drive. STAT3 mutations can accelerate this process or render CTL clones semiautonomous and not reliant on physiologic stimulation.