Rigorous benchmarking of T-cell receptor repertoire profiling methods for cancer RNA sequencing.

The ability to identify and track T-cell receptor (TCR) sequences from patient samples is becoming central to the field of cancer research and immunotherapy. Tracking genetically engineered T cells expressing TCRs that target specific tumor antigens is important to determine the persistence of these cells and quantify tumor responses. The available high-throughput method to profile TCR repertoires is generally referred to as TCR sequencing (TCR-Seq). However, the available TCR-Seq data are limited compared with RNA sequencing (RNA-Seq). In this paper, we have benchmarked the ability of RNA-Seq-based methods to profile TCR repertoires by examining 19 bulk RNA-Seq samples across 4 cancer cohorts including both T-cell-rich and T-cell-poor tissue types. We have performed a comprehensive evaluation of the existing RNA-Seq-based repertoire profiling methods using targeted TCR-Seq as the gold standard. We also highlighted scenarios under which the RNA-Seq approach is suitable and can provide comparable accuracy to the TCR-Seq approach. Our results show that RNA-Seq-based methods are able to effectively capture the clonotypes and estimate the diversity of TCR repertoires, as well as provide relative frequencies of clonotypes in T-cell-rich tissues and low-diversity repertoires. However, RNA-Seq-based TCR profiling methods have limited power in T-cell-poor tissues, especially in highly diverse repertoires of T-cell-poor tissues. The results of our benchmarking provide an additional appealing argument to incorporate RNA-Seq into the immune repertoire screening of cancer patients as it offers broader knowledge into the transcriptomic changes that exceed the limited information provided by TCR-Seq.

Deletion of CD36 exhibits limited impact on normal hematopoiesis and the leukemia microenvironment.

BACKGROUND: CD36 has been identified as a potential therapeutic target both in leukemic cells and in the tumor immune microenvironment. In acute myeloid leukemia (AML), we found that APOC2 acts with CD36 to promote leukemia growth by activating the LYN-ERK signaling. CD36 also plays a role in lipid metabolism of cancer associated T-cells leading to impaired cytotoxic CD8+ T-cell and enhanced Treg cell function. To establish CD36 as a viable therapeutic target in AML, we investigated whether targeting CD36 has any detrimental impact on normal hematopoietic cells. METHODS: Differential expression data of CD36 during human and mouse normal hematopoiesis were examined and compared. Cd36 knockout (Cd36-KO) mice were evaluated for blood analysis, hematopoietic stem cells and progenitors (HSPCs) function and phenotype analyses, and T cells in vitro expansion and phenotypes in comparison with wild type (WT) mice. In addition, MLL-PTD/FLT3-ITD leukemic cells were engrafted into Cd36-KO and WT mice, and leukemia burden was compared between groups. RESULTS: RNA-Seq data showed that Cd36 expression was low in HSPCs and increased as cells matured. Phenotypic analysis revealed limited changes in blood count except for a slight yet significantly lower red blood cell count and hemoglobin and hematocrit levels in Cd36-KO mice compared with WT mice (P < 0.05). In vitro cell proliferation assays of splenocytes and HSPCs from Cd36-KO mice showed a similar pattern of expansion to that of cells from WT mice. Characterization of HSPCs showed similar percentages of the different progenitor cell populations between Cd36-KO with WT mice. However, Cd36-KO mice exhibited ~ 40% reduction of the number of colonies developed from HSPCs cells compared with WT mice (P < 0.001). Cd36-KO and WT mice presented comparably healthy BM transplant in non-competitive models and developed similar leukemia burden. CONCLUSIONS: Although the loss of Cd36 affects the hematopoietic stem cell and erythropoiesis, limited detrimental overall impact was observed on normal Hematopoietic and leukemic microenvironments. Altogether, considering the limited impact on normal hematopoiesis, therapeutic approaches to target CD36 in cancer are unlikely to result in toxicity to normal blood cells.

Profiling of T Cell Repertoire in SARS-CoV-2-Infected COVID-19 Patients Between Mild Disease and Pneumonia.

The coronavirus disease 2019 (COVID-19) pandemic, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has become a public health emergency. The most common symptoms of COVID-19 are fever, cough, and fatigue. While most patients with COVID-19 present with mild illness, some patients develop pneumonia, an important risk factor for mortality, at early stage of viral infection, putting these patients at increased risk of death. So far, little has been known about differences in the T cell repertoires between COVID-19 patients with and without pneumonia during SARS-CoV-2 infection. Herein, we aimed to investigate T cell receptor (TCR) repertoire profiles and patient-specific SARS-CoV-2-associated TCR clusters between COVID-19 patients with mild disease (no sign of pneumonia) and pneumonia. The TCR sequencing was conducted to characterize the peripheral TCR repertoire profile and diversity. The TCR clustering and CDR3 annotation were exploited to further discover groups of patient-specific TCR clonotypes with potential SARS-CoV-2 antigen specificities. Our study indicated a slight decrease in the TCR repertoire diversity and a skewed CDR3 length usage in patients with pneumonia compared to those with mild disease. The SARS-CoV-2-associated TCR clusters enriched in patients with mild disease exhibited significantly higher TCR generation probabilities and most of which were highly shared among patients, compared with those from pneumonia patients. Importantly, using similarity network-based clustering followed by the sequence conservation analysis, we found different patterns of CDR3 sequence motifs between mild disease- and pneumonia-specific SARS-CoV-2-associated public TCR clusters. Our results showed that characteristics of overall TCR repertoire and SARS-CoV-2-associated TCR clusters/clonotypes were divergent between COVID-19 patients with mild disease and patients with pneumonia. These findings provide important insights into the correlation between the TCR repertoire and disease severity in COVID-19 patients.

Human Granulocyte-Macrophage Colony-Stimulating Factor Fused to Elastin-Like Polypeptides Assembles Biologically-Active Nanoparticles.

Human granulocyte-macrophage colony-stimulating factor (hGMCSF) is crucial in the immune system as it stimulates survival, proliferation, differentiation, and functional activation of myeloid hematopoietic cells. hGMCSF is integral to approved therapies, including monoclonal antibodies against checkpoint inhibitors, chimeric antigen receptors, and prevention of chemotherapy-induced neutropenia. Recombinant hGMCSF can be purified from Escherichia. coli; however, it forms inclusion bodies that require solubilization and refolding. Alternatively, this manuscript describes its fusion with an elastin-like polypeptide (ELP). Previously reported as purification tags and solubility enhancers, ELPs are recombinant polypeptides that undergo reversible temperature-dependent phase separation. This report is the first to show that fusion to an ELP enables direct purification of hGMCSF fusions from the soluble fraction of bacterial lysate. Surprisingly, these ELP-fusions assemble stable, small, spherical nanoparticles that maintain pro-mitotic activity of hGMCSF. These nanoparticles exhibit ELP-mediated phase separation; however, nanoparticle assembly significantly increases the entropic and enthalpic cost of phase separation compared to ELP alone. The attachment of a high molecular weight ELP to a difficult-to-express protein, like hGMCSF, appears to be a useful strategy to stabilize bioactive, protein-based nanoparticles, which may have broad applications in medicine and biology.

Clinical and preclinical characterization of CD99 isoforms in acute myeloid leukemia.

In an effort to identify target genes in acute myeloid leukemia (AML), we compared gene expression profiles between normal and AML cells from various publicly available datasets. We identified CD99, a gene that is up-regulated in AML patients. In 186 patients from The Cancer Genome Atlas AML dataset, CD99 was over-expressed in patients with FLT3-ITD and was down-regulated in patients with TP53 mutations. CD99 is a trans-membrane protein expressed on leukocytes and plays a role in cell adhesion, trans-endothelial migration, and T-cell differentiation. The CD99 gene encodes two isoforms with distinct expression and functional profiles in both normal and malignant tissues. Here we report that, although the CD99 long isoform initially induces an increase in cell proliferation, it also induces higher levels of reactive oxygen species, DNA damage, apoptosis and a subsequent decrease in cell viability. In several leukemia murine models, the CD99 long isoform delayed disease progression and resulted in lower leukemia engraftment in the bone marrow. Furthermore, the CD99 monoclonal antibody reduced cell viability, colony formation, and cell migration, and induced cell differentiation and apoptosis in leukemia cell lines and primary blasts. Mechanistically, CD99 long isoform resulted in transient induction followed by a dramatic decrease in both ERK and SRC phosphorylation. Altogether, our study provides new insights into the role of CD99 isoforms in AML that could potentially be relevant for the preclinical development of CD99 targeted therapy.

Anti-CD99 scFv-ELP nanoworms for the treatment of acute myeloid leukemia.

CD99 is a transmembrane glycoprotein shown to be upregulated in various malignancies. We have previously reported CD99 to be highly upregulated and present a viable therapeutic target in acute myeloid leukemia (AML). Currently, no therapy against CD99 is under clinical investigation. As a surface molecule, CD99 can be targeted with an antibody-based approach. Here, we have developed a new modality to target CD99 by engineering a fusion protein composed of a single-chain variable fragment antibody (anti-CD99 scFv) conjugated with a high molecular weight elastin-like polypeptide (ELP), A192: α-CD99-A192. This fusion protein assembles into multi-valent nanoworm with optimal physicochemical properties and favorable pharmacokinetic parameters (half-life: 16 h). α-CD99-A192 nanoworms demonstrated excellent in vitro and in vivo anti-leukemic effects. α-CD99-A192 induced apoptotic cell death in AML cell lines and primary blasts and prolonged overall survival of AML xenograft mouse model.

Upregulation of the EMT marker vimentin is associated with poor clinical outcome in acute myeloid leukemia.

BACKGROUND: Vimentin (VIM) is a type III intermediate filament that maintains cell integrity, and is involved in cell migration, motility and adhesion. When overexpressed in solid cancers, vimentin drives epithelial to mesenchymal transition (EMT) and ultimately, metastasis. The effects of its overexpression in AML are unclear. METHODS: In this study, we analyzed the TCGA data of 173 AML patients for which complete clinical and expression data were available. In this analysis, we assessed the association between VIM mRNA expression and patient's clinical and molecular characteristics including clinical outcome. RESULTS: VIM overexpression was associated with higher white blood count (< p = 0.0001). Patients with high VIM expression have worse overall survival (OS) and disease-free survival (DFS) compared with patients with low VIM expression (median OS; 7.95 months vs 19.2 months; p = 0.029). After age-stratification, high VIM expression was significantly associated with worse overall survival in older patients (age ≥ 60; median OS: 5.4 vs 9.9 months: p = 0.0257) but not in younger patients (age < 60). In stratification analysis according to cytogenetic status, high VIM expression was significantly associated with shorter OS (7.95 vs 24.6 months: p = 0.0102) in cytogenetically normal, but not in cytogenetic abnormal AML. CONCLUSIONS: Collectively, the data indicate that overexpression of the EMT marker vimentin is associated with poor clinical outcome in older patients with cytogenetically normal AML; and therefore may play a role in this disease.

The era of immunogenomics/immunopharmacogenomics.

Although germline alterations and somatic mutations in disease cells have been extensively analyzed, molecular changes in immune cells associated with disease conditions have not been characterized in depth. It is clear that our immune system has a critical role in various biological and pathological conditions, such as infectious diseases, autoimmune diseases, drug-induced skin and liver toxicity, food allergy, and rejection of transplanted organs. The recent development of cancer immunotherapies, particularly drugs modulating the immune checkpoint molecules, has clearly demonstrated the importance of host immune cells in cancer treatments. However, the molecular mechanisms by which these new therapies kill tumor cells are still not fully understood. In this regard, we have begun to explore the role of newly developed tools such as next-generation sequencing in the genetic characterization of both cancer cells and host immune cells, a field that is called immunogenomics/ immunopharmacogenomics. This new field has enormous potential to help us better understand changes in our immune system during the course of various disease conditions. Here we report the potential of deep sequencing of T-cell and B-cell receptors in capturing the molecular contribution of the immune system, which we believe plays critical roles in the pathogenesis of various human diseases.

Quantitative characterization of T-cell repertoire and biomarkers in kidney transplant rejection.

BACKGROUND: T-cell-mediated rejection (TCMR) remains a major cause of kidney allograft failure. The characterization of T-cell repertoire in different immunological disorders has emerged recently as a novel tool with significant implications. We herein sought to characterize T-cell repertoire using next generation sequencing to diagnose TCMR. METHODS: In this prospective study, we analyzed samples from 50 kidney transplant recipients. We collected blood and kidney transplant biopsy samples at sequential time points before and post transplant. We used next generation sequencing to characterize T-cell receptor (TCR) repertoire by using illumina miSeq on cDNA synthesized from RNA extracted from six patients' samples. We also measured RNA expression levels of FOXP3, CD8, CD4, granzyme and perforin in blood samples from all 50 patients. RESULTS: Seven patients developed TCMR during the first three months of the study. Out of six patients who had complete sets of blood and biopsy samples two had TCMR. We found an expansion of the TCR repertoire in blood at time of rejection when compared to that at pre-transplant or one-month post transplant. Patients with TCMR (n = 7) had significantly higher RNA expression levels of FOXP3, Perforin, Granzyme, CD4 and CD8 in blood samples than those with no TCMR (n = 43) (P = 0.02, P = 0.003, P = 0.002, P = 0.017, and P = 0.01, respectively). CONCLUSIONS: Our study provides a potential utilization of TCR clone kinetics analysis in the diagnosis of TCMR. This approach may allow for the identification of the expanded T-cell clones associated with the rejection and lead to potential noninvasive diagnosis and targeted therapies of TCMR.

Lenalidomide-mediated enhanced translation of C/EBPα-p30 protein up-regulates expression of the antileukemic microRNA-181a in acute myeloid leukemia.

Recently, we showed that increased miR-181a expression was associated with improved outcomes in cytogenetically normal acute myeloid leukemia (CN-AML). Interestingly, miR-181a expression was increased in CN-AML patients harboring CEBPA mutations, which are usually biallelic and associate with better prognosis. CEBPA encodes the C/EBPα transcription factor. We demonstrate here that the presence of N-terminal CEBPA mutations and miR-181a expression are linked. Indeed, the truncated C/EBPα-p30 isoform, which is produced from the N-terminal mutant CEBPA gene or from the differential translation of wild-type CEBPA mRNA and is commonly believed to have no transactivation activity, binds to the miR-181a-1 promoter and up-regulates the microRNA expression. Furthermore, we show that lenalidomide, a drug approved for myelodysplastic syndromes and multiple myeloma, enhances translation of the C/EBPα-p30 isoform, resulting in higher miR-181a levels. In xenograft mouse models, ectopic miR-181a expression inhibits tumor growth. Similarly, lenalidomide exhibits antitumorigenic activity paralleled by increased miR-181a expression. This regulatory pathway may explain an increased sensitivity to apoptosis-inducing chemotherapy in subsets of AML patients. Altogether, our data provide a potential explanation for the improved clinical outcomes observed in CEBPA-mutated CN-AML patients, and suggest that lenalidomide treatment enhancing the C/EBPα-p30 protein levels and in turn miR-181a may sensitize AML blasts to chemotherapy.

Differences in the mutational landscape of clonal hematopoiesis of indeterminate potential among races and between male and female patients with cancer.

Clonal hematopoiesis of indeterminate potential (CHIP) has emerged as a significant precursor to hematological malignancies and is associated with several age-related diseases. We leveraged public data to explore differences in the mutational landscape of CHIP between males (Ms) and females (Fs) and across diverse racial populations. DNA (cytosine-5) methyltransferase 3 alpha (DNMT3A) mutations were substantially more prevalent in Fs than in Ms (38.94% vs. 31.37%, p-value: < 0.001, q-value: < 0.001). Additional sex combs-like 1 (ASXL1) mutations were more frequent in Ms than Fs (5.82% vs. 2.69%, p-value < 0.001, q-value < 0.001). In the racial cohorts with sufficient sample sizes, STAT5B and CSF1R mutations were most frequent in Asian populations (1.40% and 0.84%), followed by Black populations (0.98% and 0.24%) and White populations (0.29% and 0.09%) (p-value: < 0.001 , q-value: 0.023 for both genes). Several other CHIP mutations were enriched in Black: RARA, SMAD2, CDKN1B, CENPA, CTLA4, EIF1AX, ELF3, MSI1, MYC, SOX17, and AURKA. On the other hand, H3C1, H3C4, and MYCL were enriched in the Asian cohort. Our analysis highlights sex and racial differences in CHIP mutations among patients with cancer. As CHIP continues to gain recognition as a critical precursor to malignancies and other diseases, understanding how these differences contribute to CHIP's underlying mechanisms and clinical implications is critical.

The DNA methylation landscape across the TCR loci in patients with acute myeloid leukemia.

The capacity of T cells to initiate anti-leukemia immune responses is determined by the ability of their receptors (TCRs) to recognize leukemia neoantigens. Epigenetic mechanisms including DNA methylation contribute to shaping the TCR repertoire composition and diversity. The DNA hypomethylating agents (HMAs) have been widely used in the treatment of acute myeloid leukemia (AML) and myelodysplastic syndrome (MDS). Whether DNA HMAs directly influence TCR gene loci methylation patterns remains unknown. By analyzing public datasets, we compared methylation patterns across TCR loci in AML patients and healthy controls. We also explored how HMAs influence TCR loci DNA methylation in patients with AML. While methylation patterns are largely conserved across the TCR loci, certain V genes exhibit high interindividual variability. Although overall methylation levels within the TCR loci did not show significant differences, specific sites, including 32 TRAV and 12 TRBV sites exhibited distinct methylation patterns when comparing T cells from healthy donors to those from patients with AML. In leukemic cells, decitabine treatment demethylates sites across the TRAV and TRBV genes. While not as significant, a similar pattern of demethylation is observed in T cells. Pretreatment AML samples exhibit higher methylation beta values in differentially methylated positions (DMPs) compared with non-DMPs. Methylation levels of certain TRAV and TRBV genes in leukemic cells are associated with patients' risk status. The presence of disease specific TCR loci methylated signatures that are associated with clinical outcome presents an opportunity for therapeutic intervention. HMAs can modulate the TCR loci methylation patterns, yet whether they could reprogram the TCR repertoire composition remains to be explored.

FLT3/CD99 Bispecific Antibody-Based Nanoparticles for Acute Myeloid Leukemia.

Cluster of differentiation 99 (CD99) is a receptor that is significantly upregulated in acute myeloid leukemia (AML). FMS-like tyrosine kinase 3 internal tandem duplication mutation in AML (FLT3-ITD AML) exhibits even higher levels of CD99 expression. Our group previously employed a novel peptide platform technology called elastin-like polypeptides and fused it with single-chain antibodies capable of binding to FLT3 (FLT3-A192) or CD99 (CD99-A192). Targeting either FLT3 or CD99 using FLT3-A192 or CD99-A192 led to AML cell death and reduced leukemia burden in AML mouse models. Here, we report on the development of a novel Co-Assembled construct that is capable of binding to both CD99 and FLT3 and the antileukemia activity of the bispecific construct in FLT3-ITD AML preclinical models. This dual-targeting Co-Assembled formulation exhibits cytotoxic effects on AML cells (AML cell lines and primary blasts) and reduced leukemia burden and prolonged survival in FLT3-ITD AML mouse models. Altogether, this study demonstrates the potential of an innovative therapeutic strategy that targets both FLT3 and CD99 in FLT3-ITD AML. SIGNIFICANCE: This study investigates a dual-targeting strategy in acute myeloid leukemia (AML), focusing on FLT3 and CD99. The approach demonstrates enhanced therapeutic potential, presenting a novel option for AML treatment.

A rigorous benchmarking of alignment-based HLA typing algorithms for RNA-seq data.

Accurate identification of human leukocyte antigen (HLA) alleles is essential for various clinical and research applications, such as transplant matching and drug sensitivities. Recent advances in RNA-seq technology have made it possible to impute HLA types from sequencing data, spurring the development of a large number of computational HLA typing tools. However, the relative performance of these tools is unknown, limiting the ability for clinical and biomedical research to make informed choices regarding which tools to use. Here we report the study design of a comprehensive benchmarking of the performance of 12 HLA callers across 682 RNA-seq samples from 8 datasets with molecularly defined gold standard at 5 loci, HLA-A, -B, -C, -DRB1, and -DQB1. For each HLA typing tool, we will comprehensively assess their accuracy, compare default with optimized parameters, and examine for discrepancies in accuracy at the allele and loci levels. We will also evaluate the computational expense of each HLA caller measured in terms of CPU time and RAM. We also plan to evaluate the influence of read length over the HLA region on accuracy for each tool. Most notably, we will examine the performance of HLA callers across European and African groups, to determine discrepancies in accuracy associated with ancestry. We hypothesize that RNA-Seq HLA callers are capable of returning high-quality results, but the tools that offer a good balance between accuracy and computational expensiveness for all ancestry groups are yet to be developed. We believe that our study will provide clinicians and researchers with clear guidance to inform their selection of an appropriate HLA caller.

Genomics and transcriptomic alterations of the glutamate receptors in acute myeloid leukemia.

Glutamine and glutamate have been widely explored as potential therapeutic targets in acute myeloid leukemia (AML). In addition to its bioenergetic role in leukemia cell proliferation, L-glutamate is a neurotransmitter that acts on glutamate receptors. However, the role of glutamate receptors in AML is largely understudied. Here, we comprehensively analyze the genomic and transcriptomic alterations of glutamate receptor genes in AML using publicly available data. We investigated the frequency of mutations in the glutamate receptor genes and whether an association exist between the presence of these mutations and clinical and molecular characteristics or patient's clinical outcome. We also assessed the dysregulation of glutamate receptor gene expression in AML with and without mutations and whether gene dysregulation is associated with clinical outcomes. We found that 29 (14.5%) of 200 patients with AML had a mutation in at least one glutamate receptor gene. The DNMT3A mutations were significantly more frequent in patients with mutations in at least one glutamate receptor gene compared with patients without mutations (13 of 29 [44.8%] vs. 41 of 171 [23.9%], p value: 0.02). Notably, patients with mutations in at least one glutamate receptor gene survived shorter than patients without mutations; however, the results did not reach statistical significance (overall survival: 15.5 vs. 19.0 months; p value: 0.10). Mutations in the glutamate receptor genes were not associated with changes in gene expression and the transcriptomic levels of glutamate receptor genes were not associated with clinical outcome.

Features of the TCR repertoire associate with patients' clinical and molecular characteristics in acute myeloid leukemia.

Background: Allogeneic hematopoietic stem cell transplant remains the most effective strategy for patients with high-risk acute myeloid leukemia (AML). Leukemia-specific neoantigens presented by the major histocompatibility complexes (MHCs) are recognized by the T cell receptors (TCR) triggering the graft-versus-leukemia effect. A unique TCR signature is generated by a complex V(D)J rearrangement process to form TCR capable of binding to the peptide-MHC. The generated TCR repertoire undergoes dynamic changes with disease progression and treatment. Method: Here we applied two different computational tools (TRUST4 and MIXCR) to extract the TCR sequences from RNA-seq data from The Cancer Genome Atlas (TCGA) and examine the association between features of the TCR repertoire in adult patients with AML and their clinical and molecular characteristics. Results: We found that only ~30% of identified TCR CDR3s were shared by the two computational tools. Yet, patterns of TCR associations with patients' clinical and molecular characteristics based on data obtained from either tool were similar. The numbers of unique TCR clones were highly correlated with patients' white blood cell counts, bone marrow blast percentage, and peripheral blood blast percentage. Multivariable regressions of TCRA and TCRB median normalized number of unique clones with mutational status of AML patients using TRUST4 showed significant association of TCRA or TCRB with WT1 mutations, WBC count, %BM blast, and sex (adjusted in TCRB model). We observed a correlation between TCRA/B number of unique clones and the expression of T cells inhibitory signal genes (TIGIT, LAG3, CTLA-4) and foxp3, but not IL2RA, CD69 and TNFRSF9 suggestive of exhausted T cell phenotypes in AML. Conclusion: Benchmarking of computational tools is needed to increase the accuracy of the identified clones. The utilization of RNA-seq data enables identification of highly abundant TCRs and correlating these clones with patients' clinical and molecular characteristics. This study further supports the value of high-resolution TCR-Seq analyses to characterize the TCR repertoire in patients.