Search for: All records

Neoantigens are derived from tumor-specific somatic mutations. Neoantigen-based synthesized peptides have been under clinical investigation to boost cancer immunotherapy efficacy. The promising results prompt us to further elucidate the effect of neoantigen expression on patient survival in breast cancer. We applied Kaplan–Meier survival and multivariable Cox regression models to evaluate the effect of neoantigen expression and its interaction with T-cell activation on overall survival in a cohort of 729 breast cancer patients. Pearson’s chi-squared tests were used to assess the relationships between neoantigen expression and clinical pathological variables. Spearman correlation analysis was conducted to identify correlations between neoantigen expression, mutation load, and DNA repair gene expression. ERCC1, XPA, and XPC were negatively associated with neoantigen expression, while BLM, BRCA2, MSH2, XRCC2, RAD51, CHEK1, and CHEK2 were positively associated with neoantigen expression. Based on the multivariable Cox proportional hazard model, patients with a high level of neoantigen expression and activated T-cell status showed improved overall survival. Similarly, in the T-cell exhaustion and progesterone receptor (PR) positive subgroups, patients with a high level of neoantigen expression showed prolonged survival. In contrast, there was no significant difference in the T-cell activation and PR negative subgroups. In conclusion, neoantigens may serve as immunogenic agents for immunotherapy in breast cancer. 

Single-cell RNA sequencing technology provides an opportunity to study gene expression at single-cell resolution. However, prevalent dropout events result in high data sparsity and noise that may obscure downstream analyses in single-cell transcriptomic studies. We propose a new method, G2S3, that imputes dropouts by borrowing information from adjacent genes in a sparse gene graph learned from gene expression profiles across cells. We applied G2S3 and ten existing imputation methods to eight single-cell transcriptomic datasets and compared their performance. Our results demonstrated that G2S3 has superior overall performance in recovering gene expression, identifying cell subtypes, reconstructing cell trajectories, identifying differentially expressed genes, and recovering gene regulatory and correlation relationships. Moreover, G2S3 is computationally efficient for imputation in large-scale single-cell transcriptomic datasets. 

Abstract Transfer RNA-derived fragments (tRFs) are a new class of small non-coding RNAs and play important roles in biological and physiological processes. Prediction of tRF target genes and binding sites is crucial in understanding the biological functions of tRFs in the molecular mechanisms of human diseases. We developed a publicly accessible web-based database, tRFtarget (http://trftarget.net), for tRF target prediction. It contains the computationally predicted interactions between tRFs and mRNA transcripts using the two state-of-the-art prediction tools RNAhybrid and IntaRNA, including location of the binding sites on the target, the binding region, and free energy of the binding stability with graphic illustration. tRFtarget covers 936 tRFs and 135 thousand predicted targets in eight species. It allows researchers to search either target genes by tRF IDs or tRFs by gene symbols/transcript names. We also integrated the manually curated experimental evidence of the predicted interactions into the database. Furthermore, we provided a convenient link to the DAVID® web server to perform downstream functional pathway analysis and gene ontology annotation on the predicted target genes. This database provides useful information for the scientific community to experimentally validate tRF target genes and facilitate the investigation of the molecular functions and mechanisms of tRFs. 

Longitudinal phenotypes have been increasingly available in genome-wide association studies (GWAS) and electronic health record-based studies for identification of genetic variants that influence complex traits over time. For longitudinal binary data, there remain significant challenges in gene mapping, including misspecification of the model for phenotype distribution due to ascertainment. Here, we propose L-BRAT (Longitudinal Binary-trait Retrospective Association Test), a retrospective, generalized estimating equation-based method for genetic association analysis of longitudinal binary outcomes. We also develop RGMMAT, a retrospective, generalized linear mixed model-based association test. Both tests are retrospective score approaches in which genotypes are treated as random conditional on phenotype and covariates. They allow both static and time-varying covariates to be included in the analysis. Through simulations, we illustrated that retrospective association tests are robust to ascertainment and other types of phenotype model misspecification, and gain power over previous association methods. We applied L-BRAT and RGMMAT to a genome-wide association analysis of repeated measures of cocaine use in a longitudinal cohort. Pathway analysis implicated association with opioid signaling and axonal guidance signaling pathways. Lastly, we replicated important pathways in an independent cocaine dependence case-control GWAS. Our results illustrate that L-BRAT is able to detect important loci and pathways in a genome scan and to provide insights into genetic architecture of cocaine use.