Myelodysplasia after clonal hematopoiesis with APOBEC3-mediated CYBB inactivation in retroviral gene therapy for X-CGD.

Stem cell gene therapy using the MFGS-gp91phox retroviral vector was performed on a 27-year-old patient with X-linked chronic granulomatous disease (X-CGD) in 2014. The patient's refractory infections were resolved, whereas the oxidase-positive neutrophils disappeared within 6 months. Thirty-two months after gene therapy, the patient developed myelodysplastic syndrome (MDS), and vector integration into the MECOM locus was identified in blast cells. The vector integration into MECOM was detectable in most myeloid cells at 12 months after gene therapy. However, the patient exhibited normal hematopoiesis until the onset of MDS, suggesting that MECOM transactivation contributed to clonal hematopoiesis, and the blast transformation likely arose after the acquisition of additional genetic lesions. In whole-genome sequencing, the biallelic loss of the WT1 tumor suppressor gene, which occurred immediately before tumorigenesis, was identified as a potential candidate genetic alteration. The provirus CYBB cDNA in the blasts contained 108 G-to-A mutations exclusively in the coding strand, suggesting the occurrence of APOBEC3-mediated hypermutations during the transduction of CD34-positive cells. A hypermutation-mediated loss of oxidase activity may have facilitated the survival and proliferation of the clone with MECOM transactivation. Our data provide valuable insights into the complex mechanisms underlying the development of leukemia in X-CGD gene therapy.

A novel Menin-MLL1 inhibitor, DS-1594a, prevents the progression of acute leukemia with rearranged MLL1 or mutated NPM1.

BACKGROUND: Mixed lineage leukemia 1-rearranged (MLL1-r) acute leukemia patients respond poorly to currently available treatments and there is a need to develop more effective therapies directly disrupting the Menin‒MLL1 complex. Small-molecule-mediated inhibition of the protein‒protein interaction between Menin and MLL1 fusion proteins is a potential therapeutic strategy for patients with MLL1-r or mutated-nucleophosmin 1 (NPM1c) acute leukemia. In this study, we preclinically evaluated the new compound DS-1594a and its salts. METHODS: We evaluated the preclinical efficacy of DS-1594a as well as DS-1594a·HCl (the HCl salt of DS-1594a) and DS-1594a·succinate (the succinic acid salt of DS-1594a, DS-1594b) in vitro and in vivo using acute myeloid leukemia (AML)/acute lymphoblastic leukemia (ALL) models. RESULTS: Our results showed that MLL1-r or NPM1c human leukemic cell lines were selectively and highly sensitive to DS-1594a·HCl, with 50% growth inhibition values < 30 nM. Compared with cytrabine, the standard chemotherapy drug as AML therapy, both DS-1594a·HCl and DS-1594a·succinate mediated the eradication of potential leukemia-initiating cells by enhancing differentiation and reducing serial colony-forming potential in MLL1-r AML cells in vitro. The results were confirmed by flow cytometry, RNA sequencing, RT‒qPCR and chromatin immunoprecipitation sequencing analyses. DS-1594a·HCl and DS-1594a·succinate exhibited significant antitumor efficacy and survival benefit in MOLM-13 cell and patient-derived xenograft models of MLL1-r or NPM1c acute leukemia in vivo. CONCLUSION: We have generated a novel, potent, orally available small-molecule inhibitor of the Menin-MLL1 interaction, DS-1594a. Our results suggest that DS-1594a has medicinal properties distinct from those of cytarabine and that DS-1594a has the potential to be a new anticancer therapy and support oral dosing regimen for clinical studies (NCT04752163).

Quantitative assessment of copy number alterations by liquid biopsy for neuroblastoma.

Liquid biopsy, a method of detecting genomic alterations using blood specimens, has recently attracted attention as a noninvasive alternative to surgical tissue biopsy. We attempted quantitative analysis to detect amplification of MYCN (MYCNamp) and loss of heterozygosity at 11q (11qLOH), which are clinical requisites as prognostic factors of neuroblastoma (NB). In this study, cell-free DNA (cfDNA) was extracted from plasma samples from 24 NB patients at diagnosis. Copy numbers of MYCN and NAGK genes were quantitatively analyzed by droplet digital PCR (ddPCR). 11qLOH was also assessed by detecting allelic imbalances of heterozygous single nucleotide polymorphisms in the 11q region. The results obtained were compared to those of specimens from tumor tissues. The correlation coefficient of MYCN copy number of cfDNA and tumor DNA was 0.88 (p < 0.00001). 11qLOH was also accurately detected from cfDNA, except for one case with localized NB. Given the high accuracy of liquid biopsy, to investigate components of cfDNA, the proportion of tumor-derived DNA was estimated by examining the variant allele frequency of tumor-specific mutations in cfDNA. The proportion of tumor-derived DNA in cfDNA was 42.5% (range, 16.9%-55.9%), suggesting sufficient sensitivity of liquid biopsy for NB. In conclusion, MYCN copy number and 11qLOH could be quantitatively analyzed in plasma cfDNA by ddPCR assay. These results suggest that plasma cfDNA can be substituted for tumor DNA and can also be applied for comprehensive genomic profiling analysis.

De Novo Mutations Activating Germline TP53 in an Inherited Bone-Marrow-Failure Syndrome.

Inherited bone-marrow-failure syndromes (IBMFSs) include heterogeneous genetic disorders characterized by bone-marrow failure, congenital anomalies, and an increased risk of malignancy. Many lines of evidence have suggested that p53 activation might be central to the pathogenesis of IBMFSs, including Diamond-Blackfan anemia (DBA) and dyskeratosis congenita (DC). However, the exact role of p53 activation in each clinical feature remains unknown. Here, we report unique de novo TP53 germline variants found in two individuals with an IBMFS accompanied by hypogammaglobulinemia, growth retardation, and microcephaly mimicking DBA and DC. TP53 is a tumor-suppressor gene most frequently mutated in human cancers, and occasional germline variants occur in Li-Fraumeni cancer-predisposition syndrome. Most of these mutations affect the core DNA-binding domain, leading to compromised transcriptional activities. In contrast, the variants found in the two individuals studied here caused the same truncation of the protein, resulting in the loss of 32 residues from the C-terminal domain (CTD). Unexpectedly, the p53 mutant had augmented transcriptional activities, an observation not previously described in humans. When we expressed this mutant in zebrafish and human-induced pluripotent stem cells, we observed impaired erythrocyte production. These findings together with close similarities to published knock-in mouse models of TP53 lacking the CTD demonstrate that the CTD-truncation mutations of TP53 cause IBMFS, providing important insights into the previously postulated connection between p53 and IBMFSs.

Copy number alteration analysis for neuroblastoma using droplet digital polymerase chain reaction.

BACKGROUND: Neuroblastoma (NB) is a malignant tumor derived from the neural crest. MYCN amplification is a well-known adverse molecular prognostic factor for NB. Genome copy number alterations (CNAs) such as chromosome (Chr) 11q deletion, 1p deletion, and 17q gain are associated with a poor prognosis. Fluorescence in situ hybridization (FISH) and Southern blotting analysis are frequently used to detect MYCN amplification. Although comparative genomic hybridization (CGH) and single-nucleotide polymorphism (SNP) chip arrays can easily detect CNAs, these methods are impractical for clinical use due to their cost and run time. Consequently, genome copy number analysis using digital droplet PCR has become widely used to monitor CNAs. METHODS: In this study, we used digital droplet polymerase chain reaction to detect MYCN amplification and Chr 11q CNA, which was used for risk stratification according to the International Neuroblastoma Risk Group classification system. We compared the results with data from SNP chip arrays in seven NB cell lines and eight primary NB samples. RESULTS: Digital droplet PCR assays successfully detected MYCN amplification and 11q CNA. The results were very consistent with those obtained by SNP chip assay. CONCLUSIONS: Digital droplet PCR can be conducted more rapidly than FISH or Southern blotting. Accordingly, it should be useful for on-site clinical applications aimed at detecting CNAs in NB and performing risk stratification promptly after diagnosis.

Evolutionary basis of HLA-DPB1 alleles affects acute GVHD in unrelated donor stem cell transplantation.

HLA-DPB1 T-cell epitope (TCE) mismatching algorithm and rs9277534 SNP at the 3' untranslated region (3'UTR) in the HLA-DPB1 gene are key factors for transplant-related events in unrelated hematopoietic cell transplantation (UR-HCT). However, the association of these 2 mechanisms has not been elucidated. We analyzed 19 frequent HLA-DPB1 alleles derived from Japanese healthy subjects by next-generation sequencing of the entire HLA-DPB1 gene region and multi-SNP data of the HLA region in 1589 UR-HCT pairs. The risk of acute graft-versus-host disease (aGVHD) was analyzed in 1286 patients with single HLA-DPB1 mismatch UR-HCT. The phylogenetic tree constructed using the entire gene region demonstrated that HLA-DPB1 alleles were divided into 2 groups, HLA-DP2 and HLA-DP5. Although a phylogenetic relationship in the genomic region from exon 3 to 3'UTR (Ex3-3'UTR) obviously supported the division of HLA-DP2 and HLA-DP5 groups, which in exon 2 showed intermingling of HLA-DPB1 alleles in a non-HLA-DP2 and non-HLA-DP5-group manner. Multi-SNP data also showed 2 discriminative HLA-DPB1 groups according to Ex3-3'UTR. Risk of grade 2-4 aGVHD was significantly higher in patient HLA-DP5 group mismatch than patient HLA-DP2 group mismatch (hazard ratio, 1.28; P = .005), regardless of donor mismatch HLA-DP group. Regarding TCE mismatch, increasing risk of aGVHD in patient HLA-DP5 group mismatch and TCE-nonpermissive mismatch were observed only in patients with TCE-permissive mismatch and patient HLA-DP2 group mismatch, respectively. Evolutionary analysis revealed that rs9277534 represented a highly conserved HLA-DPB1 Ex3-3'UTR region and may provoke aGVHD differently to TCE mismatching algorithm, reflecting exon 2 polymorphisms. These findings enrich our understanding of the mechanism of aGVHD in HLA-DPB1 mismatch UR-HCT.

Prognostic relevance of integrated genetic profiling in adult T-cell leukemia/lymphoma.

Adult T-cell leukemia/lymphoma (ATL) is a heterogeneous group of peripheral T-cell malignancies characterized by human T-cell leukemia virus type-1 infection, whose genetic profile has recently been fully investigated. However, it is still poorly understood how these alterations affect clinical features and prognosis. We investigated the effects of genetic alterations commonly found in ATL on disease phenotypes and clinical outcomes, based on genotyping data obtained from 414 and 463 ATL patients using targeted-capture sequencing and single nucleotide polymorphism array karyotyping, respectively. Aggressive (acute/lymphoma) subtypes were associated with an increased burden of genetic and epigenetic alterations, higher frequencies of TP53 and IRF4 mutations, and many copy number alterations (CNAs), including PD-L1 amplifications and CDKN2A deletions, compared with indolent (chronic/smoldering) subtypes. By contrast, STAT3 mutations were more characteristic of indolent ATL. Higher numbers of somatic mutations and CNAs significantly correlated with worse survival. In a multivariate analysis incorporating both clinical factors and genetic alterations, the Japan Clinical Oncology Group prognostic index high-risk, older age, PRKCB mutations, and PD-L1 amplifications were independent poor prognostic factors in aggressive ATL. In indolent ATL, IRF4 mutations, PD-L1 amplifications, and CDKN2A deletions were significantly associated with shorter survival, although the chronic subtype with unfavorable clinical factors was only marginally significant. Thus, somatic alterations characterizing aggressive diseases predict worse prognosis in indolent ATL, among which PD-L1 amplifications are a strong genetic predictor in both aggressive and indolent ATL. ATL subtypes are further classified into molecularly distinct subsets with different prognosis. Genetic profiling might contribute to improved prognostication and management of ATL patients.

Gene expression and risk of leukemic transformation in myelodysplasia.

Myelodysplastic syndromes (MDSs) are a heterogeneous group of clonal hematopoietic disorders with a highly variable prognosis. To identify a gene expression-based classification of myelodysplasia with biological and clinical relevance, we performed a comprehensive transcriptomic analysis of myeloid neoplasms with dysplasia using transcriptome sequencing. Unsupervised clustering of gene expression data of bone marrow CD34+ cells from 100 patients identified 2 subgroups. The first subtype was characterized by increased expression of genes related to erythroid/megakaryocytic (EMK) lineages, whereas the second subtype showed upregulation of genes related to immature progenitor (IMP) cells. Compared with the first so-called EMK subtype, the IMP subtype showed upregulation of many signaling pathways and downregulation of several pathways related to metabolism and DNA repair. The IMP subgroup was associated with a significantly shorter survival in both univariate (hazard ratio [HR], 5.0; 95% confidence interval [CI], 1.8-14; P = .002) and multivariate analysis (HR, 4.9; 95% CI, 1.3-19; P = .02). Leukemic transformation was limited to the IMP subgroup. The prognostic significance of our classification was validated in an independent cohort of 183 patients. We also constructed a model to predict the subgroups using gene expression profiles of unfractionated bone marrow mononuclear cells (BMMNCs). The model successfully predicted clinical outcomes in a test set of 114 patients with BMMNC samples. The addition of our classification to the clinical model improved prediction of patient outcomes. These results indicated biological and clinical relevance of our gene expression-based classification, which will improve risk prediction and treatment stratification of MDS.

Somatic Mutations and Clonal Hematopoiesis in Aplastic Anemia.

BACKGROUND: In patients with acquired aplastic anemia, destruction of hematopoietic cells by the immune system leads to pancytopenia. Patients have a response to immunosuppressive therapy, but myelodysplastic syndromes and acute myeloid leukemia develop in about 15% of the patients, usually many months to years after the diagnosis of aplastic anemia. METHODS: We performed next-generation sequencing and array-based karyotyping using 668 blood samples obtained from 439 patients with aplastic anemia. We analyzed serial samples obtained from 82 patients. RESULTS: Somatic mutations in myeloid cancer candidate genes were present in one third of the patients, in a limited number of genes and at low initial variant allele frequency. Clonal hematopoiesis was detected in 47% of the patients, most frequently as acquired mutations. The prevalence of the mutations increased with age, and mutations had an age-related signature. DNMT3A-mutated and ASXL1-mutated clones tended to increase in size over time; the size of BCOR- and BCORL1-mutated and PIGA-mutated clones decreased or remained stable. Mutations in PIGA and BCOR and BCORL1 correlated with a better response to immunosuppressive therapy and longer and a higher rate of overall and progression-free survival; mutations in a subgroup of genes that included DNMT3A and ASXL1 were associated with worse outcomes. However, clonal dynamics were highly variable and might not necessarily have predicted the response to therapy and long-term survival among individual patients. CONCLUSIONS: Clonal hematopoiesis was prevalent in aplastic anemia. Some mutations were related to clinical outcomes. A highly biased set of mutations is evidence of Darwinian selection in the failed bone marrow environment. The pattern of somatic clones in individual patients over time was variable and frequently unpredictable. (Funded by Grant-in-Aid for Scientific Research and others.).

Variegated RHOA mutations in adult T-cell leukemia/lymphoma.

Adult T-cell leukemia/lymphoma (ATLL) is a distinct form of peripheral T-cell lymphoma with poor prognosis, which is caused by the human T-lymphotropic virus type 1 (HTLV-1). In contrast to the unequivocal importance of HTLV-1 infection in the pathogenesis of ATLL, the role of acquired mutations in HTLV-1 infected T cells has not been fully elucidated, with a handful of genes known to be recurrently mutated. In this study, we identified unique RHOA mutations in ATLL through whole genome sequencing of an index case, followed by deep sequencing of 203 ATLL samples. RHOA mutations showed distinct distribution and function from those found in other cancers. Involving 15% (30/203) of ATLL cases, RHOA mutations were widely distributed across the entire coding sequence but almost invariably located at the guanosine triphosphate (GTP)-binding pocket, with Cys16Arg being most frequently observed. Unexpectedly, depending on mutation types and positions, these RHOA mutants showed different or even opposite functional consequences in terms of GTP/guanosine diphosphate (GDP)-binding kinetics, regulation of actin fibers, and transcriptional activation. The Gly17Val mutant did not bind GTP/GDP and act as a dominant negative molecule, whereas other mutants (Cys16Arg and Ala161Pro) showed fast GTP/GDP cycling with enhanced transcriptional activation. These findings suggest that both loss- and gain-of-RHOA functions could be involved in ATLL leukemogenesis. In summary, our study not only provides a novel insight into the molecular pathogenesis of ATLL but also highlights a unique role of variegation of heterologous RHOA mutations in human cancers.

ASXL2 regulates hematopoiesis in mice and its deficiency promotes myeloid expansion.

Chromosomal translocation t(8;21)(q22;q22) which leads to the generation of oncogenic RUNX1-RUNX1T1 (AML1-ETO) fusion is observed in approximately 10% of acute myelogenous leukemia (AML). To identify somatic mutations that co-operate with t(8;21)-driven leukemia, we performed whole and targeted exome sequencing of an Asian cohort at diagnosis and relapse. We identified high frequency of truncating alterations in ASXL2 along with recurrent mutations of KIT, TET2, MGA, FLT3, and DHX15 in this subtype of AML. To investigate in depth the role of ASXL2 in normal hematopoiesis, we utilized a mouse model of ASXL2 deficiency. Loss of ASXL2 caused progressive hematopoietic defects characterized by myeloid hyperplasia, splenomegaly, extramedullary hematopoiesis, and poor reconstitution ability in transplantation models. Parallel analyses of young and >1-year old Asxl2-deficient mice revealed age-dependent perturbations affecting, not only myeloid and erythroid differentiation, but also maturation of lymphoid cells. Overall, these findings establish a critical role for ASXL2 in maintaining steady state hematopoiesis, and provide insights into how its loss primes the expansion of myeloid cells.

BRCA1 alterations with additional defects in DNA damage response genes may confer chemoresistance to BRCA-like breast cancers treated with neoadjuvant chemotherapy.

The BRCA-like phenotype is a feature that some sporadic breast cancers share with those occurring in BRCA1 or BRCA2 mutation carriers. As tumors with the phenotype have defects in the DNA damage response pathway, which may increase sensitivity to drugs such as DNA cross-linking agents and PARP inhibitors, a method to identify this phenotype is important. The prediction of chemoresistance, which frequently develops in these tumors, is also crucial for improving therapy. We examined genomic aberrations and BRCA1 promoter methylation in tumors of 73 breast cancer (20 HR-/HER2- and 53 HR+/HER2-) patients, who received neoadjuvant chemotherapy with anthracycline, cyclophosphamide, and taxane, using SNP array CGH and quantitative PCR. The methylation and/or loss or uniparental disomy (UPD) of BRCA1 (BRCA1 alterations) and the loss or UPD of BRCA2 (BRCA2 alterations) were detected in 27 (37%) and 21 (29%), respectively, of the 73 tumors. Tumors with BRCA1 or BRCA2 alterations were associated with a higher number of genomic aberrations (P < 0.001 and P < 0.001) and higher percentage of TP53 alterations (P < 0.001 and P < 0.001) than those without. Overall survival (OS) rates were similar between patients with or without BRCA1 or BRCA2 alterations. However, when 27 patients with BRCA1-altered tumors were classified into those with or without the loss or UPD of PALB2, PAGR1, RAD51B, FANCM, MLL4, or ERCC1/2 in tumors, patients with additional defects in DNA damage response genes had worse OS (P = 0.037, 0.045, 0.038, 0.044, 0.041, or 0.019) than those without. These defects may confer chemoresistance and predict poor outcomes in patients with BRCA1-altered breast cancer.

Genome-wide surveillance of mismatched alleles for graft-versus-host disease in stem cell transplantation.

Acute graft-versus-host disease (aGVHD) represents one of the major complications in allogeneic stem cell transplantation and is primarily caused by genetic disparity between the donor and recipient. In HLA-matched transplants, the disparity is thought to be determined by loci encoding minor histocompatibility antigens (minor H antigens), which are presented by specific HLA molecules. We performed a genome-wide association study (GWAS) to identify minor H antigen loci associated with aGVHD. A total of 500 568 single nucleotide polymorphisms (SNPs) were genotyped for donors and recipients from 1589 unrelated bone marrow transplants matched for HLA-A, -B, -C, -DRB1, and -DQB1, followed by the imputation of unobserved SNPs. We interrogated SNPs whose disparity between the donor and recipient was significantly associated with aGVHD development. Without assuming HLA unrestriction, we successfully captured a known association between HLA-DPB1 disparity (P = 4.50 × 10(-9)) and grade II-IV aGVHD development, providing proof of concept for the GWAS design aimed at discovering genetic disparity associated with aGVHD. In HLA-restricted analyses, whereby association tests were confined to major subgroups sharing common HLA alleles to identify putative minor H antigen loci, we identified 3 novel loci significantly associated with grade III-IV aGVHD. Among these, rs17473423 (P = 1.20 × 10(-11)) at 12p12.1 within the KRAS locus showed the most significant association in the subgroup, sharing HLA-DQB1*06:01. Our result suggested that a GWAS can be successfully applied to identify allele mismatch associated with aGVHD development, contributing to the understanding of the genetic basis of aGVHD.

Frequent pathway mutations of splicing machinery in myelodysplasia.

Myelodysplastic syndromes and related disorders (myelodysplasia) are a heterogeneous group of myeloid neoplasms showing deregulated blood cell production with evidence of myeloid dysplasia and a predisposition to acute myeloid leukaemia, whose pathogenesis is only incompletely understood. Here we report whole-exome sequencing of 29 myelodysplasia specimens, which unexpectedly revealed novel pathway mutations involving multiple components of the RNA splicing machinery, including U2AF35, ZRSR2, SRSF2 and SF3B1. In a large series analysis, these splicing pathway mutations were frequent (∼45 to ∼85%) in, and highly specific to, myeloid neoplasms showing features of myelodysplasia. Conspicuously, most of the mutations, which occurred in a mutually exclusive manner, affected genes involved in the 3'-splice site recognition during pre-mRNA processing, inducing abnormal RNA splicing and compromised haematopoiesis. Our results provide the first evidence indicating that genetic alterations of the major splicing components could be involved in human pathogenesis, also implicating a novel therapeutic possibility for myelodysplasia.

Two susceptibility loci to Takayasu arteritis reveal a synergistic role of the IL12B and HLA-B regions in a Japanese population.

Takayasu arteritis (TAK) is an autoimmune systemic vasculitis of unknown etiology. Although previous studies have revealed that HLA-B*52:01 has an effect on TAK susceptibility, no other genetic determinants have been established so far. Here, we performed genome scanning of 167 TAK cases and 663 healthy controls via Illumina Infinium Human Exome BeadChip arrays, followed by a replication study consisting of 212 TAK cases and 1,322 controls. As a result, we found that the IL12B region on chromosome 5 (rs6871626, overall p = 1.7 × 10(-13), OR = 1.75, 95% CI 1.42-2.16) and the MLX region on chromosome 17 (rs665268, overall p = 5.2 × 10(-7), OR = 1.50, 95% CI 1.28-1.76) as well as the HLA-B region (rs9263739, a proxy of HLA-B*52:01, overall p = 2.8 × 10(-21), OR = 2.44, 95% CI 2.03-2.93) exhibited significant associations. A significant synergistic effect of rs6871626 and rs9263739 was found with a relative excess risk of 3.45, attributable proportion of 0.58, and synergy index of 3.24 (p ≤ 0.00028) in addition to a suggestive synergistic effect between rs665268 and rs926379 (p ≤ 0.027). We also found that rs6871626 showed a significant association with clinical manifestations of TAK, including increased risk and severity of aortic regurgitation, a representative severe complication of TAK. Detection of these susceptibility loci will provide new insights to the basic mechanisms of TAK pathogenesis. Our findings indicate that IL12B plays a fundamental role on the pathophysiology of TAK in combination with HLA-B(∗)52:01 and that common autoimmune mechanisms underlie the pathology of TAK and other autoimmune disorders such as psoriasis and inflammatory bowel diseases in which IL12B is involved as a genetic predisposing factor.

High-risk HLA alleles for severe acute graft-versus-host disease and mortality in unrelated donor bone marrow transplantation.

HLA molecules play an important role for immunoreactivity in allogeneic hematopoietic stem cell transplantation. To elucidate the effect of specific HLA alleles on acute graft-versus-host disease, we conducted a retrospective analysis using 6967 Japanese patients transplanted with T-cell-replete marrow from an unrelated donor. Using unbiased searches of patient and donor HLA alleles, patient and/or donor HLA-B*51:01 (patient: HR, 1.37,P<0.001; donor: HR, 1.35,P<0.001) and patient HLA-C*14:02 (HR, 1.35,P<0.001) were significantly associated with an increased risk of severe acute graft-versus-host disease. The finding that donor HLA-C*14:02 was not associated with severe acute graft-versus-host disease prompted us to elucidate the relation of these high-risk HLA alleles with patient and donor HLA-C allele mismatches. In comparison to HLA-C allele match, patient mismatched HLA-C*14:02 showed the highest risk of severe acute graft-versus-host disease (HR, 3.61,P<0.001) and transplant-related mortality (HR, 2.53,P<0.001) among all patient mismatched HLA-C alleles. Although patient HLA-C*14:02 and donor HLA-C*15:02 mismatch was usually KIR2DL-ligand mismatch in the graft-versus-host direction, the risk of patient mismatched HLA-C*14:02 for severe acute graft-versus-host disease was obvious regardless of KIR2DL-ligand matching. The effect of patient and/or donor HLA-B*51:01 on acute graft-versus-host disease was attributed not only to strong linkage disequilibrium of HLA-C*14:02 and -B*51:01, but also to the effect of HLA-B*51:01 itself. With regard to clinical implications, patient mismatched HLA-C*14:02 proved to be a potent risk factor for severe acute graft-versus-host disease and mortality, and should be considered a non-permissive HLA-C mismatch in donor selection for unrelated donor hematopoietic stem cell transplantation.