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.

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.

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.

Genetic abnormalities in myelodysplasia and secondary acute myeloid leukemia: impact on outcome of stem cell transplantation.

Genetic alterations, including mutations and copy-number alterations, are central to the pathogenesis of myelodysplastic syndromes and related diseases (myelodysplasia), but their roles in allogeneic stem cell transplantation have not fully been studied in a large cohort of patients. We enrolled 797 patients who had been diagnosed with myelodysplasia at initial presentation and received transplantation via the Japan Marrow Donor Program. Targeted-capture sequencing was performed to identify mutations in 69 genes, together with copy-number alterations, whose effects on transplantation outcomes were investigated. We identified 1776 mutations and 927 abnormal copy segments among 617 patients (77.4%). In multivariate modeling using Cox proportional-hazards regression, genetic factors explained 30% of the total hazards for overall survival; clinical characteristics accounted for 70% of risk. TP53 and RAS-pathway mutations, together with complex karyotype (CK) as detected by conventional cytogenetics and/or sequencing-based analysis, negatively affected posttransplant survival independently of clinical factors. Regardless of disease subtype, TP53-mutated patients with CK were characterized by unique genetic features and associated with an extremely poor survival with frequent early relapse, whereas outcomes were substantially better in TP53-mutated patients without CK. By contrast, the effects of RAS-pathway mutations depended on disease subtype and were confined to myelodysplastic/myeloproliferative neoplasms (MDS/MPNs). Our results suggest that TP53 and RAS-pathway mutations predicted a dismal prognosis, when associated with CK and MDS/MPNs, respectively. However, for patients with mutated TP53 or CK alone, long-term survival could be obtained with transplantation. Clinical sequencing provides vital information for accurate prognostication in transplantation.

Transcriptomic rationale for synthetic lethality-targeting ERCC1 and CDKN1A in chronic myelomonocytic leukaemia.

Despite the absence of mutations in the DNA repair machinery in myeloid malignancies, the advent of high-throughput sequencing and discovery of splicing and epigenetics defects in chronic myelomonocytic leukaemia (CMML) prompted us to revisit a pathogenic role for genes involved in DNA damage response. We screened for misregulated DNA repair genes by enhanced RNA-sequencing on bone marrow from a discovery cohort of 27 CMML patients and 9 controls. We validated 4 differentially expressed candidates in CMML CD34+ bone marrow selected cells and in an independent cohort of 74 CMML patients, mutationally contextualized by targeted sequencing, and assessed their transcriptional behavior in 70 myelodysplastic syndrome, 66 acute myeloid leukaemia and 25 chronic myeloid leukaemia cases. We found BAP1 and PARP1 down-regulation to be specific to CMML compared with other related disorders. Chromatin-regulator mutated cases showed decreased BAP1 dosage. We validated a significant over-expression of the double strand break-fidelity genes CDKN1A and ERCC1, independent of promoter methylation and associated with chemorefractoriness. In addition, patients bearing mutations in the splicing component SRSF2 displayed numerous aberrant splicing events in DNA repair genes, with a quantitative predominance in the single strand break pathway. Our results highlight potential targets in this disease, which currently has few therapeutic options.

Rational management approach to pure red cell aplasia.

Pure red cell aplasia is an orphan disease, and as such lacks rationally established standard therapies. Most cases are idiopathic; a subset is antibody-mediated. There is overlap between idiopathic cases and those with T-cell large granular lymphocytic leukemia, hypogammaglobulinemia, and low-grade lymphomas. In each of the aforementioned, the pathogenetic mechanisms may involve autoreactive cytotoxic responses. We selected 62 uniformly diagnosed pure red cell aplasia patients and analyzed their pathophysiologic features and responsiveness to rationally applied first-line and salvage therapies in order to propose diagnostic and therapeutic algorithms that may be helpful in guiding the management of prospective patients, 52% of whom were idiopathic, while the others involved large granular lymphocytic leukemia, thymoma, and B-cell dyscrasia. T-cell-mediated responses ranged between a continuum from polyclonal to monoclonal (as seen in large granular lymphocytic leukemia). During a median observation period of 40 months, patients received a median of two different therapies to achieve remission. Frequently used therapy included calcineurin-inhibitors with a steroid taper yielding a first-line overall response rate of 76% (53/70). Oral cyclophosphamide showed activity, albeit lower than that produced by cyclosporine. Intravenous immunoglobulins were effective both in parvovirus patients and in hypogammaglobulinemia cases. In salvage settings, alemtuzumab is active, particularly in large granular lymphocytic leukemia-associated cases. Other potentially useful salvage options include rituximab, anti-thymocyte globulin and bortezomib. The workup of acquired pure red cell aplasia should include investigations of common pathological associations. Most effective therapies are directed against T-cell-mediated immunity, and therapeutic choices need to account for associated conditions that may help in choosing alternative salvage agents, such as intravenous immunoglobulin, alemtuzumab and bortezomib.

Molecular features of early onset adult myelodysplastic syndrome.

Myelodysplastic syndromes are typically diseases of older adults. Patients in whom the onset is early may have distinct molecular and clinical features or reflect a demographic continuum. The identification of differences between "early onset" patients and those diagnosed at a traditional age has the potential to advance understanding of the pathogenesis of myelodysplasia and may lead to formation of distinct morphological subcategories. We studied a cohort of 634 patients with various subcategories of myelodysplastic syndrome and secondary acute myeloid leukemia, stratifying them based on age at presentation and clinical parameters. We then characterized molecular abnormalities detected by next-generation deep sequencing of 60 genes that are commonly mutated in myeloid malignancies. The number of mutations increased linearly with age and on average, patients >50 years of age had more mutations. TET2, SRSF2, and DNMT3A were more commonly mutated in patients >50 years old compared to patients ≤50 years old. In general, patients >50 years of age also had more mutations in spliceosomal, epigenetic modifier, and RAS gene families. Although there are age-related differences in molecular features among patients with myelodysplasia, most notably in the incidence of SRSF2 mutations, our results suggest that patients ≤50 years old belong to a disease continuum with a distinct pattern of early onset ancestral events.

Genetic alterations of the cohesin complex genes in myeloid malignancies.

Somatic cohesin mutations have been reported in myelodysplastic syndrome (MDS) and acute myeloid leukemia (AML). To account for the morphologic and cytogenetic diversity of these neoplasms, a well-annotated cohort of 1060 patients with myeloid malignancies including MDS (n = 386), myeloproliferative neoplasms (MPNs) (n = 55), MDS/MPNs (n = 169), and AML (n = 450) were analyzed for cohesin gene mutational status, gene expression, and therapeutic and survival outcomes. Somatic cohesin defects were detected in 12% of patients with myeloid malignancies, whereas low expression of these genes was present in an additional 15% of patients. Mutations of cohesin genes were mutually exclusive and mostly resulted in predicted loss of function. Patients with low cohesin gene expression showed similar expression signatures as those with somatic cohesin mutations. Cross-sectional deep-sequencing analysis for clonal hierarchy demonstrated STAG2, SMC3, and RAD21 mutations to be ancestral in 18%, 18%, and 47% of cases, respectively, and each expanded to clonal dominance concordant with disease transformation. Cohesin mutations were significantly associated with RUNX1, Ras-family oncogenes, and BCOR and ASXL1 mutations and were most prevalent in high-risk MDS and secondary AML. Cohesin defects were associated with poor overall survival (27.2 vs 40 months; P = .023), especially in STAG2 mutant MDS patients surviving >12 months (median survival 35 vs 50 months; P = .017).

Patterns of missplicing due to somatic U2AF1 mutations in myeloid neoplasms.

Recently, recurrent mutations of spliceosomal genes were frequently identified in myeloid malignancies, as well as other types of cancers. One of these spliceosomal genes, U2AF1, was affected by canonical somatic mutations in aggressive type of myeloid malignancies. We hypothesized that U2AF1 mutations causes defects of splicing (missplicing) in specific genes and that such misspliced genes might be important in leukemogenesis. We analyzed RNA deep sequencing to compare splicing patterns of 201 837 exons between the cases with U2AF1 mutations (n = 6) and wild type (n = 14). We identified different alternative splicing patterns in 35 genes comparing cells with mutant and wild-type U2AF1. U2AF1 mutations are associated with abnormal splicing of genes involved in functionally important pathways, such as cell cycle progression and RNA processing. In addition, many of these genes are somatically mutated or deleted in various cancers. Of note is that the alternative splicing patterns associated with U2AF1 mutations were associated with specific sequence signals at the affected splice sites. These novel observations support the hypothesis that U2AF1 mutations play a significant role in myeloid leukemogenesis due to selective missplicing of tumor-associated genes.

Deep sequencing of the T-cell receptor repertoire in CD8+ T-large granular lymphocyte leukemia identifies signature landscapes.

New massively parallel sequencing technology enables, through deep sequencing of rearranged T-cell receptor (TCR) Vß complementarity-determining region 3 (CDR3) regions, a previously inaccessible level of TCR repertoire analysis. The CDR3 repertoire diversity reflects clonal composition, the potential antigenic recognition spectrum, and the quantity of available T-cell responses. In this context, T-large granular lymphocyte (T-LGL) leukemia is a chronic clonal lymphoproliferation of cytotoxic T cells often associated with autoimmune diseases and various cytopenias. Using CD8(+) T-LGL leukemia as a model disease, we set out to evaluate and compare the TCR deep-sequencing spectra of both patients and healthy controls to better understand how TCR deep sequencing could be used in the diagnosis and monitoring of not only T-LGL leukemia but also reactive processes such as autoimmune disease and infection. Our data demonstrate, with high resolution, significantly decreased diversity of the T-cell repertoire in CD8(+) T-LGL leukemia and suggest that many T-LGL clonotypes may be private to the disease and may not be present in the general public, even at the basal level.

STAT3 mutations indicate the presence of subclinical T-cell clones in a subset of aplastic anemia and myelodysplastic syndrome patients.

Large granular lymphocyte leukemia (LGL) is often associated with immune cytopenias and can cooccur in the context of aplastic anemia (AA) and myelodysplastic syndromes (MDS). We took advantage of the recent description of signal transducer and activator of transcription 3 (STAT3) mutations in LGL clonal expansions to test, using sensitive methods, for the presence of these mutations in a large cohort of 367 MDS and 140 AA cases. STAT3 clones can be found not only in known LGL concomitant cases, but in a small proportion of unsuspected ones (7% AA and 2.5% MDS). In STAT3-mutated AA patients, an interesting trend toward better responses of immunosuppressive therapy and an association with the presence of human leukocyte antigen-DR15 were found. MDSs harboring a STAT3 mutant clone showed a lower degree of bone marrow cellularity and a higher frequency of developing chromosome 7 abnormalities. STAT3-mutant LGL clones may facilitate a persistently dysregulated autoimmune activation, responsible for the primary induction of bone marrow failure in a subset of AA and MDS patients.

BRCC3 mutations in myeloid neoplasms.

Next generation sequencing technologies have provided insights into the molecular heterogeneity of various myeloid neoplasms, revealing previously unknown somatic genetic events. In our cohort of 1444 cases analyzed by next generation sequencing, somatic mutations in the gene BRCA1-BRCA2-containing complex 3 (BRCC3) were identified in 28 cases (1.9%). BRCC3 is a member of the JAMM/MPN+ family of zinc metalloproteases capable of cleaving Lys-63 linked polyubiquitin chains, and is implicated in DNA repair. The mutations were located throughout its coding region. The average variant allelic frequency of BRCC3 mutations was 30.1%, and by a serial sample analysis at two different time points a BRCC3 mutation was already identified in the initial stage of a myelodysplastic syndrome. BRCC3 mutations commonly occurred in nonsense (n=12), frameshift (n=4), and splice site (n=5) configurations. Due to the marginal male dominance (odds ratio; 2.00, 0.84-4.73) of BRCC3 mutations, the majority of mutations (n=23; 82%) were hemizygous. Phenotypically, BRCC3 mutations were frequently observed in myelodysplastic syndromes and myelodysplastic/myeloproliferative neoplasms and associated with -Y abnormality (odds ratio; 3.70, 1.25-11.0). Clinically, BRCC3 mutations were also related to higher age (P=0.01), although prognosis was not affected. Knockdown of Brcc3 gene expression in murine bone marrow lineage negative, Sca1 positive, c-kit positive cells resulted in 2-fold more colony formation and modest differentiation defect. Thus, BRCC3 likely plays a role as tumor-associated gene in myelodysplastic syndromes and myelodysplastic/myeloproliferative neoplasms.

Whole-exome sequencing enhances prognostic classification of myeloid malignancies.

PURPOSE: To date the standard nosology and prognostic schemes for myeloid neoplasms have been based on morphologic and cytogenetic criteria. We sought to test the hypothesis that a comprehensive, unbiased analysis of somatic mutations may allow for an improved classification of these diseases to predict outcome (overall survival). EXPERIMENTAL DESIGN: We performed whole-exome sequencing (WES) of 274 myeloid neoplasms, including myelodysplastic syndrome (MDS, N=75), myelodysplastic/myeloproliferative neoplasia (MDS/MPN, N=33), and acute myeloid leukemia (AML, N=22), augmenting the resulting mutational data with public WES results from AML (N=144). We fit random survival forests (RSFs) to the patient survival and clinical/cytogenetic data, with and without gene mutation information, to build prognostic classifiers. A targeted sequencing assay was used to sequence predictor genes in an independent cohort of 507 patients, whose accompanying data were used to evaluate performance of the risk classifiers. RESULTS: We show that gene mutations modify the impact of standard clinical variables on patient outcome, and therefore their incorporation hones the accuracy of prediction. The mutation-based classification scheme robustly predicted patient outcome in the validation set (log rank P=6.77 × 10(-21); poor prognosis vs. good prognosis categories HR 10.4, 95% CI 3.21-33.6). The RSF-based approach also compares favorably with recently-published efforts to incorporate mutational information for MDS prognosis. CONCLUSION: The results presented here support the inclusion of mutational information in prognostic classification of myeloid malignancies. Our classification scheme is implemented in a publicly available web-based tool (http://myeloid-risk. CASE: edu/).

Mutations in the spliceosome machinery, a novel and ubiquitous pathway in leukemogenesis.

Myelodysplastic syndromes (MDSs) are chronic and often progressive myeloid neoplasms associated with remarkable heterogeneity in the histomorphology and clinical course. Various somatic mutations are involved in the pathogenesis of MDS. Recently, mutations in a gene encoding a spliceosomal protein, SF3B1, were discovered in a distinct form of MDS with ring sideroblasts. Whole exome sequencing of 15 patients with myeloid neoplasms was performed, and somatic mutations in spliceosomal genes were identified. Sanger sequencing of 310 patients was performed to assess phenotype/genotype associations. To determine the functional effect of spliceosomal mutations, we evaluated pre-mRNA splicing profiles by RNA deep sequencing. We identified additional somatic mutations in spliceosomal genes, including SF3B1, U2AF1, and SRSF2. These mutations alter pre-mRNA splicing patterns. SF3B1 mutations are prevalent in low-risk MDS with ring sideroblasts, whereas U2AF1 and SRSF2 mutations are frequent in chronic myelomonocytic leukemia and advanced forms of MDS. SF3B1 mutations are associated with a favorable prognosis, whereas U2AF1 and SRSF2 mutations are predictive for shorter survival. Mutations affecting spliceosomal genes that result in defective splicing are a new leukemogenic pathway. Spliceosomal genes are probably tumor suppressors, and their mutations may constitute diagnostic biomarkers that could potentially serve as therapeutic targets.