Leveraging functional genomic annotations and genome coverage to improve polygenic prediction of complex traits within and between ancestries.

We develop a method, SBayesRC, that integrates genome-wide association study (GWAS) summary statistics with functional genomic annotations to improve polygenic prediction of complex traits. Our method is scalable to whole-genome variant analysis and refines signals from functional annotations by allowing them to affect both causal variant probability and causal effect distribution. We analyze 50 complex traits and diseases using ∼7 million common single-nucleotide polymorphisms (SNPs) and 96 annotations. SBayesRC improves prediction accuracy by 14% in European ancestry and up to 34% in cross-ancestry prediction compared to the baseline method SBayesR, which does not use annotations, and outperforms other methods, including LDpred2, LDpred-funct, MegaPRS, PolyPred-S and PRS-CSx. Investigation of factors affecting prediction accuracy identifies a significant interaction between SNP density and annotation information, suggesting whole-genome sequence variants with annotations may further improve prediction. Functional partitioning analysis highlights a major contribution of evolutionary constrained regions to prediction accuracy and the largest per-SNP contribution from nonsynonymous SNPs.

Genome-wide analysis in over 1 million individuals of European ancestry yields improved polygenic risk scores for blood pressure traits.

Hypertension affects more than one billion people worldwide. Here we identify 113 novel loci, reporting a total of 2,103 independent genetic signals (P < 5 × 10-8) from the largest single-stage blood pressure (BP) genome-wide association study to date (n = 1,028,980 European individuals). These associations explain more than 60% of single nucleotide polymorphism-based BP heritability. Comparing top versus bottom deciles of polygenic risk scores (PRSs) reveals clinically meaningful differences in BP (16.9 mmHg systolic BP, 95% CI, 15.5-18.2 mmHg, P = 2.22 × 10-126) and more than a sevenfold higher odds of hypertension risk (odds ratio, 7.33; 95% CI, 5.54-9.70; P = 4.13 × 10-44) in an independent dataset. Adding PRS into hypertension-prediction models increased the area under the receiver operating characteristic curve (AUROC) from 0.791 (95% CI, 0.781-0.801) to 0.826 (95% CI, 0.817-0.836, ∆AUROC, 0.035, P = 1.98 × 10-34). We compare the 2,103 loci results in non-European ancestries and show significant PRS associations in a large African-American sample. Secondary analyses implicate 500 genes previously unreported for BP. Our study highlights the role of increasingly large genomic studies for precision health research.

Genetic control of DNA methylation is largely shared across European and East Asian populations.

DNA methylation is an ideal trait to study the extent of the shared genetic control across ancestries, effectively providing hundreds of thousands of model molecular traits with large QTL effect sizes. We investigate cis DNAm QTLs in three European (n = 3701) and two East Asian (n = 2099) cohorts to quantify the similarities and differences in the genetic architecture across populations. We observe 80,394 associated mQTLs (62.2% of DNAm probes with significant mQTL) to be significant in both ancestries, while 28,925 mQTLs (22.4%) are identified in only a single ancestry. mQTL effect sizes are highly conserved across populations, with differences in mQTL discovery likely due to differences in allele frequency of associated variants and differing linkage disequilibrium between causal variants and assayed SNPs. This study highlights the overall similarity of genetic control across ancestries and the value of ancestral diversity in increasing the power to detect associations and enhancing fine mapping resolution.

MXenes as emerging adsorbents for removal of environmental pollutants.

MXenes are a recently emerging class of two-dimensional nanomaterials that have gained considerable interest in the field of environmental protection. Owing to their high surface area, abundant terminal groups, and unique two-dimensional layered structures, MXenes have demonstrated high efficacy as adsorbents for various pollutants. Here we focused on the latest developments in the field of MXene-based adsorbents, including the structure and properties of MXenes, their synthesis and modification methods, and their adsorption performance and mechanisms for various pollutants. Among the pollutants that have been reported to be adsorbed by MXenes are radionuclides (U(VI), Sr(II), Cs(I), Eu(III), Ba(II), Th(IV), and Tc(VII)/Re(VII)), heavy metals (Hg(II), Cu(II), Cr(VI), and Pb(II)), dyes, per- and polyfluoroalkyl substances (PFAS), antibiotics (tetracycline, ciprofloxacin, and sulfonamides), antibiotic resistance genes (ARGs), and other contaminates. Moreover, future directions in MXene research are also suggested in this review.

Rapid anti-myeloma activity by T cells expressing an anti-BCMA CAR with a human heavy-chain-only antigen-binding domain.

Multiple myeloma (MM) is a rarely curable malignancy of plasma cells. MM expresses B cell maturation antigen (BCMA). We developed a fully human anti-BCMA chimeric antigen receptor (CAR) with a heavy-chain-only antigen-recognition domain, a 4-1BB domain, and a CD3ζ domain. The CAR was designated FHVH33-CD8BBZ. We conducted the first-in-humans clinical trial of T cells expressing FHVH33-CD8BBZ (FHVH-T). Twenty-five patients with relapsed MM were treated. The stringent complete response rate (sCR) was 52%. Median progression-free survival (PFS) was 78 weeks. Of 24 evaluable patients, 6 (25%) had a maximum cytokine-release syndrome (CRS) grade of 3; no patients had CRS of greater than grade 3. Most anti-MM activity occurred within 2-4 weeks of FHVH-T infusion as shown by decreases in the rapidly changing MM markers serum free light chains, urine light chains, and bone marrow plasma cells. Blood CAR+ cell levels peaked during the time that MM elimination was occurring, between 7 and 15 days after FHVH-T infusion. C-C chemokine receptor type 7 (CCR7) expression on infusion CD4+ FHVH-T correlated with peak blood FHVH-T levels. Single-cell RNA sequencing revealed a shift toward more differentiated FHVH-T after infusion. Anti-CAR antibody responses were detected in 4 of 12 patients assessed. FHVH-T has powerful, rapid, and durable anti-MM activity.

MiR-20 promotes cartilage repair in knee osteoarthritis rats via BMP2/Smad1 signaling pathway.

To explore the effect of micro ribonucleic acid (miR)-20b on knee osteoarthritis rats by regulating the bone morphogenetic protein 2 (BMP2)/Smad1 pathway, a total of 36 SD rats were randomly divided into normal group (n=12), model group (n=12) and miR-20b mimics group (n=12). The rats in normal group were fed normally, while those in model group and miR-20b mimics group were used to establish knee osteoarthritis models. After modeling, model group was not given any intervention, but miR-20b mimics group received intra-articular injection of miR-20b mimics once a day for 2 weeks. Basso, Beattie and Bresnahan (BBB) limb motor function scoring was performed at 1, 5, 7 and 14 days after the modeling, and samples were obtained after 2 weeks of intervention. Next, hematoxylin and eosin (H&E) staining was applied to observe tissue morphology, Markin's scoring was utilized to evaluate articular cartilage degeneration, and immunohistochemistry was employed to detect the expressions of BMP2 and Smad1. Thereafter, the expression of miR-20b was detected via qPCR, the content of cartilage oligomeric matrix protein (COMP) and C-telopeptide of type II collagen (CTX-II) was measured via enzyme-linked immunosorbent assay (ELISA), and the expressions of BMP2 and Smad1 proteins were examined via Western blotting (WB). BBB limb motor function scoring showed that compared with that in normal group, the BBB limb motor function score of rats in the other two groups was reduced (P<0.05). In comparison with that in model group, the BBB limb motor function score in miR-20b mimics group was increased from the 7th day after intervention (P<0.05). In addition, H&E staining results manifested that the articular surface in normal group was smooth and flat, with normal morphology, clear structure and no obvious damage. In model group, the articular surface was not smooth and uneven, and more articular cartilage fractures, morphological disorders and structural damages could be observed. Moreover, the articular surface in miR-20b mimics group was slightly damaged and smoother, and its morphology and structure were markedly improved in contrast to that in model group. The Markin's score in normal group was lower than that in model group and miR-20b mimics group (P<0.05), and it was overtly decreased in miR-20b mimics group in comparison with that in model group (P<0.05). Next, immunohistochemistry demonstrated that compared with normal group, the other two groups had lowered positive expressions of BMP2 and Smad1 (P<0.05). In comparison with model group, miR-20b mimics group exhibited notably raised positive expressions of BMP2 and Smad1 (P<0.05). Then it was found from qPCR results that the expression level of miR-20b in the other two groups was overtly reduced compared with that in normal group (P<0.05), and it was prominently elevated in miR-20b mimics group in contrast to that in model group (P<0.05). Besides, ELISA illustrated that the content of COMP and CTX-II in the cartilage tissues in the other two groups was evidently reduced compared with that in normal group (P<0.05), and it was increased prominently in miR-20b mimics group compared with that in model group (P<0.05). Finally, it was revealed by WB examination that the relative expression levels of BMP2 and Smad1 proteins in the other two groups markedly declined in comparison with those in normal group (P<0.05), and they were elevated in contrast to those in model group (P<0.05). MiR-20b can promote cartilage repair and improve articular function in knee osteoarthritis rats by up-regulating the BMP2/Smad1 signaling pathway.

Impact of humid climate on rheumatoid arthritis faecal microbiome and metabolites.

Studies have shown that high humidity is a condition that aggravates the pain of rheumatoid arthritis (RA), but the relevant mechanism is controversial. Currently, there is a lack of experimental animal studies on high humidity as an adverse factor related to the pathogenesis of RA. We used healthy SD rats and collagen-induced arthritis (CIA) rats to investigate the effects of high humidity on arthritis. Integrated metabolomics analyses of faeces and 16S rRNA sequencing of the faecal microbiota were performed to comprehensively assess the diversity of the faecal microbiota and metabolites in healthy and CIA rats. In this study, high humidity aggravated arthritis in CIA rats, which manifested as articular cartilage lesions, increased arthritis scores, and an increase in proinflammatory cytokines. High humidity had a certain effect on the articular cartilage extent, arthritis score and proinflammatory cytokines of healthy rats as well. Furthermore, high humidity caused significant changes in faecal microbes and metabolites in both healthy and CIA rats. 16S rRNA sequencing of faecal samples showed that high humidity increased the amount of inflammation-related bacteria in healthy and CIA rats. Faecal metabolomics results showed that high humidity significantly altered the level of faecal metabolites in healthy rats and CIA rats, and the changes in biological functions were mainly related to the inflammatory response and oxidative stress. Combined analysis showed that there was a strong correlation between the faecal microbiota and faecal metabolites. High humidity is an adverse factor for the onset and development of RA, and its mechanism is related to the inflammatory response and oxidative stress. However, the question of how high humidity impacts RA pathogenesis needs to be further investigated.

Learning functional conservation between human and pig to decipher evolutionary mechanisms underlying gene expression and complex traits.

Assessment of genomic conservation between humans and pigs at the functional level can improve the potential of pigs as a human biomedical model. To address this, we developed a deep learning-based approach to learn the genomic conservation at the functional level (DeepGCF) between species by integrating 386 and 374 functional profiles from humans and pigs, respectively. DeepGCF demonstrated better prediction performance compared with the previous method. In addition, the resulting DeepGCF score captures the functional conservation between humans and pigs by examining chromatin states, sequence ontologies, and regulatory variants. We identified a core set of genomic regions as functionally conserved that plays key roles in gene regulation and is enriched for the heritability of complex traits and diseases in humans. Our results highlight the importance of cross-species functional comparison in illustrating the genetic and evolutionary basis of complex phenotypes.

Leveraging base-pair mammalian constraint to understand genetic variation and human disease.

Thousands of genomic regions have been associated with heritable human diseases, but attempts to elucidate biological mechanisms are impeded by an inability to discern which genomic positions are functionally important. Evolutionary constraint is a powerful predictor of function, agnostic to cell type or disease mechanism. Single-base phyloP scores from 240 mammals identified 3.3% of the human genome as significantly constrained and likely functional. We compared phyloP scores to genome annotation, association studies, copy-number variation, clinical genetics findings, and cancer data. Constrained positions are enriched for variants that explain common disease heritability more than other functional annotations. Our results improve variant annotation but also highlight that the regulatory landscape of the human genome still needs to be further explored and linked to disease.

Leveraging Base Pair Mammalian Constraint to Understand Genetic Variation and Human Disease.

Although thousands of genomic regions have been associated with heritable human diseases, attempts to elucidate biological mechanisms are impeded by a general inability to discern which genomic positions are functionally important. Evolutionary constraint is a powerful predictor of function that is agnostic to cell type or disease mechanism. Here, single base phyloP scores from the whole genome alignment of 240 placental mammals identified 3.5% of the human genome as significantly constrained, and likely functional. We compared these scores to large-scale genome annotation, genome-wide association studies (GWAS), copy number variation, clinical genetics findings, and cancer data sets. Evolutionarily constrained positions are enriched for variants explaining common disease heritability (more than any other functional annotation). Our results improve variant annotation but also highlight that the regulatory landscape of the human genome still needs to be further explored and linked to disease.

mBAT-combo: A more powerful test to detect gene-trait associations from GWAS data.

Gene-based association tests aggregate multiple SNP-trait associations into sets defined by gene boundaries and are widely used in post-GWAS analysis. A common approach for gene-based tests is to combine SNPs associations by computing the sum of χ2 statistics. However, this strategy ignores the directions of SNP effects, which could result in a loss of power for SNPs with masking effects, e.g., when the product of two SNP effects and the linkage disequilibrium (LD) correlation is negative. Here, we introduce "mBAT-combo," a set-based test that is better powered than other methods to detect multi-SNP associations in the context of masking effects. We validate the method through simulations and applications to real data. We find that of 35 blood and urine biomarker traits in the UK Biobank, 34 traits show evidence for masking effects in a total of 4,273 gene-trait pairs, indicating that masking effects is common in complex traits. We further validate the improved power of our method in height, body mass index, and schizophrenia with different GWAS sample sizes and show that on average 95.7% of the genes detected only by mBAT-combo with smaller sample sizes can be identified by the single-SNP approach with a 1.7-fold increase in sample sizes. Eleven genes significant only in mBAT-combo for schizophrenia are confirmed by functionally informed fine-mapping or Mendelian randomization integrating gene expression data. The framework of mBAT-combo can be applied to any set of SNPs to refine trait-association signals hidden in genomic regions with complex LD structures.

Mechanism of interventional effect and targets of Zhuyu pill in regulating and suppressing colitis and cholestasis.

Zhuyu pill (ZYP) is a traditional Chinese medicine prescription composed of two drugs, Coptis chinensis Franch. and Tetradium ruticarpum (A. Jussieu) T. G. Hartley, and is commonly used in the clinical treatment of diseases of the digestive system. However, the mechanism underlying the effect of ZYP on colitis remains unclear. In this study, a colitis rat model was induced with 2,4,6-trinitro-benzenesulfonic acid (TNBS, 100 mg/kg) and treated with ZYP (low dose: 0.6 g/kg, high dose: 1.2 g/kg). Disease activity index, colonic weight index, and weight change ratio were used to evaluate the model and efficacy. LC-MS and 16S rRNA gene sequencing were used to measure differences in fecal metabolism and microorganism population among the control, model, low-dose ZYP, and high-dose ZYP groups. To elucidate the mechanism of interventional effect of ZYP, Spearman correlation analysis was used to analyze the correlation between fecal metabolism and fecal microbial number. High-dose and low-dose ZYP both exhibited significant interventional effects on colitis rat models, and high-dose ZYP produced a better interventional effect compared with low-dose ZYP. Based on a metabolomics test of fecal samples, significantly altered metabolites in the model and high-dose ZYP treatment groups were identified. In total, 492 metabolites were differentially expressed. Additionally, sequencing of the 16S rRNA gene in fecal samples revealed that the high-dose ZYP could improve TNBS-induced fecal microbiota dysbiosis. Ultimately, changes in tryptophan metabolism and Firmicutes and Gammaproteobacteria populations were detected after ZYP treatment in both colitis and cholestasis. Therefore, we conclude that tryptophan metabolism and Firmicutes and Gammaproteobacteria populations are the core targets of the anti-inflammatory effect of ZYP. These findings provide a scientific basis for further investigation of the anti-inflammatory mechanism of ZYP in the future.

Adoptive Cellular Therapy with Autologous Tumor-Infiltrating Lymphocytes and T-cell Receptor-Engineered T Cells Targeting Common p53 Neoantigens in Human Solid Tumors.

Adoptive cellular therapy (ACT) targeting neoantigens can achieve durable clinical responses in patients with cancer. Most neoantigens arise from patient-specific mutations, requiring highly individualized treatments. To broaden the applicability of ACT targeting neoantigens, we focused on TP53 mutations commonly shared across different cancer types. We performed whole-exome sequencing on 163 patients with metastatic solid cancers, identified 78 who had TP53 missense mutations, and through immunologic screening, identified 21 unique T-cell reactivities. Here, we report a library of 39 T-cell receptors (TCR) targeting TP53 mutations shared among 7.3% of patients with solid tumors. These TCRs recognized tumor cells in a TP53 mutation- and human leucocyte antigen (HLA)-specific manner in vitro and in vivo. Twelve patients with chemorefractory epithelial cancers were treated with ex vivo-expanded autologous tumor-infiltrating lymphocytes (TIL) that were naturally reactive against TP53 mutations. However, limited clinical responses (2 partial responses among 12 patients) were seen. These infusions contained low frequencies of mutant p53-reactive TILs that had exhausted phenotypes and showed poor persistence. We also treated one patient who had chemorefractory breast cancer with ACT comprising autologous peripheral blood lymphocytes transduced with an allogeneic HLA-A*02-restricted TCR specific for p53R175H. The infused cells exhibited an improved immunophenotype and prolonged persistence compared with TIL ACT and the patient experienced an objective tumor regression (-55%) that lasted 6 months. Collectively, these proof-of-concept data suggest that the library of TCRs targeting shared p53 neoantigens should be further evaluated for the treatment of patients with advanced human cancers. See related Spotlight by Klebanoff, p. 919.

Assessing the contribution of rare variants to complex trait heritability from whole-genome sequence data.

Analyses of data from genome-wide association studies on unrelated individuals have shown that, for human traits and diseases, approximately one-third to two-thirds of heritability is captured by common SNPs. However, it is not known whether the remaining heritability is due to the imperfect tagging of causal variants by common SNPs, in particular whether the causal variants are rare, or whether it is overestimated due to bias in inference from pedigree data. Here we estimated heritability for height and body mass index (BMI) from whole-genome sequence data on 25,465 unrelated individuals of European ancestry. The estimated heritability was 0.68 (standard error 0.10) for height and 0.30 (standard error 0.10) for body mass index. Low minor allele frequency variants in low linkage disequilibrium (LD) with neighboring variants were enriched for heritability, to a greater extent for protein-altering variants, consistent with negative selection. Our results imply that rare variants, in particular those in regions of low linkage disequilibrium, are a major source of the still missing heritability of complex traits and disease.

Molecular signatures of antitumor neoantigen-reactive T cells from metastatic human cancers.

The accurate identification of antitumor T cell receptors (TCRs) represents a major challenge for the engineering of cell-based cancer immunotherapies. By mapping 55 neoantigen-specific TCR clonotypes (NeoTCRs) from 10 metastatic human tumors to their single-cell transcriptomes, we identified signatures of CD8+ and CD4+ neoantigen-reactive tumor-infiltrating lymphocytes (TILs). Neoantigen-specific TILs exhibited tumor-specific expansion with dysfunctional phenotypes, distinct from blood-emigrant bystanders and regulatory TILs. Prospective prediction and testing of 73 NeoTCR signature-derived clonotypes demonstrated that half of the tested TCRs recognized tumor antigens or autologous tumors. NeoTCR signatures identified TCRs that target driver neoantigens and nonmutated viral or tumor-associated antigens, suggesting a common metastatic TIL exhaustion program. NeoTCR signatures delineate the landscape of TILs across metastatic tumors, enabling successful TCR prediction based purely on TIL transcriptomic states for use in cancer immunotherapy.

Improved analyses of GWAS summary statistics by reducing data heterogeneity and errors.

statistics from genome-wide association studies (GWAS) have facilitated the development of various summary data-based methods, which typically require a reference sample for linkage disequilibrium (LD) estimation. Analyses using these methods may be biased by errors in GWAS summary data or LD reference or heterogeneity between GWAS and LD reference. Here we propose a quality control method, DENTIST, that leverages LD among genetic variants to detect and eliminate errors in GWAS or LD reference and heterogeneity between the two. Through simulations, we demonstrate that DENTIST substantially reduces false-positive rate in detecting secondary signals in the summary-data-based conditional and joint association analysis, especially for imputed rare variants (false-positive rate reduced from >28% to <2% in the presence of heterogeneity between GWAS and LD reference). We further show that DENTIST can improve other summary-data-based analyses such as fine-mapping analysis.

A generalized linear mixed model association tool for biobank-scale data.

Compared with linear mixed model-based genome-wide association (GWA) methods, generalized linear mixed model (GLMM)-based methods have better statistical properties when applied to binary traits but are computationally much slower. In the present study, leveraging efficient sparse matrix-based algorithms, we developed a GLMM-based GWA tool, fastGWA-GLMM, that is severalfold to orders of magnitude faster than the state-of-the-art tools when applied to the UK Biobank (UKB) data and scalable to cohorts with millions of individuals. We show by simulation that the fastGWA-GLMM test statistics of both common and rare variants are well calibrated under the null, even for traits with extreme case-control ratios. We applied fastGWA-GLMM to the UKB data of 456,348 individuals, 11,842,647 variants and 2,989 binary traits (full summary statistics available at http://fastgwa.info/ukbimpbin ), and identified 259 rare variants associated with 75 traits, demonstrating the use of imputed genotype data in a large cohort to discover rare variants for binary complex traits.

iTRAQ-based quantitative proteomics analysis of the hepatoprotective effect of melatonin on ANIT-induced cholestasis in rats.

The therapeutic effects of melatonin on cholestatic liver injury have received widespread attention recently. The aim of the present study was to investigate the mechanisms of the anti-cholestatic effects of melatonin against α-naphthyl isothiocyanate (ANIT)-induced liver injury in rats and to screen for potential biomarkers of cholestasis through isobaric tags for relative and absolute quantitation (iTRAQ) proteomics. Rats orally received melatonin (100 mg/kg body weight) or an equivalent volume of 0.25% carboxymethyl cellulose sodium salt 12 h after intraperitoneal injection of ANIT (75 mg/kg) and were subsequently sacrificed at 36 h after injection. Liver biochemical indices were determined and liver tissue samples were stained using hematoxylin-eosin staining, followed by iTRAQ quantitative proteomics to identify potential underlying therapeutic mechanisms and biomarkers. The results suggested that the expression levels of alanine transaminase, aspartate aminotransferase, total bilirubin and direct bilirubin were reduced in the rats treated with melatonin. Histopathological observation indicated that melatonin was effective in the treatment of ANIT-induced cholestasis. iTRAQ proteomics results suggested that melatonin-mediated reduction in ANIT-induced cholestasis may be associated with enhanced antioxidant function and relieving abnormal fatty acid metabolism. According to pathway enrichment analysis using the Kyoto Encyclopedia of Genes and Genomes, the major metabolic pathways for the metabolism of melatonin are fatty acid degradation, the peroxisome proliferator-activated receptor signaling pathway, fatty acid metabolism, chemical carcinogenesis, carbon metabolism, pyruvate metabolism, fatty acid biosynthesis and retinol metabolism, as well as drug metabolism via cytochrome P450. Malate dehydrogenase 1 and glutathione S-transferase Yb-3 may serve as potential targets in the treatment of ANIT-induced cholestasis with melatonin.

Identification of neoantigen-reactive T lymphocytes in the peripheral blood of a patient with glioblastoma.

The adoptive transfer of naturally occurring T cells that recognize cancer neoantigens has led to durable tumor regressions in select patients with cancer. However, it remains unknown whether such T cells can be isolated from and used to treat patients with glioblastoma, a cancer that is refractory to currently available therapies. To answer this question, we stimulated patient blood-derived memory T cells in vitro using peptides and minigenes that represented point mutations unique to patients' tumors (ie, candidate neoantigens) and then tested their ability to specifically recognize these mutations. In a cohort of five patients with glioblastoma, we found that circulating CD4+ memory T cells from one patient recognized a cancer neoantigen harboring a mutation in the EED gene (EEDH189N) that was unique to that patient's tumor. This finding suggests that neoantigen-reactive T cells could indeed be isolated from patients with glioblastoma, thereby providing a rationale for further efforts to develop neoantigen-directed adoptive T cell therapy for this disease.

Direct identification of neoantigen-specific TCRs from tumor specimens by high-throughput single-cell sequencing.

BACKGROUND: Recognition of neoantigens by T cells plays a major role in cancer immunotherapy. Identification of neoantigen-specific T-cell receptors (TCRs) has become a critical research tool for studying T cell-mediated responses after immunotherapy. In addition, neoantigen-specific TCRs can be used to modify the specificity of T cells for T cell-based therapies targeting tumor-specific mutations. Although several techniques have been developed to identify TCR sequences, these techniques still require a significant amount of labor, making them impractical in the clinical setting. METHODS: Thanks to the availability of high-throughput single-cell sequencing, we developed a new process to isolate neoantigen-specific TCR sequences. This process included the isolation of tumor-infiltrating T cells from a tumor specimen and the stimulation of T cells by neoantigen-loaded dendritic cells, followed by single-cell sequencing for TCR and T-cell activation markers, interferon-γ and interleukin-2. RESULTS: In this study, potential neoantigen-specific TCRs were isolated from three melanoma and three colorectal tumor specimens. These TCRs were then synthesized and transduced into autologous T cells, followed by testing the recognition of neoantigens. A total of 28 neoantigen-specific TCRs were identified by this process. If identical TCR sequences were detected from two or more single cells, this approach was highly reliable (100%, 19 out of 19 TCRs). CONCLUSION: This single-cell approach provides an efficient process to isolate antigen-specific TCRs for research and clinical applications.